Translational control of one-carbon metabolism underpins ribosomal protein phenotypes in cell division and longevity
Figures

Doubling time and normal Rpl levels are not associated with the replicative lifespan of single rpl deletion mutants.
(A) Scatterplot between the generation time (x-axis; from Tables S2 in Steffen et al., 2012) and replicative lifespan (y-axis; as percentage of the wild type lifespan, from Table 3 in Steffen et al., 2012 and Table S2 in McCormick et al., 2015). (B) Scatterplot between the abundance of each deleted Rpl typically found in wild type cells (x-axis; the median number of molecules per cell, from Table S4 from Ho et al., 2018) and replicative lifespan (y-axis; as percentage of the wild type lifespan, from Table 3 in Steffen et al., 2012 and Table S2 in McCormick et al., 2015). The Spearman correlation coefficients (ρ) shown in each case were calculated with the rcorr function of the Hmisc R language package, comparing the pair shown in each panel. All the values used as input for this figure and analyses are in Figure 1—source data 1.
-
Figure 1—source data 1
Lifespan, doubling time, and RP levels in ribosomal protein mutants.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig1-data1-v2.xlsx

Loss of Rpl22 does not affect the relative abundance of other ribosomal proteins in ribosomes.
(A) Schematic overview of the approach to query ribosomal protein abundances in rpl22 deletion mutants (see Materials and methods). All strains were in the haploid BY4742 background. The doubling time values for each strain are shown on top (Td; in min, from Table S2 in Steffen et al., 2012). (B) Bar plot displaying on the x-axis the levels of Rpl22Ap and Rpl22Bp in wild type and rpl22 deletion mutants, as indicated on the y-axis. The abundances shown correspond to the peak area intensities obtained from SWATH-MS (see Materials and methods and Figure 2—source data 1). (C) The Log2-transformed ratios of the intensities corresponding to the ribosomal proteins detected are indicated for each pairwise comparison among the 4 strains we analyzed. Changes higher than 2-fold are outside the dashed red lines in each panel. Note that not all the changes were statistically significant, as indicated. The red arrows indicate the only cases where the differences were significant (Log2FC ≥ 1 and p-value<0.05; bootstrap-based ANOVA; see Materials and methods) The Spearman correlation coefficients (ρ) shown in each case were calculated with the rcorr function of the Hmisc R language package, comparing the pair shown in each panel.
-
Figure 2—source data 1
SWATH-mass spectrometry measurements of ribosomal protein abundances.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig2-data1-v2.xlsx

Reduced protein synthesis in cells lacking Rpl22Ap.
(A) Schematic of the approach to measure incorporation of a methionine analog (HPG) in newly synthesized proteins (see Materials and methods). (B) The strains shown on the x-axis (all in the BY4742 background, the same as in Figure 1) were incubated with HPG for 30 m, processed for CLICK reaction with the Alexa dye, and then the fluorescence per cell was quantified with flow cytometry. At least 8,000–10,000 cells were quantified in each sample, from four independent experiments in each case. On the y-axis is the mean fluorescence per cell for each sample divided by the average of the four WT samples. Significant differences in pair-wise comparisons were indicated by the non-parametric Kruskal-Wallis test (p=0.004848), while the p-values shown were based on the posthoc Nemenyi test (performed with the PMCMR R language package). (C) Representative images of HPG-Alexa-labeled cells, which were also stained with DAPI to visualize the nuclei.

Querying synchronous, dividing cells lacking ribosomal protein paralogs.
(A) Schematic of our experimental approach to identify gene expression changes in rpl22 and rpl34 paralog mutants at the transcriptional or translational level, from cells at several different stages of the cell cycle. (B) Cell size (in fL, x-axis) and the percentage of budded cells (%Budded, y-axis) from each pool of the indicated strains. The values shown are the weighted averages, from the different elutriated samples in each pool. (C) On the y-axis, are the Log2-transformed TPM values from representative transcripts known to increase in abundance in S (histone H4; HHF1), G2 phase (cyclin CLB2), or are constitutively expressed (ACT1). Cell size is on the x-axis.
-
Figure 3—source data 1
Raw read sequencing data.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig3-data1-v2.xlsx
-
Figure 3—source data 2
TPM-normalized sequencing data.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig3-data2-v2.xlsx

