Advertisement
Research Article
Transparent Process
Open Access

Mitochondrial ubiquinone–mediated longevity is marked by reduced cytoplasmic mRNA translation

Marte Molenaars, Georges E Janssens, Toon Santermans, Marco Lezzerini, Rob Jelier, Alyson W MacInnes  Correspondence email, Riekelt H Houtkooper  Correspondence email
Published 31 August 2018. DOI: 10.26508/lsa.201800082

Abstract

Mutations in the clk-1 gene impair mitochondrial ubiquinone biosynthesis and extend the lifespan in Caenorhabditis elegans. We demonstrate here that this life extension is linked to the repression of cytoplasmic mRNA translation, independent of the alleged nuclear form of CLK-1. Clk-1 mutations inhibit polyribosome formation similarly to daf-2 mutations that dampen insulin signaling. Comparisons of total versus polysomal RNAs in clk-1(qm30) mutants reveal a reduction in the translational efficiencies of mRNAs coding for elements of the translation machinery and an increase in those coding for the oxidative phosphorylation and autophagy pathways. Knocking down the transcription initiation factor TATA-binding protein-associated factor 4, a protein that becomes sequestered in the cytoplasm during early embryogenesis to induce transcriptional silencing, ameliorates the clk-1 inhibition of polyribosome formation. These results underscore a prominent role for the repression of cytoplasmic protein synthesis in eukaryotic lifespan extension and suggest that mutations impairing mitochondrial function are able to exploit this repression similarly to reductions of insulin signaling. Moreover, this report reveals an unexpected role for TATA-binding protein-associated factor 4 as a repressor of polyribosome formation when ubiquinone biosynthesis is compromised.

Introduction

Slowing down mitochondrial metabolism is a well-known means of increasing the lifespan of multiple species (Dillin et al, 2002; Copeland et al, 2009; Houtkooper et al, 2013; Breitenbach et al, 2014). In C. elegans, one such model is the clk-1(qm30) mutant that harbors a deletion in the gene encoding a ubiquinone (UQ) biosynthesis enzyme (Felkai et al, 1999). UQ is a redox-active lipid that plays a central role in the electron transport chain and mitochondrial oxidative phosphorylation (OXPHOS) by carrying electrons from complexes I and II to complex III (Stefely & Pagliarini, 2017). In addition to lifespan extension and decreased respiration, clk-1(qm30) mutants are developmentally delayed and have reduced progeny (Felkai et al, 1999). Mammalian Clk1, known as Mclk-1 or Coq7, is also involved in the lifespan (Liu et al, 2005; Stepanyan et al, 2006; Takahashi et al, 2014). Transgenic mice, rescued from embryonic lethality via the transgenic expression of mouse clk-1, with reduced Clk1 expression live longer and have smaller bodies than WT mice, demonstrating a conserved role for clk-1 in longevity across species (Takahashi et al, 2014). However, it is not fully understood how the reduction of UQ biosynthesis extends the lifespan of these animals.

Until recently, CLK-1 protein was thought to reside exclusively in the mitochondria. However, a distinct nuclear form of CLK-1 has been reported to independently regulate lifespan by mitochondrial–nuclear retrograde signaling (Monaghan et al, 2015). Adding back a copy of the clk-1 gene that lacks the mitochondrial targeting signal (MTS) to the clk-1 deletion mutant led to a subtle but significant reduction in the lifespan extension phenotype observed in the full clk-1 deletion mutant (Monaghan et al, 2015). However, the existence and function of a nuclear form of CLK-1 remains controversial. Other studies fail to detect any nuclear localization of the MCLK1/CLK-1 proteins or any biological activity of a C. elegans CLK-1 protein devoid of an MTS (Liu et al, 2017).

The importance of communication between stressed mitochondria and the rest of the cell is becoming increasingly appreciated (Quiros et al, 2016). One such reaction to stress is carried out by the mitochondrial unfolded protein response (UPRmt) (Durieux et al, 2011; Baker et al, 2012; Houtkooper et al, 2013). GCN-2, an eIF2α kinase that modulates cytosolic protein synthesis, is phosphorylated in response to the activated UPRmt in clk-1(qm30) mutants in a manner required for lifespan extension (Baker et al, 2012). Another nonmitochondrial factor, TATA-binding protein-associated factor 4 (TAF-4), is a transcription factor that has been implicated in the extension of clk-1(qm30) mutant lifespan (Khan et al, 2013). TAF-4 was identified in an RNAi screen for transcription factors required for the lifespan extension of clk-1(qm30) mutants (Khan et al, 2013). Beyond clk-1, TAF-4 was also required for the lifespan extension phenotype of two mutants of the mitochondrial electron transport chain (isp-1 and tpk-1) (Khan et al, 2013). TAF-4 is a component of the transcription factor II D mRNA transcription complex and is best known for its role in transcriptional silencing in early embryogenesis by becoming sequestered in the cytoplasm because of phosphorylation by OMA-1 (Oocyte MAturation defective) (Walker et al, 2001). The loss of TAF-4 in C. elegans suppresses the lifespan extension phenotype induced by clk-1 mutations; however, the mechanism through which this happens is not known.

Reducing protein synthesis is another well-known means by which eukaryotes extend the lifespan (Hansen et al, 2007; Pan et al, 2007). Caloric restriction, amino acid reduction, and knocking down myriad factors involved in mRNA translation (such as ribosomal proteins, elongation/initiation factors, or tRNA synthetases) all increase the lifespan (Masoro, 2000; Min & Tatar, 2006; Hansen et al, 2007; Pan et al, 2007). The genetic or pharmacological inhibition of the mechanistic target of rapamycin (mTOR)/nutrient sensing pathway or insulin receptor signaling pathway is also marked by reduced mRNA translation rates and a suppression of polyribosome formation (Genolet et al, 2008; Stout et al, 2013). Proteomic analysis of the daf-2 mutants previously revealed a substantial reduction of ribosomal proteins, translation factors, and protein metabolism components coupled to a repression of polyribosome formation (Stout et al, 2013). Moreover, our previous work demonstrated that similarly to daf-2 mutants, mutations in clk-1 worms also result in a strong repression of polyribosome formation (Essers et al, 2015). This led us to hypothesize the existence of unexplored regulatory links between dysfunctional mitochondria and cytoplasmic protein synthesis that are contributing to eukaryotic lifespan extension. In particular, we investigated which subsets of mRNAs are preferentially translated in the longer lived C. elegans with impaired UQ biosynthesis.

Results

RNAseq of clk-1(qm30) strains shows major changes in transcription

To uncover how the transcriptome of C. elegans changes on deletion of clk-1 and the influence of the clk-1 mitochondrial localization signal, we performed next-generation sequencing of total RNAs (depleted of ribosomal RNA [rRNA]) isolated from worms at the L4 stage that precedes young adulthood. Total RNA pools depleted of rRNA isolated from clk-1(qm30) mutants were compared with (1) clk-1(qm30)+WT mutants rescued with the WT clk-1 and (2) clk-1(qm30)+nuc mutants with predominantly the nuclear form of clk-1 (Fig 1A) that was described in Monaghan et al (2015). The clk-1(qm30)+WT mutants have WT phenotypes and are rescued with fully functioning clk-1, whereas the clk-1(qm30)+nuc mutants have clk-1 only residing in the nucleus and not in the mitochondria (Monaghan et al, 2015). The principal component analysis (PCA) plot and correlation matrix show the expected clustering of the biological triplicates of each strain (Fig 1B–C). Major changes were observed in the RNA pool of the clk-1(qm30)+WT mutants compared with either clk-1(qm30) or clk-1(qm30)+nuc, with 9,626 and 8,685 genes being differentially expressed, respectively. Far fewer differentially expressed genes were measured when comparing clk-1(qm30) with clk-1(qm30)+nuc, with only 247 genes up-regulated and 571 genes down-regulated. The subtle changes in RNA expression are consistent with the observation that introducing the nuclear form of clk-1 in clk-1(qm30) results in only small phenotypic changes compared with the clk-1(qm30) mutants (Monaghan et al, 2015; Liu et al, 2017).

