ABSTRACT
Pumilio domain RNA-binding proteins are known mainly as posttranscriptional repressors of gene expression that reduce mRNA translation and stability. Trypanosoma brucei has 11 PUF proteins. We show here that PUF2 is in the cytosol, with roughly the same number of molecules per cell as there are mRNAs. Although PUF2 exhibits a low level of in vivo RNA binding, it is not associated with polysomes. PUF2 also decreased reporter mRNA levels in a tethering assay, consistent with a repressive role. Depletion of PUF2 inhibited growth of bloodstream-form trypanosomes, causing selective loss of mRNAs with long open reading frames and increases in mRNAs with shorter open reading frames. Reexamination of published RNASeq data revealed the same trend in cells depleted of some other proteins. We speculate that these length effects could be caused by inhibition of the elongation phase of transcription or by an influence of translation status or polysomal conformation on mRNA decay.
INTRODUCTION
Pumilio domain proteins are involved in various aspects of RNA metabolism and are found throughout eukaryotic evolution (1–3). They usually contain approximately 8 Puf repeats, each of which binds a single nucleotide with a consensus target sequence of UGUA (4). Pumilio proteins were first discovered as regulators of differentiation in Drosophila melanogaster, and subsequently they have been found to function in differentiation of other animal species (3, 5). Most pumilio proteins studied so far have been found to be cytosolic, but there are also conserved pumilio proteins that are located in the nucleolus and are involved in rRNA processing (1, 6).
The yeast Saccharomyces cerevisiae has 6 cytosolic pumilio domain proteins, none of which is essential: even cells lacking Pufs1, 2, 3, 4, and 5p are still viable, although they have extensive changes in gene expression (7). The yeast cytosolic Puf proteins assist in mRNA localization, suppress translation, and promote mRNA decay. Puf4p and Puf5p bind to adjacent sites in the 3′-untranslated region (3′-UTR) of the HO mRNA, repressing its expression and causing deadenylation (8); they seem to act via interaction with the eukaryotic initiation factor 4E (eIF4E)-binding proteins Caf20p and Eap1p (9). Puf1p acts in conjunction with Puf3p (10), which is involved in the localization of mRNAs that encode mitochondrial proteins to the mitochondrial membrane (11). Meanwhile, Puf6p is involved in translational repression and localization of the ASH1 mRNA (12). Pumilio domain proteins in malaria parasites effect stage-specific translational repression (13, 14). In contrast, some pumilio proteins have also been shown to increase expression: for example, in Caenorhabditis elegans olfactory sensory neurons, Fas-binding factor 1 (FBF-1) promotes translation of egl-4 mRNA by binding to its 3′-UTR (15).
African trypanosomes are excellent models for the study of posttranscriptional regulation: transcription is polycistronic, with no apparent control of RNA polymerase II initiation. Individual mRNAs are produced by trans splicing and polyadenylation (16). Two life cycle stages are routinely cultured in the laboratory: the bloodstream form, which is similar to the forms that multiply in mammalian blood and tissue fluids, and the procyclic form, which is similar to the parasites in the tsetse fly midgut. Different mRNA levels—and also regulation during differentiation—are achieved by modulation of mRNA decay (17) and trans-splicing (18–20). Trypanosome mRNA decay pathways are similar to those of other eukaryotes: mRNAs are deadenylated by a CAF1-NOT complex (21), then the mRNA body is attacked from the 3′ end by the exosome (22–24) and in the 5′-3′ direction by XRNA (25). Depletion of the deadenylase CAF1 decreases deadenylation and therefore inhibits degradation of most mRNAs (22, 26); depletion of the core exosome subunit RRP45 decreases 3′-5′ degradation by the exosome (22, 23). Meanwhile, depletion of XRNA stabilizes mRNAs that are normally very unstable (17). Numerous trypanosome RNA-binding proteins have been implicated in the control of mRNA fate (27–29).
Trypanosoma brucei has 11 pumilio domain proteins. Three of them, PUF7, PUF8, and PUF10, are in the nucleolus and have been implicated in aspects of rRNA metabolism (30), although PUF10 also has functions in expression of mRNAs that are synthesized by RNA polymerase I (31). Neither PUF1 (32) nor PUF5 (33) is essential in procyclic forms, and they may also be dispensable in bloodstream forms; their roles are not known. PUF9 binds to, and stabilizes, a few mRNAs in the S phase of the cell cycle (34). We describe here studies of PUF2.
MATERIALS AND METHODS
Trypanosome culture.Culturing and transfection of trypanosomes (Lister 427 strain expressing the tet repressor from pHD1313 [1313 trypanosomes]) were performed as described in reference 35. For the expression of tetracycline-inducible genes, 100 ng ml−1 of tetracycline (Tet) (36) was used. For monitoring the growth, procyclic cells were grown in minimal essential medium-Pros medium at 27°C and were diluted to a density of 5 × 105 cells ml−1. Bloodstream-form trypanosomes were grown in HMI-9 medium at 37°C with 5% CO2 and diluted to 1 × 105 cells ml−1. To monitor growth after RNA interference (RNAi), trypanosome cultures were initially diluted and grown without selecting antibiotics, then diluted to 2 × 105 cells ml−1 for bloodstream forms and 5 × 105 cells ml−1 for procylic forms before tetracycline addition. Thereafter, cells were similarly diluted daily as required. Cells with MKT1 RNAi have been described elsewhere (37).
Plasmid constructs.Plasmids and oligonucleotides are listed in Tables S1A and S1B, respectively, in the supplemental material.
To express V5-tagged PUF2 from the PUF2 locus, 400-bp-long fragments from the 5′-UTR and the open reading frame (ORF) of PUF2 were amplified by using primers CZ3539 (SacII) and CZ3540 (XbaI) and CZ3537 (XhoI) and CZ3538 (ApaI), respectively. The PCR products were ligated into the plasmid Bla/V5 (38). The resulting plasmid, pHD2180, was digested with SacII and ApaI preceding transfection into bloodstream-form and procyclic-form cells. For in situ tagging using eYFP (enhanced yellow fluorescent protein), a 708-bp fragment from the C terminus of PUF2 was amplified using the primers CZ4024 and CZ4025. The resulting PCR fragment was ligated into the KpnI and BamHI sites of p3043 (39) (kind gift from Mark Carrington). The resulting plasmid, pHD2230, was linearized with EcoNI before use in transfecting bloodstream-form trypanosomes.
For tandem affinity purification (TAP) tagging, the PUF2 ORF was amplified using the primers CZ3942 (HindIII) and CZ3943 (HpaI), digested with the corresponding restriction enzymes, and ligated into the linearized pHD918 (24) to generate pHD2236. After linearizing with NotI, pHD2236 was transfected into 1313 bloodstream-form trypanosomes.
RNA interference against PUF2 was done by making a stem-loop RNAi plasmid as described in reference 36. A 432-bp region (bp 399 to 831 of the ORF) of the PUF2 gene was selected for RNAi (the sequence was verified for its specificity by using the RNAit program [http://trypanofan.path.cam.ac.uk/software/RNAit.html]). The stem-loop PUF2 RNAi vector was referred to as pHD2183, and after digestion with NotI the vector was used to transfect 1313 bloodstream- and procyclic-form cells.
For tethering, the PUF2 ORF was PCR amplified using the primers CZ3913 (ApaI) and CZ4025 (BamHI) and ligated into vector pHD2202 (40) to give pHD2304. The plasmid was cotransfected for inducible expression of the fusion protein, λN-PUF2-myc, into cell lines constitutively expressing a chloramphenicol acetyltransferase (CAT) reporter with an actin 3′-UTR (pHD1991) and box B actin 3′-UTR (pHD2277). CAT assays were done as described in reference 41.
Antibody production and affinity purification.The first 942 bp of the PUF2 ORF from the start codon (excluding the PUF domain) were cloned in the pQTEV-based (Addgene, Cambridge, MA) vector pHD1746 (10×His tags), to give pHD2181. Escherichia coli BL21 cells were used to express the recombinant protein. The PUF2 fragment was purified from the soluble fraction under native conditions by using a Ni-nitrilotriacetic acid superflow apparatus (Qiagen, Germany). Antibody against the purified PUF2 protein fragment was raised in a rabbit (Charles River Laboratories, Kisslegg, Germany). Antibodies were affinity purified using His-tagged PUF2 bound to nitrocellulose membrane (42). To check the sensitivity and specificity, 0.1 to 5.0 ng of the recombinant protein was loaded on a 12% SDS-PAGE, and the blot was probed with antibody diluted to 1:10,000. The antibody was able to detect a minimum of 0.2 ng of protein by Western blotting with little cross-reactivity with other trypanosome proteins, and cleavage of the His tag did not affect the signal (data not shown).
Indirect immunofluorescent staining and Western blotting.For immunofluorescent staining, the steps of slide preparation and immunofluorescent staining were done as described in reference 43, with incubation at room temperature with gentle shaking unless otherwise specified. The z-stacks of the cells were obtained and the images were deconvoluted (Wiener algorithm) using the Olympus Cell-R microscope.
Antibodies used for immunofluorescence and/or Western blotting were the following: anti-PUF2 (this study), anti-myc (Santa Cruz Biotechnology), anti-aldolase (44), anti-peroxiredoxin (kind gift from Luisa Krauth-Siegel), anti-SCD6 (kind gift from Mark Carrington), and anti-V5 (AbD SeroTec).
Northern blotting.The 1313 bloodstream-form trypanosomes in the logarithmic growth phase (0.8 × 106 to 1.2 × 106 cells ml−1) were harvested. RNA was prepared using Trifast reagent (Peqlab, Germany). Ten micrograms of total RNA was run on a denaturing formaldehyde agarose gel and then blotted on a nylon membrane (Hybond-N; Amersham, GE Healthcare). Radioactive DNA probes (Prime-IT RmT random primer labeling kit; Stratagene) were used for hybridization, and the signals were detected by using a phosphorimager (FLA 7000; Fujifilm).
Cell fractionation and polysome preparation.Approximately 5 × 108 1313 bloodstream-form cells (density, ≈1 × 106 cells ml−1) were fractionated as described in reference 45. Polysomes were purified from the same number of cells according to the method described in reference 46. The polysomal samples were trichloroacetic acid precipitated and resuspended in 4× Laemmli sample buffer. Blots were incubated with antibodies against PUF2, trypanothione reductase (TR; kind gift from L. Krauth-Siegel's lab), and small ribosomal subunit protein S9.
TAP and coimmunoprecipitation.A total of 1 × 1010 bloodstream-form cells expressing the C-terminally TAP-tagged PUF2 (TAP-PUF2) were harvested after 24 h of induction with 100 ng ml−1 tetracycline. Tandem affinity purification was done as described previously (47). The proteins were separated by SDS-PAGE, and the entire lane from the TAP-PUF2 purification was examined by mass spectrometry. Coimmunoprecipitation was done as described in reference 48. Proteins were detected by Western blotting using antibodies against V5 tag (1:2,000; AbD seroTec), the Myc tag (1:2,000; Santa Cruz Biotechnology), and PUF2 (1:2,000; present study).
Cross-linking and immunoprecipitation.PUF2 cross-linking to RNA followed by immunoprecipitation (CLIP) was performed. The method used was a combination of the CLIP (49) and the PAR-CLIP (photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation) (50, 51) methods with minor changes and without in vivo labeling with 4-thio-uracil (52). Approximately 2 × 109 bloodstream-form trypanosomes expressing V5-PUF2 were used for the pulldown assay. An equal number of 1313 bloodstream-form cells without any tagged protein expression was taken as a negative control. Around 2 × 108 bloodstream-form cells expressing V5-UBP1 (53) were used as a positive control. All bloodstream-form cells were harvested at a density of approximately 1.0 × 106 cells ml−1. All cultures were cross-linked two times with 0.4 J cm−2 254-nm UV light in a Stratalinker apparatus with mixing in between. Immunoprecipitation was done using anti-V5 antibody-coupled agarose beads (Bethyl Laboratories, TX). RNase T1 (10 U μl−1 of the original bead volume; Fermentas, Germany) used to digest the RNA-protein complexes bound to the beads. The remaining RNAs were dephosphorylated, followed by 5′-end labeling with 32P. A 4 to 12% Novex NuPAGE bis-Tris gel (Novex, Germany) was used to separate the RNA-protein complexes. Aliquots of around 5% of the lysate and unbound and bound fractions were used for Western blotting to check the efficiency of the immunoprecipitation.
Transcriptome analysis.For the RNAi experiment, trypanosomes were grown for 1 day in the presence of tetracycline, and poly(A)+ mRNA was sheared and prepared for sequencing according to the standard Illumina protocol. One PUF2 RNAi batch was sequenced on the GAII and the second with HiSeq. Data analysis, including motif searches and enrichment calculations, were performed as described previously (22, 54). All data were log transformed before the correlation coefficient calculations, in order to reduce the influence of outliers and correct for a nonnormal distribution of the data.
RESULTS
PUF2 is a cytoplasmic protein that is essential for normal growth of bloodstream-form trypanosomes.PUF2 expression was examined in two ways. First, a sequence encoding a V5 tag was integrated in frame to the 5′ end of the open reading frame. Later, we expressed the N-terminal 314 amino acids (without the PUF domains) as a His-tagged fusion in E. coli and raised an antibody. Figure 1A (lane 1) shows that the antibody recognized a protein from bloodstream-form trypanosomes that migrated slightly faster than the 100-kDa marker. In cells in which one PUF2 gene was N-terminally tagged with a sequence encoding the V5 epitope, the total level of PUF2 remained similar (Fig. 1A, lane 2). (The tag is too short to affect mobility.) Cells bearing a plasmid carrying a version with a tandem affinity tag showed the appearance of an additional, slower-migrating band after tetracycline addition (Fig. 1A, lanes 3 to 6). The specificity of the antibody was demonstrated by the absence of a signal after RNAi targeting of PUF2 (Fig. 1A, compare lanes 7 and 8). To obtain an approximate estimate of the abundance of PUF2 in bloodstream-form cells, we loaded different amounts of cell extract and the recombinant protein onto gels and compared the signals by Western blotting. In several experiments, 5 ng (1.5 × 104 molecules) of the recombinant protein fragment gave a signal of similar intensity to that from 5 × 106 cells (an example is shown in Fig. 1A), which means that there are probably 1 × 104 to 2 ×104 molecules of PUF2 per bloodstream-form cell or 0.5 to 1 molecules of PUF2 per mRNA (55). Unfortunately, we were unable to obtain any soluble PUF2 for in vitro studies.
Expression and localization of PUF2. (A) Amounts of PUF2 in cells. A 1:10,000 dilution of the antibody was used to probe 5 × 106 bloodstream-form wild-type cells, along with cells overexpressing TAP-tagged PUF2 and PUF2 RNAi cells. Different amounts of purified PUF2 fragment were used: 5.0 to 200 ng was loaded to check the sensitivity and specificity of the antibody. The antibody was also used to estimate the expression level of the protein (in terms of number of molecules per cell). (B) Intracellular location of V5-PUF2. Bloodstream-form cells expressing V5-PUF2 were used. Anti-V5 antibody was used to detect the localization of V5-PUF2 protein, while peroxiredoxin was used as a control for cytoplasmic localization. Cells without a V5-tagged protein were included as negative control.
An immunofluorescence assay using the bloodstream-form cells with constitutive expression of V5-PUF2 showed that the tagged protein was localized in discrete foci in the cytoplasm (Fig. 1B). The cytoplasmic localization was confirmed in a immunofluorescence using the affinity-purified PUF2 antibody (see Fig. S1 in the supplemental material) and PUF2 tagged at the C terminus with eYFP (see Fig. S2A and B in the supplemental material). This was consistent with a role in mRNA metabolism. The pattern for PUF2-YFP was clearly distinct from that seen with a soluble protein, peroxiredoxin (see Fig. S2C), but neither V5-PUF2 nor PUF2-YFP formed granules after heat shock and neither colocalized perceptibly with the stress granule marker SCD6 (see Fig. S2D). The levels of PUF2 in bloodstream forms were not affected by 1-h incubations at 16°C or 41°C or by treatment with the protease inhibitor MG132 (10 μg/ml) (data not shown).
To examine the role of PUF2, we decreased its expression by RNAi. PUF2 knock-down lead to a growth defect in the bloodstream form of the parasite after 2 days of induction. (Fig. 2; see also Fig. S3 in the supplemental material). No effect on procyclic growth was seen. These results are consistent with those previously seen in a high-throughput RNAi screen (56).
PUF2 RNAi in bloodstream-form cells. (A) RNAi against PUF2 in bloodstream-form cells (without any tagged protein) led to a growth defect. The cumulative cell numbers are plotted. Results are means and standard deviations for 10 independent cell lines. (B) Western blot probed with a 1:10,000 dilution of PUF2 antibody for one of the cell lines, showing that RNAi led to a decrease in the PUF2 protein level. Aldolase was used as a loading control (1:50,000 dilution of the antibody).
To find out whether PUF2 interacts with other proteins, we subjected a tagged version to tandem affinity purification. In an initial experiment, we examined only those bands that were absent in the control (TAP-only) purification. In the second experiment, all copurifying proteins were analyzed. To find out whether proteins were specifically associated with PUF2, we compared the new results with those for other purifications (unpublished results) (57). Numerous proteins copurified specifically with PUF2, but none of them was found in both purifications (see Table S2 in the supplemental material). Of those found in the second experiment, which was likely to yield more comprehensive results, several were associated with RNA metabolism. They included several other proteins with RNA-binding domains: ZC3H41, RBP20, RBP33, ALBA4, and another protein with two weak RNA recognition motifs. It was tempting to speculate that the Nudix hydrolase encoded by Tb927.11.9810 could be the decapping enzyme, but this sequence matches S. cerevisiae Ysa1p (E-values of 7.5e−19 and 3.3e−18 for BLASTP and reciprocal BLASTP, respectively). Results from mutagenesis experiments indicate that Ysa1p is involved in ADP ribose and NAD+ metabolism (58). Its closest human homologue is Nudt5, which lacks decapping activity (59).
Some PUF2 is bound to RNAs in vivo.In order to verify the RNA-binding ability of PUF2, a modified in vivo CLIP analysis was performed. The cells were irradiated at 254 nm, which leads to covalent bond formation between RNA and bound proteins. After immunoprecipitation, all unbound RNA except that protected by the protein was digested. Remaining RNAs were end labeled with 32P, and RNA-protein complexes were detected after denaturing gel electrophoresis. As a negative control, we took extracts from cells with no V5-tagged proteins and immunoprecipitated the extracts by using anti-V5 antibody-coated beads. A strong background smear was seen (Fig. 3A). As a positive control, we used extracts from cells expressing V5-UBP1 (53, 60). The control Western blot showed that most of the V5-UBP1 was attached to the beads (Fig. 3B), and the autoradiogram revealed a clear band of 32P-labeled material running at about 40 kDa that was absent from the negative control (Fig. 3A). Using cells expressing V5-PUF2, the anti-V5 immunoprecipitation was apparently not quite as efficient as for V5-UBP1 (Fig. 3B), but there was nevertheless a clear signal at the expected position of 100 kDa that was absent in both the negative and positive controls. Thus, PUF2 is clearly able to bind mRNAs.
PUF2 is bound to some RNA in vivo, as shown via a CLIP assay of V5-PUF2 and V5-UBP1 from bloodstream-form trypanosomes. Wild-type bloodstream-form cells were used as a negative control. Cells were UV irradiated, and extracts were subjected to anti-V5 immunoprecipitation and limiting RNA digestion. The RNA stubs were labeled with 32P, and proteins were separated by denaturing gel electrophoresis. (A) Autoradiogram revealing 32P-labeled RNA cross-linked to the V5-tagged proteins. The unique bands corresponding to the RNA associated with V5-PUF2 (*) and V5-UBP1 (**) are indicated. Results are for 1 × 109 V5-PUF2 or wild-type cells and 2 × 108 V5-UBP1 cells. (B) Western blot showing the efficiency of the immunoprecipitations for the experiment shown in panel A. Blots were probed with anti-V5 antibody. V5-PUF2 (*) and V5-UBP1 (**) are indicated alongside the blot. I, input; Ft, flowthrough; E, eluate.
To determine the nature of the bound mRNAs, we first tried to amplify RNA from the CLIP experiment, but we were unable to obtain sufficient products. Therefore, we repeated the immunoprecipitation of V5-tagged PUF2 and sent the RNA for sequencing without prior RNase treatment or gel electrophoresis. Cells expressing Myc-tagged and V5-tagged UBP1 served as controls, since we expected UBP1 to bind to different sets of mRNAs. After UV cross-linking and immunoprecipitation with appropriate antibody-coupled beads, the bound RNAs were purified and sequenced. For each cell line, RNA in the unbound fraction served as the control. Results are shown in Table S3 of the supplemental material. Although RNAs were specifically selected by the procedures, there was almost no difference between the mRNAs selected by V5-PUF2 and those selected by V5-UBP1. Either these two proteins have similar specificities or the binding is independent of V5-tagged protein. Surveys of the selected 3′-UTRs using MEME suggested enrichment of the sequence AUAUAU (E-value, 6.3e−1,220), and Motifsampler revealed UUGUUUGU (consensus score, 1.76; information content, 1.68).
PUF2 is mostly not associated with polysomes and can decrease levels of bound mRNA.To find out whether PUF2 was associated with polysomes in bloodstream-form cells, we subjected extracts to sucrose gradient centrifugation (Fig. 4A) and examined protein distributions. Ribosomal protein S9 (from the small ribosomal subunit [positive control]) was enriched in the polysomal fractions (Fig. 4B and C). Although PUF2 signal was also observed in the polysomal fractions, the distribution was similar to that of the negative control, TR (Fig. 4B and C). Therefore, only a small proportion, if any, of the PUF2 pool is associated with polysomes.
PUF2 is mainly not associated with polysomes in bloodstream-form trypanosomes. (A) Polyribosomes were separated on a sucrose gradient, and the absorbance at 260 nM was monitored (relative values are graphed). (B) PUF2 was detected on Western blots when using the anti-PUF2 antibody. TR was used as a negative control, and ribosomal protein S9 was the positive control. The amounts of extract loaded (in terms of relative cell equivalents) are shown beneath the blot. C, is unfractionated extract (1/100 of input). (C) The bars show relative amounts of PUF2, TR, and S9 in the gradient. Quantitated Western blot signals were corrected for loading for plotting, such that each bar represents the same number of cell equivalents.
We next assayed the ability of PUF2 to increase or decrease the abundance and translation of associated mRNA. To do this, we used a bloodstream-form cell line that expresses a CAT mRNA bearing 5 copies of the boxB recognition sequence. We also inducibly expressed PUF2 with the λN peptide at the N terminus and a myc tag on its C terminus (40). The λN peptide binds to the boxB sequence and therefore “tethers” any protein to which it is fused to the mRNA (61, 62) (Fig. 5A). A bloodstream-form cell line that included the reporter construct without the λN recognition motif (no boxB) was used as a negative control. Soon after the cell lines were obtained, expression of the λN-PUF2-myc protein was detectable in the absence of tetracycline but clearly increased after tetracycline addition (Fig. 5B). Induction of λN-PUF2-myc expression did not affect the level of CAT protein or RNA if there was no boxB present on the RNA but decreased levels of both if the mRNA had the boxB sequences (Fig. 5C). We later reexamined the cell lines and found that expression of the λN-PUF2-myc protein was now easily detected in the absence of tetracycline (Fig. 5D). We therefore compared CAT expression with that for lines carrying no inducible λN plasmid. Again, the presence of λN-PUF2-myc had no effect in the absence of boxB (Fig. 5E); with boxB, the λN-PUF2-myc caused a decrease in CAT activity to 30%. Previous results showed that tethering of other, neutral proteins has no significant effect on reporter expression (37, 54).
Tethering of PUF2 leads to a decrease in reporter the mRNA level. (A) Schematic representation of the RNA tethering strategy. The left panel shows that the fusion protein binds to the reporter mRNA through the interaction between λN (in the fusion protein) and the stem-loop structure (boxB sequence, in the reporter). The reporter has five stem-loop structures, but for simplicity only one is shown in the figure. When the boxB is not present in the mRNA (3′-UTR), the fusion protein is unable to bind. (B, upper two panels) Expression levels of the λN-PUF2-myc fusion proteins, measured by Western blotting with an anti-myc antibody. Tetracycline was included in the cultures for 24 h where indicated. Aldolase was used as a loading control. (Lower two panels) A Northern blot showing expression of CAT mRNA in the same experiment, with methylene blue-stained rRNA as a loading control. (C) Expression of CAT activity in cells with or without expression of λN-PUF2-myc, after tetracycline induction for 1 day. Results are from three different experiments, with standard deviations shown. A typical λN-PUF2-myc expression result for this experiment is shown in panel B. (D) Expression of λN-PUF2-myc fusion proteins, in the same cell lines as used for panel B (the details are described above in the panel B legend), but after storage in liquid nitrogen. (E) Expression of CAT activity and mRNA in cells with expression of λN-PUF2-myc, with tetracycline induction for 1 day. A typical λN-PUF2-myc expression result for this experiment is shown in panel D. Results are from three different experiments, and standard deviations are shown.
Although it is clearly possible that the artificial tethering, or the presence of the N- and C-terminal tags, could affect PUF2 function, the tethering result suggests that PUF2 may decrease the abundance of mRNAs to which it is bound. The decrease in CAT protein was similar to the decrease in the RNA, suggesting that the CAT protein was reduced because there was less CAT RNA. We cannot, however, rule out the possibility that PUF2 first inhibits translation; in this case, the loss of binding of the mRNA by ribosomes or translation factors could lead to destruction of the mRNA as a secondary effect. Both mechanisms would be consistent with the fact that PUF2 is mainly not associated with polysomes. If PUF2 binding leads to RNA destruction, the bound RNA will not be translated anymore, and if PUF2 inhibits translation initiation, the RNA will not be associated with polysomes.
PUF2 depletion results in specific increases and decreases in mRNAs.To find out whether PUF2 expression influences steady-state mRNA levels, we examined the transcriptomes of PUF2-depleted cells (see Table S4 in the supplemental material). The experiment was done in duplicate. We found 65 open reading frames that were more than 2-fold increased in both experiments, and 107 that were reproducibly decreased. The proteins encoded by the increased mRNAs (Fig. 6A) were preferentially located in the mitochondrion (28%, P value 4.2e−5); moreover, five of the proteins that were of unknown function and location had strongly predicted mitochondrial targeting signals. Meanwhile, the downregulated genes (Fig. 6B) were enriched in those encoding cytoskeletal proteins, including several that are known to be in the flagellum (19%; P = 1e−4).
Functions of genes regulated by PUF2. Genes showing at least a 2-fold increase or decrease in RNA expression after PUF2 RNAi were counted according to protein function. The likelihoods of having these proportions by chance are indicated for those categories that were significantly enriched. (A) Genes with increased expression; (B) genes with decreased expression.
In RNASeq experiments, the amounts of RNAs are relative. That means that by default, a large increase in some mRNAs leads to apparent decreases in others, and vice versa. We therefore used Northern blotting, with stable RNA standards, to examine the decreases in mRNAs coding for paraflagellar rod proteins 1 and 2 (see Fig. S4A and B in the supplemental material), two intraflagellar transport proteins (see Fig. S4C and D), CCCH zinc finger protein ZC3H41 (see Fig. S4E), RAD50 (see Fig. S4F), the RNA recognition motif protein RBP20 (see Fig. S4G), and cyclin 6 (see Fig. S4H). In most cases, the results were confirmed, although the effects were modest. An increase in the zinc finger protein ZC3H20 mRNA was observed (see Fig. S4I), but two other RNA-binding protein mRNAs that supposedly increased were not much affected (see Fig. S4J and K). The results showed that mRNAs both increase and decrease after PUF2 RNAi.
The protein coding region or mRNA length affects the response to RNAi.While examining the results, we noticed that the proteins encoded by the mRNAs that increased were all rather large, while those that decreased seemed small. For example, mRNAs encoding dynein heavy chains were increased, while those encoding the intermediate and light chains were unaffected. We therefore plotted the effect of PUF2 RNAi against the coding sequence length. Remarkably, there was a strong negative correlation: short open reading frame mRNAs increased, and long ones decreased (Fig. 7A; see also Table S4 in the supplemental material). The apparent differential effects of PUF2 RNAi on mRNAs for mitochondrial and cytoskeletal proteins could be artifactual, since the former have shorter protein coding regions than the latter (see Fig. S5A). In an attempt to find mRNAs whose regulation was not just related to open reading frame length, we calculated a correction factor (see Fig. S5B and C and Table S3, sheet 3, in the supplemental material). A total of 75 mRNAs were reproducibly at least 2-fold increased after the length correction, of which 18% encoded mitochondrial proteins. For the 177 reproducibly decreased mRNAs, 8 (4.5%) were cytoskeletal and 9 (5%) were mitochondrial. The expression differences are small, however, and could reflect stochastic variations in the length effect.
Correlations between RNAi effects and coding region lengths. In all cases except for the results shown in panel B, the effect of RNAi is plotted on the vertical axis and the coding region length is on the horizontal axis, using logarithmic scales. The effects of RNAi targeting of different gene products are plotted as follows: (level of mRNA after RNAi)/(level in wild-type cells). Correlation coefficients (R) are shown for the log-transformed data. Linear regression lines (plotted using Microsoft Excel on log2-transformed data) are included where relevant; The formulas were also for log2-transformed data. The numbers in parentheses are for the processed mRNA length, as annotated in TriTypDB. (A) PUF2 RNAi effect versus CDS length. (B) Comparison of PUF2 and MKT1 RNAi effects. (C) CTR9 RNAi effect versus CDS length. (D) ZC3H11 RNAi effect versus CDS length. (E) CAF1 RNAi effect versus CDS length. (F) XRNA RNAi effect versus CDS length.
We wondered whether protein coding region length would always influence expression in cells experiencing growth inhibition, so we compared our results with RNASeq or microarray data from other cell lines with RNAi that targeted expression of essential proteins (see Table S4, sheet 4, in the supplemental material). MKT1 is partially associated with polysomes and may be required for translation (37). We were startled to find that there was a correlation between the effects of RNAi targeting PUF2 and those for MKT1 (Fig. 7B) and that the MKT1 RNAi effect also showed an influence of coding region length (see Fig. S6A in the supplemental material). Intriguingly, RNAi targeting CTR9, which might be involved in transcription elongation (57), showed an even stronger effect (Fig. 7C). There was no correlation with 3′-UTR length (R2 value always under 0.12). ZC3H11 is required mainly for stabilization of chaperone transcripts in procyclic forms, but RNAi causes rapid death of bloodstream forms (54), while RBP10 levels correlate with expression of a subset of bloodstream-form-specific mRNAs (52). The effects of depleting ZC3H11 showed only a minor length influence (Fig. 7D), while for RBP10 there was none at all (see Fig. S6B). Neither stumpy forms (63) (see Fig. S6C) nor bloodstream forms grown to high density (63, 64) (data not shown) showed preferential loss of long coding sequence transcripts. Thus, the length effect is not simply a result of growth inhibition.
We also examined data obtained after depletion of mRNA degradation enzymes. After depletion of the CAF1 deadenylase or the core exosome subunit RRP45 (22), the steady-state levels of short open reading frame mRNAs were preferentially increased in the same way as for PUF2 (Fig. 7E; see also Fig. S6C in the supplemental material). In clear contrast, after depletion of the 5′-3′ exoribonuclease XRNA (17), longer-coding-region mRNAs were preferentially retained (Fig. 7E). We determined whether these effects depended on the coding region or the overall mRNA length by repeating the analysis using the annotated “processed mRNA” length. The correlation coefficients are compared in Fig. 7, and scatter plots are provided in Fig. S7 in the supplemental material. The results were similar, which was not really surprising since the mRNA length is strongly influenced by the coding region length. It therefore is not possible to distinguish between the two possibilities. There was no correlation with 3′-UTR length (the lengths used are provided in Table S4 in the supplemental material).
Finally, in an attempt to find PUF2-specific targets that were not related to the 3′-5′ decay pathway, we searched for mRNAs that were increased after PUF2 depletion but not after depletion of PAN2, CAF1, or RRP45. Only 41 genes qualified, of which 39% encoded mitochondrial proteins. However, many of these also showed increased mRNA after RNAi targeting MKT1, RBP10, or CTR9 (see Table S4, sheet 5, in the supplemental material). It is therefore difficult to tell which, if any, of the PUF2 RNAi effects on the transcriptome are direct.
DISCUSSION
Role of PUF2 in bloodstream-form trypanosomes.Trypanosome PUF2 is required for normal growth of bloodstream forms, can bind to RNA, is not associated with polysomes, and can suppress the level of both mRNA and protein in a tethering assay. The number of PUF2 molecules per cell, 10 to 20,000, is of the same order of magnitude as the number of mRNA molecules (55). The results from tethering experiments suggested that PUF2 may target a subset of mRNAs for degradation.
Identification of PUF2 targets through analysis of RNA that coprecipitated with V5-PUF2 yielded results that were difficult to interpret, since the bound mRNAs looked very similar to those selected when we used V5-tagged UBP1. Motif searches in the 3′-UTRs of putative target mRNAs did not reveal any canonical pumilio-binding domains (UGUA) but did select some sequences with U alternating with A or G. The significance of this is dubious, since the specificity of the precipitation was highly questionable. Both the inability to find mRNAs that are specifically associated with PUF2 and the lack of polysome association of PUF2 are consistent with a role for PUF2 in stimulating decay of its targets, since its binding will lead to destruction of the RNA target.
An alternative method to find PUF2 targets was to examine the transcriptome after PUF2 RNAi, in the expectation that the targets would be increased. We initially thought that PUF2 depletion caused specific increases in mRNAs encoding cytoskeletal proteins and decreases in mRNAs encoding mitochondrial proteins. Upon closer scrutiny, though, we observed a remarkable correlation between the effects of the RNAi and the coding sequence length. Long-coding-region mRNAs were depleted, while mRNAs with short coding regions increased in abundance. Some of the losses and increases were confirmed by Northern blotting. Intriguingly, the same phenomenon was seen in cells in which deadenylation pathways (CAF1, CNOT10, PAN2), or the exosome (RRP45), were depleted, and also after RNAi targeting both MKT1 and CTR9. MKT1 and CTR9 do not have the same function as PUF2—indeed, MKT1 increased mRNA levels in a tethering assay (37), the opposite effect to that of PUF2. One possibility is that the length-related effects of PUF2 depletion are indirect, through reductions in mRNAs encoding degradation machinery components. Indeed, PUF2 RNAi reduced the NOT1 mRNA by 70%, but it is difficult to distinguish cause and effect in this case, since the NOT1 coding region is 7 kb long.
One possibility is that trypanosome PUF2 depletion triggers a more general pathway that results in length-dependent changes in mRNA abundance. However, at present, this conclusion is very speculative and the mechanism is unknown.
The influence of coding region length on RNAi responses.The effects of RNAi targeting of 3′ degradation pathways, PUF2, CTR9, and MKT1 on mRNA abundance were all inversely correlated with coding region or mRNA length. Since the mRNA length mostly correlates with the coding region length, we do not know which is the determining factor. We initially suspected that the length effect was an artifact from RNASeq, since similar phenomena have not been observed among trypanosome microarray data. However, several different control data sets were used, over a period of more than 2 years, using both GAII and HiSeq technologies, and the most recent PUF2 experimental samples (with and without tetracycline) were sequenced together in the same run. Also, depletion of XRNA, which inhibits 5′-3′ mRNA degradation, had the opposite effect (an increase in the longer mRNAs), and the effects of ZC3H11 depletion did not correlate with coding sequence length. In all experiments, the RNA was fragmented before library building. It is therefore unclear how a technical artifact of this sort could arise.
In wild-type trypanosomes, there is no relationship between coding region length and RNA decay rate (see Fig. S5E in the supplemental material), and it is difficult to envisage a mechanism by which either polyadenylation or trans-splicing could be affected by the mRNA length. The strongest length-dependent effect, as judged by both correlation coefficients and the slope of the regression line, was from depletion of CTR9, which likely functions as a transcription elongation factor (57); the next strongest were for the exosome subunit RRP45 and the deadenylase CAF1. (We appreciate that variations in the effectiveness of RNAi could influence this.) Therefore, one possibility is that the length dependence is linked to transcription elongation. Under this scenario, a slowing in elongation would lead to a decrease in the abundances of longer mRNAs, and the effects of depleting other proteins would be caused by some sort of feedback inhibition that decreased the transcription rate. It is already known that RNAi against CAF1 does not result in an overall increase in mRNA abundance (22), even though it almost abolishes RNA decay. This suggests some sort of cross talk between factors required for mRNA decay and transcription, as has already been documented for yeast (65–68). Defects in the CAF1-NOT complex might inhibit mRNA turnover sufficiently to reduce the supply of nucleotides; alternatively, more specific feedback mechanisms could exist. XRNA depletion has different consequences, perhaps because it affects only a small subset of low-abundance, unstable mRNAs (17).
An alternative explanation for the length effects is that there are CDS-length-dependent differences in translation initiation or elongation. It has long been known that mRNAs can circularize during translation, through interaction between poly(A)-binding protein and eIF4G (69, 70). The interaction is thought to inhibit mRNA decay, by preventing deadenylation and/or decapping. It is not clear how many mRNAs show this interaction in vivo, or whether it is influenced by mRNA length. An alternative is spontaneous folding of the RNA, independent of the eIF4G–poly(A)-binding protin interaction (71). In yeast, mRNAs with longer coding regions have shorter poly(A) tails than those with shorter coding regions (72), and mRNAs with short coding regions are more dependent on eIF4G for initiation (73). Such differences in polysome conformation might influence the choice of mRNA decay pathway in cis.
The results shown in Fig. 7 revealed the importance of comparing transcriptome results with other available data sets and of considering unexpected correlations. In doing so, we have revealed a novel paradigm in the control of trypanosome gene expression.
ACKNOWLEDGMENTS
We thank Stuart Archer (Victor Chang Cardiac Research Institute, Darlington, Australia) for hints regarding translation of long coding regions. We also thank L. Krauth-Siegel for the anti-peroxiredoxin antibody and Mark Carrington (University of Cambridge) for ant-SCD6.
This work was supported by Deutsche Forschungsgemeinschaft grants Cl112/14-1 and Cl112/17-1 and by the Bundesministerium für Bildung und Forschung (Sysmo grant “The silicon trypanosome”).
FOOTNOTES
- Received 16 January 2014.
- Accepted 20 March 2014.
- Accepted manuscript posted online 28 March 2014.
Supplemental material for this article may be found at http://dx.doi.org/10.1128/EC.00018-14.
- Copyright © 2014, American Society for Microbiology. All Rights Reserved.