Submit an Article
Online ISSN: 1945-4589
Impact Journals, LLC. is the publisher of Aging: www.impactjournals.com.
Impact Journals is a member of the Wellcome Trust List of Compliant Publishers.
Impact Journals is a member of the Society for Scholarly Publishing.
Learn about our FREE Post-Publication Promotion Services
Research Paper Volume 6, Issue 10 pp 835—851

PHA-4/FOXA-regulated microRNA feed forward loops during Caenorhabditis elegans dietary restriction

Awadhesh Pandit1, , Vaibhav Jain , Neeraj Kumar , Arnab Mukhopadhyay ,

  • 1 Molecular Aging Laboratory, National Institute of Immunology, Aruna Asaf Ali Marg, New Delhi 110067, India

Received: September 1, 2014       Accepted: October 27, 2014       Published: October 31, 2014      

https://doi.org/10.18632/aging.100697
How to Cite

Copyright: © 2014 Pandit et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Dietary restriction (DR) increases life span and delays the onset of age-related diseases across species. However, the molecular mechanisms have remained relatively unexplored in terms of gene regulation. In C. elegans, a popular model for aging studies, the FOXA transcription factor PHA-4 is a robust genetic regulator of DR, although little is known about how it regulates gene expression. We profiled the transcriptome and miRNAome of an eat-2 mutant, a genetic surrogate of DR, by Next Generation sequencing and find that most of the miRNAs are upregulated in the young-adult worms, none significantly downregulated. Interestingly, PHA-4 can potentially regulate the expression of most of these miRNA genes. Remarkably, many of the PHA-4-regulated genes that are induced during DR are also targets of the PHA-4-upregulated miRNAs, forming a large feed-forward gene regulatory network. The genes targeted by the feed-forward loops (FFLs) are enriched for functions related to ubiquitin-mediated decay, lysosomal autophagy, cellular signalling, protein folding etc., processes that play critical roles in DR and longevity. Together our data provides a framework for understanding the complex and unique regulatory network employed during DR, suggesting that PHA-4 employs such FFLs to fine-tune gene expression and instil robustness in the system during energy crisis.

Introduction

Dietary restriction (DR) is the only intervention that can consistently increase life span and delay onset of age-related diseases in most of the model systems tested [1]. Since the initial discovery in 1935 [2] that animals fed less food lived substantially longer, DR has been shown to increase life span in many vertebrates [1, 3]. In several non-vertebrates that are important models for ageing research, like yeast, Drosophila and Caenorhabditis elegans, DR increases life and health span. In mammals, DR significantly protects against age-related disorders like type II diabetes, atherosclerosis and Alzheimer's disease [1]. Although a large body of literature exists about the physiological consequences of DR, the molecular mechanisms by which DR increases longevity are still unclear [4].

Invertebrate models like C. elegans have provided tremendous insight into the mechanisms of aging and longevity assurance in general, and DR in particular. DR is achieved through multiple genetic and extrinsic interventions in worms [5-8]. Mutations in an acetylcholine receptor gene, eat-2 lead to decreased pharyngeal pumping resulting in lowered food intake [6]. The eat-2 model of DR remains the most-studied genetic surrogate used to decipher the mechanisms associated with the process, although DR-like phenotypes can be obtained by other gene perturbations [9, 10]. DR can also be initiated by diluting the bacterial food that the worms ingest [5, 11]. Interestingly, life span extension on DR brought about by both genetic and non-genetic interventions depends upon the function of a conserved FOXA transcription factor, PHA-4 [11]. The molecular mechanisms explaining how PHA-4 mediates DR-induced longevity are currently less clear.

MicroRNAs (miRNAs) have emerged as important and conserved gene expression modulators across phylogeny [12]. Primary miRNAs are transcribed in the nucleus and processed by the Drosha complex to produce precursor miRNA (pre-miRNA, ~70 nucleotides long) [13, 14]. These pre-miRNAs are exported out of the nucleus, by the Exportin complex, into the cytoplasm where they are processed by the Dicer enzyme to generate mature miRNA [13, 14]. The mature miRNA associates with specific Argonaute proteins leading either to the degradation of target mRNA or attenuation of translation, depending on the homology of its seed region with the 3′-UTR of the mRNA [13-16]. MiRNA-mRNA interactions have been proposed to shield protein expression from various oscillations in transcription during development, thereby installing robustness in the system that prevents phenotypic fluctuations- a process called canalization [17-22]. Although contributions of miRNA have been documented in many biological processes in C. elegans and other organisms, including aging [23-28], little is known about their role in dietary restriction.

In an effort to study the transcriptional and post-transcriptional response to DR in details, we profiled the mRNA and miRNA of eat-2(ad1116) by Next Generation Sequencing (NGS). Here we show that many of the miRNAs (82 out of 368) are upregulated in the day 1 young-adult mutant worms. Remarkably none are significantly downregulated. This pattern was unique to eat-2 mutation-mediated DR as it was significantly different from the profile of low insulin signalling mutant, daf-2(e1370). Using published ChIP-seq data for PHA-4 [29, 30] and qRT-PCR validation, we show that 65 of the 82 miRNAs are potential direct transcriptional targets of PHA-4. The transcriptomics followed by bioinformatics analysis and qRT-PCR validation showed that PHA-4 may regulate a large portion of the transcripts that are upregulated in eat-2(ad1116). Interestingly, the PHA-4-regulated miRNAs target a major portion of the PHA-4-regulated transcripts during DR, thereby forming transcription factor incoherent feed forward loops (FFLs). The FOXO transcription factor DAF-16 does not appear to extensively utilize such FFLs during conditions of low insulin signalling. We hypothesize that these FFLs modulate the expression of PHA-4-dependent genes during energy-depleted conditions, building robustness into the system that reduce phenotypic fluctuations. These genes are involved in ubiquitin-mediated protein degradation, lysosomal autophagy, protein folding, signalling as well as metabolism, processes that need to be strongly regulated during energy crisis. Thus, our study elucidates the complexity of gene regulation following initiation of DR in eat-2(ad1116) and defines the central role of FOXA/PHA-4 in this process, justifying its position as a robust genetic regulator of DR-induced longevity.

Results

Next generation sequencing reveals eat-2(ad1116) miRNA profile

In order to gain an insight into the complex interplay of miRNAs during DR, we performed NGS analysis to compare miRNA profiles of wild-type (WT) N2 Bristol and eat-2(ad1116). The eat-2(ad1116) strain lives considerably longer than WT; however, none of the strains have significant mortality at day 1 (equivalent to young-adults) or day 8 of adulthood (aging worms) (Figure S1). The eat-2 mutants have defective pharyngeal pumping and consume less food since hatching. We hypothesized that the system in eat-2(ad1116) is already geared up for long life and would have significantly modified post-transcriptional cellular environment on day 1 of adulthood i.e., at the young-adult stage. Additionally, we wanted to study the changes that occur later in their life, as compared to WT, in order to understand the late-life cellular response to DR. Synchronized worms were harvested on day 1 or 8 of adulthood and small RNA sequencing libraries prepared (see experimental procedures). Each sample was run in a single lane of a flow cell in the Genome Analyzer IIx (GAIIx) platform (Illumina Inc., USA) that led to the generation of around 20 million reads per lane with an average trimmed read size of ~22 nucleotides (Figure S2A). All reads mapped to the C. elegans genome while 81-95% readily mapped to the miRBase release 19 that has a collection of 368 worm miRNAs (Figure S2B, 1). We detected 185 miRNA in day 1/young-adult worms and 187 in day 8/aging WT samples, while 224 and 161 miRNA were detected in eat-2(ad1116) samples in day 1 and day 8 samples, respectively, all miRNA having greater than 10 read counts (Table S1). We were also able to detect between 0.5-1.2 million reads that mapped to known 21U-RNAs (Table S1). After eliminating rRNA, tRNA and ncRNA reads, the remaining un-annotated reads were processed for novel miRNA discovery, as discussed below.

