BACKGROUND
methionine is one of two sulfur-containing amino acids that are incorporated into proteins during translation.
among twenty amino acids, methionine plays a special role in the biosynthesis of proteins because its codon aug is also the most common translation initiation codon.
in eukaryotes, the binding of the anticodon of the initiator met-trna to the initiation codon aug is required for initiating translation  <cit> .
this interaction is highly conserved across species.
met-trna is required for the assembly of 40s ribosome and thus may regulate the mechanism of ribosome scanning and entry, potentially serving as an important control point for translation .
since translational regulation is a key step in gene regulation, sensing intracellular methionine level and adjusting the global gene expression program through translational control may be an important strategy to coordinate cell’s metabolic state with its growth.
methionine has also been known to play important roles in a wild range of biological phenomena including growth, development, fertility, cancer and aging .
it has been widely reported that methionine intervention can effectively regulate the lifespan of numerous model organisms.
in particular, methionine restriction  has been shown to extend the lifespan of a range of species, including yeast, worm, fly and mouse .
it has also been suggested that the lifespan extension by caloric restriction, defined as reduced caloric intake without malnutrition, can be attributed to methionine restriction  <cit> .
in addition to the effect on lifespan, methionine restriction also slows or reduces many characteristics associated with senescence, such as immune and lens aging, increased igf-i and insulin levels, and cumulated oxidative damages  <cit> .
methionine restriction has also been studied extensively in anticancer therapies, either alone or in association with the other treatments, and is considered as a useful therapeutic strategy for treating various cancers  <cit> .
thus, characterizing the global gene expression program induced by metr and understanding the mechanisms by which metr regulates gene expression are important not only for understanding the basic principles of gene regulation but also for promoting human health.
translational regulation by general amino acid starvation has been extensively studied and the pathway involved has been elucidated before  <cit> .
in the canonical model, amino acid starvation leads to the accumulation of uncharged trna, activating the gcn <dig> kinase, which phosphorylates eif <dig> , ultimately affecting the translation  <cit> .
as a general strategy for sensing amino acid depletion, this may also be the mechanism to sense and respond to metr.
methionine may also work through other mechanisms to affect translation.
it has been reported that intracellular methionine availability can regulate cellular translational capacity and metabolic homeostasis by controlling the thiolation status of the wobble-uridine  nucleotides on lysine, glutamine, or glutamate trnas  <cit> .
methionine may also affect gene expression by converting to s-adenosyl methionine  <cit> , which serves as the predominant methyl donor for rrna-trna modifications and dna/protein methylations.
although there has been significant progress in understanding the roles methionine may play in gene regulation, a systematic study on the global gene expression program controlled by methionine, especially at the translational level, is still lacking.
in this work, we use ribosomal profiling and rna-seq to compare the translational and transcriptional profiles of cells growing in the normal and methionine restricted media.
we systematically characterize the translational and transcriptional programs induced by methionine restriction and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast s. cerevisiae as the model organism.
methods
yeast strains and media
yeast strains used for the ribosomal profiling and rna-seq experiments were by <dig> .
synthetic dextrose  medium and methionine restriction  medium was used in the ribosome profiling/rna-seq experiments, the sd medium contained 2%  glucose,  <dig>  g/l yeast nitrogen base  without amino acid, 20 mg/l adenine, 20 mg/l l-arginine hcl, 100 mg/l l-aspartic acid, 20 mg/l l-histidine hcl, 100 mg/l l-leucine, 30 mg/l l-isoleucine, 30 mg/l l-lysine hcl, 20 mg/l l-methionine, 50 mg/l l-phenylalanine, 200 mg/l l-threonine, 20 mg/l l-tryptophan, 30 mg/l l-tyrosine, 20 mg/l uracil, 150 mg/l l-valine, 100 mg/l glutamic acid and 4 g/l serine.
the metr media has the similar ingredients as the sd medium except for 4 mg/l l-methionine concentration.
media were freshly made before the experiments.
all nutrients were purchased from sigma-aldrich corporation.
ribosome profiling and rna-seq of cells growing in sd vs. metr media
the initial cell culture was incubated in 300 ml sd medium overnight to an od <dig>  <dig>  ~  <dig> , then diluted by five fold using fresh sd media and incubated for another 4 h under  <dig> °c to an od <dig>  <dig>  ~  <dig> .
the sample was then divided equally into two aliquots.
cells were separated from the media by spin-down at 3000 g for 5 min and re-suspended in sd and metr media respectively.
all samples were incubated for another hour before harvesting.
all the steps were carried out at 30 °c.
ribosomal profiling and rna-seq experiments were carried out using the protocol developed by ingolia et al.
<cit>  raw sequences were obtained from illumina hiseq  <dig>
sequence analysis and quantification of differential gene expression
sequence reads were aligned to the most recent s. cerevisiae genome using soapaligner/soap <dig>  with default setting  <cit> .
after trimming off the adapters, reads aligned to rrna and trna sequences were filtered out.
the rest of the reads were then aligned to the genome sequence.
finally, reads that did not align to the genome were aligned to all the cdss to retain those covering the splicing junctions.
after the alignment, we counted the number of reads starting at each position across the whole genome; for ribosomal footprinting data, the starting position of each read was shifted by 15 bps towards the 3′ end, to adjust for the offset due to ribosome protection.
to get the abundance of reads covering each gene, we sum all reads with starting position from the start to the stop codon, excluding the first 50 bps from the transcription start site , to alleviate the effect of the biased distribution of reads around the tss  <cit> .
for a reliable gene expression comparison between two conditions, we excluded genes with less than  <dig> total raw reads   <cit> .
we computed the fold change of mrna or footprint as the ratio of the corresponding reads from the two conditions, with total reads normalized to adjust the median fold change to  <dig>  to estimate the statistical significance of the fold change, we did the following analysis: for each gene g
i, we collected the other  <dig> genes with the most similar number of reads and compute the standard deviation σ
i from the log  of these  <dig> genes.
we then computed the zi=logfiσi, where f
i is the fold change of g
i comparing metr and sd.
then a p-value was calculated from the z
i to measure the significance of the fold change.
quantification of translational efficiency changes
the translational efficiency changes were calculated as the ratio of ribosomal footprints fold change to mrna fold changes for each gene.
translational efficiency change z-score is calculated by normalized the efficiency change with the standard deviation.
calculation of the tf module z-scores and kegg pathway z-scores
we used the transcription factor targets from the analysis by mcissac et al.
based on the systematic chip-chip data  <cit> , using  <dig>  as the p-value cutoff and the strongest conservation between species  <cit> .
we computed the rank sum test z-scores comparing the fold changes of the target genes vs. none-target genes for each transcription factor.
the sign of the z-score reflects the overall direction of the gene expression change in the modules; positive z-score indicates overall induction and negative z-score indicates overall repression.
we used tf modules with at least  <dig> targets for this analysis.
we used a similar method to compute the kegg  <cit>  pathway z-scores by grouping the genes from the same kegg pathway.
flow cytometer measurement of the protein abundance changes upon methionine restriction
yeast gfp-tag strains were selected from the yeast gfp library  <cit> .
for each gfp strain, we picked three single clones from the plate and cultured them overnight to saturation.
we then diluted each cell culture on the second day and grew them to od <dig>  <dig> – <dig>  in a  <dig> well plate.
cells were collected  and re-suspended in metr or sd media.
we then used flow cytometer to measure the gfp signals  in each sample after 4 h .
the cellular gfp concentration was computed by normalizing the gfp signal with the cell size using forward scattering signal  for each individual cell.
the difference of gfp concentration between metr and sd was computed as the ratio of the medians of the normalized gfp under the two conditions.
then the mean gfp fold change was calculated from the three biological replicates.
RESULTS
ribosomal profiling and rna-seq revealed a broad spectrum of transcriptional and translational changes induced by metr
we performed global gene expression profiling by rna-seq and ribosomal profiling  <cit> , comparing cells growing in normal vs. metr  conditions.
the sequencing result is of high quality with at least 50× coverage per sample .
ribosome profiling quantifies ribosome protected rna  for all the genes in the genome and thus can measure the translational efficiency when combined with the total amount rna from the rna-seq data.
we observed a broad spectrum of gene expression changes in response to metr, both at the transcriptional and the translational levels .
for the majority of the genes that changed expression, the regulation is at the transcriptional level, as the fold change of the footprint is proportional to the fold change of the transcript level .
there is a subset of genes whose translational efficiency  are increased or decreased compared to all the genes , indicating that they are under translational regulation.
we observed that a subset of genes with decreased transcriptional level tend to have a decreased translational efficiency, suggesting that they are under both transcriptional and translational control .
overall there are 110/ <dig> genes whose footprints went up/down by more than four-fold under metr, and 149/ <dig> genes increased/decreased their translational efficiency by more than two-fold.fig.
<dig> rna-seq and ribosomal profiling revealed global transcriptional and translational regulation by methionine restriction.
log <dig> of the fold change  of ribosomal-protected rna is plotted against log <dig> fold change of the mrna, for all the genes.
the size of the dots represents the number of reads  for the gene.
translational efficiency is measured by the ratio of ribosomal footprints fold-changes to the mrna fold-changes, quantified by a z-score
we found that genes with increased expression  are enriched for those involved in the amino acid biosynthetic process, including genes coding for enzymes for methionine biosynthetic pathway and sulfate assimilation.
those down-regulated genes  are enriched for protein synthesis  and rna methylation .
we have selected a few top repressed/induced genes and measured the corresponding protein level change upon methionine restriction using flow cytometry and gfp reporter strains.
the results are consistent with the footprint measurements
to identify the functional themes of the gene expression program, we analyzed the gene expression changes by organizing genes into functional groups with shared transcriptional regulators.
genes were grouped into transcription modules -- genes co-regulated by the same transcription factor, using the tf – target relationships previously identified by a genomic chip-chip analysis  <cit> .
we then analyzed the expression change of genes in the tf modules collectively by calculating a z-score for the whole module .
this approach allows a simpler functional organization of the transcriptome and improves the statistical power when the targets of a tf have small but coherent fold changes.
the analysis revealed that a number of tf modules are significantly up/down regulated , with clear functional themes.
tf modules down-regulated are involved in protein synthesis , cell cycle progression  and glycolysis .
tf modules up-regulated are involved in methionine biosynthesis  and general amino acid starvation response , general stress response , cellular maintenance , respiration , and iron utilization .
these observations suggest that in response to metr, cell slows down its growth and cell cycle progression and increases its stress resistance and cellular maintenance, in addition to the obvious increase of methionine pathway genes.fig.
<dig> transcription factors that play important roles in the regulation of gene expression by methionine restriction, identified by the transcription module analysis.
module z-scores measures the collective change of all the targets of the transcription factor relative to other genes .
transcription modules with z-scores >  <dig>  or < − <dig>  based on mrna changes or footprint changes were shown
interestingly, under metr cell also decreases glycolysis and increases respiration, and increased respiration was linked to lifespan extension caused by caloric restriction, suggesting that metr may also require increased respiration to extend lifespan  <cit> .
we also performed a similar analysis using kegg  <cit>  pathways.
overall the results are consistent with the go and tf module analyses.
the pathway analysis also reveals a few interesting pathways missed by the go analysis, including the induction of autophagy and ubiquitin-mediated proteolysis , suggesting that the cell tries to recycle amino acids under methionine restriction.
because of methionine’s special role in regulating translation, we are particularly interested in the subset of genes that were subjected to translational regulation by metr.
no enrichment of gene ontology categories was found in genes that increased their translational efficiency more than two-fold.
genes with decreased translational efficiency are enriched for carboxylic acid metabolic process, urea metabolic process, amino acid biosynthetic pathway  and protein synthesis  .
ribosomal genes are suppressed at the mrna level and suppressed even further at the footprint level.
similar results were also obtained from kegg pathway analysis .
potential mechanisms for translational regulation by metr
to investigate the mechanisms for translational regulation by metr, we analyzed the potential regulatory role of 5′utrs in the change of translational efficiency induced by metr.
we analyzed ribosome loading patterns by calculating the ratio of reads on the 5′ utr or 3′ utr over the coding region  for the group of genes with increased, decreased or no efficiency changes .
the ratio of reads on 5′utr vs. coding region increase dramatically for most of the genes except for those with increased translational efficiency .
the ratio for genes with increased translational efficiency  was high in the sd condition , and does not change much under metr.
the translation of gcn <dig> is known to be regulated by the 5′uorfs  <cit>  and its 5′ to coding ratio decreased drastically under metr condition.
this is consistent with the canonical model of gcn <dig> regulation  <cit>  and with the previous ribosome profiling experiments under general amino acid starvation condition  <cit> .fig.
<dig> ribosomal occupancy pattern of the 5′ utr and the coding regions for translationally regulated genes.
a. the ratio of ribosomal footprint reads in 5′utr to open reading frame under sd and metr condition, genes were grouped by their translational efficiency changes: up: efficiency change z-score > 2; down: efficiency change z-score < − <dig> and other genes.
b. scatter plot showing the fold change of reads in 5′utr vs. the fold change of reads in the open reading frame for the genes with efficiency change >  <dig>  dot size indicating the length of the 5′utr.
c. distribution of the ribosomal footprint reads on 5′utr and coding regions of gcn <dig>  xbp <dig>  and nit1
to further analyze the change of ribosome loading pattern under metr, we compared the fold change of the footprint in the 5′utr and the coding regions for each gene whose translational efficiency increased under metr .
there are three classes of genes with distinct patterns of the change of 5′utr reads vs. that of the coding sequences.
class one genes increase the footprints significantly more in their coding region compared to the 5′utr region , indicating significantly more loading of ribosome at the canonical start site.
this group includes gcn <dig> and several other genes such as nit <dig>  suggesting that they might be regulated in the similar fashion as gcn <dig>  class two genes increase the footprints more in their 5′utr region than the coding region , suggesting that the increased translational efficiency was due to ribosome loading at the non-canonical start in the 5′utr region.
thus the mechanism can be quite distinct from that for the regulation of gcn <dig>  class  <dig> genes show uniform changes in the 5′utr region and the coding region .
these results suggest that even for the genes with increased translational efficiency, there are potentially distinct regulatory mechanisms for different genes.
there is no obvious correlation between the length of the 5′utr region and the class .
translational repression of ribosome biogenesis genes by metr strongly correlated with the higher codon frequency in lysine, glutamine and glutamate suggesting translational regulation through trna thiolation
one potential mechanism through which methionine may directly regulate translation is through modulation of trna thiolation which is important for efficient translation of gene enriched in lysine , glutamine  and glutamate  codons  <cit> .
under sulfur starvation, trna thiolation will be downregulated.
if this mechanism operates under metr condition, we expect that genes enriched with kqe codons will have lower translation efficiency.
we calculated the pearson correlation between the frequency of k, q, e individually or combined with the translational efficiency changes under metr .
there is a negative correlation  between the kqe frequency and translational efficiency.
the negative correlation becomes stronger  when considering only the frequency of k. this correlation becomes even more pronounced when considering specific gene categories.
in the gene ontology analysis, we identified several groups of genes whose translational efficiency are significantly down-regulated by metr, including the ribosome biogenesis genes .
this group of genes also have a higher frequency of lysine, glutamine and glutamate codon , showing a significant negative correlation with the translational efficiency change .
this suggested that the repression of the translational efficiency of ribosome biogenesis genes may be controlled by the thiolation pathway.
for genes in other categories that are translational downregulated, there is no bias in the k, q, e codon frequency.
in addition, when excluding all the genes in the ribosome biogenesis category, there is no correlation between the codon frequency of kqe and the translational efficiency changes, indicating that the trna thiolation pathway can only explain part of the translational efficiency changes.
there is no significant correlation between the translational efficiency changes and the frequency of methionine codon.fig.
<dig> correlation between translational efficiency and the codon frequency of lysine, glutamine and glutamate  for the ribosome biogenesis genes.
a. top gene ontology categories enriched in genes with decreased translational efficiency .
b. histogram of the codon frequency of kqe in ribosome biogenesis genes compared with other genes.
c. histogram of the translational efficiency fold-change comparing the ribosome biogenesis genes with other genes.
d. ribosome biogenesis genes showed a higher kqe frequency and decreased translational efficiency in the scatter plot
discussion
translational regulation is a key step in gene regulation and plays an important role in cellular response to changing environment.
so far translational regulation has been much less well studied compared to transcriptional regulation.
the recent development of the ribosome profiling technique made it possible to study translational regulation at a fine resolution  <cit> .
as a special amino acid, methionine is coded by the translation initiation codon and methionyl trnai  is required for the assembly of the translation initiation complex  <cit> , thus sensing the cellular level of methionine may be an important mechanism for controlling translation and for coordinating the metabolic state of a cell with its growth.
in this work, we have quantified the global transcriptional and the translational programs induced by methionine restriction, using ribosome profiling and rna-seq.
we have identified hundreds of genes whose transcript level and/or translational efficiency changed significantly.
analysis of transcriptional changes based on transcription modules revealed clear functional themes.
while ribosomal genes and genes responsible for carbohydrate metabolism and cell cycle progression  are repressed, genes responsible for methionine and general amino acid synthesis, stress response, and cellular maintenance  are induced, indicating that cell slows down its growth and increases its stress resistance and maintenance/repair in response to methionine depletion.
interestingly, metr seems to induce respiration and decrease glycolysis, suggesting that the intra-cellular methionine level is coordinated with carbohydrate metabolism.
since methionine plays a key role in translational regulation, this suggests that regulation of translation is coordinated with metabolic state of the cell.
it is also worth noticing that iron utilization is increased under metr.
since methionine contains sulfur, this suggests that the cell’s response is to coordinate sulfur with iron, perhaps in making sulfur-iron clusters shown to be important for regulating lifespan  <cit> .
our analysis of genes whose translational efficiency is significantly changed by metr suggested a few mechanisms for translational regulation through methionine.
one well-studied mechanism for translational regulation is the regulation of gcn <dig> under general amino acid starvation, which involves a pathway triggered by uncharged trna.
gcn <dig> translation is regulated by several upstream uorfs.
under normal condition, only the 5′ uorfs are translated.
under amino acid starvation condition, uncharged trna activates the gcn <dig> kinase which phosphorylates the translation initiation factor eif2-alpha, leading to the translation of gcn <dig>  previous ribosomal profiling of general amino acid starvation showed high ribosome occupancy in the 5′uorf region of gcn <dig> which significantly decreases upon aa starvation, and at the mean time aa starvation induces a drastic increase of ribosome occupancy in the coding region  <cit> .
consistent with this observation, we found a similar pattern of ribosome loading at gcn <dig>  in addition, we found several other genes with ribosome loading pattern similar to gcn <dig>  suggesting that they might be regulated in the similar fashion.
interestingly, we also found a group of genes with significantly increased translational efficiency, but with opposite ribosome loading patterns .
these genes have much increased loading at their 5′utr compared to their coding regions upon metr, suggesting that the potential mechanism can be quite different from that for regulating gcn <dig>  our study provided good candidate genes/reporters for the detailed mechanistic study of translational regulation.
previously, laxman et al.
suggested that methionine can regulate translation through modulation of trna thiolation.
their study indicated that the intracellular methionine level directly controls the thiolation status of wobbleuridine  nucleotides present on lysine , glutamine , or glutamate  trnas, and that thiolated trnas lead to more efficient translation of genes enriched for kqe codons  <cit> .
our analysis of genes whose translational efficiency significantly decreased under metr lent additional support to this model.
using gene ontology  analysis, we found the translational efficiency of rrna processing genes are significantly downregulated by metr, and that this group of genes is significantly enriched for kqe codons .
while laxman et al.
study employed analysis of the proteomes of thiolation mutants, our study directly measured the translational efficiency of all genes under metr condition, providing complementary evidence supporting the thiolation model.fig.
<dig> comparison of the transcriptional and translational changes by general amino acid starvation and methionine restriction.
showing are scatter plots of the transcriptional changes , ribosomal footprint changes , translational efficiency changes , and ribosome occupancy of the 5′utr and the coding regions under amino acid starvation condition compared with the rich media
translational regulation by general amino acid starvation has been studied previously by ribosome profiling  <cit> .
we compare the gene expression profile of metr and amino acid starvation  <cit> .
overall, there is significant overlap between the transcriptional and translational changes induced by metr and amino acid starvation , which is not surprising as methionine is also restricted in the amino acid starvation.
amino acid starvation induced a stronger translational efficiency change , while only a few genes show more efficiency change in metr.
the footprint read coverage changes in the 5′utr region in amino acid starvation condition is similar to metr, showing a strongly increased ratio of reads in 5′utr over coding sequences for most of the genes.
genes with increased translational efficiency also start with a higher 5′utr read ratio which increased only marginally under amino acid starvation  similar to metr.
although these patterns are similar, the specificity of metr allowed us to infer potential regulatory mechanisms directly related to methionine.
metr is known to be able to extend the lifespan of a wide range of species.
our study identified a number genes with changed transcription and translational efficiency under metr; these genes can be good candidates for analyzing the downstream effectors of lifespan extension by metr.
for example, increased autophagy and respiration have been linked to the lifespan extension by caloric restriction, which is another well-known regimen that extends lifespan across species.
future studies based on the genes we identified should provide new insight into the mechanism of lifespan extension by metr.
CONCLUSIONS
in this work, we characterize the translational and transcriptional programs induced by metr and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast s. cerevisiae as the model organism.
using ribosomal profiling and rna-seq, we systematically compared the translational and transcriptional profiles of cells growing in the normal and methionine restricted media.
we observed a broad spectrum of gene expression changes in response to metr, including hundreds of genes whose transcript level and/or translational efficiency changed significantly.
these genes fall into specific functional classes that are informative of the physiological state of the cell under metr.
analysis of ribosome loading patterns of genes with increased translational efficiency suggested mechanisms both similar and different from the canonical model of translational regulation by general amino acid starvation.
analysis of the genes with decreased translational efficiency added support to the thiolation model of translational regulation by methionine.
since metr extends the lifespan of many species, the list of genes we identified in this study can be good candidates for studying the downstream effectors of lifespan extension.
additional files
additional file 1: sequencing reads statistics.
reads statistics for the rnaseq and ribosome footprint
additional file 2: mrna and footprint fold changes with p-values and translational efficiency under methionine restriction.
there are  <dig> spread sheets in this file:  <dig>  “metr foldchanges”: fold changes in mrna and footprint.
<dig>  “mrna_foldchange with p-val”: fold changes in mrna with p-values.
<dig>  “footpring_foldchange with p-val”: fold changes in footprint with p-values.
<dig>  “metr efficiency changes”: translational efficiency changes under methionine restriction.
additional file 3: figure s <dig>  volcano plot of the fold change and p-values.
.
transcription change under metr and the associated p-value computed from mrna data as described in the method, .
translation change under metr and the associated p-value computed from footprint data.
the p-values are provided in additional file  <dig>
additional file 4: gene ontology enrichment analysis of the genes with high footprint changes or translational efficiency changes.
there are  <dig> spread sheets in this file:  <dig>  metr footprint up  <dig> fold: enriched go biological process categories in the genes up-regulated under methionine restriction by more than four-fold in footprint level.
<dig>  metr footprint down  <dig> fold: enriched go biological process categories in the genes down-regulated under methionine restriction by more than four-fold in footprint level.
<dig>  metr efficiency up  <dig> fold: enriched go biological process categories in the genes with increased translational efficiency by at least two-fold.
<dig>  metr efficiency down  <dig> fold: enriched go biological process categories in the genes with decreased translational efficiency by at least two-fold.
additional file 5: figure s <dig>  validation of protein level changes under metr by flow cytometer.
red bars are genes with increased footprint reads while blue bars represent genes with decreased footprint reads.
the mean gfp fold changes and error bars are computed from three biological replicates.
additional file 6: module scores and kegg pathway scores under methionine restriction and amino acids starvation.
there are two spread sheets:  <dig>  module_scores: transcription factor module scores and corresponding p-values from the fold change data under methionine restriction or amino acids starvation.
<dig>  kegg pathway_scores: kegg pathway scores and corresponding p-values from the fold change data under methionine restriction or amino acids starvation.
additional file 7: translational efficiency changes under methionine restriction and the codon frequency of lysine, glutamine, glutamic acid and methionine.
this file shows the translational efficiency change for each gene and its codon frequency of lysine, glutamine, glutamic acid and methionine.
in addition, the length of the gene and the 3′utr and 5′utr length are also listed.
additional file 8: figure s <dig>  similar transcriptional and translational regulation under metr and general amino acid starvation revealed by transcription factor modules:  module scores from mrna data.
module scores from footprint data.
abbreviations
chipchromatin immunoprecipitation
kqe codonslysine , glutamine  and glutamate  codons
metrmethionine restriction
soapaligner/soap2a member of the short oligonucleotide analysis package
tsstranscription start site

 methionine is one of two sulfur-containing amino acids that are incorporated into proteins during translation. among twenty amino acids, methionine plays a special role in the biosynthesis of proteins because its codon aug is also the most common translation initiation codon. in eukaryotes, the binding of the anticodon of the initiator met-trna to the initiation codon aug is required for initiating translation. this interaction is highly conserved across species. met-trna is required for the assembly of 40s ribosome and thus may regulate the mechanism of ribosome scanning and entry, potentially serving as an important control point for translation. since translational regulation is a key step in gene regulation, sensing intracellular methionine level and adjusting the global gene expression program through translational control may be an important strategy to coordinate cell’s metabolic state with its growth. methionine has also been known to play important roles in a wild range of biological phenomena including growth, development, fertility, cancer and aging. it has been widely reported that methionine intervention can effectively regulate the lifespan of numerous model organisms. in particular, methionine restriction has been shown to extend the lifespan of a range of species, including yeast, worm, fly and mouse. it has also been suggested that the lifespan extension by caloric restriction, defined as reduced caloric intake without malnutrition, can be attributed to methionine restriction. in addition to the effect on lifespan, methionine restriction also slows or reduces many characteristics associated with senescence, such as immune and lens aging, increased igf-i and insulin levels, and cumulated oxidative damages. methionine restriction has also been studied extensively in anticancer therapies, either alone or in association with the other treatments, and is considered as a useful therapeutic strategy for treating various cancers. thus, characterizing the global gene expression program induced by metr and understanding the mechanisms by which metr regulates gene expression are important not only for understanding the basic principles of gene regulation but also for promoting human health. translational regulation by general amino acid starvation has been extensively studied and the pathway involved has been elucidated before. in the canonical model, amino acid starvation leads to the accumulation of uncharged trna, activating the gcn kinase, which phosphorylates eif, ultimately affecting the translation. as a general strategy for sensing amino acid depletion, this may also be the mechanism to sense and respond to metr. methionine may also work through other mechanisms to affect translation. it has been reported that intracellular methionine availability can regulate cellular translational capacity and metabolic homeostasis by controlling the thiolation status of the wobble-uridine nucleotides on lysine, glutamine, or glutamate trnas. methionine may also affect gene expression by converting to s-adenosyl methionine, which serves as the predominant methyl donor for rrna-trna modifications and dna/protein methylations. although there has been significant progress in understanding the roles methionine may play in gene regulation, a systematic study on the global gene expression program controlled by methionine, especially at the translational level, is still lacking. in this work, we use ribosomal profiling and rna-seq to compare the translational and transcriptional profiles of cells growing in the normal and methionine restricted media. we systematically characterize the translational and transcriptional programs induced by methionine restriction and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast s. cerevisiae as the model organism. methods yeast strains and media yeast strains used for the ribosomal profiling and rna-seq experiments were by. synthetic dextrose medium and methionine restriction medium was used in the ribosome profiling/rna-seq experiments, the sd medium contained 2% glucose,  g/l yeast nitrogen base without amino acid, 20 mg/l adenine, 20 mg/l l-arginine hcl, 100 mg/l l-aspartic acid, 20 mg/l l-histidine hcl, 100 mg/l l-leucine, 30 mg/l l-isoleucine, 30 mg/l l-lysine hcl, 20 mg/l l-methionine, 50 mg/l l-phenylalanine, 200 mg/l l-threonine, 20 mg/l l-tryptophan, 30 mg/l l-tyrosine, 20 mg/l uracil, 150 mg/l l-valine, 100 mg/l glutamic acid and 4 g/l serine. the metr media has the similar ingredients as the sd medium except for 4 mg/l l-methionine concentration. media were freshly made before the experiments. all nutrients were purchased from sigma-aldrich corporation. ribosome profiling and rna-seq of cells growing in sd vs. metr media the initial cell culture was incubated in 300 ml sd medium overnight to an od  ~ , then diluted by five fold using fresh sd media and incubated for another 4 h under °c to an od  ~ . the sample was then divided equally into two aliquots. cells were separated from the media by spin-down at 3000 g for 5 min and re-suspended in sd and metr media respectively. all samples were incubated for another hour before harvesting. all the steps were carried out at 30 °c. ribosomal profiling and rna-seq experiments were carried out using the protocol developed by ingolia et al. raw sequences were obtained from illumina hiseq sequence analysis and quantification of differential gene expression sequence reads were aligned to the most recent s. cerevisiae genome using soapaligner/soap with default setting. after trimming off the adapters, reads aligned to rrna and trna sequences were filtered out. the rest of the reads were then aligned to the genome sequence. finally, reads that did not align to the genome were aligned to all the cdss to retain those covering the splicing junctions. after the alignment, we counted the number of reads starting at each position across the whole genome; for ribosomal footprinting data, the starting position of each read was shifted by 15 bps towards the 3′ end, to adjust for the offset due to ribosome protection. to get the abundance of reads covering each gene, we sum all reads with starting position from the start to the stop codon, excluding the first 50 bps from the transcription start site, to alleviate the effect of the biased distribution of reads around the tss. for a reliable gene expression comparison between two conditions, we excluded genes with less than total raw reads. we computed the fold change of mrna or footprint as the ratio of the corresponding reads from the two conditions, with total reads normalized to adjust the median fold change to to estimate the statistical significance of the fold change, we did the following analysis: for each gene g i, we collected the other genes with the most similar number of reads and compute the standard deviation σ i from the log of these genes. we then computed the zi=logfiσi, where f i is the fold change of g i comparing metr and sd. then a p-value was calculated from the z i to measure the significance of the fold change. quantification of translational efficiency changes the translational efficiency changes were calculated as the ratio of ribosomal footprints fold change to mrna fold changes for each gene. translational efficiency change z-score is calculated by normalized the efficiency change with the standard deviation. calculation of the tf module z-scores and kegg pathway z-scores we used the transcription factor targets from the analysis by mcissac et al. based on the systematic chip-chip data, using as the p-value cutoff and the strongest conservation between species. we computed the rank sum test z-scores comparing the fold changes of the target genes vs. none-target genes for each transcription factor. the sign of the z-score reflects the overall direction of the gene expression change in the modules; positive z-score indicates overall induction and negative z-score indicates overall repression. we used tf modules with at least targets for this analysis. we used a similar method to compute the kegg pathway z-scores by grouping the genes from the same kegg pathway. flow cytometer measurement of the protein abundance changes upon methionine restriction yeast gfp-tag strains were selected from the yeast gfp library. for each gfp strain, we picked three single clones from the plate and cultured them overnight to saturation. we then diluted each cell culture on the second day and grew them to od – in a well plate. cells were collected and re-suspended in metr or sd media. we then used flow cytometer to measure the gfp signals in each sample after 4 h. the cellular gfp concentration was computed by normalizing the gfp signal with the cell size using forward scattering signal for each individual cell. the difference of gfp concentration between metr and sd was computed as the ratio of the medians of the normalized gfp under the two conditions. then the mean gfp fold change was calculated from the three biological replicates.  ribosomal profiling and rna-seq revealed a broad spectrum of transcriptional and translational changes induced by metr we performed global gene expression profiling by rna-seq and ribosomal profiling, comparing cells growing in normal vs. metr conditions. the sequencing result is of high quality with at least 50× coverage per sample. ribosome profiling quantifies ribosome protected rna for all the genes in the genome and thus can measure the translational efficiency when combined with the total amount rna from the rna-seq data. we observed a broad spectrum of gene expression changes in response to metr, both at the transcriptional and the translational levels. for the majority of the genes that changed expression, the regulation is at the transcriptional level, as the fold change of the footprint is proportional to the fold change of the transcript level. there is a subset of genes whose translational efficiency are increased or decreased compared to all the genes, indicating that they are under translational regulation. we observed that a subset of genes with decreased transcriptional level tend to have a decreased translational efficiency, suggesting that they are under both transcriptional and translational control. overall there are 110/ genes whose footprints went up/down by more than four-fold under metr, and 149/ genes increased/decreased their translational efficiency by more than two-fold.fig. rna-seq and ribosomal profiling revealed global transcriptional and translational regulation by methionine restriction. log of the fold change of ribosomal-protected rna is plotted against log fold change of the mrna, for all the genes. the size of the dots represents the number of reads for the gene. translational efficiency is measured by the ratio of ribosomal footprints fold-changes to the mrna fold-changes, quantified by a z-score we found that genes with increased expression are enriched for those involved in the amino acid biosynthetic process, including genes coding for enzymes for methionine biosynthetic pathway and sulfate assimilation. those down-regulated genes are enriched for protein synthesis and rna methylation. we have selected a few top repressed/induced genes and measured the corresponding protein level change upon methionine restriction using flow cytometry and gfp reporter strains. the results are consistent with the footprint measurements to identify the functional themes of the gene expression program, we analyzed the gene expression changes by organizing genes into functional groups with shared transcriptional regulators. genes were grouped into transcription modules -- genes co-regulated by the same transcription factor, using the tf – target relationships previously identified by a genomic chip-chip analysis. we then analyzed the expression change of genes in the tf modules collectively by calculating a z-score for the whole module. this approach allows a simpler functional organization of the transcriptome and improves the statistical power when the targets of a tf have small but coherent fold changes. the analysis revealed that a number of tf modules are significantly up/down regulated, with clear functional themes. tf modules down-regulated are involved in protein synthesis, cell cycle progression and glycolysis. tf modules up-regulated are involved in methionine biosynthesis and general amino acid starvation response, general stress response, cellular maintenance, respiration, and iron utilization. these observations suggest that in response to metr, cell slows down its growth and cell cycle progression and increases its stress resistance and cellular maintenance, in addition to the obvious increase of methionine pathway genes.fig. transcription factors that play important roles in the regulation of gene expression by methionine restriction, identified by the transcription module analysis. module z-scores measures the collective change of all the targets of the transcription factor relative to other genes. transcription modules with z-scores >  or < − based on mrna changes or footprint changes were shown interestingly, under metr cell also decreases glycolysis and increases respiration, and increased respiration was linked to lifespan extension caused by caloric restriction, suggesting that metr may also require increased respiration to extend lifespan. we also performed a similar analysis using kegg pathways. overall the results are consistent with the go and tf module analyses. the pathway analysis also reveals a few interesting pathways missed by the go analysis, including the induction of autophagy and ubiquitin-mediated proteolysis, suggesting that the cell tries to recycle amino acids under methionine restriction. because of methionine’s special role in regulating translation, we are particularly interested in the subset of genes that were subjected to translational regulation by metr. no enrichment of gene ontology categories was found in genes that increased their translational efficiency more than two-fold. genes with decreased translational efficiency are enriched for carboxylic acid metabolic process, urea metabolic process, amino acid biosynthetic pathway and protein synthesis. ribosomal genes are suppressed at the mrna level and suppressed even further at the footprint level. similar results were also obtained from kegg pathway analysis. potential mechanisms for translational regulation by metr to investigate the mechanisms for translational regulation by metr, we analyzed the potential regulatory role of 5′utrs in the change of translational efficiency induced by metr. we analyzed ribosome loading patterns by calculating the ratio of reads on the 5′ utr or 3′ utr over the coding region for the group of genes with increased, decreased or no efficiency changes. the ratio of reads on 5′utr vs. coding region increase dramatically for most of the genes except for those with increased translational efficiency. the ratio for genes with increased translational efficiency was high in the sd condition, and does not change much under metr. the translation of gcn is known to be regulated by the 5′uorfs and its 5′ to coding ratio decreased drastically under metr condition. this is consistent with the canonical model of gcn regulation and with the previous ribosome profiling experiments under general amino acid starvation condition.fig. ribosomal occupancy pattern of the 5′ utr and the coding regions for translationally regulated genes. a. the ratio of ribosomal footprint reads in 5′utr to open reading frame under sd and metr condition, genes were grouped by their translational efficiency changes: up: efficiency change z-score > 2; down: efficiency change z-score < − and other genes. b. scatter plot showing the fold change of reads in 5′utr vs. the fold change of reads in the open reading frame for the genes with efficiency change >  dot size indicating the length of the 5′utr. c. distribution of the ribosomal footprint reads on 5′utr and coding regions of gcn xbp and nit1 to further analyze the change of ribosome loading pattern under metr, we compared the fold change of the footprint in the 5′utr and the coding regions for each gene whose translational efficiency increased under metr. there are three classes of genes with distinct patterns of the change of 5′utr reads vs. that of the coding sequences. class one genes increase the footprints significantly more in their coding region compared to the 5′utr region, indicating significantly more loading of ribosome at the canonical start site. this group includes gcn and several other genes such as nit suggesting that they might be regulated in the similar fashion as gcn class two genes increase the footprints more in their 5′utr region than the coding region, suggesting that the increased translational efficiency was due to ribosome loading at the non-canonical start in the 5′utr region. thus the mechanism can be quite distinct from that for the regulation of gcn class genes show uniform changes in the 5′utr region and the coding region. these results suggest that even for the genes with increased translational efficiency, there are potentially distinct regulatory mechanisms for different genes. there is no obvious correlation between the length of the 5′utr region and the class. translational repression of ribosome biogenesis genes by metr strongly correlated with the higher codon frequency in lysine, glutamine and glutamate suggesting translational regulation through trna thiolation one potential mechanism through which methionine may directly regulate translation is through modulation of trna thiolation which is important for efficient translation of gene enriched in lysine, glutamine and glutamate codons. under sulfur starvation, trna thiolation will be downregulated. if this mechanism operates under metr condition, we expect that genes enriched with kqe codons will have lower translation efficiency. we calculated the pearson correlation between the frequency of k, q, e individually or combined with the translational efficiency changes under metr. there is a negative correlation between the kqe frequency and translational efficiency. the negative correlation becomes stronger when considering only the frequency of k. this correlation becomes even more pronounced when considering specific gene categories. in the gene ontology analysis, we identified several groups of genes whose translational efficiency are significantly down-regulated by metr, including the ribosome biogenesis genes. this group of genes also have a higher frequency of lysine, glutamine and glutamate codon, showing a significant negative correlation with the translational efficiency change. this suggested that the repression of the translational efficiency of ribosome biogenesis genes may be controlled by the thiolation pathway. for genes in other categories that are translational downregulated, there is no bias in the k, q, e codon frequency. in addition, when excluding all the genes in the ribosome biogenesis category, there is no correlation between the codon frequency of kqe and the translational efficiency changes, indicating that the trna thiolation pathway can only explain part of the translational efficiency changes. there is no significant correlation between the translational efficiency changes and the frequency of methionine codon.fig. correlation between translational efficiency and the codon frequency of lysine, glutamine and glutamate for the ribosome biogenesis genes. a. top gene ontology categories enriched in genes with decreased translational efficiency. b. histogram of the codon frequency of kqe in ribosome biogenesis genes compared with other genes. c. histogram of the translational efficiency fold-change comparing the ribosome biogenesis genes with other genes. d. ribosome biogenesis genes showed a higher kqe frequency and decreased translational efficiency in the scatter plot discussion translational regulation is a key step in gene regulation and plays an important role in cellular response to changing environment. so far translational regulation has been much less well studied compared to transcriptional regulation. the recent development of the ribosome profiling technique made it possible to study translational regulation at a fine resolution. as a special amino acid, methionine is coded by the translation initiation codon and methionyl trnai is required for the assembly of the translation initiation complex, thus sensing the cellular level of methionine may be an important mechanism for controlling translation and for coordinating the metabolic state of a cell with its growth. in this work, we have quantified the global transcriptional and the translational programs induced by methionine restriction, using ribosome profiling and rna-seq. we have identified hundreds of genes whose transcript level and/or translational efficiency changed significantly. analysis of transcriptional changes based on transcription modules revealed clear functional themes. while ribosomal genes and genes responsible for carbohydrate metabolism and cell cycle progression are repressed, genes responsible for methionine and general amino acid synthesis, stress response, and cellular maintenance are induced, indicating that cell slows down its growth and increases its stress resistance and maintenance/repair in response to methionine depletion. interestingly, metr seems to induce respiration and decrease glycolysis, suggesting that the intra-cellular methionine level is coordinated with carbohydrate metabolism. since methionine plays a key role in translational regulation, this suggests that regulation of translation is coordinated with metabolic state of the cell. it is also worth noticing that iron utilization is increased under metr. since methionine contains sulfur, this suggests that the cell’s response is to coordinate sulfur with iron, perhaps in making sulfur-iron clusters shown to be important for regulating lifespan. our analysis of genes whose translational efficiency is significantly changed by metr suggested a few mechanisms for translational regulation through methionine. one well-studied mechanism for translational regulation is the regulation of gcn under general amino acid starvation, which involves a pathway triggered by uncharged trna. gcn translation is regulated by several upstream uorfs. under normal condition, only the 5′ uorfs are translated. under amino acid starvation condition, uncharged trna activates the gcn kinase which phosphorylates the translation initiation factor eif2-alpha, leading to the translation of gcn previous ribosomal profiling of general amino acid starvation showed high ribosome occupancy in the 5′uorf region of gcn which significantly decreases upon aa starvation, and at the mean time aa starvation induces a drastic increase of ribosome occupancy in the coding region. consistent with this observation, we found a similar pattern of ribosome loading at gcn in addition, we found several other genes with ribosome loading pattern similar to gcn suggesting that they might be regulated in the similar fashion. interestingly, we also found a group of genes with significantly increased translational efficiency, but with opposite ribosome loading patterns. these genes have much increased loading at their 5′utr compared to their coding regions upon metr, suggesting that the potential mechanism can be quite different from that for regulating gcn our study provided good candidate genes/reporters for the detailed mechanistic study of translational regulation. previously, laxman et al. suggested that methionine can regulate translation through modulation of trna thiolation. their study indicated that the intracellular methionine level directly controls the thiolation status of wobbleuridine nucleotides present on lysine, glutamine, or glutamate trnas, and that thiolated trnas lead to more efficient translation of genes enriched for kqe codons. our analysis of genes whose translational efficiency significantly decreased under metr lent additional support to this model. using gene ontology analysis, we found the translational efficiency of rrna processing genes are significantly downregulated by metr, and that this group of genes is significantly enriched for kqe codons. while laxman et al. study employed analysis of the proteomes of thiolation mutants, our study directly measured the translational efficiency of all genes under metr condition, providing complementary evidence supporting the thiolation model.fig. comparison of the transcriptional and translational changes by general amino acid starvation and methionine restriction. showing are scatter plots of the transcriptional changes, ribosomal footprint changes, translational efficiency changes, and ribosome occupancy of the 5′utr and the coding regions under amino acid starvation condition compared with the rich media translational regulation by general amino acid starvation has been studied previously by ribosome profiling. we compare the gene expression profile of metr and amino acid starvation. overall, there is significant overlap between the transcriptional and translational changes induced by metr and amino acid starvation, which is not surprising as methionine is also restricted in the amino acid starvation. amino acid starvation induced a stronger translational efficiency change, while only a few genes show more efficiency change in metr. the footprint read coverage changes in the 5′utr region in amino acid starvation condition is similar to metr, showing a strongly increased ratio of reads in 5′utr over coding sequences for most of the genes. genes with increased translational efficiency also start with a higher 5′utr read ratio which increased only marginally under amino acid starvation similar to metr. although these patterns are similar, the specificity of metr allowed us to infer potential regulatory mechanisms directly related to methionine. metr is known to be able to extend the lifespan of a wide range of species. our study identified a number genes with changed transcription and translational efficiency under metr; these genes can be good candidates for analyzing the downstream effectors of lifespan extension by metr. for example, increased autophagy and respiration have been linked to the lifespan extension by caloric restriction, which is another well-known regimen that extends lifespan across species. future studies based on the genes we identified should provide new insight into the mechanism of lifespan extension by metr.  in this work, we characterize the translational and transcriptional programs induced by metr and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast s. cerevisiae as the model organism. using ribosomal profiling and rna-seq, we systematically compared the translational and transcriptional profiles of cells growing in the normal and methionine restricted media. we observed a broad spectrum of gene expression changes in response to metr, including hundreds of genes whose transcript level and/or translational efficiency changed significantly. these genes fall into specific functional classes that are informative of the physiological state of the cell under metr. analysis of ribosome loading patterns of genes with increased translational efficiency suggested mechanisms both similar and different from the canonical model of translational regulation by general amino acid starvation. analysis of the genes with decreased translational efficiency added support to the thiolation model of translational regulation by methionine. since metr extends the lifespan of many species, the list of genes we identified in this study can be good candidates for studying the downstream effectors of lifespan extension. additional files additional file 1: sequencing reads statistics. reads statistics for the rnaseq and ribosome footprint additional file 2: mrna and footprint fold changes with p-values and translational efficiency under methionine restriction. there are spread sheets in this file: “metr foldchanges”: fold changes in mrna and footprint. “mrna_foldchange with p-val”: fold changes in mrna with p-values. “footpring_foldchange with p-val”: fold changes in footprint with p-values. “metr efficiency changes”: translational efficiency changes under methionine restriction. additional file 3: figure s volcano plot of the fold change and p-values.. transcription change under metr and the associated p-value computed from mrna data as described in the method,. translation change under metr and the associated p-value computed from footprint data. the p-values are provided in additional file additional file 4: gene ontology enrichment analysis of the genes with high footprint changes or translational efficiency changes. there are spread sheets in this file: metr footprint up fold: enriched go biological process categories in the genes up-regulated under methionine restriction by more than four-fold in footprint level. metr footprint down fold: enriched go biological process categories in the genes down-regulated under methionine restriction by more than four-fold in footprint level. metr efficiency up fold: enriched go biological process categories in the genes with increased translational efficiency by at least two-fold. metr efficiency down fold: enriched go biological process categories in the genes with decreased translational efficiency by at least two-fold. additional file 5: figure s validation of protein level changes under metr by flow cytometer. red bars are genes with increased footprint reads while blue bars represent genes with decreased footprint reads. the mean gfp fold changes and error bars are computed from three biological replicates. additional file 6: module scores and kegg pathway scores under methionine restriction and amino acids starvation. there are two spread sheets: module_scores: transcription factor module scores and corresponding p-values from the fold change data under methionine restriction or amino acids starvation. kegg pathway_scores: kegg pathway scores and corresponding p-values from the fold change data under methionine restriction or amino acids starvation. additional file 7: translational efficiency changes under methionine restriction and the codon frequency of lysine, glutamine, glutamic acid and methionine. this file shows the translational efficiency change for each gene and its codon frequency of lysine, glutamine, glutamic acid and methionine. in addition, the length of the gene and the 3′utr and 5′utr length are also listed. additional file 8: figure s similar transcriptional and translational regulation under metr and general amino acid starvation revealed by transcription factor modules: module scores from mrna data. module scores from footprint data. abbreviations chipchromatin immunoprecipitation kqe codonslysine, glutamine and glutamate codons metrmethionine restriction soapaligner/soap2a member of the short oligonucleotide analysis package tsstranscription start site 
