Although first discovered in viruses, previous studies have identified operational −1 ribosomal frameshifting (−1 RF) signals in eukaryotic genomic sequences, and suggested a role in mRNA stability. Here, four yeast −1 RF signals are shown to promote significant mRNA destabilization through the nonsense mediated mRNA decay pathway (NMD), and genetic evidence is presented suggesting that they may also operate through the no-go decay pathway (NGD) as well. Yeast EST2 mRNA is highly unstable and contains up to five −1 RF signals. Ablation of the −1 RF signals or of NMD stabilizes this mRNA, and changes in −1 RF efficiency have opposing effects on the steady-state abundance of the EST2 mRNA. These results demonstrate that endogenous −1 RF signals function as mRNA destabilizing elements through at least two molecular pathways in yeast. Consistent with current evolutionary theory, phylogenetic analyses suggest that −1 RF signals are rapidly evolving cis-acting regulatory elements. Identification of high confidence −1 RF signals in ∼10% of genes in all eukaryotic genomes surveyed suggests that −1 RF is a broadly used post-transcriptional regulator of gene expression.
Programmed ribosomal frameshifting (PRF) is has historically been associated with the study of viruses. PRF signals stochastically redirect ribosomes into new reading frames and viral PRF promotes synthesis of C-terminally extended fusion proteins. The most well defined PRF signals direct ribosomes to slip by one nucleotide in the 5′ (−1) direction. −1 PRF signals typically contain three elements: a ‘slippery site’ composed of seven nucleotides (X XXY YYZ, incoming zero-frame indicated by spaces) where shifting occurs; a short spacer sequence and a downstream stimulatory structure, typically an mRNA pseudoknot (1–3). Current models posit that the pseudoknot directs ribosomes to pause with their aminoacyl- (aa-) and peptidyl-tRNAs positioned over the slippery sequence, where re-pairing of the non-wobble bases of both tRNAs with the −1 frame codons occurs (4–7).
It is now clear that PRF is employed by organisms representing every branch in the tree of life, suggesting an ancient and possibly universal mechanism for controlling the expression of actively translated mRNAs (8). The past few years have witnessed several reports describing in silico identification of recoding signals using a variety of computational approaches (9–16). While the methodologies of each study covered a broad range of bioinformatics techniques, the general goal with the exceptions of (9,15) was to first find out-of-frame ORFs followed by the identification of PRF signals in the overlapping region between them. While this can identify new classes of PRF signals, it is based on the assumption that PRF outcomes should mimic those observed in viral genomes and thus cannot identify new operational outcomes of frameshifting.
While ‘outcome-neutral’ approaches using mRNA motifs known to promote efficient PRF cannot identify new classes of frameshift signals, they enable an expansion of our understanding of operational uses for PRF. The seminal study in this field searched the yeast genome for operational −1 ribosomal frameshift (−1 RF) promoting motifs resembling well characterized examples of viral −1 RF signals, identifying ∼260 putative such elements (9). This work was limited by incomplete annotation of the yeast genome and insufficient computational resources available at the time. New bioinformatics tools were subsequently developed and applied using faster and more robust computational platforms. The results showed that: pattern matching approaches coupled with a predictive method for folding RNA sequences provided a dramatic improvement in the results; −1 RF motifs are widespread in the genome of Saccharomyces cerevisiae and many have predicted secondary structures with statistically significant measures of free energy (15). This analysis showed that ∼11% of yeast genes contain at least one high probability −1 RF signal. Furthermore, we demonstrated that nine putative −1 RF signals selected from a variety of S. cerevisiae genes/genome promoted efficient recoding in vivo. More recently, this bioinformatics protocol has been applied to additional genomes. Currently, more than 25 genomes have been analyzed and it appears that 8–12% of genes contain at least one potential −1 RF signal (see PRFdB at http://prfdb.umd.edu/) (17).
A key finding was that the outcome and function of −1 RF differs significantly between the viral and ‘cellular’ contexts. In viruses, PRF controls the stoichiometries of structural versus enzymatic proteins (18). In contrast, ‘cellular’ RF events redirect elongating ribosomes to premature termination codons, suggesting that −1 RF is used to control cellular mRNA abundance and stability through the nonsense-mediated mRNA decay (NMD) pathway. While PRF is required for the production of functional products, in the cellular context RF appears to operate in a different manner. Thus, in the current work, PRF is used to connote frameshift signals whose function is to produce C-terminally extended proteins with novel functions, while RF is used to refer to frameshift signals that operate to direct ribosomes to premature termination codons. A proof-of-principle experiment demonstrated that a viral −1 PRF signal can function as an mRNA destabilizing element and that mRNA destabilization required NMD (19). Here, rapid degradation of a reporter mRNA through NMD is demonstrated for four cellular yeast −1 RF signals. Further, genetic evidence suggests that the presence of the RF-stimulating pseudoknot may promote mRNA destabilization through the no-go decay (NGD) pathway (20). The EST2 gene, encoding the catalytic subunit of telomerase (21), was used to delve deeper into the relationships between −1 RF and mRNA stability. The EST2 mRNA is destabilized by −1 RF primarily via NMD. Ablation of its five −1 RF signals resulted in stabilization of the EST2 mRNA, and an inverse correlation between −1 RF efficiency and EST2 mRNA steady-steady state abundance was observed.
MATERIALS AND METHODS
Strains, genetic manipulations and media
Escherichia coli DH5α was used to amplify plasmid DNA. Transformations of E. coli were performed as described previously using the calcium chloride method (22). Yeast cells were transformed using the alkali cation method (23). Yeast strains used in this study are shown in Supplementary Table S1. Yeast were grown on YPAD and synthetic complete media (H−) (24). yRP2056, yRP2077 were kind gifts from R. Parker. YJB2659 (generously provided by Judith Berman) was sporulated and strains JD1276, JD1281, JD1287 and JD1288 were obtained by tetrad dissection.
Generation of mRNA stability vectors
Dual luciferase and mRNA stability plasmids have been previously described (15). Oligonucleotide primers were purchased from Integrated DNA Technologies (Coralville, IA, USA) and are shown in Supplementary Table S2. Computationally identified putative −1 RF signals were amplified from yeast genomic DNA using PCR using Oligonucleotide primers which terminated in a SalI restriction site at the 5′ and BamHI at the 3′. The zero-frame dual-luciferase reporter plasmid (pJD375) along with the -1 RF signal containing dsDNA fragments were digested using these restriction enzymes and ligated together to generate endogenous -1 RF signal containing dual-luciferase vectors. Oligonucleotide primers were chosen to terminate in KpnI restriction sites and amplify 41 and 30 bases of Renilla and firefly luciferase derived sequences respectively. The resulting amplicons were cloned into the KpnI site 492 bases into the PGK1 open reading frame of the unmodified PGK1 containing vector (pJD741). A premature termination codon vector (pJD828) was generated by cutting the readthrough (pJD753) with BamHI and backfilling with Klenow fragment. Plasmids so generated are described in Supplementary Table S3.
Generation of EST2 open reading frame mutants
Full length EST2 in a centromeric plasmid and the diploid S. cerevisiae EST2 deletion strain were generously provided by the Berman lab and have been previously described (25). Individual mutant strains were obtained by tetrad dissection. Five potentially significant −1 RF signals were identified in the EST2 open reading frame using the Predicted Ribosomal Frameshift Database (17). The wobble bases of five slippery heptamers were mutagenized to synonymous codons by oligonucleotide site-directed mutagenesis using the QuickChange II XL Site-Directed Mutagenesis Kit (Stratagene). Oligonucleotide design and reaction conditions were performed as recommended by the manufacturer with minor modifications. All mutations were confirmed by sequencing. Plasmids so generated are described in Supplementary Table S3.
Steady state and time course RNA blot analyses of PGK1 harboring endogenous −1 RF signals
mRNA stability vectors were transformed into wild-type yeast (JD1158), upf1Δ or upf2Δ (JD1181 or JD1367), xrn1Δ (JD1170), dcp1Δ (JD1122), ski2Δ (JD1345), ski3Δ (JD19) and dom34Δ (JD1363) cells. The EST2 mRNA stability vector (pJD754) was transformed into rpb1-1 (JD977) and rpb1-1/Upf− (JD978) cells and time courses were performed as described previously (26). Total RNA was extracted with acid phenol/chloroform (pH = 4.5) from mid-logarithmic cell cultures (27), or with Trizole© Reagent following the manufacturer’s directions (Invitrogen, Carlsbad, CA, USA). RNA (northern) blotting was performed as previously described (19). Equal amounts of RNA (1, 2 or 4 μg) were separated through 1% agarose-formaldehyde gels. RNA samples were transferred and UV cross linked to Hybond-N-membranes (Amersham). Blots were hybridized with γ[32P] 5′-end-labeled oligonucleotides specific for U3 snoRNA (loading control) and the exogenous Renilla fragment (experimental). Messenger RNAs were identified using a GeneStorm phosphoimager (Bio-Rad) and quantified using QuantifyOne (Bio-Rad). Each experiment was repeated three or more times and averaged to generate graphs. Error bars for calculations including ratios of ratios may be approximated using either of the two following calculations. ‘Average Ratio’ is defined as the ratio of the two calculated values. Valuectrl is the value of the control (the denominator of the ratio) while Stdevctrl is the calculated standard deviation of the control. Similarly, valueexp and stdevexp are the value and standard deviations of the experimental (the numerator of the ratio). Repctrl and repexp are the number of replicates performed for the control and experimental respectively. The error bars in the graphs of ratios of ratios use the approximated standard error.
Quantitative Real Time PCR
Full length EST2 expression vectors (pJD641), EST2 mutant vectors (pJD796) and null plasmids (pJD315) were transformed into WT (JD1281), EST2 deletion (JD1287), UPF2 (JD1288) and EST2/UPF2 (JD1276) deletion strains. Total RNA was extracted with acid phenol/chloroform (pH = 4.5) from mid-logarithmic cell cultures. In parallel, total RNAs were extracted from isogenic rpl3Δ strains expressing wild-type RPL3 (JD1228), the down-frameshifting rpl3-R247A allele (AM-L3R247A), or the up-frameshifting rpl3-W255C/P257S allele (JD1229). To prevent amplification from contaminating cellular DNA, RNA was treated with DNase I before reverse transcription using Turbo DNase (Ambion). cDNA was generated using the Bio-Rad iScript cDNA synthesis kit and used in the LightCycler real-time PCR system. PCR reactions were performed with 2 µl of cDNA in 20-µl reaction mixtures containing ∼10 nM each sense and antisense primer, and 1x LightCycler 480 SYBR Green I Master Mix (Roche). PCR cycles were run as follows: 1 cycle of 95°C for 10 min; 40 cycles of 95°C for 10 s, 54°C for 20 s and 72°C for 20 s. U3 snoRNA was chosen as a reference gene.
The SPR6, EST2, BUB3 and TBF1 orthologs from the genomes of S. paradoxus, S. mikatae, S. bayanus, S. castellii, S. kudriavzevii and S. kluyveri were extracted from the Yeast Gene Order Browser (http://wolfe.gen.tcd.ie/ygob/) (28). Orthologs were identified for all genes. The nucleotide sequences were analyzed for the presence of potential −1 RF signals as previously described (15,17). Results are compiled in Supplementary Table S4.
Cellular −1 RF signals are mRNA destabilizing elements
Four operational yeast cellular −1 RF signals derived from the BUB3, EST2, SPR6 and TBF1 genes were employed to test the hypothesis that −1 RF signals function as mRNA destabilization elements. The slippery heptamers for these −1 RF signals begin at nucleotides 858, 1653, 279 and 1521 of their respective ORFs. These were cloned into a yeast PGK1 reporter gene so that frameshifted ribosomes are directed to PTCs. All inserts were flanked by sequences derived from Renilla and firefly luciferase genes, providing unique exogenous sequences for specific detection of the reporter mRNAs. Two additional PGK1 reporters without −1 RF signals were used as controls: a readthrough reporter encoded a continuous ORF, while a PTC control contained an in-frame UAA termination codon (Figure 1). Reporters were introduced into wild-type yeast cells; their steady state mRNA abundances were determined by RNA blot analysis and normalized to U3 snoRNA controls (Figure 2). A minimum of three independent blots were performed for all experiments. We note that the blots shown in Figure 2 are simplified for the purpose of publication and that the strains are not all isogenic with one another. In contrast, the bar graphs shown in Figure 2 represent data summarized from multiple blots using isogenic strains.
In wild-type cells, all four of the cellular −1 RF signals and the in-frame PTC containing control were less abundant than the PGK1 reporter mRNA (Figure 2A). The decrease in mRNA steady-state abundance varied from ∼100-fold of the readthrough control (EST2) to ∼0.19-fold of wild-type (TBF1). Experiments were also performed in upf1Δ and dom34Δ strains, and the U3 snoRNA-normalized signal intensities were compared among the same signals between wild-type and mutant strains to determine the relative contributions of NMD and NGD on steady-state abundance of the −1 RF signal-containing reporters (Figure 2B and C). The PTC containing mRNA was only affected through the NMD pathway: 28-fold increased abundance in upf1Δ cells relative to wild-type cells, but no change in dom34Δ cells. The TBF1 −1 RF signal similarly affected the reporter signal only through NMD (∼4-fold). In contrast, the steady-state abundance of the EST2 and BUB3 −1 RF signal-containing reporter mRNAs were increased in both the upf1Δ and dom34Δ mutants: the EST2 signal was 35-fold less effective in decreasing mRNA abundance in upf1Δ cells when compared to the WT strain and ∼14-fold less effective in dom34Δ cells, while the values for the BUB3 signal were ∼7-fold and 8-fold, respectively. The steady-state abundance of the SPR6 −1 RF signal containing reporter mRNA was primarily increased in dom34Δ cells (∼2.0-fold). Deletion of DCP1, XRN1 and SKI3, all of which function downstream of UPF1 or DOM34, also generally increased the abundance of the reporter mRNAs (Figure 2D–F). We note however that, in the case of the dcp1Δ cells, the reporter mRNAs were relatively abundant, most likely because of the presence of the 3′ → 5′ exonuclease activity of the exosome, and/or due to decapping activity contributed by other factors, e.g. by the presence of the L-A virus (29). These results establish that endogenous cellular −1 RF signals can decrease mRNA steady-state abundance in yeast through the NMD pathway. In addition, the data are consistent with the hypothesis that a subset of these signals may also affect mRNA abundance through NGD, although substantiation of this claim requires further studies, e.g. to monitor the abundance of the endonucleolytic cleavage products and mRNA stability assays in NGD− strains.
The EST2 −1 RF signal at nucleotide 1653 is primarily destabilized by −1 RF induced NMD
Figure 2 suggests that −1 RF induced NMD is the major cause of decreased mRNA steady-state abundance by the EST2 −1 RF signal beginning at nucleotide 1653. To confirm this, a series of time course mRNA decay assays were performed employing the PGK1-EST2 −1 RF reporter, the readthrough control, and the PTC containing construct in cells harboring the temperature sensitive rpb1-1 allele of RNA polymerase II. At the zero time point, cells were shifted to the non-permissive temperature (42°C) to arrest transcription of mRNAs, total cellular mRNAs were extracted at 0, 1, 2, 4, 8 and 16 min. subsequent to the temperature shift, and RNA blots were hybridized with the firefly luciferase and U3 snoRNA probes. While the readthrough control was stable in wild-type cells (Figure 3A and D), both the PTC containing control and the reporter containing the EST2 −1 RF signal promoted rapid exponential decay of the reporter mRNA, thus demonstrating that this −1 RF signal can operate as an mRNA destabilizing element (Figure 3B–D). In a parallel experiment using rpb1-1 upf1Δ cells, all of the reporter mRNAs remained stable (Figure 3E–H). The rapid decay kinetic profile of the EST2 −1 RF containing reporter, and its stabilization in NMD-deficient cells are consistent with NMD being the major decay pathway triggered by this element (19). To independently test of this, the A AAA AAT slippery site was partially inactivated by mutating it to G AAG AAC. This silent mutation stabilized the reporter mRNA ∼19-fold compared to the wild-type slippery site in wild-type, i.e. DOM34 cells (Figure 3I). Interestingly, this is less than the 35-fold stabilization in upf1Δ cells. One would expect that, since NMD is dependent of −1 RF, then inactivation of −1 RF should be quantitatively the same as inactivation of NMD. To address this, the steady state abundance of the G AAG AAC slippery site containing PGK1 reporter was assayed in an isogenic dom34Δ strain (Figure 3I, dom34Δ lane). This combination increased the steady-state abundance of the reporter mRNA to near wild-type levels.
Ablation of −1 RF signals increases the steady-state abundance of the yeast EST2 mRNA, and −1 RF efficiency inversely correlates with EST2 mRNA abundance
The EST family of yeast genes is named after their ‘Ever Shortening Telomere’ phenotype (30). EST2 encodes the catalytic subunit of telomerase and the other three EST genes either encode protein subunits of telomerase (EST1 and EST3) or a telomere-associated regulator of telomerase (CDC13/EST4) (31). Telomere elongation occurs in late S phase, although Est2p is associated to varying extents with telomeric chromatin throughout the cell cycle, and telomerase defects result in chromosome instability and rapid senescence (32). The very low abundance EST2 mRNA is stabilized in NMD-deficient cells (33,34).
Computational analyses revealed that EST2 contains four additional high confidence −1 RF signals beginning at positions 72, 1215, 1326 and 1995 (Supplementary Figure S1). The positions of the five predicted −1 RF signals in the EST2 ORF are shown in Figure 4A. Silent protein coding changes were introduced into the slippery sites of all 5 of the −1 RF signals in a full-length EST2 clone expressed from a low copy vector (pEST2ssΔ, Figure 4A). Clones expressing either wild-type EST2 (pEST2) or pEST2ssΔ were introduced into isogenic est2Δ or est2Δ upf1Δ cells, and qRT-PCR analyses were performed. These silent mutations resulted ∼8.5-fold increase in the abundance of the full-length EST2ssΔ mRNA relative to wild-type EST2 mRNA (Figure 4B). Similarly, abrogation of NMD increased the abundance of the wild-type EST2 and EST2ssΔ mRNAs ∼5.8-fold and ∼7.0-fold, respectively. To independently monitor the influence of −1 RF on mRNA abundance, the steady-state abundance of the EST2 mRNA was monitored in isogenic cells expressing up- and down-frameshift promoting alleles of RPL3 (which encodes ribosomal protein L3) by qRT PCR. EST2 mRNA abundances were normalized to U3 snoRNA in cells expressing wild-type RPL3, the rpl3-R247A allele which decreases −1 RF from the L-A frameshift signal to ∼55% of wild-type levels (35), and the rpl3-W255C/P247S allele which increases −1 RF by ∼1.6-fold (36). Relative to wild-type cells, steady-state abundance of the EST2 mRNA was increased by 1.29 ± 0.04 fold in cells expressing rpl3-R247A, and decreased to 0.55 ± 0.02 in cells expressing rpl3-W255C/P247S (Figure 4C). Taken together, these experiments demonstrate that −1 RF induced NMD plays a significant role in destabilizing EST2 mRNA.
Programmed −1 ribosomal frameshifting, but not specific −1 RF signals appears to be conserved and rapidly evolving in budding yeasts
If regulation of gene expression through −1 RF is biologically significant, then −1 RF signals should be present in orthologous mRNAs from other budding yeast species. To address this, the BUB3, EST2, SPR6 and TBF1 orthologs were identified in S. paradoxus, S. mikatae, S. bayanus, S. castellii, S. kudriavzevii and S. kluyveri, and analyzed for potentially significant −1 RF signals as previously described (15). At first glance, these analyses reveal that no single −1 RF signal is completely conserved among the budding yeasts (Supplementary Table S4). However, closer analysis shows that strong candidate −1 RF signals can be identified in the orthologs of all of these genes, although not in every species. For example, as noted above, the S. cerevisiae EST2 mRNA contains five potential −1 RF signals. Similarly, the S. paradoxus ortholog contains five potential −1 RF signals, although none are identical to the S. cerevisiae elements. S. mikatae EST2 appears to harbor two potential −1 RF signals, S. bayanus has three, and S. castelli contains two, and S. kluyveri has three. None were identified in the S. kudriavzevii EST2 ortholog. Turning to SPR6, the S. cerevisiae mRNA contains a second potential −1 RF signal beginning at nucleotide 348 in addition to that identified beginning at nucleotide 279 (see http://cbmgintra.umd.edu/prfdb/index.cgi/detail?id=1755&accession=SGDID:S0000917&slipstart=348). Both the S. paradoxus and S. kudriavzevii SPR6 orthologs contain three potential −1 RF signals, but none were identified in the S. mikatae S. bayanus or S. castelli orthologs. Interestingly, the S. kluyveri SPR6 ortholog contains a slippery site followed by a strong stem-loop structure; while this may or may not constitute a −1 RF signal, it does suggest the presence of a rapidly evolving cis-acting element (see Discussion below). S. cerevisiae BUB3 contains the operational −1 RF signal at nucleotide 858, plus potential −1 RF signals beginning at nucleotides 27 and 732. The orthologous mRNAs in S. paradoxus, S. bayanus, S. castelli and S. kudriavzevii each appear to have one potential −1 RF signal, but none were identified in S. mikatae or S. kluyveri. Lastly, the S. cerevisiae TBF1 mRNA has only the single confirmed −1 RF signal. The S. castelli, S. kudriavzevii and S. kluyveri orthologs contain two each, and the S. mikatae has one. No potential −1 RF signals were identified in either S. paradoxus or S. bayanus. As a control, six S. cerevisiae genes lacking predicted −1 RF signals were selected (PGK1, HHT1, TEF2, MIC14, CMD1 and GRX1), orthologs from the six other yeast species identified, and these were in turn queried for the presence of putative −1 RF signals. These analyses revealed that none of the orthologs of these six genes contain predicted −1 RF signals (Supplementary Table S5, and hyperlinked data therein). The potential evolutionary significance these observations are discussed below.
In a prior proof-of-principle experiment, we utilized the well characterized −1 PRF signal from the yeast L-A dsRNA virus to demonstrate that these elements can generally function as mRNA destabilizing elements through the NMD pathway (19). Subsequently, a bioinformatics approach was used to determine that potential −1 RF signals are widely found in all genomes examined, and that the great majority of these are predicted to direct elongating ribosomes to premature termination codons (15,17). Here, we show that these chromosomally encoded, endogenous −1 RF signals can also function as cis-acting mRNA destabilizing elements, both in the context of a reporter mRNA, and also in one case in a natural context. Further, we demonstrated that −1 RF signals can differentially affect mRNA abundance through the NMD pathway, and the data are also consistent with destabilization through NGD. These are modeled in Figure 5. In support of this idea, the EST2, BUB3 and SPR6 mRNAs were all stabilized in upf1Δ, upd2Δ/nmd2Δ, upf3Δ, dcp1Δ and xrn1Δ cells (37), and the half-lives of these mRNAs were less than the mean in wild-type cells (38). Interestingly, TBF1 is not represented in these databases. In the case of a ribosome shifting reading frame into a PTC, the surveillance complex led by the Upf proteins signals rapid decapping by Dcp1p/Dcp2p, followed by deadenylation and exonucleolytic decay via Xrn1p and the exosome. In parallel, the NGD pathway can be activated by ribosomes that are stalled at strong secondary structures in mRNAs. Stalled ribosomes are freed from mRNAs by Dom34p/Hbs1p, promoting exonucleolytic cleavage at unpaired nucleotides near the pause, thus resulting in two mRNA fragments which become substrates for decapping and exonucleolytic decay [reviewed in (39)]. The findings presented here suggest that cells are not only well equipped to deal with aberrant messages which contain premature termination codons and to clear stalled ribosomes from mRNAs, but have also evolved to capitalize upon these functions to post-transcriptionally regulate gene expression.
The strength of these signals to function as mRNA destabilizing elements should be equal to a combination of (i) their strengths as −1 RF signals and (ii) their abilities to block ribosome progression, i.e. their thermodynamic stability. The EST2 signal is both highly efficient at promoting −1 RF [∼55%, see (15)], and is predicted to be quite stable (∼−27 to −24 kcal/mol depending on the particular folding solution, see http://prfdb.umd.edu/). It is important to note however that the software used to predict mRNA pseudoknots can neither identify base triples, which make major contributions to frameshifting (40–44), let alone calculate their contributions to thermodynamic stability. Regardless, this combination of high frameshifting and thermodynamic stability results in very strong destabilization via NMD (Figure 2B), and perhaps NGD as well (Figure 2C). As discussed previously (19), the exponential decay profile suggests that NMD can occur beyond the ‘pioneer round’ of translation.
In contrast to EST2, the TBF1 signal promoted ∼5% frameshifting (15), but is not predicted to be highly stable (−9.5 kcal/mol). Thus, all of its mRNA destabilization activity was through NMD (compare Figure 2B with C). The thermodynamic stability of the BUB3 signal is predicted to have an intermediate value to EST2 and TBF1 (∼−12 kcal/mol), and hence the potential contribution of NGD to the stability of its reporter was significant. Interestingly, this signal only promoted ∼1% frameshifting (15), yet the contribution of NMD to its destabilization was greater than observed for TBF1. One possible explanation for this apparent discrepancy may stem from the fact that, in order to measure frameshifting, one base had to be deleted from the spacer region between the slippery site and the stimulatory pseudoknot. Changes in the length and composition of this spacer are known to affect rates of −1 PRF (45), and thus the −1 RF values so determined cannot be taken as absolute. In contrast, the reporters used to monitor mRNA stability contained the native sequences. In light of this, it is likely that the native BUB3 −1 RF signal promotes more frameshifting than the TBF1 signal. Lastly, the SPR6 −1 RF signal is predicted to be quite stable (∼−20 kcal/mol), yet promoted very low levels of frameshifting (∼0.5%) (15). Accordingly, destabilization via NMD was negligible for this element, while NGD appeared to be the major contributor.
Beyond the pro forma demonstration that −1 RF signals can decrease cellular mRNA abundance, it is important to begin to understand the biological function of this phenomenon. As a first step in this direction, we showed that silently mutating the slippery sites in 5 predicted −1 RF signals within a full-length clone of EST2 significantly stabilized its encoded mRNA (Figure 5B). Similarly, abrogation of NMD stabilized this message, while changes in −1 RF efficiency inversely correlated with EST2 steady-state mRNA abundance. Est2p is the reverse transcriptase subunit of the telomerase holoenzyme (21). Interestingly, prior studies have demonstrated that this mRNA, along with other mRNAs encoding proteins having telomere-associated functions, are stabilized in NMD− yeast cells (33,37). Analysis of the Programmed Ribosomal Frameshift Database (http://prfdb.umd.edu/) reveals that, along with the other four putative −1 RF signals in the EST2 mRNA, the mRNAs encoding Est1p, Stn1p, Cdc13p and Orc5p, all components or regulators of telomerase that are stabilized in NMD− cells, also contain high confidence −1 RF signals (Supplementary Figure S2). In addition, the EST3 mRNA contains a +1 PRF signal (46). Intriguingly, telomerase is limiting in cells: while a yeast cell contains 64 chromosome ends, there are only ∼29 telomerase molecules per cell, and telomerase is preferentially recruited to short telomeres (47). Additionally, Tbf1p is a telobox containing general regulatory factor that binds to TTAGGG repeats within subtelomeric anti-silencing regions (48). Intriguingly, ablation of NMD (49) or overexpression of single components of telomerase-associated proteins, i.e. the TEL1 RNA, Est2p, Stn1p or Cdc13p resulted in changes in telomere length (47,50,51).
We hypothesize that yeast cells use −1 RF to limit the expression of these proteins in order to maintain the correct stoichiometric balance among telomere associated components. Corollary to this, mutations that alter −1 RF and/or NMD should affect telomere function, and should thus show phenotypic defects similar to those observed in telomerase mutants, e.g. cell cycle progression defects. Indeed, we have isolated numerous such mutants [reviewed in (52)], and have reported that the mof2-1 and mof5-1 mutants, which affect both NMD and −1 RF tend to accumulate large mother-daughter cells, and/or multiply budded cells, typical of G2/M cell cycle defects (53). Similarly, upf1Δ cells have abnormally elongated buds, and decreased telomere lengths (54,55). Intriguingly, mof6-1 mutants, which only affect −1 RF, arrest as large, unbudded cells, typical of M-phase exit defects (53). These observations suggest that stabilization of the mRNAs encoding multiple telomere-associated proteins may have dominant negative effects on telomere homeostasis, and that NMD and −1 RF may regulate different aspects of the cell cycle. Additionally, the central role of Bub3p at the mitotic cell cycle spindle assembly checkpoint and the progeroid phenotypes caused by Bub3p deficiency (56) suggest a more general role for −1 RF in control of cell growth and division. Lastly, the expression of Spr6p during sporulation (57) suggests a role for −1 RF in this developmental process as well. Future studies will dissect the roles of the −1 RF signals in these mRNAs.
Finally, if −1 RF is widely used to regulate gene expression, then it should be well conserved. The major problem associated with attempting a phylogenetic analysis of −1 RF signals is the inherent limitations of the software used to predict them. In short, it is not well enough developed to automatically identify matching motifs. In an attempt to begin to address this issue, the orthologous Bub3p, Est2p, Spr6p and Tbf1p’s in six closely related yeast species were identified, their nucleotide sequences extracted from the Yeast Gene Order Browser (28), and analyzed for the presence of potential −1 RF signals. These analyses revealed that while specific −1 RF signals do not appear to be evolutionarily conserved, −1 RF itself may be relatively well-enough conserved as a mechanism to post-transcriptionally regulate the expression of these genes across many but not all species examined (see Supplementary Table S4). Importantly, the control experiment showed that putative −1 RF signals were not detected in any of the orthologs of six S. cerevisiae genes that themselves were not predicted to contain −1 RF signals. These observations are in agreement with current evolutionary theory based on analyses of 5′ UTR sequences of Drosophila species proposing that rapid rates of mutation in cis-acting regulatory elements drives speciation because they confer very specific effects on gene expression, as opposed to mutations affecting protein structure, the pleiotrophic effects of which impose very high penalties on fitness (58–60). The findings presented in this work, i.e. that while −1 RF signals are not conserved in orthologs from other yeasts, these ORFs contain generically contain −1 RF signals, and that the absence of −1 RF signals in other orthologous ORFs, agrees very well with this theory of molecular evolution.
Supplementary Data are available at NAR Online.
National Institutes of Health (R01GM058859, R21GM068123 to J.D.D.); National Science Foundation (MCB-0084559); NIH Virology Training Grant (T32 AI51967 to A.T.B., partially); University of Maryland College of Chemical and Life Sciences Doctoral Fellowship (to A.T.B., partially). Funding for open access charge: NIH 5R01 GM058859.
Conflict of interest statement. None declared.
We would like to thank Judith Berman and Roy Parker for the gifts of yeast strains and plasmids.