Nothing Special   »   [go: up one dir, main page]

Skip to main content
RNA logoLink to RNA
. 2008 Dec;14(12):2580–2596. doi: 10.1261/rna.1351608

Molecular characterization of human Argonaute-containing ribonucleoprotein complexes and their bound target mRNAs

Markus Landthaler 1, Dimos Gaidatzis 2,3, Andrea Rothballer 1, Po Yu Chen 1, Steven Joseph Soll 1, Lana Dinic 1, Tolulope Ojo 1, Markus Hafner 1, Mihaela Zavolan 2, Thomas Tuschl 1
PMCID: PMC2590962  PMID: 18978028

Abstract

microRNAs (miRNAs) regulate the expression of mRNAs in animals and plants through miRNA-containing ribonucleoprotein particles (RNPs). At the core of these miRNA silencing effector complexes are the Argonaute (AGO) proteins that bind miRNAs and mediate target mRNA recognition. We generated HEK293 cell lines stably expressing epitope-tagged human AGO proteins and other RNA silencing-related proteins and used these cells to purify miRNA-containing RNPs. Mass spectrometric analyses of the proteins associated with different AGO proteins revealed a common set of helicases and mRNA-binding proteins, among them the three trinucleotide repeat containing proteins 6 (TNRC6A,-B,-C). mRNA microarray analyses of these miRNA-associated RNPs revealed that AGO and TNRC6 proteins bind highly similar sets of transcripts enriched in binding sites for highly expressed endogenous miRNAs, indicating that the TNRC6 proteins are a component of the mRNA-targeting miRNA silencing complex. Together with the very similar proteomic composition of each AGO complex, this result suggests substantial functional redundancy within families of human AGO and TNRC6 proteins. Our results further demonstrate that we have developed an effective biochemical approach to identify physiologically relevant human miRNA targets.

Keywords: microRNA, Argonaute, GW182/TNRC6, microRNA targets, immunoprecipitation

INTRODUCTION

miRNAs comprise a large family of evolutionary conserved small regulatory RNAs that have important regulatory functions in a wide array of biological processes, including developmental timing, cell differentiation, cell proliferation, apoptosis, and patterning of the nervous system (Ambros 2004; Jones-Rhoades et al. 2006; Bushati and Cohen 2007).

miRNA biogenesis has been reviewed in detail (Lee et al. 2006a; Chapman and Carrington 2007). Briefly, miRNAs are processed from long primary RNA polymerase II transcripts (pri-miRNAs) by the RNase III Drosha and the associated double-stranded RNA-binding protein DGCR8 into hairpin pre-miRNAs (Denli et al. 2004; Gregory et al. 2004; Han et al. 2004, 2006; Landthaler et al. 2004). The pre-miRNAs are exported to the cytoplasm by exportin-5 (Yi et al. 2003; Lund et al. 2004) and are further processed by the RNase III nuclease DICER (Bernstein et al. 2001; Grishok et al. 2001; Hutvágner et al. 2001; Ketting et al. 2001). Typically, only one strand of the miRNA duplex is incorporated into the miRNA silencing effector protein (miRNPs) complexes (Khvorova et al. 2003; Schwarz et al. 2003). The mature miRNA then guides miRNPs to messenger RNAs (mRNAs).

The recognition of mRNAs by miRNPs results in mRNA cleavage of targeted transcripts that are perfectly complementary to the miRNA (Hutvágner and Zamore 2002; Yekta et al. 2004; Barth et al. 2008) and inhibition of translation in case of limited base-pairing interactions between the 5′-end of the miRNA and complementary sequences in the 3′-untranslated regions (UTRs) of the target mRNAs (Wightman et al. 1993). Experimental evidence suggests that the inhibition of protein synthesis mediated by miRNPs takes place at the stage of translational initiation (Humphreys et al. 2005; Pillai et al. 2005; Thermann and Hentze 2007), as well as during the later elongation step (Olsen and Ambros 1999; Seggerson et al. 2002; Zeng et al. 2002; Petersen et al. 2006). In addition, miRNAs can also induce degradation of mRNA targets despite imperfect miRNA–mRNA base pairing (Bagga et al. 2005; Jing et al. 2005; Krutzfeldt et al. 2005; Lim et al. 2005; Rehwinkel et al. 2005; Giraldez et al. 2006; Schmitter et al. 2006; Wu et al. 2006; Grimson et al. 2007; Linsley et al. 2007; Rodriguez et al. 2007).

In animal cells, miRNAs recognize their target mRNAs through base-pairing interactions between 6 and 8 nucleotides (nt) at the 5′ end of the miRNA, the so-called “seed,” and complementary nucleotides in the 3′ UTR of mRNAs. Computational methods to predict miRNA binding sites in several organisms mostly rely on the interaction of the miRNA seed region with its target mRNA 3′ UTR and evolutionary conservation of the target site (John et al. 2004; Grun et al. 2005; Krek et al. 2005; Lewis et al. 2005; Gaidatzis et al. 2007; Grimson et al. 2007). Different experimental approaches have been used to identify miRNA targets and to study the physiological interactions of miRNAs and target mRNAs (Krutzfeldt et al. 2005; Lim et al. 2005; Rehwinkel et al. 2006; Schmitter et al. 2006; Vinther et al. 2006; Beitzinger et al. 2007; Easow et al. 2007; Karginov et al. 2007; Orom and Lund 2007; Zhang et al. 2007; Zhao et al. 2007; Baek et al. 2008; Hendrickson et al. 2008; Selbach et al. 2008).

Biochemical and genetic studies have shown that Argonaute (AGO) proteins are the core component of miRNP silencing complexes (Hammond et al. 2001; Martinez et al. 2002; Mourelatos et al. 2002; for reviews, see Meister and Tuschl 2004; Patel et al. 2006; Tolia and Joshua-Tor 2007). AGO proteins share a conserved PAZ and PIWI domain, and both domains are involved in the interaction with small RNAs (Song et al. 2003; Yan et al. 2003; Lingel et al. 2004; Ma et al. 2004). Structural studies further revealed that the PIWI domain adopts an RNase H-like fold, which forms the catalytic center of the AGO endonuclease in the RNA-induced silencing complex (Song et al. 2004; Ma et al. 2005; Parker et al. 2005).

A number of proteins were found to be associated with AGO proteins and the miRNP complexes (for reviewed, see Peters and Meister 2007). More recently, biochemical purifications of human EIF2C1/AGO1 and EIF2C2/AGO2 revealed additional factors, including helicases, heterogeneous nuclear ribonucleoproteins, mRNA-binding proteins and proteins involved in RNA metabolism (Hock et al. 2007). Several of these proteins function in the biogenesis and activity of miRNAs. DICER and the double-stranded RNA binding proteins, TARBP2 and PRKRA/PACT, are involved in processing of the miRNA precursor (Chendrimada et al. 2005; Haase et al. 2005; Lee et al. 2006b). RNA helicase A (DHX9) interacts with RNA-induced silencing complex (RISC) in human cells and functions in RISC loading (Robb and Rana 2007). TNRC6A/GW182 and TNRC6B, markers of mRNA-degrading cytoplasmic processing bodies (P- or GW-bodies), are required for RNA interference in human cells (Jakymiw et al. 2005; Liu et al. 2005; Meister et al. 2005). The TNRC6 proteins seem to function downstream of Argonaute in the process of mRNA degradation with the final steps likely to occur in P-bodies (Behm-Ansmant et al. 2006). Another P-body component, the helicase RCK/DDX6, is also required for translation regulation by miRNAs highlighting the importance of the cytoplasmic P-body structures in RNA silencing (Chu and Rana 2006).

We performed a comprehensive biochemical analysis to understand the molecular composition of the human RNA silencing machinery by proteomic and ribonomic profiling of the human AGO RNPs. We observed a common set of protein interactors for EIF2C1,-2,-3,-4, including DICER, TNRC6A,-B,-C (or GW182 proteins), heat-shock proteins, helicases, and various mRNA-binding proteins. The mRNA profiles of immunoprecipitates of the four AGO proteins and their associated TRNC6 proteins were strikingly similar and characterized by the enrichment of transcripts with seed sequence complementarity to expressed miRNAs. The high similarity of AGO- and TNRC6-immunopurified mRNAs suggests a substantial degree of redundancy in the regulatory function within the AGO and TNRC6 protein families, and further provides evidence that these proteins are components of the miRNA-containing mRNA-targeting complexes.

RESULTS

Proteomic analysis of human argonaute protein interactors

To characterize the molecular composition of human EIF2C1/AGO1 through EIF2C4/AGO4 (collectively referred to as AGO proteins) containing protein complexes, we generated HEK293 cells stably expressing FLAG/HA-tagged AGO proteins. FLAG/HA-AGO protein complexes were affinity purified from HEK293 total cell extracts by incubation with anti-FLAG agarose beads. The immunoprecipitated FLAG–protein complexes were eluted using FLAG peptide, followed by a second immunoprecipitation (IP) using anti-HA magnetic beads. The immunoprecipitated proteins were then separated by SDS-PAGE (Supplemental Fig. 1), excised, and identified by mass spectrometry (MS). Table 1 lists the names of the identified proteins in FLAG/HA-EIF2C1,-2,-3,-4 IPs. Surprisingly, all four AGO proteins coimmunoprecipitated a highly similar set of proteins. We identified the previously described interactors DICER, MOV10, TNRC6A, and TNRC6B (as well as the paralogous TNRC6C as a new interactor). Reciprocal purifications using tagged DICER, TNRC6A, and TNRC6B-expressing cells lines confirm these interactions (Table 1). In addition IPs of FLAG/HA-DICER coprecipitated among other proteins the know interactors TARBP2 and PRKRA/PACT (Chendrimada et al. 2005; Haase et al. 2005; Lee et al. 2006b).

TABLE 1.

Summary of proteins identified by double-affinity purifications

graphic file with name 2580tbl1a.jpg

graphic file with name 2580tbl1b.jpg

Furthermore, we encountered a number of additional proteins in AGO IPs. These copurified proteins fell within several functional groups. The largest group included the mRNA-binding proteins, in particular the heterologous nuclear ribonucleoproteins: HNRNPAB, HNRNPC, HNRNPD, HNRNPF, HNRNPH1, HNRNPK, HNRNPL, HNRNPM, SYNCRIP/HNRPQ, HNRNPR, and HNRNPU. In addition, we detected several mRNA-binding proteins with functions in mRNA transport, stabilization and translation, most prominently IGF2BP1–IGF2BP3, YBX1, ELAVL1, FUBP3, and RALY. Moreover, the polyadenylate binding proteins 1 and 4 (PABPC1 and PABPC4) were present in the IPs. In addition the RNA helicases DDX5, DDX17, and DHX9, but not DDX6, which is required for miRNA function (Chu and Rana 2006), were identified as well as a series of 40S and 60S ribosomal proteins.