Figure 1.
Figure 1. RNA-seq of the clk-1(qm30) (± nuclear or WT clk-1) mutants.

(A) Graphical illustration of clk-1(qm30), clk-1(qm30)+nuc, and clk-1(qm30)+WT lines used for RNA-seq. mRNA was isolated from three biological replicates of each mutant, harvested at the L4 stage. (B) PCA of RNA libraries showing clear distinction between clk-1(qm30), clk-1(qm30)+nuc, and clk-1(qm30)+WT (N = 3 biological replicates for each strain). (C) Correlation matrix of RNA-seq samples shows the expected clustering of the biological triplicates of each strain. (D) Top 50 down-regulated GO terms (left graph) and up-regulated GO terms (right graph) in clk-1(qm30) versus clk-1(qm30)+WT (left column) and clk-1(qm30)+nuc versus clk-1(qm30)+WT (right column). Proportions of GO terms associated with translation and mRNA processing are depicted in yellow, development and reproduction in red, metabolism in green, and ion transport in blue.

Transcripts encoding the translation machinery are strongly reduced in clk-1(qm30)

We next identified the biological pathways that are altered in clk-1(qm30) worms using The Database for Annotation, Visualization and Integrated Discovery (DAVID) and Gene Ontology (GO) analysis (Huang da et al, 2009a; Huang da et al, 2009b). Among the most up-regulated processes in clk-1(qm30) are metabolic pathways (e.g., oxidation–reduction pathway) and ion transport (Fig 1D). A substantial proportion (25%) of significantly down-regulated processes in clk-1(qm30) mutants is involved in reproduction and development (Fig 1D). Developmental delay and reduced progeny have been extensively reported as linked to the longevity phenotype of clk-1(qm30) mutants (Larsen et al, 1995; Tissenbaum & Ruvkun, 1998). Remarkably, another 20% of down-regulated GO terms are involved in translation and mRNA processing (Fig 1D). In total, 38 genes encoding ribosomal proteins, initiation, and elongation factors are down-regulated in the clk-1(qm30) mutants compared with clk-1(qm30)+WT. In addition, 28 genes specifically involved in biogenesis and processing of ribosomes were down-regulated in the clk-1(qm30) worms. Similar down-regulation of ribosomal proteins and initiation/elongation factors was observed in the clk-1(qm30)+nuc mutants. These observations suggest that mRNA translation is repressed in the clk-1(qm30)/+nuc mutants and are in line with the repressed polysome profiles we observed before (Essers et al, 2015).

There is a minor role for nuclear clk-1 in transcriptional changes

Because the nuclear form of clk-1 was suggested to associate with the chromatin in the nucleus (Monaghan et al, 2015), we expected changes in the transcriptome in the mutant expressing exclusively the nuclear form of clk-1. However, differentially expressed genes in clk-1(qm30)+nuc were very similar to the differentially expressed genes clk-1(qm30) (Fig 2A). There were 3,720 genes similarly up-regulated in the clk-1(qm30)+nuc as in the total clk-1(qm30) mutants compared with clk-1(qm30)+WT and 4,397 genes similarly down-regulated in both strains compared with the WT control (Fig 2B). The 292 up-regulated genes and 275 down-regulated genes in the clk-1(qm30)+nuc strain that were not up- or down-regulated in the clk-1(qm30) strain did not lead to significantly clustered GO terms. The 704 genes up-regulated in specifically the clk-1(qm30) strain and not in the clk-1(qm30)+nuc strain were involved in biological processes such as ion transport and oxidation–reduction (Fig 2C). The 805 genes down-regulated in the clk-1(qm30) strain but not in clk-1(qm30)+nuc were involved in biological processes such as reproduction and response to DNA damage stimuli (Fig 2C). These processes also appeared when comparing the clk-1(qm30) mutants with WT worms.

Figure 2.
Figure 2. Differentially expressed genes in clk-1(qm30)+nuc compared with clk-1(qm30).

(A) Fold changes of differentially expressed genes in clk-1(qm30) versus clk-1(qm30)+WT plotted against fold changes of clk-1(qm30)+nuc versus clk-1(qm30)+WT shows a minor role for nuclear clk-1 in transcriptional changes. Colors of the individual data points correspond to the colors of the groups of genes in the Venn diagram in (B). (B) Venn diagram shows that most genes are similarly up-regulated (3720) and down-regulated (4397) in the clk-1(qm30)+nuc and clk-1(qm30) worms when comparing with clk-1+WT. There were 292 genes exclusively up-regulated in clk-1(qm30)+nuc (blue) and 704 genes exclusively in clk-1(qm30) (green) compared with clk-1+WT. Furthermore, 175 genes were down-regulated exclusively in clk-1(qm30)+nuc (red) and 805 genes down-regulated exclusively in clk-1(qm30) (orange) compared with clk-1+WT. Differential expressions in (A) and (B) are with a threshold-adjusted P value < 0.01. (C) Significant Cluster GO Enrichments (threshold Enrichment Score > 3) associated with the 704 genes specifically up-regulated (green) and 805 down-regulated (orange) in clk-1(qm30) strain. There were no significant Cluster GO Enrichments for genes exclusively up or down-regulated in clk-1(qm30)+nuc.

Translation rates are reduced in clk-1(qm30) worms

To provide a snapshot of the translational status in the clk-1(qm30) strains, we performed polysome profiling of C. elegans at the L4 stage. Polysome profiling fractionates total cell lysate over a sucrose gradient by density and enables the distinct measurement of 40S small ribosomal subunits, 60S large ribosomal subunits, 80S monosomes (an mRNA molecule with a single ribosome attached), and multiple polysomes (mRNA molecules with multiple ribosomes attached). Polysome profiles of clk-1(qm30) reveal a strong repression of both monosomal and polysomal peaks compared with clk-1(qm30)+WT, suggesting a global repression of mRNA translation (Fig 3). Again, the clk-1(qm30)+nuc mutants show similar yet slightly less dramatic differences compared with WT versus the total clk-1(qm30) mutant compared with WT (Fig 3). Because the clk-1(qm30)+nuc does not show striking differences in transcription compared with the total clk-1(qm30) mutant, we focused our remaining studies on comparing clk-1(qm30) with clk-1(qm30)+WT. The polysome profiles, together with the RNAseq results, show that clk-1(qm30) mutants are marked by repressed mRNA translation. Further validation of translational repression in clk-1(qm30) mutants is shown by quantitative PCR (qPCR) analysis, demonstrating that several mRNAs encoding initiation factors and ribosomal proteins are down-regulated in clk-1(qm30) mutants compared with WT worms (Fig S1).

Figure 3.
Figure 3. Repressed polysome profiles in the clk-1(qm30) (±nuc) mutants.

(A) Representative traces of polysome profiles of clk-1(qm30) strains harvested at the L4 stage, when lysate is normalized to total protein levels of 500 μg. The monosomal peak and polysomal peaks (P1–P6) are indicated. (B) Quantification of polysome peak sizes (AUC). The fold change is represented compared with P1 of the clk-1(qm30)+WT. All peaks of clk-1(qm30)+WT are significantly different from both clk-1(qm30)+nuc and clk-1(qm30). No peaks were significantly different between clk-1(qm30) and clk-1(qm30)+nuc. Error bars represent mean ± SD. Significance was tested with t test and P-values were adjusted to correct for multiple testing using the Holm–Sidak method, with α = 0.05.

Figure S1.
Figure S1. qPCR analysis of encoding initiation factors and RPs.

Relative expression levels determined by qRT–PCR (WT versus clk-1(qm30)) normalized by expression of pmp-3 as reference gene. Significance was tested with t test and P-values were adjusted to correct for multiple testing using the Holm–Sidak method, with α = 0.05.

Total RNA does not mirror highly translated polysomal RNA

