Original Research ARTICLE
The Arabidopsis TOR Kinase Specifically Regulates the Expression of Nuclear Genes Coding for Plastidic Ribosomal Proteins and the Phosphorylation of the Cytosolic Ribosomal Protein S6
- 1Institut Jean-Pierre Bourgin, Institut National de la Recherche Agronomique, AgroParisTech, Centre National de la Recherche Scientifique, Université Paris-Saclay, Versailles, France
- 2Université Paris-Sud–Université Paris-Saclay, Orsay, France
- 3Umeå Plant Science Center, Department of Plant Physiology, Umeå University, Umeå, Sweden
- 4Institut de Biologie Moléculaire des Plantes, UPR 2357 CNRS, Université de Strasbourg, Strasbourg, France
- 5Plateforme PAPPSO, UMR GQE-Le Moulon, Gif sur Yvette, France
- 6Laboratoire de Génétique et Biophysique des Plantes, UMR 7265, DSV, IBEB, SBVME, CEA, CNRS, Aix-Marseille Université, Faculté des Sciences de Luminy, Marseille, France
- 7Plateforme Protéomique Strasbourg-Esplanade, CNRS FRC1589, Institut de Biologie Moléculaire et Cellulaire, Strasbourg, France
Protein translation is an energy consuming process that has to be fine-tuned at both the cell and organism levels to match the availability of resources. The target of rapamycin kinase (TOR) is a key regulator of a large range of biological processes in response to environmental cues. In this study, we have investigated the effects of TOR inactivation on the expression and regulation of Arabidopsis ribosomal proteins at different levels of analysis, namely from transcriptomic to phosphoproteomic. TOR inactivation resulted in a coordinated down-regulation of the transcription and translation of nuclear-encoded mRNAs coding for plastidic ribosomal proteins, which could explain the chlorotic phenotype of the TOR silenced plants. We have identified in the 5′ untranslated regions (UTRs) of this set of genes a conserved sequence related to the 5′ terminal oligopyrimidine motif, which is known to confer translational regulation by the TOR kinase in other eukaryotes. Furthermore, the phosphoproteomic analysis of the ribosomal fraction following TOR inactivation revealed a lower phosphorylation of the conserved Ser240 residue in the C-terminal region of the 40S ribosomal protein S6 (RPS6). These results were confirmed by Western blot analysis using an antibody that specifically recognizes phosphorylated Ser240 in RPS6. Finally, this antibody was used to follow TOR activity in plants. Our results thus uncover a multi-level regulation of plant ribosomal genes and proteins by the TOR kinase.
During their life, living organisms have to adapt their growth and development to exogenous factors such as stresses and nutrient availability. Therefore, they have evolved different regulatory pathways to increase the perception of environmental cues and to fasten the required metabolic modifications. These pathways employ conserved key players that link energy depletion, which is often the result of stresses and nutrient limitation, to anabolic and catabolic cellular activities. One of the most important pathway that is found in all eukaryotes is the one related to the target of rapamycin (TOR) protein kinase. TOR is a large kinase, which operates in at least two multi-protein complexes (TORC1 and TORC2; for reviews, see: Wullschleger et al., 2006; Laplante and Sabatini, 2012; Albert and Hall, 2015) and controls a wealth of biological outputs. In animals and yeast, it is well known that TOR positively regulates protein synthesis and anabolic activities when the growth conditions are favorable, while repressing the mechanisms implicated in recycling and catabolism (Laplante and Sabatini, 2012; Shimobayashi and Hall, 2014). Indeed, the production of proteins is particularly energy consuming since it requires ribosome biogenesis as well as mRNA translation (Warner, 1999).
In plants, there is so far only evidence for the presence of the TORC1 complex which comprises the conserved Regulatory-associated protein of TOR (RAPTOR) and the Lethal with Sec 13 (LST8) proteins, (for reviews, see: Robaglia et al., 2012; Henriques et al., 2014; Xiong and Sheen, 2014; Rexin et al., 2015; Dobrenel et al., 2016). TOR has already been shown to control a vast array of biological processes in plants (Deprost et al., 2007; Caldana et al., 2013; Xiong and Sheen, 2014; Dong et al., 2015) and a link between the TOR complex and protein translation has been evidenced. We have indeed shown earlier that TOR inactivation, either after silencing (Deprost et al., 2007) or by using a TOR inhibitor (Sormani et al., 2007) leads to a decrease in polysome abundance. It has also been shown that the translation reinitiation after a long upstream open reading frame by the plant viral reinitiation factor transactivator-viroplasmin is mediated by its physical association with the TOR protein (Schepetilnikov et al., 2011). Furthermore, TOR activity appears essential for translation reinitiation of cellular mRNA containing short-ORFs in the 5′ UTR (Schepetilnikov et al., 2013).
The biochemical analysis of the plant cytoplasmic ribosome showed that it contains 81 different proteins, 33 for the small 40S subunit and 48 for the large 80S subunit (Giavalisco et al., 2005; Carroll et al., 2008; Carroll, 2013; Browning and Bailey-Serres, 2015; Hummel et al., 2015). In animal cells, TOR regulates cap-dependent translation by phosphorylating and stimulating the activity of the ribosomal protein S6 kinase (S6K), a conserved target of TOR which phosphorylates the 40S ribosomal protein S6 (RPS6, Barbet et al., 1996; Holz et al., 2005), and by repressing the inhibitory effect of eIF4E-binding protein (Ma and Blenis, 2009; Albert and Hall, 2015). Consistently, in yeast, TOR inhibition by rapamycin leads to an 80% reduction in overall translation (Barbet et al., 1996). It has been shown in mammals that S6K is activated in a TOR-dependent manner by phosphorylation of Thr389 and Thr229 (Ma and Blenis, 2009), resulting subsequently in the phosphorylation of serine residues in the C-terminal extremity of the ribosomal protein RPS6. Early on RPS6, which is located at the right foot of the 40S subunit, was identified as the only phosphorylated protein in the ribosome small subunit (Gressner and Wool, 1976).
It has been postulated that TOR regulates the translation of a particular sub-set of mRNAs containing a 5′ terminal tract oligopyrimidine (TOP) motif (for a review, see Meyuhas and Kahan, 2015). Canonical TOP mRNAs harbor a C residue on position 1 followed by a stretch of 4–15 pyrimidines. The first evidence suggested that TOR activates TOP mRNA translation through phosphorylation of S6K and RPS6 but this hypothesis was later questioned since TOP mRNA are normally translated in S6K-deficient mice (Meyuhas and Kahan, 2015). The Arabidopsis genome contains two tandem-repeated S6K genes and the proteins encoded by these genes are directly phosphorylated by TOR (Mahfouz et al., 2006; Schepetilnikov et al., 2011, 2013; Xiong et al., 2013). The phosphorylation level of S6K proteins is positively correlated to their capacity to phosphorylate RPS6 (Turck et al., 1998, 2004; Mahfouz et al., 2006; Browning and Bailey-Serres, 2015). In yeast, untargeted phosphoproteomic analyses pointed RPS6 to be the main phosphorylation target of TOR in the ribosome (Huber et al., 2009). In plants, RPS6 was found to be the major phosphorylated ribosomal protein in tomato cells and this phosphorylation was found to be reduced after heat stress (Scharf and Nover, 1982) or in oxygen deprived maize roots (Bailey-Serres and Freeling, 1990). Later a survey of post-translational modifications of the Arabidopsis ribosomal proteins only identified Ser240 in RPS6 and Ser137 in RPL13, together with acidic proteins, as being phosphorylated (Carroll et al., 2008). Multiple phosphorylation sites were detected in the RPS6 C-terminal region including Ser238 and Ser241 for maize (Williams et al., 2003) and Ser237 and Ser240 for Arabidopsis (Chang et al., 2005). Moreover Ser240 phosphorylation was found to be induced by light and high CO2 conditions (Turkina et al., 2011; Boex-Fontvieille et al., 2013). More recently, Nukarinen et al. (2016) have observed a strong induction of RPS6 Ser240 phosphorylation when the activity of the Sucrose non-fermenting 1-Related Kinase 1 (SnRK1) is decreased. However, despite the accumulation of data showing variations in plant RPS6 phosphorylation in response to several stresses, the precise role of this C-terminal phosphorylation in the regulation of translation remains largely elusive.
Since TOR was found to affect translation in plants, we undertook a global phosphoproteomic, transcriptomic and translatomic analysis of the ribosomal fraction after TOR inactivation. Interestingly, we observed a strong effect of TOR inactivation on the expression of nuclear-encoded plastidic ribosomal proteins (pRPs) and the main phosphorylation site controlled by TOR activity was found to be Ser240 in the cytoplasmic ribosomal protein (cRP) RPS6. Finally, we made use of this specific phosphorylation site to design a robust Western-based method for quantifying TOR activity in plant extracts.
Materials and Methods
Plant Materials and Growth Conditions
Seeds of two independent ethanol-inducible TOR RNAi lines (5.2 and 6.3, described in Deprost et al., 2007) as well as an ethanol-inducible GUS overexpressing line (as a control) (Deprost et al., 2007) were grown in vitro under long day conditions (16 h light/8 h night) for 7 days on solid 1/5 Murashige and Skoog medium supplemented with sucrose 0.3% (w/v) at a constant temperature of 25°C and a light intensity of 75 μE.m-2.s-1. The plants were subsequently treated with ethanol vapor for either 3 or 10 days. Whole plantlets from two independent biological replicates of each condition were then harvested in the middle of the light period and directly snap frozen in liquid nitrogen, grinded and subjected immediately to the ribosome enrichment protocol.
Ribosomal subunits (40S and 60S), monoribosomes (80S) and polyribosomes were isolated from the plantlet powder according to Bailey-Serres and Freeling (1990) with minor modifications. Freshly harvested and grinded plantlets were homogenized at a final concentration of 10% (w/v) in the ice-cold extraction buffer (0.2 M Tris-HCl [pH 9], 0.4 M KCl, 0.025 M EGTA, 0.035 M MgCl2, 0.2 M sucrose) supplemented with 2% (v/v) Triton X-100, 2% (v/v) Tween 20, 2% (v/v) NP-40 and 1% (w/v) sodium deoxycholate. The extracts were incubated on ice for 10 min to solubilize membrane-bound ribosomes and centrifuged at 2880 × g for 15 min at 4°C. The supernatants were layered over a sucrose cushion (0.04 M Tris-HCl [pH 9], 0.2 M KCl, 0.005 M EGTA, 0.03 M MgCl2, 1.75 M sucrose) and ultracentrifuged at 225 000 × g for 14 h. The ribosome enriched pellet was resuspended in 300 μl of Laemmli buffer (Laemmli, 1970) and denatured at 100°C for 10 min.
For the proteomic characterization, ribosome enriched fractions were first submitted to a short migration through the stacking gel of a SDS-PAGE, in order to remove the rRNA and the possible chemical contaminant, including detergents. After a Coomassie staining, the unique band of proteins, for each sample, was cut and divided into five pieces that were submitted, in gel, to the tryptic digestion, reduction and alkylation. Peptide containing fractions were then analyzed by nano LC-MS/MS as previously described (Boex-Fontvieille et al., 2013). Briefly on-line liquid chromatography was performed on a NanoLC-Ultra system (Eksigent). Eluted peptides were analyzed with a Q-Exactive mass spectrometer (Thermo Electron) using a nano-electrospray interface (non-coated capillary probe, 10 μ i.d; New Objective). Peptides and the corresponding proteins were identified and grouped with X!TandemPipeline using the X!Tandem Piledriver (2015.04.01) release (Craig and Beavis, 2004) and the TAIR10 protein library with the phosphorylation of serine, threonine and tyrosine as a potential peptide modification. Precursor mass tolerance was 10 ppm and fragment mass tolerance was 0.02 Th. Identified proteins were filtered and grouped using the X!TandemPipeline v3.3.41. Data filtering was achieved according to a peptide E-value lower than 0.01. The false discovery rate (FDR) was estimated to 0.92%. Relative quantification was performed using the MassChroQ software (Valot et al., 2011) by peak area integration on extracted ion chromatograms (XICs) within a 10 ppm window, after LC-MS/MS chromatogram alignment and spike filtering.
Arabidopsis seedlings grown on MS agar plates in standard 16/8 h and 21/17°C day/night conditions were transferred to liquid MS media supplemented with 10 μM NAA (Sigma-Aldrich). Total protein extracts were precipitated with 0.1 M ammonium acetate in 100% methanol, reduced, alkylated and digested overnight with trypsin (Promega, Madison, WI, USA) in 50 mM ammonium bicarbonate. Resulting peptides were vacuum-dried and re-suspended in 250 mM acetic acid with 30% acetonitrile for phosphopeptide enrichment with Phos-Select Iron Affinity Gel (Sigma-Aldrich) according to the protocol from Thingholm et al. (2008). Eluted phosphopeptides were desalted and analyzed by nano LC-MS/MS on a TripleTOF 5600 (Sciex, Canada) coupled a NanoLC-2DPlus system with nanoFlex ChiP module (Eksigent, Sciex).
Transcriptome and Translatome Analysis
Transcriptomic and translatomic analyses were performed on two biological replicates using 7-day-old plantlets from the two independent TOR RNAi and GUS control lines grown in vitro and treated with ethanol for 24 h. Transcriptome analyses using CATMA arrays were performed on total RNA preparations as previously described (Moreau et al., 2012). For translatomic analyses total RNA was extracted and polysomal fractions were purified on sucrose gradients after ultracentrifugation as previously described (Deprost et al., 2007; Sormani et al., 2011). Polysome-bound RNAs were extracted using guanidinium hydroxychloride and precipitated by isopropanol and linear acrylamide as a carrier. Subsequently, RNAs were reverse transcribed and hybridized on CATMA arrays as described above for the determination of differentially translated mRNAs (Sormani et al., 2011). Statistical analysis of each comparison was based on two dye swaps and followed by the analysis described by Gagnot et al. (2008) and Moreau et al. (2012). Briefly, an array-by-array normalization was performed to remove systematic biases. To determine differentially expressed genes, we performed a paired t-test on the log ratios averaged on the dye swap. The raw p-values were adjusted by the Bonferroni method, which controls the family-wise error rate to keep a strong control of the false positives in a multiple-comparison context. We considered as being differentially expressed the probes with a Bonferroni P-value ≤0.05, as described by Gagnot et al. (2008). The results are available online in the CatDB database2.
Antibodies directed against phosphorylated RPS6A were obtained by conjugating to keyhole limpet hemocyanin the SRLpSSAAAKPSVTA (phosphoSer240) peptide (produced by Proteogenix, Schiltigheim, France) and injecting two New Zealand White female rabbits (performed by Proteogenix). Seven injections were performed over a period of 56 days. Then a preliminary ELISA test was performed at day 63 to evaluate the titer of the antibodies and the rabbits were bled at day 70. The obtained antisera were first depleted against immobilized non-phosphorylated peptide then specific antibodies were purified using the phosphorylated SRLpSSAAAKPSVTA peptide. The specificity of the purified antibodies was evaluated using an ELISA test with the phosphorylated and non-phosphorylated peptides (produced by Proteogenix) (Supplementary Figure S3).
SDS-PAGE and Western Blotting
Primary antibodies used in this study are directed against mammalian RPS6 (Cell Signaling Technology #2317S), rapeseed RPL13 (Sáez-Vásquez et al., 2000) and Arabidopsis RPS14 (Agrisera AS09 477).
For detection of phosphorylated RPS6, total proteins were extracted from either wild-type (Col-0), TOR RNAi or control Arabidopsis lines using the Laemmli buffer and blotted with our RPS6 phospho-specific antibody. Bradford assay (Bio-Rad) was performed to quantify total protein concentrations. Ten micrograms of proteins were separated by SDS-PAGE gels and transferred to polyvinylidene difluoride membranes (PVDF, Bio-Rad) by electroblotting. Membranes were probed with either antiphospho-RPS6 (P-RPS6) rabbit polyclonal antiserum (dilution 1:5000) or with anti-RPS6 mouse monoclonal IgG (dilution 1:1000). Goat anti-rabbit IgG-HRP (horseradish peroxidase 1:2000, Santa Cruz Biotechnologies) and Goat anti-mouse IgG-HRP (1:2000, Santa Cruz Biotechnologies) were used as secondary antibodies. Immunodetection was performed by using enhanced chemiluminescent (ECL) substrates for HRP as recommended by the manufacturer (Clarity Western ECL blotting substrate Bio-Rad). Transferred proteins on PVDF membranes were visualized by Ponceau S staining.
The 5′ UTR sequences of the ribosomal proteins mRNAs were obtained from the TAIR10 database3. The identification of the motifs was performed with the online MEME software4 (Bailey et al., 2009) with a motif recognition size comprised between 6 and 50 nt. Only the representative isoforms of the genes were used.
Arabidopsis RPS6A: At4g31700 and RPS6B: At5g10360. Proteomic raw data are available in the Protic database under the following accession name: tor_inactivation5.
Ribosome Enrichment by Density Ultracentrifugation
The aim of this study was to identify, by an untargeted proteomic analysis, modifications in the Arabidopsis ribosome fraction in response to TOR inactivation. First, we evaluated the suitability of our ribosome extraction method for LC-MS/MS analysis. To do so, 7-day-old seedlings of Arabidopsis (Col-0) were harvested and, to prevent protease, phosphatase, or kinase activities, immediately submitted to the ribosome purification protocol (see, Materials and Methods). Finally, high molecular weight particles, including polysomes, were pelleted by ultracentrifugation through a sucrose cushion.
We then submitted the ultracentrifugated fraction to SDS-PAGE and resolved proteins were stained using silver nitrate (Figure 1A). The obtained protein profile is typical of purified plant ribosomal fractions (Carroll et al., 2008). The presence of ribosomal proteins in this fraction was confirmed by Western blot analysis which allowed to normalize the protein fractions by diluting the samples according to Western blot quantifications (Figure 1B).
FIGURE 1. Determination of ribosomal protein amounts prior to LC-MS/MS analyses. Plant extracts were submitted to ultracentrifugation through a sucrose cushion to obtain a ribosome-enriched fraction. Pellets were resuspended in Laemmli buffer and analyzed for the abundance of ribosomal proteins by SDS-PAGE and Western blot. (A) Silver nitrate stained gel after SDS-PAGE. (B) Western blot against the RPS6, RPL13, and RPS14 ribosomal proteins. Exp1 and Exp2 correspond to two independent biological replicates. GUS is the control line, RNAi1 and RNAi2 are the two independent TOR RNAi lines. All the lines were induced with ethanol. MM, molecular marker.
To inactivate TOR, we used two independent ethanol-inducible TOR RNAi lines (based on the AlcR/AlcA operon) that we previously obtained and characterized (Deprost et al., 2007) and compared them to a control line expressing an ethanol-inducible GUS gene (GUS control, Deprost et al., 2007; Dobrenel et al., 2011).
Identification of Ribosomal Proteins by LC-MS/MS Analysis
Seedlings of two independent RNAi lines and of the GUS control were grown for 7 days in vitro and then treated with ethanol to induce TOR inactivation. We repeated the same experiment twice and then analyzed the six samples together. In order to remove eventual contaminating rRNA as well as the chemicals, the samples were first submitted to a short migration through a SDS-PAGE stacking gel. Proteomic as well as phosphoproteomic analyses by LC-MS/MS identified a total of 5936 spectra, corresponding to 1508 unique peptide sequences (raw data are available in the Protic database). Peptides were matched to protein sequences from TAIR10 and grouped according to sequence homology using the X!TandemPipeline with the phosphorylation of serine, threonine and tyrosine as a potential peptide modification. By this method, 361 different proteins were potentially identified by at least two peptide sequences, belonging to 217 groups (corresponding presumably to protein families, based on sequence homologies). Among these 361 proteins, 210 were identified by the presence of at least two proteotypic (i.e., specific of a given protein) peptides. Based on the previous annotations of the ribosomal proteins (Sormani et al., 2011; Tiller and Bock, 2014; Hummel et al., 2015), we found that more than half of the 361 proteins were ribosomal proteins (147 correspond to cytosolic ribosomes and 46 to organelle ribosomes) (Figure 2A). Using an in silico analysis, Sormani et al. (2011) refined the annotation of ribosomal proteins (including plastidic and mitochondrial ones). We used this list to identify the proteins (and therefore the peptides) corresponding to the mitochondrial and plastidic ribosomes. All identified organellar peptides corresponded to plastidic ribosomes.
FIGURE 2. Effect of TOR inactivation on the cytosol and organelle ribosomal proteins. (A) Distribution of the proteins identified in at least one of the samples depending on whether they are part of the cytosolic ribosome, the organelle ribosome or non-ribosomal. (B) Venn diagram showing the number of peptides corresponding to the cytosolic ribosomal proteins identified in the GUS control and in the TOR RNAi lines. The intersections of two peptide sets are shown in orange and the intersection of the three sets in red. The total numbers of peptides are shown and the areas are representative of the number of common or specific peptides. (C), same as (B) for the organelle ribosome except that intersections are shown in medium and dark green for two or three peptides sets, respectively. (D) Boxplot representing the relative abundance (log2 fold change) of the peptides derived from ribosomal proteins in the TOR RNAi compared to the GUS control lines depending of their localization in the cytosolic or plastidic ribosomes. For each peptide a mean abundance was obtained from the four repetitions (two independent experiments using the two RNAi lines) and the abundance ratios were compared.
Target of Rapamycin inactivation in the RNAi lines did not largely affect the number of detected peptides originating from the cytosolic ribosomes and thus most peptides were found both in the RNAi and in the control lines (between 419 and 449 peptides for RNAi2 and GUS lines, respectively; Figure 2B). Conversely the pool of peptides coming from pRPs was specifically depleted in the TOR RNAi lines with only 116 and 111 peptides detected for the TOR RNAi1 and RNAi2 lines, respectively, compared to 201 peptides in the control GUS line (Figure 2C). One third of the pRPs peptides were only present in the control GUS line (74 out of 218) whereas a much smaller number of peptides (17) were specifically found in the RNAi lines.
Quantitative Analysis of the Expression of Ribosomal Protein Encoding Genes
In order to exclude biases that could be caused by potential technical issues, like the mass spectrometer being occupied by some abundant peptides eluting near the peptides of interest, we quantified the traces of the peaks corresponding to the identified peptides of the plastidic and cytoplasmic RPs. These peak areas were then compared between the RNAi and GUS control lines for each peptide (Figure 2D). Such a quantitative approach confirmed that most of peptides resulting from the fragmentation of the pRPs are indeed less abundant in the TOR RNAi samples. This also confirmed that TOR inactivation did not have any strong effect on the accumulation of the cRPs peptides (Figure 2D). Thus these results suggest a global decrease in the abundance of pRPs resulting in a lower number and amount of peptides detected in the LC-MS/MS analysis.
We then used the number of spectra (or peptide hits) per protein to estimate the protein abundance for the pRPs and cRPs. Indeed it has been shown that there is a proportional relationship between these two values (Allet et al., 2004; Liu et al., 2004). We exclusively used the spectra corresponding to proteotypic peptides in order to avoid a bias that would be caused by the very high sequence homology within ribosomal protein families. By this method, we showed a coordinated down-regulation of the pRPs while the cRPs have a much less coordinated profile and are globally only slightly affected by the TOR inactivation (Figure 3). To better understand the role of TOR in regulating plant gene expression, we performed transcriptomic and translatomic analyses after TOR inactivation, in which we monitored the mRNA levels on total and polysome-bound RNA samples. The variations in the abundance of polysome-bound mRNA were determined after purification of polysomes on a sucrose gradient, extraction of RNA and microarray hybridization. To better identify the primary TOR-regulated mRNA targets, we decided to shorten the time of ethanol-mediated RNAi induction. Two biological repetitions were performed using each time the two independent RNAi lines. The resulting four transcriptomic and translatomic experiments were submitted to a statistical analysis to identify common differentially expressed genes when compared to the GUS control lines (see Materials and Methods). When focusing on the ribosomal proteins, we found a large difference in the expression profiles of the corresponding nuclear genes depending on whether the gene product is part of the cytoplasmic or the plastidic ribosome (Figure 3; Supplementary Table S1). First, almost all detected nuclear-encoded mRNA coding for pRPs, whether located in the small or in the large subunit, showed either no change in abundance or were down-regulated after TOR inactivation when compared to the GUS control line treated with ethanol (Figure 3A). On the opposite mRNAs coding for cRPs were mostly up-regulated (Figures 3B,C) except for RPL18a-1 which was the only transcript showing a reproducible decrease in abundance. The translatomic analysis mostly mirrored the transcriptomic variations but for some genes coding for cytosolic ribosomes, like RPL40 and RPL41, the total mRNA abundance did not vary significantly whereas they were more engaged in polysomes compared to the GUS control, suggesting a higher translation of these genes (Figure 3C). Taken together, these data suggest that there is a coordinated down regulation of the nuclear genes coding for the pRPs at the transcriptional, translational (polysomal loading), and protein levels in response to TOR inactivation (Figure 3A).
FIGURE 3. Effect of TOR inactivation on the transcription, translation, and abundance of ribosomal proteins. (A) Expression of the nuclear genes coding for chloroplastic ribosomal proteins at the total mRNA level (transcriptomic analysis: Tx), at the polysome-bound mRNA level (translatomic analysis: Tl) and at the protein level (proteomic analysis). (B,C) Expression level of the genes coding for the cytosolic ribosomal proteins and localized in the small or large ribosome subunit, respectively. For the transcriptome and translatome experiments, mRNA abundances in the TOR RNAi lines were compared to the GUS control line treated with the same 24 h ethanol induction time. The results represent the mean of the two TOR RNAi lines in the two independent biological replicates. Intensity ratios are shown in log2 scale according to the color scale shown in the figure on the bottom right for transcriptomic and translatomic experiments (Transcr./Transl.). For the proteomic experiment, two independent biological replicates (X1 and X2) using the two TOR RNAi lines are presented. Orange and pale green represent a quantitative up- and down-regulation, respectively, and red and dark green represent a qualitative regulation, respectively (peptides either present or absent). Only the proteotypic spectra were used. Black boxes indicate no change and white boxes represent missing data.
Next we examined whether this co-regulation involves the recognition of a conserved motif in the 5′ UTR sequences of the nuclear-encoded mRNA coding for pRPs. These sequences were analyzed for enriched motifs using the MEME software (Bailey and Elkan, 1994). A strongly significant motif composed of a stretch of pyrimidines was identified in this set of 5′ UTRs (Figure 4A). This sequence is reminiscent of the animal TOP motif found within the 5′UTR of mammalian ribosomal protein encoding mRNAs, which confers translational regulation by TOR (Meyuhas and Kahan, 2015). The MEME analysis also revealed the presence of a second A/G-enriched motif (Figure 4A). On the contrary, the 5′ UTRs of the mRNAs coding for cRPs were significantly enriched for a TTTAGGGTTT motif (Figure 4B), which is similar to the telo-box consensus sequence already identified in the promoters of these genes (Figure 4B) (McIntosh et al., 2011; Wang et al., 2011). Next we analyzed the 5′ UTRs of the nuclear genes coding for pRPs, which were less engaged in polysomes after TOR inactivation (Figure 3A). MEME analysis identified a shorter pyrimidine-rich motif that is more similar to canonical TOP motifs, but even more to the pyrimidine-rich translational element (PRTE), a motif identified in animal genes which translation is controlled by TOR (Figure 4C) (Hsieh et al., 2012). This motif is found in the majority of animal TOR targets and, unlike conventional TOP motifs, does not reside at the start of the mRNA sequence (Hsieh et al., 2012). Consistently, the motif we identified is also rarely present at the start of the analyzed mRNAs (Figure 4D; Supplementary Figure S1) but is significantly enriched in pRP transcripts translation of which is affected by TOR inactivation. A purine-rich motif was also identified in this subset of genes and was often found 3′ to the PRTE-like motif (Figure 4D).
FIGURE 4. Motifs identified in the 5′ UTRs of genes coding for ribosomal proteins. (A) Motifs identified as being significantly enriched in the 5′ UTRs of the nuclear genes coding for the plastidic RPs after analysis by the MEME software. (B) Same for genes coding for the cytosolic RPs. The telo-box motif found in cRPs promoters is shown. (C) Same as for (A) but only the 5′ UTRs of the genes showing a down-regulation in the translatome experiments following TOR inactivation were kept. (D) Positions of the motifs identified in (C) in the 5′ UTRs of the analyzed sequences. See Supplementary Figure S1 for details.
We then mined the public Genevestigator transcriptome database (Hruz et al., 2008) for information about the expression profile of plastidic ribosomal genes encoded by the nuclear genome. Interestingly these genes were found to be down-regulated in estradiol-inducible TOR RNAi lines (Xiong et al., 2013; Supplementary Figure S2) which suggests that TOR inactivation reproducibly reduces their expression. Nuclear genes coding for pRPs are also strongly induced during germination or after light treatment, which is consistent with their role in the formation of the photosynthetic machinery. Conversely they were repressed in response to various stresses such as pathogen infection, drought, hypoxia, or increased temperature as well as in response to extended night. Finally these genes were slightly induced in one experiment of sucrose feeding but globally repressed in several experiments of nitrogen starvation (Supplementary Figure S2).
TOR Dependent Phosphorylation of RPS6 on Ser240
Even if most of the cRPs were not affected by TOR inactivation (Figure 2D) we found 30 peptides which were significantly down regulated (t-test p-value < 0.05; Supplementary Table S2). In most of the case, these peptides are not proteotypic, thus making the conclusions more complicated. The two cytosolic RPS6 paralogs were identified with three different proteotypic peptides that are significantly down-regulated in a quantitative manner following TOR inactivation (Figure 3B). In total, 15 peptides were detected for this protein family. Among these 15 peptides, we identified a phosphorylated peptide corresponding to the RPS6B C-terminal extremity. The X!Tandem analysis predicted a phosphorylation site on Ser240 (SRLpSSAPAKPVAA: Figure 5; Supplementary Figure S3). This serine seems to be conserved in the eukaryotic RPS6 proteins, including plants, and was previously identified as being phosphorylated following TOR activation in yeast and animals (Pende, 2006; Meyuhas, 2008; Yerlikaya et al., 2016) (Figure 6). In order to clarify the actual position of the phosphorylated residue in this peptide, the phosphoproteomic results were submitted to a Phoscalc statistical analysis (MacLean et al., 2008). This analysis failed to clearly identify the phosphorylated site between Ser237, 240 or 241 (Supplementary Figures S3C,D). Nevertheless it showed that this RPS6B C-terminal peptide carries only one phosphate group. An independent phosphoproteomic analysis of Arabidopsis seedlings treated with auxin was also performed (see Materials and Methods). Indeed previous studies suggested that auxin mediates TOR, and thus S6K1, activation in Arabidopsis as shown by phosphorylation of TOR at Ser2424 and S6K1 at Thr449 (Schepetilnikov et al., 2013). Since S6K1 is directly involved in the phosphorylation of RPS6, we asked whether this ribosomal protein is phosphorylated at Ser240 in Arabidopsis extracts treated by auxin. After purification of the phosphopeptides by immobilized metal affinity chromatography (IMAC), Ser240 was again identified as a potential phosphorylation site both in RPS6A and B (Supplementary Table S3; Supplementary Figures S3E,F). In some cases RPS6A Ser240 was identified together with Ser237 (pSRLpSSAPAKPVAA) but again the discrimination between phosphorylated Ser240 and Ser241 was not always possible (Supplementary Table S3).
FIGURE 5. Downregulation of RPS6A and RPS6B phosphorylation in response to TOR inactivation. (A) Phosphorylation level of the RPS6B C-terminal peptide after 72 h of ethanol treatment. (B) Same for the RPS6A and (C) for the RPS6B proteins after 10 days of ethanol treatment. Each bar corresponds to the peak intensity of the phosphorylated peptide normalized by a proteotypic peptide from the same protein (GENDLPGLTDTEKPR for RPS6A and GVSDLPGLTDTEKPR for RPS6B). For each experiment, the results of the GUS control (GUS), the TOR RNAi1 (R1) and the TOR RNAi2 (R2) lines are presented. Exp1 and Exp2 correspond to two independent biological replicates with 72 h or 10 days of TOR inactivation.
FIGURE 6. Comparison of RPS6 amino acid sequences. Arabidopsis RPS6s (AtRPS6A, At4g31700.1, and AtRPS6B, At5g10360.1), Nicotiana tabacum RPS6 (NtRPS6, P29345.2), Zea mays RPS6s (ZmRPS6-1, NP_001105544.1 and ZmRPS6-2, NP_001105634.1), Homo sapiens RPS6 (HsRPS6, NP_001001.2), Mus musculus RPS6 (MmRPS6, NP_033122.1), Drosophila melanogaster RPS6s (DmRPS6A, NP_727213.1 and DmRPS6B, NP_511073.1) and Saccharomyces cerevisiae RPS6 (ScRPS6, NP_015235.1) protein sequences were aligned by T-Coffee (http://tcoffee.crg.cat). The level of conservation is represented by a color code in which the most conserved residues are in black, those with intermediate conservation are gray and the least conserved are in white. RPS6 peptides identified in Arabidopsis by the proteomic analysis are represented by lines over the corresponding amino acid sequences (black for the peptides common to the two isoforms and gray for peptides specific to one isoform). An asterisk shows the conserved phosphorylated Ser240 residue.
Even if the precise localization of the phosphorylated serine residue in the RPS6B protein C-terminus could not be determined without ambiguities, we compared the abundance of the C-terminal monophospho-SRLSSAPAKPVAA peptide between the TOR RNAi lines and the GUS controls. After 72 h of exposure to ethanol we observed a modest decrease in the RPS6B C-terminal phosphorylated peptide (Figure 5A). This decrease in phosphorylation was only observed in one of two experiments. Since the abundance of the phosphorylated peptide was already low in the control line for the second experiment, it could be that the decrease in phosphorylation level was less obvious in this case. The corresponding RPS6A peptide was not found in this analysis. Thus to confirm this TOR-dependent decrease in RPS6 C-terminus phosphorylation, we performed a longer silencing induction by ethanol for 10 days. As for the 72 h treatment, biological duplicates have been used. The abundance of the monophospho C-terminal peptides was compared for both RPS6A and RPS6B proteins between the RNAi and the control GUS lines (Figures 5B,C). The amount of phosphorylated C-terminal peptide was lower (with a mean of 40% decrease) for both the RPS6A and the RPS6B proteins in the TOR RNAi lines when treated with ethanol for 10 days.
To confirm this phosphoproteomic analysis, a phospho-specific antibody was raised against a synthetic SRLpSSAAAKPSVTA peptide in which only Ser240 was phosphorylated. The obtained polyclonal antibody was purified against this peptide and the eluted antibody fraction detected a single band corresponding in size to RPS6 proteins (Supplementary Figure S4A; Figure 7). Furthermore this antibody was found to be highly specific for the phosphorylated SRLpSSAAAKPSVTA peptide, showing no reaction with the control non-phosphorylated peptide in an ELISA test (Supplementary Figure S4B).
FIGURE 7. Detection of RPS6 phosphorylation by Western blot assay. Total protein extracts obtained from seedlings were separated by SDS-PAGE and blotted onto a membrane. After incubation with the phospho-specific antibody against RPS6 Ser240 (P-RPS6) or a monoclonal antibody against total mammalian RPS6, blots were revealed by a secondary antibody linked to HRP activity and imaged with a CCD camera (see Materials and Methods for details and Supplementary Figure S4). (A) AZD treatments inhibit Ser240 RPS6 phosphorylation. Six day-old seedlings were either mock or AZD treated for 24 h. NT, non-treated plants at time 0 and 24 h. (B) Sucrose treatment induces RPS6 phosphorylation. Six day-old seedlings were transferred to sugar-free medium for 24 h and then either mock (0–24 h) or sucrose (0,5%) treated for 24 h. 0: plants before sucrose induction at time 0. (C) Silencing of TOR decreases RPS6 phosphorylation. The control (GUS) or TOR RNAi lines were grown for 7 days and induced with 5% (v/v) ethanol (EtOH). Bottom panels show Ponceau Red staining of the membranes.
This antibody was used in an optimized Western blot assay and the obtained signal was very strong even when the antibody was diluted 1/5000. This band co-migrated with the band decorated by a specific monoclonal antibody against mammalian RPS6 (Figure 7) which was subsequently used to quantify the amount of total RPS6 in protein extracts. Both treating Arabidopsis seedlings with AZD-8055, a strong and specific second generation TOR inhibitor (Montané and Menand, 2013), or silencing the expression of TOR by a 7-day ethanol treatment, resulted in a significant and dose-dependent decrease in the signal obtained with the RPS6 phospho-specific antibody, whereas the total amount of RPS6 protein only decreased slightly (Figures 7A,C). Since this antibody is highly specific for the phosphorylated RPS6 protein, this confirms that there is a reproducible decrease in Ser240 phosphorylation level following TOR inactivation. Accordingly it was previously shown in animals that AZD inhibits TOR, and as a consequence S6K activity and RPS6 phosphorylation (Chresta et al., 2010). TOR and S6K were previously shown to be activated by sugars like sucrose and glucose (Xiong et al., 2013). Consistently RPS6 phosphorylation was augmented by the addition of sucrose when supplied to sugar-starved Arabidopsis seedlings (Figure 7B). Next we examined the kinetic of changes in RPS6 phosphorylation after either AZD-8055 or sucrose treatments (Figure 8). As soon as 1 h after AZD addition, a significant decrease in Ser240 phosphorylation was observed (Figure 8A). For sucrose, an increase in phosphorylation was detected 2 h after supply. However it is difficult at this stage to discriminate between a direct signaling effect of sucrose and an indirect consequence of sugar metabolism (Figure 8B). Altogether these data suggest a strong positive correlation between the RPS6 phosphorylation level, as detected by this specific polyclonal antibody, and TOR activity. Therefore, it seems that in plants, like in animals and yeast, RPS6 phosphorylation can be used as a robust and sensitive readout for TOR activity.
FIGURE 8. Kinetics of variations in RPS6 phosphorylation following AZD-8055 or sucrose addition. Total protein extracts obtained from seedlings were separated by SDS-PAGE and blotted onto a membrane. After incubation with the phospho-specific antibody against RPS6 Ser240 (P-RPS6) blots were revealed by a secondary antibody linked to HRP activity and imaged with a CCD camera. (A) Kinetics of the inhibition of RPS6 Ser240 phosphorylation by AZD 8055. Seven day-old seedlings were either mock (DMSO) or AZD-8055 (2 μM) treated and then harvested at the indicated time. (B) Kinetics of the induction of RPS6 Ser240 phosphorylation by sucrose. Seven day-old seedlings were transferred in sugar-free medium for 24 h and then either mock (0) or sucrose (0.5%) treated. Seedlings were harvested at the indicated time. Bottom panels show Ponceau Red staining of the membranes.
Several studies have investigated the impact of TOR inhibition in Arabidopsis on transcript or metabolite levels (Deprost et al., 2007; Moreau et al., 2012; Caldana et al., 2013; Xiong et al., 2013; Dong et al., 2015) but hitherto the global proteome has not been examined. In this paper, we investigated the expression of ribosomal proteins and genes using transcriptome, translatome, proteome and phosphoproteome analyses following silencing of the TOR gene. Concerning the cytoplasmic ribosome, we identified in our proteomic experiments 65 families of ribosomal proteins, corresponding to 69 ribosomal protein isoforms identified with at least two proteotypic peptides. Hummel et al. (2015) identified 165 cRPs isoforms by LC-MS/MS after a tryptic digestion. We were able to find more than one third of these protein paralogs and only 16 families were not found in this analysis. Seven of these 16 families cannot be identified by LC-MS after a tryptic digestion mainly because of the small size of the resulting peptides due to their high content in lysine and arginine. Thus, only nine out of the 81 ribosomal protein families were missing in our analysis (RPL22, RPL35a, RPL38, RPP3, RPS27, RPS28, RPS29, RPS30, and RACK1). Moreover, this work has been focusing solely on young seedlings while Hummel et al. (2015) were also using rosettes in their analysis and some specific paralogs may be expressed only in specific tissues or at some specific developmental stages (Weijers et al., 2001; Sormani et al., 2011). Since we performed our proteomic analyses using plants silenced for TOR expression by a long ethanol treatment (3 and 10 days), our analysis of changes in protein abundance may reveal steady-state long-term, and sometimes indirect, adaptation to a decrease in TOR activity.
In plastids translation occurs on 70S bacterial-type ribosomes. About half of the plastid small ribosome subunit (30S) proteins and most of the large subunit (50S) proteins are encoded by the nuclear genome (Yamaguchi and Subramanian, 2000). TOR silencing (Deprost et al., 2007; Xiong and Sheen, 2012; Caldana et al., 2013), inhibition by rapamycin (Sormani et al., 2007; Ren et al., 2011) or by AZD-8055 (Montané and Menand, 2013; Li et al., 2015) as well as mutations affecting the TORC1 complex (Moreau et al., 2012; Kravchenko et al., 2015) consistently result in leaf chlorosis and yellowing. We show here that TOR inhibition results in a coordinated decrease in pRP expression, at the level of protein abundance but also at the total and translated mRNA levels, which could explain these chlorotic phenotypes (Figure 3). Whether this is the result of a decreased synthesis or an increased degradation of the chloroplast components by autophagy, which is induced after TOR inactivation (Liu and Bassham, 2010), remains to be determined. Interestingly the expression of nuclear genes coding for cytosolic proteins was found to be mostly induced whereas the level of proteins often decreased. The same trends were observed in N-limited Chlamydomonas where the levels of pRPs as well as their corresponding mRNAs decreased in response to N starvation whereas only protein levels decreased for cRPs (Schmollinger et al., 2014). The same effect of TOR inactivation on the expression profile of the pRPs was confirmed by the comparison with transcriptomic data obtained using estradiol-inducible TOR RNAi line (Supplementary Figure S2) (Xiong et al., 2013). Interestingly, nuclear genes coding for pRPs were down-regulated in response to several abiotic or biotic stresses suggesting that they may play an important role in the adaptation to stresses (Supplementary Figure S2). These genes were also repressed by ABA but strongly induced by the application of brassinolide. This is consistent with the inhibitory effects of ABA on the growth-promoting hormones like brassinosteroid and with the role of TOR in brassinosteroid (Zhang et al., 2009, 2016) and ABA signaling (Kravchenko et al., 2015; Li et al., 2015), which could result in TOR inhibition. As reported earlier by Pal et al. (2013) we found diverse and uncoordinated responses to variations in sugar supply for nuclear genes coding for pRPS while their expression was repressed by nitrogen starvation as observed in Chlamydomonas (Schmollinger et al., 2014). Altogether, these data suggest that nuclear genes coding for the pRPs are controlled by TOR at multiple levels to integrate environmental cues for the regulation of chloroplastic translation.
In animals TOR is known to regulate the translation of TOP-containing mRNAs (Hsieh et al., 2012; Thoreen et al., 2012; Meyuhas and Kahan, 2015). This motif is particularly present in the 5′ UTR of genes coding for ribosomal proteins or components of the translation machinery (Meyuhas and Kahan, 2015). Using a MEME analysis we did not detect any specific enrichment for TOP motif in cRPs. The most abundant motif was related to the telo-box (Figure 4A). This DNA motif is found in the 5′ regions, often close to the start codon, of nuclear genes coding for both mitochondrial and cytosolic RPs, but not in the genes encoding pRPs (Wang et al., 2011). Conversely, we found a highly significant occurrence of a pyrimidine-rich motif in the 5′ UTR of nuclear genes coding for pRPs. This motif is reminiscent of TOP motifs (Figure 4A). A TOP-like motif was also previously identified in mRNA coding for ribosomal proteins in maize embryonic axes (Jiménez-López et al., 2011) but these motifs are so far poorly described in plants. Canonical TOP motifs are located at the start of the mRNA, which is not the case in our analysis. Nevertheless, a previous study has demonstrated the presence of several transcription start sites in genes coding for pRPs (Lagrange et al., 1993). Interestingly for the plastidic RPL21 gene, the start site specifically used in leaves produced a mRNA starting with a canonical TOP motif. The 5′ UTRs of cRP genes which are less translated after TOR inactivation were enriched in a motif that is strikingly similar to the PRTE motif found within the 5′ UTRs of animal genes controlled by TOR at the level of translation (Figures 4C,D) (Hsieh et al., 2012). It is thus tempting to hypothesize that TOR has been recruited in plants to regulate specifically in leaves the translation of nuclear mRNAs coding for chloroplastic ribosomal proteins. It was previously shown that translation of animal TOP-containing mRNAs can be differentially regulated in vitro in a wheat germ extract (Shama and Meyuhas, 1996) and that auxin stimulates S6 ribosomal protein phosphorylation on maize ribosomes and the recruitment of TOP-like mRNAs for translation (Beltrán-Peña et al., 2002). Since TOR is also activated by auxin (Schepetilnikov et al., 2013) these data suggest that plant mRNAs containing TOP-related motifs could also be regulated in a TOR- and phosphorylated S6-dependent manner.
Only phosphorylation of the RPS6 protein could be identified in a reproducible manner in previous unbiased phosphoproteomic analyses of the ribosomes performed in eukaryotes (Huber et al., 2009; Hsu et al., 2011; Yu et al., 2011). Nevertheless, despite a wealth of studies and several hypotheses, the precise biological role of these conserved phosphorylation events remains disputed and unclear. For example, expression of human RPS6 containing alanine at all phosphorylated serine residues did not modify the overall translation rate, even for TOP mRNAs (Ruvinsky et al., 2005). Instead cell growth and size as well as ribosome biogenesis were affected (Ruvinsky et al., 2005; Chauvin et al., 2014). The same result was observed in yeast expressing non-phoshorylatable RPS6 (Yerlikaya et al., 2016). However, phosphorylation of the C-terminal Ser residues of RPS6 has been used as a robust and recognized readout for TOR activity in animals and yeast (Meyuhas, 2008, 2015; Yerlikaya et al., 2016). In this study, a decrease in RPS6 phosphorylation in response to TOR inactivation was observed (Figures 5 and 7). The phosphoproteomic analysis identified a C-terminal phosphorylation site in both the RPS6A and RPS6B proteins that is TOR activity-dependent without unambiguously determining which of the C-terminal serine residues is modified. In Arabidopsis Ser237 was previously identified by MALDI-TOF as being phosphorylated but the absence of fragmentation in the C-terminal region hindered the precise localization of the other modification sites by MS/MS analysis (Chang et al., 2005). Several phosphoproteomic studies of the plant ribosome have already shown the presence of a phosphorylation site in the C-terminal peptide of the RPS6 and have suggested that the Ser240 could be one of the modified residues together with Ser229, 231, or 237 (Carroll et al., 2008; Turkina et al., 2011; Boex-Fontvieille et al., 2013). A global phosphoproteome analysis of Arabidopsis identified Ser237 and 240 as being phosphorylated together with Ser247 and Thr249 (Reiland et al., 2009). Ser240 is conserved in all plant RPS6 sequences whereas Ser241 is missing in the maize and tobacco sequences (Figure 6). Conversely Ser237 is found in all plant sequences but only Ser240 can be aligned with one of the known phosphorylated serine in the yeast (Ser232) or human (Ser235) RPS6 sequence (Figure 6; Meyuhas, 2015). It is well known in yeast and animals that TOR activity controls RPS6 phosphorylation through activation of S6K (Wullschleger et al., 2006; Biever et al., 2015; Meyuhas, 2015) and Mahfouz et al. (2006) have established that TOR interacts with S6K through RAPTOR to activate RPS6 phosphorylation. Nevertheless it should be noted that S6K is also activated by the 3-phosphoinositide-dependent protein kinase 1 (PDK1) which operated after S6K phosphorylation by TOR (Mahfouz et al., 2006; Otterhag et al., 2006). Interestingly the SnRK1 kinase, which probably acts antagonistically to TOR (Dobrenel et al., 2016), was recently shown to interact with and phosphorylate RAPTOR in Arabidopsis (Nukarinen et al., 2016). Moreover, a strong increase in RPS6 Ser240 phosphorylation was also observed after SnRK1 inactivation. This is coherent with the hypothesis that SnRK1 inhibits TOR activity, and hence RPS6 phosphorylation, presumably through RAPTOR phosphorylation (Nukarinen et al., 2016).
Western blot assays using the RPS6 Ser240 phospho-specific antibody demonstrated that phosphorylation of this residue decreased following TOR inactivation either by silencing or by using a specific inhibitor (Figure 7). Therefore this assay could be used as a TOR readout in plants. Previous assays for TOR activity in plants were based on the detection of Thr449 phosphorylation in S6K by commercial antibodies directed against phosphorylated Thr389 in animal S6K. However, these antibodies produce many non-specific bands in Western blot assays (Schepetilnikov et al., 2011; Xiong and Sheen, 2012). Moreover the abundance of plant S6K is low in plants whereas RPS6 is present in large amounts.
Taken together these data show that the TOR-dependent C-terminal RPS6 phosphorylation is conserved in plants like in other eukaryotes. We have taken advantage of this conserved phosphorylation to design a sensitive and specific assay to monitor TOR activity in plants. The question that remains open is the biological role of RPS6 phosphorylation. Structural studies have shown that RPS6 is accessible to the solvent, and hence to kinases, but the disordered C-terminal region is unfortunately absent from the resolved ribosome structure (Khatter et al., 2015). Nevertheless the charge modifications produced by phosphorylation of RPS6 probably have important biological roles either within the ribosome or for extra-ribosomal functions of RPS6. Indeed it was recently reported that RPS6 affects ribosomal RNA production and interacts with the HD2B histone deacetylase (Kim et al., 2014). More work is therefore needed to elucidate the role of TOR in regulating translation or development through the conserved RPS6 phosphorylation.
TD, MS, MZ, CR, LR, JH, and CM conceived the research plan and supervised the experiments, TD, EM-M, CF, MA, MD, MM, JC, and OL performed the experiments and analyzed the data, TD, CF, JH, CR, and CM wrote the article with contributions of all the authors.
This work was supported by ANR grants (ANR14-CE19-007 and ANR11- SV6-01002) to CF, CR, LR, MS, and CM. MM was supported by a joint Ph.D. grant from INRA (Plant Biology Department) and DSV CEA.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Annemarie Matthes for fruitful discussions and Rodnay Sormani for his help in purifying ribosomes and for his advises. We particularly wish to thank Proteogenix (France) for excellent service and for producing peptides and the antibodies against RPS6.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01611/full#supplementary-material
- ^ http://pappso.inra.fr/bioinfo/xtandempipeline/
- ^ http://urgv.evry.inra.fr/cgi-bin/projects/CATdb/consult_expce.pl?experiment_id=302
- ^ http://www.arabidopsis.org/
- ^ http://meme-suite.org/
- ^ moulon.inra.fr/protic/tor_inactivation
Allet, N., Barrillat, N., Baussant, T., Boiteau, C., Botti, P., Bougueleret, L., et al. (2004). In vitro and in silico processes to identify differentially expressed proteins. Proteomics 4, 2333–2351. doi: 10.1002/pmic.200300840
Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208. doi: 10.1093/nar/gkp335
Barbet, N. C., Schneider, U., Helliwell, S. B., Stansfield, I., Tuite, M. F., and Hall, M. N. (1996). TOR controls translation initiation and early G1 progression in yeast. Mol. Biol. Cell 7, 25–42. doi: 10.1091/mbc.7.1.25
Beltrán-Peña, E., Aguilar, R., Ortíz-López, A., Dinkova, T. D., and De Jiménez, E. S. (2002). Auxin stimulates S6 ribosomal protein phosphorylation in maize thereby affecting protein synthesis regulation. Physiol. Plant. 115, 291–297. doi: 10.1034/j.1399-3054.2002.1150216.x
Biever, A., Puighermanal, E., Nishi, A., David, A., Panciatici, C., Longueville, S., et al. (2015). PKA-dependent phosphorylation of ribosomal protein S6 does not correlate with translation efficiency in striatonigral and striatopallidal medium-sized spiny neurons. J. Neurosci. 35, 4113–4130. doi: 10.1523/JNEUROSCI.3288-14.2015
Boex-Fontvieille, E., Daventure, M., Jossier, M., Zivy, M., Hodges, M., and Tcherkez, G. (2013). Photosynthetic control of Arabidopsis leaf cytoplasmic translation initiation by protein phosphorylation. PLoS ONE 8:e70692. doi: 10.1371/journal.pone.0070692
Caldana, C., Li, Y., Leisse, A., Zhang, Y., Bartholomaeus, L., Fernie, A. R., et al. (2013). Systemic analysis of inducible target of rapamycin mutants reveal a general metabolic switch controlling growth in Arabidopsis thaliana. Plant J. 73, 897–909.
Carroll, A. J., Heazlewood, J. L., Ito, J., and Millar, A. H. (2008). Analysis of the Arabidopsis cytosolic ribosome proteome provides detailed insights into its components and their post-translational modification. Mol. Cell. Proteom. 7, 347–369. doi: 10.1074/mcp.M700052-MCP200
Chang, I. F., Szick-Miranda, K., Pan, S., and Bailey-Serres, J. (2005). Proteomic characterization of evolutionarily conserved and variable proteins of Arabidopsis cytosolic ribosomes. Plant Physiol. 137, 848–862. doi: 10.1104/pp.104.053637
Chauvin, C., Koka, V., Nouschi, A., Mieulet, V., Hoareau-Aveilla, C., Dreazen, A., et al. (2014). Ribosomal protein S6 kinase activity controls the ribosome biogenesis transcriptional program. Oncogene 33, 474–483. doi: 10.1038/onc.2012.606
Chresta, C. M., Davies, B. R., Hickson, I., Harding, T., Cosulich, S., Critchlow, S. E., et al. (2010). AZD8055 is a potent, selective, and orally bioavailable ATP-competitive mammalian target of rapamycin kinase inhibitor with in vitro and in vivo antitumor activity. Cancer Res. 70, 288–298. doi: 10.1158/0008-5472.CAN-09-1751
Deprost, D., Yao, L., Sormani, R., Moreau, M., Leterreux, G., Nicolaï, M., et al. (2007). The Arabidopsis TOR kinase links plant growth, yield, stress resistance and mRNA translation. EMBO Rep. 8, 864–870. doi: 10.1038/sj.embor.7401043
Dobrenel, T., Caldana, C., Hanson, J., Robaglia, C., Vincentz, M., Veit, B., et al. (2016). TOR signaling and nutrient sensing. Annu. Rev. Plant Biol. 67, 261–285. doi: 10.1146/annurev-arplant-043014-114648
Dobrenel, T., Marchive, C., Sormani, R., Moreau, M., Mozzo, M., Montané, M. H., et al. (2011). Regulation of plant growth and metabolism by the TOR kinase. Biochem. Soc. Trans. 39, 477–481. doi: 10.1042/BST0390477
Dong, P., Xiong, F., Que, Y., Wang, K., Yu, L., Li, Z., et al. (2015). Expression profiling and functional analysis reveals that TOR is a key player in regulating photosynthesis and phytohormone signaling pathways in Arabidopsis. Front. Plant Sci. 6:677. doi: 10.3389/fpls.2015.00677
Gagnot, S., Tamby, J. P., Martin-Magniette, M. L., Bitton, F., Taconnat, L., Balzergue, S., et al. (2008). CATdb: a public access to Arabidopsis transcriptome data from the URGV-CATMA platform. Nucleic Acids Res. 36, D986–D990. doi: 10.1093/nar/gkm757
Giavalisco, P., Wilson, D., Kreitler, T., Lehrach, H., Klose, J., Gobom, J., et al. (2005). High heterogeneity within the ribosomal proteins of the Arabidopsis thaliana 80S ribosome. Plant Mol. Biol. 57, 577–591. doi: 10.1007/s11103-005-0699-3
Holz, M. K., Ballif, B. A., Gygi, S. P., and Blenis, J. (2005). mTOR and S6K1 mediate assembly of the translation preinitiation complex through dynamic protein interchange and ordered phosphorylation events. Cell 123, 569–580. doi: 10.1016/j.cell.2005.10.024
Hruz, T., Laule, O., Szabo, G., Wessendorp, F., Bleuler, S., Oertle, L., et al. (2008). Genevestigator v3: a reference expression database for the meta-analysis of transcriptomes. Adv. Bioinformatics 2008:420747. doi: 10.1155/2008/420747
Hsieh, A. C., Liu, Y., Edlind, M. P., Ingolia, N. T., Janes, M. R., Sher, A., et al. (2012). The translational landscape of mTOR signalling steers cancer initiation and metastasis. Nature 485, 55–61. doi: 10.1038/nature10912
Hsu, P. P., Kang, S. A., Rameseder, J., Zhang, Y., Ottina, K. A., Lim, D., et al. (2011). The mTOR-regulated phosphoproteome reveals a mechanism of mTORC1-mediated inhibition of growth factor signaling. Science 332, 1317–1322. doi: 10.1126/science.1199498
Huber, A., Bodenmiller, B., Uotila, A., Stahl, M., Wanka, S., Gerrits, B., et al. (2009). Characterization of the rapamycin-sensitive phosphoproteome reveals that Sch9 is a central coordinator of protein synthesis. Genes Dev. 23, 1929–1943. doi: 10.1101/gad.532109
Hummel, M., Dobrenel, T., Cordewener, J. J., Davanture, M., Meyer, C., Smeekens, S. J., et al. (2015). Proteomic LC-MS analysis of Arabidopsis cytosolic ribosomes: identification of ribosomal protein paralogs and re-annotation of the ribosomal protein genes. J. Proteom. 128, 436–449. doi: 10.1016/j.jprot.2015.07.004
Jiménez-López, S., Mancera-Martínez, E., Donayre-Torres, A., Rangel, C., Uribe, L., March, S., et al. (2011). Expression profile of maize (Zea mays L.). embryonic axes during germination: translational regulation of ribosomal protein mRNAs. Plant Cell Physiol. 52, 1719–1733. doi: 10.1093/pcp/pcr114
Kim, Y. K., Kim, S., Shin, Y. J., Hur, Y. S., Kim, W. Y., Lee, M. S., et al. (2014). Ribosomal protein S6, a target of rapamycin, is involved in the regulation of rRNA genes by possible epigenetic changes in Arabidopsis. J. Biol. Chem. 289, 3901–3912. doi: 10.1074/jbc.M113.515015
Kravchenko, A., Citerne, S., Jéhanno, I., Bersimbaev, R. I., Veit, B., Meyer, C., et al. (2015). Mutations in the Arabidopsis Lst8 and raptor genes encoding partners of the TOR complex, or inhibition of TOR activity decrease abscisic acid (ABA). synthesis. Biochem. Biophys. Res. Commun. 467, 992–997. doi: 10.1016/j.bbrc.2015.10.028
Lagrange, T., Franzetti, B., Axelos, M., Mache, R., and Lerbs-Mache, S. (1993). Structure and expression of the nuclear gene coding for the chloroplast ribosomal protein L21: developmental regulation of a housekeeping gene by alternative promoters. Mol. Cell. Biol. 13, 2614–2622. doi: 10.1128/MCB.13.4.2614
Li, L., Song, Y., Wang, K., Dong, P., Zhang, X., Li, F., et al. (2015). TOR-inhibitor insensitive-1 (TRIN1). regulates cotyledons greening in Arabidopsis. Front. Plant Sci. 6:861. doi: 10.3389/fpls.2015.00861
MacLean, D., Burrell, M. A., Studholme, D. J., and Jones, A. M. (2008). PhosCalc: a tool for evaluating the sites of peptide phosphorylation from mass spectrometer data. BMC Res. Notes 1:30. doi: 10.1186/1756-0500-1-30
Mahfouz, M. M., Kim, S., Delauney, A. J., and Verma, D. P. (2006). Arabidopsis TARGET OF RAPAMYCIN interacts with RAPTOR, which regulates the activity of S6 kinase in response to osmotic stress signals. Plant Cell 18, 477–490. doi: 10.1105/tpc.105.035931
McIntosh, K. B., Degenhardt, R. F., and Bonham-Smith, P. C. (2011). Sequence context for transcription and translation of the Arabidopsis RPL23aA and RPL23aB paralogs. Genome 54, 738–751. doi: 10.1139/G11-029
Montané, M. H., and Menand, B. (2013). ATP-competitive mTOR kinase inhibitors delay plant growth by triggering early differentiation of meristematic cells but no developmental patterning change. J. Exp. Bot. 64, 4361–4374. doi: 10.1093/jxb/ert242
Moreau, M., Azzopardi, M., Clément, G., Dobrenel, T., Marchive, C., Renne, C., et al. (2012). Mutations in the Arabidopsis homolog of LST8/GβL, a partner of the target of rapamycin kinase, impair plant growth, flowering, and metabolic adaptation to long days. Plant Cell 24, 463–481. doi: 10.1105/tpc.111.091306
Nukarinen, E., Nägele, T., Pedrotti, L., Wurzinger, B., Mair, A., Landgraf, R., et al. (2016). Quantitative phosphoproteomics reveals the role of the AMPK plant ortholog SnRK1 as a metabolic master regulator under energy deprivation. Sci. Rep. 6, 31697. doi: 10.1038/srep31697
Otterhag, L., Gustavsson, N., Alsterfjord, M., Pical, C., Lehrach, H., Gobom, J., et al. (2006). Arabidopsis PDK1: identification of sites important for activity and downstream phosphorylation of S6 kinase. Biochimie 88, 11–21. doi: 10.1016/j.biochi.2005.07.005
Pal, S. K., Liput, M., Piques, M., Ishihara, H., Obata, T., Martins, M. C., et al. (2013). Diurnal changes of polysome loading track sucrose content in the rosette of wild-type arabidopsis and the starchless pgm mutant. Plant Physiol. 162, 1246–1265. doi: 10.1104/pp.112.212258
Reiland, S., Messerli, G., Baerenfaller, K., Gerrits, B., Endler, A., Grossmann, J., et al. (2009). Large-scale Arabidopsis phosphoproteome profiling reveals novel chloroplast kinase substrates and phosphorylation networks. Plant Physiol. 150, 889–903. doi: 10.1104/pp.109.138677
Ren, M., Qiu, S., Venglat, P., Xiang, D., Feng, L., Selvaraj, G., et al. (2011). Target of rapamycin regulates development and ribosomal RNA expression through kinase domain in Arabidopsis. Plant Physiol. 155, 1367–1382. doi: 10.1104/pp.110.169045
Ruvinsky, I., Sharon, N., Lerer, T., Cohen, H., Stolovich-Rain, M., Nir, T., et al. (2005). Ribosomal protein S6 phosphorylation is a determinant of cell size and glucose homeostasis. Genes Dev. 19, 2199–2211. doi: 10.1101/gad.351605
Sáez-Vásquez, J., Gallois, P., and Delseny, M. (2000). Accumulation and nuclear targeting of BnC24, a Brassica napus ribosomal protein corresponding to a mRNA accumulating in response to cold treatment. Plant Sci. 156, 35–46. doi: 10.1016/S0168-9452(00)00229-6
Schepetilnikov, M., Dimitrova, M., Mancera-Martínez, E., Geldreich, A., Keller, M., and Ryabova, L. A. (2013). TOR and S6K1 promote translation reinitiation of uORF-containing mRNAs via phosphorylation of eIF3h. EMBO J. 32, 1087–1102. doi: 10.1038/emboj.2013.61
Schepetilnikov, M., Kobayashi, K., Geldreich, A., Caranta, C., Robaglia, C., Keller, M., et al. (2011). Viral factor TAV recruits TOR/S6K1 signalling to activate reinitiation after long ORF translation. EMBO J. 30, 1343–1356. doi: 10.1038/emboj.2011.39
Schmollinger, S., Mühlhaus, T., Boyle, N. R., Blaby, I. K., Casero, D., Mettler, T., et al. (2014). Nitrogen-sparing mechanisms in chlamydomonas affect the transcriptome, the proteome, and photosynthetic metabolism. Plant Cell 26, 1410–1435. doi: 10.1105/tpc.113.122523
Shama, S., and Meyuhas, O. (1996). The translational cis-regulatory element of mammalian ribosomal protein mRNAs is recognized by the plant translational apparatus. Eur. J. Biochem. 236, 383–388. doi: 10.1111/j.1432-1033.1996.00383.x
Sormani, R., Masclaux-Daubresse, C., Daniel-Vedele, F., Daniele-Vedele, F., and Chardon, F. (2011). Transcriptional regulation of ribosome components are determined by stress according to cellular compartments in Arabidopsis thaliana. PLoS ONE 6:e28070. doi: 10.1371/journal.pone.0028070
Sormani, R., Yao, L., Menand, B., Ennar, N., Lecampion, C., Meyer, C., et al. (2007). Saccharomyces cerevisiae FKBP12 binds Arabidopsis thaliana TOR and its expression in plants leads to rapamycin susceptibility. BMC Plant Biol. 7:26. doi: 10.1186/1471-2229-7-26
Thingholm, T. E., Jensen, O. N., Robinson, P. J., and Larsen, M. R. (2008). SIMAC (sequential elution from IMAC)., a phosphoproteomics strategy for the rapid separation of monophosphorylated from multiply phosphorylated peptides. Mol. Cell. Proteomics 7, 661–671. doi: 10.1074/mcp.M700362-MCP200
Thoreen, C. C., Chantranupong, L., Keys, H. R., Wang, T., Gray, N. S., and Sabatini, D. M. (2012). A unifying model for mTORC1-mediated regulation of mRNA translation. Nature 485, 109–113. doi: 10.1038/nature11083
Turck, F., Kozma, S. C., Thomas, G., and Nagy, F. (1998). A heat-sensitive Arabidopsis thaliana kinase substitutes for human p70s6k function in vivo. Mol. Cell. Biol. 18, 2038–2044. doi: 10.1128/MCB.18.4.2038
Turck, F., Zilbermann, F., Kozma, S. C., Thomas, G., and Nagy, F. (2004). Phytohormones participate in an S6 kinase signal transduction pathway in Arabidopsis. Plant Physiol. 134, 1527–1535. doi: 10.1104/pp.103.035873
Turkina, M. V., Klang Årstrand, H., and Vener, A. V. (2011). Differential phosphorylation of ribosomal proteins in Arabidopsis thaliana plants during day and night. PLoS ONE 6:e29307. doi: 10.1371/journal.pone.0029307
Wang, J., Wang, Y., Wang, Z., Liu, L., Zhu, X. G., and Ma, X. (2011). Synchronization of cytoplasmic and transferred mitochondrial ribosomal protein gene expression in land plants is linked to Telo-box motif enrichment. BMC Evol. Biol. 11:161. doi: 10.1186/1471-2148-11-161
Weijers, D., Franke-van Dijk, M., Vencken, R. J., Quint, A., Hooykaas, P., and Offringa, R. (2001). An Arabidopsis minute-like phenotype caused by a semi-dominant mutation in a RIBOSOMAL PROTEIN S5 gene. Development 128, 4289–4299.
Williams, A. J., Werner-Fraczek, J., Chang, I. F., and Bailey-Serres, J. (2003). Regulated phosphorylation of 40S ribosomal protein S6 in root tips of maize. Plant Physiol. 132, 2086–2097. doi: 10.1104/pp.103.022749
Yamaguchi, K., and Subramanian, A. R. (2000). The plastid ribosomal proteins. Identification of all the proteins in the 50 S subunit of an organelle ribosome (chloroplast). J. Biol. Chem. 275, 28466–28482. doi: 10.1074/jbc.M004350200
Yerlikaya, S., Meusburger, M., Kumari, R., Huber, A., Anrather, D., Costanzo, M., et al. (2016). TORC1 and TORC2 work together to regulate ribosomal protein S6 phosphorylation in Saccharomyces cerevisiae. Mol. Biol. Cell 27, 397–409. doi: 10.1091/mbc.E15-08-0594
Yu, Y., Yoon, S. O., Poulogiannis, G., Yang, Q., Ma, X. M., Villén, J., et al. (2011). Phosphoproteomic analysis identifies Grb10 as an mTORC1 substrate that negatively regulates insulin signaling. Science 332, 1322–1326. doi: 10.1126/science.1199484
Zhang, L. Y., Bai, M. Y., Wu, J., Zhu, J. Y., Wang, H., Zhang, Z., et al. (2009). Antagonistic HLH/bHLH transcription factors mediate brassinosteroid regulation of cell elongation and plant development in rice and Arabidopsis. Plant Cell 21, 3767–3780. doi: 10.1105/tpc.109.070441
Keywords: phosphorylation, plastid, proteomic, ribosome, RPS6, TOR kinase, transcriptomic, translatomic
Citation: Dobrenel T, Mancera-Martínez E, Forzani C, Azzopardi M, Davanture M, Moreau M, Schepetilnikov M, Chicher J, Langella O, Zivy M, Robaglia C, Ryabova LA, Hanson J and Meyer C (2016) The Arabidopsis TOR Kinase Specifically Regulates the Expression of Nuclear Genes Coding for Plastidic Ribosomal Proteins and the Phosphorylation of the Cytosolic Ribosomal Protein S6. Front. Plant Sci. 7:1611. doi: 10.3389/fpls.2016.01611
Received: 12 July 2016; Accepted: 12 October 2016;
Published: 07 November 2016.
Edited by:Wim Van den Ende, KU Leuven, Belgium
Reviewed by:Matthew Paul, Rothamsted Research (BBSRC), UK
Yan Xiong, Shanghai Center for Plant Stress Biology, China
Copyright © 2016 Dobrenel, Mancera-Martínez, Forzani, Azzopardi, Davanture, Moreau, Schepetilnikov, Chicher, Langella, Zivy, Robaglia, Ryabova, Hanson and Meyer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Christian Meyer, firstname.lastname@example.org