The interaction of selected proteins with the AGO proteins was confirmed by co-IP experiments (Fig. 1). Cells stably expressing FLAG/HA-EIF2C1 were transfected with plasmids encoding the myc-tagged version of several proteins that were identified by MS. Extracts from transfected cells were immunoprecipitated with anti-FLAG agarose and FLAG eluted. The eluted proteins were separated by SDS-PAGE, blotted, and probed with an anti-myc antibody. The myc-tagged proteins EIF2C1, EIF2C2, DICER, TNRC6C, MOV10, DDX17, PABPC4, YBX1, IGF2BP3, and ELAVL1 could be coimmunoprecipitated, whereas the mRNA-binding protein QKI, which was not identified in the MS analysis, but is expressed in HEK293 cells (data not shown), could not be detected. Treatment of the lysate with RNase T1 prior to the IP abolished the interaction with most proteins except DICER and TNRC6C, suggesting that the interactions between FLAG/HA-EIF2C1 and the other myc-tagged proteins were RNA mediated.

FIGURE 1.

FIGURE 1.

Proteins identified by mass spectrometry interact with EIF2C1. HEK293 cells stably expressing FLAG/HA-EIF2C1 were transfected with constructs encoding myc-tagged variants of the proteins identified by MS. Transfected cells were lysed and cytoplasmic extracts prepared in the absence or presence (+RNase T1) of ribonuclease T1. FLAG/HA-EIF2C1 was immunoprecipitated using anti-FLAG beads and eluted from the beads using FLAG peptide. Cell extracts and FLAG-eluted proteins (FLAG-IP) were separated by SDS-PAGE, blotted, and probed with anti-HA and anti-myc antibodies.

AGO proteins associate with ribosomes

miRNAs and miRNA-repressed mRNAs, as well as AGO proteins, have been described to be associated with polysomes (Olsen and Ambros 1999; Seggerson et al. 2002; Kim et al. 2004; Nelson et al. 2004; Maroney et al. 2006; Nottrott et al. 2006; Balzer and Moss 2007; Vasudevan and Steitz 2007). We re-examined the distribution of EIF2C2, exemplary for the AGO family, using polysome gradients. Cytoplasmic extracts of FLAG/HA-tagged cells were prepared and fractionated on sucrose gradients (Fig. 2). The absorbance profiles at 254 nm showed a pattern of ribosomes of exponentially growing cells with ribosomal subunits, monosomes, and polysomes. The protein fractions of the gradient were separated by SDS-PAGE, blotted, and probed with an antibody against ribosomal protein L7 (RPL7), a component of the 80S monosome and the 60S large subunit, to identify the ribosome-containing fractions. Consistent with previous findings (Maroney et al. 2006; Nottrott et al. 2006; Balzer and Moss 2007), probing the blot with an anti-HA antibody indicated the presence of FLAG/HA-EIF2C2 throughout the gradient, including the the polysomal fractions. In polysome gradients of a cellular extract that was prepared in the presence of 30 mM EDTA, to dissociate the polysome into the ribosomal subunits and monosomes, FLAG/HA-EIF2C2 shifted from the denser fractions to the lighter fractions of the gradient, suggesting that EIF2C2 cofractionates with polysomes as observed previously for miRNAs (Maroney et al. 2006; Nottrott et al. 2006).

FIGURE 2.

FIGURE 2.

EIF2C2 is associated with polysomes in HEK293 cells. Cytoplasmic extracts were prepared by cell lysis with buffer containing 5 mM MgCl2 or 30 mM EDTA, centrifuged through sucrose gradient, and fractionated. Proteins of each fraction were TCA-precipitated and equal aliquots were separated by SDS-PAGE, blotted, and probed with antibodies as indicated. Absorbance of fractions at 254 nm is shown.

Next, we immunoprecipitated FLAG/HA-EIF2C2 with anti-FLAG beads from cytoplasmic extracts prepared with 2 mM EDTA to partially dissociate the polysomes. The precipitated protein complexes were eluted with FLAG peptide and loaded on a 10%–45% sucrose gradient and fractionated by centrifugation. Probing for FLAG/HA-EIF2C2 with an anti-HA antibody indicated that the epitope-tagged AGO protein was present in at least three major complexes, which we designated complexes I, II, and III (Fig. 3A), similar to observations made by Hock et al. (2007). The coimmunoprecipitated ribonuclease DICER sedimented with complex I, likely representing the RISC loading complex previously described (Pellino et al. 2005). Interestingly, in complex III we were able to detect the small and large ribosomal subunit proteins, RPS6 and RPL7, respectively. The presence of both of these proteins suggests that FLAG/HA-EIF2C2 is associated with the 80S ribosome, whereas complex II constitutes the small ribosomal subunit. Northern hybridization using a probe for miR-16, a highly abundant miRNA in HEK293 cells, showed the association of this miRNA with each of the three complexes (Fig. 3B).

FIGURE 3.

FIGURE 3.

EIF2C2 is associated with ribosomes. (A) FLAG/HA-EIF2C2 was immunoprecipitated from a cytoplasmic extract prepared with lysis buffer containing 2 mM EDTA. The epitope-tagged protein was eluted using FLAG peptide and separated on a 10%–45% sucrose gradient for 8.5 h. Proteins of each fraction were TCA precipitated and separated by SDS-PAGE (silver-stained gel is shown on top). Western analysis of the separated fractions was performed using antibodies as indicated. (B) miR-16 is associated with three FLAG/HA-EIF2C2-containing complexes. FLAG/HA-EIF2C2 was purified and displayed on a sucrose gradient as described in A. The distribution of FLAG/HA-EIF2C2 in the gradient was examined by Western analysis using an anti-HA antibody and compared with the presence of RPS6 and RPL7 with antibodies against these ribosomal proteins. RNA was extracted from the material loaded on to the gradient (L), and each fraction of the gradient, pooled as indicated, separated by 15% denaturing PAGE, blotted, and probed for miRNA miR-16. A total of 20 μg of HEK293 cellular RNA (HEK) was loaded along with 10 and 4 fmol of synthetic miR-16 as control.

miRNA target identification by AGO immunoprecipitation

To confirm that the characterized complexes indeed represent miRNPs involved in silencing, we developed a method to characterize the mRNAs present in these complexes. To identify the mRNAs associated with FLAG/HA-tagged AGO protein complexes, we isolated total RNA from the FLAG-immunoprecipitates (IPs) and performed an mRNA microarray analysis. IPs of AGO-containing protein complexes was confirmed by Western analysis (Supplemental Fig. 2A) and by Northern hybridization for miR-16 (Supplemental Fig. 2B).

mRNA profiles of the cleared lysates prior to the IPs were highly similar with correlation coefficient values of 0.94–0.99. The high similarities of the expression profiles between cells expressing various tagged proteins indicate that the expression of the tagged AGO proteins had little effect on mRNA expression profiles. The normalized mRNA profiles of IPs and cleared lysates were compared to determine the degree of enrichment of each individual transcript in the IP versus the lysate. The distribution of the fold changes of the enrichment and underrepresention of transcripts that were immunopurified in different experiments are summarized in Supplemental Figure 3.

For a subset of experiments, we also compared the immunoprecipitated mRNAs with the mRNA profile of the supernatants after the IP. Since contrasts generated by comparing the profile of the IP mRNA versus lysate mRNA or supernatant mRNA were qualitatively similar, only the IP/lysate contrasts were used for further analyses.

To determine whether or not different AGO proteins bind to the same mRNA transcripts, we selected the top 2.5% immunoprecipitated transcripts in each individual experiment indicated in Figure 4A, and used the union of these transcripts to perform two-dimensional hierarchical clustering. We found a striking similarity between the mRNAs purified by anti-FLAG-agarose IPs of epitope-tagged EIF2C1,-2,-3,-4. This set of enriched transcripts was a result of specific binding of mRNAs directly or indirectly to the immunoprecipitated proteins, since it was obtained independently of the antibody or bead matrix used in these experiments. IPs of the FLAG/HA-tagged EIF2C2 expressing cells using an anti-FLAG, an anti-HA, and a monoclonal antibody against EIF2C2 bound by Protein G magnetic beads, as well as an IP of FLAG/HA-QKI-expressing cells using the anti-EIF2C2 antibody, immunopurified a set of transcripts highly similar to that obtained from FLAG-agarose IPs. The high similarity of the mRNAs bound by EIF2C1,-2,-3,-4 is mirrored by the striking resemblance in the proteomic composition, suggesting that AGO proteins likely have redundant regulatory functions.

FIGURE 4.

FIGURE 4.

AGO and TNRC6 proteins immunopurify a similar set of transcripts enriched for sequences complementary to miRNA seeds (A) Two-dimensonal hierarchical clustering for the enriched transcripts in IPs. IPs were performed from extracts of cells stably expressing the FLAG/HA-tagged protein as indicated using anti-FLAG agarose beads (FA) or Protein G Dynabeads bound anti-FLAG (FD), anti-HA (HD), and monoclonal anti-Ago2 (mAD) antibodies with the number of the biological replicates indicated in parentheses. Immunoprecipitated transcripts were defined as transcripts that have a high rank (top 2.5%) in each of the indicated IP experiments. Number of transcripts is shown on the left. The left panel additionally shows seed enrichment (SG1, NSG3) for the transcripts clustered in the right panel. (B) Seed complement enrichment for HEK293 miRNAs in immunoprecipitated transcripts compared with transcripts that were not enriched in the IP. The top 5% most enriched transcripts and random subsets containing the same number of not enriched transcripts obtained from each IPs were scanned for the presence of sequences complementary to the three HEK293 seed groups (SG1–SG3) and the control nonseed group 1–3 (NSG1–NSG3), representing control seed groups of miRNAs not expressed in HEK293 cells. IPs are described in A. (C) Enrichment for conserved seed complements for HEK293 miRNAs in immunoprecipitated transcripts. CSG1 (conserved seed group 1) refers to conserved seed hits of the top five expressed miRNAs. SG1 refers to seed hits of the top five expressed miRNA seed groups without conservation information, and NSG1 and NSG2 represent control miRNA seed groups. (D) Predicting miRNA profiles from AGO-immunoprecipitated transcripts. In the left panel a linear regression model was used to infer the activities of each miRNA seed family (nucleotides 1–8, 133 seed families in total; see Landgraf et al. 2007) based on the representation of transcripts carrying seed matches in AGO-IP. Seed families whose reverse complements are frequent in the immunoprecipitated transcripts receive high activity scores, whereas seed families that do not correlate with the IP data receive low scores. The right panel shows a histogram of the Pearson correlation coefficients resulting from the comparison of the miRNA “activity” predicted by the linear model with each of the measured miRNA expression profiles. The profiles that went into the construction of the histogram are from data of the 177 samples of Landgraf et al. (2007) and the HEK293 profile determined by 454 sequencing in this study (red triangle).