In addition to providing a global overview of mRNA translation rates, polysome profiling allows the specific isolation of the RNAs associated with the polyribosomes for comparison with the total RNAs in different clk-1(qm30) strains (Fig 4A). The PCA plot and correlation matrix show clustering of polysomal and total RNA of the biological triplicates of each strain (Fig S2). We observed that changes in the polysomal RNA of clk-1(qm30) compared with the polysomal RNA of clk-1(qm30)+WT were not simply mirroring the changes observed in the total mRNA pools of these mutants (Fig 4B). For instance, there are 434 genes specifically up-regulated in the polysomal RNA pool carrying highly translated messages that are not up-regulated in the total RNA pool of clk-1(qm30) mutants (Fig 4C). A larger set of 1,657 genes is up-regulated in both polysomal and total RNA, whereas the largest part (2,767 genes) is exclusively up-regulated in the total RNA pool (Fig 4C). A similar pattern is observed with down-regulated genes in clk-1(qm30) mutants, with 466 genes detected in the polysomal RNA pool, 2,574 in both polysomal and total RNA pools, and 2,628 in the strictly total RNA pool (Fig 4C).

Figure 4.
Figure 4. Analysis of total and polysomal RNA clk-1(qm30) and clk-1(qm30)+WT.

(A) Schematic representation of polysomal fraction (blue) used for RNA isolation and total RNA (green) from C. elegans mutants clk-1(qm30) and clk-1(qm30)+WT that were used for RNAseq. (B) Fold changes of differentially expressed genes in the total RNA of clk-1(qm30) versus clk-1(qm30)+WT plotted against fold changes in polysomal RNA of clk-1(qm30) versus clk-1(qm30)+WT. Colors of the individual data points correspond to the colors of the groups of genes in the Venn diagram in (C). (C) Venn diagram showing 1657 genes were similarly up-regulated and 2574 genes down-regulated in both the polysomal as the total RNA in the clk-1(qm30) strain compared with clk-1(qm30)+WT. A total of 434 genes were exclusively up-regulated in the polysomal RNA (blue), whereas 2767 genes were exclusively up-regulated in the total RNA (green). Furthermore, 477 genes were exclusively down-regulated in the polysomal RNA (red) and 2628 genes were exclusively down-regulated in the total RNA. Differential expressions (B, C) are with a threshold-adjusted P value < 0.01. (D) Significant Cluster GO Enrichments (threshold Enrichment Score > 3) associated with the genes specifically up and down-regulated in groups of genes in (C) (colors of the bars correspond again to the color of the groups in (C)). (E) Fib-1 is reduced in both the total and polysomal RNA of clk-1(qm30) compared with clk-1(qm30)+WT. (F, G) FIB-1 levels are reduced in clk-1(qm30) worms compared with N2 worms. Error bars represent mean ± SD, significance was tested with t test, ***P < 0.0005.

Figure S2.
Figure S2. Global analysis of polysomal and total RNA-seq.

(A) Correlation matrix of RNA-seq samples shows the expected clustering of total and polysomal RNA of the biological triplicates of each strain. (B) PCA of RNA libraries showing clear distinction between polysomal and total RNA of each of the strains (N = 3 biological replicates for each strain).

OXPHOS transcripts are highly enriched in polysomal RNA of clk-1(qm30) strains

To determine the functional relevance of the differentially expressed RNAs in total and polysomal RNA pools, we performed DAVID GO analysis. For instance, ion transport processes are up-regulated in the total mRNA of clk-1(qm30) mutants. However, the polysomal expression data indicate that despite this higher abundance in the total pool, mRNAs coding for ion transport are not more highly translated in clk-1(qm30) mutants (Fig 4D). In contrast, mRNAs involved in metabolic processes such as oxidation–reduction, glycolysis, and lipid metabolism are substantially up-regulated in both the total and the polysomal RNA pools (Fig 4D). This suggests a metabolic compensation for the reduction of OXPHOS capacity by glycolysis and lipid metabolism in the clk-1(qm30) mutants. Processes that were exclusively up-regulated in the polysomal RNA pool (and thus exclusively highly translated) in clk-1(qm30) mutants are involved in ATP hydrolysis, protein import into mitochondria, and mitochondrial electron transport (Fig 4D). More specifically, 50 transcripts encoding proteins involved in mitochondrial OXPHOS were enriched in clk-1(qm30) polysomes (Table 1). This up-regulation does not seem complex-specific because mRNAs coding for proteins of all OXPHOS complexes (I–V) were enriched in clk-1(qm30) polysomes. These data suggest a model whereby proteins involved in ATP production and the import of these proteins into the mitochondria are highly translated despite lower overall expression levels in clk-1(qm30) mutants compared with WT to attempt to overcome their energy deficit.

View this table:
Table 1.

Identified polysomal transcripts involved in OXPHOS that, compared with clk-1(qm30)+WT, are enriched in clk-1(qm30).

RNAs coding for ribosomal proteins are substantially reduced in clk-1(qm30) polysomes

An even more substantial down-regulation of transcripts involved in translation and ribosome biogenesis processes was revealed in the polysomal RNA pool of clk-1(qm30) deficient worms compared with clk-1(qm30)+WT, confirming the repression of translational machinery proteins (Table 2). In the polysomal RNA pool of clk-1(qm30) mutants, 62 ribosomal proteins, 16 elongation/initiation factors, 14 amino-acyl tRNA synthetases, and 10 other proteins involved in mRNA processing were down-regulated (Table 2). These results are reminiscent of the proteomic analysis demonstrating a significant reduction of ribosomal proteins and translation factors in the daf-2 insulin receptor mutants (Stout et al, 2013). Furthermore, a subset of these transcripts coding for proteins involved in translation (21 ribosomal proteins, 11 amino-acyl tRNA synthetases, and eight other elongation/initiation factors) was even down-regulated in the total pool of RNA. This indicates that the down-regulation of translation machinery is partially regulated in a transcriptional manner. Interestingly, we observed strong down-regulation of fib-1, a gene coding for a protein involved in the regulation and maturation of rRNA, in both total and polysomal RNA (Fig 4E). The FIB-1 protein resides in the nucleoli where ribosome biogenesis initiates. It was previously shown that small nucleoli are a cellular hallmark of longevity and metabolic health (Buchwalter & Hetzer, 2017; Tiku et al, 2017). A FIB-1::GFP reporter strain for nucleoli revealed reduced FIB-1 levels in other long-lived C. elegans strains (Tiku et al, 2017). In line with this and the RNA-seq data presented here, the clk-1(qm30) worms also show reduced FIB-1 levels compared with WT N2 worms that is accompanied by reduced nucleolar size (Fig 4F and G). Taken together, it appears that an attempt to overcome an energy deficit in clk-1(qm30) mutants is coupled to the suppression of ribosome biogenesis and mRNA translation, which would conserve flagging energy for other critical cell processes.

View this table:
Table 2.

Identified polysomal transcripts involved in mRNA translation that, compared with clk-1(qm30)+WT, are reduced in clk-1(qm30).

Differential translational efficiency (TE) of mRNA reveals a role for TOR signaling

Because there are mRNAs that are selectively being translated in the clk-1(qm30) mutant or the WT independently of their abundance in the total mRNA pool, we calculated the TE of individual mRNA transcripts by taking the ratio between total RNA and polysomal RNA (Fig 5A). DAVID GO analysis on transcripts with high or low TE revealed similar biological processes as observed in the groups of genes specifically up- and down-regulated in the polysomal RNA presented in Fig 4D (Fig 5B). The list of genes with significantly different TEs was then analyzed using the “worm longevity pathway” in Kyoto Encyclopedia of Genes and Genomes (Ogata et al, 1999) using Pathview (Luo & Brouwer, 2013), to investigate the overlap in longevity pathways known in C. elegans. Besides hits in the UPRmt and stress resistance, both already described in clk-1(qm30) mutants (Durieux et al, 2011), the results also revealed possible involvement of the mTOR signaling pathway in the clk-1(qm30) mutants. Hence, we also checked the mTOR signaling pathway in Kyoto Encyclopedia of Genes and Genomes Pathview. In C. elegans, homologs of TOR pathway genes all extend the lifespan upon siRNA knockdown (Korta et al, 2012). Interestingly, our RNA-seq data show that let-363 (homologue of TOR) is preferentially translated in WT but not in clk-1(qm30) mutants (Fig 5C and D). Other members of the TOR complex 1 (TORC1), including daf-15 (RAPTOR), clk-2 (TEL2), and R10H10.7 (TTI1), also showed decreased TE in clk-1(qm30) mutants (Fig 5C and D). Moreover, our results show that a key player in autophagy, lgg-1, is preferentially translated in clk-1(qm30) mutants but not in the clk-1(qm30)+WT (Fig 5C and D).

