Suppression of transcriptional drift extends C. elegans lifespan by postponing the onset of mortality
Abstract
Longevity mechanisms increase lifespan by counteracting the effects of aging. However, whether longevity mechanisms counteract the effects of aging continually throughout life, or whether they act during specific periods of life, preventing changes that precede mortality is unclear. Here, we uncover transcriptional drift, a phenomenon that describes how aging causes genes within functional groups to change expression in opposing directions. These changes cause a transcriptome-wide loss in mRNA stoichiometry and loss of co-expression patterns in aging animals, as compared to young adults. Using Caenorhabditis elegans as a model, we show that extending lifespan by inhibiting serotonergic signals by the antidepressant mianserin attenuates transcriptional drift, allowing the preservation of a younger transcriptome into an older age. Our data are consistent with a model in which inhibition of serotonergic signals slows age-dependent physiological decline and the associated rise in mortality levels exclusively in young adults, thereby postponing the onset of major mortality.
https://doi.org/10.7554/eLife.08833.001eLife digest
All organisms age, leading to gradual declines in the body’s systems and eventually death. How certain genetic mutations and drugs delay the effects of aging and promote survival to an older age is a question many researchers are exploring. One way this problem is investigated is by looking at how the activity – or expression – of different genes changes during aging.
Scientists interested in understanding aging and longevity often study a simple worm called Caenorhabditis elegans. This worm normally lives for about three weeks, and young C. elegans are able to produce offspring within days of hatching. This accelerated life cycle allows scientists to observe the entire lifespan of the worms. Over time, experiments have shown that DNA damage, changes in behavior and changes to gene expression are all markers of aging in the worms.
Now, Rangaraju et al. describe how changes in gene expression patterns that begin early in the lives of C. elegans shorten their lifespan. Specifically, in groups of genes that work together, some genes increase expression, while others decrease expression with age. This phenomenon is called “transcriptional drift” and leads to an age-associated loss of coordination among groups of genes that help orchestrate specific tasks.
Rangaraju et al. show that an antidepressant called mianserin prevents transcriptional drift in many of C. elegans’ genes: young worms treated with the drug resist the effects of aging on the transcriptome and maintain coordinated patterns of gene expression for longer. Maintaining coordinated patterns of gene expression postpones the onset of age-related bodily declines and extends the life of treated worms by extending the duration of young adulthood and postponing the onset of age-associated death. The drug also appears to protect against stress-induced changes in gene expression. This suggests that some of the age-related shifts in gene expression occur when cells fail to recover normal gene expression patterns after a stressful event.
Questions that remain to be investigated in future studies are whether other longevity mechanisms also extend lifespan by preserving coordinated gene expression patterns, and whether other longevity mechanisms act by extending specific periods of life.
https://doi.org/10.7554/eLife.08833.002Introduction
The most widely used standard to measure aging of an organism is the quantification of lifespan (Partridge and Gems, 2007). Lifespan relates to aging, as the latter causes the degeneration of tissues and organs, thereby increasing mortality due to systemic functional tissue failure (Balch et al., 2008; Bishop et al., 2010; David et al., 2010; Haigis and Sweet-Cordero, 2011; Taylor and Dillin, 2011; Gladyshev, 2013; Burkewitz et al., 2015; Currais, 2015). Several genetic and pharmacological strategies have been shown to prolong the lifespan of various organisms, including C. elegans (Kenyon et al., 1993; Kaeberlein et al., 1999; Curran and Ruvkun, 2007; Evason et al., 2008; Onken and Driscoll, 2010; Alavez et al., 2011; Chin et al., 2014; Ye et al., 2014; Tatum et al., 2015). Mutations in age-1 or daf-2, for example, slow degenerative processes occurring throughout life, thereby constantly lowering mortality rates (Johnson, 1990; Kenyon et al., 1993; Taylor et al., 2014). Age-associated degenerative processes such as a decline in proteostatic capacity are not necessarily restricted to older organisms but can also be observed in young adults (Labbadia and Morimoto, 2015a; 2015b). This raises the possibility of degenerative processes that occur only in young adults and thus specifically contribute to the rise of mortality during young adulthood. Any longevity mechanisms preventing such a degenerative process would specifically slow mortality rates during the period of young adulthood, effectively prolonging its duration to postpone the onset of major age-associated mortality around midlife (Bartke, 2015). However, to identify such mechanisms would require mortality-independent metrics of age-associated change, as age-associated mortality rates during young adulthood are difficult to determine by demographic analysis against the back drop of non-aging-related death events (Partridge and Gems, 2007; Beltran-Sancheza et al., 2012).
In C. elegans, mortality-independent metrics of aging include age-associated decline of various behaviors or physiological parameters such as movement or stress resistance (Huang et al., 2004; Bansal et al., 2015). Molecular markers of aging include sets of genes whose expression change with age, such as micro-RNAs, electron transport chain (ETC) components, or genes involved in posttranslational modifications such as methylation (Budovskaya et al., 2008; de Magalhaes et al., 2009; Pincus et al., 2011; Horvath et al., 2015). However, aging also increases DNA damage, affects nuclear architecture, chromatin complexes, chromatin modifications, and the transcriptional machinery (Mostoslavsky et al., 2006; Scaffidi and Misteli, 2006; Feser et al., 2010; Greer et al., 2011; Maures et al., 2011; Fushan et al., 2015). Therefore, an emerging alternative approach to measure specific gene expression changes with age is to quantify the progressive imbalance in gene expression patterns as a function of age. Two such approaches, one measuring transcriptional noise, the cell-to-cell variation in gene expression, and the other measuring decreasing correlation in the expression of genetic modules, showed a loss of co-expression patterns with age (Bahar et al., 2006; Southworth et al., 2009). These studies suggest that age-associated changes can be measured independently from mortality by tracking the loss of gene expression patterns that are observed in young animals.
In the present study, we set out to investigate the mechanisms by which the atypical antidepressant mianserin extends lifespan by recording the transcriptional dynamics of mianserin-treated and untreated C. elegans across different ages. These studies revealed that aging causes transcriptional drift, an evolutionarily conserved phenomenon in which the expression of genes change in opposing directions within functional groups. These changes cause a transcriptome-wide loss in mRNA stoichiometry and loss of co-expression patterns in aging animals, as compared to young adults. Mianserin treatment reduced age-associated transcriptional drift across ~80% of the transcriptome, preserving many characteristics of transcriptomes of younger animals. We used transcriptional drift along with mortality analysis as metrics to monitor aging and find that mianserin treatment extended lifespan by exclusively slowing age-associated changes in young adults, thereby postponing the onset of mortality.
Results
Aging causes a loss of co-expression patterns observed in young adults
To better understand how aging changes gene expression patterns in a eukaryotic organism, and how these changes are affected by longevity, we measured gene expression changes in mianserin-treated or untreated C. elegans by RNA-sequencing (RNA-seq; Figure 1a). Cohort #1 was a time series to study how gene expression patterns change over time in control (water) animals or in animals treated with mianserin on day 1 of adulthood (24 hr after L4 stage). Cohort #2 was designed to study dosage effects of increasing concentrations of mianserin with aging, and cohort #3 was designed to study the effects of delayed mianserin-treatment of worms treated at day 3 or 5 of adulthood (Figure 1a). Lifespan of a sub-population of each cohort was simultaneously assessed to ensure the effect of mianserin.
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.08833.004
-
Figure 1—source data 2
- https://doi.org/10.7554/eLife.08833.005
-
Figure 1—source data 3
- https://doi.org/10.7554/eLife.08833.006
-
Figure 1—source data 4
- https://doi.org/10.7554/eLife.08833.007
-
Figure 1—source data 5
- https://doi.org/10.7554/eLife.08833.008
Comparison of gene expression profiles of age-matched mianserin-treated and untreated controls, showed that approximately 3,000–6,000 genes changed with age in response to mianserin treatment (FDR<0.1, Figure 1—source data 1) (Robinson and Oshlack, 2010; Kim et al., 2013; Lawrence et al., 2013). We separated genes into sets that showed increased or decreased expression in response to mianserin, to conduct gene-set enrichment analysis. This revealed hundreds of gene ontologies (GO) that changed in response to mianserin (Figure 1—source data 2) (Ashburner et al., 2000; Mi et al., 2005). We observed that many GOs were enriched for both, genes that increased as well as decreased as a consequence of aging. This observation complicated any interpretation on whether pathways were activated or inhibited in response to mianserin, and how the associated function (GO) relates to mianserin-induced lifespan extension (Figure 1b).
We observed a similar scenario by conducting gene-set enrichment analysis for gene expression changes in response to age in untreated animals. As seen with mianserin, many GO annotations were enriched for both up- as well as downregulated genes at any given age (Figure 1c; Figure 1—source data 3, 4), making it difficult to interpret whether those pathways are being activated or inhibited with age. We generated 50 representative pie charts out of the 249 GO annotations that contained genes that increased or decreased in expression by day 10 due to aging. These charts suggested that as animals age and become older, genes change expression in opposing directions, disrupting relative mRNA ratios within the GO, when compared to young adults (Figure 1—figure supplement 1). Thus, aging changed the stoichiometric relationship between mRNAs belonging to the same functional group (GO). In many cases, the fractions of genes that increased, decreased or did not change in expression showed no consistent pattern, nor provided any insight into the pathway activity (Figure 1—figure supplement 1).
Because the expression patterns observed in many GOs were difficult to interpret in terms of functional change, we turned to investigate expression changes in the superoxide detoxification pathway, a well-defined cellular function that declines with age (Ashburner et al., 2000; Mi et al., 2005; Kumsta et al., 2011; Bansal et al., 2015; Rangaraju et al., 2015a). As expected from our previous studies (Rangaraju et al., 2015a), the expression levels of some superoxide detoxification genes were higher in mianserin-treated animals compared to age-matched controls (Figure 1d). Exceptions were the expression levels of sod-4 and sod-5, which were lowered upon mianserin treatment (Figure 1d). However, plotting expression changes of superoxide detoxification genes as a function of age (Figure 1e, left panel) revealed again a scenario in which genes changed in opposing directions as seen in the pie charts for many GOs before (Figure 1—figure supplement 1). Some mRNAs including those of sod-4, -5 increased with age, while some decreased (sod-1, -2, prdx-2, 3, 6) and some did not change (ctl-1, 2, 3), leading to an overall 5-10-fold change in stoichiometric balance among superoxide detoxification-associated mRNAs by day 5 (Figure 1e, left panel). More interestingly, if the expression of an sod increased with age, mianserin treatment prevented the increase and if the expression of an sod decreased with age, mianserin prevented the decrease (Figure 1e, right panel). Thus, when we took the mRNA expression levels of young animals into account, the emerging picture suggested that mianserin treatment attenuated age-associated gene expression changes.
We therefore asked whether the complex gene-set enrichment patterns observed comparing mianserin-treated and untreated samples (Figure 1b,c) could be explained by mianserin preventing expression changes due to age. Indeed, many GO annotations that increased expression with age were decreased by mianserin treatment and vice versa (Figure 1f). This attenuation of age-associated changes by mianserin treatment was even more pronounced for individual genes (Figure 1g). Analyzing cohort #1 showed a significant change in expression levels of 3,367 genes, as the animals aged from day 1 to day 3, and a change in 5,947 genes from day 1 to day 10 (FDR < 0.1) (Figure 1g, significant genes only). Mianserin treatment reduced these age-associated expression changes in over 90% of cases. Including all age-associated expression changes for the 19,196 different transcripts present in our data-set, we found that mianserin treatment attenuated age-associated changes in transcription in 15,095 out of 19,169 genes (80%, binomial P < 10–100). Thus, most of the changes observed between mianserin-treated and untreated animals are due to mianserin preventing transcriptional changes with age.
When we excluded all genes that changed due to age and were attenuated by mianserin, we obtained a much smaller gene-set consisting of mianserin-induced changes that was enriched for GOs related to stress, xenobiotic and immune-responses, as well as genes associated with aging and the determination of lifespan (Table 1, Figure 1—source data 5). These GOs have been previously shown to be regulated by serotonin in C. elegans with the exception of the xenobiotic response (Zhang et al., 2005; Petrascheck et al., 2007; Rangaraju et al., 2015a). Thus, accounting for age-associated transcriptional changes dramatically simplified a seemingly very complex gene-expression pattern (Figure 1b,c). It revealed that mianserin affected expression of a small set of physiological functions that are known to be regulated by serotonin and have been shown to be required for mianserin-induced lifespan extension or for aging in general (Garsin et al., 2003; Rangaraju, et al., 2015; Petrascheck, et al., 2007) (Table 1; Figure 1f; Figure 1—source data 5).
Based on these observations, we classified gene expression changes for groups of genes into two types. Type I changes describe whether the overall expression across an entire functional group/pathway increases or decreases i.e. whether the pathway is up or down regulated with age. Type II changes describe the relative changes in gene expression among genes within functional groups with respect to each other. We named the type II change transcriptional drift. As animals age, genes within functional groups change expression levels in opposing directions resulting in the disruption of the co-expression patterns seen in young adults.
To analyze the effects of aging on transcriptional drift (type II), we designed graphs that plot the log-fold changes (log [old/young reference day1]) in gene expression as a function of age. Such a plot can be constructed for whole transcriptomes as well as for any functional subset of genes, for example, genes involved in oxidative phosphorylation or lysosome biology (Figure 1h,i). In young adults, the log-fold change is 0 and values close to 0 therefore suggest gene expression as seen in young adults (Figure 1h,i). To quantify transcriptional drift changes with age (type II), we calculated the variance of the log-fold change for genes involved in each pathway. For the purpose of this study, we will refer to this variance as drift-variance (see Materials and methods). If gene expression ratios within a pathway stay constant with age, drift-variance will stay small. If a majority of genes within a pathway change expression in opposing directions or if the rates by which they change differ dramatically, drift-variance will increase. Note that “transcriptional drift” is different from “transcriptional noise” in that the former analyzes variance among genes within the same biological replicates, whereas the latter analyzes variance of the same genes among biological replicates. Hence, how far the aging transcriptome deviates away from the transcriptome seen in young adults can be graphed in a Tukey-style box plot, which plots the drift-variance as a function of age (Figure 1h,i). We will refer to these plots as drift-plots (Figure 1h; Figure 2—figure supplement 1a–d).
Longevity mechanisms attenuate transcriptional drift-variance
We constructed drift-plots for all 19,196 genes in the data of cohort #1, which revealed a dramatic increase in drift-variance with age, showing a progressive loss of mRNA stoichiometries and co-expression patterns observed in young-adults (Figure 2a, shaded region encompassing the whiskers of Tukey-plot). This effect was also seen in other publicly available data-sets of aging C. elegans transcriptomes and drift-variance continued to increase with age at least until day 20 (Figure 2—figure supplement 1e). Mianserin treatment attenuated the effect of aging across the whole transcriptome and preserved the co-expression patterns observed in young-adults into later age. To test whether transcriptional drift is driven by a small subset of mRNAs or a transcriptome-wide phenomenon, we randomly divided the transcriptome into subsamples of ~1,000 genes. Each subsample showed identical increases in drift-variance with age, confirming a transcriptome-wide effect (Figure 2—figure supplement 1f).
We previously showed, that the effect of mianserin to extend lifespan is dose-dependent (Petrascheck et al., 2007). To explore a possible quantitative relationship between longevity and drift-variance, we generated drift-plots for transcriptomes of animals treated with increasing doses of mianserin (Figure 1a, cohort #2). Increasing doses of mianserin progressively increased longevity and decreased drift-variance as measured in 5-day-old animals (Figure 2b,c; Table 2). Thus, remarkably, by varying the dose of a single molecule, it was possible to control the degree to which aging drives the loss of transcriptional co-expression away from patterns observed in young adults. These results suggested a quantitative relationship between mianserin-induced longevity and its effect on drift-variance.
Our previous studies had also shown that mianserin does not extend lifespan when added to 5-day-old post-reproductive adult animals (Petrascheck et al., 2007). Thus, we next tested whether mianserin attenuates transcriptional drift-variance independently of longevity by treating older animals. Mianserin did not attenuate transcriptional drift-variance when added on day 5 (Figure 2d). Adding mianserin on day 3 of adulthood caused a small extension of lifespan and a corresponding small attenuation of drift-variance, further supporting a quantitative relationship between suppression of drift-variance and extension of lifespan (Figure 1a, cohort #3, Figure 2d,e; Table 2). However, mianserin fully induced the xenobiotic response by up to 1,000-fold irrespective of whether added on day 1 or day 5 (Figure 2f). Therefore, the lack of an effect of mianserin when added to day 5 adults cannot be attributed to reduced drug uptake. Taken together, these results show that mianserin does not attenuate drift-variance when it does not extend lifespan.
We next asked whether the attenuation of drift-variance is unique to mianserin or whether it is observed in other lifespan-extension paradigms (Figure 2g). We asked whether reduced insulin signaling also attenuates drift-variance by analyzing the previously published gene expression data-sets of long-lived C. elegans daf-2 RNAi-treated and vector control animals (Murphy et al., 2003). Analyses of drift-variance for these data-sets showed that treatment with daf-2 RNAi attenuated drift-variance (Figure 2g). Moreover, mianserin and daf-2 RNAi attenuated age-associated drift of overlapping sets of genes. Of the 6,958 genes for which expression levels were detected at all ages in both data-sets, 58% (4,078 genes, binomial P= 6.3e-47) were attenuated by both longevity-extending mechanisms. This overlap is consistent with experiments showing that these two longevity mechanisms partially overlap, potentially explaining why mianserin only causes a +11% lifespan extension in daf-2(e1370) mutant animals instead of 31% seen in the parallel wild-type experiments (Petrascheck et al., 2007). Thus, lifespan extension by mianserin or daf-2 RNAi attenuates transcriptional drift in overlapping sets of genes.
Conversely, suppressing longevity by daf-16(RNAi) prevented the attenuation of drift-variance by daf-2(RNAi) and increased it beyond what was seen in control animals (Figure 2g). Thus, the activation of DAF-16 target genes leads to the attenuation of transcriptional drift in thousands of genes across the transcriptome. Taken together, these results show that drift-variances increase with age in C. elegans and are attenuated in two different longevity paradigms (Figure 2a,g).
From a technical perspective, the comparison between the mianserin data and the Murphy data (Murphy et al., 2003) also shows that the phenomenon of transcriptional drift is robust enough not to be influenced by the presence of eggs in the animals or the method of sterilization, as our study used FUDR and the Murphy et al. (2003) study used sterile mutants (Figure 2a,g; Figure 2—figure supplement 2).
Attenuating drift-variances in redox-pathways preserves homeostatic capacity
The results above suggested that preserving low drift-variance in transcriptomes preserves longevity. We therefore asked whether attenuating drift-variance in specific pathways preserves homeostatic capacity, the ability of pathways to appropriately respond to a stimulus or stress. Throughout life, organisms respond to stimuli by activating or repressing transcriptional programs, an ability that is lost with age. We hypothesized that one way by which regulatory ability may be lost could be due to a failure to return to their precise steady-state transcriptional levels after stimulation. This would give rise to increases in drift-variance (Figure 3a), as seen in the drift plots for oxidative phosphorylation or lysosome biology (Figure 1h,i). In this model, slight initial deviations in gene expression levels would be compounded over time resulting in imbalanced stoichiometries between pathway components resulting in functional decline with age (Figure 3a).
Our previous studies showed that mianserin protected C. elegans from oxidative stress by a neuronal mechanism that modulated peripheral stress response genes (NEUROX) (Rangaraju et al., 2015a). We therefore constructed drift plots for redox-associated pathways that showed that mianserin indeed increased the overall expression of oxidative stress response genes (type I) relative to age-matched controls but also attenuated transcriptional drift (type II) (Figure 3b; Table 3).
We therefore asked whether mianserin treatment increased resistance to oxidative stress by either directly activating the oxidative stress response or whether attenuating transcriptional drift would preserve homeostatic capacity into older age (Rahman et al., 2013). Animals were treated with water or mianserin on day 1 of adulthood, followed by treatment with the reactive oxygen species (ROS) generator paraquat on day 3, 5, or 10 (Figure 3c). On day 3 of adulthood, no difference in stress resistance between mianserin-treated and untreated animals was observed. As animals grew older (day 5 and day 10), mianserin treatment greatly improved stress resistance (Figure 3c; Table 4). Again, as with lifespan, delaying the start of mianserin treatment to day 3 and day 5 progressively reduced its protective effect on stress resistance, this time measured in animals subjected to paraquat on day 10 of adulthood (Figure 3d; Table 5). Thus, mianserin treatment specifically improves stress resistance in older (day 5 and day 10) but not in younger (day 3) animals consistent with a model in which it preserves the homeostatic capacity of redox function.
To further distinguish between a model in which mianserin directly activates an oxidative stress response from one that preserves the homeostatic capacity by attenuating drift-variance, we asked whether mianserin enhanced (direct activation) or attenuated (preserving capacity) genes that change in response to oxidative stress (Figure 3e). Oliveira et al. identified 252 genes that were upregulated and 88 genes that were downregulated in young C. elegans in response to oxidative stress, and can therefore be considered an experimentally determined oxidative stress signature (Oliveira et al., 2009). We hypothesized that a direct activation of the oxidative stress response by mianserin would mimic the increase in expression of the 252 genes and the decrease in the expression of the 88 genes as seen in response to oxidative stress. However, we observed an attenuation rather than an activation of the oxidative stress signatures, consistent with preserving homeostatic capacity rather than a direct activation. Genes that increased in response to oxidative stress (252) showed a lower expression while genes that decreased (88) in response to oxidative stress showed a higher expression in age-matched mianserin-treated animals (Figure 3e). Consistent with the functional data, differences in the oxidative stress signature were only observed in older animals (day 5, 10), but not in younger day 3 animals. These results are consistent with a model in which mianserin treatment preserves the redox system from age-associated decline, thus improving redox capacity in older age.
Mianserin requires the serotonin receptor SER-5 to preserve low drift-variances
In mammals, mianserin antagonizes serotonergic signals sent by 5-HT2A/C receptors (Gillman, 2006). We next asked whether preservation of redox capacity and reducing drift-variance in redox pathways by mianserin depends on serotonergic signaling. To identify the serotonergic receptor, we treated multiple mutants, each deficient in signaling by a single G-protein coupled receptor (GPCR) with mianserin on day 1, followed by increasing concentrations of paraquat on day 5 to induce oxidative stress (Figure 4a,b; Table 6). Mianserin was unable to protect multiple ser-5 mutant alleles (ok3087, tm2647, tm2654) from oxidative stress (Figure 4a,b; Figure 4—figure supplement 1a; Table 6). In addition, seven structurally distinct serotonergic antagonists/inverse agonists also protect from oxidative stress in a ser-5 dependent manner (Figure 4—figure supplement 1b; Table 7). Furthermore, mianserin did not protect animals unable to synthesize serotonin (tph-1(mg280)) (Figure 4a; Table 6) (Sze et al., 2000).
We next asked whether SER-5 was also required for mianserin to preserve low transcriptional drift-variances in redox-related genes. We measured redox gene expression levels by qRT-PCR in wild-type 5-day-old N2 and ser-5(ok3087) animals that were treated with mianserin or water on day 1 (Figure 5a,b; Figure 1—figure supplement 1a). In N2 samples, mianserin increased the expression of stress response genes that drift down with age (sod-1, sod-2, prdx-2, -3, -6) and decreased the expression of stress response genes that drift up with age (sod-4, sod-5, all hsp-16s), an effect that was not observed in ser-5(ok3087) mutants. In contrast, SER-3 and SER-4, two receptors we previously showed to be required for lifespan extension by mianserin, were dispensable for stress protection (Figure 4a,b) (Petrascheck et al., 2007), as well as for the attenuation of drift-variance in redox-associated genes (Figure 4—figure supplement 1c). Thus, in wild-type animals, mianserin treatment preserved low drift-variances in redox-related genes into older age (day 5), in a ser-5 dependent manner (Figure 5a,b).
Importantly, ser-5 mutants were specifically defective in their response to mianserin, but showed no defect in their response to oxidative stress. Young (day 1) wild-type N2 animals and ser-5(ok3087) mutants showed a nearly identical response to oxidative stress (Figure 5c). The age-specific effects of ser-5 could not be attributed to expression changes, as ser-5 expression remained constant from day 1 to day 10 in our RNA-seq experiment.
To test the hypothesis that mianserin preserved the homeostatic capacity of the redox system, as suggested by Figure 3e, we asked whether the treatment with mianserin on day 1 of adulthood led to an enhanced redox gene expression in response to the stressor paraquat in older animals (day 5). We therefore challenged older mianserin-treated or control animals (day 5) with paraquat for 8 hr and measured redox-gene expression by qRT-PCR (Figure 5d). Mianserin treatment led to an enhanced transcription of redox genes in response to paraquat as compared to age-matched control animals. The enhanced response was ser-5 dependent (Figure 5d). Thus, SER-5 is required for mianserin to attenuate age-associated increases in drift-variance in redox genes, and to preserve the homeostatic capacity of the redox system into older age.
Furthermore, lifespan-extension by mianserin was strongly reduced or abrogated in ser-5, snt-1 and unc-26 mutant animals (Figure 5e, f; Figure 5—figure supplement 1a; Table 8). Seven additional serotonergic antagonists/inverse agonists also extended lifespan in a manner that was partially or fully dependent on ser-5 (Figure 5—figure supplement 1b). Thus, these results show that inhibiting serotonergic signals via SER-5 extends lifespan, attenuates age-associated drift-variance in the redox system and preserves the homeostatic capacity of the redox system.
Mianserin prolongs lifespan by slowing age-associated change in young adults
We next asked whether drift-variance could be used as a metric to monitor age-associated change in young adults. Comparing drift-variances between mianserin-treated and untreated animals, we noticed that by day 10, mianserin-treated animals exhibited a drift-variance slightly lower than that of 3-day-old control animal (P=0.37). This suggested that mianserin-treated animals showed a ~7–8 day delay in age-associated transcriptional change compared to age-matched controls (Figure 2a).
Principle component analysis (PCA), a different statistical method to analyze differences between transcriptomes, confirmed this observation (Figure 6a). PCA showed that control samples aligned on the x-axis (dimension 1) according to age and that 10 day-old mianserin-treated animals aligned closer to 3-day-old than to 10-day-old control animals. These results suggested that the physiological shift that results in the 7–8 day lifespan extension observed in mianserin-treated animals at the end of a lifespan assay was already observable by day 10.
We therefore asked whether mianserin slowed age-associated physiological change specifically in early adulthood causing a 7–8 day delay by day 10. If so, mianserin would be expected to specifically lower the mortality rate in young but not in old adults. However, the number of age-associated death events in young adults is too low to directly determine changes in age-associated mortality rates before the age of day 10. As we are comparing mortality in animals either treated with water or mianserin that is added to the same population of worms on day 1 of adulthood, we can confidently state that mortality levels are identical between mianserin-treated and untreated adults at the start of the experiment. Any difference in mortality levels observed from day 1 onwards must therefore be the result of a change in mortality rate by mianserin.
Plotting a mortality curve for over 3,000 mianserin-treated or untreated animals showed a significantly lower mortality level for mianserin-treated animals by day 12 (Figure 6b, Figure 6—figure supplement 1a). Therefore, mianserin treatment decelerated the rise in mortality levels between day 1 and 12 of adulthood. From then on, the mortality curves were parallel showing a 7–8 day shift in mortality across the remaining lifespan. The parallel nature suggested that mianserin did not affect mortality rates past day 12 and that its effect on lifespan was restricted to the period of early adulthood (Figure 6b, Figure 6—figure supplement 1a) (Mair et al., 2003; Vaupel, 2010). Power calculations confirmed that these mortality curves were sufficiently powered to detect a one day difference in lifespan in over 90% of the experiments (α=0.01) (Figure 6—figure supplement 1b) (Ye et al., 2014). These results further supported a model in which mianserin treatment specifically lowered age-associated change in early adulthood, causing a shift in physiology and mortality that can be observed in transcriptomes by day 10.
We reasoned that if the effect of mianserin on lifespan precedes the onset of mortality and is completed by day 10, mianserin treatment beyond day 10 should be dispensable. Alternatively, if mianserin still influenced mortality later in life, shorter exposures would lead to a shorter lifespan extension compared to a lifelong exposure. We therefore limited mianserin exposure to 8 hr, 1, 5, 10 and 15 days and compared their lifespan with animals treated for the entire life (Figure 6c,d). Exposing the animals for 5 or 10 days was sufficient to extend lifespan to the same extent as lifelong exposure (Figure 6c,d). Shorter exposures (8 hr, 1 day) also extended lifespan, but not by as much, showing that removing mianserin from the culture is an effective means to restrict its action (Figure 6c,d). Taken together, these results are most consistent with a model in which mianserin specifically lowers the rate of age-associated change during the first few days of adulthood, thereby extending their longevity (Figure 6e) and postponing the onset of mortality. While the change in age-associated mortality rate during early adulthood is too small to be accuratly determined, when we measured drift-variance, it allowed us to monitor the age-associated change in the transcriptome during early adulthood (Figure 6e,f).
Since the effect of mianserin in early adulthood overlapped with the reproductive period (first 5 days of adulthood), we asked whether mianserin treatment increased reproductive lifespan as has been observed in tph-1(mg280) mutants (Sze et al., 2000). Mianserin treatment blocks serotonin-induced egg-laying (Petrascheck et al., 2007), but had a minor effect on amount or timing of spontaneous egg-laying and brood size (Figure 6g). Most importantly, mianserin did not increase reproductive longevity (Figure 6g).
We further considered the possibility that mianserin acted by a mechanism similar to lifespan extension by germline ablation (Figure 6h). Two previous findings suggested otherwise: i) Lifespan extension by germline ablation depends on daf-16, while mianserin does not (Arantes-Oliveira et al., 2002; Petrascheck et al., 2007); ii) germline ablation increases lifespan of eat-2(ad1116) mutants while mianserin does not (Crawford et al., 2007). We measured whether mianserin treatment mimicked the increased proteasome activity observed in glp-1 mutants (Vilchez et al., 2012) (Figure 6h). A 24 hr mianserin treatment did not increase the proteasome activity, as measured by a fluorescence-based assay for chymotrypsin-like activity. On day 5, mianserin slightly decreased proteasome activity, consistent with a slight increase in drift-variance in proteasome-related genes (Figure 6h; Figure 6—figure supplement 1c). We concluded that mianserin specifically lowers the rate of age-associated change in somatic tissues and does not involve a mechanism directly related to the germline.
Transcriptional drift-variance increases with age in mice and humans
Our data demonstrate that changes in drift-variance provide a metric for aging that correlates with mortality in C. elegans. To test whether drift-variance also increases with age in mammals, we re-analyzed published gene expression data-sets obtained from aging mouse tissues, aging human brains, and from fibroblasts derived from Hutchinson-Gilford progeria syndrome patients (Figure 7) (Lu et al., 2004; Liu et al., 2011; Jonker et al., 2013). We calculated drift-variances from brain, kidney, liver, lung, and spleen based on gene expression data-sets from mice aged 13, 26, 52, 78, 104 and 130 weeks. We calculated drift-variances using 13-week-old mice as a young reference (see Methods) and pooled mice into age-bins of 30, 60 and 100 weeks to reduce variability. Drift-variance increased in all tissues with age (Figure 7a). Compared to the drift-variance changes observed in C. elegans (Figure 2a), these changes however were small.
Because the 13-week-old mice were used as reference for young age (see methods), the drift-variance in the 30-week-old group including the 13-week-old sample is artificially low (Figure 6a, see material and methods). To better reflect the actual variance of the 30-week-old group, we set aside the data of one 13-week-old mouse to use as a young reference and recalculated drift-variances for all samples (Figure 7b). This strategy has the advantage that we can observe the real drift-variance for the 30-week-old group by excluding the reference data-set, but has the disadvantage that the results are less robust as they all depend on a single reference sample. Plotting drift-variance for each organ as a function of age confirmed that as mice age, drift-variance increases in all organs (Figure 7b). It will be interesting to learn if the different rates by which drift-variance increases in different organs will also be observed in other data-sets.
We re-analyzed the data from Lu et al. that recorded gene expression profiles from 32 human brains aged 26 to 106 years of age (frontal cortex) (Figure 7c) (Lu et al., 2004). For the first plot, we binned the data into 20-year bins and calculated the overall drift-variance for each 20-year bin. As a young age reference, we used the mean gene expression of adults below 30 (26, 26, 27, and 29) (see Materials and methods). This analysis shows that over the entire population, drift-variance remains relatively stable until the age of sixty, and then starts to rise (Figure 6c). We also plotted the drift-variance of each individual as a function of age. This revealed a significant correlation (Spearman, rho=0.6, P=0.0014) between age and drift-variance in the human brain.
Irrespective of the age of the mother, the aging process starts afresh for each new generation. We therefore hypothesized that aging must be reversed with each new generation and asked whether it is possible to reverse increases in drift-variances. To address this question, we re-analyzed the data-set generated by Liu et al. who derived induced pluripotent stem cells (iPSCs) from fibroblasts of healthy controls (BJ) and patients suffering from Hutchinson-Gilford progeria syndrome (HPGS), an accelerated aging syndrome (Figure 7e) (Liu et al., 2011). As a young-reference to calculate drift-variance, we used human H9 embryonic stem cells (ESC). As expected for a premature aging syndrome, fibroblasts from HGPS patients showed increased drift-variance relative to BJ control fibroblasts (Figure 7e). Furthermore, nuclear reprogramming reduced drift-variance in iPSCs to levels closer to those seen in H9 embryonic stem cells. Thus, increases in drift-variance are reversed by nuclear reprogramming in vitro
Discussion
In this study, we have analyzed the dynamics of aging C. elegans transcriptomes and how these dynamics are affected by mianserin treatment. We separate transcriptional changes across groups into those that characterize activation or inhibition of entire pathways (type I) and those that characterize the relative expression levels among genes (type II, transcriptional drift, Figure 1h,i). In C. elegans, transcriptional drift continuously increases with age across the transcriptome, substantially altering stoichiometric balances observed in young animals (Figure 2a). Longevity mechanisms induced by either pharmacologically blocking serotonergic signaling or by blocking insulin signaling by daf-2 RNAi attenuate transcriptional drift (Figure 2a,g). Abolishing lifespan extension by these mechanisms by either blocking serotonergic signaling too late (mianserin, day 5) or by addition of daf-16 RNAi (daf-2) abolished the attenuation of drift-variance (Figure 2).
Detailed analysis of redox-related pathways showed that mianserin-reduced drift-variances are associated with improved stress resistance in older age (Figure 3). Mutations in the serotonin receptor SER-5 that abolish the effect of mianserin on drift-variance also abolished its effect on stress resistance and lifespan (Figure 4, 5).
Using transcriptome-wide drift-variance values as a metric for age showed that mianserin treatment attenuated the age-associated increase of drift-variance, thereby preserving the characteristics of a much younger (~3 days-old) transcriptome up to chronological day 10 (Figure 2a, 6a). These results showed that mianserin caused a 7–8 days delay in age-associated transcriptional change and suggested that the physiological changes leading to a lifespan extension were already completed by day 10.
Measuring mortality levels supported this conclusion. By day 12, the entire mortality curve was shifted parallel by 7–8 days (Figure 6b) showing that the physiological delay leading to a lifespan extension was already completed. Experiments in which animals were exposed to mianserin for limited periods of time confirmed that mianserin exposure for the first 5–10 days of adulthood was necessary and sufficient to fully extend lifespan (Figure 6c,d). The most parsimonious explanation that accounts for all these results is that mianserin treatment slows degenerative processes specifically between day 1 and 10, extending the duration of the period of young adulthood thereby postponing the onset of major mortality around mid-life (Figure 6e,f).
Biological interpretation of transcriptional drift-variance
Aging has been shown to cause DNA damage, degeneration of the nuclear architecture, loss of histones, loss of histone modification (Kaeberlein et al., 1999; Scaffidi and Misteli, 2006; Burgess et al., 2012). These changes contribute to the degenerative phenotypes observed with aging (Mostoslavsky et al., 2006; Feser et al., 2010; Peleg et al., 2010). In the present study, we used expression patterns of young adults as a reference to monitor the aging process across the transcriptome. We found that aging causes the expression of genes within functional groups to drift apart, causing a loss of co-expression patterns as observed in young adults. We quantified this phenomenon using drift-variance, defined as the variance in gene expression among genes. It is important to distinguish transcriptional noise, which measures the variance of the same genes among samples (Bahar et al., 2006), from transcriptional drift, which measures variance among genes within the same samples. At present it is unclear whether transcriptional drift is the consequence of a regulated program or of degenerative changes in the nucleus that lead to a loss of transcriptional control. Consistent with a regulated program are recent findings that the germline actively represses the activation of heat shock promoters via histone methylation, causing a decline in heat shock capacity (Labbadia and Morimoto, 2015b). Consistent with degenerative changes are recent findings that show the loss of histone methylation to cause aberrant gene expression that increases with age leading to a transcriptional drift-like effect (Somel et al., 2006; Mercken et al., 2013; Pu et al., 2015; Sen et al., 2015).
Irrespective of whether transcriptional drift is the consequence of a regulated program or a degenerative change, its effect on pathway function is likely to be detrimental. Many physiological processes depend on appropriate stoichiometry of their components. Large and persistent deviations in mRNA balance as measured by drift-variance are likely to result in stoichiometric imbalances in protein complexes, negatively affecting proteostasis as has been recently observed (Houtkooper et al., 2013; Walther et al., 2015). Our results modulating drift-variance for redox genes via mianserin and SER-5 certainly suggest that the age-associated increases in drift-variance are associated with regulatory decline (Figures 3, 5). Attenuation of transcriptional drift in the redox system was associated with an improved homoestatic capacity, i.e. an improved ability of the redox system to appropriately respond to outward stimuli.
Transcriptional drift also provided a useful concept to analyze aging transcriptomes. Accounting for its effects dramatically simplified what was an initially excessively complex expression pattern (Figure 1). Excluding gene expression changes due to drift left a set of genes that changed expression in response to mianserin treatment that was enriched for genes related to stress, innate immunity, aging and the xenobiotic response. With the exception of the xenobiotic response, which is expected to be triggered by addition of a foreign substance such as mianserin (Figure 2f), all other functions have been linked to serotonin signaling (Table 1) (Zahn et al., 2006; Petrascheck et al., 2007; Rangaraju et al., 2015a).
Further, in accordance with the hypothesis that increases in drift-variance are a signature of aging in the transcriptome, we find that drift-variance is attenuated by two longevity mechanisms (mianserin and daf-2 RNAi) across large sections of the transcriptome. Many of the age-associated changes that were reversed by mianserin were also reversed by daf-2 RNAi (58%). This overlap is consistent with chemical epistasis experiments. Treating daf-2(e1370) mutants with mianserin causes only a partial extension of lifespan (11% instead of 31%) (Petrascheck et al., 2007) consistent with the idea that many of the genes attenuated by mianserin treatment are already attenuated in daf-2(e1370) mutants and thus do not further contribute to a lifespan extension. It should be noted that age-associated increases in drift-variance do not contradict the idea that transcription factors regulate longevity. Activation of DAF-16 target genes by daf-2RNAi prevent age-associated drift of thousands of genes, thus resulting in a net decrease of drift, even though a transcriptional program has been induced (Figure 2g). Our experiment did not address the questions whether increasing drift-variance beyond what occurs naturally with age accelerates aging and whether attenuation of transcriptional drift-variance is universal to all longevity mechanisms.
At this point, it is prudent to mention possible pitfalls associated with transcriptional drift analysis. Drift-variance calculations require data-sets that include multiple ages (3 or more) as direct statistical comparisons to the young-reference are not permissible. Furthermore, in the context of GO annotations, it is important to realize that if a given GO annotation contains significant numbers of mis-annotated genes, these genes may change expression in a different direction giving the erroneous impression of transcriptional drift. To account for these effects in our study, we i) used the experimentally determined oxidative stress signature derived from Olivera et al (Figure 3e), and ii) used a robust Levene’s test to determine statistical differences. The robust Levene’s test uses a 10% trimmed mean, which removes large outliers such as those that would be expected by mis-annotation. These safeguards, however, are only effective if the number of mis-annotated genes is small relative to the total number of genes.
Conceptually, transcriptional drift is not a biomarker for aging. It is a metric for aging similar to lifespan measurements that can be used to monitor age-associated physiological changes on the molecular level within groups of genes. Lifespan measurements record the fraction of organisms alive in different cohorts at any given time to compare rates of aging, while drift-variance allows a similar comparison based on transcriptional drift-variance. What made drift-variance measures essential for the present study was that it allowed us to monitor age-associated physiological changes in young animals, at a time when age-associated mortality levels are too low to be accurately determined (see below).
Period-specific lifespan extension
Measuring lifespan of mianserin-treated and untreated C. elegans revealed a mean lifespan extension of 7–8 days (Figure 2). Lifespan measurements detect differences after the majority of the animals have died and make no statements about the period during which the relevant physiological events that lead to an increase in lifespan occur (Figure 2c,e) (Mair et al., 2003; Partridge and Gems, 2007). The finding that transcriptional drift values in mianserin-treated animals already showed a 7–8 day delay in physiological change as early as day 10 suggested a model in which the physiological events responsible for the 7–8 days lifespan extension take place (and conclude) prior to day 10 (Figure 2a, 6a,e).
Determining mortality levels at different ages confirmed this model. Mianserin or water is added on day 1 of adulthood to the same preparation of N2 animals. The mortality levels of both cohorts (water, mianserin) are therefore identical at the start of the experiment. Thus, the lower mortality level observed on day 12 in mianserin-treated animals is the result of a lower mortality rate prior to day 12 (Figure 6b). Furthermore, mianserin ceases to affect mortality rates past day 12 as evident by highly parallel mortality curves (Figure 6b). As with the results obtained with drift measurements, the most plausible explanation is that mianserin treatment specifically decelerates the rise in mortality in young adults leading to a lower mortality level sometime between day 10 to day 12 that persists throughout life, ultimately revealing itself in a 7–8 day lifespan extension (~30–40% increase in lifespan) (Figure 6b).
Analysis of drift-variance, PCA, mortality and survivorship independently arrive at the same 7–8 days delay in physiology, either measured as a feature of transcriptomes or by recording death times. All methods suggest that the delay is completed before day 10 or 12 and therefore occurs during early adulthood. We further experimentally confirmed this suggestion by showing that treatment for the first five or ten days of life was necessary and sufficient to achieve the same lifespan extension observed with lifelong treatment (Figure 6c,d).
Even though this period exactly overlaps with the reproductive period, the effect of mianserin appears to be specific to somatic tissue (Figure 6g,h). In contrast to germline ablation, mianserin extends lifespan of daf-16 mutants but not of eat-2 mutants (Crawford et al., 2007; Petrascheck et al., 2007; Vilchez et al., 2012) and does not increase proteasome activity as observed in glp-1 mutants (Figure 6h). It is still possible that the mianserin-induced lifespan extension interacts or depends on the germline, but if it does, the connection is more indirect potentially similar to what has been observed for dietary restriction (Crawford et al., 2007).
Lifespan extension mechanisms that decelerate the rate of mortality are generally interpreted as slowing the aging process, while a parallel shift as the one we observe with mianserin is interpreted as a constant risk factor that causes a proportional shift in the overall risk of death (Mair et al., 2003; Harrison et al., 2009; Vaupel, 2010; Kirkwood, 2015). Our data do not challenge any of these prior interpretations, but add a further possibility. Parallel shifts may also be brought about by a period extension in which the rate of age-associated physiological change is specifically lowered in young adults. Age-associated mortality in young adults is very low compared to extrinsic mortality factors and thus changes in age-associated mortality rates are difficult to reliably determine (Partridge and Gems, 2007; Beltran-Sancheza et al., 2012). Specific changes in mortality rates during early adulthood therefore can go unnoticed but manifest themselves later as parallel shifts at the time when age-associated mortality levels are sufficiently high to be reliably determined. Whether the attenuation of physiological changes specific to young adults that affects later mortality, as seen for mianserin, is the equivalent of slowing aging in young adults is a debate for the general aging community.
In summary, this work describes the phenomenon of transcriptional drift and how it can be used as a metric for aging. Using this metric, we show that blocking serotonergic signals by mianserin delays age-associated physiological changes such as transcriptional drift and mortality exclusively during early adulthood, thus extending the duration of this period and postponing the onset of age-associated mortality.
Materials and methods
Measurement of transcriptional drift and drift-variance
View detailed protocolAnalyzing the RNA-seq data in aging C. elegans, we observed dramatic changes in the transcriptome with age. We simply termed these changes ‘transcriptional drift’, to emphasize the ambiguity of these changes. These changes could either be the result of regulated changes as part of a biological program, or caused by a progressive loss of transcriptional control with age. Note that a progressive loss of transcriptional control does not necessarily have to result in random changes. A gene that is continuously activated in young animals may be less activated in older animals due to a progressive functional decline in the transcriptional machinery. Thus, a gradual loss of transcriptional control would cause an age-associated decline in expression of that gene in a non-random fashion. Conversely, repressive chromatin is lost with age leading to increases in transcription that are repressed in young animals. As most physiological processes depend at least to some degree on transcriptional regulation, we propose that expression changes of genes within the same pathway that go into opposing directions (drift-variance increases) are detrimental for the functionality of the pathway (as seen for redox pathways in Figure 3b). These changes may also allow us to indirectly track the functional decline by measuring transcriptional drift.
Calculating transcriptional drift and drift-variance
Request a detailed protocolTranscriptional drift (td) is the change in transcript level of a gene at a given age from its level in young animals (“young reference”). As all the subsequent calculations depend on the age chosen for “young reference” we made sure to indicate the age used as a “young reference” for each plot (see below). For all the C. elegans work, the “young reference” age was day 1, at the onset of reproductive maturity in adulthood.
For any gene x, transcriptional drift (td) is defined as (Equation 1).
or, which is the same as
where, ‘cpm’ stands for counts per million; ‘t’ stands for time in days, weeks or years, dependent on the organism.
Equation 1 normalizes the level of transcription for all genes to 0 for a young animal. Note: If several biological replicates are available for the age of the young reference, a variance for the young age can be calculated (see the section below titled ‘Variance for “the young reference”’).
To evaluate changes in co-expression, we calculated the drift-variance (dv) (Equation 3) over a group of n genes with transcriptional drift-values ranging from tdi=1 to tdn.
Thus, if genes maintain a youthful co-expression pattern, drift-variance stays relatively small. If large fractions of genes within a GO or an entire transcriptome change expression in opposing directions, the drift-variance increases, suggesting a loss of youthful co-expression patterns as shown in Figure 1h,i.
Variance for the “young reference”
Request a detailed protocolIf multiple replicate data-sets for the “young reference” age are available, it is possible to plot drift-variance for the young reference as well. There are two ways to incorporate multiple “young reference” data-sets, each of which has its advantages or disadvantages.
Method #1 uses all “young reference” samples to calculate a mean gene expression level for each individual gene to generate the “young reference” values for Equation 1. Method #1 will result in a drift-variance for the “young reference” age as well, but this drift-variance is too small and should not be used for statistical comparisons due to circular referencing. The advantage of method #1 is that the results for all subsequent ages are more robust as the inclusion of several “young reference” samples thereby reducing the overall noise (used in Figures 2a,g, 3b, 7a,c,e).
Method #2 allows calculating a real drift-variance value for young animals by setting aside one or several samples as the “young reference.” These samples are only used as references and therefore do not contribute to the drift-variance in each plot. For the remaining experimental replicates of the same age, transcriptional drift is then calculated using Equation 1 without including any of the “young reference” samples.” This will result in a drift-variance greater than 0 for the youngest age and show how much drift varies between young animals. Method #2 has the disadvantage that if there are only few young reference samples are available, and only one is used as a young reference, all values of the graph depend on a single reference sample. We used this method #2 to calculate the variances for Figure 7b,d. The case of 7d was ideal as there were 4 samples less than 30 years of age which were set aside as reference and that allowed us to calculate the “young reference”-mean over all 4 samples. As drift-variances for these 4 samples are artificially low due to self referencing they were excluded from the plot. Ideally, an experiment would have 4–6 gene expression replicates for the “young reference” age, in which case, half of them could be used as references, the others as experimental samples.
How transcriptional drift and variance relate to measures like fold-changes in transcription is shown in Supplementary Figure 2a–d. To determine whether the differences in variance were statistically different, we used the Brown-Forsythe version of the Levene’s test, as implemented in STATA software.
Calculations for drift-plots in Figures
Request a detailed protocolFigure: 1g: Volcano plot used mean cpm values from all three biological replicates.
The 0 line (young reference, day 1 expression, yellow line) indicates the expected expression level for young day 1 adult animals.
Black: Each dot represents one of for the 3,367 genes that significantly change expression with age between day 1 and day 3. The -log10(P-value) of the P-value comparing day3 water vs day 1 water is shown as a function of the the log2(cmp day 3 water / cpm day1 water).
Blue: Same 3,367 genes as above. However the -log10(P-value) comparing day3 mianserin vs day 1water is shown as a function of the the log2(cmps day 3 mianserin / cpm day1 water). Note: both data-sets (black and blue) use identical y- coordinates to demonstrate the reduction in age-associated changes upon mianserin-treatment. (cpm stands for: counts per million).
Young Reference: To obtain a ‘young reference’ value for each individual gene the mean expression level across all three biological replicates of young day 1 old water-treated C. elegans animals was calculated.
Figure 1h, i: Drift plots for genes involved in oxidative phosphorylation (KEGG pathway: cel 00190) and the lysosome (KEGG pathway: cel 04142). Only one out of three replicates was used to generate these plots. Transcriptional drift for oxidative phosphorylation and lysosomal genes (line graphs) was calculated using Equation 1 and plotted as a function of C. elegans age (gray lines). At each age, the transcriptional drift-variance across all genes within the pathway was calculated using Equation 2 and plotted as Tukey-style box plots omitting outliers. Tukey plots were superimposed over the line graphs. See Equation 1, 3. Outliers were only omitted for graphical purposes but not for statistical testing (robust Levene’s test). The lines for each gene were included in these two plots, superimposed on the Tukey-style box plot to illustrate the significance and utility of the box plots in visualizing transcriptional drift.
Young reference: As a “young reference” value for each individual gene, the expression level of young day 1 old water-treated C. elegans animals was used. Only replicate #1 of our data-set was used.
Figure 2a: Drift plots for all 19,196 genes in our data-set of water-treated control and mianserin-treated animals. Tukey plots show drift-variance calculated for the entire transcriptome (Equation 3). See Equation 1, 3. Outliers were only omitted for graphical purposes, but not for statistical testing (robust Levene’s test).
Young reference: To obtain a “young reference” value for each individual gene, the mean expression level across all three biological replicates of young day 1 old water-treated C. elegans animals was calculated.
Figure 2b: Drift plots show transcriptional drift on day 5 for 19,196 genes as a function of mianserin concentration. For each concentration, drift-variances were calculated for 5-day-old animals that were treated with increasing concentrations of mianserin on day 1, and plotted as Tukey-style box plots as a function of mianserin concentrations, excluding outliers. Outliers were only removed for graphical purposes but not for statistical testing (robust Levene’s test).
Young reference: To obtain a “young reference” value for each individual gene, the mean expression level across all three biological replicates of young day 1 old water-treated C. elegans animals was calculated.
Figure 2d: Drift plots show transcriptional drift on day 10 of adulthood for 19,196 genes as a function of age when mianserin-treatment was started. Tukey plots show drift-variance calculated for the entire transcriptome on day 10 (Equation 3) as a function of age at which mianserin-treatment was initiated.
Young reference: To obtain a “young reference” value for each individual gene, the mean expression levels across all three biological replicates of young day 1 old water-treated C. elegans animals was calculated.
Figure 2f: Log2 fold changes in expression for each gene shown in the y-axis were calculated by the formula: y = log2(cpm treatment day 10/cpm water day 1).
Figure 2g: The data from Murphy et al. were dowloaded from the Princeton Puma database. Expression values were calculated using the following variables in the data-set: expression value = ch1netmean/ch2normalizednetmean. Drift plots for control- RNAi, daf-2(RNAi) treated and daf-16(RNAi); daf-2(RNAi) treated animals were plotted as transcriptional drift-variance as a function of C. elegans age. To plot drift-variance for the entire transcriptome as function of age in days, we binned the data as follows. Day 0 (8 hr), day 1 (24 hr), day 2 (28 hr, 40 hr, 52 hr), day 4 (72 hr, 96 hr), day 6 (144 hr, 196 hr).
Young reference: As a “young reference” value for each individual gene we used the expression level at 8 hr of age. The young reference was determined for each RNAi treatement specifically (control RNAi, daf-16(RNAi); daf-2(RNAi), daf-2(RNAi).
Figure 3e: The log fold gene expression with age was calculated for each of the 252 genes that are known to be upregulated in response to oxidative stress and for each of the 88 genes known to be downregulated in response to oxidative stress. We then performed a linear fit for each set of genes for water-treated (gray) and mianserin-treated (blue) samples. Shaded region shows the 95% confidence interval.
Figure 7a, b: 7a) Drift plots showing transcriptional drift and drift-variance in different tissues across different mouse ages. For each age, the drift-variance was calculated across the entire transcriptome (Equation 3) and plotted as Tukey-style box plots omitting outliers. As only three mice were available for each age, we pooled two ages for each age bin.
7b) Drift-variance for each tissue as a function of age.
Young Reference: 7a: To obtain a “young reference” value for each individual gene, the mean expression level across all three biological replicates of young 13-week-old mice was calculated for each tissue.
Young Reference 7b: To obtain “young reference” values for each individual gene, we used one single 13-week-old replicate as a “young reference” from each tissue. The data from the “young reference” did not contribute to the graph and thus show a real transcriptional drift-variance.
Figure 7c, d: 7c). Drift plots showing transcriptional drift-variance in human gene expression data from frontal cortices as a function of age. For 7c, the data were pooled into 20 year bins.
7d) Plots drift-variance calculated based on Equation 3 as a function of age for each sample individually.
Young Reference: To obtain “young reference” values for each individual gene, the mean gene expression levels was calculated averaging expression levels from 4 samples aged 25 to 29 years and used as the “young reference” value in Equation 1.
Figure 2—figure supplement 1: e) The transcriptional drift plots were constructed by using the GEO data-sets GSE21784 and GSE46051, which are independent publicly available data-sets for aging C. elegans.
f) The transcriptional drift plots were constructed by sub-sampling the data from our RNA-seq. We randomly assigned half of all genes (out of 19,196) to one of 10 gene-sets each containing ~1000 genes (5%) and plotted the drift-variance for each set. All 10 sets look nearly indistinguishable to Figure 2a.
Figure 2—figure supplement 2: f) The drift plot was constructed by removing all the genes from our data-set that were not detected in the sterile CF512 strain, thereby removing genes likely resulting from eggs and germline.
g) The drift plot was constructed by removing all genes from our data-set that were detected by RNA-seq in isolated C. elegans eggs.
k) Gene-sets enriched in AFD neurons (left plot), ASE neurons (middle plot) and NSM neurons (right plot) were used to construct drift plots based on their expression in our data-set.
Principle component analysis
Request a detailed protocolPrincipal components analysis plot (Figure 6a) was generated from the counts table using multidimensional scaling as implemented by the plotMDS function in the edgeR package, which computes inter-sample distances as the root-mean-square of the 500 genes with the largest log2 fold-changes between each pair of sample (the 'leading log fold-change").
Chemicals
Solvents used to prepare stock solutions: Paraquat was dissolved in water; mianserin was dissolved either in water or DMSO as mentioned; Mirtazapine, Dihydroergotamine, LY-165,163/PAPP, Mirtazapine, Metergoline, Ketanserin, Methiothepin, and Amperozide were dissolved in DMSO; FUDR was dissolved in S-complete (Table 9).
Strains
Detailed descriptions of all strains used in this study are tabulated below. All strains were backcrossed at least 4 times with the N2 Bristol strain. All strains were maintained as described in (Brenner, 1974). The strains with name starting with VV were generated by outcrossing to N2 Bristol strain in our lab (Table 10).
Lifespan assay and analysis
Request a detailed protocolLifespan assays were conducted in 96-well plates as described in (Solis and Petrascheck, 2011; Rangaraju et al., 2015b). Briefly, age-synchronized animals were cultured in S-complete media containing E. coli OP50 as feeding bacteria (~2 × 109 bacteria mL−1) in 96-well plates, such that 5–15 worms are in each well. At the L4 stage, FUDR was added to prevent animals from producing offspring. Solvent (water or DMSO) or small molecules were added on day 1 of adulthood, exposing the worms to control or compound treatment until the end of the assay. When used, DMSO was kept to a final concentration of 0.33% v/v. Live animals were scored visually, based on movement induced by shaking and application of light to each well. Animals were scored three times a week, until 95% of animals were dead in all the tested conditions. Statistical analysis was performed using the Mantel–Haenszel version of the log-rank test.
Stress resistance assays
Request a detailed protocolResistance to oxidative stress was determined by measuring survival of mianserin-treated and untreated worms after a 24 hr exposure to the ROS-generator paraquat (Methyl viologen). Experimental worm cultures were set up as described in Lifespan assays. For dose response assays, paraquat was added to a final concentration of 0, 25, 50, 75, 100 mM on day 5 of adulthood. For paraquat time-course experiment (Figure 3c), paraquat was added 3 days, 5 days, or 10 days after addition of mianserin on day 1 of adulthood. For mianserin time-course experiment (Figure 3d), 50 µM mianserin was added on day 1, day 3 or 5 of adulthood, followed by 100 mM paraquat on day 10. For all experiments, survival of worms was assessed 24 hr after paraquat addition and expressed as the percentage of live versus total animals.
RNA-sequencing (RNA-seq) transcriptional studies and data analysis
Request a detailed protocolMianserin-induced changes in transcription were determined by RNA-seq. A total of 12 conditions were tested each run in three biological replicates. N2 worms were cultured in 96-well plates as described in (Solis and Petrascheck, 2011). Animals in cohort #1 were treated on day 1 with water (solvent) or 50 µM mianserin, and harvested on day 3, 5, and 10 of adulthood. Animals in cohort #2 were treated with water (solvent control) or mianserin (2, 10, or 50 µM) on day 1 of adulthood and harvested on day 5. Animals in cohort #3 were treated with water (solvent) or 50 µM mianserin on day 1, day 3 and day 5 and harvested on day 10 (See Figure 1a). RNA was also harvested from untreated day 1 adults, to obtain the “young reference”. Harvested animals were washed three times in ice cold Dulbecco’s phosphate buffer saline and frozen in liquid nitrogen. A parallel lifespan assay was conducted for all cohorts to ensure mianserin action. Three biological replicates were harvested for every cohort. To extract RNA, frozen worms were re-suspended in ice-cold Trizol, zirconium beads, and glass beads (cat # 03961-1-103 and cat # 03961-1-104) in the ratio of 5:1:1 respectively, and disrupted in Precellys lysing system (6500 rpm, 3 x 10 s cycles) followed by chloroform extraction. For RNA-seq, the extracted RNA was precipitated and purified further using Qiagen RNAeasy Mini kit columns (cat # 74104). RNA was precipitated using isopropanol and washed once with 75% ethanol. Integrity of the RNA was confirmed with a Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). To prepare the library, 100 ng of total RNA per sample was processed using NuGEN Encore Complete DR RNA-seq Prep Kit (NuGEN; San Carlos; CA, USA), as per manufacturer’s instructions. The libraries were sequenced using v2 sequencing chemistry in a HiSeq2000 platform (Illumina, San Diego, CA, USA). A single-read sequencing approach was used with 100 cycles, resulting in reads with a length of 100 nucleotides each. Libraries containing their own index sequences were sequenced in a multiplex manner by pooling six libraries per lane. Resulting sequences were obtained after 20–30 million reads per sample. Sequence data were extracted in FASTQ format and used for data analysis.
RNA-seq data analysis
Request a detailed protocolRNA-seq data were analyzed by aligning the reads to the C. elegans reference genome and transcriptome from WormBase using Tophat 2 (Kim et al., 2013), and unambiguously mapped reads were counted for each annotated gene in each sample (Lawrence et al., 2013). Data were normalized for sequencing depths (counts per million, cpm) but not for gene length as no comparisons between genes within the same sample were made. The quasi-likelihood F-test from the edgeR package (Robinson and Oshlack, 2010; Lund et al., 2012) was used to test these counts for statistically significant differential gene expression between water- and mianserin-treated samples, while controlling for expression differences between the 3 biological replicates. We performed multiple testing correction by using the Benjamini-Hochberg procedure to compute a false discovery rate (FDR) value for each gene, and we considered an FDR less than 10% to be significant (Benjamini and Hochberg, 1995; Zhang et al., 2009).
Quantitative real-time PCR (qRT-PCR) and data analysis
Request a detailed protocolAll qRT-PCR experiments were conducted according to the MIQE guidelines (Bustin et al., 2009), except that samples were not tested in a bio-analyzer, but photometrically quantified using a Nanodrop. All strains were cultured in 96-well plates as described in (Solis and Petrascheck, 2011). Water (solvent) or mianserin were added on day 1 of adulthood and worms were harvested on day 5. RNA was extracted as described above, followed by DNAse (Sigma, cat # AMPD1-1KT) treatment and reverse transcription using iScript RT-Supermix (BIO-RAD, cat # 170–8841) at 42ºC for 30 min. Quantitative PCR reactions were set up in 384-well plates (BIO-RAD, cat # HSP3901), which included 2.5 µl Bio-Rad SsoAdvanced SYBR Green Supermix (cat # 172–5264) or Kapa SYBR Fast master mix (cat # KK4602), 1 µl cDNA template (2.5 ng/µl, to final of 0.5 ng/µl in 5 µl PCR reaction), 1 µl water, and 0.5 µl of forward and reverse primers (150 nM final concentration for BIO-RAD SYBR mix and 75 nM final for Kapa SYBR mix) (see Table below for oligo pairs used for qRT-PCR of genes tested). Quantitative PCR was carried out using a BIO-RAD CFX384 Real-Time thermocycler (95ºC, 3 min; 40 cycles of 95ºC 10 s, 60ºC 30 s; Melting curve: 95ºC 5 s, 60ºC- 95ºC at 0.5ºC increment, 10 s). Gene expression was normalized to three reference genes, rcq-5, crn-3 and rpl-6, using the BIO-RAD CFX Manager software. Statistical significance was determined using Student’s t-test (Table 11).
Measurement of 26S proteasome activity
Request a detailed protocolWild-type N2 worms were cultured as described (Solis and Petrascheck, 2011). Water or Mianserin 50 µM were added on day 1 and 26S proteasome activity was assayed on day 2 and day 5 using the Millipore Proteasome activity kit (cat# APT280), following manufacturer’s protocol. Equal number of worms per condition were washed off culture media using ice cold Dulbecco’s phosphate buffer saline and freshly lysed using Precellys system (6500 rpm, 3 x 10 s cycles) in assay buffer (25 mM HEPES, pH 7.5, 0.5mM EDTA, 0.05% NP-40, and 0.001% SDS (w/v)). Chymotrypsin-like proteasome activity in the lysates were assessed using the Suc-LLVY-AMC substrate and fluorogenic AMC substrate cleavage was measured in 20 min intervals for 120 min. A subset of lysates were pre-incubated with Lactacystin (12.5 µM final) to ensure specificity of AMC cleavage by 26S proteasome. The amount of cleaved AMC fragments were quantified using TECAN xfluor safire II system at excitation of 360 nm and emission of 480 nm. The resulting readings were normalized to the total protein content in the samples measured using Bradford assay.
Mortality curve and probability of detection
Request a detailed protocolMortality curves were generated based on the life table provided in Figure 6—figure supplement 1, tabulating death times of 15 independent experiments performed over 5 years. Each experiment consisted of 2 cohorts (water or 50 µM mianserin) and each cohort consisted of ~100 worms each amounting to ~1500 worms per condition. Power of detection was determined by Monte-Carlo simulations using a parametric model with parameters derived from our survival data of a cohort of over 5,026 N2 animals. The power of detection plot (Figure 6—figure supplement 1) shows the probability to detect a true lifespan extension with a significance level α=0.01 as a function of percent increase in lifespan for an experiment consisting of n animals. An accuracy of 1 day is the equivalent of a 5% increase in lifespan.
Data availability
-
Suppression of Transcriptional Drift Extends C. elegans Lifespan by Postponing the Onset of MortalityPublicly available at the NCBI Gene Expression Omnibus (Accession no: GSE63528).
References
-
Gene ontology: tool for the unification of biologyNature Genetics 25:25–29.https://doi.org/10.1038/75556
-
Uncoupling lifespan and healthspan in caenorhabditis elegans longevity mutantsProceedings of the National Academy of Sciences 112:E277–E286.https://doi.org/10.1073/pnas.1412192112
-
Early life events can shape aging and longevityCurrent Aging Science 8:11–13.https://doi.org/10.2174/1874609808666150422123355
-
Early cohort mortality predicts the rate of aging in the cohort: a historical analysisJournal of Developmental Origins of Health and Disease 3:380–386.https://doi.org/10.1017/S2040174412000281
-
Controlling the false discovery rate - a practical and powerful approach to multiple testingJournal of the Royal Statistical Society Series B-Methodological 57:289–300.
-
DNA damage, chromatin, and transcription: the trinity of agingCurrent Opinion in Cell Biology 24:724–730.https://doi.org/10.1016/j.ceb.2012.07.005
-
The metabolite alpha-ketoglutarate extends lifespan by inhibiting ATP synthase and TORNature 510:397–401.
-
Ageing and inflammation – a central role for mitochondria in brain health and diseaseAgeing Research Reviews 21:30–42.https://doi.org/10.1016/j.arr.2015.02.001
-
The molecular signature and cis-regulatory architecture of a c. elegans gustatory neuronGenes & Development 21:1653–1674.https://doi.org/10.1101/gad.1560107
-
Elevated histone expression promotes life span extensionMolecular Cell 39:724–735.https://doi.org/10.1016/j.molcel.2010.08.015
-
A systematic review of the serotonergic effects of mirtazapine in humans: implications for its dual action statusHuman Psychopharmacology: Clinical and Experimental 21:117–125.https://doi.org/10.1002/hup.750
-
Accelerated epigenetic aging in down syndromeAging Cell 14:491–495.https://doi.org/10.1111/acel.12325
-
Measurements of age-related changes of physiological processes that predict lifespan of caenorhabditis elegansProceedings of the National Academy of Sciences 101:8084–8089.https://doi.org/10.1073/pnas.0400848101
-
Deciphering death: a commentary on gompertz (1825) 'on the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies'Philosophical Transactions of the Royal Society B: Biological Sciences 370:20140379.https://doi.org/10.1098/rstb.2014.0379
-
Points of significance: visualizing samples with box plotsNature Methods 11:119–120.https://doi.org/10.1038/nmeth.2813
-
Effects of oxidative stress on behavior, physiology, and the redox thiol proteome of caenorhabditis elegansAntioxidants & Redox Signaling 14:1023–1037.https://doi.org/10.1089/ars.2010.3203
-
The biology of proteostasis in aging and diseaseAnnual Review of Biochemistry 84:435–464.https://doi.org/10.1146/annurev-biochem-060614-033955
-
Software for computing and annotating genomic rangesPLoS Computational Biology 9:e1003118.https://doi.org/10.1371/journal.pcbi.1003118
-
Detecting differential expression in RNA-sequence data using quasi-likelihood with shrunken dispersion estimatesStatistical Applications in Genetics and Molecular Biology 11:.https://doi.org/10.1515/1544-6115.1826
-
The PANTHER database of protein families, subfamilies, functions and pathwaysNucleic Acids Research 33:D284–D288.https://doi.org/10.1093/nar/gki078
-
High-throughput small-molecule screening in caenorhabditis elegansMethods in Molecular Biology 1263:139–155.https://doi.org/10.1007/978-1-4939-2269-7_11
-
H3K36 methylation promotes longevity by enhancing transcriptional fidelityGenes & Development 29:1362–1376.https://doi.org/10.1101/gad.263707.115
-
Measuring caenorhabditis elegans life span in 96 well microtiter platesJournal of Visualized Experiments.https://doi.org/10.3791/2496
-
Gene expression becomes heterogeneous with ageCurrent Biology 16:R359–R360.https://doi.org/10.1016/j.cub.2006.04.024
-
The embryonic cell lineage of the nematode caenorhabditis elegansDevelopmental Biology 100:64–119.https://doi.org/10.1016/0012-1606(83)90201-4
-
Aging as an event of proteostasis collapseCold Spring Harbor Perspectives in Biology 3:a004440.https://doi.org/10.1101/cshperspect.a004440
-
Systemic stress signalling: understanding the cell non-autonomous control of proteostasisNature Reviews Molecular Cell Biology 15:211–217.https://doi.org/10.1038/nrm3752
Article and author information
Author details
Funding
NIH Office of the Director (New Innovator Award)
- Michael Petrascheck
Ellison Medical Foundation (New Scholar in Aging Award)
- Michael Petrascheck
Muscular Dystrophy Association (Development Grant/ Postdoctoral Fellowship)
- Sunitha Rangaraju
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was funded by grants to MP, from the NIH (DP2 OD008398), a grant from The Ellison Medical foundation (AG-NS-0928-12), an MDA Development Grant for SR, and an NSF GRFP Fellowship for GMS. SEE is supported by The Ellison Medical Foundation (AG-NS-0950-1), and by a Baxter Foundation Young Faculty Award. Some strains were provided by Shigen-Japan or the CGC, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440). We thank Jim Priess (U. Washington), and Bruce Bowerman (U. Oregon) for advice, and Dr. Veena Prahlad, Dr. Eros Lazzerini Denchi, Dr. Maria Carretero, Dr. Bruno Conti, Dr. Andrew Chisholm and Caroline Broaddus for critical reading of the manuscript.
Copyright
© 2015, Rangaraju et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 7,691
- views
-
- 1,267
- downloads
-
- 82
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
B-cell repertoires are characterized by a diverse set of receptors of distinct specificities generated through two processes of somatic diversification: V(D)J recombination and somatic hypermutations. B-cell clonal families stem from the same V(D)J recombination event, but differ in their hypermutations. Clonal families identification is key to understanding B-cell repertoire function, evolution, and dynamics. We present HILARy (high-precision inference of lineages in antibody repertoires), an efficient, fast, and precise method to identify clonal families from single- or paired-chain repertoire sequencing datasets. HILARy combines probabilistic models that capture the receptor generation and selection statistics with adapted clustering methods to achieve consistently high inference accuracy. It automatically leverages the phylogenetic signal of shared mutations in difficult repertoire subsets. Exploiting the high sensitivity of the method, we find the statistics of evolutionary properties such as the site frequency spectrum and dN/dS ratio do not depend on the junction length. We also identify a broad range of selection pressures spanning two orders of magnitude.