To demonstrate that the mRNP complexes isolated by the FLAG/HA-EIF2C1,-2,-3,-4 IPs indeed contained mRNAs targeted by miRNAs, the 3′ UTRs of the enriched transcripts were scanned for the presence of sequences complementary to the seed sequences of miRNAs expressed in HEK293. For this purpose, we clustered the seed sequences into seed groups (SG) according to the expression level of the corresponding miRNA family members in HEK293 cells (Table 2). SG1 comprises the seed sequences of the five most highly expressed miRNA families (S1–S5), whereas SG2 and SG3 contain the miRNA families that rank S6–S10 and S11–S15, respectively (Table 2). The density of seed-complementary motifs in the 3′ UTR of the enriched mRNAs was compared with the density of these motifs in a set of size-matched 3′ UTRs that were not enriched in the IP. A value of 1 indicates that the density of miRNA seed-complementary motifs is the same between immunoprecipitated transcripts and mRNAs that are expressed but not immunoprecipitated (control set). By performing repeated random selections of control transcripts, we obtained an estimate of the variance in the calculated enrichment in seed-complementary sites. Moreover, we compared the enrichment obtained for the most highly expressed miRNAs to that calculated for random subsets of miRNAs that were not expressed in HEK293 cells (NSG1–NSG3; see Supplemental Table 1). Figure 4B shows that the immunopurified mRNAs are enriched in sites complementary to SG1 through SG3, respectively, whereas no enrichment could be observed for control seed groups (NSG1–NSG3). The seed complement enrichment generally decreased with the relative expression of the corresponding miRNAs.

TABLE 2.

Expression of HEK293 miRNA seed families

graphic file with name 2580tbl2.jpg

To provide further evidence that the immunoprecipitated transcripts are indeed functional miRNA targets, we computed the enrichment of conserved miRNA seed complements using a method for estimating the probability that a miRNA binding site is under evolutionary selection (Gaidatzis et al. 2007). We found a significantly stronger enrichment in these sites, indicating that the IP captures functional miRNA targets, which were selected during evolution (Fig. 4C).

Yet another indication that our approach identifies bona fide miRNA targets is provided by a linear regression model, with which we attempt to predict the transcript enrichment in the IP in terms of the number of miRNA seed matches in the 3′ UTR of transcripts and a vector of miRNA-dependent weights, representing the concentrations of miRNAs in a cell. By fitting these miRNA-specific weights we derive a relative miRNA expression profile in the cell. More precisely, given that the prediction of miRNA target sites is based solely on the seed sequence of miRNAs, we cannot distinguish between different miRNAs carrying the same seed sequence; therefore, what we predict is the seed expression profile. We compared the predicted profile with the 177 profiles obtained from different tissues (Landgraf et al. 2007), and found that the predicted profile most strongly correlates with the experimentally determined HEK293 miRNA profile (Fig. 4D). Thus, the AGO-immunoprecipitated mRNAs enable one to reconstruct, to some extent, the miRNA expression profile of the sample. Figure 4D additionally shows a positive correlation coefficient between the predicted HEK293 miRNA expression profile and the profiles obtained experimentally from other types of cells, likely due to the presence of a subset of miRNAs with broad tissue distribution.

miRNA transfection and immunopurification of miR-122 targets

Another method to define miRNA targets relied on miRNA transfection experiments to identify miRNA-regulated mRNAs by measuring the reduction of transcripts caused by the addition of the exogenous miRNA (Lim et al. 2005; Grimson et al. 2007). To provide additional experimental support for our approach to identify miRNA-targeted mRNAs, we transfected FLAG/HA-EIF2C2-expressing cells with miR-122, a liver-specific miRNA, and compared the cellular RNA and RNA isolated from FLAG IPs with that of mock-transfected cells, similar to experiments performed previously (Karginov et al. 2007; Hendrickson et al. 2008). Northern analysis showed that miR-122 was transfected and recovered in the immunoprecipitate (Fig. 5A). Furthermore, mRNAs immunoprecipitated from the miR-122-transfected cells exhibit a strong enrichment in miR-122 seed-complementary sites (Fig. 5B), whereas the transcripts immunoprecipitated from mock-transfected cells were enriched only in sequences complementary to the HEK293 endogenous SG1–SG3, but not for the miR-122 seed.

FIGURE 5.

FIGURE 5.

miR-122 transfection and immunoprecipitation of miR-122 mRNA targets. (A) Cells stably expressing FLAG/HA-EIF2C2 were mock transfected and transfected with a miR-122 duplex. Fifteen hours after transfection cells were lysed and the epitope-tagged protein was immunoprecipitated from cytoplasmic extracts with anti-FLAG antibody. RNA was extracted from the cleared cell lysate and the IP. RNA was separated by 15% denaturing PAGE and probed for miR-16 and miR-122. A total of 10 and 4 fmol of synthetic miR-16 and miR-122, respectively, were loaded as standards. The miR-16 blot was reprobed with an oligonucleotide antisense to valine tRNA to ensure equal loading of lysate RNA. (B) Analysis of enriched mRNAs in Ago-IP. The EIF2C2-IP-enriched transcripts were scanned for the presence of sequences complementary to the three HEK293 seed groups (SG1–SG3), the control nonseed groups 1–3 (NSG1–NSG3), and the miR-122 seed in the 3′ UTR. The first and second bar group show the enrichment of complements in immunoprecipitated transcripts of the mock-transfected and mir-122-transfected cells, respectively.

In summary, these analyses demonstrated that our IP protocol was able to capture targets of both endogenously expressed, as well as transfected miRNAs, and that the immunoprecipitated mRNA population reflects the miRNA population present in the cell.

AGO and TNRC6 protein IPs enrich similar mRNAs with miRNA seed complements

Next, we compared the mRNA populations obtained by IP of proteins associated with AGO, including DICER, TNRC6A-C, IGF2BP1, PABPC4, and QKI as control (summarized in Supplemental Table 2). We found a striking similarity between the mRNAs purified by anti-FLAG-agarose IPs of epitope-tagged EIF2C1,-2,-3,-4 and TNRC6A-C (Fig. 4A), indicating that AGO and TNRC6 proteins are part of identical or largely overlapping mRNP complexes. This conclusion was also supported by the observation that the interaction between AGO and TNRC6 proteins is RNase resistant (Fig. 1). Furthermore, a set of about 600 transcripts was shown to be reproducibly associated with EIF2C1,-2,-3,-4 and TNRC6A,-B,-C proteins (Fig. 4A).

AGO and DICER IPs had some mRNAs in common. However, the overlap was far less significant than the one observed between AGO and the TNRC6-immunoprecipitated transcripts (Supplemental Fig. 4), although we could show a stable protein–protein interaction between DICER and EIF2C1 (Fig. 1).

The sets of transcripts immunoprecipitated by the IGF2BP and QKI proteins overlapped to some extend with the sets of transcripts immunoprecipitated by AGO and TNRC6. IGF2BP1 was identified in the double affinity purifications of AGO proteins and was shown to interact with AGO proteins in co-IP experiments in an RNA-dependent manner (data not shown).

In contrast, the PABPC4-IPed RNA did not reveal any significant overlap, although an RNA-mediated interaction with EIF2C1 could be observed (Fig. 1). In fact, some of the most enriched transcripts in PABPC4-immunoprecipitated are only present at low levels in the AGO-IPs (Supplemental Fig. 4).

As expected, TNRC6-immunoprecipitated transcripts showed a clear enrichment in miRNA seed-complementary sites analogous to AGO-immunoprecipitated transcripts (Fig. 4B,C), whereas no seed enrichment was observed for mRNAs in DICER, IGF2BP1, QKI, and PABPC4 IPs. These results indicated that AGO and TNRC6 protein-containing mRNP complexes associated with miRNA targets and are components of the miRNP silencing effector complexes. The seed enrichment analysis further suggested that neither QKI nor IGF2BP1 are part of the miRNA effector complexes defined by AGO and TNRC6 proteins. Interestingly, the PABPC4-IPed transcripts showed a significant reduction in evolutionarily selected seed-complementary sites, suggesting that the poly(A) tails of miRNA-targeted mRNAs might be less accessible for the epitope-tagged PABPC4, consistent with the observation of shorter poly(A) tails on miRNA-targeted mRNAs (Behm-Ansmant et al. 2006; Giraldez et al. 2006; Wu et al. 2006; Wakiyama et al. 2007).

Experimental validation of miRNA targets

To verify experimentally if the enriched mRNAs are indeed negatively regulated by HEK293 miRNAs, we generated dual luciferase reporter constructs containing the 3′ UTRs of 16 of the 20 most enriched mRNAs in EIF2C1,-2,-3,-4 IPs (Supplemental Table 2). We were unable to amplify the other four unusually long 3′ UTRs. The luciferase activity increased reproducibly for 14 of the 16 reporter constructs when cotransfected with a pool of antisense 2′-O-methyl oligoribonucleotides to block the activity of the 12 most highly expressed HEK293 miRNAs, compared with a transfection with an unspecific 2′-O-methyl control oligoribonucleotide (Fig. 6). The binding of the respective mRNAs to AGO protein complexes and repression by miRNAs confirmed that the IPed mRNAs are regulated by miRNAs in vivo.

FIGURE 6.

FIGURE 6.

miRNA-mediated repression of 3′ UTRs fused to luciferase reporter genes. Luciferase activity from HEK293 cells cotransfected with each reporter/control psiCHECK construct and 2′-O-methyl oligoribonucleotides antisense to 12 HEK293 miRNAs (red bars) was normalized to that from cotransfection of each reporter construct with 2′-O-methyl oligoribonucleotide antisense to Drosophila miRNA bantam (blue bars). Transfection of parental psiCHECK vector without inserted 3′ UTR (psiCHECK) and 3′ UTR with two artificial binding sites each for miR-16 and miR-196b (psiCHECK T) are indicated. Error bars represent standard deviation of experiment performed with six replicates. Normalized Renilla luciferase (RL) versus firefly luciferase (FL) activities are indicated.

DISCUSSION

In this present study, we provide the results of a comprehensive proteomic and ribonomic profiling study to elucidate the protein, the small RNA, and target mRNA composition of the four human AGO protein and other RNA silencing-related protein complexes.

Proteomic analysis

In contrast to other studies, our survey included the proteomic characterization of all four human AGO family members (EIF2C1,-2,-3,-4) and revealed that all of them associate with a similar set of helicases, hnRNPs, and mRNA-binding proteins, as well as DICER and members of the TNRC6 protein family. The composition of EIF2C1 and EIF2C2 RNP complexes showed a significant overlap with interactors recently identified (Hock et al. 2007). We were unable to detect several proteins that were identified by a single peptide in the MS analysis by Hock et al. (2007). Among these were several DEAD/DEXH helicases (DDX39, DDX46, DHX36) and the RNA-binding protein RBM4, which was shown to be involved in miRNA-mediated regulation (Hock et al. 2007). This suggests that these proteins might only be transiently associated with AGO protein-containing complexes.