Figure 5.
Figure 5. Translation efficiency (TE) of transcripts in clk-1(qm30) worms.

(A) Volcano plot of log2 fold change of TE (total RNA:polysomal RNA) of transcripts clk-1(qm30):clk-1(qm30)+WT. The blue data points represent transcripts with high TE in clk-1(qm30) worms being shifted from the total RNA to the highly translated polysomal RNA. The red data points represent transcripts that have low TE in clk-1(qm30) being shifted from polysomal RNA to the total RNA. Differentially translationally regulated genes in red and blue with threshold-adjusted P value < 0.01. (B) Significant Cluster GO Enrichments (threshold Enrichment Score > 3) associated with significantly different TEs (colors of the bars correspond again to the data points in (A)). (C) TE of individual transcripts involved in TOR pathway in clk-1(qm30)+WT (left bar) clk-1(qm30) (middle bar) and their ratio (right bar). Colors of the ratio bars correspond again to the data points in (A). Clk-1(qm30) worms slow reduced TE of let-363, daf-15, clk-2, and R10H10.7 and an increased TE for lgg-1. (D) Schematic overview of transcripts involved in TOR pathway that have altered TE in clk-1(qm30) versus clk-1(qm30)+WT represented in (C) and their involvement in longevity.

taf-4 is required for repressed polysome profiles in clk-1(qm30)

TAF-4 is one of few transcription factors reported to be critical for the clk-1(qm30) longevity phenotype (Khan et al, 2013). To test if TAF-4 has a role in polysome formation in clk-1(qm30) worms, we performed RNAi of taf-4 in the clk-1(qm30) mutants. Here, we show that RNAi of taf-4 significantly increases the monosomal and polysomal peaks of clk-1(qm30) compared with the control empty vector HT115 (Fig 6A and C), suggesting this transcription factor plays a direct role in the suppression of polyribosome formation upon the loss of UQ biosynthesis. This increase is not observed in the polysome profiles upon taf-4 RNAi in the clk-1(qm30)+WT (Fig 6B). These data suggest that TAF-4 loss partially restores polyribosome formation in clk-1(qm30) mutants similar to the previously reported attenuation of the longevity phenotype (Khan et al, 2013). Because TAF-4 is a transcription factor, we looked at changes in expression of specific genes involved in translation upon taf-4 RNAi in the clk-1(qm30) mutants. We observed significant up-regulation of the translation initiation factor eif-1 when taf-4 is knocked down in clk-1(qm30) worms, suggesting that TAF-4 has a direct role in repressing eif-1 expression (Fig 6D). Eif-1 is important for translation initiation, and yeast eif-1 mutants (sui1−) show a reduction of polysomes compared with WT (Yoon & Donahue, 1992). Moreover, the reduction of eif-1 extends lifespan in N2 worms (Wang et al, 2010). Taken all together, these results suggest that the reduction of the translation machinery observed in clk-1(qm30) mutants may at least in part be driven by TAF-4 repression, and may help to explain the previously described link between TAF-4– and clk-1(qm30)– mediated lifespan extension (Khan et al, 2013).

Figure 6.
Figure 6. RNAi of taf-4 partially restores repressed peaks in clk-1(qm30) mutants.

(A) Representative traces of polysome profiles of clk-1(qm30) worms fed HT115 bacteria transformed with the empty vector (EV) or expressing taf-4 RNAi when harvested at the L4 stage, when lysate is normalized to total protein levels of 500 μg. The monosomal peak and polysomal peaks (P1–P5) are indicated. (B) Representative traces of polysome profiles of clk-1(qm30)+WT worms fed HT115 bacteria transformed with the EV or expressing taf-4 RNAi, when lysate is normalized to total protein levels of 500 μg. The monosomal peak and polysomal peaks are indicated. (C) Quantification of polysome peak sizes (AUC). The fold change is represented compared with P1 of the clk-1(qm30)+WT fed with HT115 bacteria. Polysomal peaks P1 and P2 were significantly different between clk-1(qm30)+WT fed with HT115 and taf-4 RNAi bacteria. No peaks were significantly different between clk-1(qm30)+WT fed with HT115 and taf-4 RNAi bacteria. Error bars represent mean ± SD. Significance was tested with t test and P-values were adjusted to correct for multiple testing using the Holm–Sidak method, with α = 0.05. (D) Relative expression levels determined by qRT–PCR in clk-1(qm30)(+WT) fed with either the EV or taf-4 RNAi bacteria. Expression levels were normalized using the geometrical mean of reference genes Y45F10D.4, tba-1, and csq-1. Significance was tested using one-way ANOVA with Sidak's multiple comparisons test. Error bars represent mean ± SEM. **P < 0.01, significance indicated only between clk-1(qm30) EV and clk-1(qm30) taf-4 RNAi.

Discussion

This study was initiated following our previous observation that the polysome profiles of clk-1(qm30) mutants revealed the same repression of polysomes as in the daf-2 mutants (Essers et al, 2015). We continued this line of inquiry by isolating and sequencing polysomal and total RNAs from clk-1(qm30) mutants and compared them with clk-1(qm30)+WT worms. These results demonstrated that compared with clk-1(qm30)+WT, very few mRNAs coding for components of the translation machinery such as initiation factors, elongation factors, and ribosomal proteins are present on clk-1(qm30) polyribosomes. In contrast, the clk-1(qm30) polyribosomes reveal a strong preference for translating RNAs coding for OXPHOS proteins. The reduction in the TE of mRNAs coding for the translation machinery is highly reminiscent of the proteomic landscape we described in the daf-2 mutants, where there is a scarcity of ribosomal proteins (RPs) yet also no appreciable reduction in total RP mRNA expression (Depuydt et al, 2013; Essers et al, 2015). In addition, the small nucleoli in clk-1(qm30) worms confirming reduced ribosome biogenesis is a hallmark of cellular aging that was previously observed in other longevity mutants (Tiku et al, 2017). Taken all together, the results suggest that reduced ribosome biogenesis and shutting down protein synthesis by altering the TE of mRNAs coding for the translation machinery is a common mechanism used by several different systems that drive lifespan extension.

Our results go on to suggest that the proposed nuclear form of CLK-1 has a relatively small effect on alleviating the polyribosome reduction observed in the clk-1(qm30) mutants. The subtle differences in lifespan and development previously described (Monaghan et al, 2015) might be explained by reestablishing expression of some reproduction genes in the clk-1(qm30)+nuc mutant, which are well known to associate with lifespan (Hsin & Kenyon, 1999; Monaghan et al, 2015). Because it has been shown that under certain conditions proteins may be imported into the mitochondria without any MTS (Ruan et al, 2017), it is possible that a small proportion of the nuclear CLK-1 enters the mitochondria without an MTS. The present study does not disprove the existence of nuclear form of CLK-1; however, an obvious function of a nuclear CLK-1 is not apparent with these data.