Transcripts with altered relative abundance in rpl22 and rpl34 paralog deletion mutants.
(A) Decision diagram for identifying transcripts that were differentially expressed in paralog deletion mutants. (B) The number of transcripts with significantly different levels (adjusted p-value/FDR < 0.05, Log2FC ≥ 2, identified both from the babel and DESeq2 R language packages) at any one cell size pool, between any pairwise comparison between WT, rpl22, and rpl34 mutants is shown. The data for wild type diploid cells were from Blank et al., 2017. The Gene Ontology terms (The Gene Ontology Consortium, 2019) that were most enriched in each case are shown, based on the PANTHER (Mi et al., 2019) platform classification, incorporating the Holm-Bonferroni correction. All the gene names of loci with significant changes in transcript levels are in Figure 4—source data 1 (sheet ‘cumulative’). The ‘UP’ and ‘DOWN’ groupings correspond to each pairwise comparison shown, with the ‘UP’, or ‘DOWN’, group higher, or lower, for the strain in the numerator of the ratio, respectively. (C) Reduced levels of transcripts encoding glycolytic/gluconeogenic enzymes in rpl22aΔ compared to rpl22bΔ cells. The Log2-transformed ratio of the corresponding TPM values is shown in each case, across the different cell size pools. The data were hierarchically clustered without supervision and displayed with the pheatmap R package. Each row corresponds to a separate 5 fL cell size interval, from 40 to 75 fL, from top to bottom. The arrows indicate transcripts for key enzymes of glycolysis and gluconeogenesis. mRNAs with missing values across the cell sizes we analyzed were not included in the heatmap.
-
Figure 4—source data 1
mRNAs with significantly different abundances between any two strains.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig4-data1-v2.xlsx