Of the other identified interactors, DICER and the members of the TNRC6 protein family have been shown to be essential in miRNA precursor processing (Bernstein et al. 2001; Grishok et al. 2001; Hutvágner et al. 2001) and miRNA-mediated translational repression and mRNA decay, respectively (Jakymiw et al. 2005; Liu et al. 2005; Meister et al. 2005; Rehwinkel et al. 2005; Behm-Ansmant et al. 2006; Eulalio et al. 2008). We provide further biochemical evidence that these proteins form stable complexes with AGO proteins.

DICER cosediments with EIF2C2 in the lighter density complex I, which likely represents the RISC loading complex (RLC). Our analysis of the DICER immunoprecipitated mRNAs suggests that DICER and its associated proteins are not part of the miRNA effector complexes, since the DICER-enriched transcripts shared no significant similarity to AGO-immunoprecipitated miRNA targets. Recent in vitro reconstitution experiments of the human RLC indicated that EIF2C2 dissociates from this complex once it is loaded with miRNAs, suggesting that AGO proteins disengage from the RLC before targeting mRNAs (MacRae et al. 2008).

We showed previously that TNRC6B interacts with EIF2C1 and EIF2C2, and that the protein localizes to P-bodies and is required for miRNA-mediated mRNA cleavage (Meister et al. 2005). Here, we provide evidence that the human TNRC6 proteins interact stably with the mRNA-bound miRNP effector complexes, since miRNA-targeted mRNAs are highly enriched in TNRC6 proteins complexes. Zhang et al. (2007) made similar observations for AIN-1 and AIN-2, two C. elegans proteins homologous to the human TNRC6 protein family. As observed for the two TNRC6 proteins in the nematode, the mRNAs immunoprecipitated by the three human TNRC6 members overlapped extensively, supporting the idea that the different proteins may act redundantly and regulate a highly similar set of mRNAs.

Our proteomic analyses further revealed that AGO-containing complexes (as well as TNRC6-containing complexes) contained components of the translational machinery. We identified ribosomal proteins belonging to the small and large subunits, whereas other proteins commonly associated with the P-bodies, like proteins of the decapping machinery (Jakymiw et al. 2005; Liu et al. 2005; Meister et al. 2005) and the helicase DDX6 (Chu and Rana 2006), were not detected. In addition, we observed cosedimentation of EIF2C2 with polysomes, and association with the small subunit of the ribosome (complex II) as well as the 80S ribosome (complex III). Our results contrast previous observations (Hock et al. 2007), which described sedimentation coefficients for AGO mRNP complexes II and III of ∼19S and 25–30S, respectively.

The observed localization of EIF2C2 to polysomes is consistent with several previous studies suggesting that translation is regulated by miRNPs at a step after initiation. In C. elegans, miRNA-targeted mRNAs were discovered in polysomes (Olsen and Ambros 1999; Seggerson et al. 2002). Additional studies demonstrated that siRNAs and miRNAs associate with polysomes (Kim et al. 2004; Nelson et al. 2004; Maroney et al. 2006; Nottrott et al. 2006; Petersen et al. 2006). Furthermore, the D. melanogaster 80S RISC complex was found to contain ribosomes (Pham et al. 2004).

miRNA target identification

The isolation of mRNAs bound to EIF2C1,-2,-3,-4 and TNRC6A,-B,-C-containing protein complexes resulted in the identification of a highly similar set of transcripts, whose 3′ UTRs are enriched for sequences that are complementary to miRNA seed sequences (positions 1–8).

Our biochemical approach identified miRNA targets based on the association of these transcripts with human AGO and TNRC6 protein complexes. These results are consistent with the previously described physical interactions and functional relationship of TNRC6A and B with AGO proteins (Jakymiw et al. 2005; Liu et al. 2005; Meister et al. 2005; Wakiyama et al. 2007; Eulalio et al. 2008). TNRC6C is a novel component of this interaction network. Moreover, the nearly identical sets of mRNAs coimmunoprecipitated by the different AGO and TNRC6 proteins, respectively, suggest a high degree of redundancy between the members of each protein family with regard to their mRNA targets, and may also exhibit redundancy on a functional level. The targeting redundancy of human AGO proteins is reflected in their ability to bind identical miRNAs (Liu et al. 2004; Meister et al. 2004).

Our approach focused on the identification of physiologically relevant targets of endogenously expressed miRNAs in human cells, as similar efforts have in Drosophila Schneider S2 cells and in C. elegans (Easow et al. 2007; Zhang et al. 2007). Other IP approaches in human cells focused on the identification of mRNA targets of single, transfected miRNAs that are bound to EIF2C2 (Karginov et al. 2007; Hendrickson et al. 2008).

We identified about 600 transcripts to be regulated in HEK293 cells, a number that is substantially lower than the number of predicted targets for miRNAs expressed in HEK293 cells and the about 1200 and 2500 enriched transcripts recently identified by IP of EIF2C2 in HEK293 cells (Karginov et al. 2007; Hendrickson et al. 2008). Our seed complement density analysis revealed that only the highly expressed miRNAs leave detectable signatures in the immunoprecipitated mRNAs, suggesting that the highly expressed miRNAs have a biologically significant function or are amenable to experimental analysis of miRNA regulation in biological systems.

The dual luciferase reporter assay results demonstrated that the absolute degree of regulation mediated by miRNAs varies between 1.5- and 3.5-fold. The magnitude of regulation is likely dependent on the number of miRNP effector complexes associated with the target mRNA, since the level of repression achieved is dependent on both the amount of mRNA and the number of available miRNP complexes (Doench and Sharp 2004). In addition, miRNPs function in the context of other post-transcriptional regulatory networks. AGO and TNRC6 protein complexes contain a number of mRNA-binding proteins (mRBPs) with functions in mRNA transport, stabilization, and translation, such as IGF2BP1 (Yisraeli 2005). Recent observations indicate that the interplay of miRNPs and trans-acting mRBPs regulate the abundance of mRNAs and the rate of their translation in response to stress and environmental signals (Leung and Sharp 2007). The following examples illustrate the competition or synergism of various mRNA-binding factors and the modular architecture of post-transcriptional regulatory elements. HuR/ELAVL1, an AU-rich element-binding protein also identified in our AGO IPs, relieves a miRNA-induced inhibition of translation in cells subjected to different stress conditions (Bhattacharyya et al. 2006). The RNA-binding protein Dnd1 inhibits miRNA access to target mRNA (Kedde et al. 2007). Furthermore, EIF2C2 and fragile X mental retardation-related protein 1 were shown to activate the translation of tumor necrosis factor-α upon cell cycle arrest (Vasudevan et al. 2007).

miRNA target predictions have been valuable tools for defining putative miRNA-regulated transcripts. However, prediction methods have several kinds of limitations. First, the rules for base-pairing of the miRNA seed with the evolutionary conserved target transcript sites were shown not to be exclusive indicators of miRNA/mRNA interactions (Krutzfeldt et al. 2005; Didiano and Hobert 2006). Annotations of mRNA sequences, in particular 3′ UTR sequences, might be incomplete or inaccurate (Rajewsky 2006). Prediction methods have mainly focused on target identification in 3′ UTR sequences and the impact of potential binding sites in other locations of the transcript remains to be determined. Finally, the information about expression levels and patterns of miRNAs and mRNAs were not taken into consideration in any of the computational predictions. Experimental data sets that allow for the testing of these predictions have recently been reported (Baek et al. 2008; Selbach et al. 2008).

The first experimental approaches to identify miRNA-regulated mRNAs on a transcriptome-wide level involved recording the stability of transcripts upon miRNA expression and inhibition (Krutzfeldt et al. 2005; Lim et al. 2005). These approaches have their limitations, as only target mRNAs were revealed that show changes in steady-state levels. Furthermore, these methods, as well as recently reported quantitative proteomic approaches (Baek et al. 2008; Selbach et al. 2008), rely on the introduction of miRNAs or miRNA-inhibitory “antagomirs,” which could have an unintended impact on physiological regulatory interactions between miRNAs and their target mRNAs and are limited to the study of the altered miRNA.

Our experimental approach, as well as other IP-based methods, which identify miRNA-associated transcripts independent of the effect on mRNA stability, (Beitzinger et al. 2007; Easow et al. 2007; Karginov et al. 2007; Zhang et al. 2007; Hendrickson et al. 2008) allow the discovery of miRNA-regulated mRNAs without prior knowledge of the rules of interaction between miRNA and target mRNAs and the evolutionary conservation of this regulatory relationship. In our case, the target identification by IP indicates the feasibility of the method to define physiologically relevant miRNA-regulated transcripts without overexpression or depletion of miRNAs. We expect that a more detailed analysis of these experimentally determined miRNA-targeted transcripts will provide further insight into the sequence and structural features that affect the regulation of these transcripts in their natural environment, as well as the combinatorial regulation by different miRNA families.

In summary, we have shown that IPs of AGO and TNRC6 proteins resulted in an enrichment of miRNA targets, which correlated with the miRNA expression profile in HEK293 cells and was only detectable for the 15 most strongly expressed miRNA seed families. Furthermore, our IP experiments using a monoclonal antibody against human EIF2C2 demonstrate the possibility to identify miRNA-targeted mRNAs in human cells and tissues, which allows one to uncover the physiological miRNA regulatory networks in humans. Ultimately, additional biochemical approaches are necessary to determine the exact in vivo binding sites of mRBPs and miRNPs, thus establishing a comprehensive high-resolution map of the mRNA-RBP/miRNP interaction landscape. Defining this “messenger ribonucleoprotein code” and its sequence requirements will provide a better understanding of the mechanisms mediating post-transcriptional control of gene expression and the complexity of this regulatory network.

MATERIAL AND METHODS

Plasmids and antibodies

Plasmids and antibodies are described in the Supplemental Materials. Plasmids are made available through Addgene (www.addgene.com).

Double-affinity purification of FLAG/HA-tagged protein complexes

Cells were harvested from 8–10 L of suspension culture by centrifugation, washed in PBS, and frozen after resuspension of the cell pellet (about 15 mL) in an equal volume of freezing buffer (50 mM HEPES-KOH at pH 7.4, 1 mM EGTA, 1 mM MgCl2, 100 mM KCl, and 10% glycerol). Frozen cells were pulverized with a mortar and lysed in 3 mL/g frozen cell suspension (75 mM HEPES-KOH at pH 7.4, 1.5 mM EGTA, 1.5 mM MgCl2, 150 mM KCl, and 15% glycerol, 0.075% NP-40, EDTA-free Protease Inhibitor Cocktail [Roche]) by sonication. The cell lysate was cleared by centrifugation at 25.000g for 45 min at 4°C. IP was performed for 3 h at 4°C using the FLAG-specific monoclonal antibody M2 covalently coupled to agarose beads (Sigma). The immunoprecipitate was washed six times in Wash300 buffer (50 mM HEPES-KOH at pH 7.4, 1 mM EGTA, 1 mM MgCl2, 3 mM KCl, and 10% glycerol, 0.05% NP-40, 1 mM DTT, EDTA-free Protease Inhibitor Cocktail) and twice in Wash100 buffer (as Wash300 buffer but with 100 mM KCl). FLAG-tagged proteins were eluted by addition of 500 μg of synthetic 1×FLAG peptide per 1 mL of Wash100 buffer, followed by rotation for 1 h at 4°C. FLAG-eluted proteins were further immunoprecipitated using μMACS HA-beads (Miltenyi) for 1 h and eluted in 1× SDS-loading buffer. Eluted proteins were separated on 4%–12% NUPAGE SDS-PAGE (Invitrogen) and stained with Gelcode Blue (Pierce).

