- Research article
- Open Access
A metabolic signature of long life in Caenorhabditis elegans
BMC Biologyvolume 8, Article number: 14 (2010)
Many Caenorhabditis elegans mutations increase longevity and much evidence suggests that they do so at least partly via changes in metabolism. However, up until now there has been no systematic investigation of how the metabolic networks of long-lived mutants differ from those of normal worms. Metabolomic technologies, that permit the analysis of many untargeted metabolites in parallel, now make this possible. Here we use one of these, 1H nuclear magnetic resonance spectroscopy, to investigate what makes long-lived worms metabolically distinctive.
We examined three classes of long-lived worms: dauer larvae, adult Insulin/IGF-1 signalling (IIS)-defective mutants, and a translation-defective mutant. Surprisingly, these ostensibly different long-lived worms share a common metabolic signature, dominated by shifts in carbohydrate and amino acid metabolism. In addition the dauer larvae, uniquely, had elevated levels of modified amino acids (hydroxyproline and phosphoserine). We interrogated existing gene expression data in order to integrate functional (metabolite-level) changes with transcriptional changes at a pathway level.
The observed metabolic responses could be explained to a large degree by upregulation of gluconeogenesis and the glyoxylate shunt as well as changes in amino acid catabolism. These responses point to new possible mechanisms of longevity assurance in worms. The metabolic changes observed in dauer larvae can be explained by the existence of high levels of autophagy leading to recycling of cellular components.
See associated minireview: http://jbiol.com/content/9/1/7
The nematode Caenorhabditis elegans normally has a life-span of about three weeks. The dauer larva, however, lives for up to eight times longer . In addition, mutations in scores of genes have been identified that increase longevity. These genes have been grouped into several pathways including the Insulin/Insulin-Like signalling pathway (IIS) [2–6], the dietary restriction pathway [7, 8] and the translation control pathway , but how they regulate ageing individually and together is still obscure. What is certain, however, is that each of them influences the metabolism of the worm in some fashion. This has been shown by the discovery that particular longevity pathways control, or at least interact with, key regulators of metabolism [10–13] as well as many metabolic enzymes [14–22].
Despite these advances, our understanding of how altered metabolism influences longevity in worms, indeed, if it does so at all, remains very incomplete. In part, this is because attention has focused almost exclusively on the genes that control metabolism rather than metabolites themselves. Yet gene activity can only give a very dim outline of the activity of a metabolic network since much regulation occurs at the post-transcriptional, or even post-translational level, for example, by allosteric interactions among metabolites and the enzymes that catalyse them [23, 24]. One way to investigate the activity of metabolic networks in a more direct fashion is metabolite profiling (also sometimes called metabolomics or metabonomics). Metabolomics has been previously combined with functional genomics to study a variety of biological problems and species [25–28] including, recently, C. elegans [29–31]. Here, we apply it to investigating the metabolic networks of a series of worms that are, for one reason or another, long-lived. Most of our long-lived worms are defective for components of the IIS pathway and one is translation defective; but we also study the dauer stage that forms when larvae are grown under stressful conditions. We show that all these long-lived worms have metabolic profiles that are not only very different from normal worms but also very similar to each other; in other words, that there is a metabolic signature for long-life in worms. The existence of this signature is surprising since the IIS and translation pathways are, at least, thought to influence longevity by quite distinct mechanisms [9, 12]. This signature is composed of metabolites that function in several distinct parts of the network, including carbohydrate, amino acid and choline metabolism. Since our ultimate goal is an integrated model of worm metabolism, we also interrogate existing global gene expression data from daf-2 mutant worms  to give a general account of how the metabolic networks of long-lived worms differ from those with normal life-spans.
Results and discussion
Long-lived mutants have distinctive metabolic profiles
Of the various pathways known to regulate longevity in worms, the best known is the Insulin/Insulin-Like signalling (IIS) pathway [32, 33]. Many mutations that disrupt components of this pathway affect the ability of larval worms to enter and leave the dauer stage, but they also increase the longevity and stress resistance of adults as well as reduce their fecundity [2–6, 34].
We began by studying m41, a hypomorph mutation that disrupts daf-2 which encodes a tyrosine kinase that is expressed throughout the worm and is thought to act as a receptor for many of the 37 insulin-like ligands present in the C. elegans genome [35, 36]. daf-2(m41) hermaphrodites are 10 to 90% longer-lived than wild-type worms [3, 17, 37, 38] (our data not shown). Since m41 is a dauer-constitutive temperature-sensitive mutation we grew these worms at the permissive temperature, 15°C, until L4, transferred them to 22.5°C, and assayed their metabolites as old adults (240 hours). We did this by freezing the worms instantly in liquid nitrogen, extracting polar metabolites, and then acquiring 1H NMR spectra. The spectra showed a range of resonances from small molecule metabolites, typical of tissue extracts. We then divided the spectra into bins each chosen to represent as far as possible a single metabolite resonance. Principal components analysis (PCA) and hierarchical cluster analysis (HCA) of the reduced data showed that daf-2(m41) and wild-type samples have distinct metabolic profiles with little overlap between the two groups of samples in PC1 (Figure 1A). The loadings along this axis showed that many NMR-detectable metabolites contribute to the difference between the genotypes (Figure 1B).
Many IIS mutations exist and they differ in the severity and kind of their phenotypic effects [3, 34]. So, in a separate experiment we simultaneously studied three daf-2 mutations: m41, e1370 and m596 as well as daf-28(sa191) which disrupts an insulin-like ligand thought to bind DAF-2. DAF-28 is thought to activate DAF-2 and so promote normal, reproductive growth and longevity, but sa191 is a dominant negative gain-of-function allele . Like daf-2 hypomorphs, daf-28(sa191) is partly dauer constitutive, has long-lived adults, and can be repressed by mutations in daf-16 . In this experiment, we used fewer samples of worms per genotype than in previous experiments, but sampled three ages, L1 (15 hours), middle-age (144 hours) and old age (240 hours), and raised them first at 15°C and then at 25°C. Considering just the old-age samples we found that all of these mutants have distinctive metabolic profiles, with e1370 and m41 having respectively the least and most distinctive metabolic phenotypes relative to wild-type (Figure 2A). The rank order of allele strength for longevity at 22.5°C and 25°C is e1370 <m41 <m596 (our data not shown; ), but daf-2 alleles have a variety of phenotypic effects which do not all show the same rank order of severity [3, 34, 39]. The three daf-2 mutations affect different parts of the receptor: e1370 disrupts the intracellular tyrosine kinase domain while m41 and m596 disrupt, respectively, the Cysteine Rich and Leucine Rich L2 extracellular domains [34, 39]. Some L2 domain mutations in the human Insulin receptors have very low ligand binding affinity ; a similar property of daf-2(m596) may explain why its metabolome resembles that of daf-28(sa191) which disrupts a putative ligand.
When surveying the four daf-2(-) mutants, we also looked at another kind of long-lived mutant, ife-2(ok306), which disrupts a gene encoding an isoform of the eukaryotic translation initiation factor, eIF4E . Since this mutation does not require DAF-16 to confer increased longevity, it is thought that IFE-2 works either downstream or in parallel to DAF-2 to regulate longevity. We found that the metabolic profile of ife-2(ok306) is very similar to that of the ILS mutants: cluster analysis and PCA do not clearly separate ife-2(ok306) worms from ILS mutants (Figure 2A).
Since we sampled all of our mutants at three ages we were also able to study, at least crudely, when the mutant worms acquired their distinctive metabolic profiles. PCA and cluster analysis shows that all mutants had distinctive metabolic profiles even as larvae, but in all cases the metabolic profiles became increasingly different from wild-type with age (Figure 2B). Since reproduction has not yet begun in L1 larvae, which do not even have gonads, the distinctive profiles of the long-lived mutants cannot be entirely due to a decrease in metabolic resources allocated to reproduction.
Finally, in this same experiment, we also studied dauer larvae. Dauers form when L2 worms are crowded or deprived of food. They do not feed, have very distinctive transcriptional profiles, physiologies and morphologies, are very stress resistant and also do not age . We raised dauers at two temperatures, 20°C (n = 3) and 25°C (n = 5), and compared them to old adults (240 hours) raised at the same temperatures (n = 5 for both). Clustering and PCA showed that dauers and adults have unambiguously distinct metabolic profiles as do worms raised at different temperatures, with temperature nested within the two stages (Figure 2C). Comparing dauers to L1s or young adults gave very similar results (data not shown). We found that metabolite levels showed strong stage × temperature interactions. This is reflected in the reversal of the relative positions of dauers and adults along the PC 2 axis depending on the temperature and also in analysis of variance on individual bins (data not shown).
The metabolic signature of long life in worms
What are the metabolic features of long-life? In C. elegans, many mutants and environmental treatments confer increased longevity, but the devices by which they do so, or whether or not they are the same, remains unclear. One reason for this is that few studies assay more than one class of long-lived worm using the same technology and experimental conditions. Since we have studied three classes of long-lived worm, dauers, IIS defective and a translation-defective mutant, using the same metabolic profiling technique, we can directly compare them and ask what, if anything, do they have in common?
To do this we determined the relative concentrations of 26 metabolites (Figure 3A; Additional file 1). Then, to identify the metabolic signatures of long-life we ranked them by consistency and direction of response in long-lived worms relative to wild-type (Figure 3B). Although this figure shows the results of all of our experiments, here we focus on the single experiment in which worms were raised at 25°C and sampled at 10 days after hatching.
We found that the metabolic responses of our long-lived worms were strikingly similar. More than half of the examined metabolites show qualitatively similar changes in dauers, IIS mutants, and ife-2 mutant worms. This result was surprising since IIS and ife-2 mutations ostensibly influence very different aspects of the worm's physiology. We propose that these metabolites constitute a minimal metabolic signature of long-life in worms.
One of the signature metabolites was the disaccharide trehalose. An important carbohydrate storage molecule in nematodes, trehalose is thought to confer stress resistance in many invertebrates [41–43]. Previous studies have shown that the expression of genes involved in its synthesis are elevated in dauers and IIS mutants so our finding that trehalose pool sizes are elevated in dauers and IIS mutants was expected; indeed, trehalose has been proposed as a longevity assurance sugar [19, 44]. ife-2(ok306) worms show high trehalose levels as well implying that a deficiency in protein synthesis can affect carbohydrate metabolism as well.
This similarity across different classes of long-lived worms is also seen in amino acid levels. Of the 12 amino acids we studied, 11 are regulated in IIS mutants. Of these, 10 are regulated in the same way in ife-2(oK306) and 5 are in dauers. We also found the following metabolites consistently regulated across dauer, IIS, and ife-2(0) mutants: choline, phosphocholine, and glycerophosphocholine (GPC), which are associated with lipid metabolism; acetate, malate and succinate, which are associated with carbohydrate metabolism; propanoate and NAD+. Some of these longevity-signature metabolites are expected from previous studies of long-lived worms. This is particularly true of metabolites that have a role in carbohydrate metabolism and we consider them in greater detail below. Others, such as altered pool sizes of amino acids, choline, phosphocholine, propanoate and NAD+, could point to new mechanisms of longevity assurance in C. elegans. We also note that, although we have ignored metabolic responses peculiar to particular worm strains, we cannot exclude the possibility that they influence longevity as well.
Autophagy and the dauer metabolome
Although the dauer metabolome resembles that of long-lived mutants in many ways, we also found that it has some unique features. We detected pools of two post-translationally modified amino acids, phosphoserine and hydroxyproline, in dauers and only dauers. Since hydroxylation and phosphorylation generally take place on peptides rather than free amino acids, these pools are likely the result of protein degradation. The most obvious source of free hydroxyproline is collagen: the C. elegans genome contains approximately 175 genes encoding collagens which are used in the basement membranes and cuticle, two prolyl 4-hydroxylases, and at least five peptidases that are required for the processing or turn-over of cuticle collagens [45, 46]. The source of free phosphoserine is less obvious, however, given that we studied high concentration metabolites, they probably are not derived from the phosphorylated serines found in signalling pathways but rather represent structural components. We suspect that the phosphoserine pool seen in dauers is derived from serine phosphoglyceride lipids freed by turnover of membrane phospholipids. This is supported by the observation that choline compounds (choline, phosphocholine, and GPC) are also strongly increased in dauers, and altered choline compound concentrations are frequently observed in mammalian tumours where they mark the membrane turnover characteristic of rapidly proliferating cells . Since dauers do not feed they rely on energy stores such as fats and glycogen to survive, and much evidence shows that the beta-oxidation and glyoxylate pathways that metabolize fatty acids are upregulated in dauers [19, 48–50]. One explanation, then, for the elevated levels of these modified amino acids is that dauers are utilizing spare extracellular matrix, and other proteins, as yet another energy store. Consistent with this idea, dauer morphogenesis requires extensive autophagy . A mutation that abolishes autophagy also decreases longevity specifically in daf-2(e1370) which implies that high levels of autophagy promote long-life ; if so we did not detect any obvious signatures of this.
Some, but not all, signature metabolite responses require DAF-16
Much evidence shows that the longevity prolonging effects of IIS mutants are mediated by the FOXO transcription factor, DAF-16 [2, 52, 53]. Down-regulation of DAF-2 signalling results in nuclear localization and hence activation of DAF-16 which, in turn, activates or represses many genes which contribute to longevity [14, 54, 55]. Since the longevity prolonging effect of daf-2 mutants is repressed by null mutations of daf-16, one way of disentangling the phenotypes of IIS mutants that contribute causally to longevity from those that do not, is to ask whether they, too, depend on DAF-16 activity. The reasoning, previously applied to transcriptomic and proteomic data, is that any molecular phenotype that contributes to longevity should be abolished by inhibiting DAF-16 activity just as longevity itself [15, 21].
To find out whether IIS control over metabolism was also DAF-16 dependent we compared the metabolomes of 144 hour-old wild-type worms to worms carrying either another daf-2 allele, e1370, or the null daf-16(m26) mutation, or both. Unsupervised methods (PCA and HCA) divide these samples into two major groups: a cluster which contains 7/8 daf-2(-) samples and a cluster that contains the rest (Figure 4A). Thus the metabolic phenotype of daf-2 is at least partially DAF-16 dependent; however the double mutants form a sub-cluster clearly distinct from wild-type implying that not all of the daf-2 phenotype is so.
By measuring metabolite levels directly in single mutant daf-2(e1370) and daf-16(m26) and double mutant daf-2(e1370);daf-16(m26) worms we were able to apply this test to 11 metabolites (Figure 4B). Of these only four, phosphocholine, and the closely related amino acids isoleucine, valine, and leucine, showed the classic pattern of DAF-16 dependence: concentrations of each were substantially elevated or repressed in daf-2(-), but not daf-16(0) or double mutant samples. Several other metabolites showed more complex patterns of epistasis. For example, lysine, lactate and glycerol concentrations were all significantly different in either daf-16(0) or daf-2(-);daf-16(0) worms or both compared to wild-type worms (P = 0.05; two-tailed t-test) and so fail the classical test of DAF-16 dependence. Surprisingly, trehalose showed no sign of DAF-16 dependence: daf-16(0) samples have wild-type levels of the sugar, but double mutants are not different from daf-2(-). This is in contrast to previous results shown by transcriptomic studies [15, 19], which may perhaps reflect the fact that changes in gene expression levels alone do not necessarily equate to functional differences [56, 57].
Classical DAF-16 dependence make isoleucine, valine, leucine and phosphocholine strong candidates for having a causal role in long life, particularly as all four are signature metabolites as defined above. Conversely, the absence of DAF-16 dependence in trehalose suggests that DAF-2 regulates it via another transcription factor parallel to DAF-16 and that it may not contribute to long life. The interpretation of non-classical daf-16 epistasis is less clear. The classical test supposes that DAF-16 is fully repressed in normal worms, but activated in the absence of DAF-2 signalling by translocation of the transcription factor from cytoplasm to nuclei [55, 58]. This model is certainly too simple since daf-16(0) mutants have a variety of subtle phenotypes such as rapid growth, early reproduction and a slightly reduced lifespan, and normal worms have at least some DAF-16 visible in their nuclei [59, 60]. Some metabolites are, then, also apparently sensitive to low levels of DAF-16 activity.
Metabolic targets of DAF-2 signalling
We have shown that the pool sizes of many metabolites differ between long-lived and normal worms and that some of these differences are shared by various long-lived mutants, in particular the several daf-2 alleles that we studied. But what genes does daf-2 regulate that result in these changes? In order to investigate this, we mapped some of our signature metabolites onto a standard metazoan metabolic network so that we could identify those parts of the network that are altered in daf-2 worms. We then asked whether the genes encoding metabolic enzymes that work in the same parts of the network are regulated as well, and if so, whether the two sources of data could be used to give an economical account of how the metabolism of C. elegans is altered in daf-2 worms. As an initial guide we used NEMAPATH  to identify, for each metabolite, the pathways in which they might work and the C. elegans genes that might act in them. We then interrogated a previously published global expression dataset based on daf-2(e1370) worms raised under conditions similar to ours  and examined the expression patterns of these genes for patterns of co-regulation.
We found that at least five of the metabolites regulated in daf-2(-) worms, malate, acetate, succinate, glucose and trehalose, were linked by three major pathways: the glyoxylate shunt, gluconeogenesis, and starch metabolism. Consistent with previous studies [15, 19, 21] we found that genes encoding enzymes in these pathways are up-regulated in daf-2(-) worms. In all, the C. elegans genome contains 38 genes encoding components of 20 enzymes that work in these pathways. Of these 38 genes, 14 are significantly up-regulated, 5 are significantly down-regulated, 18 are not regulated, and 1 has no data (Figure 5, Additional file 2).
While the observed shifts in carbohydrate metabolism might have been expected from the results of gene expression studies, the changes in amino acid metabolism shown by our data were not. Most amino acid pools are upregulated in long-lived worms (8/12 in the IIS mutants; 10/12 in ife-2). One possible explanation for this is that protein synthesis is generally repressed in IIS mutants as it is in ife-2  and that the amino acid pool represents a surplus. daf-2(e1370) worms may have reduced protein synthesis since the expression of their t-RNA synthetases are generally repressed (Additional file 2). However, there is no correlation between amino acid pool size and tRNA synthetase expression (data not shown). It is more likely then that amino acid pool sizes are dictated by catabolic pathways that direct them to energy production or other uses. Consistent with this, many genes that encode components of phenylalanine and tyrosine catabolic pathways, among them tyrosinase and phenylalanine hydroxylase, are regulated in daf-2 (Additional file 2) [21, 62]. Intriguingly, melanin, a tyrosine metabolite, has recently been discovered in worm cuticles where it is thought to have a protective function .
The most striking response among the amino acids, however, is the upregulation of the branched chain amino-acids (BCAAs) isoleucine, leucine and valine. The pool sizes of these amino acids are positively correlated across long-lived mutants (Figure 3A). Furthermore, unlike most other metabolites, their upregulation in daf-2 is entirely DAF-16 dependent, making them strong candidates for being causally involved in longevity (Figure 4B). Like other animals, C. elegans cannot synthesize these amino acids , and so any difference in their relative concentrations must be due to a change in either protein turnover or their catabolism. In fact, BCAA pool sizes are co-regulated in many circumstances such as growth in worms  or obesity in humans [65, 66]. This co-regulation is a consequence of them sharing the first two steps in their catabolic pathways: transamination by BCAT and oxidative carboxylation by the mitochondrial BCKD enzyme complex . In daf-2(e1370) worms, BCAT expression is wild-type, but all four genes encoding components of the BCKD complex are strongly downregulated (Additional file 2). We hypothesize that downregulation of the BCKD complex is responsible for the increased BCAA pool sizes of daf-2 worms. This hypothesis also suggests a way to manipulate BCAA pool sizes to test their contribution to long life. Strong inactivation of BCKD-complex genes in worms causes severe embryonic and larval phenotypes (; wormbase.org) and, in humans, maple syrup urine disease, a metabolic disorder resulting in encephalopathy and death ; however, it remains possible that more subtle elevation of BCAA levels by diet or partial downregulation of the BCKD complex will confer long life.
By studying the metabolic profiles of a range of long-lived worms we have identified a metabolic signature of long life common to dauers, IIS mutants and a translation defective mutant. Some of the metabolites that comprise this signature, such as those involved in carbohydrate metabolism, are expected from studies of global gene expression; others, such as those involving amino acid metabolism, are new. The existence of a common metabolic signature for long life suggests that longevity pathways that have been previously considered independent may, in fact, regulate the same regions of the metabolic network. By interrogating an existing global gene expression profile dataset on daf-2 worms, we have identified some of those regions. We find that the changes in carbohydrate metabolism can be explained by upregulation of the glyoxylate shunt and gluconeogenesis. We also find that the general elevation of amino acid pool sizes in long-lived worms is likely due to regulation of catabolic pathways that divert amino acids away from protein synthesis and to energy metabolism or other survival functions. Dauers show many of the features of long-lived mutants, but they also show elevated levels of phosphoserine, hydroxyproline, and choline compounds; we suggest that these are probably the consequence of extensive autophagy.
Although we have shown how the study of metabolite concentrations and gene expression can be used together to give a consistent account of the metabolism of long-lived worms, we recognize that we have considered only a few of the many metabolic differences that may give these worms their remarkable life-spans. Ultimately, we would like to provide a general metabolic model for long-life - one that quantitatively explains longevity in terms of the expression of genes that encode metabolic enzymes, the activity of those enzymes, and the metabolic concentrations and fluxes that result from them. That goal however, requires a well-founded native metabolic network for C. elegans as well as a general account of its kinetics and how they are regulated, none of which currently exist. Nevertheless, the recent development of system-wide metabolic models, particularly of microbes, suggests that such a global model of C. elegans metabolism may be possible in the near future [69–72]. Our study begins to provide the empirical basis upon which it will depend.
We used the following strains: wild-type N2 Bristol, DR26 daf-16(m26), CB1370 daf-2(e1370), DR1564 daf-2(m41), DR1565 daf-2(m596), DR1309 daf-16;daf-2(m26;e1370), JT191 daf-28(sa191), KX15 ife-2(ok306) and RB579 ife-2(ok306), the outcrossed version of KX15. They were provided by the Caenorhabditis elegans Center (CGC) at the University of Minnesota.
Worms were grown using standard techniques. daf-2(-) alleles are temperature-sensitive dauer constitutive. Synchronized populations were grown at the permissive temperature, 20°C, until L4 and then transferred to 22.5°C or 25°C. Controls and non-ts strains were treated in the same way. Dauers were induced by growing worms as above but at high density. Worms were sampled for NMR at either 15 h (L1), 144 h (adult, six days), or 240 h (adult, 10 days) post hatch. To ensure that only old worms were sampled at 10 days, their offspring were removed by repeated filtering using 50 μm Nitex membranes (Sefar Ltd., Bury, UK).
Nuclear magnetic resonance (NMR) spectroscopy
For NMR analysis worms were washed off plates and collected in 0.5 ml of M9 buffer , quick-frozen in liquid nitrogen and stored at -80°C. We then ground the tissue at liquid nitrogen temperatures in a mortar and pestle, and added 2 ml ice-cold methanol directly to the ground tissue to give a final concentration of 80% methanol. We transferred the extract to microcentrifuge tubes, and rinsed the mortars with an additional 2 ml of 80% methanol. We then centrifuged the extracts (10 minutes, 16,000 g) and dried the supernatants in a rotary vacuum concentrator. We rehydrated each sample in 650 μl of NMR buffer (100% 2H2O, 0.97 mM sodium trimethylsilyl-2,2,3,3-2H4-propionate (TSP), 0.1 M phosphate buffer pH 7.0), centrifuged again to remove any particulate matter, and transferred 600 μl to 5 mm NMR tubes.
NMR spectra were acquired essentially as described by Beckonert et al.  using a Bruker Avance DRX600 spectrometer (Bruker BioSpin, Rheinstetten, Germany) with a field strength of 14.1 T and consequent 600 MHz 1H resonance frequency, equipped with a 5 mm cryogenically-cooled inverse geometry probe. A 1D NOESY pulse sequence was used for water suppression, with an acquisition time of 1.36 s, and an additional relaxation delay of 3.5 s, with presaturation during the relaxation recovery and 0.1 s mixing time giving a 5 s recycle time; we collected 160 transients per sample, following four dummy scans to allow the system to approach a steady state. The data were acquired into 32 K points over a 12 kHz spectral width.
NMR data processing and analysis
We carried out initial data processing in iNMR v.2. The summed transients were multiplied by an exponential apodization function equivalent to 0.5 Hz line broadening and zero-filled by 50%, followed by Fourier transformation. The spectra were referenced to the TSP resonance at 0 ppm, and phase correction and first-order baseline correction carried out using the software's proprietary algorithms. We visually identified peaks in the spectra and divided them manually into bins (integrals); compared to equal-interval binning of the entire spectrum, this has the effect of reducing the total number of variables, aligning each bin more closely with an individual resonance, and excluding spectral regions that contain only noise across all samples. Around 40 detectable metabolites could be readily identified in routine 1D spectra of the worm extracts (Additional file 3). There were also a number of resonances from as-yet unassigned metabolites (for example, singlets at δ 1.993 and δ 1.986 ppm represent probable N-acetyl groups). In addition, we re-processed all spectra in Chenomx NMR Suite 4.6 (Chenomx, Edmonton, AB, Canada) and quantified metabolite concentrations for selected metabolites by computer-assisted manual fitting of metabolites. This software fits idealized spectra made up of combinations of Lorentzian peaks, based on authentic standards . We assigned metabolite resonances by comparing their multiplicity and chemical shift to compounds found in the Chenomx database. This was supplemented by 2D NMR experiments (COSY and HSQC) acquired for typical samples, and additional comparisons to our own in-house standards data and other online databases. All of the metabolites fitted were present in the Chenomx proprietary database, except for trehalose, which we added to the database.
We then normalized the data by dividing each profile by a single normalization factor, the median fold change across all compounds relative to a reference profile (a median of all profiles), as described by Dieterle et al. , and log-transformed them by log10(ni + x). The transformation reduced the dominating effect of the high intensities of a few metabolites, so that intensity of otherwise weak and insignificant peaks is increased. The constant x was chosen such that the dependence between standard deviation and intensity was removed for a series of technical replicates, that is, increasing homoscedasticity (the principle is discussed by Purohit et al. , although a different transform was used by these authors). Multivariate analyses (principal components analysis, PCA, and hierarchical cluster analysis, HCA) were performed in either Aabel (Gigawiz, Tulsa, OK, USA) or JMP (SAS UK, Marlow, UK) as appropriate. Using the fitted data (that is, from Chenomx), we then tested the level of each metabolite in each experiment against its appropriate control (Additional file 1). In order to identify those metabolites that are generally regulated across multiple longevity treatments, we used a Fisher's combined probability test (Additional file 1, Part 4). As an additional test of the overall consistency in the pattern of regulation, we applied a sample randomization test and found as a whole the data set is highly structured (P < 1e-10).
Gene expression analysis
We downloaded Shaw et al.'s  gene expression data from PUMAdb http://puma.princeton.edu/, and interrogated it for the expression of genes putatively encoding all metabolic enzymes. In these experiments, daf-2(e1370) worms were grown at 25°C and so the results are also relevant to our data. Gene identities and metabolic pathways were based on NemaPath , a version of KEGG [78–80] native for nematode sequences http://www.nematode.net/cgi-bin/keggview.cgi.
branched chain amino acid
branched chain ketoacid dehydrogenase
hierarchical cluster analysis
heteronuclear single quantum coherence
nuclear magnetic resonance
nuclear Overhauser effect spectroscopy
principal component analysis
Klass M, Hirsh D: Non-ageing developmental variant of Caenorhabditis elegans. Nature. 1976, 260: 523-525. 10.1038/260523a0.
Kenyon C, Chang J, Gensch E, Rudner A, Tabtiang R: A C. elegans mutant that lives twice as long as wild type. Nature. 1993, 366: 461-464. 10.1038/366461a0.
Gems D, Sutton AJ, Sundermeyer ML, Albert PS, King KV, Edgley ML, Larsen PL, Riddle DL: Two pleiotropic classes of daf-2 mutation affect larval arrest, adult behavior, reproduction and longevity in Caenorhabditis elegans. Genetics. 1998, 150: 129-155.
Morris JZ, Tissenbaum HA, Ruvkun G: A phosphatidylinositol-3-OH kinase family member regulating longevity and diapause in Caenorhabditis elegans. Nature. 1996, 382: 536-539. 10.1038/382536a0.
Paradis S, Ailion M, Toker A, Thomas JH, Ruvkun G: A PDK1 homolog is necessary and sufficient to transduce AGE-1 PI3 kinase signals that regulate diapause in Caenorhabditis elegans. Genes Dev. 1999, 13: 1438-1452. 10.1101/gad.13.11.1438.
Paradis S, Ruvkun G: Caenorhabditis elegans Akt/PKB transduces insulin receptor-like signals from AGE-1 PI3 kinase to the DAF-16 transcription factor. Genes Dev. 1998, 12: 2488-2498. 10.1101/gad.12.16.2488.
Lakowski B, Hekimi S: The genetics of caloric restriction in Caenorhabditis elegans. Proc Natl Acad Sci USA. 1998, 95: 13091-13096. 10.1073/pnas.95.22.13091.
Panowski SH, Wolff S, Aguilaniu H, Durieux J, Dillin A: PHA-4/Foxa mediates diet-restriction-induced longevity of C-elegans. Nature. 2007, 447: 550-10.1038/nature05837.
Syntichaki P, Troulinaki K, Tavernarakis N: eIF4E function in somatic cells modulates ageing in Caenorhabditis elegans. Nature. 2007, 445: 922-926. 10.1038/nature05603.
Apfeld J, O'Connor G, McDonagh T, DiStefano PS, Curtis R: The AMP-activated protein kinase AAK-2 links energy levels and insulin-like signals to lifespan in C-elegans. Genes Dev. 2004, 18: 3004-3009. 10.1101/gad.1255404.
Jia K, Chen D, Riddle DL: The TOR pathway interacts with the insulin signaling pathway to regulate C. elegans larval development, metabolism and life span. Development. 2004, 131: 3897-3906. 10.1242/dev.01255.
Sheaffer KL, Updike DL, Mango SE: The Target of Rapamycin pathway antagonizes pha-4/FoxA to control development and aging. Curr Biol. 2008, 18: 1355-1364. 10.1016/j.cub.2008.07.097.
Wang MC, O'Rourke EJ, Ruvkun G: Fat Metabolism Links Germline Stem Cells and Longevity in C. elegans. Science. 2008, 322: 957-960. 10.1126/science.1162011.
Lee SS, Kennedy S, Tolonen AC, Ruvkun G: DAF-16 target genes that control C-elegans life-span and metabolism. Science. 2003, 300: 644-647. 10.1126/science.1083614.
Murphy CT, McCarroll SA, Bargmann CI, Fraser A, Kamath RS, Ahringer J, Li H, Kenyon C: Genes that act downstream of DAF-16 to influence the lifespan of Caenorhabditis elegans. Nature. 2003, 424: 277-284. 10.1038/nature01789.
McElwee J, Bubb K, Thomas JH: Transcriptional outputs of the Caenorhabditis elegans forkhead protein DAF-16 (vol 2, pg 111, 2003). Aging Cell. 2003, 2: 341-341. 10.1046/j.1474-9728.2003.00043.x.
Halaschek-Wiener J, Khattra JS, McKay S, Pouzyrev A, Stott JM, Yang GS, Holt RA, Jones SJM, Marra MA, Brooks-Wilson AR, Riddle DL: Analysis of long-lived C. elegans daf-2 mutants using serial analysis of gene expression. Genome Res. 2005, 15: 603-615. 10.1101/gr.3274805.
Hamilton B, Doug YQ, Shindo M, Liu WY, Odell I, Ruvkun G, Lee SS: A systematic RNAi screen for longevity genes in C-elegans. Genes Dev. 2005, 19: 1544-1555. 10.1101/gad.1308205.
McElwee JJ, Schuster E, Blanc E, Thornton J, Gems D: Diapause-associated metabolic traits reiterated in long-lived daf-2 mutants in the nematode Caenorhabditis elegans. Mech Ageing Dev. 2006, 127: 458-472. 10.1016/j.mad.2006.10.002.
Shaw WM, Luo S, Landis J, Ashraf J, Murphy CT: The C-elegans TGF-beta dauer pathway regulates longevity via insulin signaling. Curr Biol. 2007, 17: 1635-1645. 10.1016/j.cub.2007.08.058.
Dong MQ, Venable JD, Au N, Xu T, Park SK, Cociorva D, Johnson JR, Dillin A, Yates JR: Quantitative mass spectrometry identifies insulin signaling targets in C-elegans. Science. 2007, 317: 660-663. 10.1126/science.1139952.
Hansen M, Hsu AL, Dillin A, Kenyon C: New genes tied to endocrine, metabolic, and dietary regulation of lifespan from a Caenorhabditis elegans genomic RNAi screen. PLoS Genet. 2005, 1: 119-128. 10.1371/journal.pgen.0010017.
Kell DB: Metabolomics and systems biology: making sense of the soup. Curr Opin Microbiol. 2004, 7: 296-307. 10.1016/j.mib.2004.04.012.
ter Kuile BH, Westerhoff HV: Transcriptome meets metabolome: hierarchical and metabolic regulation of the glycolytic pathway. FEBS Lett. 2001, 500: 169-171. 10.1016/S0014-5793(01)02613-8.
Fiehn O, Kopka J, Dormann P, Altmann T, Trethewey RN, Willmitzer L: Metabolite profiling for plant functional genomics (vol 187, pg 2000). Nature Biotechnology. 1157, 19: 173-173.
Raamsdonk LM, Teusink B, Broadhurst D, Zhang NS, Hayes A, Walsh MC, Berden JA, Brindle KM, Kell DB, Rowland JJ, Westerhoff HV, van Dam K, Oliver SG: A functional genomics strategy that uses metabolome data to reveal the phenotype of silent mutations. Nat Biotechnol. 2001, 19: 45-50. 10.1038/83496.
Castrillo JI, Zeef LA, Hoyle DC, Zhang N, Hayes A, Gardner DC, Cornell MJ, Petty J, Hakes L, Wardleworth L, Rash B, Brown M, Dunn WB, Broadhurst D, O'Donoghue K, Hester SS, Dunkley TP, Hart SR, Swainston N, Li P, Gaskell SJ, Paton NW, Lilley KS, Kell DB, Oliver SG: Growth control of the eukaryote cell: a systems biology study in yeast. J Biol. 2007, 6: 4-10.1186/jbiol54.
Vijayendran C, Barsch A, Friehs K, Niehaus K, Becker A, Flaschel E: Perceiving molecular evolution processes in Escherichia coli by comprehensive metabolite and gene expression profiling. Genome Biol. 2008, 9: R72-10.1186/gb-2008-9-4-r72.
Blaise BJ, Giacomotto J, Elena B, Dumas ME, Toulhoat P, Segalat L, Emsley L: Metabotyping of Caenorhabditis elegans reveals latent phenotypes. Proc Natl Acad Sci USA. 2007, 104: 19808-19812. 10.1073/pnas.0707393104.
Atherton HJ, Jones OA, Malik S, Miska EA, Griffin JL: A comparative metabolomic study of NHR-49 in Caenorhabditis elegans and PPAR-alpha in the mouse. FEBS Lett. 2008, 582: 1661-1666. 10.1016/j.febslet.2008.04.020.
Swire J, Fuchs S, Bundy JG, Leroi AM: The cellular geometry of growth drives the amino acid economy of Caenorhabditis elegans. Proc Biol Sci. 2009, 276: 2747-2754. 10.1098/rspb.2009.0354.
Antebi A: Ageing - When less is more. Nature. 2007, 447: 536-537. 10.1038/447536a.
Kenyon C: The plasticity of aging: Insights from long-lived mutants. Cell. 2005, 120: 449-460. 10.1016/j.cell.2005.02.002.
Patel DS, Garza-Garcia A, Nanji M, McElwee JJ, Ackerman D, Driscoll PC, Gems D: Clustering of genetically defined allele classes in the Caenorhabditis elegans DAF-2 insulin/IGF-1 receptor. Genetics. 2008, 178: 931-946. 10.1534/genetics.107.070813.
Li WQ, Kennedy SG, Ruvkun G: daf-28 encodes a C-elegans insulin superfamily member that is regulated by environmental cues and acts in the DAF-2 signaling pathway. Genes Dev. 2003, 17: 844-858. 10.1101/gad.1066503.
Pierce SB, Costa M, Wisotzkey R, Devadhar S, Homburger SA, Buchman AR, Ferguson KC, Heller J, Platt DM, Pasquinelli AA, Liu LX, Doberstein SK, Ruvkun G: Regulation of DAF-2 receptor signaling by human insulin and ins-1, a member of the unusually large and diverse C-elegans insulin gene family. Genes Dev. 2001, 15: 672-686. 10.1101/gad.867301.
Larsen PL, Albert PS, Riddle DL: Genes That Regulate Both Development and Longevity in Caenorhabditis-Elegans. Genetics. 1995, 139: 1567-1583.
McCulloch D, Gems D: Sex-specific effects of the DAF-12 steroid receptor on aging in Caenorhabditis elegans. Ann N Y Acad Sci. 2007, 1119: 253-259. 10.1196/annals.1404.018.
Scott BA, Avidan MS, Crowder CM: Regulation of hypoxic death in C-elegans by the insulin/IGF receptor homolog DAF-2. Science. 2002, 296: 2388-2391. 10.1126/science.1072302.
Adams TE, Epa VC, Garrett TPJ, Ward CW: Structure and function of the type 1 insulin-like growth factor receptor. Cell Mol Life Sci. 2000, 57: 1050-1093. 10.1007/PL00000744.
Behm CA: The role of trehalose in the physiology of nematodes. Int J Parasitol. 1997, 27: 215-229. 10.1016/S0020-7519(96)00151-8.
Pellerone FI, Archer SK, Behm CA, Grant WN, Lacey MJ, Somerville AC: Trehalose metabolism genes in Caenorhabditis elegans and filarial nematodes. Int J Parasitol. 2003, 33: 1195-1206. 10.1016/S0020-7519(03)00173-5.
Singer MA, Lindquist S: Thermotolerance in Saccharomyces cerevisiae: the Yin and Yang of trehalose. Trends Biotechnol. 1998, 16: 460-468. 10.1016/S0167-7799(98)01251-7.
Lamitina ST, Strange K: Transcriptional targets of DAF-16 insulin signaling pathway protect C-elegans from extreme hypertonic stress. American Journal of Physiology-Cell Physiology. 2005, 288: C467-C474. 10.1152/ajpcell.00451.2004.
Craig H, Isaac RE, Brooks DR: Unravelling the moulting degradome: new opportunities for chemotherapy?. Trends in Parasitol. 2007, 23: 248-253. 10.1016/j.pt.2007.04.003.
Myllyharju J, Kivirikko KI: Collagens, modifying enzymes and their mutations in humans, flies and worms. Trends Genet. 2004, 20: 33-43. 10.1016/j.tig.2003.11.004.
Ackerstaff E, Glunde K, Bhujwalla ZM: Choline phospholipid metabolism: A target in cancer cells?. J Cell Biochem. 2003, 90: 525-533. 10.1002/jcb.10659.
Burnell AM, Houthoofd K, O'Hanlon K, Vanfleteren JR: Alternate metabolism during the dauer stage of the nematode Caenorhabditis elegans. Experimental Gerontol. 2005, 40: 850-856. 10.1016/j.exger.2005.09.006.
Oriordan VB, Burnell AM: Intermediary Metabolism in the Dauer Larva of the Nematode Caenorhabditis-Elegans .1. Glycolysis, Gluconeogenesis, Oxidative-Phosphorylation and the Tricarboxylic-Acid Cycle. Comp Biochem Physiol Biochem Mol Biol. 1989, 92: 233-238. 10.1016/0305-0491(89)90271-X.
Wang J, Kim SK: Global analysis of dauer gene expression in Caenorhabditis elegans. Development. 2003, 130: 1621-1634. 10.1242/dev.00363.
Melendez A, Talloczy Z, Seaman M, Eskelinen EL, Hall DH, Levine B: Autophagy genes are essential for dauer development and life-span extension in C-elegans. Science. 2003, 301: 1387-1391. 10.1126/science.1087782.
Lin K, Dorman JB, Rodan A, Kenyon C: daf-16: An HNF-3/forkhead family member that can function to double the life-span of Caenorhabditis elegans. Science. 1997, 278: 1319-1322. 10.1126/science.278.5341.1319.
Ogg S, Paradis S, Gottlieb S, Patterson GI, Lee L, Tissenbaum HA, Ruvkun G: The Fork head transcription factor DAF-16 transduces insulin-like metabolic and longevity signals in C-elegans. Nature. 1997, 389: 994-999. 10.1038/40194.
Libina N, Berman JR, Kenyon C: Tissue-specific activities of C-elegans DAF-16 in the regulation of lifespan. Cell. 2003, 115: 489-502. 10.1016/S0092-8674(03)00889-4.
Lin Y, Huang X: Motor neuron axon guidance in Caenorhabditis elegans. Dev Biol. 2001, 235: 231-232.
Griffin TJ, Gygi SP, Ideker T, Rist B, Eng J, Hood L, Aebersold R: Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae. Mol Cell Proteomics. 2002, 1: 323-333. 10.1074/mcp.M200001-MCP200.
Gygi SP, Rochon Y, Franza BR, Aebersold R: Correlation between protein and mRNA abundance in yeast. Mol Cell Biol. 1999, 19: 1720-1730.
Lee MH, Park H, Shim G, Lee J, Koo HS: Regulation of gene expression, cellular localization, and in vivo function of Caenorhabditis elegans DNA topoisomerase I. Genes to Cells. 2001, 6: 303-312. 10.1046/j.1365-2443.2001.00423.x.
Cypser JR, Johnson TE: Hormesis in Caenorhabditis elegans dauer-defective mutants. Biogerontology. 2003, 4: 203-214. 10.1023/A:1025138800672.
Weinkove D, Halstead JR, Gems D, Divecha N: Long-term starvation and ageing induce AGE-1/PI 3-kinase-dependent translocation of DAF-16/FOXO to the cytoplasm. Bmc Biology. 2006, 4: 1-10.1186/1741-7007-4-1.
Wylie T, Martin J, Abubucker S, Yin Y, Messina D, Wang ZY, McCarter JP, Mitreva M: NemaPath: online exploration of KEGG-based metabolic pathways for nematodes. Bmc Genomics. 2008, 9: 525-10.1186/1471-2164-9-525.
Fisher AL, Page KE, Lithgow GJ, Nash L: The Caenorhabditis elegans K10C2.4 gene encodes a member of the fumarylacetoacetate hydrolase family - A caenorhabditis elegans model of type i tyrosinemia. Journal of Biological Chemistry. 2008, 283: 9127-9135. 10.1074/jbc.M708341200.
Calvo AC, Pey AL, Ying M, Loer CM, Martinez A: Anabolic function of phenylalanine hydroxylase in Caenorhabditis elegans. Faseb Journal. 2008, 22: 3046-3058. 10.1096/fj.08-108522.
Payne SH, Loomis WF: Retention and loss of amino acid biosynthetic pathways based on analysis of whole-genome sequences. Eukaryot Cell. 2006, 5: 272-276. 10.1128/EC.5.2.272-276.2006.
Brosnan JT, Brosnan ME: Branched-chain amino acids: enzyme and substrate regulation. J Nutr. 2006, 136: 207S-211S.
Newgard CB, An J, Bain JR, Muehlbauer MJ, Stevens RD, Lien LF, Haqq AM, Shah SH, Arlotto M, Slentz CA, Rochon J, Gallup D, Ilkayeva O, Wenner BR, Yancy WS, Eisenson H, Musante G, Surwit RS, Millington DS, Butler MD, Svetkey LP: A branched-chain amino acid-related metabolic signature that differentiates obese and lean humans and contributes to insulin resistance. Cell Metab. 2009, 9: 311-326. 10.1016/j.cmet.2009.02.002.
Fraser AG, Kamath RS, Zipperlen P, Martinez-Campos M, Sohrmann M, Ahringer J: Functional genomic analysis of C. elegans chromosome I by systematic RNA interference. Nature. 2000, 408: 325-330. 10.1038/35042517.
Chuang DT, Shih VE: Maple syrup urine disease (branched-chain ketoaciduria). The Metabolic and Molecular Bases of Inherited Disease. Edited by: Scriver CR, Beaudet AL, Sly WS, Valle D. 2001, New York: McGraw-Hill, 2: 1971-2005, 8
Lafaye A, Junot C, Pereira Y, Lagniel G, Tabet JC, Ezan E, Labarre J: Combined proteome and metabolite-profiling analyses reveal surprising insights into yeast sulfur metabolism. J Biol Chem. 2005, 280: 24723-24730. 10.1074/jbc.M502285200.
Coquin L, Feala JD, McCulloch AD, Paternostro G: Metabolomic and flux-balance analysis of age-related decline of hypoxia tolerance in Drosophila muscle tissue. Mol Syst Biol. 2008, 4: 233-10.1038/msb.2008.71.
Feala JD, Coquin L, McCulloch AD, Paternostro G: Flexibility in energy metabolism supports hypoxia tolerance in Drosophila flight muscle: metabolomic and computational systems analysis. Mol Syst Biol. 2007, 3: 99-10.1038/msb4100139.
Steuer R, Gross T, Selbig J, Blasius B: Structural kinetic modeling of metabolic networks. Proc Natl Acad Sci USA. 2006, 103: 11868-11873. 10.1073/pnas.0600013103.
Brenner S: The genetics of Caenorhabditis elegans. Genetics. 1974, 77: 71-94.
Beckonert O, Keun HC, Ebbels TM, Bundy J, Holmes E, Lindon JC, Nicholson JK: Metabolic profiling, metabolomic and metabonomic procedures for NMR spectroscopy of urine, plasma, serum and tissue extracts. Nat Protoc. 2007, 2: 2692-2703. 10.1038/nprot.2007.376.
Weljie AM, Newton J, Mercier P, Carlson E, Slupsky CM: Targeted profiling: quantitative analysis of 1H NMR metabolomics data. Anal Chem. 2006, 78: 4430-4442. 10.1021/ac060209g.
Dieterle F, Ross A, Schlotterbeck G, Senn H: Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in H-1 NMR metabonomics. Analytical Chemistry. 2006, 78: 4281-4290. 10.1021/ac051632c.
Purohit PV, Rocke DM, Viant MR, Woodruff DL: Discrimination models using variance-stabilizing transformation of metabolomic NMR data. Omics-a Journal of Integrative Biology. 2004, 8: 118-130. 10.1089/1536231041388348.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-484. 10.1093/nar/gkm882.
Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006, 34: D354-357. 10.1093/nar/gkj102.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.
SD and JGB were funded by the Natural Environment Research Council (UK).
SF and SKD carried out metabolomic studies and primary analysis of data and helped write the paper. JMV and JSS carried out further analyses. JGB and AML conceived and coordinated the study, constructed the figures, and wrote the paper. All authors read and approved the final manuscript.
Silke Fuchs, Jacob G Bundy contributed equally to this work.