We compared the miRNA profiles of WT worms collected at days 1 or 8 with a previously published study on changes in miRNA profile during aging in WT nematode [31] (Table S2). Out of a total of 21 miRNA reported to be upregulated in day 8 compared to day 1 in their study, 16 matched to our data. We found an additional 20 miRNA upregulated in day 8 samples. When we compared the downregulated miRNAs from the two studies, 7 out of 18 matched our data while we obtained an additional 24 miRNAs, probably because of increased depth of sequencing. It may also be due to the less number of miRNAs (174) present in the miRbase release 14, to which the earlier data was mapped. Together, the sequencing depth and coverage attained in our study includes most of the known miRNAs and has the potential to discover new and low expressing ones.

Unique changes in miRNA profile in eat-2(ad1116) worms

We compared the normalized miRNA expression (TPM, see experimental procedures) of WT with eat-2(ad1116) either in day 1/young-adults or in day 8/aging worms and found that most of the miRNAs were distinctly upregulated in young-adult eat-2(ad1116) (Figure 1A, Figure S3 for analysis flow chart, Table 1, Table S3). For this, we considered all miRNAs that had read counts over 10 in both the strains on a particular day, with minimum fold change of ±1.5 having a p<0.05. Out of a total of 184 miRNA common between young-adult WT and eat-2(ad1116) on day 1, 105 were significantly upregulated in eat-2(ad1116). Interestingly, none were significantly downregulated. Seventy nine miRNA had no significant change in expression pattern. We found that 40 miRNA were exclusively expressed in eat-2(ad1116) (Table S4); out of them cel-miR-41-5p, cel-miR-40-5p, cel-miR-39-5p, cel-miR-4813-5p had read counts more than 100. For further bioinformatic analysis using these upregulated miRNAs in day 1/young-adult worms, we did not consider the strandedness of the miRNAs and only analysed 82 out of the 105 as uniquely upregulated miRNAs.