The inhibition of respiration in the clk-1(qm30) mutants triggers a retrograde response by up-regulating cell-protective and metabolic genes (Cristina et al, 2009). Our work is in line with these findings and in this study, we confirm the importance of retrograde communication between stressed mitochondria and the rest of the cell. In Drosophila, it was shown that dietary restriction–induced longevity was regulated by activation of the translational repressor 4E-BP, which is the eukaryotic translation initiation factor 4E binding protein (Zid et al, 2009). This resulted in increased ribosomal loading of nuclear-encoded OXPHOS mRNAs. This high TE of mitochondrial proteins was crucial for the dietary restriction–induced lifespan extension because inhibition of individual mitochondrial subunits of OXPHOS complexes diminished the lifespan extension (Zid et al, 2009). Another study showed that reduction of cytosolic protein synthesis could suppress ageing-related mitochondrial degeneration in yeast (Wang et al, 2008). This postponed mitochondrial degeneration upon reducing cytoplasmic protein synthesis together with our results suggesting a robust strategy where mitochondria communicate with the cytoplasmic translation machinery to equilibrate cellular energy generated by one compartment and used in another. Interestingly, the mitochondrial dysfunction observed in young Mclk1+/− mice was shown to paradoxically result in an almost complete protection from the age-dependent loss of mitochondrial function later in life (Lapointe et al, 2009). The increased TE of mitochondrial mRNAs we observed here in the mitochondrial clk-1(qm30) mutants could induce this physiological state that ultimately develops long-term beneficial effects for aging.

Our results found that mRNAs coding for OXPHOS complexes I–V substantially enriched on polysomes of clk-1(qm30) worms. Given the marked reduction of total polysomes in these mutants, we can assume that most of the few remaining polysomes are present solely to maintain OXPHOS complexes at viable levels. This translational up-regulation of OXPHOS mRNAs on ribosomes was also recently observed by ribosome profiling in the context of host shutoff during vaccinia virus infection in human cells (Dai et al, 2017). This host shutoff facilitates resource availability for viral evasion and is coupled to a global inhibition of protein synthesis in the host (Bercovich-Kinori et al, 2016). In addition, the authors propose that this translational up-regulation of OXPHOS transcripts is due to their relatively short 5′ UTRs, which are known to play a regulatory role in TE (Araujo et al, 2012; Dai et al, 2017). A similar mechanism coupling energy metabolism and translational control is observed in the present study in the context of longevity.

Besides OXPHOS transcripts, this study also identified an increased TE for the autophagy gene lgg-1, the C. elegans orthologue of LC3 in humans. Autophagy is an indispensable function in organisms undergoing lifespan extension. In the absence of bec-1, the C. elegans orthologue of the autophagy gene APG6/VPS30/beclin1, daf-2–mutant worms fail to undergo lifespan extension (Melendez et al, 2003). Reducing bec-1 in worms subject to caloric restriction (eat-2 mutants) or inhibition of the mTOR pathway (reducing let-363 [TOR] or daf-15 [RAPTOR] expression) also has the same effect of preventing lifespan extension (Hansen et al, 2008). These findings are in line with our observation that mTOR is translationally repressed in clk-1(qm30) mutants (Fig 5D) and fits well with the notion that a reduction of growth can often lead to lifespan extension, in line with current theories proposing hyperfunction as a driving force in aging (Blagosklonny, 2006; Blagosklonny & Hall, 2009). Furthermore, our results suggest a novel mechanism for inducing autophagy in systems undergoing lifespan extension by affecting the TE of mRNAs coding for crucial components of the mTOR-autophagy regulatory axis.

Finally, we found an unexpected role for TAF-4 as a repressor of polyribosome formation when UQ biosynthesis is compromised. Although TAF-4 is essential for all RNA transcriptions in the early embryo in C. elegans (Walker et al, 2001), this requirement is different in later developmental stages, and worms exposed to taf-4 RNAi starting from L1 throughout their remaining life grow normally (Khan et al, 2013). The authors propose a role for TAF-4 in the retrograde response model of longevity in mutants with dysfunctional mitochondria. In this model, dysfunctional mitochondria signal via cAMP response element-binding proteins to activate TAF-4 in a way that increases life-extending gene expression. Here, we show that one of the consequences of TAF-4 loss in clk-1(qm30) mutants, and perhaps also in isp-1 and tpk-1 mutants, is a failure of the system that shuts down cytosolic mRNA translation in response to mitochondria impairment. This unexpected ability of TAF-4 to control protein synthesis, by regulating the expression of eif-1, will be very interesting for further study.

This study demonstrates that molecular pathways of long-lived clk-1(qm30) mutants with dysfunctional mitochondria converge with those of other distinctly different long-lived mutants to reduce cytoplasmic mRNA translation. It also reveals a striking change in the transcripts that the polyribosomes prefer to translate when mitochondria are compromised. Ribosomal proteins and other protein parts of the translation machinery are very marginally translated, whereas the mitochondrial OXPHOS genes are still robustly translated despite a substantial reduction in the number of polysomes. We propose that this mechanism is likely recapitulated in many more models of lifespan extension and reaffirms studies that aim to extend human lifespan by reducing protein synthesis.

Materials and Methods

Nematode growing and conditions

C. elegans strains used MQ130 clk-1(qm30), OL0100 ukSi1[pclk-1::clk-1::gfp, cb-unc-119(+)], and OL0120 ukSi2[pclk-1::clk-1 ΔMTS (13-187)::gfp; cb-unc-119(+)] was kindly provided by AJ Whitmarsh and G Poulin (Monaghan et al, 2015) and cguIs001 (FIB-1::GFP) was kindly provided by SJ Lo and T Ma (Lee et al, 2010). Worms were cultured at 20°C on nematode growth medium agar plates seeded with OP50 strain Escherichia coli. Synchronized worms were harvested and snap-frozen at the L4 larval stage for either total mRNA isolation or polysome profiling. For RNAi knockdown experiments, synchronized eggs were plated on nematode growth media containing IPTG (containing 2 mM IPTG) with a bacterial lawn of either E. coli HT115 (RNAi control strain, containing an empty vector) or taf-4 RNAi bacteria.

Polysome profiling

Gradients of 17–50% sucrose (11 ml) in gradient buffer (110 mM KAc, 20 mM MgAc2, and 10 mM Hepes, pH 7.6) were prepared on the day before use in thin-walled, 13.2-ml, polyallomer 14 89-mm centrifuge tubes (Beckman Coulter). Nematodes were lysed in 500 ml polysome lysis buffer (gradient buffer containing 100 mM KCl, 10 mM MgCl2, 0.1% NP-40, 2 mM DTT, and RNaseOUT [Thermo Fisher Scientific]) using a Dounce homogenizer. The samples were centrifuged at 3,500 g for 10 min to remove debris and the supernatant was subjected to bicinchoninic acid assay (BCA) protein assay. In all, 500 mg of total protein for each sample was loaded on top of the sucrose gradients. The gradients were ultracentrifuged for 2 h at 40,000 rpm in a SW41Ti rotor (Beckman-Coulter). The gradients were displaced into a UA6 absorbance reader (Teledyne ISCO) using a syringe pump (Brandel) containing 60% sucrose. Absorbance was recorded at an optical density of 254 nm. All steps of this assay were performed at 4°C or on ice and all chemicals came from Sigma-Aldrich unless stated otherwise. Polysome peaks were quantified by measuring the area under the curve (AUC) in ImageJ. Significance was tested with t test, and P-values were adjusted to correct for multiple testing using the Holm–Sidak method, with α = 0.05.

RNA sequencing

Total/polysomal RNA isolation

For isolation of total RNA, worms were homogenized with a 5-mm steel bead using a TissueLyser II (Qiagen) for 5 min at a frequency of 30 times/s in the presence of TRIzol (Invitrogen), and the isolation was continued according to the manufacturer's protocol. Polysomal fractions from two experiments were pooled and mRNA was extracted using TRIzol LS (Invitrogen) according to the manufacturer's protocol. Contaminating genomic DNA was removed using RNase-Free DNase (Qiagen) and samples were cleaned up with the RNeasy MinElute Cleanup kit (Qiagen) and sequenced by GenomeScan BV at a 20 M read depth.

Library preparation