among twenty amino acids, methionine has a special role as it is coded by the translation initiation codon and methionyl-trnai is required for the assembly of the translation initiation complex. thus methionine may play a special role in global gene regulation. these genes show clear functional themes, suggesting that cell slows down its growth and cell cycle progression and increases its stress resistance and maintenance in response to metr. using ribosomal profiling and rna-seq, we observed a broad spectrum of gene expression changes in response to metr and identified hundreds of genes whose transcript level and/or translational efficiency changed significantly. analysis of genes whose translational efficiency changed significantly under metr revealed different modes of translational regulation: 1) ribosome loading patterns in the 5′utr and coding regions of genes with increased translational efficiency suggested mechanisms both similar and different from that for the translational regulation of gcn under general amino acid starvation condition; 2) genes with decreased translational efficiency showed strong enrichment of lysine, glutamine, and glutamate codons, supporting the model that methionine can regulate translation by controlling trna thiolation. in this work, we characterize the translational and transcriptional programs induced by methionine restriction and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast s. cerevisiae as the model organism. metr induced a broad spectrum of gene expression changes at both the transcriptional and translational levels, with clear functional themes indicative of the physiological state of the cell under metr. different modes of translational regulation were induced by metr, including the regulation of the ribosome loading at 5′utr and regulation by trna thiolation. since metr extends the lifespan of many species, the list of genes we identified in this study can be good candidates for studying the mechanisms of lifespan extension. interestingly, under metr cell also decreases glycolysis and increases respiration, and increased respiration was linked to lifespan extension caused by caloric restriction. methionine has also been known to play important roles in cell growth, development, cancer, and aging.
 among twenty amino acids, methionine has a special role as it is coded by the translation initiation codon and methionyl-trnai is required for the assembly of the translation initiation complex. thus methionine may play a special role in global gene regulation. methionine has also been known to play important roles in cell growth, development, cancer, and aging. in this work, we characterize the translational and transcriptional programs induced by methionine restriction and investigate the potential mechanisms through which methionine regulates gene expression, using the budding yeast s. cerevisiae as the model organism.  using ribosomal profiling and rna-seq, we observed a broad spectrum of gene expression changes in response to metr and identified hundreds of genes whose transcript level and/or translational efficiency changed significantly. these genes show clear functional themes, suggesting that cell slows down its growth and cell cycle progression and increases its stress resistance and maintenance in response to metr. interestingly, under metr cell also decreases glycolysis and increases respiration, and increased respiration was linked to lifespan extension caused by caloric restriction. analysis of genes whose translational efficiency changed significantly under metr revealed different modes of translational regulation: 1) ribosome loading patterns in the 5′utr and coding regions of genes with increased translational efficiency suggested mechanisms both similar and different from that for the translational regulation of gcn under general amino acid starvation condition; 2) genes with decreased translational efficiency showed strong enrichment of lysine, glutamine, and glutamate codons, supporting the model that methionine can regulate translation by controlling trna thiolation.  metr induced a broad spectrum of gene expression changes at both the transcriptional and translational levels, with clear functional themes indicative of the physiological state of the cell under metr. different modes of translational regulation were induced by metr, including the regulation of the ribosome loading at 5′utr and regulation by trna thiolation. since metr extends the lifespan of many species, the list of genes we identified in this study can be good candidates for studying the mechanisms of lifespan extension. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. http://dx.doi.org/ /100000002national institutes of healthag043080li hao issue-copyright-statement© the author 2017 