Reduced translational efficiency of transcripts encoding enzymes of the methionine and 1C metabolic pathways in rpl22aΔ cells.
(A) Heatmap of the transcripts with significantly different translational efficiencies (see Figure 3) between the two Rpl22 paralog deletion strains. The ratio of the TPM values for the ribosome footprints ((ribo) against the corresponding values for the mRNA reads (rna) define the translational efficiency of each locus in each paralog mutant. The data were hierarchically clustered and displayed with the pheatmap R package. Each row corresponds to a separate cell size interval, from 40 to 75 fL, from top to bottom. The Gene Ontology terms highlighted were enriched, based on the PANTHER platform classification, incorporating the Holm-Bonferroni correction. (B) Diagram of the metabolic pathways involving enzymes whose mRNAs have altered translational efficiency in rpl22aΔ vs rpl22bΔ cells. Proteins whose mRNAs have lower translational efficiency in rpl22aΔ cells are shown in blue. (C) Immunoblots from WT (RPL22+), rpl22aΔ, or rpl22bΔ cells, carrying a MET3-TAP allele expressed from its endogenous chromosomal location. All the strains were otherwise isogenic. The signal corresponding to Met3p-TAP was detected with the PAP reagent, while the Pgk1p signal represents loading. The band intensities were quantified using the ImageJ software package, and the relative abundance of Met3p-TAP in rpl22aΔ: rpl22bΔ cells is shown at the bottom, from three independent such experiments.
-
Figure 5—source data 1
mRNAs with significantly different translational efficiency between any two strains.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig5-data1-v2.xlsx

Transcripts with altered translational efficiency (TE) in rpl22 and rpl34 paralog deletion mutants.
(A) Decision diagram for identifying transcripts that had altered translational efficiency in paralog deletion mutants. (B) The number of transcripts with significantly different levels (adjusted p-value/FDR < 0.05, Log2FC ≥ 1, identified from all three statistical packages) at any one cell size pool, between any pairwise comparison between WT, rpl22, and rpl34 mutants is shown. The data for wild type diploid cells were from Blank et al., 2017. The Gene Ontology terms that were most enriched in each case are shown, based on the PANTHER platform classification, incorporating the Holm-Bonferroni correction. All the gene names of loci with significant changes in translational efficiency are also in Figure 5—source data 1 (sheet ‘cumulative’). The ‘UP’ and ‘DOWN’ groupings correspond to each pairwise comparison shown, with the ‘UP’, or ‘DOWN’, group higher, or lower, for the strain in the numerator of the ratio, respectively.

The translational efficiency of GCN4 is de-repressed in cells lacking Rpl22Ap.
The Log2-tranformed ratio of the TPM values for the ribosome footprints ((ribo) against the corresponding values for the mRNA reads (rna) define the translational efficiency of GCN4 (y-axis), for each of the cell sizes we queried (x-axis), in the strains shown. Asterisks indicate that the translational efficiency was significantly reduced (FDR < 0.05) compared to other mRNAs in that strain (based on the babel ‘combined’ output).

Dysregulation of the translational efficiency of transcripts encoding gene products involved in cytoplasmic translation in cells lacking Rpl22 or Rpl34 paralogs.
Heatmaps of the transcripts with significantly different translational efficiencies (see Figure 2) between wild type diploid cells (from Blank et al., 2017) and rpl22aΔ/rpl22aΔ (A), rpl22bΔ/rpl22bΔ (B), rpl34aΔ/rpl34aΔ (C), rpl34bΔ/rpl34bΔ (D) cells. The ratio of the TPM values for the ribosome footprints ((ribo) against the corresponding values for the mRNA reads (rna) define the translational efficiency of each locus in each strain. The data were hierarchically clustered and displayed with the pheatmap R package. Each column corresponds to a separate cell size interval, from 40 to 75 fL, from left to right. The Gene Ontology terms highlighted were enriched, based on the PANTHER platform classification, incorporating the Holm-Bonferroni correction. mRNAs with missing values across the cell sizes we analyzed were not included in the heatmaps. All the gene names of loci with significant changes in translational efficiency are also in Figure 5—source data 1 (sheet ‘cumulative’).

Features of mRNAs with altered translational efficiency in rpl mutants.
Top, All the mRNAs shown in Figure 5—figure supplement 3 that had altered translational efficiency between wild type and any of the rpl deletion strains were placed in the following groups: UP, if their median of their log2[(ribo/rna)WT/(ribo/rna) rplΔ] values from all cell cycle points was positive; DOWN, if their median of their log2[(ribo/rna)WT/(rio/rna)rplΔ] values from all cell cycle points was negative; OTHER if they were not identified as being under translational control between wild type and any of the rpl deletion strains. The values for ORF, 5’-leader (5’-UTR), and 3’-UTR length for each of the mRNAs were from Lin and Li, 2012. The protein abundance values were the median number of molecules per cell, from the studies analyzed in Ho et al., 2018. The violin plots were generated with the lattice R language package. Significant differences in pair-wise comparisons were indicated by the non-parametric Kruskal-Wallis test, while the p-values shown were for all significant differences based on the posthoc Nemenyi test (performed with the PMCMR R language package). Bottom, The same analysis as in the top panels for the mRNAs with different translational efficiency between the rpl22aΔ and rpl22bΔ mutants, based on their log2[(ribo/rna)rpl22aΔ/(ribo/rna)rpl22bΔ] values.

Metabolic profiling indicates reduced flux through central metabolic pathways and the folate cycle in rpl22aΔ cells.
(A) Eleven metabolites shown at the bottom had significantly reduced levels in rpl22aΔ cells (Log2FC ≥ 1, p<0.05; based on bootstrapped ANOVA; see Materials and methods) and they were significantly enriched for the metabolic pathways shown to the right (FDR < 0.05). Pathway enrichment analysis was done with the MetaboAnalyst R language package. The metabolites were identified with untargeted, MS-based profiling of primary metabolites and biogenic amines, and targeted amino acid analysis. Metabolites indicated with gray in the Table are part of the pathways shown to the right. (B) The Log2-transformed peak intensities from the MS-based profiling of the metabolites shown in A (except Glycine) are on the y-axis. The strains used in the analysis are on the x-axis. (C) The Log2-transformed levels (in nmoles) of amino acids, after PTH-derivatization, Edman degradation and HPLC detection, are shown on the y-axis. The red arrows indicate the only amino acids (Gly and Trp) whose levels were significantly lower in rpl22aΔ cells (Log2FC ≥ 1, p<0.05; based on bootstrapped ANOVA; see Materials and Methods). The strains used in the analysis are on the x-axis, and they were in the BY4742 background.
-
Figure 6—source data 1
Metabolite and amino acid abundances.
- https://cdn.elifesciences.org/articles/53127/elife-53127-fig6-data1-v2.xlsx

Exogenous addition of metabolites from pathways affected in rpl22aΔ cells does not suppress the slower growth of these cells.
Cultures from the indicated strains shown to the left (in the diploid BY4743 background in A and B; in the haploid BY4742 background in C), were spotted onto solid plates with standard undefined media (YPD; see Materials and Methods. In each case, aliquots from saturated cultures were spotted in 10-fold serial dilutions. The plates contained the compounds shown on top in each case (GSH is reduced glutathione and GSSG is oxidized glutathione), at the indicated concentrations. The plates were photographed after they were incubated for 3 days at 30°C.

Nobiletin is accumulated in rpl22aΔ cells.
(A) Structure of nobiletin. (B) Peak intensities assigned to nobiletin from mass spectrometry are on the y-axis. The strains used in the analysis are on the x-axis, and they were in the BY4742 background.

Deletion of enzymes of one-carbon metabolic pathways extends replicative lifespan in yeast.
Survival curves for MATα (BY4742) cells (shown in black), compared to experiment-matched cells (shown in red) lacking SHM1 (A), SHM2 (B), ADE17 (C), ADE2 (D), or ADE3 (E). Mean lifespans are shown in parentheses, along with the number of cells assayed in each case. In the case of shm1Δ, shm2Δ, ade17Δ, and ade3Δ cells, the lifespan extension was significant (p<0.0001; based on the log-rank test).

Deletion of enzymes of one-carbon metabolic enzymes does not further extend the replicative longevity of rpl22aΔ cells.
Survival curves for wild type MATα (BW885) cells (shown in black), compared to experiment-matched cells of the indicated genotype. Mean lifespans are shown in parentheses. The lifespans were determined from 60 cells in each case.

Loss of Shm2p and Ade3p impinges on multiple cell cycle phases but in distinct ways.
(A) Double shm2Δ,ade3Δ deletion mutants are slow-growing. Two representative tetrad dissections from shm2Δ x ade3Δ crosses are shown. The yellow diamond indicates the shm2Δ,ade3Δ segregants. All the strains used in B-E were segregants from the same shm2Δ x ade3Δ crosses and, except as indicated, isogenic otherwise. (B) Violin plots of the mean and birth size of the indicated strains, calculated from ≥12 asynchronous cultures in each case. The plots were generated with the lattice R language package. Significant differences in pair-wise comparisons were indicated by the non-parametric Kruskal-Wallis test, while the p-values shown were for all significant differences based on the posthoc Nemenyi test (performed with the PMCMR R language package). (C) Representative DNA content histograms from the indicated strains, from at least 10,000 cells and ≥5 independent asynchronous cultures in each case. On the x-axis is fluorescence per cell, while the cell number is on the y-axis. The average and sd of the percentage of cells with unreplicated DNA (%G1) is shown. There were no statistically significant differences among the strains, based on the non-parametric Kruskal-Wallis test. (D) The rate of cell size increase (k, in h−1, shown on the y-axis) was calculated as described in Soma et al., 2014, for the strains shown on the x-axis. Each data point is the average of two technical replicates, from synchronous elutriated cultures. (E) The critical size (y-axis) for the strains shown on the x-axis was calculated as described in Soma et al., 2014, from the same cultures as in D.

The translational efficiency of SHM2 is cell cycle-regulated.
Cell size (in fL, x-axis) and the translational efficiency (TE) of SHM2, shown as Log2(expressed ratios) across the cell size series (y-axis), in the indicated strains. The values shown are the average from the three independent pools, for each 5 fL interval across.
Tables
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Strain, strain background (S. cerevisiae) | BY4743 | Giaever et al., 2002 | RRID:SCR_003093 | MATa/α his3Δ1/his3Δ1 leu2Δ0/leu2Δ0 LYS2/lys2Δ0 met15Δ0/MET15 ura3Δ0/ura3Δ0 |
Strain, strain background (S. cerevisiae) | 32672 | Giaever et al., 2002 | RRID:SCR_003093 | rpl22aΔ::KanMX/rpl22aΔ::KanMX, BY4743 otherwise |
Strain, strain background (S. cerevisiae) | 35844 | Giaever et al., 2002 | RRID:SCR_003093 | rpl22bΔ::KanMX/rpl22bΔ::KanMX, BY4743 otherwise |
Strain, strain background (S. cerevisiae) | 30192 | Giaever et al., 2002 | RRID:SCR_003093 | rpl34aΔ::KanMX/rpl34aΔ::KanMX, BY4743 otherwise |
Strain, strain background (S. cerevisiae) | 31445 | Giaever et al., 2002 | RRID:SCR_003093 | rpl34bΔ::KanMX/rpl34bΔ::KanMX, BY4743 otherwise |
Strain, strain background (S. cerevisiae) | BY4742 | Steffen et al., 2012 | MATα his3Δ1 leu2Δ0 lys2Δ0 ura3Δ0 | |
Strain, strain background (S. cerevisiae) | KS976 | Steffen et al., 2012 | rpl22aΔ::URA3, BY4742 otherwise | |
Strain, strain background (S. cerevisiae) | KS979 | Steffen et al., 2012 | rpl22bΔ::URA3, BY4742 otherwise | |
Strain, strain background (S. cerevisiae) | KS999 | Steffen et al., 2012 | rpl22a,bΔ::URA3, BY4742 otherwise | |
Strain, strain background (S. cerevisiae) | BY4741 | Giaever et al., 2002 | RRID:SCR_003093 | MATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0 |
Strain, strain background (S. cerevisiae) | MET3-TAP | Dharmacon | YSC1178-202231887 | MET3-TAP::HIS3M × 6, BY4741 otherwise |
Strain, strain background (S. cerevisiae) | HB147 | This study | rpl22aΔ::URA3, MET3-TAP::HIS3M × 6, BY4741 otherwise | |
Strain, strain background (S. cerevisiae) | HB171 | This study | rpl22bΔ::URA3, MET3-TAP::HIS3M × 6, BY4741 otherwise | |
Strain, strain background (S. cerevisiae) | 13403 | Giaever et al., 2002 | RRID:SCR_003093 | shm1Δ::KanMX, BY4742 otherwise |
Strain, strain background (S. cerevisiae) | 12669 | Giaever et al., 2002 | RRID:SCR_003093 | shm2Δ::KanMX, BY4742 otherwise |
Strain, strain background (S. cerevisiae) | 16561 | Giaever et al., 2002 | RRID:SCR_003093 | ade17Δ::KanMX, BY4742 otherwise |
Strain, strain background (S. cerevisiae) | 12384 | Giaever et al., 2002 | RRID:SCR_003093 | ade2Δ::KanMX, BY4742 otherwise |
Strain, strain background (S. cerevisiae) | 6591 | Giaever et al., 2002 | RRID:SCR_003093 | ade3Δ::KanMX, BY4741 otherwise |
Strain, strain background (S. cerevisiae) | NM64 | This study | MATα, shm2Δ::KanMX, ade3Δ::KanMX, met -, lys- | |
Strain, strain background (S. cerevisiae) | NM65 | This study | MATa, shm2Δ::KanMX, ade3Δ::KanMX, met -, lys- | |
Strain, strain background (S. cerevisiae) | NM66 | This study | MATα, ade17Δ::KanMX, rpl22aΔ::URA, his-, lys-, leu- | |
Strain, strain background (S. cerevisiae) | NM67 | This study | MATα, shm2Δ::KanMX, rpl22aΔ::URA, his-, leu-, met- | |
Strain, strain background (S. cerevisiae) | NM68 | This study | MATα, shm1Δ::KanMX, rpl22aΔ::URA, his-, leu-, met- | |
Strain, strain background (S. cerevisiae) | BW885 | This study | MATα his3Δ1 leu2Δ0 lys2Δ0 ura3Δ0 | |
Other | Yeast extract | Sigma-Aldrich | Y1625 | |
Other | Peptone | Sigma-Aldrich | P5905 | |
Chemical compound, drug | Dextrose | Sigma-Aldrich | D9434 | |
Chemical compound, drug | Cycloheximide | Calbiochem | 239763 M | |
Chemical compound, drug | Sodium azide | Sigma-Aldrich | S2002 | |
Chemical compound, drug | Tris(hydroxymethyl)aminomethane | Sigma-Aldrich | 252859 | |
Chemical compound, drug | Tris base | Roche | TRIS-RO | |
Chemical compound, drug | Sodium chloride | Sigma-Aldrich | S7653 | |
Chemical compound, drug | Magnesium chloride hexahydrate | USP | 1374248 | |
Chemical compound, drug | DTT | Sigma-Aldrich | D0632 | |
Chemical compound, drug | Triton X-100 | Sigma-Aldrich | T8787 | |
Peptide, recombinant protein | Turbo DNase I | ThermoFisher | AM2238 | |
Other | Glass beads | Scientific Industries | SI-BG05 | |
Other | 13 × 51 mm polycarbonate centrifuge tubes | Beckman Coulter | 349622 | |
Chemical compound, drug | Sucrose | Sigma-Aldrich | S0389 | |
Chemical compound, drug | Phosphate buffered saline (PBS) | Sigma-Aldrich | P4417 | |
Commercial assay or kit | Click-iT HPG Alexa Fluor 488 Protein Synthesis Assay Kit | ThermoFisher | C10428 | |
Chemical compound, drug | DAPI (4',6-Diamidino-2-Phenylindole, Dihydrochloride) | ThermoFisher | D1306 | |
Commercial assay or kit | Ribo-Zero Magnetic Gold Kit (Yeast) | Epicentre | MRZY1324 | |
Commercial assay or kit | SciptSeq v2 RNA-Seq Library Preparation Kit | Epicentre | SSV21124 | |
Antibody | Peroxidase Anti-Peroxidase (PAP) Soluble Complex | Sigma-Aldrich | P1291 | (1:1000) |
Antibody | Anti-Pgk1p antibody, rabbit polyclonal | abcam | ab38007 | (1:1000) |
Other | NovexWedgeWell4–12% Tris-Glycine gels | ThermoFisher | XP04125 | |
Software, algorithm | MetaboAnalyst | https://www.metaboanalyst.ca/ | RRID:SCR_015539 | Web server for statistical, functional and integrative analysis of metabolomics data |
Software, algorithm | AccuComp Z2 | Beckman Coulter | 383550 | Software to monitor number and size of cells with Z2 cell counter |
Software, algorithm | NIS-Elements | https://www.nikoninstruments.com/Products/Software | RRID:SCR_014329 | Microscope imaging software suite used with Nikon products |
Software, algorithm | ImageJ | https://imagej.net/ | RRID:SCR_003070 | Image processing software |
Software, algorithm | Adobe Photoshop | https://www.adobe.com/products/photoshop.html | RRID:SCR_014199 | Image processing software |
Software, algorithm | RStudio | http://www.rstudio.com/ | RRID:SCR_000432 | Software for the R statistical computing environment |
Software, algorithm | SGD | http://www.yeastgenome.org/ | RRID:SCR_004694 | Saccharomyces Genome Database |
Software, algorithm | R | https://www.r-project.org | v3.5.2 RRID:SCR_001905 | Statistical Computing Environment |
Software, algorithm | PANTHER | http://www.geneontology.org/ | RRID:SCR_002811 | Gene ontology enrichment analysis |