The samples were processed for Illumina using the NEBNext Ultra Directional RNA Library Prep kit (NEB: New England Biolabs [NEB] #E7420) according to the manufacturer's protocol. Briefly, rRNA was depleted from total RNA using the rRNA depletion kit (NEB #E6310). After fragmentation of the rRNA-reduced RNA, a cDNA synthesis was performed to ligate with the sequencing adapters and PCR amplification of the resulting product. Quality and yield after sample preparation was measured with the Fragment Analyzer. Size of the resulting products was consistent with the expected size distribution (a broad peak between 300 and 500 bp). Clustering and DNA sequencing using the Illumina cBot and HiSeq 4000 was performed according to the manufacturer's protocol with a concentration of 3.0 nM of DNA. HiSeq control software HCS v3.4.0, image analysis, base calling, and quality check were performed with the Illumina data analysis pipeline RTA v2.7.7 and Bcl2fastq v2.17. RNA-seq data have been deposited in the ArrayExpress database at European Molecular Biology Laboratory: European Bioinformatics Institute under accession number: E-MTAB-7166.

Bioinformatics analysis

Mapping of the reads

The quality of the reads in the fastq files was confirmed using FastQC version 0.11.4 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc). FastQC performs multiple quality tests and provides the user with warnings about possible problems with the raw data. Paired end reads were mapped to the C. elegans genome using HISAT2 version 2.1.0(Kim et al, 2015). The WBcel235 genome assembly was downloaded from WormBase and used as reference genome for the mapping. After successful mapping, counts per gene were extracted from the .SAM files with HTSeq version 0.9.1 (Anders et al, 2015). The -stranded=reverse setting was used, conforming with the NEB Ultra Directional RNA Library Prep kit procedure.

Statistical analysis

The statistical analysis was performed in R using DESeq2 version 1.16.1 (Love et al, 2014). DESeq2 models the raw counts with the negative binomial distribution and normalizes between samples by calculating size factors using the median ratio method (Anders & Huber, 2010). The statistical test is based on a generalized linear model (GLM) using the negative binomial model as an error distribution model. To more accurately model dispersion parameters, the trend of dispersion to abundance is taken into account using an empirical Bayes procedure. The significance of the GLM coefficients is determined using a Wald test. The P-values are adjusted to correct for multiple testing. In the differential translational regulation analysis, the fitted GLM took the following form:log(μi)=β0+β1Xsample+β2Xtype+β3Xsample Xtype.(1)

Here, the Wald test was performed on the interaction parameter 3 to infer significance. The log2 fold change calculated in this case corresponds to the following:log2 fold change=log 2(strain Apolysomal/strain Atotalstrain Bpolysomal/strain Btotal).(2)where in Equation (1), the sample refers to the strain and type to polysomal or total RNA. Experimental and statistical normalizations were performed for each condition individually, so counts for a specific gene in the polysomal dataset can be possibly higher than in the total RNA dataset, even though the first is technically a subset of the latter. Still, we can test if genes are highly translated in one strain when it is not in the other.

Fluorescence microscopy imaging

For imaging, synchronized L4 stage N2 and clk-1(qm30) worms carrying FIB-1::GFP were placed on top of a 2% agarose pad on slides and anaesthetized using 1 mM levamisole. Nucleoli in the pharynx were imaged using 100× magnification (PLAN APO 100×/1.40/OIL Phaco3) with a Leica DM6 microscope and Leica LAS-X software. Total and average nucleolar area per worm was quantified using ImageJ.

qPCR of total and polysomal RNA

cDNA synthesis was performed with 1 μg of total RNA using the QuantiTect Reverse Transcription Kit (Qiagen). LightCycler 480 SYBR Green I Master (Roche) was used for qPCR analysis, primers are listed in Table 3. Data were analyzed with LightCycler 480 software release 1.5 and LinRegPCR version 2015.3. Relative gene expression was normalized to the geometric mean of reference genes Y45F10D.4, tba-1, and csq-1.

View this table:
Table 3.

Primers for qPCR.

Acknowledgements

The authors thank AJ Whitmarsh and G Poulin (University of Manchester) for providing the MQ130, OL0100, and OL0120 C. elegans clk-1(qm30) strains. We also thank SJ Lo and T Ma (Chang Gung University) for providing the FIB-1::GFP strain. Work in the Houtkooper group is financially supported by an European Research Council Starting grant (638290) and a VIDI grant from ZonMw (91715305) and the Velux Stiftung. Work in the Jelier group is financially supported by a Katholieke Universiteit Leuven grant C14/16/060. AW MacInnes and M Lezzerini are supported by E-Rare-2, the ERA-Net for Research on Rare Diseases (ZonMW 40-44000-98-1008). GE Janssens is supported by a 2017 Federation of European Biochemical Society long-term fellowship. We thank the financial support of the European Cooperation in Science and Technology initiative called “Group of C. elegans New Investigators in Europe (GENiE; BM1408).”

Author Contributions

  • M Molenaars: conceptualization, formal analysis, investigation, visualization, methodology, and writing—original draft.

  • GE Janssens: data curation, formal analysis, funding acquisition, investigation, visualization, and writing—review and editing.

  • T Santermans: data curation, formal analysis, investigation, and writing—review and editing.

  • M Lezzerini: formal analysis, investigation, and writing—review and editing.

  • R Jelier: data curation, formal analysis, supervision, funding acquisition, investigation, and writing—review and editing.

  • AW MacInnes: conceptualization, formal analysis, supervision, funding acquisition, investigation, and writing—original draft.

  • RH Houtkooper: conceptualization, resources, formal analysis, supervision, funding acquisition, investigation, and writing—original draft.

Conflict of Interest Statement

The authors declare that they have no conflict of interest.

  • Received April 30, 2018.
  • Revision received August 23, 2018.
  • Accepted August 23, 2018.
Creative Commons logoCreative Commons logohttps://creativecommons.org/licenses/by/4.0/

