The present invention relates to RNA interference.
RNA interference (“RNAi”) refers to the silencing of the expression of a gene through the introduction of an RNA duplex into a cell. In RNAi, the RNA duplex is designed such that one strand (the antisense strand) has a region (the antisense region) that is complementary to a region of a target sequence, and the other strand (the sense strand) has a region (the sense region) that is complementary to the antisense strand. In mammals, RNAi requires the use of a small interfering RNA molecule (“siRNA”) that contains both an antisense region and a sense region. Use of longer molecules in mammals results in the undesirable interferon response.
One problem with applying RNAi techniques is that an siRNA that is directed against one particular target may silence another gene. This is referred to as an “off-target effect,” which has been observed to result in 1.5 to 5-fold changes in the expression of dozens to hundreds of genes by either transcript degradation or translation attenuation mechanisms. Off-target effects can occur from either the sense strand or the antisense strand and can occur when as few as eleven base pairs of complementarity exists between the siRNA and target. Jackson et al., (2003) “Expression profiling reveals off-target gene regulation by RNA,” Nat. Biotechnol. 21, 635-7.
Off-target gene silencing can present a significant challenge in the interpretation of large-scale RNAi screens for gene function and the identification and the use of optimal lead components for therapeutic applications. At one time, it was believed that off-target effects were due to overall identity of either strand of an siRNA duplex and a sequence other than the target. However, the inventors have determined that overall identity, i.e., based on all or most of the nucleotides in either the sense and/or antisense region being the same as or complementary to a region of a gene that is not being targeted, cannot very well predict off-target effects, except for near perfect matches.
One solution known to persons of ordinary skill for reducing off-target effects has been to use modifications of nucleotides at select positions within the duplex. Examples of these modifications are described in PCT application, PCT/US2005/011008, publication number WO 2005/097992 A2. However, modifications are not effective on all siRNA, can be expensive, and are not applicable to DNA-based RNAi (i.e. vector driven RNAi).
Further, when running an experiment with a given or candidate siRNA there is a challenge of determining whether any particular phenotype that is observed is due to silencing of a target gene or to an off-target effect. The present invention is directed to this challenge.
The present invention is directed toward determining whether a phenotype is due to an off-target effect in RNAi mediated gene-silencing applications. Additionally, through the use of the methods, libraries and computer program products of the present invention, a person of ordinary skill can reduce the likelihood that an siRNA that is selected will have undesirable levels of off-target effects and determine whether an siRNA induced phenotype is due to an off-target effect or silencing of a target gene.
According to a first embodiment, the present invention provides a method for selecting an siRNA for gene silencing in humans, said method comprising: (a) selecting a target gene, wherein the target gene comprises a target sequence; (b) selecting a candidate siRNA, wherein said candidate siRNA comprises 18-25 nucleotide base pairs that form a duplex comprised of an antisense region and a sense region and said antisense region of said candidate siRNA is at least 80% complementary to a region of said target sequence; (c) comparing a sequence of the nucleotides at positions 2-7 of said antisense region of said candidate siRNA to a dataset wherein said dataset comprises the nucleotide sequences of the 3′ UTR regions (3′ untranslated regions) of a set of human RNA sequences; (d) optionally comparing a sequence of the nucleotides at positions 2-7 of said sense region of said candidate siRNA to said dataset; and (e) selecting said candidate siRNA as an siRNA for gene silencing, if said sequence of the nucleotides at positions 2-7 of said antisense region are 100% complementary to sequences within fewer than 2000 distinct 3′ UTRs of mRNA within said dataset and optionally the nucleotides at positions 2-7 of said sense region are 100% complementary to sequences within fewer than 2000 distinct 3′ UTR regions of mRNA within the dataset.
Two thousand distinct 3′ UTRs represents approximately 8.5% of the 23,500 known human NM 3′ UTR sequences (in Refseq 15). As databases change in size and differ across organisms it may be useful to set the limit as 5%-15% of the known sequences in a given dataset. Preferably for any organism considered, there are at least 5,000, and more preferably at least 10,000 known sequences in a dataset when the method is applied. For humans it was observed that based on the known number of sequences, the set of seeds that appear in fewer than 2000 distinct 3′ UTRs excludes essentially all of the seed sequences that do not contain the CG dinucleotide. Accordingly, although there may be more than 2000 distinct 3′ UTRs that contain certain seeds with the CG dinucleotide, there are substantially no seeds that appear in fewer than 2000 distinct 3′ UTRs that do not contain this dinucleotide.
Positions 2-7 may be referred to as a hexamer sequence. Alternatively, one may focus on positions 2-8, which may be referred to as a heptamer sequence. The nucleotide sequence of the siRNA that is complementary to the 3′ UTR may be referred to as a “seed sequence,” regardless of whether positions 2-7 or 2-8 of the sense or antisense strand. The siRNA that is selected for gene silencing may be introduced into a cell and used to silence the target gene while causing a relatively low level of off-target effects. When performing the above-described method, one may start with one candidate siRNA, a plurality of siRNAs, or all possible siRNAs that contain antisense regions that are complementary to a region of a target sequence. Preferably the antisense region is at least 80% complementary to a region of the target sequence. In some embodiments it is at least 90% complementary to a region of the target sequence. In some embodiments it is 100% complementary to a region of the target sequence.
In a second embodiment, the present invention provides a method for converting an siRNA having desirable silencing properties, yet undesirable off-targeting effects, into an siRNA that retains the silencing properties (or has a functionality that is decreased by no more than 10%, more preferably no more than 5% and most preferably no more than 3%), yet has the lower levels of off-target effects. The method comprises comparing the sequence of the seed of the siRNA with a database comprising low frequency seed complements (or 3′ UTRs that may be searched according to the frequency of sequences that are six or seven bases in length) and identifying one or more single nucleotide changes that could be incorporated into the seed sequence of the siRNA such that the seed sequence is converted to a sequence with a low seed frequency complement without losing silencing activity. Unless otherwise specified, a low frequency seed complement is a sequence that appears in fewer than 2000 distinct human 3′ UTRs. A sequence that appears more than one time in a 3′ UTR for a given mRNA sequence is counted as only a single occurrence for the purpose of the present invention. The aforementioned silencing activity could be determined empirically and/or predicted through rational design criteria as described below.
In a third embodiment, the present invention provides a method of designing a library of siRNA sequences. The method comprises collecting siRNA sequences of at least 100 siRNAs that target at least 25 different genes, wherein said siRNA sequences comprise 18-25 bases, and at least 25% of the siRNA sequences have a hexamer sequence at positions 2-7 of an antisense sequence selected from reverse complement of the sequences of the group consisting of the sequences in Table V below.
The library could in its simplest form be created by identifying a set of candidate siRNA for a plurality of target sequences, and manually typing them into a computer database such that on average at least one of every four siRNAs that are input contains a seed sequence that is the reverse complement of a sequence identified in Table V. Preferably the siRNA within the library all have a selected level of functionality, which may for example be determined by trial and error or may be predicted to be among the most functional through bioinformatics techniques such as those described in U.S. Ser. No. 10/714,333 or PCT/US04/14885. When the library contains both siRNA with seed sequences that are the reverse complement of those within Table V and siRNA with seed sequences that are not the reverse complement of those within Table V, preferably the siRNA that have seed sequences that are the reverse complement of the hexamers in Table V are denoted or otherwise tagged as containing such a sequence for easy identification by a user or computer program.
In a fourth embodiment, the present invention provides a library of siRNA sequences, said library comprising a collection of siRNA sequences of at least 100 siRNAs that target at least 25 different genes, wherein said siRNA sequences comprise 18-25 bases, and at least 25% of the siRNA sequences have a hexamer sequence at positions 2-7 of an antisense sequence selected from the group consisting of the reverse complement of the sequences in Table V below. This library may be populated through the entry of data into an appropriate computer program. As persons of ordinary skill are aware, the computer program will include code for receiving data corresponding to nucleic acid sequences and for searching among this type of data. Preferably, the library also contains a means to differentiate between ORF, and untranslated sequences, (e.g., 5′ UTR and 3′ UTR). Further, although positions 2-7 of the antisense strand are referenced above, this information is understood to refer implicitly to positions 13-18 of the opposite strand in a 19-mer (or corresponding positions in a strand of a different length e.g., positions 17-22 in a 23-mer, positions 19-24 in a 25-mer).
In a fifth embodiment, the preset invention provides a computer program product for use in conjunction with a computer system, the computer program product comprising a computer readable storage medium and a computer program mechanism embedded therein, the computer program mechanism comprising: (a) an input module, wherein said input module permits a user to identify a target sequence; (b) a database mining module, wherein said database mining module is coupled to said input module and is capable of searching an siRNA database comprised of at least 100 siRNA sequences that target at least 25 different genes, wherein each of said siRNA sequences comprises 18-25 bases; and (c) an output module, wherein said output module is coupled to said database mining module and said output module is capable of providing to said user an identification of one or more siRNA sequences from said database where each siRNA that is identified comprises an antisense sequence that is at least 80% complementary to a region of said target sequence and at least 25% of the siRNA sequences identified from said database have a hexamer sequence at positions 2-7 of said antisense sequence selected from the group consisting of the reverse complement of sequences in Table V below. In some embodiments, at least 25% of the siRNA also have a hexamer sequence at positions 2-7 of the sense sequence selected from the group consisting of the reverse complement of sequences in Table V.
In a sixth embodiment, the present invention provides a method of determining whether a phenotype observed with a given siRNA for a target gene in an RNA interference experiment is target specific or is a false positive result. The method comprises: (a) introducing the given siRNA into a first target cell, wherein said given siRNA comprises a sense region and an antisense region, each of which is 18-25 nucleotides in length; (b) measuring said phenotype in said first target cell; (c) introducing a control siRNA into a second target cell, wherein said control siRNA comprises a sense region and an antisense region, each of which is 18-25 nucleotides in length, wherein positions 2-7 of the antisense region of the control siRNA form the same nucleotide sequence as that of positions 2-7 of the antisense region of the given siRNA, wherein the positions 2-7 are counted relative to the 5′ terminus of the antisense regions of the given siRNA and control siRNA, and the rest of the control sequence is scaffold; (d) measuring said phenotype in said second target cell after (c); and (e) comparing the phenotype in said first target cell with the phenotype in said second target cell, whereby, if the phenotype in said first target cell is similar (i.e., both results score as “positive” for a given phenotype in an assay as judged by any one or a number of art accepted statistical and non-statistical methods) to that observed in said second target cell, the phenotype observed in said first target cell is determined to be a false positive result.
In a seventh embodiment, the present invention provides a library of siRNA molecules (this is also referred to as a control siRNA library or seed library), wherein said library comprises a collection of at least 25 siRNAs, wherein each siRNA comprises and antisense region that is 18-25 nucleotides in length, wherein positions 2-7 or 2-8 of the antisense region of each of said siRNA sequences comprises a unique sequence of six or seven contiguous nucleotides and a constant sequence at all other positions of the antisense region.
In an eighth embodiment, the present invention provides a method for constructing a control siRNA library, wherein said library comprises a collection of at least twenty-five siRNAs, wherein each siRNA comprises a sense region and an antisense region, and each of the sense and antisense region is 18-25 nucleotides in length. The method comprises: creating a list of said at least twenty-five siRNA sequences, wherein each of said at least twenty-five sequences comprises a unique sequence of six contiguous nucleotides at positions 2-7 of said antisense region and a constant sequence at all other positions other than the 2-7 positions, wherein the constant sequence forms a neutral scaffolding sequence. A library is to comprise both sense and antisense regions even if only one is recited, because through standard Watson-Crick bases pairing, information about one strand (or region) will provide information about the other. If only one strand is recited, in some embodiments one will assume 100% complementarity between the antisense and sense regions.
The present invention provides methods for reducing off-target effects during gene silencing and methods for selecting siRNA for use in these applications. The present invention also provides libraries and computer program products that assist in increasing the likelihood that an siRNA will have reduced off-target effects and/or provide means for determining whether an observed phenotype is due to an off-target effect.
The inventors have discovered that the number of off-targets generated by an siRNA can be limited by choosing an siRNA that has a sense and/or antisense strand with seed sequences that is/are complementary to the 3′UTR of a limited number of genes in the target genome. As the frequency at which a seed match appears in the population of 3′ UTRs of a genome is predictive of the number of off-targets, it is possible to select for siRNAs that have fewer off-targets based on their seed region.
To that end, according to a first embodiment the present invention comprises a method for selecting an siRNA for gene silencing in a human cell. The method comprises: (a) selecting a target gene, wherein the target gene comprises a target sequence; (b) selecting a candidate siRNA, wherein said candidate siRNA comprises 18-25 nucleotide base pairs that form a duplex comprised of an antisense region and a sense region and said antisense region of said candidate siRNA is at least 80% complementary to said target sequence; (c) comparing a sequence of the nucleotides at positions 2-7 of said antisense region of said candidate siRNA to a dataset wherein said dataset comprises the nucleotide sequences of the 3′ UTRs of a set of human RNA sequences or a data set that is comprised of the frequency of all of the hexamers in the 3′UTR transcriptome; (d) optionally, comparing a sequence of the nucleotides at positions 2-7 of said sense region of said candidate siRNA to said dataset; and (e) selecting said candidate siRNA as an siRNA for gene silencing, if said sequence of the nucleotides at positions 2-7 of said antisense region (and optionally of said sense region) are complementary to sequences that appear in the 3′ UTRs of fewer than 2000 distinct mRNA. Once selected, the sequence may be displayed to a user in for example printed form or displayed on a computer screen. The sequence may also be stored in an electronic memory device. Additionally, the sequence may also be synthesized, including by either enzymatic or chemical means to form an siRNA duplex.
A similar method can be devised based on the frequency of heptamer sequences. However, because there are four times as many possible heptamer sequences, each heptamer sequence will occur on average less frequently than each hexamer sequence. Accordingly, one could look to select siRNA that have heptamer sequences at positions 2-8 of the antisense region and optionally the sense seed region that appears in fewer than 500 distinct 3′ UTRs of human mRNA.
One may omit step (d) when employing this method, in which case during step (e), one would only compare the seed sequence within the antisense region to the 3′ UTR regions (i.e., determine the presence of the reverse complement of the seed sequence). Preferably, step (d) is not omitted unless the duplex will be modified (e.g. through chemical modifications) or contain another cause of strand bias that reduces the likelihood that the sense strand can induce RNAi and thus is rendered essentially incapable of generating undesirable levels of off-target effects. Alternatively, as most rational design algorithms select for siRNA that preferentially introduce the antisense strand into RISC, this method can also be used to minimize the contributions that the sense strand seed makes to off-target effects.
The number of distinct 3′ UTRs in which the reverse complement of seed sequences appear that is selected as the cut off for an organism is selected based on the discovery that the appearance of the complement of seed sequences in 3′ UTRs forms a bimodal distribution. As described more fully in example 4 below and
When the 4096 possible hexamer seeds are binned by the number of distinct human NM 3′ UTRs in which their complements appear, the resulting histogram shows a clear bimodal distribution. The sharp secondary peak at the left of the histogram represents a distinct population of 3′ UTRs with low frequency seed complement. This low frequency may be due to the ubiquitous presence of the CG dinucleotide in these seeds, as the CG dinucletoide is rare in mammals. For humans, the cut off frequency between the two nodes is located at approximately 2000 distinct 3′ UTRs (see
For the rat, this point is approximately 600 for known sequences (see
With respect to implementing the present invention, and as persons skilled in the art are aware, if one assumes 100% complementarity between the sense and antisense strands and one knows the length of the duplex, by examining one strand, information is implicitly provided about the other strand. Thus in a 20-mer duplex, information about positions 2-7 of the antisense strand may be learned by focusing on positions 14-19 of the sense strand.
The Datasets
The terms “dataset” and “database” are used interchangeably and refer to sets or libraries of sequences. The sequences of a database can represent the total collection of e.g., 3‘UTRs of an organism’s genome, or expressed 3′ UTRs for e.g. a particular cell type. Accordingly, databases include but are not limited to those that contain the complete or cell specific mRNA sequences or 3′ UTR sequences e.g., GenBank or Pacdb (http://harlequin.jax.org/pacdb/), or datasets that comprise the frequency of all complements of hexamers or heptamers in the 3′UTR of the transcriptome of the target cell or organism. Such databases can be used to select targets and candidate siRNAs. Additionally, cDNA databases preferably generated using poly-dT primers can be used to select targets and candidate siRNAs. Alternatively or additionally, databases may compromise siRNA sequences. These sequences may be defined by parameters that include but are not limited to length, target sequences, species and predicted or empirical functionality. The siRNA sequences may also have data associated with them that identify gene(s) that they target.
The data may be stored on relational databases or file based databases. Examples of relational databases include but are not limited to Sequel Server, Oracle, and MySeql. An example of a file-based database includes but is not limited to File Maker Pro.
The Target Gene
A “target gene” is any gene that one wishes to silence. As persons skilled in the art are aware, typically siRNAs silence a target gene by becoming associated with RISC (the RNA Induced Silencing Complex) and then cleaving or inhibiting the translation of the target gene messenger RNA (“mRNA”). The mRNA comprises both a coding sequence, which will be translated into a protein or polypeptide, and a 3′ UTR (3′ untranslated region). The mRNA may contain other areas as well, including a 5′ UTR, and/or a tail (e.g., poly A tail). The target gene may be selected based on the desire to study or to knockdown (i.e., reduce expression of) that gene. The “target sequence” is, unless otherwise specified, a portion of the mRNA that codes for a protein. The phrases “target specific effect,” “target-specific gene knockdown” and “target specific” as used herein mean a measurable effect (e.g., a decrease in target mRNA levels, protein levels, or particular phenotype) that is associated with RISC-mediated cleavage of said mRNA. This is to be distinguished from an off-target effect, which is generally: (1) unintended; and (2) mediated by complementarity between the seed region of an siRNA and e.g., a sequence in the 3′UTR of the unintended target gene.
The siRNA
After a gene is selected, at least one candidate (also referred to as a “given”) siRNA is examined, and preferably a plurality of candidate siRNAs is examined. An siRNA is a short interfering ribonucleic acid, that unless otherwise specified contains a sense region of 18-25 and antisense region of 18-25. The antisense region and the sense region may be at least 80% complementary to each other. The antisense region and the sense region may be at least 90% complementary to each other. Unless otherwise specified, they are assumed to be 100% complementary to each other. In addition to an antisense region and a sense region, an siRNA may have one or more overhangs of up to six bases on any, a plurality, all or none of the 3′ and 5′ ends of the sense and antisense regions. Further, unless otherwise specified, within the definition of an siRNA are shRNAs.
When working in mammals such as humans, chimpanzees, rats, mice, horses, sheep, goats, cows, dogs, cats, fugu, etc., preferably each of the antisense region and the sense region of the siRNA comprises 18-25 bases, more preferably 19-25 bases, even more preferably 19-24 bases and most preferably 19-23 bases. Preferably the antisense region is at least 80% complementary to a region of the target sequence. In some embodiments, it is at least 90% complementary to a region of the target sequence. In some embodiments, it is at least 95% complementary to a region of the target sequence. In other embodiment it is 100% complementary to a region of the target sequence. Unless otherwise specified, the antisense region and the region of the target sequence are presumed to be 100% complementary to each other.
The base pairs of an siRNA will form a duplex comprised of an antisense region and a sense region. A candidate siRNA may be comprised of either two separate strands, one of which comprises the antisense region (which may form the entire or be part of the antisense strand) and the other of which comprises the sense region (which may form the entire or be part of the sense strand). The candidate siRNA may also comprise one long strand, such as a hairpin siRNA. Alternatively, the candidate siRNA may comprise a fractured or nicked hairpin that is a duplex comprised of two strands, one of which contains all of the sense region and part of the antisense region, while the other strand contains part of the antisense region. Similarly, a fractured or nicked hairpin may be a duplex comprised of two strands, one of which contains all of the antisense region and part of the sense region, while the other strand comprises part of the sense region. These types of hairpin molecules are also described in pending U.S. patent application Ser. No. 11/390,829, which was filed on Mar. 28, 2006 and published as US 2006-0223777 A1 on Oct. 5, 2006.
The candidate siRNA may have blunt ends or overhangs on either or both of the 5′ or 3′ ends on either or both strands. If any overhangs are present, preferably they will be 1-6 base pairs in length and on the 3′ end of either or both of the antisense strand or sense strand. More preferably, the overhangs will be 2 base pairs in length on the 3′ end of the antisense or sense strand. If the siRNA is a hairpin or fractured hairpin molecule, it will also contain a loop structure.
The candidate siRNA may have modifications, such as 5′ phosphate groups, modifications of the 2′ carbon of the ribose sugars, and internucleotide modifications. Exemplary modifications include 2′-O-alkyl modifications (e.g., 2′-O-methyl, 2′-O-ethyl, 2′-O-propyl, 2′-O-isoproyl, 2′-O-butyl), 2′fluoro modifications, 2′ orthoester modifications, and internucleotide thio modifications. The modifications may be included to increase stability and/or specificity.
Modifications can be added to siRNA to enable users: (1) to apply the invention to one strand; or (2) to enhance the efficiency of the invention. As described in USPTO patent application Ser. No. 11/019,831, publication no. US2005-0223427A1 chemical modifications can be added to enhance specificity. Thus, for example, addition of a 5′ phosphate group on the first antisense nucleotide, and 2′ O-alkyl modifications (e.g., 2′ O-methyl) on the first sense nucleotide and the second sense nucleotide eliminate the ability of the sense strand to enter RISC, and thus would allow users to confine the method of the invention to the antisense strand.
Alternatively, the method of the invention can be applied to both strands to identify siRNA with desirable traits, and subsequently modifications can be added to both strands (e.g., (1) a 5′ phosphate group on the first antisense nucleotide, and 2′ O-alkyl modifications (e.g., 2′ O-methyl) on the first 5′ sense nucleotide, the second 5′ sense nucleotide, the first 5′ antisense nucleotide and the second 5′ antisense nucleotide; or (2) a 5′ phosphate group on the first 5′ antisense nucleotide, and 2′ O-alkyl modifications (e.g., 2′ O-methyl) of the first 5′ sense nucleotide, the second 5′ sense nucleotide and the second 5′ antisense nucleotide) to minimize off-targets further. When modifications are present, all nucleotides that are not specifically identified as having a modification are preferably unmodified, i.e., they have 2′OH groups on their ribose sugars. Thus, the presence of modifications such as 2′ modifications on one or both strands does not preclude application of the current invention. In fact, because certain modifications may reduce off-target effects, but not to the degree desired, in some instances it is advantageous to apply the current invention to both strands of a duplex regardless of whether there are any chemical modifications or other bases for strand bias.
The phrase “first 5′ sense nucleotide” refers to the 5′ most nucleotide of the sense region, and thus this nucleotide would be part of the duplex formed with the antisense region. The phrase “second 5′ sense nucleotide” refers to the next 5′ most nucleotide of the sense region. The second 5′ sense nucleotide is immediately adjacent to and downstream (i.e. 3′) of the first 5′ sense nucleotide, and thus would also be part of the duplex formed. The phrase “first 5′ antisense nucleotide” refers to the 5′ most nucleotide of the antisense region. The phrase “second 5′ antisense nucleotide” refers to the next 5′ most nucleotide of the antisense region. The second 5′ antisense nucleotide is immediately adjacent to and downstream of the first 5′ antisense nucleotide. The first 5′ antisense nucleotide and second 5′ antisense nucleotide are also each part of the duplex formed with the sense region. Thus, any 5′ overhangs do not affect the definition of the aforementioned first or second 5′ nucleotides.
The nucleotides within each region may also be referred to by their positions relative to the 5′ terminus of that region. Thus, the first 5′ antisense nucleotide is located at position 1 of the antisense region, the second 5′ antisense nucleotide is located at position 2 of that region, the third 5′ antisense nucleotide is located at position 3 of that region, the fourth 5′ antisense nucleotide is located at position 4 of that region, the fifth 5′ antisense nucleotide is located at position 5 of that region, etc. A similar convention can be used to identify the nucleotides of the sense region; however, note that in a duplex of 19 base pairs, position 1 of the sense region will appear opposite position 19 of the antisense region. Unless otherwise specified the hexamer and heptamer sequences that are examined in the context of the present invention refer to positions 2-7 and 2-8, respectively of the antisense and/or sense regions of the siRNA.
Previous investigations known to persons of ordinary skill in the art have suggested that off-target effects could be eliminated by minimizing the overall levels of complementarity between an siRNA and unintended targets in the genome of interest. The inventors have demonstrated that this technique is not viable (see Birmingham et al., (2006) “3′ UTR seed matches, but not overall identity, are associated with RNAi off-targets” Nature Methods 3:199-204) and instead, have identified key parameters that allow RNAi users to minimize off-target effects. First, as shown in Example 1, it was observed that the 3′ UTR of off-targeted genes frequently have one or more sequences that are the reverse complement of the seed sequence of an siRNA. Second, as shown in Example 2, the inventors observed that the frequency at which all hexamers and/or heptamers appear in the 3′ UTR sequences of any given genome (e.g. human, mouse, and rat genomes) varies considerably. It was also observed that an association exists between the number of off-targets generated by a particular siRNA, and the frequency at which the reverse complement of the seed sequence of the siRNA appears in the 3′ UTRs of the genome. Based on these observations, the present inventors developed a method for minimizing off-target effects described herein and methods for distinguishing whether a phenotype is due to silencing of a targeted gene or an off-target effect.
When seeking to reduce off-target effects, preferably one focuses on positions 2-7 of the antisense region and/or sense region or positions 2-8 of the antisense region and/or sense region of a candidate siRNA. In some embodiments, it is preferable to consider both strands because either strand could in theory generate an off-target effect. Focusing on a smaller number of positions may lead to false positive matches and focusing on a greater number of positions may lead to false negative results.
As noted above, according to one embodiment of the present invention, one examines positions 2-7 or 2-8 of the antisense region and/or positions 2-7 or 2-8 of the sense region of a candidate siRNA and compares the sequence of the nucleotides located at those positions to the dataset containing sequences from the 3′ UTRs of mRNA of for example, a genome (e.g. a human genome 3′ UTR dataset or other mammalian or other organism's 3′ UTR dataset) to determine whether complementary exists in one or more instances. In some embodiments, preferably, the dataset comprises the 3′ UTRs of at least distinct 1500 mRNA sequences, more preferably of at least 2000 distinct mRNA sequences, and even more preferably of at least 3000 distinct mRNA sequences. In some embodiments, the 3′ UTR regions of all known mRNAs for a species or cell type are within the dataset (e.g. HeLa cells, or MCF7 cells). Preferably, the dataset is also species specific. In some embodiments, when trying to reduce off-target effects in cells expressing human genes, the dataset comprises a sufficiently large set of expressed 3′ UTR regions of human mRNA, if not all known such regions. Alternatively, the data set might be composed of all of the seed complements for a particular cell type, tissue, or organism, and a listing of their frequencies.
After one examines positions 2-7 or positions 2-8 of the antisense region and/or the sense region of a candidate siRNA or collection of siRNA, one may select desirable siRNA based on the frequency of the seed matches in (i.e. instances of complementarity to) the distinct 3′ UTR of e.g. the mRNA dataset. siRNA, for example, can be selected on the basis of having seed sequences that are complementary to sequences in fewer than about 2000 distinct 3′ UTRs, more preferably fewer than about 1500, even more preferably, fewer than about 1000 and even most preferably, fewer than about 500 sequences in 3′ UTR regions. Note that a sequence may appear two or more times within a 3′ UTR of a given gene. In these cases each additional occurrence would not be considered an additional match.
Although not wishing to be bound by any one theory, it is postulated that the advantage of using siRNA that have low seed complement frequencies in the 3′ UTR regions is due to the relatively limited amount of RISC in a cell. RISC is an integral part of gene silencing in mammals, and RISC may be guided to a target by at least two means. First, RISC may be guided to a target when there is full complementarity of a region of the siRNA to the target sequence, typically a region of at least 18 nucleotides. Second, RISC may be guided to another RNA molecule when there is complementarity between positions 2-7 or 2-8 of the antisense region or positions 2-7 or 2-8 of the sense region of the siRNA and a sequence in the 3′ UTR of another molecule.
There are 4096 (46) different sequences for the six nucleotides from positions 2-7, and 16,384 (47) different sequences for the seven nucleotides from positions 2-8 assuming canonical bases, i.e., A, C, G, U. Thus, the method for comparing the candidate siRNA to a dataset comprising 3′ UTRs may be performed most easily by a computer algorithm. The use of computer algorithms to manipulate and to select nucleotide sequences is well known to persons of ordinary skill in the art.
The dataset could be organized by inputting all or a sufficiently large set of mRNA, including their 3′ UTRs. Then one, a plurality, or all candidate siRNAs of a given size or multiple sizes could be compared against the dataset to determine the number of times that the antisense seed sequence and/or the sense seed sequence are complementary to 3′ UTR sequences in the dataset. One could weed out siRNAs that do not have seeds with low frequency seed complements. Alternatively, one could create a dataset of distinct 3′ UTRs, search for the number of distinct 3′UTRs that contain each 6 or 7-mers repeat then develop a database that contains each hexamer or heptamer sequence and the frequency at which it appears in the 3′UTR transcriptome.
The result of the frequency of the 1081 least frequent hexamers based on human 3′ UTRs in RefSeq Version 17 from the NCBI database is identified in Table V. The seed sequences of the candidate siRNA could, for example, then be compared against this set of information to look for complementary sequences and thus determine the likelihood of off-target effects.
The datasets of the siRNAs of the present invention may be organized into specific libraries. For example, one may create a library of at least 100 different siRNAs that target at least 25 different genes (e.g., an average of four siRNA per target) where at least 25% of the siRNA have a seed sequence that is the complement of a sequence selected from Table V. Preferably there are at least 200 different siRNA, more preferably at least 500 different siRNA, even more preferably at least 1000 different siRNA, even more preferably at least 2000 different siRNA, even more preferably at least 5000 different siRNA. Further, preferably the library contains siRNA that target at least 50 different genes, more preferably at least 100 different genes, even more preferably at least 200 different genes, even more preferably at least 400 different genes, even more preferably at least 500 different genes, and even more preferably at least 1000 different genes. A more comprehensive library would contain siRNA that target the entire genome. For example, such a library may contain 100,000 siRNAs for about 25,000 different genes (four siRNAs per gene).
In some embodiments, preferably at least 40%, more preferably at least 50%, even more preferably at least 80%, even more preferably at least 90% and most preferably 100% of the siRNA in a particular collection have a seed sequence that is the reverse complement of a sequence selected from Table V.
The method for selecting siRNA of the present invention may be used in combination with methods for selecting siRNA based on rational design to increase functionality. Rational design is, in simplest terms, the application of a proven set of criteria that enhance the probability of identifying a functional or hyperfunctional siRNA. These methods are for example described in commonly owned WO 2004/045543 A2, published on Jun. 3, 2004, U.S. Patent Publication No. 2005-0255487 A1, published on Nov. 17, 2005, and WO 2006/006948 A2 published on Jan. 19, 2006 the teachings of which are incorporated by reference herein. When selecting siRNA for the aforementioned libraries, one may apply rational design criteria to a set of candidate siRNAs, and then weed out some or all sequences that do not meet the aforementioned seed criteria. Thus, in these circumstances, the seed criteria may be a filter applied to rational design criteria. Alternatively, one could weed out some or all sequences that do not satisfy the seed criteria, and then apply rational design criteria.
Combining the methods of the invention with siRNA selected by rational design as described above may allow users to simplify the application of the method by focusing on the seed sequence of the antisense strand. Rationally designed siRNA are (in part) selected on the basis that the antisense strand of the duplex (i.e. the strand that is complementary to the desired target) is preferentially loaded into RISC. For that reason, off-targets of rationally designed siRNA are predominantly the result of annealing of the seed region of the antisense strand with the sequences in the 3′ UTR of the off-targeted gene. Therefore, in cases where rationally designed siRNA having an antisense strand bias are being used, it is possible to confine the method of the invention to the antisense strand alone, and ignore possible off-target contributions by the sense strand.
The siRNA selected according to the present invention may be used in both in vitro and in vivo applications, in for example, connection with the introduction of siRNA into mammalian cells.
The siRNA used in connection with the present invention may be synthesized and introduced into a cell. Methods for synthesizing siRNA of desired sequences are well known to persons of ordinary skill in the art. These methods include but are not limited to generating duplexes of two separate strands and unimolecular molecules that form duplexes by chemical synthesis, enzymatic synthesis, or expression vectors of siRNA or shRNA.
In another embodiment, the invention provides a method for converting an siRNA having desirable silencing properties, yet undesirable off-targeting effects, into an siRNA that retains the silencing properties, yet has fewer off-targets. The method comprises comparing the sequence of the seed of the siRNA(s) with a database comprising low frequency seed complements and identifying one or more single nucleotide changes that could be incorporated into the seed sequence of the siRNA such that the frequency of the seed complement is converted from a moderate or high frequency, to a low frequency, without losing silencing activity. In one non-limiting example of this method, highly functional siRNA containing an sense seed of 5′-AGGCCG, 5′-ACCCCG, or 5′-ACGCCT (seed complement frequencies of 2376, 2198, and 2001 based on all human NM 3′ UTRs derived from NCBI RefSeq 15) can be converted to a low frequency seed complement (5′-ACGCCG, 472 appearances) by altering a single nucleotide, thus generating an siRNA with a seed that has a low frequency seed complement. A “low frequency seed complement” refers to a sequence of bases whose complement appears relatively infrequently in the 3′ UTR region of mRNAs, e.g., appears in equal to or fewer than about 2000 distinct 3′ UTR regions, more preferably fewer than about 1500 3′ UTR regions, even more preferably, fewer than about 1000 3′ UTR regions, and most preferably fewer than about 500 times in 3′ UTRs. By changing a based within the siRNA, the antisense region of the siRNA may have a lower degree of complementarity with the target. In some embodiments, when the nucleotide of the antisense region is changed, the corresponding nucleotide of the sense region is changed as well.
The present invention also provides a method for designing a library of siRNA sequences. By having a library of siRNA sequences, a person of ordinary skill has readily available a set of siRNAs that has been pre-screened to, for example, have a reduced level of off-target effects. In one embodiment the library contains sequences of at least 100 siRNAs that target at least 25 different genes. Larger databases such as those described above are also within this embodiment.
The sequences within the library may be for one or both strands of an siRNA duplex that is 18-25 base pairs in length. Because of standard AU, GC base pairing it is not necessary to have the code for both strands in the database. When a library has a plurality of siRNA for a given gene, a user may use individual sequences from the plurality or use them in a pool. Thus, by way of example, a user may select a highly functional siRNA such as that determined by Formula X of PCT/US04/14885 and filter those sequences by applying a low frequency seed complement criterion, which may for example, be any siRNA with a seed sequence that is the reverse complement of a sequence that is identified in Table V, or it may be an siRNA with the lowest seed complement frequency for the target, or it may be an siRNA with the lowest seed complement frequency that is among the siRNAs that have the two, three, four, five, six, seven, eight, nine, or ten highest predicted functionalities (or empirical functionalities, i.e., gene silencing capabilities if known). Alternatively, one may use pools of two, three, four, five, six etc., siRNAs that have low if not the lowest seed complement frequencies. Still further one could combine pools of two, three, four, five, six, etc. siRNAs for a target wherein within each pool one or more are selected based on functionality and one or more are selected based on seed complement frequency.
In Table V below is a list that represents hexamer nucleotide sequences that occur at least once in fewer than 2000 distinct known human NM 3′ UTRs. There are 1081 hexamer sequences in the list. As noted above, the 4096 possible hexamers are not uniformly distributed in human 3′ UTRs, instead showing a distinct bimodal distribution including a population of low-frequency hexamers (as defined above). The inventors have demonstrated that siRNAs whose seed complements occur infrequently in 3′ UTRs produce significantly fewer off-targets than those whose seed complements occur at higher frequencies. The use of “T” in the table is by convention in most databases. However, it is understood as referring to a Uracil in any RNA sequence, including any siRNA sequence.
Additionally, it may be desirable to create a library with a maximal percentage of siRNA sequences that have low seed frequency complements. Although it may be preferable for most or all sequences to have low seed frequency complements, that is not always practical for a given target gene, and other considerations such as functionality are important to consider. Thus, preferably on average at least one of every four siRNA sequences has a seed that has a low frequency complement, more preferably on average at least two of every four siRNAs have a seed with a low frequency seed complement, even more preferably on average at least three of every four siRNAs have a seed with a low frequency complement. In some embodiments at least one siRNA for each target contains a seed with a low frequency if not the lowest frequency seed complement. Table V identifies the 1081 seed complement sequences that occur in the fewest distinct human 3′ UTRs. Also included in the table under the heading “distinctnmutr3” is the number of 3′ UTRs in which a given low frequency seed complement sequence appears.
Given the presentation of Table V, a person of ordinary skill could create a database by comparing the seed sequences of a plurality of siRNA to the sequences on Table V and inputting those siRNA into a searchable database if those siRNA contain the seeds that have a seed complement frequency below a requisite level. The person of ordinary skill may also include information about the functionality of the siRNA as well as its targets. Preferably, the library is searchable through computer technology and contains a mechanism for linking the sequence data with e.g., target data and/or seed complement frequency.
The libraries of the present invention may, for example, be located on a user's hard drive, a LAN (local area network), a portable memory stick, a CD, the worldwide web or a remote server or otherwise, including storage and communication technologies that are developed in the future.
The computer program products of the present invention could be organized in modules including input modules, database mining modules and output modules that are coupled to one another. In some embodiments, the modules may be one or more of hardware, software or hybrid residing in or distributed among one or more local or remote computers. The modules may be physically separated or together and may each be a logic routine or part of a logic routine that carries out the embodiments disclosed herein. The modules are preferably accessible through the same user interface.
The software of the present invention may, for example, run on an operating system at least as powerful as Windows 2000.
The computer program may be written in any language that allows for the input of a sequence and searching within a dataset for an siRNA that targets the sequence based on complementarity or identity. For example, the computer program product may be in C#, Pearl or LISP. The program may be run on any standard personal computer or network system. Preferably the computer is of sufficient power to quickly mine large datasets, such as those of the present invention, e.g., 2.33 GHz, 256 RAM and 80 Gb.
The input module will thus be accessible to a user through a user interface and permit a user to select a target gene by for example, name, accession number and/or nucleotide sequence. The input module may offer the user the ability to request the format of the output, and the content of the output, e.g., request the lowest frequency seed complement to be output and/or the lowest frequency with a set of the highest functional siRNAs, e.g., the siRNA whose functionality is predicted to the highest by a set of rational design criteria.
The input module may then convert the inputted data into a standard syntax that is sent to the database mining module. The database mining module then searches a database containing a set of siRNA that are either complementary to or similar to a region of the target depending on whether sense or antisense information is input. The database mining module then transmits the result to the output module, which either saves the results and/or displays them on a user interface. The computer program product may be configured such that the database mining module searches within a database that is part of the computer program product, and/or configured to mine a stand alone database.
The computer program product, as well as the library and methods described herein may be used to assist persons of ordinary skill in the art to identify siRNA with reduced off-target effects.
The computer program product may be run on any standard personal computer that has sufficient power capabilities. As persons of ordinary skill in the art are aware, a more powerful computer may be able to manipulate larger amounts of data at a faster rate. Exemplary computers include but are not limited to personal computers currently sold by IBM, Apple, Dell and Gateway.
According to another embodiment, the present invention provides a method of determining whether a phenotype observed in RNAi experiment is target specific or the result or indicative of a false positive. As used herein, the term “phenotype” refers to a qualitative or quantitative characteristic measured by an assay in vitro or in cells such as the expression of one or more proteins and/or other molecules by a cell or cell death. Methods for measuring protein expression or cell death include but are not limited to counting viable cells, cell proliferation counted as cell numbers, translocation of a protein such as c-jun or NFKB, expression of a reporter gene, microarray analysis to obtain a profile of gene expression, Western Blots, cell differentiation, etc.
A “false positive” is when a test or assay wrongly attributes an effect or phenotype to a particular treatment. If an siRNA targeting gene A gives rise to cell death, this result is a false positive if in fact knockdown of gene A can be separated from the cell death phenotype.
The phenotype observed with a seed control siRNA is said to be ‘similar’ to that generated by the test siRNA when the seed control siRNA generates a phenotype that is positive as judged by the same statistical criteria that were used in the assay to identify the test siRNA, e.g., measurement of decreased protein production or determining cell death.
A “phenotype” is a detectable characteristic or appearance.
Under this method, one introduces a given (also referred to as a “candidate”) siRNA into a first target cell. Preferably the siRNA comprises a sense region and an antisense region, each of which is 18-25 nucleotides in length, exclusive of overhangs. Any overhangs may be 0-6 bases and located on the 5′ end and/or 3′ end of the sense and/or antisense regions. In some embodiments, no overhangs are present. Additionally, preferably the antisense region and the sense region are at least 80% complementary to each other. In some embodiments, they are at least 95% complementary to each other. In some embodiments that are 100% complementary to each other.
The target nucleotide sequence may be either a DNA sequence or RNA sequence.
The target cell may be any cell that either exhibits or has the potential to exhibit a particular characteristic such as the expression of a protein of interest. When the effect of an siRNA on a particular phenotype is being measured, a baseline level of that phenotype can be measured in the target cell, which may be referred to as the baseline target cell.
A given or candidate siRNA may be introduced into a first target cell. The first target cell is preferably the same cell type as the cell that was used to determine the baseline value of the phenotype of interest exists under the same conditions, e.g., (same cell density, temperature and protection from or exposure to environmental stimuli). After the given siRNA is introduced, the phenotype of interest is measured.
Additionally a control siRNA is introduced into a second target cell. The phrase “control siRNA” refers to an siRNA that has the same (antisense) seed sequence as the test siRNA, associated with a scaffold. Alternatively, a control siRNA can contain two seed sequences: the first at positions 2-7 or 2-8 on the antisense strand and the second at positions 2-7 or 2-8 on the sense strand. In these instances, the scaffold represents all of those nucleotides that are not associated with the two seed sequences.
As a person of ordinary skill in the art would appreciate, the second seed position reflects the portion of the sense strand that is complementary to positions 13-18 and 12-18 of the antisense strand in a 19-mer duplex. The use of the second seed region may be desirable when both strands of the test siRNA have the potential to enter RISC. Thus, when the control siRNA comprises only the first seed region, the sense region may contain the modifications identified above to prevent the strand comprising the sense region (e.g., a sense strand or the end of a hairpin molecule) from entering the RISC complex. An exemplary modification is a 2′-O-methyl group as positions 1 and 2 of the sense region.
The bases that are not within the seed region of the control siRNA and the complement in the other strand form a neutral scaffolding sequence. Preferably the scaffolding has a similarity of less than 80% to the bases at corresponding positions within the given (also referred to as candidate or test siRNA) siRNA, more preferably less than 60% similarity, even more preferably less than 50% similarity, even more preferably less than 20% similarity. The term “similarity” as used in this paragraph refers to the identity of a particular nucleotide at a particular position within the sense or antisense region. In some embodiments, within the scaffolding, the control and candidate siRNAs contain none of the same bases at the same positions.
In some embodiments the neutral scaffolding is derived from a sequence that has been empirically tested not to have undesirable levels of off-target effects.
In some embodiments, it is preferable to have position 1 of the antisense region be occupied by U.
Exemplary neutral scaffolding sequences are shown below as a sense sequence where N is A, U G, or C. The string of 6 Ns represents a hexamer of choice from the seed control library. Alternatively, 7 Ns (a heptamer) can replace the 6 Ns shown in the sequences below, with the base at sense position 12 (antisense position 8) in the 19-mer changing from A, C, G, or T to N. Thus, under these circumstances, the first nucleotide that is 5′ of the hexamer of Ns is replaced with another N to generate the complement of the seed heptamer. For instance, SEQ. ID NO. 13, 5′ UGGUUUACAUGUNNNNNNA 3′ would appear as SEQ. ID NO. 16: 5′ UGGUUUACAUGNNNNNNNA 3′;
Unless otherwise specified, the antisense sequences are assumed to be 100% complementary to these sense sequences:
SEQ. ID NO. 13 is derived from a “SMART selection” designed siRNA targeting GAPDH. This siRNA is one selected using rational design criteria such as those described in WO 2006/006948 A2. SEQ. ID NOs. 14 and 15 are derived from functional siRNAs targeting GAPDH and PPIB respectively.
The second target cell is preferably the same type of cell as the first target cell and maintained under the same conditions as the first target cell.
After introduction of the control siRNA into the second target cell, the phenotype is measured. This phenotype is compared to the phenotype measured after introduction of the given siRNA into the first target cell. If the phenotype of the first target cell is similar to the phenotype of the second target cell after the introduction of the siRNA, the phenotype observed in the first target cell is determined to be a false positive. A phenotype is considered similar if both phenotypes pass the threshold limit as defined by the assay or are scored as a “hit” as defined by any number of statistical methods that are used to assess assay outputs. Such statistical methods include, but are not limited to B scores and z score.
According to another embodiment, the present invention provides a library of sequences of at least twenty-five siRNA molecules that are 18-25 bases in length. Each duplex in the library comprises either one or two unique sequences and a scaffolding sequence. The one or two unique sequences are located at the positions of the seed sequences in the previous embodiments. When the siRNA comprises one unique region, the unique region is located at positions 2-7 or 2-8 within the antisense region. These positions are counted from the 5′ end of the antisense region. When the siRNA comprises two unique regions, the first unique region is located at positions 2-7 or 2-8 within the antisense region, and the second unique region is located at positions 2-7 or 2-8 of the sense region.
The library of this embodiment may contain at least 25 sequences, at least 50 sequences, at least 100 sequences, at least 200 sequences, at least 300 sequences, at least 500 sequences, at least 750 sequences, at least 1000 sequences, e.g., 1081 sequences, or all possible the number of sequences that correspond to all of the possible combinations of unique sequences. For example, when there is one seed region of six contiguous nucleotides, there are 4096 (46) unique sequences, when there is one region of seven contiguous nucleotides, there are 16384 (47) unique sequences, when there are two seed region of six contiguous nucleotides, there are 16,777,216 (412) unique sequences, when there are two seed regions, one of six contiguous nucleotides and one of seven contiguous nucleotides, there are 67,108,864 (413) unique sequences, and when there are two seed regions both of seven contiguous nucleotides, there are 268,435,456 (414) unique sequences.
In one embodiment the library comprises at least 1081 siRNA sequences, wherein 1081 of the siRNA sequences each comprises a unique sequence selected from the reverse complement of the sequences identified in table V at positions 2-7 of the antisense region and a neutral scaffolding at all other positions.
This library may be stored in a computer readable storage medium such as on a hard drive, CD or floppy disk.
According to another method, the present invention provides a method for constructing a control siRNA library. This library may contain any number of sequences with a unique seed region or unique seed regions as described above, e.g., at least 25 sequences, at least 50 sequences, at least 100 sequences, at least 200 sequences, at least 300 sequences, at least 500 sequences, at least 750 sequences, at least 1000 sequence, etc. The library comprises nucleotide sequences that describe antisense regions. This description may be through recitation of antisense region sequences themselves or recitation of sense region sequences with the understanding the antisense region will have a sequence that is the reverse complement of the sense region sequence. Additionally, the library may or may not identify overhang regions that are ultimately to be used with an siRNA.
Preferably the sequences in the seed control library are 18-25 nucleotides in length.
The method comprises creating a list of the desired number of siRNA sequences, wherein each of the sequences comprises a unique sequence of six contiguous nucleotides at the positions that correspond to positions 2-7 of the antisense region and a constant region at all other positions. In other embodiments, the unique sequence could occupy (i) positions 2-8 of the antisense region; (ii) positions 2-7 of the antisense region and positions 2-7 of the sense region; (iii) positions 2-8 of the antisense region and positions 2-8 of the sense region; (iv) positions 2-8 of the antisense region and positions 2-7 of the sense region; or (v) positions 2-7 of the antisense region and positions 2-8 of the sense region. The listing may for example be stored within the memory or a computer readable storage device.
Each of the elements within any of the aforementioned embodiments may be used in connection with any other embodiment, unless such use is inconsistent with that embodiment.
Having described the invention with a degree of particularity, examples will now be provided. These examples are not intended to and should not be construed to limit the scope of the claims in any way. Although the invention may be more readily understood through reference to the following examples, they are provided by way of illustration and are not intended to limit the present invention unless specified by and in the claims.
siRNA Synthesis. siRNA duplexes targeting human PPIB (NM—000942), MAP2K1 (NM—002755), GAPDH (NM—002046), and PPYLUC (U47295), were synthesized with 3′ UU overhangs using 2′-ACE chemistry. Scaringe, S. A. (2000) “Advanced 5′-silyl-2′-orthoester approach to RNA oligonucleotide synthesis,” Methods Enzymol. 317, 3-18; Scaringe, S. A. (2001) “RNA oligonucleotide synthesis via 5′-silyl-2′-orthoester chemistry,” Methods 23, 206-217; Scaringe, S, and Caruthers, M. H. (1999) U.S. Pat. No. 5,889,136; Scaringe, S, and Caruthers, M. H. (1999) U.S. Pat. No. 6,008,400; Scaringe, S. (2000) U.S. Pat. No. 6,111,086; Scaringe, S. (2003) U.S. Pat. No. 6,590,093.
Transfection. HeLa cells were obtained from ATCC (Manassas, Va.). Cells were grown at 37° C. in a humidified atmosphere with 5% CO2 in DMEM, 10% FBS, and L-Glutamine. All propagation media were further supplemented with penicillin (100 U/mL) and streptomycin (100 μg/mL). For transfection experiments, cells were seeded at 1.0-2.0×104 cells/well in a 96 well plate, 24 hours before the experiment in antibiotic-free media. Cells were transfected with siRNA (100 nM) using Lipofectamine 2000 (0.25 μL/well, Invitrogen) or DharmaFECT 1 (0.20 μL/well, Thermo Fisher, Inc.). For targeting of PPYLUC (U47295), cotransfections of plasmid and siRNA were performed using Lipofectamine 2000 at 0.5 μL/well in 293 cells at 2.5×104 cells/well in a 96 well plate and harvested at 24 hours.
Gene Knockdown and Cell Viability Assay. Twenty-four to seventy-two hours post-transfection, the level of target knockdown was assessed using a branched DNA assay (Genospectra) specific for the target of interest. In all experiments, GAPDH (a housekeeping gene) was used as a reference. When GAPDH was the target gene, PPIB was used as a reference. All experiments were performed in triplicate and error bars represent standard deviation from the mean. For viability studies, 25 μl of AlamarBlue reagent (Trek Diagnostic Systems) was added to each well, and cells were incubated 1-2 h at 37° C., 5% CO2. Absorbance was then read at 570 nm using a 600 nm subtraction. The optical density (OD) is proportional to the number of viable cells in culture when the reading is in the linear range (0.6 to 0.9). Transfections resulting in an OD of ≧80% of control were considered nontoxic.
Microarray Experiments. For each sample, 1 μg of total RNA isolated from siRNA-treated cells was amplified and Cy5-labeled (Cy-5 CTP, Perkin Elmer) using Agilent's Low Input RNA Fluorescent Linear Amplification Kit and hybridized against Cy3 labelled material derived from lipid treated (control) samples. Hybridizations were performed using Agilent's Human 1A (V2) Oligo Microarrays (˜21,000 unique probes) according to the published protocol (750 ng each of Cy-3 and Cy-5 labelled sample loaded onto each array). Slides were washed using 6× and 0.06×SSPE (each with 0.025% N-lauroylsarcosine), dried using Agilent's nonaqueous drying and stabilization solution, and scanned on an Agilent Microarray Scanner (model G2505B). The raw image was processed using Feature Extraction software (v7.5.1). Further analysis was performed using Spotfire Decision Site 7.2 software and the Spotfire Functional Genomics Module. Outlier flagging was not used. Off-targets were identified as genes that were down-regulated by two-fold or more (log ratio of more than −0.3) by a given siRNA in at least one experiment, but were not modulated by other functionally equivalent siRNA targeting the same gene.
Computational Analysis. The Smith-Waterman local algorithm was implemented in C# and augmented to extend alignments along the entire length of the shorter aligned sequence. The implementation also allowed the use of either uniform match rewards/mismatch costs or scoring matrices, and either linear or single affine gap costs.
The first stage of analysis used this implementation to align each strand of 12 siRNAs (including one non-rationally designed siRNA) against all GenBank mRNAs represented on the microarray chip. The 1000 highest percent identity alignments (on either strand) for each siRNA were archived. The archived alignments were analyzed to determine their identity distributions and discover alignments with experimentally off-targeted mRNAs, using the validated dataset of 347 off-targets, including all accession numbers that were sequence-specifically down-regulated by 2-fold or more in at least one biological replicate.
The parameter-testing studies defined twelve scoring matrixes designed to reward complementarity rather than identity. Each scoring matrix was combined with at least one linear gap penalty (designed to allow only one gap at a time) and one single affine gap penalty (designed to allow multiple-gap runs) of varying weights to generate the 30 parameter sets. The dataset of experimental off-targets was limited to include only those 180 that were sequence-specifically down-regulated by approximately 2-fold or more in two biological replicates for the 11 rationally designed siRNAs and had well-annotated coding sequences. A control set was chosen at random from those mRNAs that were not significantly down-regulated by any of the test siRNAs, and assigned to the siRNAs in equal numbers as in the off-target set. For each parameter set, the S-W implementation was used to align each strand of the siRNAs with their off-targets' reversed mRNA (due to the complementary nature of the scoring matrices) and the best 20 alignments were archived; the process was repeated for the control set. Analysis identified the highest percent identity archived alignment for each siRNA/mRNA pair (including both strands) and generated histograms of these highest identity distributions for each dataset under each parameter set. Since all distributions except those for sets 29 and 30 were approximately normal, each off-target/control distribution pair except these two was subjected to a two-tailed T-test to determine whether their means were significantly different. The remaining two were subjected to a chi-squared test for independence. The results of all tests were adjusted using the Bonferroni correction to account for multiple comparisons. The analysis was also conducted for each strand individually.
The seed analysis was performed using a stringent subset of the experimentally validated off-targets including only those 84 with well-annotated UTRs that were sequence-specifically down-regulated by at least 2-fold in both of two biological replicates for 8 siRNAs measured in a single experiment; the control set was correspondingly narrowed. The analysis counted occurrences of exact substrings (identical to positions 13-18 inclusive, hexamer, and 12-18 inclusive, heptamer) of the siRNA sense strand to the 5′ UTR, ORF, and 3′ UTRs of each off-target and control.
A database of experimentally validated off-targeted genes was generated from the expression signatures of HeLa cells transfected with one of twelve different siRNAs (100 nM) targeting three different genes, PPIB, MAP2K1, and GAPDH. Eleven rationally designed siRNAs having a strong antisense (AS) strand bias toward RISC entry and one non-rationally designed siRNA were transfected into cells. Rationally designed siRNAs were selected according to the methods disclosed in U.S. Patent Publication No. 2005/0255487 A1.
Genes that were down-regulated by two-fold or more (i.e. expression of 50% or less as compared to controls) by a given siRNA in one or more biological replicates, but were not modulated by other functionally equivalent siRNA targeting the same gene were designated as off-targets. Expression signatures of cells transfected with the 12 siRNAs identified 347 off-targeted genes. The expression signatures are shown in
Tables IA-IC provide the siRNA sequence, intended target, list of validated off-targets and subsets of sequences that were used in each analysis. Table IA identifies the sequences used. Table IB provides data for the experimental results. Table IC provides the results for use in the sw1, sw2 and the seed analyses. “sw1” identifies the group of validated off-targets that were used to generate
Using the Smith Waterman alignment algorithm, the sense and antisense strands for each siRNA were aligned against the more than 20,000 genes represented on Agilent's Human 1A (V2) Oligo Microarray. Gene Sequences that exhibited ≧79% identity with either the sense or antisense strands were designated as in silico predicted off-targets. Commonly used reward/penalty parameters (a match reward=2, a mismatch penalty=−2, and a linear gap penalty=−3) were employed and a maximum cutoff of 1000 alignments per siRNA was arbitrarily imposed. (Although multiple alignments between a given siRNA and mRNA were recorded, analyses were done using only the best alignment between each pair). Surprisingly, the number of in silico predicted off-targets typically exceeded the number identified by microarray analysis by 1-2 orders of magnitude, regardless of whether alignments of one or both strands were included in the analysis. Thus, comparison of the validated off-target dataset with in silico predicted off-targets showed that identity cutoffs failed to accurately predict off-targeted genes.
Table II demonstrates the discrepancy between the number of validated off-targets for each siRNA and the predicted number of targets using different identity cutoffs. Predicted numbers are based on identity matches between the sense and antisense strand of the siRNA against the GenBank genes represented on Agilent's Human 1A (V2) Oligo Microarray. Table II below demonstrates a false positive rate of over 99% at the 79% identity cutoff. This number of predicted off-targets represented more than one third of the number of mRNAs in the human genome. Moreover, only 23 of the 347 experimentally validated off-targets were identified by in silico methods using this cutoff, which represents a false negative rate of approximately 93%. Higher cutoffs (≧84% and ≧89%) produced similarly poor overlap between experimental and in silico target predictions (7 and 1 commonly identified targets using the 84%, and 89% identity filter, respectively), as well as gross mis-estimations of the number of off-targets (1278 and 54, respectively). Based on these observations, it was concluded that overall sequence identity was a poor predictor of the number and identity of off-targeted genes.
The inventors recognized that alignments are particularly sensitive to the weighting of matches, mismatches, and gaps. With the long term goal of creating a customized S-W parameter set that can distinguish between off-targeted and untargeted populations, individual siRNAs targeting human cyclophilin B (PPIB), firefly luciferase (PPYLUC), and secreted alkaline phosphatase (SEAP) were synthesized in their native state or with one of three base pair mismatches at each of the 19 positions of the duplex (48 variants per siRNA). Subsequently, a systematic single mismatch analysis of siRNA functionality was performed by transfecting each siRNA into HeLa cells and measuring the relative level of target silencing. The results of these experiments are presented in FIGS. 3A-C and demonstrate several points.
First, Ppyr/LUC #5 and ALPPL2#2 studies clearly show that the central region of the duplex (positions 9-12) is particularly sensitive to mismatches. In contrast, duplexes with mismatches at positions 18 and 19 exhibit consistent silencing, suggesting that the strength of base pairing in this region is less critical. Outside of positions 9-12 and 18-19, the inventors observed that identical mismatches at any position could have widely disparate impacts on siRNA performance. Thus, for instance, while an A-G mismatch at position 3 of the Ppyr/LUC #5 has little impact on overall duplex functionality, the same mismatch at the same position in the ALPPL2#2 targeting siRNA dramatically alters silencing efficiency.
Second, G-A and G-G mismatches at position 14 of the ALPPL2 #2 siRNA have little or no effect on functionality, but identical mismatches at the same position in the Ppyr/LUC #5 siRNA result in a loss of activity. These findings suggest that with the exceptions of positions 18 and 19 (which appear to be insensitive to base pair mismatches) the complete sequence plays a role in determining the impact of mismatches, thus preventing the development of clear position-dependent mismatch criteria. Nonetheless, analysis of all mismatches in a position independent manner identifies a decided bias (
These observed biases were incorporated into 30 additional S-W parameter sets to test whether changes in the rewards/costs associated with matches and mismatches could improve the ability to predict off-targeted genes by overall alignment identity. Table III below describes the thirty custom S-W scoring parameters sets tested.
As it is unclear how gaps are tolerated by RNAi, several different gap penalties (both linear and affine) were included in the scoring matrices. Two populations of siRNA/mRNA pairs (180 representing experimentally validated off-target interactions and 180 having no discernable off-target interactions) were analyzed with each of the 30 unique scoring schemes. Analysis of off-targeted and untargeted populations using each of the modified parameter sets failed to distinguish between the two datasets regardless of whether alignments for one or both strands were included. The finding that the distributions of maximum identity in the best alignment for each parameter set for off-targeted and untargeted populations are statistically indistinguishable (p>0.05 after application of Bonferroni correction for multiple comparisons,
Recent studies on microRNA (miRNA) mediated gene modulation have shown that complementary base pairing between the seed sequence and sequences in the 3′ UTR of mRNA is associated with miRNA-mediated gene knockdown. (Lim et al., Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs, Nature 433, 769-73 (2005)). As siRNAs and miRNAs are believed to share some portion of the RNAi machinery, the inventors investigated whether complementarity between the seed sequence of the siRNA and any region of the transcript was associated with off-targeting. To accomplish this, the 5′ UTR, ORF, and 3′ UTR of 84 experimentally determined off-target genes were scanned for exact complementary matches to the antisense seed sequence (hexamer, positions 2-7, and heptamer, positions 2-8) of their respective siRNA. This dataset of siRNAs and their off-targeted genes was then compared to a control group (84 siRNA/mRNAs that shared no off-target interactions) to determine whether seed matches in any of the three regions correlated with off-targeting. For 5′ UTR and ORF sequences, the frequency at which one or more hexamer seed matches were present in the experimental and control groups was statistically indistinguishable (at the p>0.05 level using the chi squared test for independence, frequencies were 2.3% and 5.9% for the 5″ UTR, 30.9% and 23.8% for ORF sequences, respectively). In contrast, the incidence at which one or more hexamer matches were found in the 3′ UTR of off-targets was nearly 5-fold higher than that observed in the untargeted populations (84.5% in the experimental group, 17.8% in the control group; significant with p<0.001,
Furthermore, the positive predictive value (defined as [true positives]/[true positives+false positives]) of the association between 3′ UTR hexamer seed matches and off-targeted genes increased when multiple matches were required (for two or more 3′ UTR matches: off-targeted genes=29.76%, untargeted genes=3.57%) as shown in Table IV below, for sensitivity, specificity, and positive predictive power of siRNA hexamer and heptamer seed matches.
When four 3′ UTR hexamer seed matches are present, no false positives were detected in this limited sample. As seed matches provide an enhancement over the predictive abilities of blastn and S-W homology based searches, a search tool has been developed to enable identification of all possible human off-targets for any given siRNA based on 3′ UTR hexamer seed matches. The 3′ UTR hexamer identification tool takes the 19 base pair siRNA sense sequence, identifies the corresponding hexamer of the target site, and displays the identity of all genes carrying at least one perfect hexamer seed match in the 3′ UTR. A second column may display a smaller subset of genes that have two or more perfect 3′ UTR seed matches.
The frequency at which heptamer seed matches were observed in the 5′ UTR, ORF, and 3′ UTR of experimental and control groups was similar to those documented for hexamers (heptamer frequency in experimental and control groups: 5′ UTR: 0% and 1.2%; ORF: 16.6% and 9.5%; 3′ UTR: 69.1% and 8.3%) suggesting that the relevant seed sequence may consist of 7 nucleotides (positions 2-8), and the method of the present invention may be applied by focusing on either size region. As was observed with hexamer seed matches, increases in the numbers of 3′ UTR heptamer seed matches were associated with improvements in the specificity of the association. The observed associations remain after 3′ UTR length is controlled for by examining paired off-targeted and non-targeted control 3′ UTRs with lengths equal to within thirty bases (
The work presented here demonstrates that with the exception of instances of near-perfect complementarity, the level of overall complementarity between an siRNA and any given mRNA is not associated with off-target identity. Both S-W and BLAST sequence alignment algorithms grossly overestimate the number of off-targeted genes when common thresholds are employed, suggesting that siRNA designed algorithms employing these methods may be discarding significant numbers of functional siRNAs due to unfounded specificity concerns. Moreover, the overlap between predicted and validated off-targets is minimal (0.2 to 5%) when identity thresholds ranging between ≧79% and ≧89% are employed. In addition, custom S-W parameters informed by base pair mismatch studies fail to produce alignments that distinguish between off-targeted and untargeted populations. These findings reveal that current protocols used to minimize off-target effects (e.g. BLAST and S-W) have little merit aside from eliminating the most obvious off-targets (i.e. sequences that have identical or near-identical target sites).
The sequences of human NM 3′ UTRs for RefSeq Version 17 were down loaded from NCBI (http://www.ncbi.nlm.nih.gov/). Subsequently, a comparison was made between these sequences and all 6 and 7 nt seeds (Lewis, B. P., C. B. Burge and D. P. Bartel. (2005) “Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets,” Cell 120(1):15-20) to determine the frequency at which each possible hexamer/heptamer seed obtain was observed. The results, presented in
When the 4096 possible hexamer seeds are binned by the number of human NM 3′ UTRs in which they appear, the resulting histogram shows a distinct bimodal distribution. The sharp peak at the left of the histogram represents a distinct population of low-frequency seeds. (As shown in
The low seed complement frequency threshold of 2000 distinct 3′ UTRs was arrived at by determining the uppermost boundaries of the rare-seed peak. In other animals (notably rat, in which the number of available NM RefSeq 3′ UTRs is only about ⅓ of that available for human) the 2000 threshold would not apply, but the bimodal distribution is still evident in
Thus, the threshold used for a particular organism (or for the human organism when designing against a later—and therefore larger—RefSeq database) should preferably be redetermined by plotting the above sort of histogram and selecting the upper limit of the rare seed peak. If this is not possible, then a percentage threshold may be applied (although it is not proven that the percentage of seeds in the low frequency peak is completely comparable between organisms); 2000 distinct 3′ UTRs represent approximately 8.5% of the currently known human transcriptome, so a reasonable percentage-based threshold would be to designate as low-frequency any seed that occurs in 8.5% or less of known transcripts for the genome in question. However, because the number of mRNAs for a given species and variability among the 3′ UTRs for those species, a cut off between 5% and 15% would generally be appropriate.
To better understand off-target signatures, a panel of 29 functional siRNAs (each providing >80% gene knockdown) targeting GAPDH or PPIB were individually transfected into HeLa cells (100 nM, 10K cells/well, Lipofectamine 2000). Included in this set were two siRNAs, GAPDH H15 and PPIB H17 that targeted different genes but had the same antisense seed region. Total RNA was collected 24 hrs later and subsequently analyzed (Agilent A1 human microarray using mock-transfected cells as a control reference) to determine whether siRNAs with similar seed regions generated similar off-target profiles.
A heatmap of the results of these experiments is provided in
Previously identified siRNA-off target pairs were used to investigate whether control siRNA (i.e. siRNA that had identical seed regions, but distinct, neutral scaffolds) could be used to confirm false positive phenotypes generated by test siRNA. Work by Lim et al. (NAR 33, 4527-4535, 2005) demonstrated that two unique siRNAs targeting GRK4 and BTK (respectively) down-regulated a reporter construct containing a HIF1 alpha 3′ UTR. As each of the two targeting siRNAs had the same seed region (see sequences below) and the HIF1 alpha 3′ UTR contained two exact seed complements (see bold, underlined sequence below), these results represent a classic example of a false positive phenotype induced by off-target effects.
To test the ability of control siRNA to mimic the false positive effect induced by the GRK4- and BTK-targeting siRNAs, the seed sequence of the targeting siRNAs was embedded into a neutral scaffold (see sequences below) and transfected into HeLa cells (100 nM, DharmaFECT1). Subsequently, the relative levels of HIF1 alpha mRNA were assessed by branched DNA assay to determine whether the control siRNA could mimic the false positive effects induced by the GRK4- and BTK-targeting duplexes. As shown in
A given (or candidate) siRNA may be identified that is thought to cause a particular phenotype such as cell death or a particular level of silencing. Researchers may wish to determine if the hit is due to knockdown of the gene that was being targeted, or if it was the result of an off-target effect by the siRNA.
An siRNA (also referred herein as a control siRNA or a seed control siRNA) that has the same seed as the candidate siRNA that induces the phenotype identified in the previous paragraph is selected from a seed control library. The region of the control siRNA that is not part of the seed region contains a neutral scaffold sequence that has less than 80% sequence similarity with the nucleotides of the candidate siRNA that induces the phenotype. If the original phenotype was the result of an off-target effect, then transfection of this seed control siRNA should induce an identical or similar phenotype as the candidate siRNA as defined by the thresholds of the assay.
In contrast, if the original effect was the result of the target specific knockdown, then this seed control siRNA should not induce the phenotype. The scaffolding may be selected to have no effect when a seed region other than that of the candidate siRNA is employed.
A portion of the highly functional siRNA targeting GAPDH (GAPDH duplex 4, GAPDH4 or G4 OT) was chosen as a scaffolding sequence because the duplex efficiently targets GAPDH but off-targets minimal numbers of genes otherwise. Duplexes representing 15 seeds were synthesized as chimeras in the context of the scaffold sequence of GAPDH4. The sense strand sequences are shown below with the inserted seed reverse complement sequence in bold; all duplexes were synthesized with chemical modification (modification of sense strand nucleotides 1 and 2 (counting from the 5′ end of the oligonucleotide) with 2′-O-methyl modifications, the 5′-most nucleotide of the antisense strand is phosphorylated) to ensure preferential entry of the antisense strand into RISC. In the sequences listed below, “L” represents control siRNA sequences that have low seed complement frequencies, “M” represents control siRNA sequences that have moderate seed complement frequencies, and “H” represents control siRNA sequences that have high seed complement frequencies.
Standard microarray off-targeting analysis demonstrated several points including: (1) that while none of these chimeric molecules could still target GAPDH, they all presented unique microarray signatures; and (2) that chimeric sequences that had seeds with low seed complement frequencies induced (overall) fewer off-target genes than those with moderate or high seed complement frequencies. No common genes were off-targeted among all 16 duplexes, indicating that this scaffold sequence contributes little to nothing to the identity of the off-targeted genes.
A seed control library of molecules can be constructed by synthesizing a set of 19-mer control siRNA with an overhang of 1-6 nucleotides (for example, with UU overhangs on the 3′ end of each strand). Each of the control siRNAs contains one of the possible 4,096 hexamers at the seed position (nucleotides 2-7 on the antisense strand). The reverse complement of each of these seeds is present at positions 13-18 of the sense strand. The duplexes may be synthesized with the chemical modification pattern described in the previous example so as to maximize the introduction of the antisense strand into RISC and to minimize the ability of the sense strand to generate off-target effects. (See US-2005-0223427A1, the contents of which are incorporated by reference.)
The sequence of the duplex that is not defined by the seed region (the scaffold-nucleotides 1-12 and 19 of the sense strand and its reverse complement on the antisense strand) should be selected so as not to interfere with seed-based targeting of this sequence, as well as not having any other undesired effects. Thus, the scaffold region should not contain stretches of homopolymer longer than three bases that could form unusual structures or sequences that could form a fold-back duplex (or hairpin) of that strand alone.
In addition, position 19 of the sense region is preferably an “A” (“U” at position 1 of the antisense region) to possibly allow some unwinding flexibility and to match many known, naturally occurring miRNA sequences. The entire 19-mer sense strand should be determined by BLAST or another identity algorithm to not have a 17-19 base identity with any human gene transcript, which would cause the control duplex to target another message for specific endonucleolytic cleavage by RISC in addition to the seed-based off-targeting mechanism
Examples of possible sense region sequences of scaffolds are provided in SEQ. ID NOs. 13-15. The antisense region may for example, be 100% complementary to the sense regions.
It should be noted that one may choose not to synthesize all 4096 different duplexes (i.e., control siRNAs) for a given scaffolding. One may first test an siRNA designed rationally to be highly functional. Next, one may examine the seed regions for these siRNAs to determine if they exhibit certain phenotypes. Next control siRNAs could be created that contain the seed sequences that correspond only to the seed sequences of those siRNA that show discernible phenotypes.
*Predicted target number based on overall percentage identity
This application is a continuation-in-part of U.S. application Ser. No. 11/724,346, filed Mar. 15, 2007, which claims the benefit of U.S. Provisional Application Ser. No. 60/782,970, filed Mar. 16, 2006. The entire disclosures of those applications are incorporated by reference as if set forth fully herein.
Number | Date | Country | |
---|---|---|---|
60782970 | Mar 2006 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11724346 | Mar 2007 | US |
Child | 11825461 | Jul 2007 | US |