- Research article
- Open Access
Genome-wide polysomal analysis of a yeast strain with mutated ribosomal protein S9
© Pnueli and Arava; licensee BioMed Central Ltd. 2007
- Received: 16 March 2007
- Accepted: 21 August 2007
- Published: 21 August 2007
The yeast ribosomal protein S9 (S9) is located at the entrance tunnel of the mRNA into the ribosome. It is known to play a role in accurate decoding and its bacterial homolog (S4) has recently been shown to be involved in opening RNA duplexes. Here we examined the effects of changing the C terminus of S9, which is rich in acidic amino acids and extends out of the ribosome surface.
We performed a genome-wide analysis to reveal effects at the transcription and translation levels of all yeast genes. While negligible relative changes were observed in steady-state mRNA levels, a significant number of mRNAs appeared to have altered ribosomal density. Notably, 40% of the genes having reliable signals changed their ribosomal association by more than one ribosome. Yet, no general correlations with physical or functional features of the mRNA were observed. Ribosome Density Mapping (RDM) along four of the mRNAs with increased association revealed an increase in ribosomal density towards the end of the coding region for at least two of them. Read-through analysis did not reveal any increase in read-through of a premature stop codon by the mutant strain.
The ribosomal protein rpS9 appears to be involved in the translation of many mRNAs, since altering its C terminus led to a significant change in ribosomal association of many mRNAs. We did not find strong correlations between these changes and several physical features of the mRNA, yet future studies with advanced tools may allow such correlations to be determined. Importantly, our results indicate an accumulation of ribosomes towards the end of the coding regions of some mRNAs. This suggests an involvement of S9 in ribosomal dissociation during translation termination.
- Mutant Strain
- Sucrose Gradient
- Peak Fraction
- Premature Termination Codon
The ribosome is a large RNA-protein complex having multiple roles in the translation process. The RNA moiety in eukaryotes includes four RNA molecules that provide the structural scaffold for the functional sites of the ribosome, and interacts with various proteins of the translation machinery. The protein moiety is composed of more than 70 proteins, located mostly at the ribosome's periphery. These proteins were shown to play various roles, including assistance in ribosome assembly, interaction with translation initiation factors, stabilization of the ribosome structure, and more [1, 2].
The yeast ribosomal protein S9 (sup46) is known to play a role in accurate decoding [3, 4]. Mutagenesis analysis yielded a general (omnipotent) suppressor phenotype, characterized by suppression of termination by all three stop codons, as well as misreading of various codons [5–8]. Some of these mutations are characterized by an increase in the isoelectric point of S9, suggesting a defect in acidic amino acids . Recent structural studies in prokaryotes and eukaryotes revealed that S9 is located at the head of the small ribosomal subunit near the entrance tunnel of the mRNA leading into the decoding center [10–12]. Part of the protein is in close contact with the rRNA, while other regions extend out of the ribosome. The position of the bacterial homolog of S9 (S4) prompted Yusofova et al. to suggest a role in opening the secondary structure of the incoming mRNA . Indeed, it has been recently shown that prokaryotic ribosomes having mutations in S4 (the E. coli homolog of S9) possess a lower ability to unwind RNA duplexes . Thus, S9 might exhibit a helicase activity, which may be related to its role in accurate decoding.
Here we characterized the effects of changing the acidic C terminus of S9. We performed a genome-wide polysomal analysis of a strain expressing S9 mutated in its last 10 amino acids. Of the ~800 genes for which we were able to compare the ribosomal association between wild-type and the mutant strain, 40% had a significant change in ribosomal association. Analysis of several physical and functional features of the affected transcripts did not reveal any feature that correlates with the extent of change. Ribosome density mapping (RDM) analysis for four mRNAs with increased association suggests an accumulation of ribosomes towards the 3' end of the mRNA for at least two of them. These results suggest that S9 is involved in ribosomal dissociation at the termination site.
Mutation of S9
S9 is an essential protein located at the entrance tunnel of the mRNA into the ribosome. It was implicated in functions related to translation elongation such as decoding accuracy and mRNA unwinding. To further explore functions of S9, we performed partial deletions at various regions of interest. Here we present data regarding a mutation of the last 10 amino acids (Fig. 1). These amino acids include multiple aspartic and glutamic acids (D and E) and form a DEAD-like sequence (yet no other motif of classical DEAD helicases appears in S9). They are located outside the region of interaction with the rRNA (the S4 domain) and extend out of the ribosome surface [4, 11, 12].
There are two copies of the S9 gene in yeast that express almost identical proteins. The differences between the encoded proteins are only arginine to lysine substitution in position 40, and four amino acids difference in the C terminus. Both genes were deleted by the mating and sporulation of haploids carrying either RPS9A deletion or RPS9B deletion. A single-copy plasmid carrying the mutated RPS9A coding region with 500 bp upstream and 300 bp downstream was inserted into these cells prior to sporulation, and was selected during the sporulation with the appropriate medium. The last 10 amino acids (188–197) were deleted from the mutated S9 and the deletion procedure led to the addition of three amino acids (Ser, Phe and Trp). Although these added amino acids are uncharged, we cannot exclude the possibility that the observed effects are due to their addition. The deletion resulted in the removal of four amino acids differing between S9A and S9B, and the only remaining difference between S9A and S9B is the substitution of basic amino acid at position 40. We therefore refer to the resulting altered protein as S9 rather than S9A and hereinafter designate it as rps9-Δ188.
Genome-wide analysis of rps9-Δ188
To determine the effects of mutating S9 on a broad spectrum of mRNA sequences, we utilized DNA microarrays carrying ~6,000 known and predicted yeast ORFs. We first describe the results of a transcriptome analysis that characterizes the relative effects on the steady-state mRNAs levels, and then the results of analysis of multiple polysomal fractions aimed at characterizing the effects on ribosomal association.
For transcriptome analysis, RNA was collected from wild-type and rps9-Δ188 strains grown on either rich media (YPD) or minimal media (SD plus Met, Lys, Ade). Each of the RNA samples was labeled with red fluorescent dye (Cy5) and hybridized to a DNA microarray together with a reference RNA sample labeled with a green fluorescent dye (Cy3). About 5,000 genes passed the quality selection procedure (see Methods section) in each of the microarrays (Additional File 1). To determine the genes that were affected by the S9 mutation, we compared the changes in expression of the mutant relative to the wild-type strain in each growth condition. The relative expression of only 52 and 97 genes were changed by more than two-fold in YPD and SD media, respectively (Additional File 1). Many of these genes are Ty elements, a repetitive sequence in the yeast genome that experimentally might represent only a single sequence. The small number of genes whose relative abundance had changed indicates that the mutation in S9 has minor effects on the transcriptome profile of the parental strain. In addition, although S9 was shown to be part of the yeast penta-snRNP complex , we did not observed any significant change in the relative abundance of intron-containing mRNAs. It should be noted that because of the experimental set-up, these changes are only relative and do not report on global changes (e.g., a strong and similar decrease in abundance for all mRNA).
The observation that more genes are affected in SD medium might be related to the higher transcriptional activity of some metabolic pathways under these conditions. Interestingly, the group of reduced genes in the minimal medium included genes involved in nucleic acid metabolism (mainly adenine synthesis), although the medium included adenine, and none of the adenine pathway genes is mutated in this strain. Consistent with this, we observed that the mutant strain exhibits slow growth phenotype in the absence of adenine in the growth media (Additional File 1). This may suggest an indirect role for S9 in the transcription of genes involved in nucleic acid metabolism, probably by affecting the translation of a transcription factor. Indeed, polysomal profile analysis of the Bas1 transcription factor, which is known to regulate the expression of genes encoding enzymes involved in purine, pyrimidine and histidine biosynthetic pathways [17, 18], revealed a decrease in its polysomal association from >80% in the wild-type strain to 57.3% +/- 3.2% in the mutant (data not shown).
To gain insight into the role of S9 during translation, we performed genome-wide analysis of mRNAs separated into multiple polysomal fractions. The rps9-Δ188 strain was grown in optimal conditions (YPD), cells were harvested, and large complexes were separated by velocity sedimentation in sucrose gradient. Figure 2A presents the OD 254 profile of the entire gradient that is indicative of the amounts of polysomal complexes. As can be observed, there is a strong signal in the region of the gradient that includes polysomal complexes, and no significant differences are observed when comparing the OD profile to the profile of a wild-type strain (Fig. 2B). This indicates that the majority of mRNAs in the cells is associated with ribosomes, and is presumably translated.
To characterize changes in translation of particular mRNAs, we collected multiple fractions along the gradient and analyzed the mRNA content in each fraction by DNA microarray. The first fractions were pooled together and included all complexes up to the size of 80S. This fraction was designated "free mRNA" sample, since it contains molecules that are free of fully assembled ribosome. The remaining part of the gradient, which included mRNAs associated with one or more ribosomes, was fractionated to multiple fractions of equal volume. Five in vitro transcribed bacterial RNAs were spiked into each of the fractions at equivalent amounts and were used to normalize for differences in all subsequent steps of the procedure, including RNA purification, labeling with fluorescent dyes, hybridization to DNA microarrays and scanning . This external normalization procedure retains the inherent differences in mRNA amounts at the different fractions, and therefore enables the construction of a sedimentation profile for each gene, which is refractory to global changes in mRNA levels . RNA from each fraction was labeled with red fluorescent dye and hybridized together with an unrelated RNA sample, labeled with green fluorescence as a reference to a DNA microarray. This experimental procedure was repeated three times, and the ratios obtained for all genes are presented in Additional File 2. There are many missing signals for many genes at various fractions. These missing values represent spots that did not pass the quality selection criteria for various reasons, either technical (e.g., bad spotting of the DNA on the microarray) or biological (e.g., no transcripts of that gene in this fraction). For further analysis we selected only those genes that, in general, had only one missing value in their polysomal region of the gradient (see Methods section for details). This selection resulted in ~4,000 genes for which we could assign a "peak fraction", i.e., the fraction having the highest red-to-green signal (Additional File 3). From this group of genes, we filtered in only genes in which the peak fraction repeated itself in two or more experiments (i.e., one or less fraction difference). About 2,800 genes passed this selection criterion, and from the peak fraction and the OD260 of the gradient (Fig. 2) we assigned for each mRNA the average number of ribosomes that it is mostly associated with (Additional File 4).
To determine whether the mutation in S9 led to changes in ribosomal density of specific mRNAs, we compared the findings to published values of an experiment that was done with an isogenic wild-type strain grown under the same medium and subjected to the same extraction procedure. Overall, the wild-type strain and the mutant strain yielded a similar polysomal profile (Fig. 2 and Fig. 1A in ), consistent with no major differences in the experimental procedure. About 800 genes had a reliable ribosomal association value in both the wild-type and the mutant strain (Additional File 5). Not surprisingly, this group includes genes whose transcripts are relatively highly abundant compared to the entire genome of yeast. The median abundance value for this group is 2.2 mRNA copies per cell, while for the entire yeast genome the value is only 1.2 copies per cell. The bias towards abundant mRNA is probably due to the stringent selection procedure, which filters out genes with weak signals or with missing signals from multiple microarrays. This bias probably led to other indirect effects, including the following: 1. An enrichment for genes containing introns (82 of the genes had introns, while in the entire yeast genome the abundance of such genes is <5%). Intron-containing genes have relatively high abundance of mRNAs in yeast. 2. Exclusion of relatively long transcripts (>4,000 nts), as almost all of them are expressed at low levels in yeast. The results and conclusions presented below were also obtained when much lower filtration criteria were employed (data not shown).
List of strains
MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0
BY4741 MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0 rps9a::kanMX4
BY4742 MAT á his3Δ1 leu2Δ0 lys2Δ0 ura3Δ0 rps9b::kanMX4
BY4741 MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0 rps9a::kanMX4::HIS3 rps9b::kanMX4 shm2::kanMX4::hisG +PA183(rps9-Δ188) +PA185
MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0 +PA27 (luciferase-stop)
BY4741 MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0 rps9a::kanMX4::HIS3 rps9b::kanMX4 shm2::kanMX4::hisG +PA183(rps9-Δ188)+ PA27 (luciferase-stop)
BY4741 MAT a his3Δ1 leu2Δ0 met15Δ0 ura3Δ0 rps9a::kanMX4::HIS3 rps9b::kanMX4 shm2::kanMX4::hisG +PA183(rps9-Δ188) +PA28 (luciferase 2 × stop)
List of plasmids
pGal1-luciferase lys445TAG TAA
Various features of groups of affected mRNAs
Number of Genes
ΔG 5UTR (kcal/mol)
GC content (%)
ΔG ORF (kcal/mol)
ORF length (nts)
WT abundance (copies/cell)
Mutant abundancec (copies/cell)
All 815 genes
We also utilized the GO mapper tool from SGD in order to determine if any enrichment for specific functional groups exists in any of the affected groups. For this analysis we combined several groups together (i.e. change of two and more ribosomes to either direction) in order to have groups large enough. The lists of genes that changed by two or more ribosomes was subjected to this analysis (Additional Files 6 and 7). We did not observe any specific enrichment in the process, the function or the component that these genes are assigned to.
S9 mutation leads to increased density at the 3' end
Summary of the three experimental repeats of the RDM analysis
Gene (cleavage site)
YHB1 (cleavage at the middle of the ORF)b
YHB1 (cleavage at the first and last thirds)c
YGR026W (cleavage at the middle of the ORF)d
We analyzed another two mRNAs (Nop15 and Aro7) for differences in sedimentation between their 3' and 5' halves. The results suggest that on these mRNAs as well there is an increased density of the 3' half in the mutant strain, yet the extent of increase is small and therefore its significance is unclear (Additional File 8).
rps9-Δ188 does not confer a read-through phenotype
Genome-wide analysis of the role of ribosomal proteins
Ribosomal proteins are known to play many roles during the translation process [1, 2]. Various methods were used to explore these roles, including sensitivity assays to antibiotics, analyses of reporter mRNAs carrying premature termination codons or rare codons, and structural assays that examined their positions on the ribosomes at different translation states. Here we took a different approach and used the information extracted from the ribosomal association of many different sequences in order to gain insight into the role of a particular ribosomal protein. We were able to compare the ribosomal association of ~800 mRNAs and divide them into several groups of affected mRNAs. Our limited analysis did not reveal any physical or functional feature that correlated with the extent of change (Table 3). Yet, our data might serve as a basis for more thorough analyses, utilizing novel computational tools. Most necessary are algorithms for calculating and aligning RNA structures, which may reveal correlations that we failed to identify. Another current limitation stems from the limited information available for various mRNA features, such as 3' UTRs length, local structure stability and sequence elements. Recently, a large-scale cDNA analysis provided important information about lengths of 5'UTRs in yeast . Although we did not find correlation between the 5' UTRs lengths and changes in ribosomal association, information other than the length (e.g., structural elements) may be involved. Better characterization of these features may assist in defining broader roles of S9.
Changes in transcript levels were shown in many cases to correlate with changes at the translational level, a phenomenon that was termed 'potentiation' . Such a process can occur, for example, if specific proteins are deposited on the transcript as it is synthesized , or if alternative transcription start site is selected . When comparing the changes in transcriptome to the changes in ribosome association of the 815 genes, we did not observe any correlation between the two (data not shown). This could be because our sample is relatively small and has small changes in relative transcript levels. Alternatively, the effects of S9 might only be at the translation step (e.g., during translation termination), and therefore the changes of transcript levels are only indirect consequences of the translational effects. Such indirect secondary effects are likely to be diverse, therefore simple correlation will not be observed.
The analysis presented here focused on those mRNAs having increased density, and based on the few mRNAs that were tested specifically, we propose a role for S9 in termination. However, a significant number of genes appeared to have a similar ribosomal association as in the wild-type, as well as decreased density compared to the wild-type strain. This suggests that the effects on the termination of S9 are conferred only to a subset of mRNAs, for example, mRNAs with a specific sequence or structural motif in their 3' UTR. We were unable to identify a common sequence motif using standard tools (e.g., MEME, ) and are currently attempting to combine structure alignment algorithms to this end. Alternatively, S9 might similarly affect termination on all mRNAs. Yet, due to differences in re-initiation rates between different mRNAs, the effect will appear to be differential and strong enough on some mRNAs to lower ribosomes' reloading. In such a model, the net balance of termination inhibition will be a decrease in ribosomal association for mRNAs with a strong affinity between their 5' and 3' ends, and an increase in density for mRNAs with weak interactions. Testing this model necessitates knowledge about affinities between the ends of different mRNAs or the rates of re-initiation on different mRNAs. Such information is yet to be determined.
Possible involvement of S9 in termination
We have herein focused on the sub-group of mRNAs having increased ribosomal association. For two of the tested mRNAs an increased ribosomal association towards the 3' end of the mRNA was observed; such an increase cannot be excluded for the other two mRNAs. The increased density towards the 3' end of the mRNA might be due to ribosomes reading through the stop codon and accumulating at the 3' UTR. However, if our results from the premature termination codon analysis (Fig. 6) are applicable to a natural stop codon, a likely explanation for the increased density is that ribosomes accumulate upstream to the stop codon. S9 might interact with one of the elongation or release factors that are necessary for efficient termination and ribosome recycling [31–33].
Possible roles of the acidic domain
The region of S9 that was changed includes multiple acidic amino acids. The high-resolution structures of both the prokaryotic and the eukaryotic ribosomes suggest that this region can interact with incoming mRNA or other proteins. There are numerous examples of protein-RNA interactions that are performed through negatively charged residues . These residues interact preferentially with the bases of the mRNA molecule . It is possible that once the ribosome approaches the termination region, the C terminus of S9 interacts with the incoming 3' UTR and thereby allows efficient termination. Another possible role for the acidic domain is through interaction with other proteins. The C terminus might interact with elongation or termination factors, or with other proteins necessary for efficient termination.
We presented here the first genome-wide study of changes in ribosomal association due to a mutation in the ribosomal protein (S9). For analysis we selected a protein that resides at the entrance of the mRNA leading into the ribosomes and which might interact with the incoming message. This analysis revealed groups of mRNAs that changed their ribosomal association compared to the wild-type strain. Limited bioinformatics analysis did not reveal a general correlation between the change in association and physical or functional properties of the mRNA (Table 3). Yet, this data might serve as a template for the development of more advanced tools that will identify a novel sequence or structure elements in the mRNA. We focused here on the group of mRNAs having increased association and found that at least for some the increased association is biased towards the 3' end of the transcript. This suggested a novel function for S9 in the termination process. Future studies aimed at identifying elements or factors that interact with S9 will provide insight into the exact mechanism by which S9 affects termination.
Strains and growth conditions
Strains are listed in Table 1. To allow construction of rps9a and rps9b double knockout (strain YA178), haploid strains deleted of either 9A (Euroscarf Y02747) or 9B (Euroscarf Y16961) were mated (to allow selection for deletion of 9B, the kanMX4 in strain Y16961 was replaced by HIS3). The diploid was transformed with a plasmid expressing rps9-Δ188 from its native promoter and with its 3' UTR. This strain was then sporulated and cells that were deleted of rpS9A, 9B and expressed rps9-Δ188 were selected on the appropriate medium.
The rps9-Δ188 strain also carried a deletion of the Shm2 gene and expressed a Flag-tagged Shm2 controlled by its own promoter from a plasmid. This modification was included to allow future studies in which modification of a reporter mRNA will be necessary. Shm2 is a highly expressed, non-essential gene that is involved in the one-carbon metabolism. It does not have any known role in translation or ribosome function.
Mating, sporulation and lithium acetate transformation were all carried out using conventional methods . Yeast was grown normally to the mid-logarithmic phase either in rich media (YPD) or selective media (SD with the relevant drop outs).
The plasmids used are listed in Table 2. RPS9 constructs were made by PCR amplification of a genomic region starting from 500 bp upstream to the RPS9A gene and ending 300 bp downstream to the stop codon. This yielded a PCR product that is 1,895 bp long and includes the entire RPS9A ORF with its intron, 500 bp upstream and 300 bp downstream. This fragment was inserted into pRS415 and pRS416 vectors, and served as a template for the construction of rps9-Δ188. rps9-Δ188 was constructed by Hind3 cleavage, Klenow filling-in and religation of the plasmid. All constructs were confirmed by sequencing.
Protein extracts were prepared from trichloroacetic acid-treated yeast cells as previously described . Immunoblot procedures were essentially as described  using mouse monoclonal antibodies directed against Flag (Sigma cat. no. F3165, 1:5000 dilution).
Luciferase plasmids were kindly provided by Dr. Daniel Reines  and inserted into the relevant strains (Table 1). Cells were grown to the mid-logarithmic phase, and the expression of luciferase genes was induced by 2 hr growth in galactose. Paromomycin treatment was done by the addition of paromomycin (200 μM final concentration) 30 min after the shift to galacose, and cells were grown for an additional 90 min. Luciferase assay was performed using the Luciferase Assay System from Promega (Cat. E1500).
Yeast total RNA was extracted by the hot phenol extraction method. Extraction procedure for polysomal preparation was essentially as described  except that the lysis was performed by two rounds of 1.5 min treatment in a bead beater. RDM analysis was performed as described . Briefly, 50 ml of BY4741 cells were grown to OD600 0.6–0.9 at 30°C in YPD medium (1% yeast extract, 2% Bacto peptone, 2% glucose). Cells were harvested following cycloheximide (CHX) addition to the medium to a final concentration of 0.1 mg/ml, and lysed as previously described . Lysate was separated through a sucrose gradient and the fractions containing the majority of the tested mRNA were collected. Standard fraction volume was ~600 μl, into which DTT and Ribonuclease Inhibitor (Porcine liver, Takara) were added to a final concentration of 0.15 mM and 400 units/ml, respectively. Oligodeoxynucleotides (ODN) complementary to the regions of interest were added to a final concentration of 0.16 μM and annealing was allowed for 20 min with a gradual decrease in temperature from 37°C to room temperature. Cleavage reactions were initiated by the addition of RNase H (10 units per reaction) and 5× RNase H buffer (0.1 M Tris pH 7.4, 0.5 M KCl, 0.1 M MgCl2, 0.5 mM DTT, 2.5 mg/ml CHX). Reactions were allowed for 30 min at 37°C and terminated by the addition of 400 μl of LMD (20 mM Tris pH 7.4, 140 mM KCl, 1.5 mM MgCl2, 0.5 mM DTT, 0.1 mg/ml CHX and 1 mg/ml Heparin). Cleavage products were then separated by velocity sedimentation on a 10–50% sucrose gradient, and fractions were collected into tubes containing 1.5 fraction volumes of 8 M Guanidium HCl and 2.5 fraction volume of 100% ethanol. Pallets were resuspended and the sedimentation position of each cleavage product was determined by northern analysis . Radioactive probes were prepared by random priming incorporation of 32P-labled nucleotides into a PCR fragment homologous to the tested mRNA.
Microarrays production and analysis
DNA microarrays were produced using a contact spotter (built according to the MGuide, Stanford). All known and predicted yeast ORFs (~6,000) were PCR-amplified from start to stop codon using specific primers. Samples of each reaction were analyzed by gel electrophoresis and the remainder was precipitated and used for spotting on Amino silane coated slides (GAPS II from Corning). Following hybridization, microarrays were scanned by an Axon 4000B scanner, and spots were analyzed using the GenePix 5.1 program (Axon Instruments). Low-quality spots were filtered out if they were not distinguishable from the local background (difference less than 300), or had large variations in pixels' ratio (correlation coefficient less than 0.5). Polysomal analysis normalization was done by the addition of equal amounts of five different in vitro transcribed B. subtilis RNAs to each of the fractions immediately after collection . DNA complementary to these RNAs was spotted at least 50 times on the DNA microarray, and the signals obtained from these spots were used to normalize any global differences in the signals of different microarrays . For transcriptome analysis, normalization was done according to the median ratio of all spots on the microarray.
Peak fraction and ribosomal association assignment
Three independent experimental repeats carrying out the entire protocol from yeast growth through polysomal fractions and microarray hybridization were performed. The fraction collection parameters led to 13 or 14 fractions per collection (variation due to small differences in gradient volume). In each repeat, all of the fractions corresponding to complexes smaller than 80S (i.e., mRNA free of ribosomes) were combined and hybridized together to the DNA microarray. The rest of the fractions were not combined and were hybridized individually to DNA microarrays. Every experiment therefore includes 8 to 10 microarrays, each representing a different region of the gradient. Following normalization of the DNA microarrays signals, red-to-green ratios were obtained for each gene from the different microarrays to generate its distribution profile along the gradient. In many cases, the distribution profile included missing values, namely a signal that did not pass the quality selection in a specific microarray. For further analysis, we selected only those genes that had only a single missing value, unless it appeared in one of the extreme fractions, in which case it was not counted as missing but assumed to contain very low amounts of mRNA (Additional File 3). An exception to this single missing value rule was done in the first repeat, where fraction number 10 gave a very poor signal and the entire fraction was excluded from the analysis. This fraction included mRNAs associated with two ribosomes, and analysis of the two other repeats revealed that almost none of the mRNAs exhibited their highest association with two ribosomes.
To determine the number of ribosomes that a transcript is mostly associated with, we first defined the peak fraction of sedimentation, which is the polysomal fraction where the highest number of each gene's mRNAs sediments. In each repeat for each gene, the fraction (excluding the free fraction) with the highest ratio was defined as the peak fraction. The number of ribosomes sedimenting in the peak fraction was deduced directly from peaks in the OD254 profile of the gradient. For fractions derived from the limited resolution region of the gradient (fractions 13 and 14), the number of ribosomes was extrapolated from the high-resolution region (see Additional File 2). Next, the number of ribosomes obtained from the three repeats was averaged to yield the average number of ribosomes on a transcript (Additional File 4). Averaging was done only if two or more experimental repeats appeared to deviate by one or less fractions. Cases where the peak fraction appeared to be a single ribosome were not included in the averaging unless they were either one or two ribosomes in all three repeats. The exclusion of repeats in which the peak appeared in the single ribosome fraction was based on the assumption that mRNAs with a single ribosome represented a different translational status than mRNAs in polysomes.
We wish to thank Dr. M. Shmoish for help with the statistical analysis, Dr. C. Zilberstein for extracting the ΔG values, and Dr. D. Reines for the luciferase plasmids. This work was supported by grants from the Israeli Science Foundation and the Zimmer Cancer Research Fund.
- Dresios J, Panopoulos P, Synetos D: Eukaryotic ribosomal proteins lacking a eubacterial counterpart: important players in ribosomal function. Mol Microbiol. 2006, 59 (6): 1651-1663. 10.1111/j.1365-2958.2006.05054.x.PubMedView ArticleGoogle Scholar
- Dinos G, Wilson DN, Teraoka Y, Szaflarski W, Fucini P, Kalpaxis D, Nierhaus KH: Dissecting the ribosomal inhibition mechanisms of edeine and pactamycin: the universally conserved residues G693 and C795 regulate P-site RNA binding. Mol Cell. 2004, 13 (1): 113-124. 10.1016/S1097-2765(04)00002-4.PubMedView ArticleGoogle Scholar
- Chattoo BB, Palmer E, Ono BI, Sherman F: Patterns of Genetic and Phenotypic Suppression of lys2 Mutations in the Yeast SACCHAROMYCES CEREVISIAE. Genetics. 1979, 93 (1): 67-79.PubMed CentralPubMedGoogle Scholar
- Vincent A, Liebman SW: The yeast omnipotent suppressor SUP46 encodes a ribosomal protein which is a functional and structural homolog of the Escherichia coli S4 ram protein. Genetics. 1992, 132 (2): 375-386.PubMed CentralPubMedGoogle Scholar
- Masurekar M, Palmer E, Ono BI, Wilhelm JM, Sherman F: Misreading of the ribosomal suppressor SUP46 due to an altered 40 S subunit in yeast. J Mol Biol. 1981, 147 (3): 381-390. 10.1016/0022-2836(81)90490-3.PubMedView ArticleGoogle Scholar
- Eustice DC, Wakem LP, Wilhelm JM, Sherman F: Altered 40 S ribosomal subunits in omnipotent suppressors of yeast. J Mol Biol. 1986, 188 (2): 207-214. 10.1016/0022-2836(86)90305-0.PubMedView ArticleGoogle Scholar
- Ono BI, Stewart JW, Sherman F: Serine insertion caused by the ribosomal suppressor SUP46 in yeast. J Mol Biol. 1981, 147 (3): 373-379. 10.1016/0022-2836(81)90489-7.PubMedView ArticleGoogle Scholar
- Stansfield I, Jones KM, Herbert P, Lewendon A, Shaw WV, Tuite MF: Missense translation errors in Saccharomyces cerevisiae. J Mol Biol. 1998, 282 (1): 13-24. 10.1006/jmbi.1998.1976.PubMedView ArticleGoogle Scholar
- Ishiguro J, Ono BI, Masurekar M, McLaughlin CS, Sherman F: Altered ribosomal protein S11 from the SUP46 suppressor of yeast. J Mol Biol. 1981, 147 (3): 391-397. 10.1016/0022-2836(81)90491-5.PubMedView ArticleGoogle Scholar
- Yusupov MM, Yusupova GZ, Baucom A, Lieberman K, Earnest TN, Cate JH, Noller HF: Crystal structure of the ribosome at 5.5 A resolution. Science. 2001, 292 (5518): 883-896. 10.1126/science.1060089.PubMedView ArticleGoogle Scholar
- Yusupova GZ, Yusupov MM, Cate JH, Noller HF: The path of messenger RNA through the ribosome. Cell. 2001, 106 (2): 233-241. 10.1016/S0092-8674(01)00435-4.PubMedView ArticleGoogle Scholar
- Spahn CM, Beckmann R, Eswar N, Penczek PA, Sali A, Blobel G, Frank J: Structure of the 80S ribosome from Saccharomyces cerevisiae--tRNA-ribosome and subunit-subunit interactions. Cell. 2001, 107 (3): 373-386. 10.1016/S0092-8674(01)00539-6.PubMedView ArticleGoogle Scholar
- Takyar S, Hickerson RP, Noller HF: mRNA helicase activity of the ribosome. Cell. 2005, 120 (1): 49-58. 10.1016/j.cell.2004.11.042.PubMedView ArticleGoogle Scholar
- Mager WH, Planta RJ, Ballesta JG, Lee JC, Mizuta K, Suzuki K, Warner JR, Woolford J: A new nomenclature for the cytoplasmic ribosomal proteins of Saccharomyces cerevisiae. Nucleic Acids Res. 1997, 25 (24): 4872-4875. 10.1093/nar/25.24.4872.PubMed CentralPubMedView ArticleGoogle Scholar
- Bernstein KA, Gallagher JE, Mitchell BM, Granneman S, Baserga SJ: The small-subunit processome is a ribosome assembly intermediate. Eukaryot Cell. 2004, 3 (6): 1619-1626. 10.1128/EC.3.6.1619-1626.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Stevens SW, Ryan DE, Ge HY, Moore RE, Young MK, Lee TD, Abelson J: Composition and functional characterization of the yeast spliceosomal penta-snRNP. Mol Cell. 2002, 9 (1): 31-44. 10.1016/S1097-2765(02)00436-7.PubMedView ArticleGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431 (7004): 99-104. 10.1038/nature02800.PubMed CentralPubMedView ArticleGoogle Scholar
- Mieczkowski PA, Dominska M, Buck MJ, Gerton JL, Lieb JD, Petes TD: Global analysis of the relationship between the binding of the Bas1p transcription factor and meiosis-specific double-strand DNA breaks in Saccharomyces cerevisiae. Mol Cell Biol. 2006, 26 (3): 1014-1027. 10.1128/MCB.26.3.1014-1027.2006.PubMed CentralPubMedView ArticleGoogle Scholar
- Arava Y, Wang Y, Storey JD, Liu CL, Brown PO, Herschlag D: Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae. Proc Natl Acad Sci U S A. 2003, 100 (7): 3889-3894. 10.1073/pnas.0635171100.PubMed CentralPubMedView ArticleGoogle Scholar
- van de Peppel J, Kemmeren P, van Bakel H, Radonjic M, van Leenen D, Holstege FC: Monitoring global messenger RNA changes in externally controlled microarray experiments. EMBO Rep. 2003, 4 (4): 387-393. 10.1038/sj.embor.embor798.PubMed CentralPubMedView ArticleGoogle Scholar
- Sharp PM, Li WH: The codon Adaptation Index--a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 1987, 15 (3): 1281-1295. 10.1093/nar/15.3.1281.PubMed CentralPubMedView ArticleGoogle Scholar
- Bennetzen JL, Hall BD: Codon selection in yeast. J Biol Chem. 1982, 257 (6): 3026-3031.PubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.PubMed CentralPubMedView ArticleGoogle Scholar
- Hurowitz EH, Brown PO: Genome-wide analysis of mRNA lengths in Saccharomyces cerevisiae. Genome Biol. 2003, 5 (1): R2-10.1186/gb-2003-5-1-r2.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang Y, Liu CL, Storey JD, Tibshirani RJ, Herschlag D, Brown PO: Precision and functional specificity in mRNA decay. Proc Natl Acad Sci U S A. 2002, 99 (9): 5860-5865. 10.1073/pnas.092538799.PubMed CentralPubMedView ArticleGoogle Scholar
- Miura F, Kawaguchi N, Sese J, Toyoda A, Hattori M, Morishita S, Ito T: A large-scale full-length cDNA analysis to explore the budding yeast transcriptome. Proc Natl Acad Sci U S A. 2006, 103 (47): 17846-17851. 10.1073/pnas.0605645103.PubMed CentralPubMedView ArticleGoogle Scholar
- Preiss T, Baron-Benhamou J, Ansorge W, Hentze MW: Homodirectional changes in transcriptome composition and mRNA translation induced by rapamycin and heat shock. Nat Struct Biol. 2003, 10 (12): 1039-1047. 10.1038/nsb1015.PubMedView ArticleGoogle Scholar
- Lotan R, Bar-On VG, Harel-Sharvit L, Duek L, Melamed D, Choder M: The RNA polymerase II subunit Rpb4p mediates decay of a specific class of mRNAs. Genes Dev. 2005, 19 (24): 3004-3016. 10.1101/gad.353205.PubMed CentralPubMedView ArticleGoogle Scholar
- MacKay VL, Li X, Flory MR, Turcott E, Law GL, Serikawa KA, Xu XL, Lee H, Goodlett DR, Aebersold R, Zhao LP, Morris DR: Gene expression analyzed by high-resolution state array analysis and quantitative proteomics: response of yeast to mating pheromone. Mol Cell Proteomics. 2004, 3 (5): 478-489. 10.1074/mcp.M300129-MCP200.PubMedView ArticleGoogle Scholar
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMedGoogle Scholar
- Peske F, Rodnina MV, Wintermeyer W: Sequence of steps in ribosome recycling as defined by kinetic analysis. Mol Cell. 2005, 18 (4): 403-412. 10.1016/j.molcel.2005.04.009.PubMedView ArticleGoogle Scholar
- Zavialov AV, Hauryliuk VV, Ehrenberg M: Splitting of the posttermination ribosome into subunits by the concerted action of RRF and EF-G. Mol Cell. 2005, 18 (6): 675-686. 10.1016/j.molcel.2005.05.016.PubMedView ArticleGoogle Scholar
- Hirokawa G, Nijman RM, Raj VS, Kaji H, Igarashi K, Kaji A: The role of ribosome recycling factor in dissociation of 70S ribosomes into subunits. Rna. 2005, 11 (8): 1317-1328. 10.1261/rna.2520405.PubMed CentralPubMedView ArticleGoogle Scholar
- Treger M, Westhof E: Statistical analysis of atomic contacts at RNA-protein interfaces. J Mol Recognit. 2001, 14 (4): 199-214. 10.1002/jmr.534.PubMedView ArticleGoogle Scholar
- Lejeune D, Delsaux N, Charloteaux B, Thomas A, Brasseur R: Protein-nucleic acid recognition: statistical analysis of atomic interactions and influence of DNA structure. Proteins. 2005, 61 (2): 258-271. 10.1002/prot.20607.PubMedView ArticleGoogle Scholar
- Guthrie C, Fink G: Guide to yeast genetics and molecular and cell biology. 2002, Academic pressGoogle Scholar
- Foiani M, Liberi G, Lucchini G, Plevani P: Cell cycle-dependent phosphorylation and dephosphorylation of the yeast DNA polymerase alpha-primase B subunit. Mol Cell Biol. 1995, 15 (2): 883-891.PubMed CentralPubMedView ArticleGoogle Scholar
- Ausubel FM, Brent R, Kingston RE, Moore DD, Smith JA, Seidman JG, Struhl K: Current protocols in molecular biology.Google Scholar
- Arava Y: Isolation of polysomal RNA for microarray analysis. Methods Mol Biol. 2003, 224: 79-87.PubMedGoogle Scholar
- Arava Y, Boas FE, Brown PO, Herschlag D: Dissecting eukaryotic translation and its control by ribosome density mapping. Nucleic Acids Res. 2005, 33 (8): 2421-2432. 10.1093/nar/gki331.PubMed CentralPubMedView ArticleGoogle Scholar
- Kuhn KM, DeRisi JL, Brown PO, Sarnow P: Global and Specific Translational Regulation in the Genomic Response of Saccharomyces cerevisiae to a Rapid Transfer from a Fermentable to a Nonfermentable Carbon Source. Mol Cell Biol. 2001, 21 (3): 916-927. 10.1128/MCB.21.3.916-927.2001.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.