This article is available under a Creative Commons License (Attribution 4.0 International, as described at https://creativecommons.org/licenses/by/4.0/).

References

  1. ↵
    1. Anders S,
    2. Huber W
    (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106. doi:10.1186/gb-2010-11-10-r106
    OpenUrlCrossRefPubMed
  2. ↵
    1. Anders S,
    2. Pyl PT,
    3. Huber W
    (2015) HTSeq: A Python framework to work with high-throughput sequencing data. Bioinformatics 31: 166–169. doi:10.1093/bioinformatics/btu638
    OpenUrlCrossRefPubMed
  3. ↵
    1. Araujo PR,
    2. Yoon K,
    3. Ko D,
    4. Smith AD,
    5. Qiao M,
    6. Suresh U,
    7. Burns SC,
    8. Penalva LO
    (2012) Before it gets started: Regulating translation at the 5′ UTR. Comp Funct Genomics 2012: 475731. doi:10.1155/2012/475731
    OpenUrlCrossRefPubMed
  4. ↵
    1. Baker BM,
    2. Nargund AM,
    3. Sun T,
    4. Haynes CM
    (2012) Protective coupling of mitochondrial function and protein synthesis via the eIF2alpha kinase GCN-2. PLoS Genet 8: e1002760. doi:10.1371/journal.pgen.1002760
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bercovich-Kinori A,
    2. Tai J,
    3. Gelbart IA,
    4. Shitrit A,
    5. Ben-Moshe S,
    6. Drori Y,
    7. Itzkovitz S,
    8. Mandelboim M,
    9. Stern-Ginossar N
    (2016) A systematic view on influenza induced host shutoff. eLife 5: e18311. doi:10.7554/elife.18311
    OpenUrlCrossRef
  6. ↵
    1. Blagosklonny MV
    (2006) Aging and immortality: Quasi-programmed senescence and its pharmacologic inhibition. Cell Cycle 5: 2087–2102. doi:10.4161/cc.5.18.3288
    OpenUrlCrossRefPubMed
  7. ↵
    1. Blagosklonny MV,
    2. Hall MN
    (2009) Growth and aging: A common molecular mechanism. Aging 1: 357–362. doi:10.18632/aging.100040
    OpenUrlCrossRefPubMed
  8. ↵
    1. Breitenbach M,
    2. Rinnerthaler M,
    3. Hartl J,
    4. Stincone A,
    5. Vowinckel J,
    6. Breitenbach-Koller H,
    7. Ralser M
    (2014) Mitochondria in ageing: There is metabolism beyond the ROS. FEMS Yeast Res 14: 198–212. doi:10.1111/1567-1364.12134
    OpenUrlCrossRefPubMed
  9. ↵
    1. Buchwalter A,
    2. Hetzer MW
    (2017) Nucleolar expansion and elevated protein translation in premature aging. Nat Commun 8: 328. doi:10.1038/s41467-017-00322-z
    OpenUrlCrossRef
  10. ↵
    1. Copeland JM,
    2. Cho J,
    3. Lo T Jr.,
    4. Hur JH,
    5. Bahadorani S,
    6. Arabyan T,
    7. Rabie J,
    8. Soh J,
    9. Walker DW
    (2009) Extension of Drosophila life span by RNAi of the mitochondrial respiratory chain. Curr Biol 19: 1591–1598. doi:10.1016/j.cub.2009.08.016
    OpenUrlCrossRefPubMed
  11. ↵
    1. Cristina D,
    2. Cary M,
    3. Lunceford A,
    4. Clarke C,
    5. Kenyon C
    (2009) A regulated response to impaired respiration slows behavioral rates and increases lifespan in Caenorhabditis elegans. PLoS Genet 5: e1000450. doi:10.1371/journal.pgen.1000450
    OpenUrlCrossRefPubMed
  12. ↵
    1. Dai A,
    2. Cao S,
    3. Dhungel P,
    4. Luan Y,
    5. Liu Y,
    6. Xie Z,
    7. Yang Z
    (2017) Ribosome profiling reveals translational upregulation of cellular oxidative phosphorylation mRNAs during Vaccinia virus-induced host shutoff. J Virol 91: e01858-16. doi:10.1128/jvi.01858-16
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. Depuydt G,
    2. Xie F,
    3. Petyuk VA,
    4. Shanmugam N,
    5. Smolders A,
    6. Dhondt I,
    7. Brewer HM,
    8. Camp DG 2nd.,
    9. Smith RD,
    10. Braeckman BP
    (2013) Reduced insulin/insulin-like growth factor-1 signaling and dietary restriction inhibit translation but preserve muscle mass in Caenorhabditis elegans. Mol Cell Proteomics 12: 3624–3639. doi:10.1074/mcp.m113.027383
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Dillin A,
    2. Hsu AL,
    3. Arantes-Oliveira N,
    4. Lehrer-Graiwer J,
    5. Hsin H,
    6. Fraser AG,
    7. Kamath RS,
    8. Ahringer J,
    9. Kenyon C
    (2002) Rates of behavior and aging specified by mitochondrial function during development. Science (New York, NY) 298: 2398–2401. doi:10.1126/science.1077780
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Durieux J,
    2. Wolff S,
    3. Dillin A
    (2011) The cell-non-autonomous nature of electron transport chain-mediated longevity. Cell 144: 79–91. doi:10.1016/j.cell.2010.12.016
    OpenUrlCrossRefPubMed
  16. ↵
    1. Essers PB,
    2. Nonnekens J,
    3. Goos YJ,
    4. Betist MC,
    5. Viester MD,
    6. Mossink B,
    7. Lansu N,
    8. Korswagen HC,
    9. Jelier R,
    10. Brenkman AB, et al.
    (2015) A long noncoding RNA on the ribosome is required for lifespan extension. Cell Rep 10: 307–452. doi:10.1016/j.celrep.2014.12.029
    OpenUrlCrossRef
  17. ↵
    1. Felkai S,
    2. Ewbank JJ,
    3. Lemieux J,
    4. Labbe JC,
    5. Brown GG,
    6. Hekimi S
    (1999) CLK-1 controls respiration, behavior and aging in the nematode Caenorhabditis elegans. EMBO J 18: 1783–1792. doi:10.1093/emboj/18.7.1783
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Genolet R,
    2. Araud T,
    3. Maillard L,
    4. Jaquier-Gubler P,
    5. Curran J
    (2008) An approach to analyse the specific impact of rapamycin on mRNA-ribosome association. BMC Med Genomics 1: 33. doi:10.1186/1755-8794-1-33
    OpenUrlCrossRefPubMed
  19. ↵
    1. Hansen M,
    2. Chandra A,
    3. Mitic LL,
    4. Onken B,
    5. Driscoll M,
    6. Kenyon C
    (2008) A role for autophagy in the extension of lifespan by dietary restriction in C. elegans. PLoS Genet 4: e24. doi:10.1371/journal.pgen.0040024
    OpenUrlCrossRefPubMed
  20. ↵
    1. Hansen M,
    2. Taubert S,
    3. Crawford D,
    4. Libina N,
    5. Lee SJ,
    6. Kenyon C
    (2007) Lifespan extension by conditions that inhibit translation in Caenorhabditis elegans. Aging Cell 6: 95–110. doi:10.1111/j.1474-9726.2006.00267.x
    OpenUrlCrossRefPubMed
  21. ↵
    1. Houtkooper RH,
    2. Mouchiroud L,
    3. Ryu D,
    4. Moullan N,
    5. Katsyuba E,
    6. Knott G,
    7. Williams RW,
    8. Auwerx J
    (2013) Mitonuclear protein imbalance as a conserved longevity mechanism. Nature 497: 451–457. doi:10.1038/nature12188
    OpenUrlCrossRefPubMed
  22. ↵
    1. Hsin H,
    2. Kenyon C
    (1999) Signals from the reproductive system regulate the lifespan of C. elegans. Nature 399: 362–366. doi:10.1038/20694
    OpenUrlCrossRefPubMed
  23. ↵
    1. Huang da W,
    2. Sherman BT,
    3. Lempicki RA
    (2009a) Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37: 1–13. doi:10.1093/nar/gkn923
    OpenUrlCrossRefPubMed
  24. ↵
    1. Huang da W,
    2. Sherman BT,
    3. Lempicki RA
    (2009b) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57. doi:10.1038/nprot.2008.211
    OpenUrlCrossRefPubMed
  25. ↵
    1. Khan MH,
    2. Ligon M,
    3. Hussey LR,
    4. Hufnal B,
    5. Farber R 2nd.,
    6. Munkacsy E,
    7. Rodriguez A,
    8. Dillow A,
    9. Kahlig E,
    10. Rea SL
    (2013) TAF-4 is required for the life extension of isp-1, clk-1 and tpk-1 Mit mutants. Aging 5: 741–758. doi:10.18632/aging.100604
    OpenUrlCrossRef
  26. ↵
    1. Kim D,
    2. Langmead B,
    3. Salzberg SL
    (2015) HISAT: A fast spliced aligner with low memory requirements. Nat Methods 12: 357–360. doi:10.1038/nmeth.3317
    OpenUrlCrossRefPubMed
  27. ↵
    1. Korta DZ,
    2. Tuck S,
    3. Hubbard EJ
    (2012) S6K links cell fate, cell cycle and nutrient response in C. elegans germline stem/progenitor cells. Development (Cambridge, England) 139: 859–870. doi:10.1242/dev.074047
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Lapointe J,
    2. Stepanyan Z,
    3. Bigras E,
    4. Hekimi S
    (2009) Reversal of the mitochondrial phenotype and slow development of oxidative biomarkers of aging in long-lived Mclk1+/− mice. J Biol Chem 284: 20364–20374. doi:10.1074/jbc.m109.006569
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Larsen PL,
    2. Albert PS,
    3. Riddle DL
    (1995) Genes that regulate both development and longevity in Caenorhabditis elegans. Genetics 139: 1567–1583.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Lee LW,
    2. Lo HW,
    3. Lo SJ
    (2010) Vectors for co-expression of two genes in Caenorhabditis elegans. Gene 455: 16–21. doi:10.1016/j.gene.2010.02.001
    OpenUrlCrossRef
  31. ↵
    1. Liu JL,
    2. Yee C,
    3. Wang Y,
    4. Hekimi S
    (2017) A single biochemical activity underlies the pleiotropy of the aging-related protein CLK-1. Scientific Rep 7: 859. doi:10.1038/s41598-017-00754-z
    OpenUrlCrossRef
  32. ↵
    1. Liu X,
    2. Jiang N,
    3. Hughes B,
    4. Bigras E,
    5. Shoubridge E,
    6. Hekimi S
    (2005) Evolutionary conservation of the clk-1-dependent mechanism of longevity: Loss of mclk1 increases cellular fitness and lifespan in mice. Genes Dev 19: 2424–2434. doi:10.1101/gad.1352905
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Love MI,
    2. Huber W,
    3. Anders S
    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15: 550. doi:10.1186/s13059-014-0550-8
    OpenUrlCrossRefPubMed
  34. ↵
    1. Luo W,
    2. Brouwer C
    (2013) Pathview: An R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics (Oxford, England) 29: 1830–1831. doi:10.1093/bioinformatics/btt285
    OpenUrlCrossRefPubMed
  35. ↵
    1. Masoro EJ
    (2000) Caloric restriction and aging: An update. Exp Gerontol 35: 299–305. doi:10.1016/s0531-5565(00)00084-x
    OpenUrlCrossRefPubMed
  36. ↵
    1. Melendez A,
    2. Talloczy Z,
    3. Seaman M,
    4. Eskelinen EL,
    5. Hall DH,
    6. Levine B
    (2003) Autophagy genes are essential for dauer development and life-span extension in C. elegans. Science (New York, NY) 301: 1387–1391. doi:10.1126/science.1087782
    OpenUrlAbstract/FREE Full Text
  37. ↵
    1. Min KJ,
    2. Tatar M
    (2006) Restriction of amino acids extends lifespan in Drosophila melanogaster. Mech Ageing Dev 127: 643–646. doi:10.1016/j.mad.2006.02.005
    OpenUrlCrossRefPubMed
  38. ↵
    1. Monaghan RM,
    2. Barnes RG,
    3. Fisher K,
    4. Andreou T,
    5. Rooney N,
    6. Poulin GB,
    7. Whitmarsh AJ
    (2015) A nuclear role for the respiratory enzyme CLK-1 in regulating mitochondrial stress responses and longevity. Nat Cell Biol 17: 782–792. doi:10.1038/ncb3170
    OpenUrlCrossRefPubMed
  39. ↵
    1. Ogata H,
    2. Goto S,
    3. Sato K,
    4. Fujibuchi W,
    5. Bono H,
    6. Kanehisa M
    (1999) KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 27: 29–34. doi:10.1093/nar/27.1.29
    OpenUrlCrossRefPubMed
  40. ↵
    1. Pan KZ,
    2. Palter JE,
    3. Rogers AN,
    4. Olsen A,
    5. Chen D,
    6. Lithgow GJ,
    7. Kapahi P
    (2007) Inhibition of mRNA translation extends lifespan in Caenorhabditis elegans. Aging Cell 6: 111–119. doi:10.1111/j.1474-9726.2006.00266.x
    OpenUrlCrossRefPubMed
  41. ↵
    1. Quiros PM,
    2. Mottis A,
    3. Auwerx J
    (2016) Mitonuclear communication in homeostasis and stress. Nat Rev Mol Cell Biol 17: 213–226. doi:10.1038/nrm.2016.23
    OpenUrlCrossRefPubMed
  42. ↵
    1. Ruan L,
    2. Zhou C,
    3. Jin E,
    4. Kucharavy A,
    5. Zhang Y,
    6. Wen Z,
    7. Florens L,
    8. Li R
    (2017) Cytosolic proteostasis through importing of misfolded proteins into mitochondria. Nature 543: 443–446. doi:10.1038/nature21695
    OpenUrlCrossRefPubMed
  43. ↵
    1. Stefely JA,
    2. Pagliarini DJ
    (2017) Biochemistry of mitochondrial coenzyme Q biosynthesis. Trends Biochem Sci 42: 824–843. doi:10.1016/j.tibs.2017.06.008
    OpenUrlCrossRefPubMed
  44. ↵
    1. Stepanyan Z,
    2. Hughes B,
    3. Cliche DO,
    4. Camp D,
    5. Hekimi S
    (2006) Genetic and molecular characterization of CLK-1/mCLK1, a conserved determinant of the rate of aging. Exp Gerontol 41: 940–951. doi:10.1016/j.exger.2006.06.041
    OpenUrlCrossRefPubMed
  45. ↵
    1. Stout GJ,
    2. Stigter EC,
    3. Essers PB,
    4. Mulder KW,
    5. Kolkman A,
    6. Snijders DS,
    7. van den Broek NJ,
    8. Betist MC,
    9. Korswagen HC,
    10. Macinnes AW, et al.
    (2013) Insulin/IGF-1-mediated longevity is marked by reduced protein metabolism. Mol Syst Biol 9: 679. doi:10.1038/msb.2013.35
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Takahashi K,
    2. Noda Y,
    3. Ohsawa I,
    4. Shirasawa T,
    5. Takahashi M
    (2014) Extended lifespan, reduced body size and leg skeletal muscle mass, and decreased mitochondrial function in clk-1 transgenic mice. Exp Gerontol 58: 146–153. doi:10.1016/j.exger.2014.08.003
    OpenUrlCrossRef
  47. ↵
    1. Tiku V,
    2. Jain C,
    3. Raz Y,
    4. Nakamura S,
    5. Heestand B,
    6. Liu W,
    7. Spath M,
    8. Suchiman HED,
    9. Muller RU,
    10. Slagboom PE, et al.
    (2017) Small nucleoli are a cellular hallmark of longevity. Nat Commun 8: 16083. doi:10.1038/ncomms16083
    OpenUrlCrossRef
  48. ↵
    1. Tissenbaum HA,
    2. Ruvkun G
    (1998) An insulin-like signaling pathway affects both longevity and reproduction in Caenorhabditis elegans. Genetics 148: 703–717.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Walker AK,
    2. Rothman JH,
    3. Shi Y,
    4. Blackwell TK
    (2001) Distinct requirements for C. elegans TAF(II)s in early embryonic transcription. EMBO J 20: 5269–5279. doi:10.1093/emboj/20.18.5269
    OpenUrlAbstract/FREE Full Text
  50. ↵
    1. Wang J,
    2. Robida-Stubbs S,
    3. Tullet JM,
    4. Rual JF,
    5. Vidal M,
    6. Blackwell TK
    (2010) RNAi screening implicates a SKN-1-dependent transcriptional response in stress resistance and longevity deriving from translation inhibition. PLoS Genet 6: e1001048. doi:10.1371/journal.pgen.1001048
    OpenUrlCrossRefPubMed
  51. ↵
    1. Wang X,
    2. Zuo X,
    3. Kucejova B,
    4. Chen XJ
    (2008) Reduced cytosolic protein synthesis suppresses mitochondrial degeneration. Nat Cell Biol 10: 1090–1097. doi:10.1038/ncb1769
    OpenUrlCrossRefPubMed
  52. ↵
    1. Yoon HJ,
    2. Donahue TF
    (1992) The suil suppressor locus in Saccharomyces cerevisiae encodes a translation factor that functions during tRNA(iMet) recognition of the start codon. Mol Cell Biol 12: 248–260. doi:10.1128/mcb.12.1.248
    OpenUrlAbstract/FREE Full Text
  53. ↵
    1. Zid BM,
    2. Rogers AN,
    3. Katewa SD,
    4. Vargas MA,
    5. Kolipinski MC,
    6. Lu TA,
    7. Benzer S,
    8. Kapahi P
    (2009) 4E-BP extends lifespan upon dietary restriction by enhancing mitochondrial activity in Drosophila. Cell 139: 149–160. doi:10.1016/j.cell.2009.07.034
    OpenUrlCrossRefPubMed
Back to top
Download PDF
Email Article
Citation Tools
Clk-1–dependent ubiquinone biosynthesis affects translation
Marte Molenaars, Georges E Janssens, Toon Santermans, Marco Lezzerini, Rob Jelier, Alyson W MacInnes, Riekelt H Houtkooper
Life Science Alliance Aug 2018, 1 (5) e201800082; DOI: 10.26508/lsa.201800082
Twitter logo Facebook logo Mendeley logo
Issue Cover

In this Issue

Volume 1, No. 5
October 2018
Advertisement

Jump to section

Subjects