Mass spectrometry

Protein bands were excised from stained SDS–polyacrylamide gels, washed, and digested with trypsin. Tryptic peptides were isolated and liquid chromatography–tandem mass spectrometry analyses were performed on a Nano HPLC coupled to an ion-trap (ThermoElectron Finnigan LTQ) mass spectrometer. Peptide fragment ion spectra were searched against the nonredundant protein database for matches using MASCOT (Matrix Science).

Transfection and co-IP

For miR-122 transfection experiments, FLAG/HA-EIF2C2 cells were transfected with miR-122 duplex or mock transfected using Lipofectamine RNAiMAX as described by manufacturer.

For transient transfection of cells, stably expressing FLAG/HA-EIF2C1, pDESTmyc constructs were delivered using Lipofectamine 2000 as described by the manufacturer. Cells were lysed in 50 mM HEPES-KOH (pH 7.4), 150 mM KCl, 2 mM EDTA, 0.5 mM DTT, 1 mM NaF, and 0.5% NP-40. The cell lysate was cleared. For RNase treatment, lysate was incubated for 20 min at room temperature with RNase T1 at 1 U/μL. IP was performed as described above with ANTI-FLAG M2-agarose beads. The immunoprecipitate was washed 3× in Wash300 buffer and 2× in Wash100 buffer. FLAG-tagged proteins were eluted by addition 1×FLAG peptide. Eluted proteins were loaded separated by SDS-PAGE, blotted, and probed with anti-HA and anti-myc antibodies. Primary antibodies were detected with HRP-conjugated secondary antibodies using the ECL Plus Western Blotting Detection System (Amersham).

Sucrose density fractionation

Cells were washed with ice-cold PBS, lysed in 10 mM HEPES-KOH (pH 7.4), 150 mM KCl, 5 mM MgCl2 or 30 mM EDTA, and 0.5% NP-40, and centrifuged for 5 min to pellet the nuclei. Polysome fractionations of 200 μL of cell extract were performed in 12 mL of gradients from 10%–45% (w/v) sucrose in 10 mM HEPES-KOH (pH 7.4), 150 mM KCl supplemented with either 5 mM MgCl2 or 30 mM EDTA. Lysates were separated by centrifugation at 38,000 rpm for 1.5 h in a SW41 rotor at 4°C. FLAG-peptide-eluted protein complexes from immunoprecipitates were separated by centrifugation at 38,000 rpm for 8.5 h in a SW41 rotor at 4°C in 10%–45% (w/v) sucrose gradient in 10 mM HEPES (pH 7.4), 150 mM KCl, and 5 mM MgCl2.

RNA isolation from cell lysate and FLAG-protein immunoprecipitates

About 15 mL of cell pellet were lysed in three cell pellet volumes of 50 mM HEPES-KOH (pH 7.4), 150 mM KCl, 2 mM EDTA, 0.5 mM DTT, 1 mM NaF, and 0.5% NP-40. RNA from the lysate was isolated by adding 3 vol of RNA extraction solution (4 M guanidinium isothiocyanate, 25 mM sodium citrate, 0.5% N-lauroylsarcosinate, 50 mM β-mercaptoethanol, and 50% acidic phenol) and 0.2 volumes of chloroform. RNA was ethanol-precipitated from the aqueous phase. FLAG/HA-tagged or endogenous protein was immunoprecipitated with anti-FLAG M2 agarose beads or antibody-bound Protein G Dynabeads (Invitrogen). Beads were washed three times with 50 mM HEPES-KOH (pH 7.4), 300 mM KCl, 2 mM EDTA, 0.5 mM DTT, 1 mM NaF, and 0.05% NP-40. RNA isolation from immunoprecipitated RNPs was performed as described previously (Meister et al. 2004). RNA for microarray analysis was further purified using RNeasy mini spin columns (QIAGEN). Quality of RNA was assessed using the Agilent Bioanalyzer.

Microarray analysis

A total of 2 μg of purified total RNA from HEK293 cell lysate or from immunoprecipitated RNPs was used in the One-Cycle Eukaryotic Target Labeling Assay (Affymetrix) according to the manufacturer's protocol. Biotinylated cRNA targets were cleaned up, fragmented, and hybridized to Human Genome U133 Plus 2.0 Array (Affymetrix).

miRNA Northern hybridization

miRNA Northern analysis was performed as described previously (Landthaler et al. 2004).

Dual luciferase assay

HEK293 cells were cotransfected in 96-well format (40.000 cells/well) with 100 ng of the respective psiCHECK vector and a total of 120 pmoles of a pool of 2′-O-methyl oligoribonucleotides (anti-miR-15a, anti-miR-16, anti-miR-17, anti-miR-18, anti-miR-19a, anti-miR-19b, anti-miR-20a, anti-miR-20b, anti-miR-32, anti-mir-92, anti-miR-93, and anti-mir-196b) antisense to the respective miRNA or 120 pmoles of 2′-O-methyl oligoribonucleotide antisense to the Drosophila-specific miRNA bantam with Lipofectamine 2000 (Invitrogen). Cells were lysed in 1×Passive Lysis Buffer (Promega) 15 h after transfection and analyzed using the Dual-Luciferase Reporter System (Promega) as described by the manufacturer on a BIO-TEK Clarity 96-well plate reader with double injectors.

Computational analysis of microarray data

