|Home | About | Journals | Submit | Contact Us | Français|
Aging is under genetic control in C. elegans but the mechanisms of lifespan regulation are not completely known. MicroRNAs (miRNAs) regulate various aspects of development and metabolism and one miRNA has been previously implicated in lifespan.
Here we show that multiple miRNAs change expression in C. elegans aging, including novel miRNAs, and that mutations in several of the most up-regulated miRNAs lead to lifespan defects. Some act to promote normal lifespan and stress resistance while others inhibit these phenomena. We find that these miRNAs genetically interact with genes in the DNA damage checkpoint response pathway and in the insulin signaling pathway.
Our findings reveal that miRNAs both positively and negatively influence lifespan. Since several miRNAs up-regulated during aging regulate genes in conserved pathways of aging and thereby influence lifespan in C. elegans, we propose that miRNAs may play important roles in stress response and aging of more complex organisms.
In nature, longevity is characterized by a remarkable degree of variability, but the factors that dictate this plasticity of aging are largely unknown. Research pioneered in C. elegans has identified multiple genetic mechanisms that play a significant role during organismal aging . The accumulated evidence demonstrates that aging is influenced by a complex interplay of regulatory mechanisms that respond to stress, nutrient availability and environment. Mutations in the daf-2 insulin/IGF receptor gene cause an extension of lifespan and stress resistance by deregulating the function of the DAF-16/FOXO transcription factor . Similarly, mutations in genes that modulate the response to caloric restriction and DNA repair have significant effects on longevity, through mechanisms independent of the insulin pathway . Moreover, these pathways are evolutionarily conserved in higher organisms. Nevertheless, the complete genetic circuits by which these genes modulate aging are not fully elucidated.
MicroRNAs (miRNAs), discovered in C. elegans , are endogenous, noncoding RNAs with important roles in controlling gene expression and development in plants and animals. Diverse functions have been attributed to miRNAs - including roles in developmental timing, neuronal development, metabolism, immunity, lifespan and cancer  - and target predictions suggest that they may regulate thousands of human genes . Still, only a small number of the thousands of known miRNAs have been implicated in a specific biological function, while the role of the vast majority of miRNAs remains unknown. Mutations in the developmental timing genes, lin-4 (a founding member of the miRNA family) and its target lin-14 significantly alter the normal aging process of C. elegans, thus revealing a degree of miRNA control of aging in C. elegans . In addition, microarray analysis in C. elegans revealed distinct miRNA expression changes during aging . These observations show that miRNAs may function in pathways that impact lifespan but little is known about the extent of their roles and their mechanisms.
To identify miRNAs with functions in lifespan, we undertook a deep sequencing survey of miRNAs from aged tissue in wild-type (N2) and in a long-lived daf-2(e1370) mutant C. elegans. We hypothesized that aging-specific miRNAs might exhibit differential expression in older animals or in mutants with aberrant lifespan, such as daf-2(e1370), which live twice as long as wild-type C. elegans . Total RNA was collected from N2 and daf-2 animals when they reached young adulthood (Day 0) and when animals reached Day 10, an inflection point when wild-type animals begin dying. We obtained a total of 889,762 sequences of small RNAs (15-25nt) from the 4 samples of C. elegans examined, of which 626,616 sequences (70.4%) aligned perfectly to the C. elegans genome (Table S1). Overall, our survey detected perfect matches to the mature forms of 120 out of the 155 miRNAs annotated in miRBase 13.0, which indicates the significant sequencing depth of our study. This allowed us to identify miRNAs, including novel species, that are enriched in aging animals.
We hypothesized that yet-unidentified miRNAs may contribute to gene regulatory networks that affect lifespan. After discarding previously annotated ncRNAs and mRNAs from the 626,616 C. elegans aligned reads, we obtained 83,439 non-redundant, unannotated sequences from C. elegans. These 83,439 sequences, representing 142,689 total reads, were tested to determine if they represent novel C. elegans miRNAs. To identify likely miRNA candidates we used the miRNA-discovery software package, miRDeep, which identifies novel miRNAs based on the stability of their putative pre-miRNA hairpins in addition to "signature reads" - sequences found with the deep sequencing data that are consistent with miRNA biogenesis . Our analysis revealed a set of 11 sequences that are predicted to fold into hairpin structures characteristic of bona-fide miRNAs (Fig. 1a) and represent likely novel miRNAs. Moreover, several of these novel miRNA candidates (mir-A, mir-B, mir-C and mir-I) share homology over their "seed" sequence (nucleotides 2–8 of the mature miRNA) with miRNAs in higher eukaryotes, including 3 human miRNAs, suggesting that they may be family members of these miRNAs (Fig. 1b). Interestingly, miR-I appears to be a new member of the miR-34 miRNA family, and miR-34 is one of the miRNAs that exhibited the greatest increase in expression with aging (see below) (Fig. S2 and Table 1).
Several of these novel candidates – such as miR-N, miR-L, miR-M, miR-G, miR-H, miR-I and miR-K - are expressed preferentially in aged C. elegans (Fig. S1). In addition, these miRNAs show reduced expression in daf-2(lf) worms in both young and old animals (Fig. S1). Thus, these novel miRNAs seem to be preferentially expressed in normal, aged worms – a result that is consistent with the rationale of the experiment. Finally, four miRNAs - miR-A, miR-B, miR-C and miR-I are expressed only in the daf-2 samples, and may represent miRNAs that are normally under the regulation of the IIS pathway.
In order to validate the expression of our candidate novel miRNAs, we designed custom TaqMan assays specific to the putative mature sequence of the four candidate miRNAs for which we had the fewest reads - miR-G, miR-H, miR-I, and miR-K. Our qRT-PCR results showed that indeed all four candidates were expressed in aged C. elegans, although very much less frequently, as expected, than previously characterized miRNAs, such as let-7 (Fig. 1c). Importantly, the expression of these tested candidate miRNAs was significantly reduced in mutant C. elegans animals with attenuated expression of the miRNA-associated factor ALG-1 (Fig. 1c). As alg-1 has been shown to be important for mature miRNA accumulation , our results are consistent with the classification of mir-G, mir-H, mir-I and mir-K as novel miRNAs.
In an effort to discover miRNAs with functions specific to post-developmental processes, we measured the expression changes of all miRNAs in our sequencing data over the course of aging and in long-lived mutants. We observed a relative overall reduction in the expression of known miRNAs between days 0 and 10 of adulthood in both wild-type and daf-2 animals. Whereas nearly half of all sequencing reads in Day 0 wild-type adult animals represent known miRNAs (47.9%), by Day 10 the frequency of known miRNAs dropped to 33.6% of all genome-matching reads (Table S1). This pattern was similar in the long-lived daf-2 mutants (45.1% at Day 0 to 38.1% by Day 10). We found that although the expression levels of many miRNAs remain unchanged, a subset of miRNAs exhibit considerable changes in expression with aging (Fig. S2). For purposes of comparing expression changes of individual miRNAs, here we considered only those 81 miRNAs for which we obtained more than 10 sequences reads between the 2 samples being compared. Using this criteria, 8.6% miRNAs increased and 28.4% decreased in expression by more than 2-fold from Day 0 to Day 10 of adulthood in wild-type N2 animals (Fig. S2, Table 1). daf-2 animals exhibited a nearly identical pattern: 9.8% increased while 25.6% miRNAs decreased (Fig. S2, Table 1).
miR-246 was the most up-regulated miRNA during aging, showing nearly a 6-fold increase in expression from Day 0 to Day 10 in wild-type. Other miRNAs that exhibited strong increases in expression during aging in wild-type include miR-71, which became the second most abundant miRNA by Day 10 of aging (>4-fold up-regulation), miR-34 (4-fold up-regulation), miR-253 (3.3-fold up-regulation), and the homologous miRNAs miR-238 (2.9-fold up-regulation) and miR-239a/b (>2-fold up-regulation) (Fig. S2 and Table 1). Conversely, let-7 was the most down-regulated miRNA in aging - exhibiting a 7.4-fold decrease in expression from day 0 to day 10 of adulthood in N2 animals (Table 1). Other miRNAs that decreased in expression with aging include miR-41 (>6-fold down-regulation), miR-70 (>5-fold down-regulation) and miR-252 (>5-fold down-regulation). Notably, two miRNAs with dynamic changes in expression with aging, let-7 and mir-34, are known to be involved in cancer [10, 11]. In a further link between miRNAs and lifespan, a recent report has implicated let-7 in aging-induced senescence in mouse neuronal stem cells .
We validated our measurements of the expression of several miRNAs by using several established measures of RNA expression. We tested five miRNAs that we cloned more frequently in old animals (miR-246, miR-238, miR-71, miR-240 and miR-67) and three that we cloned less frequently in old animals (let-7, miR-70, miR-237), and found that the changes in expression with aging detected by TaqMan qRT-PCR correlated with the expression changes we observed by cloning for both N2 (r=0.90) and daf-2 animals (r=0.90) (Fig. S3). We also used northern blot to compare the expression changes with aging of representative miRNAs, such as miR-246, the most highly up-regulated miRNA in our cloning survey, and, as a control, miR-66, whose expression remains constant. Consistent with our cloning results, we found that in both wild-type and daf-2 mutants that miR-66's expression remains unchanged while miR-246's expression increases significantly with aging (Fig. S3b).
In order to characterize the spatio-temporal expression pattern of these miRNAs during aging, we examined GFP expression driven by promoters of aging-associated miRNAs (Fig. 2). In confirmation of our deep-sequencing and qRT-PCR findings, we observe a significant increase in expression of GFP driven by the promoters of mir-71, mir-238, mir-239 and mir-246 during aging (Fig. 2). Interestingly, both pmir-71::GFP and pmir-239::GFP reach maximum GFP expression between days 5–7 of adulthood, which subsequently wanes, whereas pmir-238::GFP and pmir-246::GFP expression increases steadily until the end of life. As reported before, miR-71 and miR-238 are expressed ubiquitously and miR-246 is expressed in gonadal sheath cells . During aging, pmir-71::GFP is noticeably increased throughout the intestine and in the pharynx (Fig. 2). By contrast, pmir-239::GFP expression is seen primarily in intestine and in neurons (Fig. 2). Although the GFP expression patterns in transgenic animals may not fully reflect the endogenous expression patterns of these miRNAs (particularly as expression of genes introduced by microinjection is often silenced in the germline), the observation of increased expression of these miRNAs in intestine, neurons and somatic gonad is particularly interesting as signaling in these tissues has been previously shown to regulate lifespan in C. elegans .
The insulin/insulin-like growth factor-1 (IGF-1) receptor, DAF-2, is a key component of the insulin/insulin-like IGF-1 signaling (IIS) pathway in C. elegans that responds to environmental signals to regulate stress response and lifespan. Previous studies have identified a wide variety of mRNAs whose expression is altered in daf-2 mutant animals and found that a significant number of these genes could affect lifespan when mutated . Similarly, we hypothesized that miRNAs that exhibit significantly altered expression in daf-2 animals may be good candidates for novel lifespan determining genes. In an effort to identify miRNAs that depend on the IGF-1/daf-2 pathway, we compared the miRNA profiles of wild-type animals with those in daf-2(e1370) mutants.
We observe that the vast majority of miRNAs exhibit very similar changes in expression with aging between N2 and daf-2 animals (r=0.83) suggesting that the changes in miRNA expression that we observed are reproducible across strains and that the expression of the majority of miRNAs is not affected by a loss of DAF-2 (Fig. S2). Although the profile of miRNA expression in N2 and daf-2 animals was similar, we found a few specific miRNAs that exhibit altered expression in daf-2 mutant animals (Fig. S2, Table 1). Notably, several of these miRNAs also exhibited altered expression with aging (see below), suggesting that they may be particularly good candidates for genes that influence lifespan. Interestingly, the miRNA with the greatest increase in expression in young daf-2 animals, miR-237, is one of the most down-regulated miRNAs in older animals (Table 1). In fact, other miRNAs exhibit this pattern - miR-62 and miR-252 are among the most up-regulated miRNAs in daf-2 mutant animals, when compared to wild-type of the same age (Day 10), in addition to being among the most down-regulated miRNAs with aging. Conversely, miR-239b which appears to be down-regulated in young daf-2 animal is also among the most up-regulated microRNAs with aging (Table 1). Collectively, these daf-2-associated miRNAs reveal the intersection of the daf-2 pathway with miRNA-regulated genes that impact lifespan (Table 1).
We considered the hypothesis that miRNAs displaying dramatic changes in expression across aging and in daf-2 mutants may function in lifespan control. Specifically, we focused on miRNAs that are over-expressed with aging since the up-regulation of specific miRNAs appears to be a common response of various organisms to stress conditions such as DNA damage, hypoxia and nutrient deprivation [10, 15–18]. In order to uncover putative functions for these miRNAs in lifespan, we examined mutants that contain deletions for the miRNAs of interest (Table 1) from the C. elegans Genome Consortium (CGC). These mutant C. elegans were backcrossed 3–6 times to the control laboratory wild-type C. elegans strain (N2) before assaying lifespan. Knock-out mutants for all miRNAs that exhibit greater than 2-fold increases in expression with aging were tested: mir-246(n4636), mir-71(n4115), mir-34(gk437), mir-253(nDf64), mir-238(n4112), and mir-239(nDf62). In addition, we also tested the lifespans of mir-240(n4541) and mir-67(n4899) - knock-out mutants of miRNAs that exhibited more modest increases of expression with aging (>1.5 fold). Lifespan assays on these mutants revealed that four of the miRNAs that exhibit the largest expression changes with aging - miR-71, miR-238, miR-239, and miR-246 - also influence lifespan (Fig. 3). Specifically, deletions of miR-71, miR-238 and miR-246 decreased lifespan, while deletion of miR-239 lead to reproducible and significant lifespan extension as compared to wild-type animals (Fig. 3a). The mean lifespan of these mutants was altered by ~10–50% in relation to wild-type N2 animals (Table S3). Other miRNA mutants tested - mir-34(gk437 and n4276), mir-253(nDf64), mir-240(n4541) and mir-67(n4899) - did not exhibit significantly altered lifespan (data not shown). Although the effects observed with our aging-associated miRNA mutants are not as dramatic as the aberrant lifespan observed in such mutants as daf-2 and age-1 (which typically double the lifespan of C. elegans ), they are within the range of other well characterized effectors of C. elegans lifespan such as the 5–20% extension induced by a clk-1 mutation [19, 20], the 12–20% extension caused by sir-2.1 overexpression , the 14% lifespan extension induced by resveratrol treatment , and the significance cut-offs (>5–10% extension) used in two large-scale RNAi screens [23, 24].
Since miR-71, miR-238 and miR-246 were highly up-regulated during aging (Table 1), these results suggest that these miRNAs increase during aging and promote longevity. An alternative explanation for the shorter lifespans of mir-71, mir-238 and mir-246 mutants might lie in non-specific sickness caused by deletion of these genes. However, none of these mutants exhibited any physiological or behavioral abnormalities during growth and development ; Fig. S4 and Table S4). In assays that tested morphology, growth, development, locomotion, pharyngeal pumping, defecation, egg-laying, and dauer formation these mutants all exhibited wild-type behavior [25, 26]. Nevertheless, and in order to test directly for a role in lifespan, we generated, via microinjection, transgenic C. elegans that over-express the aging-associated miRNAs miR-71, miR-238, miR-239 and miR-246 via extrachromosomal arrays . Strikingly, we observe that over-expression of both miR-71 and miR-246 significantly and reproducibly increase longevity while miR-239 overexpression reduces lifespan (Fig. 3 b–d). Overexpression of miR-238 did not affect lifespan, suggesting that either this miRNA is not a bona fide aging gene or that its over-expression may cause confounding effects such as developmental abnormalities (Fig. S4a). Importantly, the lifespan phenotypes caused by miR-71, miR-239 or miR-246 overexpression are the opposite of what we observe when these miRNAs are deleted (Fig. 3a). Together, these observations strongly suggest a direct role for miR-71 and miR-246 to promote longevity and a function for miR-239 in pathways that antagonize longevity in C. elegans.
Since stress resistance is associated with lifespan , we considered the possibility that deletions of aging-associated miRNAs may also impact the response of C. elegans to stress. To that end, we tested the heat stress and oxidative stress response of C. elegans mutant adults bearing deletions of mir-71, mir-238, mir-239 and mir-246. Compared to the wild-type N2 strain, we observed that the long-lived mutant mir-239(nDf62) exhibited significant resistance to both heat stress (35°C) and oxidative stress (paraquat) (Fig. 4). Conversely, the short-lived mutant mir-71(n4115 exhibited a significantly increased sensitivity to both heat and oxidative stress while mir-238(n4112) and mir-246(n4636) caused increased sensitivity to oxidative and heat stress, respectively (Fig. 4). These results demonstrate that miR-71, miR-238, miR-239 and miR-246 function in pathways that regulate the response of C. elegans to conditions of environmental stress, in support of our observations that these miRNAs are necessary for normal lifespan in C. elegans.
In order to understand the biological pathways that mediate the function of aging-associated miRNAs, we tested their function in the IGF-1/insulin-like pathway, which is known to mediate longevity in C. elegans, Drosophila and mouse. We tested if the aging phenotypes of mir-71(n4115), mir-238(n4112), mir-239(nDf62) and mir-246(n4636) would be affected in daf-2 and daf-16 loss-of-function environments. We found that RNAi directed against daf-2 and daf-16 modified the lifespan phenotypes of mir-71(n4115) and mir-239(nDf62), suggesting that these miRNAs function through the insulin-like pathway in C. elegans. Indeed, daf-16(RNAi) completely abolished the long lifespan of mir-239(nDf62) mutants, showing that the long-lived phenotype of mir-239 mutants depends on the presence of daf-16 (Fig. 5a). In addition, mir-239(nDf62) mutant animals did not further enhance lifespan in a daf-2(RNAi) background, suggesting that mir-239 and daf-2 function in the same genetic pathway (Fig. 5a). Similarly, we observed that mir-71(n4115) animals grown on daf-2(RNAi) were shorter lived than wild-type animals grown on daf-2(RNAi), which demonstrates that loss of miR-71 partially suppresses the long lifespan induced by daf-2(RNAi) (Fig. 5b). In addition, daf-16(RNAi) did not further decrease the lifespan of mir-71(n4115) animals, suggesting that mir-71 and daf-16 function in the same pathway (Fig. 5c). By contrast, we found that the reduction in lifespan of mir-246(n4636) mutants in daf-2(RNAi) was comparable to the effect of loss of miR-246 in an empty vector RNAi (~15% reduction in lifespan), suggesting that miR-246 does not function through daf-2 (Fig. S6a). In addition, we found that miR-238 is not necessary for longevity induced by loss of DAF-2 as mir-238(n4112) did not affect the lifespan of daf-2(RNAi) animals suggesting that miR-238 functions upstream of daf-2 or through an independent pathway (Fig. S6b). Together, these observations strongly suggest that miR-71 and miR-239 function through the IGF-1/insulin-like pathway to mediate lifespan in C. elegans.
Since mir-71(n4115) did not fully suppress the long lifespan of animals on daf-2(RNAi), we considered whether miR-71 might also function in other pathways of aging. To that end, we tested whether the absence of a functional copy of miR-71 might affect the lifespan of mutants in the DNA damage checkpoint pathway, which has been shown to affect lifespan in C. elegans , We found that the increased lifespan observed when RNAi is directed against two checkpoint proteins – CDC-25.1 and CHK-1 – was completely suppressed in a mir-71(n4115) genetic background, suggesting that the function of CDC-25.1 and CHK-1 on lifespan depends on miR-71 (Fig. 5d and Fig. S6c).
In order to understand the molecular basis of miRNA function in IIS and in the DNA damage checkpoint pathways, we examined the 3′ UTRs of genes in these pathways and searched for putative binding sites to the aging-associated miRNAs. We found that several genes in the IGF-1/insulin like pathway such as AGE-1, PDK-1, AKT-1 and DAF-16 have putative binding sites for miR-71 (Table S6) as predicted by mirWIP . In addition, CHK-1, CDC-25.1 and CDC-25.2 are also predicted to be targets of miR-71 (Table S6). Indeed, miR-71 is predicted to have 4 binding sites in the 3’UTR of PDK-1 and 3 binding sites in CDC-25.1 (Fig. 6a). Since CDC-25.1 and PDK-1 antagonize lifespan, we reasoned that miR-71 may promote longevity by down-regulating these genes during adulthood. In order to test this hypothesis, we assayed the expression levels of these genes by qRT-PCR in a mir-71(n4115) background and found that the mRNA levels of both PDK-1 and CDC-25.1 in Day 10 adults were significantly elevated in mutants that lack miR-71, as compared to wild-type N2 animals (Fig. 6c). Interestingly, in young adults (Day 0), the levels of PDK-1 and CDC-25.1 were unchanged in a mir-71(n4115) background as compared to N2 animals suggesting that miR-71 functions to down-regulate these genes in older animals specifically (Fig. 6c). We also examined the expression levels of other predicted targets of miR-71 – AGE-1, DAF-16 and CHK-1 – in a mir-71(n4115) background and we found that the levels of these mRNAs were unchanged in Day 0 adults and modestly increased (< 2-fold) in Day 10 mir-71(n4115) adults, compared to wild-type worms.
Since we found that miR-239 function depends on DAF-16, we also searched for potential targets of miR-239 in the IIS pathway. Interestingly, however, none of the core genes in the IIS pathway are predicted to be targeted by miR-239. When we examined the expression of genes in the IIS pathway by qRT-PCR, we found a decrease in the mRNA levels of AGE-1 and PDK-1 in Day 10 mir-239(nDf62) animals (2.2-fold and 2.1-fold, respectively) (Fig. S6d), suggesting that miR-239 may act upstream of AGE-1/PDK-1 through a factor that is yet to be identified (Fig. 6d).
Our results reveal a greater involvement of miRNAs in C. elegans lifespan than previously appreciated, with some promoting and others suppressing longevity. We observed significant expression changes of multiple miRNAs during aging and also identified novel miRNAs in aged C. elegans. Analysis of the lifespan of strains of C. elegans containing mutations in miRNAs that exhibit the most significant changes of expression during aging demonstrated that three of the miRNAs that are most over-expressed in aged animals - miR-71, miR-238 and miR-246 - act to increase longevity in C. elegans, while a fourth miRNA - miR-239 - normally limits lifespan. We found, by genetic criteria, that both miR-71 and miR-239 function through the insulin-signaling pathway. In addition, we observed that miR-71 also interacts with the DNA damage response pathway, as it is required for the longevity extension caused by knock-down of CDC-25.1 and CHK-1. Given that miR-71 is predicted to target CDC-25.1 and CHK-1, this suggests a possible negative feedback loop, whereupon miR-71’s ability to promote longevity is antagonized by CDC-25.1, but miR-71, in turn, negatively regulates these genes (Fig. 6d). Consistent with this idea, we observed that in mutant animals that lack miR-71, the expression of CDC-25.1 in older animals is increased 8-fold. The fact that miR-71 functions in both the IIS and the DNA checkpoint pathways suggests that miR-71 is a possible link between these pathways (Fig. 6d).
Notably, the lifespan abnormalities reported here are the first observable phenotypes reported for these miRNAs. Indeed, an earlier survey failed to identify any phenotypes during growth and development for the vast majority of miRNA knock-out mutants in C. elegans . Given that the majority of miRNAs individually seem to not be essential for C. elegans development, our observation of dramatic effects on adult lifespan for four independent loss-of-function mutants is striking. These results suggest that miR-71, miR-239 and miR-246 may be aging-specific genes - functioning specifically during adulthood to regulate genetic pathways that promote lifespan. Consistent with this idea, we observed that the up-regulation in mRNA levels of two predicted targets of miR-71 was specific to older animals. The absence of any phenotypes beside abnormal lifespans for these miRNA mutants further suggests that they might not be genes with antagonistic pleiotrophic function - i.e. beneficial early in life but detrimental post-reproduction. However, we cannot rule out a subtle function for these miRNAs in development - indeed, these miRNAs are expressed dynamically during development, especially miR-71 [31, 32]. It is possible that the trade-offs inherent to the roles of these miRNAs may only be apparent in specific environmental conditions, such as nutrient limitation and non-optimal temperature. In this case, these miRNAs may be used to fine-tune the response of C. elegans to environmental stress. In support of this idea, we observe that the miRNA mutants that we found to exhibit abnormal lifespan also show correspondingly abnormal responses to stress. These results are in agreement with previous studies that demonstrate that miRNAs may be particularly well suited to modulate organismal response to stress . Since miRNAs have the potential to regulate multiple targets, increasing the level of one miRNA could potentially affect the action of a vast number of target genes. Indeed, the up-regulation of specific miRNAs has been noted in response to a variety of sources of stress in different organisms - DNA damage and oncogenic stress to mammalian cell lines, [10, 16], hypoxia in neoplastic cells , and nutrient deprivation in plants . In that context, it is particularly interesting that we observe aberrant lifespan and stress responses in C. elegans mutants of miRNAs that are up-regulated with aging. These observations lend credence to the idea that miRNAs may be particularly useful modulators of gene regulatory pathways that respond to stress to affect an organism's lifespan. Since miRNAs exhibit high degrees of conservation of sequence and function it is likely that miRNAs also play extensive roles in aging in more complex animals.
The wild-type C. elegans strain N2 (Bristol) was cultured using standard protocols . Mutant strains were obtained from the Caenorhabditis Genetics Center or generated in house and are listed in Supplemental Experimental Methods.
Lifespan assays were performed at 20°C as previously described . Statistical analysis was done using Graphpad Prism 5.0 and JMP software to determine survival difference (log-rank Mantel-Cox test). Summary of lifespan statistics is shown in Tables S3–S4.
Heat stress and oxidative stress (Paraquat) experiments were performed as previously described . At least 30 animals were tested for each strain. Mean and standard deviation were determined from experiments done in triplicate. p-values were calculated using Student's t-test.
Total RNA was harvested from worm pellets, size-selected (15–25nt) and miRNAs were cloned as described before [35, 31]. The resulting cDNAs were then concatemerized, as described , and sequenced by 454 deep sequencing (454 Lifesciences / Roche) .
A total of 889,762 reads corresponding to small RNAs (15–30nt) were obtained by 454 deep sequencing. The deep-sequencing reads were aligned to the C. elegans genome (version WS190, obtained from Wormbase) and the number of sequence reads that correspond to known miRNAs was assessed using perfect sequence matching to a database of known miRNAs (miRBase, v.13.0). To compare the differential expression of miRNAs across samples, the number of miRNA reads in each sample was normalized to the total number of reads in each sample that matched the C. elegans genome (Table S1).
To identify novel miRNAs we used the software package miRDeep . Beginning with all sequence reads that match the C. elegans genome perfectly, we discarded sequences that match previously annotated ncRNAs and mRNAs, using NCBI blastall. In addition, matches to known C. elegans miRNAs (miRBase v.13) and 21U-RNAs [37, 38, 32] were removed by perfect sequence matching. The remaining reads were then mined for putative novel miRNAs using miRDeep, using default settings , and secondary structures of putative pre-miRNA hairpins were generated using RNAFold .
We used the program mirWIP to identify targets of aging-associated miRNAs . mirWIP targets are enriched for characteristics of bona-fide miRNA targets, including the structural accessibility of target sequences, the total free energy of miRNA-target hybridization and the topology of base-pairing to the 5' seed region of the miRNA . As mirWIP’s target predictions originate from immuno-precipitation of the RISC factors AIN-1 and AIN-2 in C. elegans, they are particularly appropriate to our study, and are reported to have a lower false-positive rate than other methods . mirWIP predictions were obtained from http://126.96.36.199/query.php; (predictions for mir-238 were not available for mirWIP).
Raw sequencing reads from 454 deep sequencing and processed small RNA sequences were deposited in NCBI's Gene Expression Omnibus [GSE24510]. Additional experimental details available in Supplemental Experimental Methods.
We thank C. elegans Genetic Center for supplying strains. ADL was supported by a National Research Service Award Post-Doctoral Fellowship from the NIH (1F32AG030851). ZP was supported by a Jane Coffin Childs Post-Doctoral Fellowship. FJS was supported by a Breakthrough in Gerontology grant from the American Federation for Aging Research, the Ellison Medical Foundation and the NIH (AG033921).
Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.