Patent Application
20100035966
This invention generally relates to methods for identifying a subset of microRNA responsive genes that share at least one phenotype, and to methods and compositions for regulating cell cycle progression.
Many genes are related via common regulation, common functional molecular mechanisms, and common pathways. Understanding the relationship between genes is important for biological research and has extensive practical application in drug development and diagnostics.
MicroRNAs are a recently identified class of regulatory RNAs that target specific mRNAs for degradation or inhibition of translation, resulting in a decrease of the protein encoded by the target mRNA. Current estimates are that 30% or more of human mRNAs are regulated by miRNAs (Lewis et al., Cell 120:15-20 (2005). Studies investigating expression profiles of various miRNAs in normal and cancer cells reveals that miRNA expression patterns may have clinical relevance. (See, e.g., Yanaihara, N. et al., Cancer Cell 9:189-198, 2006). Application of various bioinformatics approaches have revealed that a single miRNA might bind to as many as 200 gene targets and these targets are often diverse in function, including, for example, transcription factors, secreted factors, receptors and transporters (see, e.g., Esquela-Kerscher and Slack, Nature Reviews 6:259-269 (2006); Baretl, D. P. et al., Nat Rev Genet 5(5):396-400 (2004)). Therefore, the deletion or overexpression of a particular miRNA is likely to be pleotropic.
To date, over 200 microRNAs have been described in humans, however, the current state of knowledge regarding microRNA targets and the determination of microRNA functions is incomplete. Although thousands of miRNA targets have been predicted using computational methods, relatively few predications have been experimentally validated. Computational methods are not optimal for predicting miRNA target sites. Bioinformatics approaches generally rely heavily on the detection of seed region (the encompasses the first 1-12 bases of the mature miRNA sequence) complementary motifs that are conserved in the 3′ UTR sequences of genes across divergent species (see, e.g., John, B. et al., PloS Biol 2(11):e363, 2004). Therefore, such methods are not predictive for microRNA targets sites that are not conserved across species, or for identifying target sites that are not perfectly matched with seed regions. Moreover, target prediction using different computational methods often do not agree. Since relatively few predicted microRNA: target interactions have been experimentally confirmed, it is difficult to know how accurate such predictions are. Available methods for validation are laborious and not easily amenable to high-throughput methodologies (see e.g., Bentwich, I., FEBS Lett 579:5904-5910 (2005)).
It is important to assign functions to miRNAs and to accurately identify miRNA responsive targets. Since a single miRNA can regulate hundreds of targets, understanding of biological pathways regulated by microRNAs is not obiouvs from examination of their targets. As functions are assigned to miRNAs, it is also important to determine which of their target(s) are responsible for a phenotype. It is also currently unknown whether the numerous miRNA responsive targets act individually or in concert.
In accordance with the foregoing, in one aspect, the present invention provides a method for identifying a subset of genes that are responsive to a selected microRNA species and share at least one phenotype, the method comprising: a) modulating the level of at least one microRNA species in a cell type of interest; b) measuring the level of expression of each member of a plurality of genes in the modulated cells to yield a plurality of expression values; c) measuring the level of expression of each member of the same plurality of genes in a reference cell type that is not modulated to yield a plurality of reference expression values; d) comparing the plurality of expression values from the modulated cells with the plurality of reference expression values to identify a set of genes that are responsive to the microRNA species, each identified gene having an expression value that is statistically different from the corresponding reference value; e) identifying transcripts within the set of genes that contain at least one microRNA responsive target sequence corresponding to the modulated microRNA species; f) modulating the level of expression of a plurality of members of the set of genes that are responsive to the microRNA species and that contain at least one microRNA responsive target sequence in the cell type of interest with a plurality of gene-specific agents; and g) comparing the cells modulated in step (a) and the cells modulated in step (f) for the presence of at least one shared phenotype to identify a subset of genes that are responsive to the selected microRNA species and exhibit a shared phenotype.
In another aspect, the present invention provides a method of inhibiting proliferation of a mammalian cell comprising introducing into said cell an effective amount of at least one at least one small interfering RNA agent (iRNA), wherein said iRNA comprises a nucleotide sequence of at least 15 nucleotides, wherein the nucleotide sequence comprises a seed region consisting of nucleotide positions 1 to 12, wherein position 1 represents the 5′ end of the iRNA nucleotide sequence and wherein said seed region comprises a nucleotide seqeuence of at least six contiguous nucleotides that is complementary to six contiguous nucleotides within positions 1 to 12 of a nucleotide sequence, wherein position 1 represents the 5′ end of the nucleotide sequence, wherein the nucleotide sequence is selected from the group consisting of SEQ ID NO: 1, SEQ ID NO:2, SEQ ID NO:3, SEQ ID NO:4, SEQ ID NO:5, SEQ ID NO:6, SEQ ID NO:7 and SEQ ID NO:8. In some embodiments, the method comprises introducing at least one iRNA that inhibits the expression of at least one miR-16 responsive gene selected from TABLE 5 into the mammalian cell.
In some embodiments, the at least one miR-16 response gene is selected from the group consisting of ATG9A, ATXN7L3, C10orf46, IPPK, C9orf42, C9orf91, CARD10, CBX6, CDC27, CDK6, COX10, H2AFX, KIAA0317, MFN2, PHF17, PPP1R11, RAB11FIP2, and SRPR.
In another aspect, the present invention provides a composition comprising a combination of gene-specific inhibitors directed to at least two miR-16 responsive target genes selected from TABLE 5. In some embodiments, the composition comprises a combination of selective inhibitors directed to at least two miR-16 responsive target genes selected from the group consisting of: ATG9A, ATXN7L3, C10orf46, IPPK, C9orf42, C9orf91, CARD10, CBX6, CDC27, CDK6, COX10, H2AFX, KIAA0317, MFN2, PHF17, PPP1R11, RAB11FIP2, and SRPR. In further embodiments, the composition comprises a combination of selective inhibitors directed to at least two miR-16 responsive target genes selected from the group consisting of CARD10, CDC27, CDK6, and C10orf46.
In another aspect, the present invention provides a composition comprising a at least one gene-specific inhibitor directed to at least one miR-16 responsive target gene selected from TABLE 6.
In yet another aspect, the present invention provides an isolated dsRNA molecule comprising one nucleotide strand that is substantially identical to a sequence selected from the group consisting of SEQ ID NO:236 to SEQ ID NO:361.
The methods of the invention are useful, for example, for identifying a subset of microRNA responsive genes that are associated with a common biological pathway. The compositions comprising gene-specific agents directed to the identified microRNA responsive genes are useful, for example, for modulating one or more biological pathways.
The foregoing aspects and many of the attendant advantages of this invention will become more readily appreciated as the same become better understood by reference to the following detailed description, when taken in conjunction with the accompanying drawings, wherein:
Unless specifically defined herein, all terms used herein have the same meaning as they would to one skilled in the art of the present invention. Practitioners are particularly directed to Sambrook et al. (1989), Molecular Cloning: A Laboratory Manual, 2nd ed., Cold Spring Harbor Press, Plainsview, N.Y. (1989), and Ausubel et al., Current Protocols in Molecular Biology (Supplement 47), John Wiley & Sons, New York (1999), for definitions and terms of the art.
As used herein, the term “microRNA species,” or (“microRNA”) or (“miRNA”) or “mi-R” refers to small, non-protein coding RNA molecules that are expressed in a diverse array of eukaryotes, including mammals. MicroRNA molecules typically have a length in the range of from 15 nucleotides to 120 nucleotides, the size depending upon the specific microRNA species and the degree of intracellular processing. Mature, fully processed miRNAs are about 15 to 30, 15-25, or 20 to 30 nucleotides in length, and more often between about 16 to 24, 17 to 23, 18 to 22, 19 to 21 or 21 to 24 nucleotides in length. MicroRNAs include processed sequences as well as corresponding long primary transcripts (pri-miRNAs) and processed precursors (pre-miRNAs). Some microRNA molecules function in living cells to regulate gene expression via RNA interference. A representative set of microRNA species is described in the publicly available miRBase sequence database as described in Griffith-Jones et al., Nucleic Acids Research 32:D109-D111 (2004) and Griffith-Jones et al., Nucleic Acids Research 34:D140-D144 (2006), accessible on the World Wide Web at the Wellcome Trust Sanger Institute website.
As used herein, the term “microRNA family” refers to a group of microRNA species that share identity across at least 6 consecutive nucleotides within nucleotide positions 1 to 12 of the 5′ end of the microRNA molecule, also referred to as the “seed region,” as described in Brennecke, J. et al., PloS biol 3(3):pe85 (2005).
As used herein, the term “microRNA family member” refers to a microRNA species that is a member of a microRNA family, including naturally occurring microRNA species and articifically generated microRNA molecules.
As used herein, the term “RNA interference” or “RNAi” refers to the silencing or decreasing of gene expression by iRNA agents (e.g., siRNAs, miRNAs, shRNAs), via the process of sequence-specific, post-transcriptional gene silencing in animals and plants, initiated by an iRNA agent that has a seed region sequence in the iRNA guide strand that is complementary to a sequence of the silenced gene.
As used herein, the term an “iRNA agent” (abbreviation for “interfering RNA agent”), refers to an RNA agent, or chemically modified RNA, which can down-regulate the expression of a target gene. While not wishing to be bound by theory, an iRNA agent may act by one or more of a number of mechanisms, including post-transcriptional cleavage of a target mRNA, or pre-transcriptional or pre-translational mechanisms. An iRNA agent can include a single strand (ss) or can include more than one strands, e.g. it can be a double stranded (ds) IRNA agent.
As used herein, the term “single strand iRNA agent” is an iRNA agent which consists of a single molecule. It may include a duplexed region, formed by intra-strand pairing, e.g., it may be, or include, a hairpin or panhandle structure. The ssRNA agents of the present invention include transcripts that adopt stem-loop structures, such as shRNA, that are processed into siRNA.
As used herein, the term “ds iRNA agent” is a dsRNA (double stranded RNA) agent that includes two strands that are not covalently linked, in which interchain hybridization can form a region of duplex structure. The dsRNA agents of the present invention include silencing dsRNA molecules that are sufficiently short that they do not trigger the interferon response in mammalian cells.
As used herein, the term “siRNA” refers to a small interfering RNA. siRNAs include short interfering RNA of about 15-60, 15-50, or 15-40 (duplex) nucleotides in length, more typically about 15-30, 15-25 or 19-25 (duplex) nucleotides in length, and is preferably about 20-24 or about 21-22 or 21-23 (duplex) nucleotides in length (e.g., each complementary sequence of the double stranded siRNA is 15-60, 15-50, 15-40, 15-30, 15-25 or 19-25 nucleotides in length, preferably about 20-24 or about 21-22, or 21-23 nucleotides in length, preferably 19-21 nucleotides in length, and the double stranded siRNA is about 15-60, 15-50, 15-40, 15-30, 15-25 or 19-25, preferably about 20-24, or about 21-22 or 19-21 or 21-23 base pairs in length). siRNA duplexes may comprise 3′ overhanges of about 1 to about 4 nucleotides, preferably about 2 to 3 nucleotides and 5′ phosphate termini. In some embodiments, the siRNA lacks a terminal phosphate. In some embodiments, one or both ends of siRNAs can include single-stranded 3′ overhangs that are two or three nucleotides in length, such as, for example, deoxythymidine (dTdT) or uracil (UU) that are not complementary to the target sequence. In some embodiments, siRNA molecules can include nucleotide analogs (e.g., thiophosphate or G-clamp nucleotide analogs), alternative base linkages (e.g., phosphorothioate, phosphonoacetate, or thiophosphonoacetate) and other modifications useful for enhanced nuclease resistance, enhanced duplex stability, enhanced cellular uptake, or cell targeting.
In certain embodiments, at least one of the two strands of the siRNA further comprises a 1-4, preferably a 2 nucleotide 3′ overhang. The nucleotide overhang can include any combination of a thymine, uracil, adenine, guanine, or cytosine, or derivatives or analogues thereof. The nucleotide overhang in certain aspects is a 2 nucleotide overhang, where both nucleotides are thymine. Importantly, when the dsRNA comprising the sense and antisense strands is administered, it directs target specific interference and bypasses an interferon response pathway. In order to enhance the stability of the short interfering nucleic acids, the 3′ overhangs can also be stabilized against degradation. In one embodiment, the 3′ overhangs are stabilized by including purine nucleotides, such as adenosine or guanosine nucleotides. Alternatively, substitution of pyrimidine nucleotides by modified analogues, e.g., substitution of uridine nucleotides in the 3′ overhangs with 2′-deoxythymidine, is tolerated and does not affect the efficiency of RNAi degradation. In particular, the absence of a 2′ hydroxyl in the 2′-deoxythymidine significantly enhances the nuclease resistance of the 3′ overhang in tissue culture medium.
As used herein, a “3′ overhang” refers to at least one unpaired nucleotide extending from the 3′ end of an siRNA sequence. The 3′ overhang can include ribonucleotides or deoxyribonucleotides or modified ribonucleotides or modified deoxyribonucleotides. The 3′ overhang is preferably from 1 to about 5 nucleotides in length, more preferably from 1 to about 4 nucleotides in length and most preferably from about 2 to about 4 nucleotides in length. The 3′ overhang can occur on the sense or antisense sequence, or on both sequences of an RNAi construct. The length of the overhangs can be the same or different for each strand of the duplex. Most preferably, a 3′ overhang is present on both strands of the duplex, and the overhang for each strand is 2 nucleotides in length. For example, each strand of the duplex can comprise 3′ overhangs of dithymidylic acid (“tt”) or diuridylic acid (“uu”).
As used herein, the siRNA molecules need not be limited to those molecules containing only RNA, but may further encompass chemically-modified nucleotides and non-nucleotides. WO2005/078097; WO2005/0020521 and WO2003/070918 detail various chemical modifications to RNAi molecules, and the contents of each reference is hereby incorporated by reference in its entirety. In certain embodiments, for example, the siRNA molecules may lack 2′-hydroxyl (2′-OH) containing nucleotides. The siRNA can be chemically synthesized or may be encoded by a plasmid (e.g., transcribed as sequences that automatically fold into duplexes with hairpin loops). siRNA can also be generated by cleavage of longer dsRNA (e.g., dsRNA greater than about 25 nucleotides in length) with the E. Coli RNAse III or Dicer. These enzymes process the dsRNA into biologically active siRNA (see, e.g., Yang et al., 2002 PNAS USA 99:9942-7; Calegari et al., 2002, PNAS USA 99:14236; Byrom et al., 2003, Ambion TechNotes 10(1): 4-6; Kawasaki et al., 2003, Nucleic Acids Res. 31: 981-7; Knight and Bass, 2001, Science 293: 2269-71; and Robertson et al., 1968, J. Biol. Chem. 243: 82). The long dsRNA can encode for an entire gene transcript or a partial gene transcript.
Non limiting examples of siRNA molecules of the invention may include a double-stranded polynucleotide molecule comprising self-complementary sense and antisense regions, wherein the antisense region comprises nucleotide sequence that is complementary to nucleotide sequence in a target nucleic acid molecule or a portion thereof (alternatively referred to as the guide region, or guide strand when the molecule contains two separate strands) and the sense region having nucleotide sequence corresponding to the target nucleic acid sequence or a portion thereof (also referred as the passenger region, or the passenger strand when the molecule contains two separate strands). The siRNA can be assembled from two separate oligonucleotides, where one strand is the sense strand and the other is the antisense strand, wherein the antisense and sense strands are self-complementary (i.e., each strand comprises nucleotide sequence that is complementary to nucleotide sequence in the other strand; such as where the antisense strand and sense strand form a duplex or double stranded structure, for example wherein the double stranded region is about 18 to about 30, e.g., about 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29 or 30 base pairs); the antisense strand (guide strand) comprises nucleotide sequence that is complementary to nucleotide sequence in a target nucleic acid molecule or a portion thereof and the sense strand (passenger strand) comprises nucleotide sequence corresponding to the target nucleic acid sequence or a portion thereof (e.g., about 15 to about 25 nucleotides of the siRNA molecule are complementary to the target nucleic acid or a portion thereof). Typically, a short interfering RNA (siRNA) refers to a double-stranded RNA molecule of about 17 to about 29 base pairs in length, preferably from 19-21 base pairs, one strand of which is complementary to a target mRNA, that when added to a cell having the target mRNA or produced in the cell in vivo, causes degradation of the target mRNA. Preferably the siRNA is perfectly complementary to the target mRNA. But it may have one or two mismatched base pairs.
Alternatively, the siRNA is assembled from a single oligonucleotide, where the self-complementary sense and antisense regions of the siRNA are linked by means of a nucleic acid based or non-nucleic acid-based linker(s). The siRNA can be a polynucleotide with a duplex, asymmetric duplex, hairpin or asymmetric hairpin secondary structure, having self-complementary sense and antisense regions, wherein the antisense region comprises nucleotide sequence that is complementary to nucleotide sequence in a separate target nucleic acid molecule or a portion thereof and the sense region having nucleotide sequence corresponding to the target nucleic acid sequence or a portion thereof. The siRNA can be a circular single-stranded polynucleotide having two or more loop structures and a stem comprising self-complementary sense and antisense regions, wherein the antisense region comprises nucleotide sequence that is complementary to nucleotide sequence in a target nucleic acid molecule or a portion thereof and the sense region having nucleotide sequence corresponding to the target nucleic acid sequence or a portion thereof, and wherein the circular polynucleotide can be processed either in vivo or in vitro to generate an active siRNA molecule capable of mediating RNAi. The siRNA can also comprise a single stranded polynucleotide having nucleotide sequence complementary to nucleotide sequence in a target nucleic acid molecule or a portion thereof (for example, where such siRNA molecule does not require the presence within the siRNA molecule of nucleotide sequence corresponding to the target nucleic acid sequence or a portion thereof), wherein the single stranded polynucleotide can further comprise a terminal phosphate group, such as a 5′-phosphate (see for example Martinez et al., 2002, Cell 110:563-574 and Schwarz et al., 2002, Molecular Cell, 10:537-568), or 5′,3′-diphosphate. In certain embodiments, the siRNA molecule of the invention comprises separate sense and antisense sequences or regions, wherein the sense and antisense regions are covalently linked by nucleotide or non-nucleotide linkers molecules as is known in the art, or are alternately non-covalently linked by ionic interactions, hydrogen bonding, van der waals interactions, hydrophobic interactions, and/or stacking interactions. In certain embodiments, the siRNA molecules of the invention comprise nucleotide sequence that is complementary to nucleotide sequence of a target gene. In another embodiment, the siRNA molecule of the invention interacts with nucleotide sequence of a target gene in a manner that causes inhibition of expression of the target gene.
As used herein, “percent modification” refers to the number of nucleotides in the iRNA, or each of the strand of the siRNA or to the collective dsRNA that have been modified. Thus 19% modification of the antisense strand refers to the modification of up to 4 nucleotides/bp in a 21 nucleotide sequence (21 mer). 100% refers to a fully modified dsRNA. The extent of chemical modification will depend upon various factors well known to one skilled in the art. Such, as for example, target mRNA, off-target silencing, degree of endonuclease degradation, etc.
As used herein, the term “shRNA” or “short hairpin RNAs” refers to individual transcripts that adopt stem-loop structures which are processed into siRNA by RNAi machinery. Typical shRNA molecules comprise two inverted repeats containing the sense and antisense target sequence separated by a loop sequence. The base-paired segment may vary from 17 to 29 nucleotides, wherein one strand of the base-paied stem is complementary to the mRNA of a target gene. The loop of the shRNA stem-loop structure may be any suitable length that allows inactivation of the target gene in vivo. While the loop may be from 3 to 30 nucleotides in length, typically it is 1-10 nucleotides in length. The base paired stem may be perfectly base paired or may have 1 or 2 mismatched base pairs. The duplex portion may, but typically does not, contain one or more bulges consisting of one or more unpaired nucleotides. The shRNA may have non-base-paired 5′ and 3′ sequences extending from the base-paired stem. Typically, however, there is no 5′ extension. The first nucleotide of the shRNA at the 5′ end is a G, because this is the first nucleotide transcribed by polymerase III. If G is not present as the first base in the target sequence, a G may be added before the specific target sequence. The 5′ G typically forms a portion of the base-paired stem. Typically, the 3′ end of the shRNA is a poly U segment that is a transcription termination signal and does not form a base-paired structure. As described in the application and known to one skilled in the art, shRNAs are processed into siRNAs by the conserved cellular RNAi machinery. Thus shRNAs are precursors of siRNAs and are, in general, similarly capable of inhibiting expression of a target mRNA transcript. For the purpose of description, in certain embodiments, the shRNA constructs of the invention target one or more mRNAs that are targeted by a member of the miR-16 family including miR-16, miR-15a, miR-15b, miR103, miR-107, miR-195 and kshv-miR-K12-6-5p. The strand of the shRNA that is antisense to the target gene transcript is also known as the “guide strand”.
As used herein, the term “microRNA responsive target site” refers to a nucleic acid sequence ranging in size from about 5 to about 25 nucleotides (such as 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, or 25 nucleotides) that is complementary (i.e., a direct match), or essentially complementary to at least a portion of a microRNA molecule. In some embodiments, the microRNA responsive target site comprises at least 6 consecutive nucleotides, at least 7 consecutive nucleotides, at least 8 consecutive nucleotides, or at least 9 consecutive nucleotides that are complementary to the seed region of a microRNA molecule (i.e., within positions 1 to 12 of the 5′ end of the microRNA molecule, referred to as the “seed region.” See, e.g., Brennecke, J. et al., PloS biol 3(3):pe85 (2005)).
As used herein, the term “isolated” in the context of an isolated nucleic acid molecule, is one which is altered or removed from the natural state through human intervention. For example, an RNA naturally present in a living animal is not “isolated.” A synthetic RNA or dsRNA or microRNA molecule partially or completely separated from the coexisting materials of its natural state, is “isolated.” Thus, an miRNA molecule which is deliberately delivered to or expressed in a cell is considered an “isolated” nucleic acid molecule.
As used herein, the term “complementary” refers to nucleic acid sequences that are capable of base-pairing according to the standard Watson-Crick complementary rules. That is, the larger purines will base pair with the smaller pyrimidines to form combinations of guanine paired with cytosine (G:C) and adenine paired with either thymine (A:T) in the case of DNA, or adenine paired with uracil (A:U) in the case of RNA.
As used herein, the term “essentially complementary” refers to microRNA target nucleic acid sequences that are longer than 5 nucleotides that are complementary (an exact match) to at least 5 consecutive nucleotides of the 5′ portion of a microRNA molecule from nucleotide positions 1 to 12, (also referred to as the “seed region”), and are at least 85% complementary (such as at least 88%, at least 90%, at least 95%, or at least 96% identical) across the remainder of the microRNA target nucleic acid sequence as compared to a naturally occurring microRNA species set forth in the miRBase sequence database, which is publically accessible on the World Wide Web at the Wellcome Trust Sanger Institute at http://microna.sangerac.uklsequences/, and as described in the following references: Ambros et al., RNA 9:277-279 (2003); Griffith-Jones, Nucleic Acids Res. 32:D109-D111 (2004); Griffith-Jones, Nucleic Acids Res. 34:D140-D144 (2006); Lagos-Quintana et al., Curr. Biol. 12(9):735-9 (2002); Lim LP et al., Science 299(5612):1540 (2003).
As used herein, the term “equivalent” with reference to a microRNA target sequence, is a nucleic acid sequence that is identical to at least 5 consecutive nucleotides of the 5′ portion of a particular microRNA target sequence from nucleotide positions 1 to 12 (also referred to as the “seed region”), and is at least 85% identical (such as at least 88%, at least 90%, at least 95%, or at least 96% identical) to the remainder of the microRNA target sequence. The comparison of sequences and determination of percent identity and similarity between two sequences can be accomplished using a mathematical algorithm of Karlin and Altschul (Proc. Natl. Acad. Sci. USA 87:2264-2268 (1990)), modified as in Karlin and Atlschul (Proc. Natl. Acad. Sci. USA 90:5873-5877 (1993)). Such an algorithm is incorporated into the NBLAST and XBLAST programs of Altshcul et al. (J. Mol. Biol. 215:403-410 (1990)).
As used herein, the term “gene” refers to a nucleic acid (e.g., DNA or RNA) sequence that comprises coding sequences necessary for the production of an RNA and/or a polypeptide, or its precursor as well as noncoding sequences (untranslated regions) surrounding the 5′ and 3′ ends of the coding sequences. The term “gene” encompasses both cDNA and genomic forms of a gene. A functional polypeptide can be encoded by a full length coding sequence or by any portion of the coding sequence as long as the desired activity or functional properties (e.g., enzymatic activity, ligand binding, signal transduction, antigenic presentation) of the polypeptide are retained. The sequences which are located 5′ of the coding region and which are present on the mRNA are referred to as 5′ untranslated sequences (“5′UTR”). The sequences which are located 3′ or downstream of the coding region and which are present on the mRNA are referred to as 3′ untranslated sequences, or (“3′UTR”).
As used herein, the term “expression cassette” refers to a nucleic acid molecule which comprises at least one nucleic acid sequence that is to be expressed, along with its transcription and translational control sequences. The expression cassette typically includes restriction sites engineered to be present at the 5′ and 3′ ends such that the cassette can be easily inserted, removed, or replaced in a gene delivery vector. Changing the cassette will cause the gene delivery vector into which it is incorporated to direct the expression of a different sequence.
As used herein, the term “phenotype” encompasses the meaning known to one of skill in the art, including modulation of the expression of one or more genes, as measured by gene expression analysis or protein expression analysis.
As used herein, the term “source of biological knowledge” refers to information that describes the function (e.g., at molecular, cellular and system levels), structure, pathological roles, toxicological implications, etc., of a multiplicity of genes. Various sources of biological knowledge can be used for the methods of the invention, including databases and information collected from public sources such as Locuslink, Unigene, SwissTrEMBL, etc., and organized into a relational database following the concept of the central dogma of molecular biology. In some embodiments, the annotation systems used by the Gene Ontology (GO) Consortium or similar systems are employed. GO is a dynamic controlled vocabulary for molecular biology which can be applied to all organisms as knowledge of gene is accumulating and changing, it is developed and maintained by Gene Ontology™ Consortium (Gene Ontology: tool for the unification of biology. The Gene Ontology Consortium (2000), Nature Genet. 25:25-29)).
As used herein, the phrase to “inhibit the proliferation of a mammalian cell” means to kill the cell, or permanently or temporarily arrest the growth of the cell. Inhibition of a mammalian cell can be inferred if the number of such cells, either in an in vitro culture vessel, or in a subject, remains constant or decreases after administration of the compositions of the invention. An inhibition of tumor cell proliferation can also be inferred if the absolute number of such cells increases, but the rate of tumor growth decreases.
The phrase “inhibiting expression of a target gene” refers to the ability of an RNAi agent, such as an siRNA, to silence, reduce, or inhibit expression of a target gene. Said another way, to “inhibit”, “down-regulate”, or “reduce”, it is meant that the expression of the gene, or level of RNA molecules or equivalent RNA molecules encoding one or more proteins or protein subunits, or activity of one or more proteins or protein subunits, is reduced below that observed in the absence of the RNAi agent. For example, an embodiment of the invention proposes inhibiting, down-regulating or reducing expression of one or more miR-16 responsive genes, by introduction of an miR-16-like siRNA molecule, below the level observed for that miR-16 responsive gene in a control cell to which an miR-16-like siRNA molecule has not been introduced. In another embodiment, inhibition, down-regulation, or reduction contemplates inhibition of the target mRNA below the level observed in the presence of, for example, an siRNA molecule with scrambled sequence or with mismatches. In yet another embodiment, inhibition, down-regulation, or reduction of gene expression with a siRNA molecule of the instant invention is greater in the presence of the invention siRNA e.g., siRNA that down regulates one or more miR-16 responsive gene mRNAs levels than in its absence. In one embodiment, inhibition, down regulation, or reduction of gene expression is associated with post transcriptional silencing, such as RNAi mediated cleavage of a target nucleic acid molecule (e.g. RNA) or inhibition of translation.
To examine the extent of gene silencing, a test sample (e.g., a biological sample from organism of interest expressing the target gene(s) or a sample of cells in culture expressing the target gene(s)) is contacted with an siRNA that silences, reduces, or inhibits expression of the target gene(s). Expression of the target gene in the test sample is compared to expression of the target gene in a control sample (e.g., a biological sample from organism of interest expressing the target gene or a sample of cells in culture expressing the target gene) that is not contacted with the siRNA. Control samples (i.e., samples expressing the target gene) are assigned a value of 100%. Silencing, inhibition, or reduction of expression of a target gene is achieved when the value of test the test sample relative to the control sample is about 95%, 90%, 85%, 80%, 75%, 70%, 65%, 60%, 55%, 50%, 45%, 40%, 35%, 30%, 25%, 20%, or 10%. Suitable assays include, e.g., examination of protein or mRNA levels using techniques known to those of skill in the art such as dot blots, northern blots, in situ hybridization, ELISA, microarray hybridization, immunoprecipitation, enzyme function, as well as phenotypic assays known to those of skill in the art.
An “effective amount” or “therapeutically effective amount” of an siRNA or an RNAi agent is an amount sufficient to produce the desired effect, e.g., inhibition of expression of a target sequence in comparison to the normal expression level detected in the absence of the siRNA or RNAi agent. Inhibition of expression of a target gene or target sequence by a siRNA or RNAi agent is achieved when the expression level of the target gene mRNA or protein is about 90%, 80%, 70%, 60%, 50%, 40%, 30%, 25%, 20%, 15%, 10%, 5%, or 0% relative to the expression level of the target gene mRNA or protein of a control sample.
In one aspect, the present invention provides a method for identifying a subset of genes that are responsive to a selected microRNA species and share at least one phenotype, the method comprising: a) modulating the level of at least one microRNA species in a cell type of interest; b) measuring the level of expression of each member of a plurality of genes in the modulated cells to yield a plurality of expression values; c) measuring the level of expression of each member of the same plurality of genes in a reference cell type that is not modulated to yield a plurality of reference expression values; d) comparing the plurality of expression values from the modulated cells with the plurality of reference expression values to identify a set of genes that are responsive to the microRNA species, each identified gene having an expression value that is statistically different from the corresponding reference value; e) identifying transcripts within the set of genes that contain at least one microRNA responsive target sequence corresponding to the modulated microRNA species; f) modulating the level of expression of a plurality of members of the set of genes that are responsive to the microRNA species and that contain at least one microRNA responsive target sequence in the cell type of interest with a plurality of gene-specific agents; and g) comparing the cells modulated in step (a) and the cells modulated in step (f) for the presence of at least one shared phenotype to identify a subset of genes that are responsive to the selected microRNA species and exhibit a shared phenotype.
As shown in
In accordance with the methods of this aspect of the invention, the level of at least one microRNA species is modulated (i.e., increased or decreased) in a cell type of interest. In one embodiment, the level of microRNA species is decreased in a cell type of interest. A decrease in microRNA expression may be achieved using any suitable method, such as introducing an inhibitory agent, such as an iRNA agent selected to inhibit expression of the endogenous gene encoding the microRNA.
In another embodiment, the level of microRNA species is increased in a cell type of interest. An increase in expression of a microRNA species may be achieved using any suitable method, such as by inducing expression of the endogenous microRNA species, by introducing an expression vector encoding a microRNA, or by introducing one or more microRNA duplex molecules into the cell type of interest.
The methods of this aspect of the invention can be practiced to modulate the level of any microRNA gene product in a cell type of interest. Cellular endogenous microRNAs (miRNAs) are an evolutionary conserved group of noncoding 17-24 nucleotide long single-stranded RNA molecules that are expressed in plants and animals (see, Bartel, D. P., Cell 116:281-297, 2004). MicroRNA single-strands are often referred to in the art as “guide strands” or “active strands”, while the complement to the active strand is referred to as the “passenger strand.” MicroRNAs have been shown to inhibit the expression of target genes by interacting with complementary sites in the 3′ untranslated region (UTR) of the target mRNAs (see, Olsen et al., Dev. Biol. 216:671-680 (1999); Bartel et al., Cell 116:281-297 (2004)). Although most predicted microRNA target recognition sites lie in 3′ UTR regions, other regions, including coding region recognition is also observed (Brennecke, J. et al., PloS Biol. 3(3):pe85 (2005); Lewis, B. P. et al., Cell 120(1):15-20 (2005)). Genetic and biochemical studies have indicated that microRNAs are processed to their mature forms by RNAse III family nucleases, and function through RNA-mediated interference (RNAi) and related pathways to regulate the expression of target genes (Hannon et al., Nature 418:244-251 (2002); Pasquinelli et al., Annu. Rev. Cell Dev. Biol. 18:495-513 (2002)). These microRNAs function as natural triggers of the RNAi pathway.
Computational and molecular cloning approaches have revealed hundreds of microRNAs that are expressed at various levels in a variety of organisms. Over 200 different mammalian microRNAs have been identified, as described in the “miRBase sequence database” which is publically accessible on the World Wide Web at the Wellcome Trust Sanger Institute website at http://microrna.sanger.ac.uk/sequences/. A list of exemplary microRNA species is also described in the following references: Ambros et al., RNA 9:277-279 (2003); Griffith-Jones, Nucleic Acids Res. 32:D109-D111 (2004); Griffith-Jones, Nucleic Acids Res. 34:D140-D144 (2006); Lagos-Quintana et al., Curr. Biol. 12(9):735-9 (2002); Lim, L. P., et al., Science 299(5612):1540 (2003).
MicroRNA expression studies have revealed tissue-specific expression levels of microRNA species in mouse embryos (Mansfield et al., Nature Genetics 36(10):1079-1083 (2004); Houbaviy et al., Developmental Cell 5:351-358 (2003)); human (Sempere et al., Genome Biology 5:R13 (2004); Thompson et al., Nature Methods 1(1):1-7 (2004); Sun, Y. et al., Nucl. Acid Res. 32(22):e188 (2004); Liu et al., Proc. Natl Acad. Sci. 101(5):9740-9744 (2004); Lu et al., Nature 435:834-838 (2005); Barad et al., Genome Research 14:2487-2494 (2004); Baskerville et al., RNA 11:241-247 (2005)); and zebrafish (Wienholds et al., Science 309:310-311 (2005)). The high conservation of expression of tissue specific and tissue-enriched miRNAs between species suggests that miRNAs may play a conserved role in the establishment and/or maintenance of a cell or tissue type (see, Sempere et al., Genome Biology 5:R13 (2004)).
Families of microRNAs have been identified whose members share a region of 5′ identity but differ in their 3′ ends. It has been shown that two different microRNA family members that shared a common 5′ sequence that was complementary to a single 8-mer seed site in the bagpipe 3′ UTR were capable of repressing expression of a reporter gene containing the 8-mer target, even though the 3′ ends of the microRNAs differed, indicating that the target site was responsive to both microRNAs in this family (Brennecke et al., PloS Biology 3(3):e85 (2005)). For example,
In one embodiment, the level of at least one microRNA species is increased in a cell type of interest by introducing the microRNA species into the cell. The introduced microRNA species may be encoded in an expression vector, or may be a chemically synthesized or recombinantly produced gene product. The microRNA species for use in the practice of the methods of the invention can be obtained using a number of standard techniques. For example, the gene products can be chemically synthesized or recombinantly produced as described in more detail below.
The microRNA species may be introduced into the cell using various methods such as infection with a viral vector encoding the microRNA, microinjection, or by transfection using electroporation or with the use of a transfection agent. Transfection methods for mammalian cells are well known in the art, and include direct injection of the nucleic acid into the nucleus of a cell, electroporation, liposome transfer or transfer mediated by lipophilic materials, receptor mediated nucleic acid delivery, bioballistic or particle acceleration, calcium phosphosphate precipitation and transfection mediated by viral vectors. For example, cells can be transfected with a liposomal transfer compound, e.g., DOTAP (N-[1-(2,3-dioleoyloxy)propyl]-N,N,N,-trimethyl-ammonium methylsulfate, Boehringer-Mannheim) or an equivalent, such as LIPOFECTIN. An exemplary method for transfecting miRNA into mammalian cells is described in EXAMPLE 1.
The methods of this aspect of the invention may be practiced using any cell type, such as primary cells or an established line of cultured cells may be used in the practice of the methods of the invention. For example, the methods may be used in any mammalian cell from a variety of species, such as a cow, horse, mouse, rat, dog, pig, goat, or primate, including a human. In some embodiments, the methods may be used in a mammalian cell type that has been modified, such as a cell type derived from a transgenic animal or a knockout mouse. In another embodiment, the microRNA is modulated in a cell type of interest due to infection with a virus that naturally encodes a microRNA, such as kshv-miR-K12-6-5p, which encodes a member of the miR-16 family (see Cai et al., PNAS 102:5570-5575 (2005)).
In some embodiments, the method of the invention is practiced using a cancer cell type. Representative examples of suitable cancer cell types that can be cultured in vitro and used in the practice of the present invention are colon cancer cells, such as wild type HCT116, wild-type DLD-1, HCT116Dicerex5 and DLD-1 Dicerex5 cells described in Cummins, J. M., et al., PNAS 103(10):3687-3692 (2006)). Other non-limiting examples of suitable cancer cell types include A549, MCF7, and TOV21G and are available from the American Type Culture Collection, Rockville, Md. In further embodiments, the cell type is an microRNA mediated cancer cell type. For example, it has been shown that miR-17, 18, 19, 20, 25, 92, 93 and 106 corresponds to clusters of miRNAs that have been found to be expressed in skeletal muscle and dendritic cells and upregulated by Myc (O'Donnell et al., Nature 435:828 (2005)) and to promote tumor growth in a mouse model of B-cell lymphoma (He et al., Nature 435:828 (2005)). As another example, it has been shown that the locus encoding miR-15a and miR-16 at 13q14 is deleted in more than half of B cell chronic lymphocytic leukemias (CLL) (Calin et al., PNAS 99:15524-15529 (2002)). miR-16 is also subject to mutations in CLL patients, but not control subjects, and its locus is subject to loss of heterozygosity (Calin et al., N Engl J Med 353:1793-1801 (2005)). Taken together, these findings suggest that the miR-15a/miR-16 locus behaves as a classical tumor suppressor locus. As described herein, it has been observed that the miR-16 family of microRNAs negatively regulates cell cycle progression by inducing G0/G1 cell cycle accumulation via coordinate regulation of miR-16 responsive targets. Accordingly, in some embodiments, the methods can be used to identify microRNA responsive genes in miR-15 or miR-16 mediated cancer cells, such as, for example, CLL or prostate cancer cells that may be useful as therapeutic targets, either alone, or in combination.
Identification of MicroRNA-responsive Modulated Transcripts
In the practice of this aspect of the invention, the amount of gene expression (as determined by measuring the amount of mRNA transcribed from a gene, or as represented by the amount of cDNA made from the transcribed mRNA, or as measured by the amount of protein produced that is encoded by the gene) is measured and compared to reference values to yield gene expression patterns that provide information about the effect of a microRNA (or IRNA agent) on a selected cell type. Any method of gene expression analysis may be utilized, such as RNA profiling using Northern blots, quantitative PCR, microarray hybridization analysis, or protein profiling using protein detection methods such as antibody detection methods, (e.g., immunoblots, ELISAs) or mass spectrometry, and the like.
In one embodiment of the method, a set of microRNA-responsive modulated transcripts are identified by measuring the gene expression profile in the microRNA transfected cells in comparison to the gene expression profile in a reference sample, such as mock-transfected cells, to yield an expression value for each transcript within the population of measured transcripts, and then performing at least one calculation on all of the expression values to yield an expression value that numerically represents the expression pattern of the population of genes in response to the microRNA. In some embodiments of the invention, two separate cell lines are transfected in parallel with the same microRNA(s), as described in EXAMPLE 1 and an intersection value is then obtained to provide a more accurate representation of gene expression.
In one embodiment, RNA transcript levels are measured to assess gene expression. Any suitable RNA preparation can be utilized, such as total cellular RNA, or such as cytoplasmic RNA or such as an RNA preparation that is enriched for messenger RNA (mRNA), such as RNA preparations that include greater than 70%, or greater than 80%, or greater than 90%, or greater than 95%, or greater than 99% messenger RNA. Typically, RNA preparations that are enriched for messenger RNA are utilized to provide the RNA template in the practice of the methods of this aspect of the invention. Messenger RNA can be purified in accordance with any art-recognized method, such as by the use of oligo-dT columns (see, e.g., Sambrook et al. 1989, Molecular Cloning—A Laboratory Manual (2nd Ed.), Vol. 1, Chapter 7, Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y.). In a specific embodiment, the mRNA molecules in the RNA sample comprise at least 100 different nucleotide sequences. In other embodiments, the mRNA molecules of the RNA sample comprise at least 500, 1,000, 5,000, 10,000, 20,000, 30,000, 40,000, 50,000, 60,000, 70,000, 80,000, 90,000 or 100,000 different nucleotide sequences. In another specific embodiment, the RNA sample is a mammalian RNA sample, the mRNA molecules of the mammalian RNA sample comprising about 20,000 to 30,000 different nucleotide sequences, or comprising substantially all of the different mRNA sequences that are expressed in the cell(s) from which the mRNA was extracted. In some embodiments, the RNA is isolated from the microRNA transfected cells during a time period from about 4 hours to about 24 hours. In some embodiments, the RNA is isolated within 6 hours after transfection.
In one exemplary embodiment of the method, the gene expression profile of the microRNA modulated cells is measured by microarray hybridization. In accordance with this embodiment, after transfection of microRNA molecules into a desired cell type messenger RNA is extracted (and may or may not be purified) from the transfected cells and may be used as a template to synthesize cDNA or cRNA which is then labeled (e.g., with a fluorescent dye). In some embodiments of the method the labeled cDNA or cRNA is then hybridized to nucleic acid molecules immobilized on a substrate (e.g., a DNA microarray). The immobilized nucleic acid molecules represent some, or all, of the genes that are expressed in the cells that were transfected with the at least one microRNA. The labeled cDNA or cRNA molecules that hybridize to the nucleic acid molecules immobilized on the DNA array are identified, and the level of expression of each hybridizing cDNA or cRNA is measured and compared to the level of expression of the same cDNA or cRNA species in control cells that were not transfected, thereby revealing a gene expression pattern that was caused by the microRNA.
In some embodiments of the method, cDNA molecules are synthesized that are complementary to the RNA template molecules. Each cDNA molecule is preferably sufficiently long (e.g., at least 50 nucleotides in length) to subsequently serve as a specific probe for the mRNA template from which it was synthesized, or to serve as a specific probe for a DNA sequence that is identical to the sequence of the mRNA template from which the cDNA molecule was synthesized. Individual DNA molecules can be complementary to a whole RNA template molecule, or to a portion thereof. Thus, a population of cDNA molecules is synthesized that includes individual DNA molecules that are each complementary to all, or to a portion, of a template RNA molecule. Typically, at least a portion of the complementary sequence of at least 95% (more typically at least 99%) of the template RNA molecules are represented in the population of cDNA molecules.
Any reverse transcriptase molecule can be utilized to synthesize the cDNA molecules, such as reverse transcriptase molecules derived from Moloney murine leukemia virus (MMLV-RT), avian myeloblastosis virus (AMV-RT), bovine leukemia virus (BLV-RT), Rous sarcoma virus (RSV) and human immunodeficiency virus (HIV-RT). A reverse transcriptase lacking RNaseH activity (e.g., S
The synthesis of the cDNA molecules can be primed using any suitable primer, typically an oligonucleotide in the range of ten to 60 bases in length. Oligonucleotides that are useful for priming the synthesis of the cDNA molecules can hybridize to any portion of the RNA template molecules, including the oligo-dT tail. In some embodiments, the synthesis of the cDNA molecules is primed using a mixture of primers, such as a mixture of primers having random nucleotide sequences. Typically, for oligonucleotide molecules less than 100 bases in length, hybridization conditions are 5° C. to 10° C. below the homoduplex melting temperature (Tm); (see generally, Sambrook et al. Molecular Cloning: A Laboratoy Manual, 2nd ed., Cold Spring Harbor Press, 1987; Ausubel et al., Current Protocols in Molecular Biology, Greene Publishing, 1987).
The cDNA molecules are typically labeled to facilitate the detection of the cDNA molecules when they are used as a probe in a hybridization experiment, such as a probe used to screen a DNA microarray, to identify an efficacy-related population of genes. The cDNA molecules can be labeled with any useful label, such as a radioactive atom (e.g., 32P), but typically the cDNA molecules are labeled with a dye. Examples of suitable dyes include fluorophores and chemiluminescers.
In the context of the present example, the labeled cDNA is hybridized to a DNA array that includes hundreds, or thousands, of identified nucleic acid molecules (e.g., cDNA molecules) that correspond to genes that are expressed in the type of cells wherein gene expression is being analyzed. Typically, hybridization conditions used to hybridize the labeled cDNA to a DNA array are no more than 25° C. to 30° C. (for example, 10° C.) below the melting temperature (Tm) of the native duplex of the cDNA that has the lowest melting temperature (see generally, Sambrook et al., Molecular Cloning: A Laboratory Manual, 2nd ed., Cold Spring Harbor Press, 1987; Ausubel et al., Current Protocols in Molecular Biology, Greene Publishing, 1987). Tm for nucleic acid molecules greater than about 100 bases can be calculated by the formula Tm=81.5+0.41% (G+C)−log(Na+). For oligonucleotide molecules less than 100 bases in length, exemplary hybridization conditions are 5° to 10° C. below Tm.
In one embodiment, the microarray is an array of polynucleotide probes, the array comprising a support with at least one surface and typically at least 100 different polynucleotide probes, each different polynucleotide probe comprising a different nucleotide sequence and being attached to the surface of the support in a different location on the surface. For example, the nucleotide sequence of each of the different polynucleotide probes can be in the range of 40 to 80 nucleotides in length. For example, the nucleotide sequence of each of the different polynucleotide probes can be in the range of 50 to 70 nucleotides in length. For example, the nucleotide sequence of each of the different polynucleotide probes can be in the range of 50 to 60 nucleotides in length. In specific embodiments, the array comprises polynucleotide probes of at least 2,000, 4,000, 10,000, 15,000, 20,000, 50,000, 80,000, or 100,000 different nucleotide sequences.
Thus, the array can include polynucleotide probes for most, or all, genes expressed in a cell, tissue, organ or organism. In a specific embodiment, the cell or organism is a mammalian cell or organism. In another specific embodiment, the cell or organism is a human cell or organism. In specific embodiments, the nucleotide sequences of the different polynucleotide probes of the array are specific for at least 50%, at least 75%, at least 80%, at least 85%, at least 90%, at least 95%, or at least 99% of the genes in the genome of the cell or organism. Most preferably, the nucleotide sequences of the different polynucleotide probes of the array are specific for all of the genes in the genome of the cell or organism. In specific embodiments, the polynucleotide probes of the array hybridize specifically and distinguishably to at least 10,000, to at least 20,000, to at least 50,000, to at least 80,000, or to at least 100,000 different polynucleotide sequences. In other specific embodiments, the polynucleotide probes of the array hybridize specifically and distinguishably to at least 90%, at least 95%, or at least 99% of the genes or gene transcripts of the genome of a cell or organism. In some embodiments, the polynucleotide probes of the array hybridize specifically and distinguishably to the genes or gene transcripts of the entire genome of a cell or organism. In further embodiments, the microarray may be designed to contain a set of probes corresponding to a set of transcripts known to be modulated in a selected biological pathway.
The use of a two-color fluorescence labeling and detection scheme to define alterations in gene expression has been described, e.g., in Schena et al. 1995, Science 270:467-470, which is incorporated by reference in its entirety for all purposes. An advantage of using cDNA molecules labeled with two different fluorophores is that a direct and internally controlled comparison of the mRNA levels corresponding to each arrayed gene in two cell states can be made, and variations due to minor differences in experimental conditions (e.g., hybridization conditions) will not affect subsequent analyses. However, it will be recognized that it is also possible to use cDNA molecules from a single cell, and compare, for example, the absolute amount of a particular mRNA in, e.g., a drug-treated or an untreated cell.
Exemplary microarrays and methods for their manufacture and use are set forth in Hughes, T. R., et al., Nature Biotechnology 19:342-347 (April 2001), which publication is incorporated herein by reference.
Signal detection and data analysis. When fluorescently labeled probes are used, the fluorescence emissions at each site of an array can be detected by scanning confocal laser microscopy. In one embodiment, a separate scan, using the appropriate excitation line, is carried out for each of the two fluorophores used. Alternatively, a laser can be used that allows simultaneous specimen illumination at wavelengths specific to the two fluorophores and emissions from the two fluorophores can be analyzed simultaneously (see Shalon et al. 1996, Genome Research 6:639-645, which is incorporated by reference in its entirety for all purposes). In one embodiment, the arrays are scanned with a laser fluorescent scanner with a computer controlled X-Y stage and a microscope objective. Sequential excitation of the two fluorophores is achieved with a multi-line, mixed gas laser and the emitted light is split by wavelength and detected with two photomultiplier tubes. Fluorescence laser scanning devices are described in Shalon et al. Genome Res. 6:639-645 (1996), and in other references cited herein. Alternatively, the fiber-optic bundle described by Ferguson et al. Nature Biotechnol. 14:1681-1684 (1996), may be used to monitor mRNA abundance levels at a large number of sites simultaneously.
Signals are recorded and may be analyzed by computer, e.g., using a 12 bit analog to digital board. In some embodiments the scanned image is despeckled using a graphics program (e.g., Hijaak Graphics Suite) and then analyzed using an image gridding program that creates a spreadsheet of the average hybridization at each wavelength at each site. If necessary, an experimentally determined correction for “cross talk” (or overlap) between the channels for the two fluors may be made. For any particular hybridization site on the transcript array, a ratio of the emission of the two fluorophores can be calculated. The ratio is independent of the absolute expression level of the cognate gene, but is useful for genes whose expression is significantly modulated by drug administration.
The relative abundance of an mRNA in two biological samples is scored as a perturbation and its magnitude determined (i.e., the abundance is different in the two sources of mRNA tested), or as not perturbed (i.e., the relative abundance is the same). Preferably, in addition to identifying a perturbation as positive or negative, it is advantageous to determine the magnitude of the perturbation. This can be carried out, as noted above, by calculating the ratio of the emission of the two fluorophores used for differential labeling, or by analogous methods that will be readily apparent to those of skill in the art.
By way of example, two samples, each labeled with a different fluor, are hybridized simultaneously to permit differential expression measurements. If neither sample hybridizes to a given spot in the array, no fluorescence will be seen. If only one hybridizes to a given spot, the color of the resulting fluorescence will correspond to that of the fluor used to label the hybridizing sample (for example, green if the sample was labeled with Cy3, or red, if the sample was labeled with Cy5). If both samples hybridize to the same spot, an intermediate color is produced (for example, yellow if the samples were labeled with fluorescein and rhodamine). Then, applying methods of pattern recognition and data analysis known in the art, it is possible to quantify differences in gene expression between the samples. Methods of pattern recognition and data analysis are described in e.g., International Publication WO 00/24936, which is incorporated by reference herein.
Comparison of Gene Expression Levels: In accordance with the methods of the invention, a set of genes are identified that are responsive to the microRNA species, wherein each individual gene has an expression value in response to the modulated level of microRNA that is significantly different from the corresponding reference expression value. An expression value is considered to display a statistically significant difference (either increased or decreased) in expression level (p<0.05) relative to a control, such as a mock-transfected cells.
Art-recognized statistical techniques can be used to compare the levels of expression of individual genes, or proteins, to identify genes, which exhibit significantly regulated expression levels in microRNA or siRNA transfected cells as compared to untreated cells. Thus, for example, a t-test can be used to determine whether the mean value of repeated measurements of the level of expression of a particular gene, is significantly regulated in a cell transfected with a particular microRNA or siRNA molecule compared to the same cell type that has not been transfected. Similarly, Analysis of Variance (ANOVA) can be used to compare the mean values of two or more populations (e.g., two or more populations of cultured cells treated with different amounts of a candidate drug) to determine whether the means are significantly different.
The following publications describe examples of art-recognized techniques that can be used to compare the levels of expression of individual genes, or proteins, in treated and untreated living things, or in diseased and non-diseased living things, to identify genes which exhibit significantly different expression levels: Nature Genetics, Vol. 32, pp. 461-552 (supplement December 2002); Bioinformatics 18(4):546-54 (April 2002); Dudoit et al. Technical Report 578, University of California at Berkeley; Tusher et al., Proc. Nat'l. Acad. Sci. U.S.A. 98(9):5116-5121 (April 2001); and Kerr et al., J. Comput. Biol. 7:819-837.
Representative examples of other statistical tests that are useful in the practice of the present invention include the chi squared test which can be used, for example, to test for association between two factors (e.g., transcriptional induction, or repression, by a drug molecule and positive or negative correlation with the presence of a disease state). Again by way of example, art-recognized correlation analysis techniques can be used to test whether a correlation exists between two sets of measurements (e.g., between gene expression and disease state). For example, an error model may be used to analyze microarray data as described by Weng L. et al., Bioinformatics 22(9): 111-21 (2006), incorporated herein by reference. Standard statistical techniques can be also be found in statistical texts, such as Modern Elementary Statistics, John E. Freund, 7th edition, published by Prentice-Hall; and Practical Statistics for Environmental and Biological Scientists, John Townend, published by John Wiley & Sons, Ltd.
Identification of at Least One microRNA Responsive Target Site in the Transcripts of the Set of microRNA Responsive Genes
In accordance with the method of the invention, transcripts within the set of microRNA responsive genes identified in the microRNA modulated cells are analyzed for the presence of at least one microRNA responsive target site corresponding to the modulated microRNA species. A microRNA responsive target site refers to a nucleic acid sequence ranging in size from about 5 to about 25 nucleotides that corresponds to the modulated microRNA. The microRNA responsive target site may be complementary (i.e., a direct match) or essentially complementary to at least a portion of the modulated microRNA molecule, or a family member thereof. In some embodiments, the microRNA responsive target site comprises a sequence that is complementary (a direct match) to at least 6 consecutive nucleotides (referred to as a hexamer) of the 5′ portion of the microRNA molecule from nucleotide positions 1 to 12 (referred to as the seed region). Exemplary seed regions for several microRNA families are provided in
The regulated transcripts may be analyzed for the presence of a microRNA responsive target site by querying a database of stored sequence information, or may be directly analyzed by hybridization to a probe, PCR amplification or direct sequencing.
The microRNA responsive target site may be present in any region of the downregulated transcript, such as the 5′ UTR, coding region, or 3′ UTR. It has been shown that microRNAs inhibit the expression of target genes by interacting with complementary sites in the 3′ untranslated region (UTR) of the target mRNAs (see Olsen et al., Dev Biol 216:671-680 (1999); Bartel et al., Cell 116:281-297 (2004)). However, other regions, including coding region recognition have also been observed (Brennecke et al., PloS Biol 3(3):pe85 (2005)). Accordingly, in some embodiments, the method includes the step of identifying at least one microRNA responsive target region in the 3′ UTR of the regulated transcript. In some embodiments, the method includes the step of identifying at least two sequences corresponding to at least two seed regions in the 3′ UTR of the regulated transcript. For example, as described in EXAMPLE 1 and TABLE 3, it has been observed that a significant percentage of downregulated transcripts identified in miR-16 transfected cells contained at least one, and often two, sequences corresponding to the miR-16 consensus seed sequence.
Once a plurality of regulated genes are identified that comprise at least one microRNA responsive target site, a plurality of members of the set of genes are modulated in the cell type of interest with a plurality of gene-specific agents. The gene-specific agents are designed to specifically modulate (either increase or decrease expression of) the identified micro-RNA-responsive genes to determine whether inhibition (or overexpression) of the individual genes targeted, either alone or in combination, produces the same expression pattern or phenotype as the corresponding microRNA. In other words, the set of identified micro-RNA responsive genes are individually targeted to see if when modulated, alone or in combination, one or more phenotypes observed in the microRNA modulated cells are phenocopied.
The gene-specific agents used to modulate the identified micro-RNA-responsive genes may be any suitable agonist or antagonist of the identified gene. As a non-limiting illustrative example, in the case where a selected microRNA species was transfected into a cell type of interest, a set of microRNA-responsive genes are identified that exhibit decreased expression, and a plurality of gene-specific inhibitory agents directed against the identified set of genes would be chosen to transfect into the cell type of interest. The phenotype of the microRNA transfected cells and the inhibitory gene-specific agent transfected cells would then be compared for the presence of at least one shared phenotype.
Gene-specific agents designed to inhibit a gene of interest include RNA inhibitors such as antisense oligonucleotides, iRNA agents, and protein inhibitors, such as antibodies, soluble receptors, and the like. IRNA agents encompass any RNA agent which can downregulate the expression of a target gene, including siRNA molecules and shRNA molecules. The siRNA molecules may be designed to inhibit a particular target gene by using an algorithm developed to increase efficiency of the siRNAs for silencing while minimizing their off-target effects, as described in Jackson et al., Nat Biotech 21:635-637 (2003), WO2006006948, and WO2005/042708, incorporated herein by reference. As described in Jackson et al., RNA 12: 1197-1205 (2006), incorporated herein by reference, and also in U.S. Patent Publication No. 20050223427, incorporated herein by reference, 2′-O-methyl modifications to specific positions within the siRNA seed region (e.g., at positions 1-5) reduce both the number of off-target transcripts and the magnitude of their regulation, while not significantly affecting silencing of the intended targets. In particular, it was observed by Jackson et al. (2006), that 2′-O-methyl ribosyl substitution in the siRNA at position 1, position 2, or positions 1+2 of the seed region were effective to reduce the off-target transcript down-regulation whiel regulation of the intended target was unaffected. The addition of O-methyl groups to the 2′ position of the ribosyl ring may be carried out using techniques known in the art, and is commonly used for RNA modification to alter thermodynamic and binding properties of modified duplexes. See e.g., Monia et al., J. Biol. Chem. 268: 14514-14522 (1993).
Exemplary siRNA sequences designed to target miR-16 downregulated transcripts are provided below in TABLE 9.
The microRNA, and iRNA agents (including shRNA, and siRNA molecules) for use in the practice of the methods of the invention and to produce the compositions of the invention may chemically synthesized or recombinantly produced using methods known in the art. for example, the RNA products may be chemically synthesized using appropriately protected ribonucleoside phosphoramidites and a conventional DNA/RNA synthesizer. Commercial suppliers of synthetic RNA molecules or synthesis reagents include Proligo (Hamburg Germany) and Dharmacon Research (Lafayette, Colo.). Exemplary microRNA molecules that may be used to practice various embodiments of the methods of this aspect of the invention are provided in TABLE 1 and TABLE 7.
In one embodiment, the invention provides chemically modified siRNA constructs designed to target miR-16 responsive targets. For example, the siRNA agent can include a non-nucleotide moiety. A chemical modification or other non-nucleotide moiety can stabilize the sense (guide strand) and antisense (passenger strand) sequences against nucleolytic degradation. Additionally, conjugates can be used to increase uptake and target uptake of the siRNA agent to particular cell types. Thus, in one embodiment the siRNA agent includes a duplex molecule wherein one or more sequences of the duplex molecule is chemically modified. Non-limiting examples of such chemical modifications include phosphorothioate internucleotide linkages, 2′-deoxyribonucleotides, 2′-O-methyl ribonucleotides, 2′-deoxy-2′-fluoro ribonucleotides, “universal base” nucleotides, “acyclic” nucleotides, 5′-C-methyl nucleotides, and terminal glyceryl and/or inverted deoxy abasic residue incorporation. These chemical modifications, when used in siRNA agents, can help to preserve RNAi activity of the agents in cells and can increase the serum stability of the siRNA agents.
In one embodiment, the first and optionally or preferably the first two internucleotide linkages at the 5′ end of the antisense and/or sense sequences are modified, preferably by a phosphorothioate. In another embodiment, the first, and perhaps the first two, three, or four internucleotide linkages at the 3′ end of a sense and/or antisense sequence are modified, for example, by a phosphorothioate. In another embodiment, the 5′ end of both the sense and antisense sequences, and the 3′ end of both the sense and antisense sequences are modified as described.
Alternatively, microRNA gene products and IRNA agents can be expressed from recombinant circular or linear DNA plasmids using any suitable promoter. Suitable promoters for expressing RNA from a plasmid include the U6 or H1 RNA polIII promoter sequences, or the cytomegalovirus promoters. Selection of other suitable promoters is within the skill in the art. The recombinant plasmids of the invention can also comprise inducible or regulatable promoters for expression of the microRNA or IRNA agent gene products in a desired cell type. For example, a vector may be designed to drive expression (e.g., using the PolIII promoter) of both the sense and antisense strands seperately, which hybridize in vivo to generate siRNA.
In one embodiment, the iRNA agent is an shRNA. A vector may be used to drive expression of short hairpin RNA (shRNA), which are individual transcripts that adopt stem-loop structures, which are processed into siRNAs by the RNAi machinery in the cell. Typically, the shRNA design comprises two inverted repeats containing the sense and antisense target sequence separated by a loop sequence. Typically, the loop sequence contains 8 to 9 bases. A terminator sequence consisting of 5-6 polydTs is present at the 3′ end and one or more cloning sequences may be added to the 5′ end using complementary oligonucleotides. A website is available for design of such vectors, see, http://www.genelink.com/sirna/shRNAhelp.asp.
An shRNA vector may be designed with an inducible promoter. For example, a lentiviral vector may be used expressing tTS (tetracycline-controlled transcriptional repressor, Clontech). For example, a tetracycline-inducible shRNA designed to target a gene, such as PLK1 may be driven from an H1 promoter, as described in Jackson et al., RNA 12:1-9 (2006). The cells of interest are infected with recombinant lentivirus and shRNA expression is induced by incubation of the cells in the presence of 50 ng/mL of doxycycline.
In some embodiments, the method of this aspect of the invention further comprises the step of comparing the set of identified microRNA responsive transcripts to a source of biological knowledge to identify a plurality of microRNA responsive transcripts that are enriched in a known biological pathway or process. As used herein, the term “source of biological knowledge” refers to information that describes the function (e.g., at molecular, cellular and system levels), structure, pathological roles, toxicological implications, etc., of a multiplicity of genes. Various sources of biological knowledge can be used for the methods of the invention, including databases and information collected from public sources such as Locuslink, Unigene, SwissTrEMBL, etc., and organized into a relational database following the concept of the central dogma of molecular biology.
In some embodiments, the annotation systems used by the Gene Ontology (GO) Consortium or similar systems are employed. GO is a dynamic controlled vocabulary for molecular biology which can be applied to all organisms as knowledge of gene is accumulating and changing; it is developed and maintained by Gene Ontology™ Consortium (Gene Ontology: Tool for the unification of biology. The Gene Ontology Consortium (2000) Nature Genet. 25:25-29). Gene annotations using GO terms provide an excellent resource for summarized knowledge on each gene. Genes with similar biological property are annotated with the same or similar GO terms and thus can be easily identified. Currently, there are three categories of GO terms: biological processes, molecular function, and cellular component. For example, the Degenerin gene is annotated with “peripheral nervous system development,” “monovalent inorganic cation transport,” “central nervous system development,” and “synaptic transmission” for biological process, “amiloride-sensitive sodium channel” for molecular function, and “integral plasma membrane protein” for cellular component.
GO cluster algorithm may be employed in computer implemented methods to automatically cluster genes based upon existing knowledge. One of skill in the art may determine that a correlation (positive or negative) exists between the expression pattern of the microRNA-modulated cells and the enrichment of regulated (i.e. responsive) transcripts in an identified biological pathway using GO terms.
In accordance with the methods of this aspect of the invention, the cells that are seperately modulated with a plurality of gene-specific agents directed against different identified miRNA responsive genes are analyzed for the presence of at least one phenotype of interest that is also displayed by cells modulated with the corresponding microRNA. The shared phenotype may be any measurable effect, or combination of measurable effects, including cellular phenotypes and expression phenotypes. For example, a phenotype can be determined with reference to stimulation, and/or inhibition, of one or more biological responses; and/or the absolute and/or relative magnitude of stimulation, and/or inhibition, of one, or more, biological responses; and/or the inability to affect (e.g., the inability to stimulate or inhibit) one, or more, biological responses. Non-limiting examples of measurable biological processes include biochemical pathways; physiological processes that contribute to the internal homeostasis of a living organism; developmental processes that contribute to the normal physical development of a living organism; and acute or chronic diseases. For example, cellular morphology may be analyzed using automated fluorescence microscopy to visualize actin filaments, microtubules and DNA, in order to detected morphological cellular phenotypes. In some embodiments, a phenotype may be measurable by performing gene or protein expression analysis on a set of predetermined transcripts or proteins. For example, a cell cycle phenotype may be measured by analyzing cell cycle progression as described in EXAMPLE 4.
In one embodiment of the method, the phenotype is a gene expression pattern induced by modulating the microRNA. In accordance with this embodiment, a gene expression profile for cells transfected with each of a plurality of the gene-specific agents targeting a plurality of microRNA responsive transcript are compared to the microRNA transfected cell profile. The gene expression patterns triggered by the microRNA and each gene-specific agent targeting a micro-RNA responsive transcript are statistically analyzed to identify the subset of genes that share the phenotype of triggering a similar expression profile in the cell type of interest.
In another embodiment of the method, a cellular phenotype is observed in transfected cells, such as a cell cycle regulatory phenotype. For example, the transfected cells can be analyzed for various properties such as viability and cell cycle distribution at various time points, as described in more detail in EXAMPLE 4.
In accordance with some embodiments, the step of comparing the down-regulated transcripts to a source of biological knowledge to identify a plurality of down-regulated transcripts that are enriched in a known biological pathway or process is used to select a suitable phenotype for use in the method of the invention. It will be appreciated by those of skill in the art that the method of the invention is not limited to a particular sequence of steps described herein, but may be practiced using the described steps in a different order, and/or with the addition or elimination of certain steps.
The methods of this aspect of the invention have been used to identify a subset of genes that are responsive to miR-16 and share a cell cycle phenotype, as described in more detail in EXAMPLES 1-7 and
In another aspect, the present invention provides a method of inhibiting proliferation of a mammalian cell comprising introducing into said cell an effective amount of at least one small interfering RNA agent (iRNAi), wherein said iRNA comprises a nucleotide sequence of at least 15 nucleotides, wherein the nucleotide sequence comprises a seed region consisting of nucleotide positions 1 to 12, wherein position 1 represents the 5′ end of the iRNA nucleotide sequence and wherein said seed region comprises a nucleotide seqeuence of at least six contiguous nucleotides that is complementary to six contiguous nucleotides located within positions 1 to 12 of a nucleotide sequence (including position 1 and position 12), wherein position 1 represents the 5′ end of the nucleotide sequence, wherein the nucleotide sequence is selected from the group consisting of SEQ ID NO: 1, SEQ ID NO:2, SEQ ID NO:3, SEQ ID NO:4, SEQ ID NO:5, SEQ ID NO:6, SEQ ID NO:7 and SEQ ID NO:8.
In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:1. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:2. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:3. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:4. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:5. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:6. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:7. In some embodiments, the nucleotide sequence is within positions 1 to 12 of SEQ ID NO:8.
In some embodiments, the iRNA is a gene-specific inhibitor of expression of at least one miR-16 responsive gene selected from TABLE 5. In some embodiments, the method comprises introducing an effective amount of an iRNA that inhibits expression of at least two miR-16 responsive genes selected from TABLE 5 into the mammalian cell. In some embodiments, the method comprises introducing an effective amount of at least one iRNA that inhibits expression of at least one miR-16 responsive gene selected from TABLE 6 into the mammalian cell.
In some embodiments, the at least one miR-16 response gene is selected from the group consisting of ATG9A, ATXN7L3, C10orf46, IPPK, C9orf42, C9orf91, CARD10, CBX6, CDC27, CDK6, COX10, H2AFX, KIAA0317, MFN2, PHF17, PPP1R11, RAB11FIP2, and SRPR. In some embodiments, the method comprises introducing a composition comprising an effective amount of a combination of iRNA agents, such as nucleic acid molecules that inhibit at least two or more miR-16 responsive targets selected from the group consisting CARD10, CDC27, CDK6 and C10orf46.
The methods of this aspect of the invention may be used to inhibit proliferation of a cancer cell, such as an miR-16 mediated cancer cell, such as CLL or prostate cancer cells. The methods of this aspect of the invention may also be used to inhibit proliferation of a mammalian cell infected with a virus expressing a microRNA species in the miR-16 family, such as, for example, a cell infected with kshv-miR-K12-6-5p.
In some embodiments, the iRNA agents are selected from the group consisting of siRNA molecules and shRNA molecules. Exemplary siRNA molecules useful in the practice of the method of the invention are provided in TABLE 8 and TABLE 10. In some embodiments, the siRNA molecules comprise at least one dsRNA molecule comprising one nucleotide strand that is substantially identical to a portion of the mRNA encoding a gene listed in TABLE 6, such as, for example, CDK6, CDC27, CARD10 and C10orf46. In one particular embodiment, the gene-specific agent of at least one miR-16 responsive gene is at least one dsRNA molecule comprising a double-stranded region, wherein one strand of the double-stranded region is substantially identical to 15 to 25 consecutive nucleotides of an mRNA encoding a gene set forth in TABLE 6, and the second strand is substantially complementary to the first, and wherein at least one end of the dsRNA has an overhang of 1 to 4 nucleotides.
In one embodiment, the iRNA agent comprises at least one of SEQ ID NO:272-277, SEQ ID NO:284-289, SEQ ID NO:290-295, SEQ ID NO:248-253, SEQ ID NO:344-349, SEQ ID NO:350-355, SEQ ID NO:356-361.
The siRNAs useful in the methods of the invention may be chemically synthesized and annealed before delivery to a cell or mammalian subject, as described supra. In some embodiments, the siRNAs are synthesized in vivo, such as from a plasmid expression system (see, e.g., Tuschl and Borkhardt, Molec Interventions 2:158-167 (2002)). Exemplary constructs for making dsRNAs are described, for example, in U.S. Pat. No. 6,573,099. In some embodiments, the siRNA or shRNA inhibitory molecules inhibit expression of a target gene by at least 30%, such as 50%, such as 60%, such as 80%, such 90% up to 100%.
The siRNA and shRNA molecules can be delivered into cells in culture using electroporation or lipophilic reagents. The siRNA molecules can be delivered into a mammalian subject, for example, by intravenous injection, direct injection into a target site (e.g., into tumors), or into mice or rats by high-pressure tail-vein injection. It has been demonstrated that synthetic siRNAs can silence target gene expression in mammalian models. For example, McCaffrey et al. (Nature 418:38-39 (2002)), described silencing of a reporter gene in mice when the reporter gene and siRNA were injected simultaneously by high-pressure tail vein injections. Moreover, Soutsched et al. (Nature 432:173-178 (2004)) demonstrated that a synthetic siRNA downregulated expression of an endogenous target gene following intravenous injection in mice. Similarly, Pulukuir et al. (J. Biol. Chem 280:36529-36540 (2005)), demonstrated that injection of plasmids expressing short hairpin RNAs (shRNAs) into tumors in mice downregulated expression of the target gene in the tumors and also caused a decrease in tumor weight.
In another aspect, the present invention provides compositions comprising a combination of nucleic acid molecules that are useful as inhibitors of at least two or more miR-16 responsive targets selected from TABLE 5 or TABLE 6. In some embodiments, the compositions comprise a combination of nucleic acid molecules that are useful as inhibitors of at least two or more miR-16 responsive targets selected from the group consisting of ATG9A, ATXN7L3, C10orf46, IPPK, C9orf42, C9orf91, CARD10, CBX6, CDC27, CDK6, COX10, H2AFX, KIAA0317, MFN2, PHF17, PPP1R11, RAB11FIP2, and SRPR. In some embodiments, the compositions comprise a combination of nucleic acid molecules that are useful as inhibitors of at least two or more coordinately regulated miR-16 responsive targets selected from the group consisting CARD10, CDC27, CDK6 and C10orf46.
In another aspect, the present invention provides an isolated dsRNA molecule comprising one nucleotide strand that is substantially identical to a sequence selected from the group consisting of SEQ ID NO:236 to SEQ ID NO:361. In some embodiments, the isolated dsRNA molecule comprises at least one of SEQ ID NO:236 to SEQ ID NO: 361. In some embodiments, at least one strand of the isolated dsRNA molecule consists of at least one of SEQ ID NO:236 to SEQ ID NO:361.
In another embodiment, pharmaceutical compositions comprising nucleic acid molecules that inhibit at least one miR-16 responsive target are provided. Such a composition contains from about 0.01 to 90% by weight (such as 1 to 20% or 1 to 10%) of a therapeutic agent of the invention in a pharmaceutically acceptable carrier. Solid formulations of the compositions for oral administration may contain suitable carriers or excipients, such as corn starch, gelatin, lactose, acacia, sucrose, microcrystalline cellulose, kaolin, mannitol, dicalcium phosphate, calcium carbonate, sodium chloride, or alginic acid. Liquid formulations of the compositions for oral administration prepared in water or other aqueous vehicles may contain various suspending agents such as methylcellulose, alginate, tragacnth, pectin, kelgin, carrageenan, acacia, polyvinylpyrrolidone, and polyvinyl alcohol.
Injectable formulations of the compositions may contain various carriers such as vegetable oils, dimethylacetamide, dimethylformamide, ethyl lactate, ethyl carbonate, isopropyl myristate, ethanol, polyols (glycerol, propylene glycol, liquid polyethylene glycol and the like). For intravenous injections, water soluble versions of the compounds may be administered by the drip method, whereby a pharmaceutical formulation containing an antifungal agent and a physiologically acceptable excipient is infused. Physiologically acceptable excipients may include, for example, 5% dextrose, 0.9% saline, Ringer's solution, or other suitable excipients. Intramuscular preparations, e.g., a sterile formulation of the compounds of the invention can be dissolved and administered in a pharmaceutical excipient was as water-for-injection, 0.9% saline, or 5% glucose solution.
Conventional methods, known to those of ordinary skill in the art of medicine, can be used to administer the pharmaceutical formulations to a mammalian subject. The pharmaceutical formulations can be administered via oral, subcutaneous, intrapulmonary, transmucosal, intraperitoneal, intrauterine, sublingual, intrathecal or intramuscular routes.
The following examples merely illustrate the best mode now contemplated for practicing the invention, but should not be construed to limit the invention. All literature citations are expressly incorporated by reference.
This Example describes a method of using microarray expression profiling of cells transfected with a panel of microRNAs duplex molecules to identify microRNA responsive targets that are coordinately down-regulated.
Methods and Materials:
24 microRNAs were selected for transfection, as shown below in TABLE 1.
AGCAGCAUUGUACAGGGCUAUGA
AUAGCCCUGUACAAUGCUGGUAU
AGCUGGUUGAAGGGGAUUAAAU
AACACUGUCUGGUAAAGAUGG
AUCUUUACCAGACAGUGAUAU
ACUGCAGUGAAGGCACUUGU
AAGUGCCUUCACUGCACUAU
AGUUUUGCAUAGAUUUGCAGAAU
AGUUUUGCAUGGAUUUGCAGAAU
AUCGUUACCAGACAGUGUAAAU
AACAUCAGUCUGAUAAGCAAAU
AUGACCUAUGAAUUGACAGAC
AGCUACAUUGUCUGCUGGGUUUC
AACCCAGCAGACAAUGUAGGUAU
AGACCGAGACAAGUGCAAAGAU
AGGCCGGGACAAGUGCAAUAAU
AAAGUGCUGUUCGUGCAGGUAG
ACCUGCACGAACAGCACUUUAU
Cell lines: Each of the 24 different miRNAs listed in TABLE 1 were transfected into HCT116 Dicerex5 and DLD-1 Dicerex5 colon tumor cells described in Cummins, J. M., et al. (PNAS 103(10):3687-3692 (2006)). Briefly described, these colon cancer tumor cell lines have homozygous disruption of the Dicer helicase domain and show reduced levels of many endogenous miRNAs (Cummins, supra). A preliminary set of experiments were carried out in matched Dicer wild type cells. It was observed that the Dicer homomorphs showed substantially similar expression patterns, with approximately 2-fold more intense expression changes following transfection of exogeneous siRNA/miRNAs than matched Dicer wild type cells (data not shown), therefore, the Dicerex5 cells were used for subsequent experiments. Wild type HCT116, DLD-1, HeLa, A549, MCF7, and TOV21G cells were obtained from the ATCC, Rockville, Md.
RNA Duplexes: RNA duplexes corresponding to mature miRNAs were designed as described in Lim, L. P., et al., Nature 433(7027):769-73 (2005), incorporated herein by reference. The miRNA duplexes used in this study are shown in TABLE 1. miRNA duplexes were ordered and produced by Dharmacon (Lafayette, Colo.).
Transfections: For miRNA transfections, cells were plated 24 hours prior to transfection. HCT116 cells were transfected in 6 well plates using Lipofectamine 2000 (Invitrogen, Carlsbad, Calif.). DLD-1, HeLa, TOV21G and A549 cells were transfected using SilentFect (Bio-Rad, Hercules, Calif.). MCF7 cells were transfected using Oligofectamine (Invitrogen, Carlsbad, Calif.). miRNA duplexes were annealed and were individually transfected at a final concentration of 25 nM for all cell lines except for HCT166, where concentration of 100 nM was used. Semi-automated 384 well transfections were carried out as described in Bartz et al., submitted.
Microarray Analysis: Total RNA was isolated from transfected cells 24 hours post transfection. Total RNA was purified using an RNAeasy kit (Qiagen). cRNA amplification and microarray analysis was performed as described in Hughes et al., Nat Biotech 19:342-347 (2001), incorporated herein by reference. Amplified cRNA from miRNA-transfected cells was hybridized against cRNA from mock-transfected cells (treated with transfection reagent in the absence of RNA duplex). Ratio hybridizations were performed with fluorescent label reversal to eliminate dye bias. Expression profile analysis and statistical error models were used as described in Hughes et al., supra. Data were analyzed using Rosetta Resolver™ (Rosetta Biosoftware, Seattle, Wash.).
Identification of Downregulated Transcripts: Gene expression profiles were determined at 24 hours post-transfection by competitive hybridization of amplified RNA from miRNA-transfected HCT116Dicerex5 cells and DLD-1 Dicerex5 cells versus mock treated cells. The 24 hour time point was chosen because it has been shown that by 24 hours mRNA silencing is maximal, but secondary transcription effects due to protein depletion are generally minimal (see Jackson et al., Nat. Biotech. 21:635-637 (2003)). Expression profiles were analyzed from the two individual cell lines transfected with each miRNA, and also as a union, or intersection of profiles from the two transfected cell lines. The intersection signature (p<0.01) between the HCT116Dicerex5 cells and DLD-1 Dicerex5 transfected cell lines was identified. It was determined that the intersection signatures yielded a more consistent set of targets, likely due to the cancellation of artifacts in the hybridization data. A collection of miRNA downregulated transcripts (p<0.05) was defined by the intersection of downregulated transcripts from two different Dicerex5 cell lines.
Results of Microarray analysis: Transcripts that were downregulated after transfection with miRNA were identified in microarray gene expression data as downregulated transcripts using a P value cut-off (P<0.05). It was observed that transfection of miRNAs sharing seed region identity resulted in patterns of downregulated transcripts were largely unique to a family of miRNAs sharing identity in residues 1-12 of the miRNA guide strand (seed region). For example, nearly identical expression profiles were obtained for the following miRNA species: (1) miR-15a (SEQ ID NO:1), miR-16 (SEQ ID NO:3), and miR-103 (SEQ ID NO:4), (2) miR-17-5p (SEQ ID NO:9), miR-20a (SEQ ID NO:10), (3) miR-141 (SEQ ID NO: 12) and miR-200a (SEQ ID NO: 13) and (4) miR-192 (SEQ ID NO: 14) and miR-215 (SEQ ID NO: 15). The seed region sequence for the exemplary miRNAs listed above is shown in TABLE 2 and
With regard to miR-15a (SEQ ID NO:1), miR15b (SEQ ID NO:2) and miR-16 (SEQ ID NO:3), it was determined that this set of miRNA species share an identical hexamer seed region (e.g., a stretch of 6 contiguous bases complementary to a seed region of nucleotide 1-6, 2-7, or 3-8), and were designated as an “miR-16 family.” It was interesting to note that miR-103 (SEQ ID NO:4) gave an overlapping, but distinct profile from miR-15a and miR-16 (data not shown). As shown in
miRNA Downregulated Gene Set Analysis:
The set of miRNA downregulated transcripts identified by microarray analysis were then tested for enrichment with respect to the presence of 1) one or more miRNA hexamer seed regions in the 3′ UTR (i.e., a stretch of 6 contiguous bases complementary to a seed region nt 1-6, 2-7, or 3-8) corresponding to one of the miRNA species in the set of transfected miRNAs, and/or 2) transcript annotation information regarding the biological function of the transcript, for example, annotation of a transcript in the public Gene Ontology Biological Process database. The Gene Ontology Biological Process database provides structure controlled Gene Ontologies (GO) that describe gene products in terms of their associated biological processes, cellular components, and molecular functions in a species-independent manner. The categories of GO terms include biological processes, molecular function, and cellular components. The Gene Ontology (GO) Database is maintained by the Gene Ontology Consortium (2000), Nature Genet. 25:25-29 (2000).
Analysis of miRNA Downregulated Gene Set for miRNA Targets
The transcripts that were found to be downregulated were examined for the presence of one or more seed hexamer matches in a 3′ UTR corresponding to one of the miRNA species transfected into the cells relative to a background set of genes on the microarray using a hypergeometric distribution (see, e.g., Benjamini, Y. & Hochberg, Y. (1995), J.R. Stat. Soc. B 57:289-300). 23 of the 24 miRNA downregulated intersection signatures showed enrichment for seed region hexamers, a result that is highly unlikely to occur by chance (E<1E-20). The family-specific transcripts were highly enriched for hexamer sequence motifs complementary to (matching) the seed region for that family (data not shown). It was also determined that the downregulated signatures were significantly enriched with respect to computationally predicted miRNA targets (Krek et al., Nat. Genet. 37:495-500 (2005); Xie et al., Nature 434:338-345 (2005)), confirming that the miRNA transfections were largely successful.
UAGCAGCACAUAAUGGUUUGUG
UAGCAGCACAUCAUGGUUUACA
UAGCAGCACGUAAAUAUUGGCG
AGCAGCAUUGUACAGGGCUAUGA
AGCAGCAUUGUACAGGGCUAUCA
UAGCAGCACAGAAAUAUUGGC
CCAGCAGCACCUAAUCCAUCGG
UAGCAGCA
CAAAGUGCUUACAGUGCAGGUAGU
UAAAGUGCUUAUAGUGCAGGUAG
UAAAGUGCUGACAGUGCAGAU
UAACACUGUCUGGUAAAGAUGG
UAACACUGUCUGGUAACGAUGU
CUGACCUAUGAAUUGACAGCC
AUGACCUAUGAAUUGACAGAC
Results of Seed Analysis: The downregulated transcripts were analyzed for enrichment of hexamer motifs complementary to the miRNA seed region in their 3′ UTRs. TABLE 3 shows the data for the miRNA-16 family including miR-103, miR-15a, and miR-16. Of the total 18,124 total transcripts on the microarray, 641 transcripts were identified that were downregulated and contained at least one miR-16 seed region in the 3′UTR.
The miR-16 family of transcripts was further analyzed with regard to timing of transcript downregulation and number of hexamer seed regions present in the 3′ UTR< as described in more detail in EXAMPLE 2. Transcripts in the intersection signature that were also regulated (p<0.05) at 6 hours in HCT Dicerex5 cells were defined as miR-16 consensus downregulated transcripts (N=116), as described in EXAMPLE 2.
Biological Annotation Analysis: The set of miRNA downregulated transcripts identified by microarray analysis were also analyzed to determine if they were enriched for transcripts associated with known biological pathways or processes. A search was done with the transcripts identified in the microarray analysis against the GO Biological Database for correlation between a set of transcripts downregulated with a particular miRNA species and a GO process term.
Results: It was observed that 3 of the 24 miRNAs used for transfection, miR-16, miR-15a, and miR-103 downregulated a set of transcripts that showed a significant correlation (E<1E-2) for transcripts associated with a GO Biological Process term, as shown above in TABLE 4. It is interesting to note that miR-16, miR-15a, and miR-103 are all from the miR-16 family, as shown in
The set of transcripts downregulated in transfections with miR-16, miR-15a and miR-103 were found to be significantly enriched for transcripts annotated with GO process terms “mitotic cell cycle” or “cell cycle,” as shown below in TABLE 4. It was observed that the degree of functional enrichment in miR-103 signatures was less significant than with miR-15a and miR-16. This is consistent with the finding that the gene expression signature resulting from miR-103 transfection overlaps with, but is distinct from, the miR-15a and miR-16 signature, and that the miR-103 seed region is offset by 1 nucleotide, as described above and shown in
This Example describes a kinetic analysis of transcript regulation after transfection with RNA duplexes for miR-16 family members miR-15a, miR-106b, and miR-16 and analysis of miR-16 downregulated transcripts for the presence of miRNA responsive target sites.
Rationale: The miRNA responsive target site analysis and transcript annotation described in EXAMPLE 1 suggested that transcripts regulated by miR-15a, miR-16, and to some extent miR-103 are involved in regulation of the cell cycle and/or cell growth. Both mitotic cell cycle transcripts and putative targets of miR-15a and miR-16 were found to be enriched in the 24-hr down-regulated signatures. In order to distinguish transcripts directly regulated or indirectly regulated by the miR-16 family, a kinetic analysis of transcript regulation was carried as described below.
Methods: HCT116Dicerex5 and DLD-1Dicerex5 cells were transfected with miR-106b, miR-16, and miR-15a duplexes, as described in EXAMPLE 1. RNA samples were isolated from the transfected cells at 0, 6, 10, 14, and 24 hours after transfection and were compared to RNA from mock-transfected cells. Gene expression signatures were determined at each time point using microarray expression analysis as described in EXAMPLE 1.
In a separate experiment, miR-16 was transfected into HCT116 wt (colon cancer), DLD-1wt (colon cancer), A549 (lung cancer), MCF7 (breast cancer), and TOV21G (ovarian cancer) cells. Wild type HCT116, DLD-1, HeLa, A549, MCF7, and TOV21G cells were obtained from the ATCC, Rockville, Md. Transfections were carried out as described in EXAMPLE 1. HCT116 cells were transfected in 6 well plates using Lipofectamine 2000 (Invitrogen, Carlsbad, Calif.). DLD-1, HeLa, and A549 cells were transfected using SilentFect (Bio-Rad, Hercules, Calif.). MCF7 and TOV21G cells were transfected using Oligofectamine (Invitrogen, Carlsbad, Calif.). miRNA duplexes were individually transfected at a final concentration of 25 nM for all cell lines except for HCT166, where concentration of 100 nM was used.
Results:
A second group of transcripts that were initially identified as downregulated in miR-15a/miR-16 transfected cells and annotated as “mitotic cell cycle” related transcripts were downregulated at 24 hours (E>1E-2) in comparison to a luciferase control at 24 hours (E=1.1E-80). This second group of transcripts that is downregulated more slowly (in contrast to the transcripts downregulated after 6 hours) appears to be indirectly affected by miR-16 modulation.
Analysis of miR-16 Family Consensus Transcripts (Downregulated after 6 Hours)
The individual transcripts downregulated after 6 hours (referred to as miR-16 family consensus transcripts), were tested for enrichment for annotation of the GO Biological Process term “cell cycle” and the miR-16 family seed region hexamer beginning at position 2. Of the N=116 transcripts that fell into the category of miR-16 consensus family transcripts, N=84 transcripts were selected for further analysis based on having a sequence record denoting ≧20 nucleotides of both CDS and 3′ UTR. The background set for miR-16 transcriptional target comparisons comprised a set of transcripts that were not found to be down-regulated by miR-16, but having a similar distribution of expression levels and CDS and 3′ UTR lengths (N=1,546).
Results:
The 84 miR-16 consensus transcripts were compared to the expression level-matched set (1546) to identify the presence of hexamer miR-16 seed target sites (positions 1-6, 2-7, or 3-8 of SEQ ID NO:3). Hexamer motifs matching the miR-16 seed target site were found in both coding sequence (CDS) and 3′ UTR regions. The number of transcripts identified with CDS target sites was not found to differ significantly between the miR-16 down-regulated transcripts and the background set (data not shown). However, it was observed that 95% of the miR-16 6 hr downregulated transcripts contained miR-16 seed target sites in their 3′ UTR, in comparison to 45% of controls (P<5E-13).
It was further determined that the miR-16 downregulated transcripts contain multiple copies of miR-16 target sites in their 3′ UTRs. The median number of miR-16 target sites in the 3′ UTR regions of miR-16 consensus transcripts was found to be 2, with the control group having less than 1% occurrence of more than 1 miR-16 target site. The significance the two groups was determined to have a Wilcoxon rank-sum P value<1E-10.
The increased number target sites per transcript for miR-16 down-regulated transcripts was partially attributable to longer 3′ UTRs in these transcripts (median of ˜1,150 nt for miR-16 down regulated transcripts, versus ˜660 nt for the background set, P<1E-2), but more significantly to a greater number of target sites per kilobase of 3′ UTR (2.3 per kb versus 0.2 in controls (P<4E-10)).
Furthermore, the longest target site per transcript was significantly longer for miR-16 down-regulated transcripts (P<5E-14). This was true even when transcripts without hexamer matches were excluded: the median length of the longest site per miR-16 target was 8 bases, vs. 6 bases for background transcripts with hexamer matches (P<3E-10). Thus, many miR-16 down-regulated transcripts have multiple target sites matching the miR-16 seed region, generally with at least one site showing extended complementarity to the miR-16 seed.
It was observed that the 3′ UTRs of miR-16 downregulated transcripts were longer (median of 1200 nt) than the 3′ UTR regions from background transcripts, with a median value of 600 nt. The significance of the difference between the groups was determined to have a Wilcoxon rank-sum P value<4E-3.
Finally, it was observed that the 3′ UTRs of miR-16 downregulated transcripts had a higher density of miR-16 target sites in their 3′ UTR (about 2% of the 3′ UTR) in comparison to background transcripts (about 0.5% of the 3′ UTR). The significance of the difference between the groups was determined to have a Wilcoxon rank-sum P value <1E-12.
Therefore, these results confirm that both mitotic cell cycle genes and putative targets of miR-15a and miR-16 were enriched in 24 hour down-regulated signatures. Further, these results suggest that miR-15a and miR-16 directly targets cell cycle genes at the G0/G1 phase of the cell cycle and indirectly effects cell cycle regulation during the S and G2/M phase. Therefore, mitotic cell cycle genes in miR-15a and miR-16 downregulated signatures is likely to be a downstream indirect effect of regulation of direct targets of miR-15a and miR-16.
Homo sapiens peroxisomal biogenesis
Homo sapiens exostoses (multiple) 2
Homo sapiens solute carrier family 11
Homo sapiens BCL2-interacting killer
Homo sapiens cyclin E1 (CCNE1),
Homo sapiens cell division cycle 27
Homo sapiens CHK1 checkpoint homolog
Homo sapiens COX10 homolog,
Homo sapiens
Homo sapiens H2A histone family,
Homo sapiens solute carrier family 9
Homo sapiens signal recognition particle
Homo sapiens WEE1 homolog (S. pombe)
Homo sapiens CDP-diacylglycerol
Homo sapiens TAR (HIV) RNA binding
Homo sapiens Rho GDP dissociation
Homo sapiens breakpoint cluster region
Homo sapiens ribosomal protein S6
Homo sapiens kinesin family member 3B
Homo sapiens sperm associated antigen 7
Homo sapiens RCE1 homolog, prenyl
Homo sapiens heat shock 70 kDa
Homo sapiens protein-L-isoaspartate
Homo sapiens transforming, acidic
Homo sapiens ubiquitin specific
Homo sapiens chromosome 1 open reading
Homo sapiens kinesin family member 1C
Homo sapiens WW domain containing E3
Homo sapiens lysophospholipase II
Homo sapiens tousled-like kinase 1
Homo sapiens coiled-coil domain
Homo sapiens nin one binding protein
Homo sapiens chromobox homolog 6
Homo sapiens caspase recruitment domain
Homo sapiens KIAA0317 (KIAA0317),
Homo sapiens mitofusin 2 (MFN2),
Homo sapiens RAB11 family interacting
Homo sapiens cytoplasmic linker
Homo sapiens anaphase promoting
Homo sapiens NMD3 homolog
Homo sapiens ubiquitin associated domain
Homo sapiens ring finger protein 138
Homo sapiens NCK interacting protein
Homo sapiens kelch repeat and BTB
Homo sapiens ubiquitin-conjugating
Homo sapiens zinc finger, CCHC domain
Homo sapiens ubiquitin-conjugating
Homo sapiens cell division cycle 37
Homo sapiens cell division cycle
Homo sapiens hypothetical protein
Homo sapiens riboflavin kinase (RFK),
Homo sapiens chromosome 20 open
Homo sapiens vacuolar protein sorting
Homo sapiens anillin, actin binding
Homo sapiens hypothetical protein
Homo sapiens HCV NS3-transactivated
Homo sapiens TGFB-induced factor 2
Homo sapiens protein phosphatase 1,
Homo sapiens ubiquitin-conjugating
Homo sapiens SMAD specific E3
Homo sapiens chromosome 9 open reading
Homo sapiens APG9 autophagy 9-like 1
Homo sapiens hypothetical protein
Homo sapiens ischemia/reperfusion
Homo sapiens solute carrier family 25
Homo sapiens PHD finger protein 17
Homo sapiens ubiquitin domain containing
Homo sapiens microtubule-associated
Homo sapiens chromosome 14 open
Homo sapiens hypothetical protein
Homo sapiens microtubule-actin
Homo sapiens F-box and WD-40 domain
Drosophila) (FBXW7), transcript
Homo sapiens beta-transducin repeat
Homo sapiens solute carrier family 35,
Homo sapiens chromosome 9 open reading
Homo sapiens BCL2-like 12 (proline rich)
Homo sapiens palmitoyl-protein
Homo sapiens transmembrane protein 55B
Homo sapiens chromosome 1 open reading
Homo sapiens solute carrier family 15,
Homo sapiens chromosome 9 open reading
Homo sapiens cytochrome b, ascorbate
Homo sapiens chromosome 10 open
Homo sapiens NOD9 protein (NOD9),
Homo sapiens U2AF homology motif
Homo sapiens family with sequence
Homo sapiens peroxisomal biogenesis
Homo sapiens exostoses (multiple) 2
Homo sapiens solute carrier family 11
Homo sapiens BCL2-interacting killer
Homo sapiens COX10 homolog,
Homo sapiens
Homo sapiens H2A histone family,
Homo sapiens solute carrier family 9
Homo sapiens signal recognition particle
Homo sapiens TAR (HIV) RNA binding
Homo sapiens kinesin family member 3B
Homo sapiens sperm associated antigen 7
Homo sapiens RCE1 homolog, prenyl
Homo sapiens heat shock 70 kDa
Homo sapiens chromosome 1 open
Homo sapiens coiled-coil domain
Homo sapiens nin one binding protein
Homo sapiens chromobox homolog 6
Homo sapiens KIAA0317 (KIAA0317),
Homo sapiens mitofusin 2 (MFN2),
Homo sapiens cytoplasmic linker
Homo sapiens anaphase promoting
Homo sapiens NMD3 homolog
Homo sapiens ubiquitin associated domain
Homo sapiens NCK interacting protein
Homo sapiens riboflavin kinase (RFK),
Homo sapiens chromosome 20 open
Homo sapiens vacuolar protein sorting
Homo sapiens hypothetical protein
Homo sapiens HCV NS3-transactivated
Homo sapiens chromosome 9 open
Homo sapiens APG9 autophagy 9-like 1
Homo sapiens hypothetical protein
Homo sapiens ischemia/reperfusion
Homo sapiens solute carrier family 25
Homo sapiens PHD finger protein 17
Homo sapiens ubiquitin domain containing
Homo sapiens solute carrier family 35,
Homo sapiens BCL2-like 12 (proline rich)
Homo sapiens transmembrane protein 55B
Homo sapiens chromosome 1 open
Homo sapiens solute carrier family 15,
Homo sapiens chromosome 9 open
Homo sapiens cytochrome b, ascorbate
Homo sapiens NOD9 protein (NOD9),
Homo sapiens U2AF homology motif
Homo sapiens family with sequence
Results: The set of miR-16 consensus downregulated transcripts (N=116) shown in TABLE 5 overlapped with, but were not identical to, computationally predicted miR-16 targets, as shown in
As shown in
With regard to cell cycle distribution of various cell types transfected with miR-16, it was observed that miR-16 induced the G0/G1 accumulation phenotype in HCT116 wt, DLD-1wt, and A549 cells (data not shown). However, the G0/G1 phenotype was not observed in Hela cells. While not wishing to be bound by theory, it is possible that miR-16 did not induce a G0/G1 accumulation in Hela cells because this cell type is known to be deficient in G1 checkpoint control.
This Example demonstrates that the miR-15a and miR-16 induced transcript downregulation induces cellular phenotypes.
Methods: 170 miRNA duplexes were transfected into HCT116dicerex5 and DLD-1Dicerex5 cells using the methods described in EXAMPLE 1. Control transfections were done in parallel with luciferase siRNA-transfected cells and with siRNA targeting PLK1, an essential mitotic kinase (Liu, X., et al., PNAS 100(10):5789-94 (2003)). Cell viability was measured 96 hours post transfection.
AUAGCCCUGUACAAUGCUGGUAU
AUAGCCCUGUACAAUGCUGGUAU
AGCUGGUUGAAGGGGAUUAAAU
AGUUUUGCAUGGAUUUGCAGAAU
AACCCAGCAGACAAUGUAGGUAU
AGACCGAGACAAGUGCAAAGAU
AACAUCAGUCUGAUAAGCAAAU
AUCUUUACCAGACAGUGAUAU
AAGUGCCUUCACUGCACUAU
ACCUGCACGAACAGCACUUUAU
AUCGUUACCAGACAGUGUAAAU
AGGCCGGGACAAGUGCAAUAAU
AGGCAGGCAUGCGGGCAGAUAAU
ACACACUUCCUUACAUUUUAAU
AGUUCUUCAACUGGCAGCAUAU
AGGCAUCAUAUAGGAGCUGCAAU
ACAAGUUAGGGUCUCAGGCAAU
AUCAGCUAACUACACUGCGUAU
AUCAGCUAAUGACACUGCGUAU
ACUCAAAAGAUGGCGGCUCAU
ACGCGUAUUCUUAAGCAAUAAU
AGCUACAGUGCUUCAUUUUAAU
ACGCGAGCCGAACGAACAAAAU
AACAACAGGAAACUACCAAAU
AGAAGCAUUUCCACACUCAU
AACAACAUGAAACUACCAAAU
AACCAGCUAAGACACUGCGAAU
AGCUCCUGGAGGACAGGCAAU
AAAGUUCUGUAGUGCACUCAAU
ACCAUAGGGUAAAACUAUUAU
AUGCAACUACAAUGCUCAU
AGUGAAGACACGGAGCCACAAU
AAAGUUCUGUGAUGCACUCAAU
AAGCUUUUUGCUCGUCUUUUAU
AGGUUCCACCCCAGCAGCCAU
AAACACCAUUGUCACACUCGAAU
ACAUAGGAAUAAAAAGCCAAAAU
AGCAGCACCUGGGGCAGUCGAU
AGUGGUCCUAAACAUUUCUCAU
ACGUGGAAUUUCCUCUAUCUAU
ACUUUUGUGACUAUGCUAAU
AAAAAUGUGCUAGUGCCAUAAU
AGUACAUCAUCUAUACUGAAAU
AUUAUUACUCACGGUACCAAU
ACCGAUUUCAAAUGGUGUUAAU
AGGAGUCUGAGCAUUUCAAU
AAGAGACCGGUUCACUGUCAAU
AAAUCCCUGGCAAUGUGUUAU
AUACAGCUAGAUAACCAAACAAU
AUGGUUAGAUCAAGCACUAAU
AGUUUUGCAUAGAUUUGCAGAAU
ACUGAAACAUGGAAGCAUUUAAU
ACAGGUAGUCUGAACACUGCGAU
AGGUCGACCGUGUAAUGUCCAU
AGUGUCAGAUACGGUGUCGAU
ACCAAAACAUGGAAGCACUUAAU
ACCAUGGCUGUAGACUGUAAAU
AAAUUCGGUUCUACAGGGAAAU
AACUUAGUAAUGUGCAAAAAU
AAGAGACCGGUUCACUGUCAAU
ACAGAUAGUCUAAACACUGCGAU
AAAGCGCCCCCAUUUUGACUAU
AACAAAAUCACUGAUGCUGCAAU
ACCAAUGCCCUAGGGGAUGGGAU
AGAUCAACCAUGUAUUAAUAU
AAGUUCUGUCAUGCACUCAAU
ACUAAAACAUGGAAGCACUAAAU
AAGCCCAGACCGCAAAAUGAU
AUUUUUCGUUAUUGCUCUUCAAU
AACGGAACCACUAGUGACUAGAU
ACCAUUGCUAAAGUGCAAUUAU
ACUGCCUUUCUCUUUAAU
AGGAUUGAGGGGGGGUUUUAU
ACUCAAACAUGGAAGCAUUUAAU
AGCAAGUACAUCCACGUUUAAU
AAUUGCGUUUGGACAAUGAAU
AAGGCAACACGGAUAAUUUAAU
ACCCCAAAAUCGAAGCACUUCAU
AGUGCUUUCUACUUUAAGAU
ACACGUGCACUGUACAAU
AAAGCACUUCCAUGUUAAACUAU
Results: TABLE 7 and
As shown in TABLE 7 and
It was also observed that miR-195 inhibited growth in DLD-1 Dicerex5 cells better than in HCT116 Dicerex5 cells (˜60% vs. ˜20%, respectively). In other experiments, growth inhibition by miR-195 was indistinguishable from that of miR-15a, miR-15b and miR-16 (data not shown). miR-103 and miR-107 were less effective at inhibiting cell growth in both cell lines tested. Of the other miRNAs that were analyzed by microarray profiling, shown in TABLE 7, miR-192, miR-215 and miR-133a gave stronger growth inhibitory phenotypes than miR-15a, miR-15b, and miR-16.
These results demonstrate that miR-15a, miR-15b, and miR-16 induce a measurable cellular phenotype consistent with their gene expression profiles. Because miR-15a, miR-15b, and miR-16 share seed region identity and induce identical expression profiles and cellular phenotypes, miR-16 was considered to be representative of the miR-16 family and was used for subsequent experiments.
This Example demonstrates that miR-16 regulates G0/G1 because cells transfected with miR-16 duplex had increased numbers of cells in G0/G1 and a corresponding decrease in S and G2/M.
Methods: HCT116Dicerex5 cells were transfected with miRNA duplexes for miR-106b (SEQ ID NO:18, 19), miR-16 (SEQ ID NO:30, 31), kshv-miR-K12-6-5p (SEQ ID NO:64, 65), kshv-mir-K12-3 (SEQ ID NO:66, 67), and miR-16 (SEQ ID NO:30, 31) containing mismatches at positions 2 and 3 (2, 3 mm) or 18 and 19 (18, 19 mm). Transfections were carried out as described in EXAMPLE 1.
24 hours after transfection, floating and adherent cells were harvested, combined, and processed. Alternatively, Nocodazole (100 ng/ml, Sigma-Aldrich) was added beginning at 20 hours post transfection, and cells were analyzed for cell cycle distribution 18 hours later (46 hours post transfection). Cells were collected by centrifugation, fixed with ice-cold 70% ethanol, washed with PBS, and resuspended in PBS containing Proridium Iodide (10 μg/ml) and RNAse A (1 mg/ml). After a final incubation at 37° C. for 30 minutes, cells were analyzed by flow cytometry using a FACSCalibur flow cytometer (Becton Dickinson). A total of 10,000 events were counted for each sample. Data was analyzed using FlowJo software (TreeStar, Ashland, Oreg.).
Results:
The fluorescence intensity of cells having diploid DNA content (2N), and cells having tetraploid DNA (4N) content was analyzed by FACS. The results are shown below in TABLE 8.
As shown in TABLE 8, miR-106b transfected cells gave a cell cycle profile indistinguishable from mock or luciferase-transfected cells. miR-16 and viral miRNA kshv-miR-K12-6-5p which expresses a microRNA from the miR-16 family caused accumulation of cells at G0/G1 (2N DNA content). kshv-miR-6-5p therefore mimics the activity of a normal cellular microRNA family. In contrast, another KSHV miRNA, kshv-miR-K12-3, which does not express an miR-16 microRNA, had no observable effect on cell cycle distribution. It has been shown that kshv-miR-K12-6-5p is found in latently infected cells, so it is possible that miR-16 expression disrupts cell cycle control in Kaposi's sarcoma-associated herpes virus infected cells, which may play a role in the establishment or maintenance of latency (see Cai et al., PNAS 102:5570-5575 (2005); Verma, S. C., et al., FEMS Microbiol Lett 222:155-163 (2003)).
The G0/G1 accumulation phenotype observed in miR-16 transfected cells was reversed by paired seed region mismatches in miR-16 at positions 2+3 (and 4+5, not shown), but not by matches outside the seed region at positions 18+19 (and 19+20, not shown). It was also determined that miR-16 containing seed region mismatches at 2+3 and 4+5 triggered different expression profiles than miR-16, and the down-regulated transcripts were found to be enriched for hexamers matching the mismatched seed regions (data not shown). Therefore, cell cycle and gene expression phenotypes induced by miR-16 are miRNA seed sequence-dependent.
The G0/G1 accumulation phenotype became clearer when the microtubule de-polymerizing drug, Nocodazole, was added before analysis to block cells from re-entering the cell cycle after mitosis. As shown in TABLE 8, Nocadazole treatment cause nearly all miR-106b transfected cells to accumulate in G2/M (4N DNA content), whereas a large fraction of miR-16-transfected cells remained in G0/G1. The extent of G0/G1 accumulation varied with the concentration of miR-16 transfected into cells, with clear effects seen at duplex concentrations of <1 nM (data not shown).
This Example demonstrates that transfection of siRNA pools directed to miR-16 consensus downregulated transcripts into HCT116 cells triggered G0/G1 accumulation of the transfected cells.
Methods:
Preparation of siRNA Pools
siRNA sequences were designed with an algorithm developed to increase efficiency of the siRNAs for silencing while minimizing their off-target effects, as described in Jackson et al., Nat Biotech 21:635-637 (2003), incorporated herein by reference. siRNA duplexes were ordered from Sigma-Proligo (Boulder, Colo.). The siRNA sequences designed to target miR-16 consensus downregulated transcripts are provided below in TABLE 8.
Transfections: HCT116Dicerex5 cells were transfected with siRNA pools (with a separate transfection for each pool) targeting 102 identified downregulated consensus transcripts containing matches to the miR-16 seed region (miR-16 targets) and 51 control transcripts that did not contain miR-16 seed region matches (non-miR-16 targets). For the siRNA pools, three siRNAs targeting the same transcript were pooled at equal molarity (final concentration of each siRNA (17 nM); total siRNA concentration (50 nM)). Transfections were carried out as described in EXAMPLE 1. miR-16 was transfected as a positive control. The transfected cells were then analyzed for cell cycle distribution. To control for variation in transfection efficiency between different experiments, the percentages of cells in G0/G1 were normalized so that background control (the average of mock and luciferase-transfected cells) and miR-16 transfected (100 nM) cells gave 0% and 100% cells in G0/G1, respectively ((% G1sample-% G1background)/(% G1miR-16-% G1background)*100). Before normalization, background values averaged 6.4+1.1 (mean±SD) % cells in G0/G1, and miR-16-transfected cell values averaged 41+5% cells in G0/G1 (based on seven independent experiments).
The non-miR-16 targets tested were: EEA, FAM57A, GBP3, JUNB, L3 MBTL3, MAST2, NAGK, NBL, NUBP, PLAU, TNKS2, ADAM17, DUSP11, KHSRP, NDST2, PKP4, SQLE, SRP54, STX7, TAF2, TFR2, THOP1, TRAF3, TSPAN4, TSPAN6, TTC3, TTF2, USP1, UVRAG, YARS, RUVBL1, B4GALT3, MBTPS1, TNFSF10, ADAM9, MTMR1, PEX11B, SUCLA2, DPM2, PCAF, RAB7L1, OSMR, ARRB1, CRY1, DCTN1, STOM, FEN1, GNG11, SLC1A33.
Results: The HCT116Dicerex5 transfection results of siRNA pools directed against miR-16 consensus transcripts are shown above in TABLE 9.
As shown in
This Example demonstrates that pools of siRNAs transfected into cells are capable of inducing the same phenotype as miR-16 transfected into cells.
Methods:
siRNA pools were selected for further analysis:
CARD10-1 (104.3% G1 in TABLE 9), CDC27-1 (29.7%), CDK6-1 (24.8%), C10orf46 (33.1%). MiR-16, miR-106(b) control duplex and siRNA pools targeting CARD10, CDC27-1, CDK6-1 and C10orf46, were each transfected into HCT116Dicerex5 cells (3 siRNAs per pool, each at 33 nM; 100 nM total concentration). Cell cycle analysis was carried out using FACs as described above in EXAMPLE 4.
Results of siRNA transfected pools: As shown in
Transfection of Individual siRNAs:
siRNAs to selected miR-16 downregulated targets that induced G0/G1 accumulation, shown below in TABLE 10, were individually transfected into HCT116Dicerex5 cells at 0.25 nM. The percent of cells in G1 was determined as described in EXAMPLE 4.
mRNA silencing was quantified at 48 hours by real-time PCR using an ABI PRISM 7900HT Sequence Detection System and Assays-on-Demand gene expression products (Applied Biosystems, Foster City, Calif.). The mRNA value for each gene was normalized relative to GUSB (No. 431088E) mRNA levels in each RNA sample.
Protein products of targeted transcripts CDK6 and CDC27 were also measured in siRNA transfected cells. Immunoblotting was performed as described in Jackson et al., RNA 12(7):1-8 (2006). Anti-CDK6 monoclonal antibody (DCS-90) was purchased from Abcam (Cambridge, Mass.).
A subset of miR-16 targets CDK6, CARD10, CDC27 and C10orf46 were chosen for analysis. For each of these targets, at least 2 siRNAs were identified that induced G0/G1 accumulation in greater than 20% of the cells (CDK6 (2/3); CARD10 (2/3); CDC27(2/6); and C10orf46 (2/6) individual siRNAs tested). mRNA silencing data for these siRNAs was also determined, as shown in TABLE 10. In all cases, it was determined that siRNAs inducing the strongest cell cycle phenotypes also silenced the mRNA of their targets the most, as determined by the % mRNA remaining.
Results: The 25 siRNA pools that gave >20% G0/G1 accumulation were deconvoluted and individual siRNAs (100 nM) were tested for their ability to phenocopy miR-16. miR-16 (1 nM and 100 nM) and miR-106b (100 nM) were transfected as controls. For 24 of the 25 pools, at least one individual member of the siRNA pool gave a phenotype as strong as the pool value. The results of the individually transfected siRNAs from each selected pool are shown above in TABLE 10 (Exp 1).
Individual siRNA molecules were then titrated to a concentration where most gave phenotypes only slightly above background (0.25 nM, compared to 100 nM total concentration used in
An experiment was carried out in which individual siRNAs (0.25 nM) for CARD10-4, KIAA0317-3 and C9orf91-1, the three siRNAs that gave the strongest phenotype shown in
A subset of targets, including CDK6, CARD10, CDC27, and C10orf46, were chosen for further analysis. An experiment was conducted to demonstrate that the cell cycle phenotypes triggered by siRNAs to these targets resulted from silencing of these targets, and not due to indirect effects from unintended targets (see. e.g., Jackson et al., Nat Biotech 21:635-637 (2003)). As shown in TABLE 10, the siRNAs that triggered the strongest cell cycle phenotypes also generally were found to silence their targets the most, as determined by % mRNA remaining.
Another experiment was carried out in which pooled siRNAs to two miR-16 targets (CDK6 and CDC27, each at 0.25 nM) or three miR-16 targets (CDK6, CDC27 and CARD10, each at 0.25 nM) were transfected into HCT116Dicerex5 cells. In each case, the pools gave a stronger phenotype than the individual siRNAs (at 0.25 nM) (data not shown).
These experiments suggested that the robust G0/G1 cell cycle phenotype elicited by miR-16 mediated regulation results from coordinated silencing of multiple miR-16 responsive targets.
This Example demonstrates the synergistic effect of transfecting a pool of siRNAs directed to different miR-16 responsive targets in comparison to transfection of the corresponding individual siRNAs.
Methods:
Transfection of Pooled siRNAs to Different miR-16 Responsive Targets
Individual siRNAs targeting CDK6 (siRNA=CDK6-3), CDC27 (siRNA=CDC27-3), CARD10 (siRNA=CARD10-4) and C10orf46 (siRNA=C10orf46-3) (shown in TABLE 11) were transfected into HCT116Dicerex5 cells at a concentration of 0.25 nM. These four siRNAs (CDK6-3, CDC27-3, CARD10-4 and C10orf46-3) were also transfected together as a pool (0.25 nM each siRNA, total concentration=1 nM). miR-16 and miR-106b were transfected as controls (1 nM each). Cell cycle phenotypes were determined by FACS analysis as described in EXAMPLE 4.
Results: As shown in
This example describes the cellular phenotype induced in cells transfected with miR-16 or miR-106b and in cells transfected with anti-miR-16 or anti-miR-106b.
Methods:
Cell Cycle Analysis
HCT116Dicerex5 cells were transfected with increasing concentrations of miR-16 duplex (SEQ ID NO:3) (0.5, 1, 10 and 100 nM), or plasmids carrying miR-16 expressed as an shRNA (miR-16 hairpin, Silva, et al., Nat Genet 37:1281-1288 (2005)), or from its endogenous locus on chromosome 13 (miR-16 locus). In both cases, expression of the precursor form of miR-16 was under the control of an H1 promoter. For transfections with plasmids, the DNA concentration was 1.5 μg DNA per 2×105 to 3×105 cells in a 6-well dish.
The transfected cells were treated with nocodazole and analyzed as described in EXAMPLE 4. The percentages of cells in G0/G1 in different experiments were normalized so that background control and miR-16-transfected cells gave 0% and 100% of cells in G0/G1, respectively [(the percentage of G1 cells in the sample—the percentage of background G1 cells)/(the percentage of G1 miR-16-transfected cells—the percentage of background G1 cells)×100].
For duplex transfections, mock-transfected cells were used to determine the percentage of background G1 cells. For plasmid transfections, cells transfected with an empty vector were used to determine the percentage of background G1 cells. miR-16 copy numbers (copies/20 pg of RNA) were determined by a quantitative primer extension PCR assay (Raymond, C. K., et al., RNA 11:1737-1744 (2005)). This assay preferentially detects mature miR-16 over duplex and hairpin forms.
Results:
Transfection with empty vector and an shRNA for miR-106b did not result in the accumulation of cells in G0/G1 or miR-16 levels greater than those in mock-transfected cells (data not shown).
As shown in
Inhibition of miR-16
Rationale: The following experiments were carried out in order to determine if the cellular phenotype induced by miR-16 over-expression experiments was reversed in loss-of-function experiments.
Methods:
HeLa cells, TOV21G, HCT116wild-type and HCT116Dicerex5 cells were transfected with luciferase siRNA, miR-16, miR-106b, and anti-miR-16 or anti-miR106b (LNA or 2′-O-methyl modified oligonucleotide inhibitors of miR-16 and miR-106b, Hutvagner, G. M., et al., PloS Biol 2:E98 (2004); Orom, U. A., et al., Gene 372:137-141 (2006)). Gene expression analysis was carried out on the transfected cells using microarrays as described in EXAMPLE 2.
Results:
TABLE 11 and TABLE 12 summarize the gene expression results from transfected cells with anti-miR-16 or anti-miR-106b, respectively. For controls, patterns of regulation of randomized sets of genes were compared from luciferase siRNA-treated cells. The median change or median percentage of increase in regulation was compared. Control sets gave 0.01%+0.6% change. The # of up-regulated, percentage of miR-16 or miR-106b consensus targets having a level of regulation >0. Control sets gave 50%±5% of targets with a level of regulation of >0. For the up-regulation P value, Wilcoxon signed-rank P values for the up-regulation of the indicated target sets. Mitotic cell cycle genes, transcripts down-regulated by the miR-16 duplex at 24 h and annotated with the GO biological process term “mitotic cell cycle”.
As shown above in TABLES 11 and 12, the median regulation of consensus miR-16 and miR-106b targets was slightly increased in cells transfected with the specific anti-miR. When considered as a group, miR-16 targets were more significantly up-regulated than miR-106b targets in cells treated with anti-miR-16 (see TABLE 11). Target regulation was greatly reduced in HCT116Dicerex5 cells, as expected given the reduced miRNA levels in these cells. Likewise, anti-miR-106b-treated cells showed significant regulation of miR-106b, but not miR-16 consensus targets (see TABLE 12). Taken together, these results indicate that in cells with wild-type Dicer function, many, if not most, miR-16 and miR-106b targets down-regulated by miRNA duplexes in gain-of-function experiments (as shown in
As further shown in TABLE 11, miR-16-regulated transcripts annotated with the GO biological process term “mitotic cell cycle” (see TABLE 4) were significantly up-regulated by anti-miR-16 in TOV21G and HCT116 wild-type cells. These transcripts were less regulated in HeLa cells, which do not show an miR-16 cell cycle phenotype, and in HCT116Dicerex5 cells, which have reduced endogenous levels of miR-16. As shown in TABLE 12, mitotic cell cycle transcripts were not significantly regulated in any cell line by anti-miR-106b. These results support the regulation of cell cycle progression by physiological levels of miR-16.
Analysis of the expression of cell cycle genes of anti-miR-16-treated TOV21G and HCT116 wild-type cells did not reveal any obvious differences from that of control-treated cells (data not shown). Thus, gene expression changes measured by microarray were not sufficient to identify a detectable cell cycle phenotype in these cells.
While not wishing to be bound by theory, it was hypothesized that cells having higher endogenous levels of miR-16 would be more susceptible to anti-miR-16-induced phenotypic changes detectable by flow cytometry. In order to determine endogenous miR-16 levels, a number of transfectable cell lines were screened using a quantitative primer extension PCR assay (Raymond, C. K. et al., RNA 11:1737-1744 (2005)). These screening experiments showed that the cell line SW1417 had elevated levels of endogenous miR-16 (˜1,500 copies/cell) (data not shown). The transfection of SW1417 cells with anti-miR-16 resulted in a significant decrease in numbers of G0G1 compared to mock-transfected cells [10%±1% decrease (three independent experiments); P<1E-3. In contrast, anti-miR-106b-transfected cells did not show significant differences from mock-transfected cells [5%±3% decrease (three independent experiments); P>5E-2]. Therefore, these results indicate that disruption of physiological miR-16 levels in certain cell types can alter cell cycle distribution.
This example demonstrates that the levels of miR-16 family-down-regulated transcripts negatively correlate with miR-195 levels in human tumors.
Rationale: An experiment was carried out to determine whether gene expression changes triggered by miRNA transfection in the in vitro models described in EXAMPLES 1-7 reflect the relationship between steady state levels of transcripts down-regulated by miR-16 family members and miRNA levels in human tumors.
Methods:
RNA was isolated from a series of 29 tumors and 28 adjacent uninvolved normal tissues. mRNA expression was measured using microarrays, and miR-195 levels were determined using a quantitative primer extension PCR assay (Raymond, C. K. et al., RNA 11:1737-1744 (2005)). miR-195 (an miR-16 family member) was chosen for this due to the availability of reliable tumor atlas expression data.
mRNA and miRNA expression levels in human tumors (including breast cancer, lung cancer, colon cancer, kidney cancer and gastric cancer), and adjacent normal tissues were expressed as ratios of these levels to expression levels in a pool of normal samples from each tissue type. Correlations between expression level ratios for miR-195 and transcripts down-regulated 24 h after the transfection of tissue culture cells with miR-16 were calculated. As a control, correlations were also calculated for ˜200 random permutations of expression ratios (random transcripts).
Results:
While illustrative embodiments have been illustrated and described, it will be appreciated that various changes can be made therein without departing from the spirit and scope of the invention.
| Filing Document | Filing Date | Country | Kind | 371c Date |
|---|---|---|---|---|
| PCT/US07/71239 | 6/14/2007 | WO | 00 | 9/8/2009 |
| Number | Date | Country | |
|---|---|---|---|
| 60813827 | Jun 2006 | US |