The probe intensities were normalized for all the microarray experiments using the bioconductor (see http://www.bioconductor.org and Gentleman et al. (2004) and gcRMA software (Wu et al. 2004). To quantify effects based on the frequency of seed-complementary sites in the 3′ UTRs, we selected, for each gene measured by the microarray, the transcript with median 3′ UTR length. This data set consisted of 14,658 transcripts. From all of the probe sets corresponding unambiguously to a given gene, we selected the one that responded best (exhibited the highest variance) across a large number of experiments performed on the Affymetrix platform that we used. This probe set was used to monitor the per-transcript expression level across our experiments. Based on the transcript expression levels computed this way, we calculated the Pearson correlation coefficient of mRNA expression between pairs of experiments.

Selection of transcripts for the analysis of target spectra of different RNA-binding proteins

In order to compare the target spectra of different RNA-binding proteins, we extracted from each sample the top 2.5% most enriched transcripts in a particular IP of a specific protein relative to lysate. This amounts to 195 transcripts per sample from a total of 7820, for which we had expression evidence in HEK293 cells. We used the union of these transcripts (1589) to perform hierarchical clustering. The clustering results are rather insensitive to variations in the number of transcripts that we select from each sample, but larger transcript sets are more difficult to visualize.

Construction of miRNA seed groups

The miRNA profile of HEK293 cells was first determined as described in the Supplemental Materials.. We used this profile to rank the miRNA expression. To improve the reliability of our statistics of miRNA target site enrichment, we performed this analysis for one group of miRNAs at a time as follows. One of the most reliable indicators of a miRNA target site is the perfect complementarity to the 5′ end of the miRNA (Lewis et al. 2005). We therefore divided miRNAs sets based on the sequence they carry at positions 1–8. For each of these sets, we computed the total expression in HEK293 by summing the counts of all miRNAs within a set. Such a sequence set typically contains paralogous miRNA genes. We then ranked the expression of these sets, and constructed groups of five sequence sets at a time, starting from the highest to the lowest expressed sequence set. We called these groups SG1, SG2, SG3, etc. For reference, we also created random groups of five sequence sets (NSG1–NSG3) that were not expressed in HEK293 cells (see Supplemental Table 1).

Enrichment of miRNA target sites in transcripts that co-IP with RNA-binding proteins

We assessed the enrichment in miRNA target sites of the set of transcripts bound by a given RNA-binding protein as follows. We selected the top 5% transcripts that were IPed with the given RNA-binding protein. For each miRNA seed family (nucleotides 1–8), we determined the density (sites/nucleotides) of seed (1–7, 2–8, and 1–8) complementary sites in the given set of transcripts. This quantity was compared with the density of the same motifs in a random subset of transcripts of the same size, which were not enriched in the IP (background set). The background transcripts where chosen such that their fold change in the IP compared with lysate was in the 25%–75% range. By repeating the procedure of selecting background subsets of transcripts 100 times, we obtained means and standard deviations for motif enrichment in the various samples.

Enrichment of evolutionarily selected target sites in transcripts that coimmunoprecipitate with mRNA-binding proteins

The algorithm described by Gaidatzis et al. (2007) was used to obtain for each putative miRNA target site a posterior probability that the site is under selection. Then, we applied the method described in the previous paragraph to compute the enrichment by counting each site with its posterior probability of being under selection, rather than with a uniform weight of 1 as we did in the previous paragraph.

Computational inference of miRNA expression based on the mRNA IP data

We used a linear model to predict the IP enrichment of a transcript in terms of the number of binding sites that it has for individual miRNAs, as well as of a factor reflecting the “concentration” of the miRNA in the cell. The miRNAs to be considered were determined based on published miRNA expression data (Landgraf et al. 2007). From this data set we extracted the list of miRNA seeds (nucleotides 1–8) that accumulated at least 50 sequences cloned from all the sampled tissues. We considered these 133 unique seeds for relating the enrichment in the IP to the miRNA expression.

The enrichment in the IP was measured, the number of binding sites for individual miRNAs in each 3′ UTR was calculated as the number of matches to positions 1–7, 2–8, or 1–8 of the miRNA, and using Matlab we fitted the weights wm corresponding to all miRNAs in the expressions Inline graphic. Here, Ei is the logarithm of the fold change of transcript i in the IP compared with lysate, and Ni,m is the number of sites for miRNA m in the 3′ UTR of transcript i. We then compared these weights with the relative miRNA frequencies measured by 454 sequencing, and we obtained a Pearson correlation coefficient of 0.37. Similarly, we correlated the inferred weights with the miRNA concentrations measured in a variety of other samples from Landgraf et al. (2007).

SUPPLEMENTAL DATA

Supplemental material can be found at http://www.rnajournal.org.

ACKNOWLEDGMENTS

We are grateful to Haiteng Deng and Joseph Fernandenz (Rockefeller University Proteomics Resource Facility) for the MS analysis, Wenxiang Zhang and Connie Zhao (Genomics Resource Center) for mRNA array analysis, and the Memorial Sloan Kettering Sequencing Core for 454 sequencing. We thank Jean Hausser for sharing some of his programs and Manuel Ascano, Thalia Farazi, Stefan Juranek, and other members of the laboratory for comments on the manuscript. We also thank Millipore for antibodies. D.G. was supported by Swiss National Fund grant #3100A0-114001 to M.Z. M.L. was supported by an Irma T. Hirschl Postdoctoral Fellowship. T.T. is an HHMI investigator, and work in his laboratory was supported by NIH grants GM073047 and GM068476.

Footnotes

Article published online ahead of print. Article and publication date are at http://www.rnajournal.org/cgi/doi/10.1261/rna.1351608.

REFERENCES

  1. Ambros V. The functions of animal microRNAs. Nature. 2004;431:350–355. doi: 10.1038/nature02871. [DOI] [PubMed] [Google Scholar]
  2. Baek D., Villen J., Shin C., Camargo F.D., Gygi S.P., Bartel D.P. The impact of microRNAs on protein output. Nature. 2008;455:64–71. doi: 10.1038/nature07242. [DOI] [PMC free article] [PubMed] [Google Scholar]
  3. Bagga S., Bracht J., Hunter S., Massirer K., Holtz J., Eachus R., Pasquinelli A.E. Regulation by let-7 and lin-4 miRNAs results in target mRNA degradation. Cell. 2005;122:553–563. doi: 10.1016/j.cell.2005.07.031. [DOI] [PubMed] [Google Scholar]
  4. Balzer E., Moss E.G. Localization of the developmental timing regulator Lin28 to mRNP complexes, P-bodies and stress granules. RNA Biol. 2007;4:16–25. doi: 10.4161/rna.4.1.4364. [DOI] [PubMed] [Google Scholar]
  5. Barth S., Pfuhl T., Mamiani A., Ehses C., Roemer K., Kremmer E., Jaker C., Hock J., Meister G., Grasser F.A. Epstein-Barr virus-encoded microRNA miR-BART2 down-regulates the viral DNA polymerase BALF5. Nucleic Acids Res. 2008;36:666–675. doi: 10.1093/nar/gkm1080. [DOI] [PMC free article] [PubMed] [Google Scholar]
  6. Behm-Ansmant I., Rehwinkel J., Doerks T., Stark A., Bork P., Izaurralde E. mRNA degradation by miRNAs and GW182 requires both CCR4:NOT deadenylase and DCP1:DCP2 decapping complexes. Genes & Dev. 2006;20:1885–1898. doi: 10.1101/gad.1424106. [DOI] [PMC free article] [PubMed] [Google Scholar]
  7. Beitzinger M., Peters L., Zhu J.Y., Kremmer E., Meister G. Identification of human microRNA targets from isolated argonaute protein complexes. RNA Biol. 2007;4:76–84. doi: 10.4161/rna.4.2.4640. [DOI] [PubMed] [Google Scholar]
  8. Bernstein E., Caudy A.A., Hammond S.M., Hannon G.J. Role for a bidentate ribonuclease in the initiation step of RNA interference. Nature. 2001;409:363–366. doi: 10.1038/35053110. [DOI] [PubMed] [Google Scholar]
  9. Bhattacharyya S.N., Habermacher R., Martine U., Closs E.I., Filipowicz W. Relief of microRNA-mediated translational repression in human cells subjected to stress. Cell. 2006;125:1111–1124. doi: 10.1016/j.cell.2006.04.031. [DOI] [PubMed] [Google Scholar]
  10. Bushati N., Cohen S.M. microRNA Functions. Annu. Rev. Cell Dev. Biol. 2007;23:175–205. doi: 10.1146/annurev.cellbio.23.090506.123406. [DOI] [PubMed] [Google Scholar]
  11. Chapman E.J., Carrington J.C. Specialization and evolution of endogenous small RNA pathways. Nat. Rev. Genet. 2007;8:884–896. doi: 10.1038/nrg2179. [DOI] [PubMed] [Google Scholar]
  12. Chendrimada T.P., Gregory R.I., Kumaraswamy E., Norman J., Cooch N., Nishikura K., Shiekhattar R. TRBP recruits the Dicer complex to Ago2 for microRNA processing and gene silencing. Nature. 2005;436:740–744. doi: 10.1038/nature03868. [DOI] [PMC free article] [PubMed] [Google Scholar]
  13. Chu C.Y., Rana T.M. Translation repression in human cells by microRNA-induced gene silencing requires RCK/p54. PLoS Biol. 2006;4:e210. doi: 10.1371/journal.pbio.0040210. [DOI] [PMC free article] [PubMed] [Google Scholar]
  14. Denli A.M., Tops B.B., Plasterk R.H., Ketting R.F., Hannon G.J. Protein and RNA Composition of human AGO protein complexes processing of primary microRNAs by the Microprocessor complex. Nature. 2004;432:231–235. doi: 10.1038/nature03049. [DOI] [PubMed] [Google Scholar]
  15. Didiano D., Hobert O. Perfect seed pairing is not a generally reliable predictor for miRNA-target interactions. Nat. Struct. Mol. Biol. 2006;13:849–851. doi: 10.1038/nsmb1138. [DOI] [PubMed] [Google Scholar]
  16. Doench J.G., Sharp P.A. Specificity of microRNA target selection in translational repression. Genes & Dev. 2004;18:504–511. doi: 10.1101/gad.1184404. [DOI] [PMC free article] [PubMed] [Google Scholar]
  17. Easow G., Teleman A.A., Cohen S.M. Isolation of microRNA targets by miRNP immunopurification. RNA. 2007;13:1198–1204. doi: 10.1261/rna.563707. [DOI] [PMC free article] [PubMed] [Google Scholar]
  18. Eulalio A., Huntzinger E., Izaurralde E. GW182 interaction with Argonaute is essential for miRNA-mediated translational repression and mRNA decay. Nat. Struct. Mol. Biol. 2008;15:346–353. doi: 10.1038/nsmb.1405. [DOI] [PubMed] [Google Scholar]
  19. Gaidatzis D., van Nimwegen E., Hausser J., Zavolan M. Inference of miRNA targets using evolutionary conservation and pathway analysis. BMC Bioinformatics. 2007;8:69. doi: 10.1186/1471-2105-8-69. [DOI] [PMC free article] [PubMed] [Google Scholar]
  20. Gentleman R.C., Carey V.J., Bates D.M., Bolstad B., Dettling M., Dudoit S., Ellis B., Gautier L., Ge Y., Gentry J., et al. Bioconductor: Open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. doi: 10.1186/gb-2004-5-10-r89. [DOI] [PMC free article] [PubMed] [Google Scholar]
  21. Giraldez A.J., Mishima Y., Rihel J., Grocock R.J., Van Dongen S., Inoue K., Enright A.J., Schier A.F. Zebrafish MiR-430 promotes deadenylation and clearance of maternal mRNAs. Science. 2006;312:75–79. doi: 10.1126/science.1122689. [DOI] [PubMed] [Google Scholar]
  22. Gregory R.I., Yan K.P., Amuthan G., Chendrimada T., Doratotaj B., Cooch N., Shiekhattar R. The microprocessor complex mediates the genesis of microRNAs. Nature. 2004;432:235–240. doi: 10.1038/nature03120. [DOI] [PubMed] [Google Scholar]
  23. Grimson A., Farh K.K., Johnston W.K., Garrett-Engele P., Lim L.P., Bartel D.P. MicroRNA targeting specificity in mammals: Determinants beyond seed pairing. Mol. Cell. 2007;27:91–105. doi: 10.1016/j.molcel.2007.06.017. [DOI] [PMC free article] [PubMed] [Google Scholar]
  24. Grishok A., Pasquinelli A.E., Conte D., Li N., Parrish S., Ha I., Baillie D.L., Fire A., Ruvkun G., Mello C.C. Genes and mechanisms related to RNA interference regulate expression of the small temporal RNAs that control C. elegans developmental timing. Cell. 2001;106:23–34. doi: 10.1016/s0092-8674(01)00431-7. [DOI] [PubMed] [Google Scholar]
  25. Grun D., Wang Y.L., Langenberger D., Gunsalus K.C., Rajewsky N. microRNA target predictions across seven Drosophila species and comparison to mammalian targets. PLoS Comput. Biol. 2005;1:e13. doi: 10.1371/journal.pcbi.0010013. [DOI] [PMC free article] [PubMed] [Google Scholar]
  26. Haase A.D., Jaskiewicz L., Zhang H., Laine S., Sack R., Gatignol A., Filipowicz W. TRBP, a regulator of cellular PKR and HIV-1 virus expression, interacts with Dicer and functions in RNA silencing. EMBO Rep. 2005;6:961–967. doi: 10.1038/sj.embor.7400509. [DOI] [PMC free article] [PubMed] [Google Scholar]
  27. Hammond S.M., Boettcher S., Caudy A.A., Kobayashi R., Hannon G.J. Argonaute2, a link between genetic and biochemical analyses of RNAi. Science. 2001;293:1146–1150. doi: 10.1126/science.1064023. [DOI] [PubMed] [Google Scholar]
  28. Han J., Lee Y., Yeom K.H., Kim Y.K., Jin H., Kim V.N. The Drosha–DGCR8 complex in primary microRNA processing. Genes & Dev. 2004;18:3016–3027. doi: 10.1101/gad.1262504. [DOI] [PMC free article] [PubMed] [Google Scholar]
  29. Han J., Lee Y., Yeom K.H., Nam J.W., Heo I., Rhee J.K., Sohn S.Y., Cho Y., Zhang B.T., Kim V.N. Molecular basis for the recognition of primary microRNAs by the Drosha-DGCR8 complex. Cell. 2006;125:887–901. doi: 10.1016/j.cell.2006.03.043. [DOI] [PubMed] [Google Scholar]
  30. Hendrickson D.G., Hogan D.J., Herschlag D., Ferrell J.E., Brown P.O. Systematic identification of mRNAs recruited to argonaute 2 by specific microRNAs and corresponding changes in transcript abundance. PLoS One. 2008;3:e2126. doi: 10.1371/journal.pone.0002126. [DOI] [PMC free article] [PubMed] [Google Scholar]
  31. Hock J., Weinmann L., Ender C., Rudel S., Kremmer E., Raabe M., Urlaub H., Meister G. Proteomic and functional analysis of Argonaute-containing mRNA-protein complexes in human cells. EMBO Rep. 2007;8:1052–1060. doi: 10.1038/sj.embor.7401088. [DOI] [PMC free article] [PubMed] [Google Scholar]
  32. Humphreys D.T., Westman B.J., Martin D.I., Preiss T. MicroRNAs control translation initiation by inhibiting eukaryotic initiation factor 4E/cap and poly(A) tail function. Proc. Natl. Acad. Sci. 2005;102:16961–16966. doi: 10.1073/pnas.0506482102. [DOI] [PMC free article] [PubMed] [Google Scholar]
  33. Hutvágner G., Zamore P.D. A microRNA in a multiple-turnover RNAi enzyme complex. Science. 2002;297:2056–2060. doi: 10.1126/science.1073827. [DOI] [PubMed] [Google Scholar]
  34. Hutvágner G., McLachlan J., Bálint É., Tuschl T., Zamore P.D. A cellular function for the RNA interference enzyme Dicer in small temporal RNA maturation. Science. 2001;93:834–838. doi: 10.1126/science.1062961. [DOI] [PubMed] [Google Scholar]
  35. Jakymiw A., Lian S., Eystathioy T., Li S., Satoh M., Hamel J.C., Fritzler M.J., Chan E.K. Disruption of GW bodies impairs mammalian RNA interference. Nat. Cell Biol. 2005;7:1267–1274. doi: 10.1038/ncb1334. [DOI] [PubMed] [Google Scholar]
  36. Jing Q., Huang S., Guth S., Zarubin T., Motoyama A., Chen J., Di Padova F., Lin S.C., Gram H., Han J. Involvement of microRNA in AU-rich element-mediated mRNA instability. Cell. 2005;120:623–634. doi: 10.1016/j.cell.2004.12.038. [DOI] [PubMed] [Google Scholar]
  37. John B., Enright A.J., Aravin A., Tuschl T., Sander C., Marks D. Human miRNA targets. PLoS Biol. 2004;2:e363. doi: 10.1371/journal.pbio.0020363. [DOI] [PMC free article] [PubMed] [Google Scholar]
  38. Jones-Rhoades M.W., Bartel D.P., Bartel B. MicroRNAS and their regulatory roles in plants. Annu. Rev. Plant Biol. 2006;57:19–53. doi: 10.1146/annurev.arplant.57.032905.105218. [DOI] [PubMed] [Google Scholar]
  39. Karginov F.V., Conaco C., Xuan Z., Schmidt B.H., Parker J.S., Mandel G., Hannon G.J. A biochemical approach to identifying microRNA targets. Proc. Natl. Acad. Sci. 2007;104:19291–19296. doi: 10.1073/pnas.0709971104. [DOI] [PMC free article] [PubMed] [Google Scholar]
  40. Kedde M., Strasser M.J., Boldajipour B., Vrielink J.A., Slanchev K., le Sage C., Nagel R., Voorhoeve P.M., van Duijse J., Orom U.A., et al. RNA-binding protein Dnd1 inhibits microRNA access to target mRNA. Cell. 2007;131:1273–1286. doi: 10.1016/j.cell.2007.11.034. [DOI] [PubMed] [Google Scholar]
  41. Ketting R.F., Fischer S.E., Bernstein E., Sijen T., Hannon G.J., Plasterk R.H. Dicer functions in RNA interference and in synthesis of small RNA involved in developmental timing in C. elegans . Genes & Dev. 2001;15:2654–2659. doi: 10.1101/gad.927801. [DOI] [PMC free article] [PubMed] [Google Scholar]
  42. Khvorova A., Reynolds A., Jayasena S.D. Functional siRNAs and miRNAs exhibit strand bias. Cell. 2003;115:209–216. doi: 10.1016/s0092-8674(03)00801-8. [DOI] [PubMed] [Google Scholar]
  43. Kim J., Krichevsky A., Grad Y., Hayes G.D., Kosik K.S., Church G.M., Ruvkun G. Identification of many microRNAs that copurify with polyribosomes in mammalian neurons. Proc. Natl. Acad. Sci. 2004;101:360–365. doi: 10.1073/pnas.2333854100. [DOI] [PMC free article] [PubMed] [Google Scholar]
  44. Krek A., Grun D., Poy M.N., Wolf R., Rosenberg L., Epstein E.J., MacMenamin P., da Piedade I., Gunsalus K.C., Stoffel M., et al. Combinatorial microRNA target predictions. Nat. Genet. 2005;37:495–500. doi: 10.1038/ng1536. [DOI] [PubMed] [Google Scholar]
  45. Krutzfeldt J., Rajewsky N., Braich R., Rajeev K.G., Tuschl T., Manoharan M., Stoffel M. Silencing of microRNAs in vivo with “antagomirs.”. Nature. 2005;438:685–689. doi: 10.1038/nature04303. [DOI] [PubMed] [Google Scholar]
  46. Landgraf P., Rusu M., Sheridan R., Sewer A., Iovino N., Aravin A., Pfeffer S., Rice A., Kamphorst A.O., Landthaler M., et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129:1401–1414. doi: 10.1016/j.cell.2007.04.040. [DOI] [PMC free article] [PubMed] [Google Scholar]
  47. Landthaler M., Yalcin A., Tuschl T. The human DiGeorge syndrome critical region gene 8 and Its D. melanogaster homolog are required for miRNA biogenesis. Curr. Biol. 2004;14:2162–2167. doi: 10.1016/j.cub.2004.11.001. [DOI] [PubMed] [Google Scholar]
  48. Lee Y., Han J., Yeom K.H., Jin H., Kim V.N. Drosha in primary microRNA processing. Cold Spring Harb. Symp. Quant. Biol. 2006a;71:51–57. doi: 10.1101/sqb.2006.71.041. [DOI] [PubMed] [Google Scholar]
  49. Lee Y., Hur I., Park S.Y., Kim Y.K., Suh M.R., Kim V.N. The role of PACT in the RNA silencing pathway. EMBO J. 2006b;25:522–532. doi: 10.1038/sj.emboj.7600942. [DOI] [PMC free article] [PubMed] [Google Scholar]
  50. Leung A.K., Sharp P.A. microRNAs: A safeguard against turmoil? Cell. 2007;130:581–585. doi: 10.1016/j.cell.2007.08.010. [DOI] [PubMed] [Google Scholar]
  51. Lewis B.P., Burge C.B., Bartel D.P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120:15–20. doi: 10.1016/j.cell.2004.12.035. [DOI] [PubMed] [Google Scholar]
  52. Lim L.P., Lau N.C., Garrett-Engele P., Grimson A., Schelter J.M., Castle J., Bartel D.P., Linsley P.S., Johnson J.M. Microarray analysis shows that some microRNAs down-regulate large numbers of target mRNAs. Nature. 2005;433:769–773. doi: 10.1038/nature03315. [DOI] [PubMed] [Google Scholar]
  53. Lingel A., Simon B., Izaurralde E., Sattler M. Nucleic acid 3′-end recognition by the Argonaute2 PAZ domain. Nat. Struct. Mol. Biol. 2004;11:576–577. doi: 10.1038/nsmb777. [DOI] [PubMed] [Google Scholar]
  54. Linsley P.S., Schelter J., Burchard J., Kibukawa M., Martin M.M., Bartz S.R., Johnson J.M., Cummins J.M., Raymond C.K., Dai H., et al. Transcripts targeted by the microRNA-16 family cooperatively regulate cell cycle progression. Mol. Cell. Biol. 2007;27:2240–2252. doi: 10.1128/MCB.02005-06. [DOI] [PMC free article] [PubMed] [Google Scholar]
  55. Liu J., Carmell M.A., Rivas F.V., Marsden C.G., Thomson J.M., Song J.J., Hammond S.M., Joshua-Tor L., Hannon G.J. Argonaute2 is the catalytic engine of mammalian RNAi. Science. 2004;305:1437–1441. doi: 10.1126/science.1102513. [DOI] [PubMed] [Google Scholar]
  56. Liu J., Rivas F.V., Wohlschlegel J., Yates J.R., 3rd, Parker R., Hannon G.J. A role for the P-body component GW182 in microRNA function. Nat. Cell Biol. 2005;7:1261–1266. doi: 10.1038/ncb1333. [DOI] [PMC free article] [PubMed] [Google Scholar]
  57. Lund E., Guttinger S., Calado A., Dahlberg J.E., Kutay U. Nuclear export of microRNA precursors. Science. 2004;303:95–98. doi: 10.1126/science.1090599. [DOI] [PubMed] [Google Scholar]
  58. Ma J.B., Ye K., Patel D.J. Structural basis for overhang-specific small interfering RNA recognition by the PAZ domain. Nature. 2004;429:318–322. doi: 10.1038/nature02519. [DOI] [PMC free article] [PubMed] [Google Scholar]
  59. Ma J.B., Yuan Y.R., Meister G., Pei Y., Tuschl T., Patel D.J. Structural basis for 5′-end-specific recognition of guide RNA by the A. fulgidus Piwi protein. Nature. 2005;434:666–670. doi: 10.1038/nature03514. [DOI] [PMC free article] [PubMed] [Google Scholar]
  60. MacRae I.J., Ma E., Zhou M., Robinson C.V., Doudna J.A. In vitro reconstitution of the human RISC-loading complex. Proc. Natl. Acad. Sci. 2008;105:512–517. doi: 10.1073/pnas.0710869105. [DOI] [PMC free article] [PubMed] [Google Scholar]
  61. Maroney P.A., Yu Y., Fisher J., Nilsen T.W. Evidence that microRNAs are associated with translating messenger RNAs in human cells. Nat. Struct. Mol. Biol. 2006;13:1102–1107. doi: 10.1038/nsmb1174. [DOI] [PubMed] [Google Scholar]
  62. Martinez J., Patkaniowska A., Urlaub H., Lührmann R., Tuschl T. Single-stranded antisense siRNAs guide target RNA cleavage in RNAi. Cell. 2002;110:563–574. doi: 10.1016/s0092-8674(02)00908-x. [DOI] [PubMed] [Google Scholar]
  63. Meister G., Tuschl T. Mechanisms of gene silencing by double-stranded RNA. Nature. 2004;431:343–349. doi: 10.1038/nature02873. [DOI] [PubMed] [Google Scholar]
  64. Meister G., Landthaler M., Patkaniowska A., Dorsett Y., Teng G., Tuschl T. Human Argonaute2 mediates RNA cleavage targeted by miRNAs and siRNAs. Mol. Cell. 2004;15:185–197. doi: 10.1016/j.molcel.2004.07.007. [DOI] [PubMed] [Google Scholar]
  65. Meister G., Landthaler M., Peters L., Chen P.Y., Urlaub H., Luhrmann R., Tuschl T. Identification of novel argonaute-associated proteins. Curr. Biol. 2005;15:2149–2155. doi: 10.1016/j.cub.2005.10.048. [DOI] [PubMed] [Google Scholar]
  66. Mourelatos Z., Dostie J., Paushkin S., Sharma A., Charroux B., Abel L., Rappsilber J., Mann M., Dreyfuss G. miRNPs: A novel class of ribonucleoproteins containing numerous microRNAs. Genes & Dev. 2002;16:720–728. doi: 10.1101/gad.974702. [DOI] [PMC free article] [PubMed] [Google Scholar]
  67. Nelson P.T., Hatzigeorgiou A.G., Mourelatos Z. miRNP:mRNA association in polyribosomes in a human neuronal cell line. RNA. 2004;10:387–394. doi: 10.1261/rna.5181104. [DOI] [PMC free article] [PubMed] [Google Scholar]
  68. Nottrott S., Simard M.J., Richter J.D. Human let-7a miRNA blocks protein production on actively translating polyribosomes. Nat. Struct. Mol. Biol. 2006;13:1108–1114. doi: 10.1038/nsmb1173. [DOI] [PubMed] [Google Scholar]
  69. Olsen P.H., Ambros V. The lin-4 regulatory RNA controls developmental timing in Caenorhabditis elegans by blocking LIN-14 protein synthesis after the initiation of translation. Dev. Biol. 1999;216:671–680. doi: 10.1006/dbio.1999.9523. [DOI] [PubMed] [Google Scholar]
  70. Orom U.A., Lund A.H. Isolation of microRNA targets using biotinylated synthetic microRNAs. Methods. 2007;43:162–165. doi: 10.1016/j.ymeth.2007.04.007. [DOI] [PubMed] [Google Scholar]
  71. Parker J.S., Roe S.M., Barford D. Structural insights into mRNA recognition from a PIWI domain-siRNA guide complex. Nature. 2005;434:663–666. doi: 10.1038/nature03462. [DOI] [PMC free article] [PubMed] [Google Scholar]
  72. Patel D.J., Ma J.B., Yuan Y.R., Ye K., Pei Y., Kuryavyi V., Malinina L., Meister G., Tuschl T. Structural biology of RNA silencing and its functional implications. Cold Spring Harb. Symp. Quant. Biol. 2006;71:81–93. doi: 10.1101/sqb.2006.71.053. [DOI] [PMC free article] [PubMed] [Google Scholar]
  73. Pellino J.L., Jaskiewicz L., Filipowicz W., Sontheimer E.J. ATP modulates siRNA interactions with an endogenous human Dicer complex. RNA. 2005;11:1719–1724. doi: 10.1261/rna.2102805. [DOI] [PMC free article] [PubMed] [Google Scholar]
  74. Peters L., Meister G. Argonaute proteins: Mediators of RNA silencing. Mol. Cell. 2007;26:611–623. doi: 10.1016/j.molcel.2007.05.001. [DOI] [PubMed] [Google Scholar]
  75. Petersen C.P., Bordeleau M.E., Pelletier J., Sharp P.A. Short RNAs repress translation after initiation in mammalian cells. Mol. Cell. 2006;21:533–542. doi: 10.1016/j.molcel.2006.01.031. [DOI] [PubMed] [Google Scholar]
  76. Pham J.W., Pellino J.L., Lee Y.S., Carthew R.W., Sontheimer E.J. A Dicer-2-dependent 80s complex cleaves targeted mRNAs during RNAi in Drosophila . Cell. 2004;117:83–94. doi: 10.1016/s0092-8674(04)00258-2. [DOI] [PubMed] [Google Scholar]
  77. Pillai R.S., Bhattacharyya S.N., Artus C.G., Zoller T., Cougot N., Basyuk E., Bertrand E., Filipowicz W. Inhibition of translational initiation by Let-7 microRNA in human cells. Science. 2005;309:1573–1576. doi: 10.1126/science.1115079. [DOI] [PubMed] [Google Scholar]
  78. Rajewsky N. microRNA target predictions in animals. Nat. Genet. 2006;38(Suppl):S8–S13. doi: 10.1038/ng1798. [DOI] [PubMed] [Google Scholar]
  79. Rehwinkel J., Behm-Ansmant I., Gatfield D., Izaurralde E. A crucial role for GW182 and the DCP1:DCP2 decapping complex in miRNA-mediated gene silencing. RNA. 2005;11:1640–1647. doi: 10.1261/rna.2191905. [DOI] [PMC free article] [PubMed] [Google Scholar]
  80. Rehwinkel J., Natalin P., Stark A., Brennecke J., Cohen S.M., Izaurralde E. Genome-wide analysis of mRNAs regulated by Drosha and Argonaute proteins in Drosophila melanogaster . Mol. Cell. Biol. 2006;26:2965–2975. doi: 10.1128/MCB.26.8.2965-2975.2006. [DOI] [PMC free article] [PubMed] [Google Scholar]
  81. Robb G.B., Rana T.M. RNA helicase A interacts with RISC in human cells and functions in RISC loading. Mol. Cell. 2007;26:523–537. doi: 10.1016/j.molcel.2007.04.016. [DOI] [PubMed] [Google Scholar]
  82. Rodriguez A., Vigorito E., Clare S., Warren M.V., Couttet P., Soond D.R., van Dongen S., Grocock R.J., Das P.P., Miska E.A., et al. Requirement of bic/microRNA-155 for normal immune function. Science. 2007;316:608–611. doi: 10.1126/science.1139253. [DOI] [PMC free article] [PubMed] [Google Scholar]
  83. Schmitter D., Filkowski J., Sewer A., Pillai R.S., Oakeley E.J., Zavolan M., Svoboda P., Filipowicz W. Effects of Dicer and Argonaute down-regulation on mRNA levels in human HEK293 cells. Nucleic Acids Res. 2006;34:4801–4815. doi: 10.1093/nar/gkl646. [DOI] [PMC free article] [PubMed] [Google Scholar]
  84. Schwarz D.S., Hutvágner G., Du T., Xu Z., Aronin N., Zamore P.D. Asymmetry in the assembly of the RNAi enzyme complex. Cell. 2003;115:199–208. doi: 10.1016/s0092-8674(03)00759-1. [DOI] [PubMed] [Google Scholar]
  85. Seggerson K., Tang L., Moss E.G. Two genetic circuits repress the Caenorhabditis elegans heterochronic gene lin-28 after translation initiation. Dev. Biol. 2002;243:215–225. doi: 10.1006/dbio.2001.0563. [DOI] [PubMed] [Google Scholar]
  86. Selbach M., Schwanhausser B., Thierfelder N., Fang Z., Khanin R., Rajewsky N. Widespread changes in protein synthesis induced by microRNAs. Nature. 2008;455:58–63. doi: 10.1038/nature07228. [DOI] [PubMed] [Google Scholar]
  87. Song J.J., Liu J., Tolia N.H., Schneiderman J., Smith S.K., Martienssen R.A., Hannon G.J., Joshua-Tor L. The crystal structure of the Argonaute2 PAZ domain reveals an RNA binding motif in RNAi effector complexes. Nat. Struct. Biol. 2003;10:1026–1032. doi: 10.1038/nsb1016. [DOI] [PubMed] [Google Scholar]
  88. Song J.J., Smith S.K., Hannon G.J., Joshua-Tor L. Crystal structure of Argonaute and its implications for RISC slicer activity. Science. 2004;305:1434–1437. doi: 10.1126/science.1102514. [DOI] [PubMed] [Google Scholar]
  89. Thermann R., Hentze M.W. Drosophila miR2 induces pseudo-polysomes and inhibits translation initiation. Nature. 2007;447:875–878. doi: 10.1038/nature05878. [DOI] [PubMed] [Google Scholar]
  90. Tolia N.H., Joshua-Tor L. Slicer and the argonautes. Nat. Chem. Biol. 2007;3:36–43. doi: 10.1038/nchembio848. [DOI] [PubMed] [Google Scholar]
  91. Vasudevan S., Steitz J.A. AU-rich-element-mediated upregulation of translation by FXR1 and Argonaute 2. Cell. 2007;128:1105–1118. doi: 10.1016/j.cell.2007.01.038. [DOI] [PMC free article] [PubMed] [Google Scholar]
  92. Vasudevan S., Tong Y., Steitz J.A. Switching from repression to activation: MicroRNAs can up-regulate translation. Science. 2007;318:1931–1934. doi: 10.1126/science.1149460. [DOI] [PubMed] [Google Scholar]
  93. Vinther J., Hedegaard M.M., Gardner P.P., Andersen J.S., Arctander P. Identification of miRNA targets with stable isotope labeling by amino acids in cell culture. Nucleic Acids Res. 2006;34:e107. doi: 10.1093/nar/gkl590. [DOI] [PMC free article] [PubMed] [Google Scholar]
  94. Wakiyama M., Takimoto K., Ohara O., Yokoyama S. Let-7 microRNA-mediated mRNA deadenylation and translational repression in a mammalian cell-free system. Genes & Dev. 2007;21:1857–1862. doi: 10.1101/gad.1566707. [DOI] [PMC free article] [PubMed] [Google Scholar]
  95. Wightman B., Ha I., Ruvkun G. Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans . Cell. 1993;75:855–862. doi: 10.1016/0092-8674(93)90530-4. [DOI] [PubMed] [Google Scholar]
  96. Wu L., Fan J., Belasco J.G. MicroRNAs direct rapid deadenylation of mRNA. Proc. Natl. Acad. Sci. 2006;103:4034–4039. doi: 10.1073/pnas.0510928103. [DOI] [PMC free article] [PubMed] [Google Scholar]
  97. Wu Z., Irizarry R.A., Gentleman R., Martinez-Murillo F., Spencer F. A model-based background adjustment for oligonucleotide expression arrays. J. Am. Stat. Assoc. 2004;99:909–917. [Google Scholar]
  98. Yan K.S., Yan S., Farooq A., Han A., Zeng L., Zhou M.M. Structure and conserved RNA binding of the PAZ domain. Nature. 2003;426:468–474. doi: 10.1038/nature02129. [DOI] [PubMed] [Google Scholar]
  99. Yekta S., Shih I.H., Bartel D.P. MicroRNA-directed cleavage of HOXB8 mRNA. Science. 2004;304:594–596. doi: 10.1126/science.1097434. [DOI] [PubMed] [Google Scholar]
  100. Yi R., Qin Y., Macara I.G., Cullen B.R. Exportin-5 mediates the nuclear export of pre-microRNAs and short hairpin RNAs. Genes & Dev. 2003;17:3011–3016. doi: 10.1101/gad.1158803. [DOI] [PMC free article] [PubMed] [Google Scholar]
  101. Yisraeli J.K. VICKZ proteins: A multi-talented family of regulatory RNA-binding proteins. Biol. Cell. 2005;97:87–96. doi: 10.1042/BC20040151. [DOI] [PubMed] [Google Scholar]
  102. Zeng Y., Wagner E.J., Cullen B.R. Both natural and designed microRNAs can inhibit the expression of cognate mRNAs when expressed in human cells. Mol. Cell. 2002;9:1327–1333. doi: 10.1016/s1097-2765(02)00541-5. [DOI] [PubMed] [Google Scholar]
  103. Zhang L., Ding L., Cheung T.H., Dong M.Q., Chen J., Sewell A.K., Liu X., Yates J.R., 3rd, Han M. Systematic identification of C. elegans miRISC proteins, miRNAs, and mRNA targets by their interactions with GW182 proteins AIN-1 and AIN-2. Mol. Cell. 2007;28:598–613. doi: 10.1016/j.molcel.2007.09.014. [DOI] [PMC free article] [PubMed] [Google Scholar]
  104. Zhao Y., Ransom J.F., Li A., Vedantham V., von Drehle M., Muth A.N., Tsuchihashi T., McManus M.T., Schwartz R.J., Srivastava D. Dysregulation of cardiogenesis, cardiac conduction, and cell cycle in mice lacking miRNA-1-2. Cell. 2007;129:303–317. doi: 10.1016/j.cell.2007.03.030. [DOI] [PubMed] [Google Scholar]

Articles from RNA are provided here courtesy of The RNA Society

RESOURCES