Differential regulation of miRNA population in WT and eat-2(ad1116). A scatter plot of miRNAs commonly expressed between WT and eat-2(ad1116) in (A) day 1 young-adult worms or (B) aging day 8 worms. For this, we considered 120 (day 1) or 78 (day 8) miRNAs (see Table S3) that changed expression in eat-2(ad1116) as compared to WT and had read counts >10. (C) Validation of sequencing results by quantitative real-time PCR. The miRNA expresion of young-adult eat-2(ad1116) was compare to WT worms collected on day 1 of adulthood. Error bars indicate standard deviation between independent biological replicates. The p-values were calculated byt-test (* represent p

Figure 1. Differential regulation of miRNA population in WT and eat-2(ad1116). A scatter plot of miRNAs commonly expressed between WT and eat-2(ad1116) in (A) day 1 young-adult worms or (B) aging day 8 worms. For this, we considered 120 (day 1) or 78 (day 8) miRNAs (see Table S3) that changed expression in eat-2(ad1116) as compared to WT and had read counts >10. (C) Validation of sequencing results by quantitative real-time PCR. The miRNA expresion of young-adult eat-2(ad1116) was compare to WT worms collected on day 1 of adulthood. Error bars indicate standard deviation between independent biological replicates. The p-values were calculated byt-test (* represent p < 0.05).

Table 1. List of miRNAs that were significantly upregulated in young-adult/Day 1 eat-2(ad1116) worms as compare to WT (more details in Table S3)

S. No.miRNAFold Changep-valueS. No.miRNAFold Changep-value
1cel-miR-82-5p47.46<0.000154cel-miR-253-3p3.3<0.0001
2cel-miR-38-3p35.19<0.000155cel-miR-237-3p3.290.0300
3cel-miR-74-5p23.950.000856cel-miR-783.220.0039
4cel-miR-54-5p23.41<0.000157cel-miR-784-5p3.16<0.0001
5cel-miR-44-5p20.61<0.000158cel-miR-240-3p3.15<0.0001
6cel-miR-35-5p20.59<0.000159cel-miR-73-3p3.12<0.0001
7cel-miR-81-5p20.33<0.000160cel-miR-58-5p3.07<0.0001
8cel-miR-229-3p19.88<0.000161cel-miR-77-3p3.07<0.0001
9cel-miR-35-3p18.36<0.000162cel-miR-85-3p3.05<0.0001
10cel-miR-73-5p14.88<0.000163cel-miR-43-3p3.03<0.0001
11cel-miR-36-3p14.83<0.000164cel-miR-787-3p2.930.0400
12cel-miR-41-3p14.65<0.000165cel-miR-60-3p2.88<0.0001
13cel-miR-39-3p13.03<0.000166cel-miR-235-3p2.87<0.0001
14cel-miR-40-3p12.78<0.000167cel-miR-786-3p2.780.0400
15cel-miR-255-3p12.60.000568cel-miR-80-5p2.73<0.0001
16cel-miR-79-5p12.420.003869cel-miR-7852.72<0.0001
17cel-miR-37-3p11.92<0.000170cel-miR-63-3p2.71<0.0001
18cel-miR-66-3p10.850.040071cel-miR-49-3p2.7<0.0001
19cel-miR-37-5p10.670.040072cel-miR-18172.680.0060
20cel-miR-72-5p10.45<0.000173cel-miR-231-3p2.64<0.0001
21cel-miR-64-3p10.180.010074cel-miR-66-5p2.54<0.0001
22cel-miR-71-3p9.43<0.000175cel-miR-74-3p2.45<0.0001
23cel-miR-77-5p8.72<0.000176cel-miR-51-5p2.3<0.0001
24cel-miR-252-5p8.08<0.000177cel-miR-47-3p2.29<0.0001
25cel-miR-238-5p7.07<0.000178cel-miR-3592.240.0100
26cel-miR-83-5p6.530.030079cel-miR-64-5p2.23<0.0001
27cel-miR-83-3p6.43<0.000180cel-miR-34-5p2.220.0100
28cel-miR-56-5p6.17<0.000181cel-miR-7932.19<0.0001
29cel-miR-72-3p6.070.040082cel-miR-54-3p2.190.0064
30cel-miR-52-3p5.03<0.000183cel-miR-259-5p2.160.0037
31cel-miR-42-3p5.03<0.000184cel-miR-34-3p2.080.0400
32cel-miR-84-3p4.880.000885cel-miR-3552.080.0006
33cel-miR-228-3p4.87<0.000186cel-miR-239a-5p2.01<0.0001
34cel-miR-241-5p4.74<0.000187cel-miR-10181.950.0019
35cel-miR-71-5p4.37<0.000188cel-miR-86-3p1.890.0036
36cel-miR-238-3p4.32<0.000189cel-miR-1-3p1.85<0.0001
37cel-miR-49374.070.040090cel-miR-2481.840.0006
38cel-miR-70-5p4.07<0.000191cel-miR-46-3p1.82<0.0001
39cel-miR-1829b4.040.005292cel-miR-246-3p1.81<0.0001
40cel-miR-75-3p3.96<0.000193cel-miR-65-5p1.8<0.0001
41cel-miR-45-5p3.950.020094cel-miR-229-5p1.73<0.0001
42cel-miR-75-5p3.940.000195cel-miR-56-3p1.72<0.0001
43cel-miR-47-5p3.80.000596cel-miR-87-3p1.710.0002
44cel-miR-358-3p3.790.040097cel-miR-2541.70<0.0001
45cel-miR-1-5p3.690.000198cel-miR-243-3p1.700.0010
46cel-miR-67-3p3.66<0.000199cel-miR-1830-3p1.69<0.0001
47cel-miR-240-5p3.62<0.0001100cel-miR-61-3p1.68<0.0001
48cel-miR-239b-5p3.55<0.0001101cel-miR-79-3p1.68<0.0001
49cel-miR-230-5p3.550.0002102cel-miR-236-3p1.61<0.0001
50cel-miR-51-3p3.5<0.0001103cel-miR-86-5p1.51<0.0001
51cel-miR-5592-5p3.4<0.0001104cel-miR-244-5p1.50<0.0001
52cel-miR-5592-3p3.40.0041105cel-miR-2-3p1.500.0500
53cel-miR-55-3p3.32<0.0001
MiRNAs that were used for QRT-PCR validation are in bold.
MiRNAs that did not possess any PHA-4 binding peaks are indicated in red. All the other miRNAs had one or more binding peaks within 5 kb upstream and 1 kb downstream of the transcription start site.

In order to validate the sequencing results, we chose 9 miRNAs having varying expression levels. Using quantitative real time PCR, we could validate the expression pattern of eight of these miRNAs; they were significantly upregulated in eat-2(ad1116) compared to WT on day 1 (Figure 1C).

Next we compared the miRNA profiles of aging WT and eat-2(ad1116) worms collected at day 8 of adulthood (Figure 1B, Table S3). We found 157 miRNAs common between WT and eat-2(ad1116); 15 were significantly upregulated in eat-2(ad1116) while 43 were downregulated compared to WT. Twenty miRNAs had no significant change in expression pattern. Four miRNAs were exclusively expressed in eat-2(ad1116) while 30 were present only in WT; out of them 4 had read counts more than 100 (Table S5).

From the above observations, it appeared that many miRNAs that are upregulated in young-adult eat-2(ad1116) on day 1 returned to the levels of WT or lower with progressing age. When we compared the expression patterns of miRNAs that were expressed in all the four samples with a minimum read count of 10, we found that 61 of the 152 (~40%) were upregulated in eat-2(ad1116) on day 1 but returned to near WT levels on day 8 (Cluster A in Figure 2A; Table S6). These included miRNAs cel-miR-70, cel-miR-51, cel-miR-56, cel-miR-230, cel-miR-238, cel-miR-83, etc. There were 10 others that were upregulated in day 1 of eat-2(ad1116), but remained upregulated even in day 8 (Cluster B). These included miRNAs cel-miR-237, cel-miR-238, cel-miR-229, cel-miR-791 and cel-miR-72. Levels of miRNAs cel-miR-241, cel-miR-53, cel-miR-57 and cel-miR-228 were decreased in day 8 in WT compare to day 1 levels, but failed to so in eat-2(ad1116) (Cluster E). Several miRNAs had expression patterns that were reverse of that in WT worms. For e.g., 11 miRNA, including cel-miR-63, cel-miR-53, cel-miR-36-41, cel-miR-254 and cel-miR-61 increased in expression from day 1 to 8 in WT worms; however, in eat-2(ad1116) it had an inverse expression pattern (Cluster C). As a corollary, only cel-miR-58 had the opposite expression pattern-decreased in WT but increased in eat-2(ad1116). Apart from miRNAs expressing in all samples, another 79 miRNA had more than 10 read count in any one or more of the 4 samples and were categorized separately (Table S7). Most of these miRNAs had low abundance as indicated by their read counts that were less than 1000. Only cel-miR-86 had 1258 read counts (TPM 58.72) on day 1 but had no detectable count on day 8. In this list also, we found 29 miRNAs were exclusively expressed in eat-2(ad1116) on day 1. These observations point to the profound changes in miRNA expression during eat-2(ad1116)-mediated dietary restriction, possibly to support metabolic changes associated with the long life span. MicroRNAs were, in general, found to be upregulated in eat-2(ad1116) young adult worms on day 1 of adulthood.

Age-dependent changes in miRNA expression. (A) A Heatmap showing expression patterns of 152 miRNAs that are commonly expressed in WT and eat-2(ad1116) at day 1 (young-adults) and day 8 (aging worms) of adulthood. The miRNAs are grouped into distinct clusters A-G based on the similarity in their expression patterns. (B) A Heatmap comparing the expression of miRNAs in WT, eat-2(ad1116) and daf-2(e1370).

Figure 2. Age-dependent changes in miRNA expression. (A) A Heatmap showing expression patterns of 152 miRNAs that are commonly expressed in WT and eat-2(ad1116) at day 1 (young-adults) and day 8 (aging worms) of adulthood. The miRNAs are grouped into distinct clusters A-G based on the similarity in their expression patterns. (B) A Heatmap comparing the expression of miRNAs in WT, eat-2(ad1116) and daf-2(e1370).

The miRNAome of eat-2(ad1116) worms is different from the long-lived insulin-like signalling mutant

Next we asked whether the profile obtained with eat-2(ad1116) is a general pattern for all long-lived mutants. In worms, the insulin-like signalling pathway is a major regulator of longevity and development [32]. Mutations in the insulin-like receptor daf-2 leads to dramatic increase in life span [33]. However, the mechanisms by which DR and insulin-like signalling increase life span are independent of each other [6]. We profiled the miRNAome of daf-2(e1370) at day 1 of adulthood and found that only 24 miRNAs were upregulated significantly while 33 were downregulated (Table S1, S18). All 24 upregulated miRNAs were also commonly detected in eat-2(ad1116) young adult worms, although the levels were higher in the latter in most cases (Figure 2B). In stark contrast to daf-2(e1370), none of the miRNAs were downregulated in eat-2(ad1116) as mentioned above. Together these data showed that the profile of miRNA in eat-2(ad1116)-mediated DR worms is unique and does not represent a common trend of long-lived mutants. However, some of the commonly upregulated miRNAs may still correspond to a signature for long-lived mutants, but needs to be validated in other strains that have enhanced longevity.

Novel miRNA that change expression in eat-2(ad1116)

Since DR induces a dramatic reprogramming of development and metabolism, we hypothesized that novel miRNAs may contribute towards this effect. Additionally, with increased depth of sequencing, we expected to detect novel miRNAs in our WT samples also. We focused our attention on the 877496 (in young adults on day 1) and 682280 (in day 8/aging worms) non-redundant un-annotated reads that aligned to the C. elegans genome in eat-2(ad1116) (Table S1). In order to discover novel miRNAs, we used the miRNA discovery package miRDeep2 (see Supporting information for reference). We found 48 novel miRNAs in young-adults and 9 in aging eat-2(ad1116) worms, with a read count greater than 10. Additionally following analysis of unassigned reads, in WT we identified 21 in young-adults and 35 in aging worms (8 and 9). Out of the 13 commonly predicted miRNAs between day 1/young-adult WT and eat-2(ad1116) (without read count or p value cut-off) worms, 9 were upregulated and 3 down-regulated as calculated on the basis of read count (Table S10). On the other hand, 2 out of 3 common miRNAs on day 8 were downregulated in eat-2(ad1116). Thus, similar to the known miRNA, we found distinct upregulation of novel miRNA population in young-adult eat-2(ad1116) worms.

In order to select predicted novel miRNA for experimental validation, we used one additional criterion. Apart from read count filter of ≥10 and presence of a star sequence, we used MiRDeep Score >1 as a cut-off. Using these criteria, we identified 7 predicted miRNA in WT and 15 in eat-2(ad1116) young-adult worms on day 1, and 5 in aging WT as well as 1 in aging eat-2(ad1116). From these, we selected 3 predicted novel miRNA for validation (Figure 3A, Table S11). We grew WT, alg-1(gk214) or alg-2(ok304) and checked the expression of the predicted miRNAs at young-adult stage. The ALG-1 and ALG-2 proteins are required for miRNA maturation [34]. We found that all of the candidate miRNAs expressed in an alg-2-dependent manner while 2 expressed in alg-1-dependent manner as determined by quantitative real-time PCR analysis (Figure 3B), confirming that they are indeed bona fide miRNAs. Out of these three novel miRNA, two are down-regulated in eat-2(ad1116) young-adult worms while one is upregulated when compared to WT (Figure 3C), matching their expression as determined by sequencing. Together, it appears that the profound transcriptional and post-transcriptional changes required following DR in eat-2(ad1116) involves the role of additional novel miRNAs.

Expression of novel miRNAs during DR. (A) The secondary structures of 3 predicted novel miRNA candidates as generated using the miRDeep2 program. (B) Expression of predicted novel miRNAs were examined in specific Argonaut gene mutants alg-1(gk214) or alg-2(ok304) by QRT-PCR (C) Expression changes of novel miRNAs in eat-2(ad1116) as compare to WT collected on day 1 of adulthood (young-adult worms) by qRT-PCR. Error bars indicate Standard Deviation between independent biological replicates. The p-values were calculated by t-test (* represents p

Figure 3. Expression of novel miRNAs during DR. (A) The secondary structures of 3 predicted novel miRNA candidates as generated using the miRDeep2 program. (B) Expression of predicted novel miRNAs were examined in specific Argonaut gene mutants alg-1(gk214) or alg-2(ok304) by QRT-PCR (C) Expression changes of novel miRNAs in eat-2(ad1116) as compare to WT collected on day 1 of adulthood (young-adult worms) by qRT-PCR. Error bars indicate Standard Deviation between independent biological replicates. The p-values were calculated by t-test (* represents p < 0.05).

PHA-4 regulates majority of the upregulated miRNA in eat-2(ad1116)

The FOXA transcription factor PHA-4 is an absolute requirement for DR-induced longevity in eat-2 mutants [11]. We asked whether PHA-4 could be directly regulating the expression of the miRNAs upregulated in eat-2(ad1116). For this, we used ChIPBase (See Supporting information for reference), a database that catalogues transcription factor binding maps on miRNA (as well as other small RNAs) and protein-coding gene promoters using ChIP-seq data from published work [29, 30]. We queried the database with all the miRNAs that were upregulated in young-adult eat-2(ad1116) worms collected on day 1. We found that 65 of the 82 promoters of upregulated miRNA have one or more PHA-4 binding peaks (Table 1), potentially pointing at direct transcriptional regulation of these miRNA by PHA-4. In order to validate this data, we chose 9 of these miRNAs that are induced in eat-2(ad1116) as compared to WT and checked their expression in eat-2(ad1116) mutants, both in presence or absence of PHA-4. We found that all these miRNAs, except cel-mir-47, are regulated in a PHA-4-dependent manner (Figure 4A). Thus, PHA-4 not only controls mRNA levels of coding genes, like sod-2 and sod-4 during DR [11], it possibly regulates the transcription of a large part of the miRNA population.

The miRNAs and mRNAs upregulated during DR are potential direct targets of PHA-4. (A) QRT-PCR analysis shows the expression levels of the miRNAs in eat-2(ad1116);rrf-3(pk1246) compared to rrf-3(pk1246), when grown on control or pha-4 RNAi. (B) QRT-PCR analysis showing the expression levels of predicted PHA-4 direct target genes in eat-2(ad1116);rrf-3(pk1246) compared to rrf-3(pk1246), grown either on control and pha-4 RNAi. The error bars represent Standard Deviation between independent biological samples. ‘ns’ represents p ≥ 0.05 by t-test.

Figure 4. The miRNAs and mRNAs upregulated during DR are potential direct targets of PHA-4. (A) QRT-PCR analysis shows the expression levels of the miRNAs in eat-2(ad1116);rrf-3(pk1246) compared to rrf-3(pk1246), when grown on control or pha-4 RNAi. (B) QRT-PCR analysis showing the expression levels of predicted PHA-4 direct target genes in eat-2(ad1116);rrf-3(pk1246) compared to rrf-3(pk1246), grown either on control and pha-4 RNAi. The error bars represent Standard Deviation between independent biological samples. ‘ns’ represents p ≥ 0.05 by t-test.

DR-induced miRNAs in eat-2(ad1116) may target a large part of the transcriptome

We discovered dramatic changes in miRNA expression levels with DR in eat-2(ad1116). We also identified PHA-4 as a transcription factor responsible for regulating most (65 out of 82) of these miRNAs.

MiRNAs target the 3′UTR of mRNA and lead to translational arrest or degradation depending on the extent of base pairing in the ‘seed’ sequence. To throw light on the biological processes being regulated by the PHA-4-dependent miRNAs in eat-2(ad1116), we predicted their target genes. Since the prediction of miRNA target is not always accurate, we used two different programs, miRanda as well as Targetscan (see Supporting information for references) and considered only the ones that came up with both the searches (Table S12). The two programs together predicted a total of 5145 possible unique targets of the 65 upregulated miRNAs (Table S13). We found that while miR-71 had the largest number of predicted targets (782), miR-243 has 34 targets. We also found that multiple miRNAs target the same mRNA (Table S13). For example, syg-1 is targeted by 18 miRNA, ptc-1 by 16, whereas ima-3 and ina-1 are targeted by 13 and 15 of the upregulated miRNA, respectively. SYG-1 is required cell autonomously in the HSNL neuron to determine synaptic specificity with vulval muscles and VC neurons [35]. PTC-1 is a Drosophila patched homolog, required for germline cytokinesis [36]. INA-1 is required during development for migration of neurons, coelomocyte precursor and distal tip cells of somatic gonad [37] while IMA-3 is the worm homolog of importin alpha involved oogenesis [38]. These genes may therefore be potentially important targets for miRNAs during DR.

Mature miRNAs, in most cases, bind to the 3′UTR of target mRNAs leading to their translational arrest [15]. Under this premise, we would expect to see down-regulation of protein levels corresponding to these predicted targets. In a previous study [39], the complete proteome of WT and eat-2(ad1116) were compared using quantitative proteomics. Out of the 167 proteins that were significantly downregulated in eat-2(ad1116), we found that 74 of them are predicted targets of the miRNAs (hypergeometric analysis, p=2.67e-8; Table S14). This suggests that PHA-4-regulated miRNAs may translationally inhibit many of their target mRNAs.

Many genes are upregulated in eat-2(ad1116) in a PHA-4-dependent manner

Next, in an effort to study the transcriptional response to DR, we compared the transcriptome of young-adult WT and eat-2(ad1116) worms collected on day 1 of adulthood, using NGS. We generated ~ 41 million reads for WT and ~42.2 million reads for eat-2(ad1116) using the GAIIx platform, as described in the experimental procedure. About 96.4% of WT reads and 97.2% of eat-2(ad1116) reads mapped to the C. elegans genome. We found that similar to the miRNA results, many more mRNAs were upregulated (3607 mRNAs) in eat-2(ad1116) compared to those that were downregulated (231 mRNAs)(Table S15).

Since PHA-4 plays such a central role in DR-mediated longevity, we asked whether these upregulated genes are under the direct control of PHA-4. For this, we looked for PHA-4 binding in the 5 kb promoter proximal region of the genes upregulated in eat-2(ad1116) using ChIPbase. We found that out of the 3607 upregulated transcripts, 2037 have at least one PHA-4 binding site (p=6.18e-171, by Hypergeometric test)(Table S15). To validate whether the presence of binding sites correlate with PHA-4 dependence in transcription, we randomly chose 10 genes. Using qRT-PCR, we found that 8 transcripts are dependent on PHA-4 (Figure 4B). These data suggests that PHA-4 may bind and upregulate a large number of genes in eat-2(ad1116); this may be one of the reasons why knocking down PHA-4 completely suppresses DR-induced life span.

Feed forward loops involving PHA-4 along with its target mRNA and miRNA

We have found that the FOXA transcription factor PHA-4 may be a critical requirement in upregulation of mRNAs as well as miRNA in eat-2(ad1116). Next, we asked whether there was an overlap between genes that are transcriptionally upregulated by PHA-4 and the targets of PHA-4-regulated miRNA. By comparing 2037 PHA-4-regulated genes to 5145 genes that are predicted targets of PHA-4-regulated miRNA, we found a significant overlap (1073 genes with a p = 3.16e-176 by Hypergeometric test; Table S16). Thus, PHA-4 forms a large number of incoherent feed forward loops [40] during DR in eat-2(ad1116) (Figure 5A, B). While on one hand, PHA-4 transcriptionally upregulates the expression of a many transcripts, the same set of transcripts are potentially targeted by miRNAs that it transcribes (Figure 5B). In fact, some of the target genes are controlled by upto 20 miRNAs indicating that they might potentially play important roles during DR.

PHA-4 controls a large number of genes by feed-forward loops using miRNAs. (A) A representation of transcription factor-miRNA feed forward loop involving PHA-4 that controls life span during DR. (B) A complex regulatory network involving PHA-4 and its target mRNAs that are also regulated by PHA-4-controlled miRNAs. The network was generated using Cytoscape v3.1.1. Only a subset of genes targeted by PHA-4-regulated miRNA [as determined by mirTarBase were used to generate a representative decongested network. The boxes corresponding to the miRNA or mRNA/genes are colored based on their fold induction in eat-2(ad1116) as compared to WT. Edges are colored based on the number of miRNAs regulating a particular gene. Dark red are for genes regulated by 20 or more miRNAs while those regulated by less than 10 miRNAs are in black. MiRNAs that have no target listed in mirTarBase are in red. (C) The genes targeted by PHA-4 transcriptionally as well as post-transcriptionally are highly connected by protein-protein interactions, as determined by STRING software.

Figure 5. PHA-4 controls a large number of genes by feed-forward loops using miRNAs. (A) A representation of transcription factor-miRNA feed forward loop involving PHA-4 that controls life span during DR. (B) A complex regulatory network involving PHA-4 and its target mRNAs that are also regulated by PHA-4-controlled miRNAs. The network was generated using Cytoscape v3.1.1. Only a subset of genes targeted by PHA-4-regulated miRNA [as determined by mirTarBase were used to generate a representative decongested network. The boxes corresponding to the miRNA or mRNA/genes are colored based on their fold induction in eat-2(ad1116) as compared to WT. Edges are colored based on the number of miRNAs regulating a particular gene. Dark red are for genes regulated by 20 or more miRNAs while those regulated by less than 10 miRNAs are in black. MiRNAs that have no target listed in mirTarBase are in red. (C) The genes targeted by PHA-4 transcriptionally as well as post-transcriptionally are highly connected by protein-protein interactions, as determined by STRING software.

Next we asked whether the insulin-like signalling pathway also uses such FFLs to control gene expression. The insulin-like signalling pathway components phosphorylate and inactivate the FOXO transcription factor DAF-16 that is a major output of this cascade [32]. Under low insulin signalling conditions, such as in daf-2(e1370), DAF-16 translocates into the nucleus and transactivates a large assortment of genes [32]. Using ChIPbase, we first determined that 22 of the 24 miRNA upregulated in daf-2(e1370) are potential direct targets of DAF-16. Next, using a published daf-2(e1370) transcriptomics data [41], we found that 687 transcripts are upregulated in daf-2(e1370); out of these 336 have DAF-16 binding sites within 5 kb of promoter, according to ChIPbase. Thereafter, we overlapped the predicted targets of the 22 DAF-16-regulated miRNA (2235 transcripts) with that of the DAF-16-dependent mRNA. We found that only 47 genes overlapped (13.98%; p=0.01 by hypergeometric test) in contrast to that of eat-2(ad1116) where PHA-4 regulated 52.67% (p = 3.16e-176) of the upregulated transcripts by FFLs. This data suggested that PHA-4-mediated FFLs are a unique feature of DR-induced longevity in eat-2(ad1116) and may not be used extensively by DAF-16/FOXO downstream of the insulin-like signalling pathway.

PHA-4 target genes regulated through FFLs control important biological functions

Next we determined the functions of all the genes that are potentially transcriptionally as well as post-transcriptionally regulated by PHA-4. Since these genes are under such tight regulation, we can expect them to be important for DR-mediated longevity. We categorized these 1073 genes based on their Gene Ontology (GO) terms as well as determined which molecular pathways these genes function in. We used “Gene Ontology enRIchment anaLysis and visualization” tool (GOrilla) and visualized the results using Revigo. We found that the GO terms associated with metabolism, protein folding and response to various chemical stimuli are enriched (Figure 6A). These genes are involved in ubiquitin mediated protein degradation, lysosomal autophagy, inositol phosphate metabolism, phosphatidylinositol signalling, ErbB signalling as determined by KEGG pathway analysis using DAVID (see Supporting information for references) (Figure 6B). These 1073 genes are also highly connected with each other through protein-protein interaction as determined by STRING database (Figure 5C). Further, these genes are enriched for transcription factors (~10% of the 1073 genes code for transcription factors as against 4.4% found in the genome; p=1.68e-12 by Hypergeometric test). Interestingly, when we compared this list of genes with the GeneAge database (http://genomics.senescence.info/genes/), we found that 86 genes are already known to be important in longevity determination having p=1.66e-12 (Table 2). Together, the genes that are regulated by PHA-4 directly or indirectly through miRNAs most likely are important players in DR-mediated longevity in eat-2(ad1116).

Genes regulated by PHA-4,transcriptionally as well as post-transcriptionally using miRNA are involved in important biological functions during DR. (A) The 1073 genes are enriched in GO terms pertaining to protein folding, cellular metabolism and response to external stimuli. GO analysis was performed using GORILA software and visualized with REVIGO. (B) KEGG pathway components for signaling, lysosome and ubiquitin-mediated decay are enriched in genes that are regulated by PHA-4-controlled FFLs.

Figure 6. Genes regulated by PHA-4,transcriptionally as well as post-transcriptionally using miRNA are involved in important biological functions during DR. (A) The 1073 genes are enriched in GO terms pertaining to protein folding, cellular metabolism and response to external stimuli. GO analysis was performed using GORILA software and visualized with REVIGO. (B) KEGG pathway components for signaling, lysosome and ubiquitin-mediated decay are enriched in genes that are regulated by PHA-4-controlled FFLs.

Table 2. List of aging-associated genes that are transcriptionally as well as post-transcriptionally regulated by PHA-4

Gene NameEffect on lifespanFunctions
acdh-13IncreaseAcyl CoA DeHydrogenase
akt-1IncreaseAKT kinase family
aps-1IncreaseAdaPtin, Small chain (clathrin associated complex)
bec-1DecreaseBEClin (human autophagy) homolog
C47D12.2Increase-
C56G2.1IncreaseKH domain-containing protein
ccr-4DecreaseCCR (yeast CCR4/NOT complex component) homolog
cec-3IncreaseChromo domain-containing protein
ceh-20DecreaseHomeobox containing protein
chc-1DecreaseClathrin Heavy Chain
cua-1DecreaseCU (copper) ATPase
daf-12IncreaseA member of the steroid hormone receptor superfamily
daf-14IncreasePossibly a transducer of the DAF-7/TGF-beta-mediated signal that promotes reproductive growth
daf-16DecreaseForkhead family transcription factor
daf-3DecreaseTranscriptional regulator
daf-4IncreaseTransmembrane serine/threonine kinase
din-1DecreaseDAF-12 Interacting Protein
dlk-1DecreaseDAP (Death Associated Protein kinase) Like Kinase
dpy-27DecreaseATP-binding protein that is a homolog of the SMC4 subunit of mitotic condensin
dve-1DecreaseDVE (Defective proVEntriculus in Drosophila) homolog
egl-18IncreaseGATA family transcription factor
egl-27IncreaseHomolog of human MTA1 (metastasis-associated protein), a part of a nucleosome remodeling and histone deacetylation (NURD) complex
egl-8Increasephospholipase C beta
egl-9IncreaseProline hydroxylase
epc-1IncreaseEnhancer of PolyComb-like
erm-1IncreaseEzrin/Radixin/Moesin
ets-4IncreaseETS class transcription factor
F09F7.5Increase-
F14B4.3IncreaseDNA-directed RNA polymerase
gcy-28IncreaseGuanylyl CYclase
gei-4DecreaseGEX Interacting protein
hpk-1DecreaseHomeodomain interacting Protein Kinase
ire-1DecreaseIRE1 kinase related
itr-1IncreaseInositol Triphosphate Receptor
klo-2Decrease-
kri-1Decreasehuman KRIT 1 (Krev interaction trapped/cerebral cavernous malformation 1) homolog
ldb-1DecreaseLIM domain Binding Protein
let-23Increase and DecreaseEGF-receptor-family transmembrane tyrosine kinase
let-502DecreaseRho-binding Ser/Thr kinase, orthologous to human myotonic dystrophy kinase (DM-kinase)
let-60IncreaseMember of the GTP-binding RAS protooncogene family
lin-15BDecreaseNovel protein that contains a THAP domain, a zinc-coordinating, site-specific DNA-binding domain
lin-23DecreaseF-box- and WD-repeat-containing protein, components of SCF (Skp1, Cullin, F-box) ubiquitin-ligase complexes that function in ubiquitin-mediated protein degradation
lin-3DecreaseEGF family of peptide growth factors
lin-5Increase-
max-2IncreaseMotor AXon guidance
mep-1DecreaseMog interacting, Ectopic P granules
mfn-1IncreaseMitoferrin
mnk-1DecreaseMAP kinase iNtegrating Kinase (MNK) homolog
mpz-2IncreaseTemporarily Assigned Gene name
mrp-5IncreaseMultidrug Resistance Protein family
mtm-3DecreaseMTM (myotubularin) family
nac-2IncreaseNADC (Na+-coupled dicarboxylate transporter) family
nfi-1DecreaseTranscription/replication factor
nhr-14IncreaseNuclear Hormone Receptor family
nhr-23IncreaseNuclear Hormone Receptor family
nhr-49DecreaseNuclear Hormone Receptor family
nhx-2IncreaseNa/H eXchanger
nkcc-1IncreaseNa-K-Cl Cotransporter homolog
npa-1IncreaseFatty acids and vitamin A binding
pab-2IncreasePolyA Binding protein
pup-2IncreasePoly(U) Polymerase
R08E3.3Increase-
rab-5DecreaseRAB family
rbr-2DecreaseRB (Retinoblastoma Binding protein) Related
rpa-1IncreaseReplication Protein A homolog
sel-5IncreaseSuppressor/Enhancer of Lin-12
set-2IncreaseSET (trithorax/polycomb) domain containing
set-6DecreaseSET (trithorax/polycomb) domain containing
sgk-1IncreaseSerum- and Glucocorticoid-inducible kinase homolog
smg-1IncreaseSuppressor with Morphological effect on Genitalia
smk-1DecreaseSMEK (Dictyostelium Suppressor of MEK null) homolog
syd-2DecreaseRequired cell autonomously in neurons for differentiation of presynaptic active zones
tax-6IncreaseRequired for inhibition and adaptation of several sensory neurons and for the normal regulation of egg-laying by serotonin
tes-1IncreaseTemporarily Assigned Gene name
top-1DecreaseTOPoisomerase
unc-26IncreaseEncodes synaptojanin, a polyphosphoinositide phosphatase
unc-32IncreaseVacuolar H+-ATPase V0 sector, subunit a
unc-52IncreaseBasement-membrane (HSPG) core
unc-62IncreaseMeis-class homeodomain protein
unc-76IncreaseEncodes FEZ family protein involved in axon-axon interactions
wdr-23IncreaseWD Repeat protein
wip-1DecreaseWiskott-Aldrich syndrome protein (WASP)-Interacting Protein and gene assignment
wts-1IncreaseWarTS (Drosophila) protein kinase homolog
wwp-1DecreaseWW domain Protein (E3 ubiquitin ligase)
Y48G1A.4IncreaseRibosome biogenesis
zfp-1Decreaseleucine zipper, zinc finger, and PHD/LAP domain protein

Discussion

Although the role of DR in life span regulation in C. elegans is well-established, the underlying molecular mechanisms are only now beginning to be elucidated. The role of the FOXA transcription factor PHA-4 in DR was identified several years ago, but very little is known about its functional regulation during DR [11]. Initial studies had indicated that PHA-4 levels increase during DR in eat-2(ad1116) and the transcription factor controls genes involved in oxidative stress resistance (like sod-2 and sod-4) [11] or autophagy [42]. In this study, we show for the first time the complexity of gene regulation brought about by this transcription factor. We found that PHA-4 upregulates a large contingent of genes in eat-2(ad1116); however, it also targets these genes using miRNAs that it controls, forming incoherent FFLs. There are not many examples of such large sets of TF-miRNA FFLs in biology, although most of the miRNA have been predicted to function in this manner [43]. Such intricate mechanisms of gene expression control may make PHA-4 a major factor in DR-induced life span regulation. Perturbation of PHA-4 may destabilize the system during DR, leading to complete suppression of life span, as observed earlier [11].

PHA-4 is a central part of a large number of incoherent feed forward loops (Figure 5B), as evident from our study. A miRNA would bind to the 3′UTR of its target mRNA to inhibit translation or in some cases, lead to its destabilization and degradation [14-16]. One of the main objectives of such miRNA-mRNA interaction is to shield protein expression from fluctuations in transcription [17-20]. This mechanism also has a built-in rheostat that ensures that translation of a protein will effectively proceed only when the expression of mRNA picks up, titrating away the miRNA inhibition [21]. Increased noise in gene expression would lead to phenoypic variability during development or environmental perturbations, a process that is curbed by canalization or robustness in the system [21]. Thus, PHA-4 may be the master regulator of gene expression that is responsible to instil robustness in the system, reducing phenotypic variability during DR. By controlling a large number of miRNAs and mRNAs in from of incoherent feed-forward loops, the transcription factor also ensures an appropriate DR response only when nutrition levels fall to a critical level. This may effectively save a lot of cellular energy by preventing the system to mount an all-out defence against nutritional deprivation when there are only minor fluctuations in food availability.

The genes that are intensely protected from expression fluctuations by PHA-4 may constitute important components of the DR response in eat-2(ad1116). When we analysed the 1073 genes that are regulated by PHA-4 transcriptionally as well as post-transcriptionally, we found them to be categorized in KEGG pathways for endocytosis, lysosome-mediated autophagy, cellular signalling and ubiquitin-mediated proteolysis. These genes were also enriched for GO terms pertaining to intracellular signalling and cellular response to organic compounds including lipids, hormones etc., as well as protein folding.

It appears that genes involved in protein turnover, under the transcriptional control of PHA-4, are particularly targeted by the FFLs. Ubiquitin-mediated proteosomal degradation involving ubiquitin ligases is a well-known requirement for DR-mediated longevity as well as for starvation response [44-47]. The HECT (Homologous to E6AP Carboxy Terminus) E3 ubiquitin ligase WWP-1 is required for DR-induced life span extension while components of a SCF E3 ubiquitin ligase have been shown to be important for fasting-induced longevity [44, 46]. In fact, our data suggests that wwp-1 is one of the targets of PHA-4 that is also regulated by miRNA. Similarly, Skp1-Cullin-F-box (SCF) E3 ligase complex genes like skr-1, cul-3, cul-4, cul-5, lin-23 as well as Ubiquitin conjugating enzyme ubc-25 are also regulated through the PHA-4-dependent FFLs. The WD40 repeat protein WDR-23 along with the ubiquitin ligase CUL4/DDB1 functions to regulate the SKN-1/NRF2 transcription factor [48]. SKN-1 activity in the neurons is required for DR-mediated longevity [49]. We found that the wdr-23 and cul-4 genes are regulated by PHA-4 FFLs. Additionally, components of the lysosome-mediated autophagy are also targeted by FFLs. Beclin1 homolog in worms bec-1 is the major regulator of autophagy, a process that requires PHA-4 activity [50]. We found that bec-1 as well as other autophagic genes like atg-2 (yeast atg2p and human ATG2A/ATG2B homolog) and atg-11 (fly Atg17-PC homolog) may be regulated transcriptionally as well as post-transcrip-tionally by PHA-4, Thus, protein turnover that provides building blocks for synthesis of new proteins and fuel cell with energy during energy-deprived conditions are primary contributors of longevity [47] and thus may be under tight regulatory control by PHA-4.

Signal transduction pathways are other important components of DR that seems to be insulated from fluctuations by PHA-4 FFLs. Components of MAP kinase signalling, calcium signalling and Wnt signalling are enriched. Signalling components may play critical role during DR to maintain cellular homeostasis and tight control of their protein expression is required to minimize phenotypic fluctuations. Well known signal transduction components, like akt-1and sgk-1(mammalian AKT and SGK homologs), tir-1(Toll-like receptor homolog), sek-1 (MKK4 homolog that acts in the p38 pathway), pkc-1 (mammalian protein kinase C epsilon homolog), nipi-3 (the human kinase Tribbles homolog), mbk-1 (orthologous to Drosophila MINIBRAIN and mammalian DYRK1A/MnbK), jun-1 (mammalian JUN homolog), let-60 (member of the GTP-binding RAS protooncogene family) etc. are potentially regulated by the FFLs; these genes include ones with known role in aging and stress. Interestingly, transcription factor genes are also highly enriched in the genes regulated by FFLs. The system may be protecting these important transcription factors from gene expression fluctuations in order to maintain proper cellular functions during DR, leading to increased life span. In fact, the genes that are part of the FFLs downstream of PHA-4 are significantly enriched for genes that are involved in life span regulation.

While this study was in review, Smith-Vikos et al. [51] reported PHA-4 as the most well-connected factor in a transcription factor-miRNA regulatory network involved in aging. They identified two miRNAs that are potentially directly regulated by PHA-4 and are upregulated during two different DR interventions, bacterial dilution and in eat-2(ad1116). We also found that these two miRNA were significantly upregulated in our study, although to a greater extent. This difference may be attributed to the stages at which the worms were harvested for analysis. Together, both the studies attest to the central role of PHA-4 in controlling miRNA expression during DR to regulate life span.

In conclusion, our study reveals the complexity of gene regulation during DR in eat-2(ad1116), showing PHA-4 as a central regulator of gene expression that may infuse robustness in the system using a large number of FFLs to protect against phenotypic variations during periods of low food availability. It will be interesting to study whether alternative models of DR in worms or other organisms have similar post-transcriptional regulation by FOXA transcription factors.

Methods

Complete and detailed experimental procedure is available as Supporting information.

Strain maintenance

All strains were maintained at 20°C using standard C. elegans techniques. All RNAi experiments were initiated using synchronized L1 stage worms. Strains used in the study are: N2 Bristol as WT, eat-2(ad1116), rrf-3(pk1246), eat-2(ad1116);rrf-3(pk1246), alg-1(gk214) and alg-2(ok304).

For collecting worms for NGS analysis, synchronized L1 larvae were placed on Nematode Growth Media (NGM) agar plates and allowed to grow till L4 stage. The plates were overlaid with FuDR to stop the hatching of the eggs. Worms were collected from the plates by washing with 1 X M9 the following day (day 1/young-adult worms) or on day 8 (aging worms). The worm pellet was washed three times with 1 X M9 before Trizol reagent (Invitrogen, USA) was added and RNA isolated as described below.

RNA isolation

RNA isolation was performed using Trizol. Briefly, worms grown on OP50 bacteria, vector or RNAi of interest were washed off the plates with M9 buffer. Thereafter, 0.4 ml of Trizol reagent was added and the worms lysed by vigorous vortexing. RNA was purified by phenol:chloroform:isoamylalcohol extraction and ethanol precipitation. RNA integrity were confirmed by analysis on an Agilent 2100 bioanalyzer, using the RNA 6000 nano kit (Agilent) to confirm that the RIN numbers are above 8.0. Alternatively, the quality of the ribosomal RNA 28 S and 18 S as determined on an agarose gel was used as a measure of integrity and the absorbance at 260/280 nm was used to determine quantity.

RNA library preparation, Next Generation sequencing and data analysis. The small RNA libraries were constructed using Small RNA sample preparation kit v1.5 according to the manufacturer's instructions (Illumina Inc., USA). Briefly, the total RNA (2μg) were ligated to 3′ RNA adapter using RNA ligase truncated and 5′ adapter using T4 RNA ligase 2 (New England Biolabs, USA). The ligation products were reverse transcribed using Superscript II Reverse Transcriptase (Life Technologies, USA) and amplified with 12 cycles of PCR. The PCR products constituting the small RNA cDNA libraries were resolved on 6% Novex TBE PAGE Gel (Invitrogen, USA) and ~150 bp fragments excised. The library was eluted from the acrylamide gel and analyzed on Agilent 2100 Bioanalyzer using DNA high sensitivity kit (Agilent Technologies, USA). NGS of cDNA libraries were performed using Illumina GAIIX for 36 cycles. A total of 2.8 GB of raw sequence data, comprising of WT and eat-2(ad1116) strains at both day 1 and day 8, was imported into the CLC Genomics Workbench 6.5.1 (CLC Bio, Denmark). The reads were trimmed off the adapter sequences and reads containing low quality bases eliminated. The trimmed raw sequences were mapped to the miRBase release 19, allowing for a maximum of two gaps or mismatches. Unpaired group comparisons, based on Transcript Per Million (TPM), were used as expression values to compare different samples. A fold change ±1.5 with a minimum read count of ≥10 were used to filter the differentially expressed miRNA. The p value cutoff was set at p≤0.05 based on Kal's Z test statistical. The sequencing data is available at GEO repository with Series record number GSE60155 and GSE61112.

RNA-Sequencing (RNA-seq) libraries of WT and eat-2(ad1116) Day 1 samples were prepared as recommended by the Illumina TruSeq™ RNA Sample Preparation kit using Low-Throughput (LT) Protocol (Illumina, Inc.,USA). NGS of libraries was performed using Illumina GA IIX for 78 cycles including 6 additional cycles for index read. Sequence reads were aligned using CLC Genomics Workbench 6.5.1 with default setting against C. elegans genome assembly (WS231). Unpaired group comparisons, based on RPKM (Reads Per Kilobase per Million mapped reads), were chosen as expression values for comparing the samples. A fold change ±2.0 and P value ≤0.05 (Kal's Z test) were used to filter the differentially expressed genes.

Author Contributions

AM conceived the project and wrote the ms. AP and NK established the NGS core at NII and standardized protocols. AP performed all the experiments. VJ performed the bioinformatics analysis with assistance from AP. AP, VJ and AM analyzed the data.

Acknowledgments

We apologize to the authors whose original work is not included in the references due to space limitations. We thank members of the Molecular Aging Lab as well as NII NGS core facility for their help and support. We also thank Dr. S. Narasimhan for suggestions on writing the ms. This research was supported by intramural funding from National Institute of Immunology (NII) and from Department of Biotechnology (DBT), Govt. of India (BT/PR13720/BAB/10/779/2010). We thank DBT for a generous infrastructure grant for establishment of the NII NGS core facility. Some strains were provided by the CGC, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440). AM was supported by Ramalingaswami Fellowship. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflicts of Interest

The author of this manuscript declares no conflict of interest.

References

  • 1. Speakman JR and Mitchell SE. Caloric restriction. Mol Aspects Med. 2011; 32:159-221. [PubMed]
  • 2. McCay C, Crowell MF, Maynard L. The effect of retarded growth upon the length of life span and upon the ultimate body size. J nutr. 1935; 10:63.
  • 3. John OH and Luigi F. Caloric restriction in humans. Experimental Gerontology. 2007; 42:.
  • 4. Mair W and Dillin A. Aging and survival: the genetics of life span extension by dietary restriction. Annu Rev Biochem. 2008; 77:727-754. [PubMed]
  • 5. Mair W, Panowski SH, Shaw RJ, Dillin A. Optimizing dietary restriction for genetic epistasis analysis and gene discovery in C. elegans. PLoS One. 2009; 4:e4535 [PubMed]
  • 6. Lakowski B and Hekimi S. The genetics of caloric restriction in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 1998; 95:13091-13096. [PubMed]
  • 7. Greer EL, Dowlatshahi D, Banko MR, Villen J, Hoang K, Blanchard D, Gygi SP, Brunet A. An AMPK-FOXO pathway mediates longevity induced by a novel method of dietary restriction in C. elegans. Curr Biol. 2007; 17:1646-1656. [PubMed]
  • 8. Kaeberlein TL, Smith ED, Tsuchiya M, Welton KL, Thomas JH, Fields S, Kennedy BK, Kaeberlein M. Lifespan extension in Caenorhabditis elegans by complete removal of food. Aging Cell. 2006; 5:487-494. [PubMed]
  • 9. Chamoli M, Singh A, Malik Y, Mukhopadhyay A. A novel kinase regulates dietary restriction-mediated longevity in Caenorhabditis elegans. Aging Cell. 2014;.
  • 10. Vora M, Shah M, Ostafi S, Onken B, Xue J, Ni JZ, Gu S, Driscoll M. Deletion of microRNA-80 activates dietary restriction to extend C. elegans healthspan and lifespan. PLoS Genet. 2013; 9:e1003737.
  • 11. Panowski SH, Wolff S, Aguilaniu H, Durieux J, Dillin A. PHA-4/Foxa mediates diet-restriction-induced longevity of C. elegans. Nature. 2007; 447:550-555..
  • 12. Gurtan AM and Sharp PA. The role of miRNAs in regulating gene expression networks. Journal of molecular biology. 2013; 425:3582-3600. [PubMed]
  • 13. Kim VN, Jinju H, Mikiko CS. Biogenesis of small RNAs in animals. Nature reviews Molecular cell biology. 2009; 10:126-139..
  • 14. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004; 116:281-297. [PubMed]
  • 15. Baek D, Villén J, Shin C, Camargo FD, Gygi SP, Bartel DP. The impact of microRNAs on protein output. Nature. 2008; 455:64-71. [PubMed]
  • 16. Farh KK, Grimson A, Jan C, Lewis BP, Johnston WK, Lim LP, Burge CB, Bartel DP. The widespread impact of mammalian MicroRNAs on mRNA repression and evolution. Science. 2005; 310:1817-1821. [PubMed]
  • 17. Siciliano V, Garzilli I, Fracassi C, Criscuolo S, Ventre S, di Bernardo D. MiRNAs confer phenotypic robustness to gene networks by suppressing biological noise. Nature communications. 2013; 4:2364.
  • 18. Kim Dh, Grün D, van Oudenaarden A. Dampening of expression oscillations by synchronous regulation of a microRNA and its target. Nature genetics. 2013; 45:1337-1344. [PubMed]
  • 19. Matteo O, Carla B, Davide C, Michele C. The role of incoherent microRNA-mediated feedforward loops in noise buffering. PLoS computational biology. 2011; 7:.
  • 20. Strovas TJ, Rosenberg AB, Kuypers BE, Muscat RA, Seelig G. MicroRNA-based single-gene circuits buffer protein synthesis rates against perturbations. ACS synthetic biology. 2014; 3:324-331. [PubMed]
  • 21. Posadas DM and Carthew RW. MicroRNAs and their roles in developmental canalization. Curr Opin Genet Dev. 2014; 27C:1-6. [PubMed]
  • 22. Hornstein E and Shomron N. Canalization of development by microRNAs. Nat Genet. 2006; 38:Suppl S20-24. [PubMed]
  • 23. Devika G and Stephen MC. miRNAs and aging: A genetic perspective. Ageing research reviews. 2014;.
  • 24. Kato M and Slack FJ. Ageing and the small, non-coding RNA world. Ageing Res Rev. 2012;.
  • 25. Smith-Vikos T and Slack FJ. MicroRNAs and their roles in aging. J Cell Sci. 2012; 125:7-17. [PubMed]
  • 26. Sachi I and Frank S. MicroRNAs and the genetic network in aging. Journal of molecular biology. 2013; 425:3601-3608. [PubMed]
  • 27. Ibáñez-Ventoso C, Yang M, Guo S, Robins H, Padgett RW, Driscoll M. Modulated microRNA expression during adult lifespan in Caenorhabditis elegans. Aging cell. 2006; 5:235-246. [PubMed]
  • 28. Ibanez-Ventoso C and Driscoll M. MicroRNAs in C. elegans Aging: Molecular Insurance for Robustness? Curr Genomics. 2009; 10:144-153..
  • 29. Niu W, Lu ZJ, Zhong M, Sarov M, Murray JI, Brdlik CM, Janette J, Chen C, Alves P, Preston E, Slightham C, Jiang L, Hyman AA, et al. Diverse transcription factor binding features revealed by genome-wide ChIP-seq in C. elegans. Genome research. 2011; 21:245-254..
  • 30. Zhong M, Niu W, Lu ZJ, Sarov M, Murray JI, Janette J, Raha D, Sheaffer KL, Lam HY, Preston E, Slightham C, Hillier LW, Brock T, et al. Genome-wide identification of binding sites defines distinct functions for Caenorhabditis elegans PHA-4/FOXA in development and environmental response. PLoS Genet. 2010; 6:e1000848 [PubMed]
  • 31. Kato M, Chen X, Inukai S, Zhao H, Slack FJ. Age-associated changes in expression of small, noncoding RNAs, including microRNAs, in C. elegans. Rna. 2011; 17:1804-1820. [PubMed]
  • 32. Kenyon CJ. The genetics of ageing. Nature. 2010; 464:504-512. [PubMed]
  • 33. Kenyon C, Chang J, Gensch E, Rudner A, Tabtiang R. A C. elegans mutant that lives twice as long as wild type. Nature. 1993; 366:461-464. [PubMed]
  • 34. Grishok A, Pasquinelli AE, Conte D, Li N, Parrish S, Ha I, Baillie DL, Fire A, Ruvkun G, Mello CC. Genes and mechanisms related to RNA interference regulate expression of the small temporal RNAs that control C. elegans developmental timing. Cell. 2001; 106:23-34. [PubMed]
  • 35. Shen K and Bargmann CI. The immunoglobulin superfamily protein SYG-1 determines the location of specific synapses in C. elegans. Cell. 2003; 112:619-630. [PubMed]
  • 36. Kuwabara PE, Lee MH, Schedl T, Jefferis GS. A C. elegans patched gene, ptc-1, functions in germ-line cytokinesis. Genes & development. 2000; 14:1933-1944. [PubMed]
  • 37. Baum PD and Garriga G. Neuronal migrations and axon fasciculation are disrupted in ina-1 integrin mutants. Neuron. 1997; 19:51-62. [PubMed]
  • 38. Geles KG and Adam SA. Germline and developmental roles of the nuclear transport factor importin alpha3 in C. elegans. Development (Cambridge, England). 2001; 128:1817-1830..
  • 39. Yuan Y, Kadiyala CS, Ching TT, Hakimi P, Saha S, Xu H, Yuan C, Mullangi V, Wang L, Fivenson E, Hanson RW, Ewing R, Hsu AL, et al. Enhanced energy metabolism contributes to the extended life span of calorie-restricted Caenorhabditis elegans. J Biol Chem. 2012; 287:31414-31426. [PubMed]
  • 40. Alon U. Network motifs: theory and experimental approaches. Nature reviews Genetics. 2007; 8:450.
  • 41. Riedel C, Dowen R, Lourenco G, Kirienko N, Heimbucher T, West J, Bowman S, Kingston R, Dillin A, Asara J, Ruvkun G. DAF-16 employs the chromatin remodeller SWI/SNF to promote stress resistance and longevity. Nature cell biology. 2013; 15:491.
  • 42. Lapierre LR, Gelino S, Melendez A, Hansen M. Autophagy and lipid metabolism coordinately modulate life span in germline-less C. elegans. Curr Biol. 2011; 21:1507-1514..
  • 43. Tsang J, Zhu J, van Oudenaarden A. MicroRNA-mediated feedback and feedforward loops are recurrent network motifs in mammals. Molecular cell. 2007; 26:753-767. [PubMed]
  • 44. Uno M, Honjoh S, Matsuda M, Hoshikawa H, Kishimoto S, Yamamoto T, Ebisuya M, Yamamoto T, Matsumoto K, Nishida E. A fasting-responsive signaling pathway that extends life span in C. elegans. Cell reports. 2013; 3:79-91. [PubMed]
  • 45. Lucanic M and Kapahi P. Ubiquitin ligases join the field of dietary restriction in C.elegans. Aging. 2009; 1:751-752. [PubMed]
  • 46. Carrano AC, Liu Z, Dillin A, Hunter T. A conserved ubiquitination pathway determines longevity in response to diet restriction. Nature. 2009; 460:396-399. [PubMed]
  • 47. Vellai T and Takacs-Vellai K. Regulation of protein turnover by longevity pathways. Adv Exp Med Biol. 2010; 694:69-80. [PubMed]
  • 48. Choe KP, Przybysz AJ, Strange K. The WD40 repeat protein WDR-23 functions with the CUL4/DDB1 ubiquitin ligase to regulate nuclear abundance and activity of SKN-1 in Caenorhabditis elegans. Mol Cell Biol. 2009; 29:2704-2715. [PubMed]
  • 49. Bishop NA and Guarente L. Two neurons mediate diet-restriction-induced longevity in C. elegans. Nature. 2007; 447:545-549. [PubMed]
  • 50. Hansen M, Chandra A, Mitic LL, Onken B, Driscoll M, Kenyon C. A role for autophagy in the extension of lifespan by dietary restriction in C. elegans. PLoS genetics. 2008; 4:.
  • 51. Smith-Vikos T, de Lencastre A, Inukai S, Shlomchik M, Holtrup B, Slack FJ. MicroRNAs Mediate Dietary-Restriction-Induced Longevity through PHA-4/FOXA and SKN-1/Nrf Transcription Factors. Curr Biol. 2014; 24:2238-2246. [PubMed]
Download PDF
Corresponding Author

Arnab Mukhopadhyay, PhD
arnab@nii.ac.in

Keywords
microRNA dietary restriction PHA-4/FOXA Transcription factor miRNA feed forward loops aging
Paper Sections