Roles of the Developmental Regulator Homothorax in Limiting Longevity in
The normal aging process is associated with stereotyped changes in gene expression, but the regulators responsible for these age-dependent changes are poorly understood. Using a novel genomics approach, we identified HOX co-factor unc-62 (Homothorax) as a developmental regulator that binds proximal to age-regulated genes and modulates lifespan. Although unc-62 is expressed in diverse tissues, its functions in the intestine play a particularly important role in modulating lifespan, as intestine-specific knockdown of unc-62 by RNAi increases lifespan. An alternatively-spliced, tissue-specific isoform of unc-62 is expressed exclusively in the intestine and declines with age. Through analysis of the downstream consequences of unc-62 knockdown, we identify multiple effects linked to aging. First, unc-62 RNAi decreases the expression of yolk proteins (vitellogenins) that aggregate in the body cavity in old age. Second, unc-62 RNAi results in a broad increase in expression of intestinal genes that typically decrease expression with age, suggesting that unc-62 activity balances intestinal resource allocation between yolk protein expression and fertility on the one hand and somatic functions on the other. Finally, in old age, the intestine shows increased expression of several aberrant genes; these UNC-62 targets are expressed predominantly in neuronal cells in developing animals, but surprisingly show increased expression in the intestine of old animals. Intestinal expression of some of these genes during aging is detrimental for longevity; notably, increased expression of insulin ins-7 limits lifespan by repressing activity of insulin pathway response factor DAF-16/FOXO in aged animals. These results illustrate how unc-62 regulation of intestinal gene expression is responsible for limiting lifespan during the normal aging process.
Published in the journal:
. PLoS Genet 9(2): e32767. doi:10.1371/journal.pgen.1003325
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1003325
Summary
The normal aging process is associated with stereotyped changes in gene expression, but the regulators responsible for these age-dependent changes are poorly understood. Using a novel genomics approach, we identified HOX co-factor unc-62 (Homothorax) as a developmental regulator that binds proximal to age-regulated genes and modulates lifespan. Although unc-62 is expressed in diverse tissues, its functions in the intestine play a particularly important role in modulating lifespan, as intestine-specific knockdown of unc-62 by RNAi increases lifespan. An alternatively-spliced, tissue-specific isoform of unc-62 is expressed exclusively in the intestine and declines with age. Through analysis of the downstream consequences of unc-62 knockdown, we identify multiple effects linked to aging. First, unc-62 RNAi decreases the expression of yolk proteins (vitellogenins) that aggregate in the body cavity in old age. Second, unc-62 RNAi results in a broad increase in expression of intestinal genes that typically decrease expression with age, suggesting that unc-62 activity balances intestinal resource allocation between yolk protein expression and fertility on the one hand and somatic functions on the other. Finally, in old age, the intestine shows increased expression of several aberrant genes; these UNC-62 targets are expressed predominantly in neuronal cells in developing animals, but surprisingly show increased expression in the intestine of old animals. Intestinal expression of some of these genes during aging is detrimental for longevity; notably, increased expression of insulin ins-7 limits lifespan by repressing activity of insulin pathway response factor DAF-16/FOXO in aged animals. These results illustrate how unc-62 regulation of intestinal gene expression is responsible for limiting lifespan during the normal aging process.
Introduction
The normal aging process involves the deterioration of a variety of tissues. In Caenorhabditis elegans, this includes loss of mobility and muscle cellular organization, decreased pharyngeal pumping, ectopic neuronal branching, deterioration of synapse function, and degeneration of the intestine [1], [2], [3], [4], [5]. The deterioration of the intestine is particularly striking, as the normal lumenal structure and even entire nuclei are lost in old worms [5]. One approach to understanding the intrinsic process of aging is to find mechanisms responsible for these changes that occur during normal aging [6].
Genome-wide microarray studies of aging in C. elegans have identified over a thousand transcripts that significantly increase or decrease in expression with age, covering a wide array of tissue-types and functions [7], [8], [9], [10]. Genes that are altered between young and old include not only downstream effects of aging, but more intriguingly may also reflect changes that are causal for tissue decline. Indeed, variability in the expression of a number of these age-regulated genes predicts remaining lifespan of individual worms, indicating that these genes can serve as biomarkers of physiological aging [11]. However, the genetic and regulatory mechanisms underlying these changes are poorly understood.
In an effort to link gene expression changes during aging to a regulator that mediates longevity, Budovskaya et al. identified a hypodermal developmental circuit involving GATA transcription factor genes elt-3, elt-5, and elt-6 as an example of a developmental circuit affecting lifespan [10]. Expression of elt-3 declined with age independently of tested stress and damage signals, and increased activity of elt-3 led to increased lifespan. These results suggest that mis-regulation of developmental pathways during aging, including the elt-3 pathway, may be an intrinsic property limiting lifespan of the organism [10]. There are many other transcription factors have been linked to longevity through the insulin signaling, dietary restriction, nutrient sensing, and stress-response pathways that may also be involved in orchestrating the normal aging process (reviewed in [12]).
Here we use a novel genomics approach based upon screening large-scale transcription factor target datasets generated by the modENCODE project to identify potential transcriptional regulators of age-dependent expression in C. elegans (E.V.N. and S.K.K. unpublished data). The most promising transcription factors identified by this approach include not only regulators involved in stress response pathways, but also developmental regulators. One such factor, UNC-62, is the C. elegans ortholog of Hox co-factor Homothorax (Hth/Meis). UNC-62 plays essential roles in development of the nervous system, hypodermis, and vulva [13], [14], [15], [16], but knockdown of unc-62 beginning in adulthood extends lifespan by ∼40% [17]. Focusing on the intestine, we identified that intestine-specific knockdown of UNC-62 recapitulates this extended lifespan and that UNC-62 decreases expression with age in the intestine. A number of mechanisms can be linked to this extension of lifespan, including shutting off yolk protein (vitellogenin) gene expression, preventing aberrant expression of neuronal genes in the intestine in old age, and increasing transcriptional resources for general intestinal gene expression. These mechanisms illustrate how antagonistic pleiotropy and allocation of cellular resources between somatic health and fertility can play important roles in specifying longevity.
Results
In order to uncover potential regulators of aging, we identified transcription factors with targets that are enriched for age-regulated genes. To identify targets bound by each transcription factor, we queried the data from 98 ChIP-seq datasets covering 57 transcription factors (many in multiple developmental stages) generated by the C. elegans modENCODE consortium [18]. For each ChIP seq experiment, we obtained DNA regions bound by the transcription factor identified by the PeakSeq algorithm (q-value<10−5) [19]. We noted that for each transcription factor, some binding sites are factor-specific (bound by only a few other transcription factors) whereas others are general (including the extreme case of HOT regions bound by more than >65% of transcription factors assayed)(Figure S1A–S1B) [20]. Compared to general binding sites, factor-specific binding sites are associated with genes whose expression and biological function are more likely to be shared with the transcription factor (E.V.N. and S.K.K. unpublished data; [20]). Thus, to improve the association between the presence of a binding site in ChIP-seq data and regulation of expression of nearby transcripts, we narrowed the targets to those that are factor-specific, defined as regions significantly enriched for a given transcription factor and up to 8 other transcription factors out of the 57 total transcription factors assayed by the modENCODE consortium(E.V.N. and S.K.K. unpublished data; [20]).
Each set of factor-specific targets was then independently compared to a set of 1106 age-dependent genes obtained from DNA microarray experiments [10], with significant enrichment in the overlap between the two sets identified using Fisher's exact test. At an enrichment cutoff of p<10−5, we identified 9 transcription factors as potential aging regulators (E.V.N. and S.K.K. unpublished data). These include transcription factors that have been previously associated with aging and longevity, such as oxidative damage response factor skn-1 [21] and developmental regulator elt-3 [10].
One transcription factor with known roles during embryonic and larval development, unc-62, has been shown to modulate lifespan [17] but was not previously linked to changes that occur during the normal aging process. UNC-62 is the C. elegans ortholog of Drosophila Homothorax and mammalian Meis, which are co-factors for HOX transcriptional regulators [22]. During development, knockdown of unc-62 activity yields phenotypes in maturation of the vulva, hypodermis, and the nervous system [13], [14], [15], [16]. However, during adulthood, reduction of unc-62 activity by RNAi extends lifespan by ∼30–40% (Figure 1A) [17]. Using the 1272 UNC-62 binding sites identified by the modENCODE consortium in day 4 young adult worms, we identified 399 factor-specific binding sites associated with 310 target genes. This set of direct targets includes 52 genes that show altered expression with age (2.9-fold enriched, p<10−15)(Figure 1B). In this work we characterize a new role for unc-62 in adults in order to explore the connection between an essential developmental regulator and aging.
The intestine-specific unc-62(7a) isoform decreases with age in the intestine
To investigate the role of unc-62 in specifying lifespan, we first set out to define the tissues in which it is expressed. Previous experiments using a 2.9 kb proximal promoter for unc-62 identified expression in vulval precursor cells, neurons, and some hypodermal cells [15]. Additionally, SAGE-tag sequencing from dissected intestines showed that unc-62 is expressed in the intestine [23]. In order to characterize the expression of unc-62 further, we made use of a GFP reporter for unc-62 in its full genomic context. This strain was generated by creating a modified fosmid with GFP inserted at the C-terminus of UNC-62, and includes regulatory elements contained within introns as well as distal promoter elements (Figure 1C). Biolistic bombardment was then used to obtain a strain expressing a low-copy, integrated UNC-62:GFP transgene. This fusion was sufficient to rescue embryonic and larval lethality phenotypes of the unc-62 deletion allele (s472)(Figure 1D).
Using this translational fusion, we observed strong UNC-62:GFP expression in a variety of tissues in the hermaphrodite, including the vulval precursor cells, the ventral cord motorneurons and other neurons, the hypodermis, and the intestine (Figure 1E). Expression in neurons and hypodermis is visible throughout development and into adulthood; expression in vulval precursor cells is only seen when that lineage emerges in the L3 and L4 stages (data not shown). For the intestine, we observed strong expression beginning in the L3 stage and continuing into adulthood (Figure 1E).
Alternative RNA splicing of unc-62 produces isoforms that differ in usage of the seventh exon (7a and 7b), which encodes the N-terminal region of the UNC-62 TALE homeodomain [13]. Transcripts including exon 7a increase in abundance more than 100-fold between embryos and adults, whereas transcripts containing 7b are expressed roughly equivalently throughout development [13]. To query the tissue-specific expression of these isoforms, we generated UNC-62:GFP reporters that express only UNC-62(7a) or UNC-62(7b). These reporters were made by inserting a stop codon in exon 7b into the UNC-62:GFP translational reporter such that it can only express UNC-62(7a):GFP, or a stop codon in exon 7a such that the reporter can only express UNC-62(7b):GFP (Figure 2A). Strains stably expressing either of these isoform-specific GFP reporters were then generated by biolistic bombardment. Fluorescent imaging of these strains indicated that UNC-62(7a) was predominantly expressed in the intestine starting in L3 and continuing through adulthood. In contrast, UNC-62(7b) was expressed in neurons, the ventral nerve cord, vulval precursor cells, and hypodermis beginning in embryos and continuing through adulthood (Figure 2A).
We next determined whether expression of the unc-62 isoforms change with age in the intestine. We measured expression of the intestinal unc-62(7a) isoform during normal aging, and observed that it decreases in the intestine by day 8 of adulthood (∼26% decrease, p<10−10)(Figure 2B–2C, Figure S2A). A similar decrease was observed in whole-worm unc-62(7a) mRNA levels between day 2 and day 8 of adulthood when measured by qPCR (Figure S2B). In contrast, unc-62(7b) was not observed in the intestine at any age. Thus, the decrease in expression of intestinal unc-62(7a) with age is not caused by changes in alternative splicing, but rather by a decrease in the level of transcription.
unc-62 acts in the intestine and hypodermis to modulate lifespan
As unc-62 is expressed in a variety of tissues, we wanted to identify the tissues in which unc-62 functions to affect lifespan of the entire organism. To do this, we performed RNAi knockdown of unc-62 in strains that generate an RNAi response only in specific tissues. These strains contain a mutation in the RNAi pathway (rde-1) that leads to inactivation of the RNAi response, as well as a transgene expressing rde-1 under a tissue-specific promoter to restore the RNAi response only in that tissue (Figure 3A) [24], [25]. We found that knockdown of unc-62 either only in intestine or only in the hypodermis was sufficient to cause a ∼30% extension of median lifespan (Figure 3C–3D). As a control, we crossed in the unc-62:GFP transgene and observed that GFP RNAi reduced GFP expression only in the proper tissue for the intestinal- and the hypodermal-specific RNAi strains (Figure S3). However, unc-62 RNAi did not have a significant effect on longevity in muscle tissue, neuronal tissue or the uterus (Figure 3E–3G). These results indicate that wild-type unc-62 activity in the intestine and hypodermis is critical for its role in limiting wild-type lifespan. In this work, we primarily focus on intestinal activity of unc-62 and the intestine-specific unc-62(7a) isoform.
Next, we asked at what time unc-62 acts to limit lifespan. Loss of unc-62 during development leads to severe developmental defects, indicating that unc-62 is beneficial up until adulthood. In contrast, unc-62 RNAi started at the first day of adulthood extends lifespan (Figure 1A). We found that unc-62 RNAi beginning at day five of adulthood could still significantly extend lifespan after the end of self-fertility (Figure S4A). This indicates that even after day five, unc-62 still performs functions that are detrimental and that limit the lifespan of wild-type worms.
Stage-specific binding of targets by the UNC-62 Homothorax transcription factor
In addition to the adult ChIP-seq targets described earlier, we analyzed UNC-62 ChIP-seq targets identified in the earlier L3 larval stage, when unc-62 appears predominantly as the 7b isoform (E.V.N. and S.K.K. unpublished data). In the L3 stage, the ChIP seq experiments identified 1193 UNC-62 binding sites, of which 251 are factor-specific. However, many young adult UNC-62 binding sites showed little enrichment in the L3 stage ChIP-seq experiment (Figure 4A and E.V.N. and S.K.K. unpublished data). We compared the 151 genes associated with these UNC-62 L3 binding sites to the set of aging-regulated genes and found only 10 UNC-62 L3 targets that are also age-dependent, which is not significantly enriched (.98-fold, p>.1). Thus, although the adult UNC-62 targets tend to be differentially expressed with age, the UNC-62 targets in L3 larvae are not. The genes bound by UNC-62 also have differential tissue-specificity: L3 larval targets are enriched for neuronal-enriched expression, whereas young adult targets are instead enriched for intestine-specific expression (E.V.N. and S.K.K. unpublished data). These distinctions corresponds to the pattern of alternative splicing of UNC-62; i.e., up until the L3 stage unc-62 is expressed predominantly as the unc-62(7b) isoform in neurons and the hypodermis, whereas in adults unc-62 appears as the unc-62(7a) isoform in the intestine and as the unc-62(7b) isoform elsewhere (E.V.N. and S.K.K. unpublished data). We focused our analysis on the young adult targets as these targets are enriched for genes with age-dependent expression.
UNC-62 is a direct and necessary activator of intestinal yolk protein genes
To understand the molecular effects of unc-62 knockdown in adults, we identified the targets of UNC-62 that are altered when unc-62 expression is knocked down. Poly(A)+ RNA was prepared from replicates of young adult hermaphrodites grown on either unc-62 RNAi or empty vector RNAi. We then prepared 3′ end-enriched RNA-seq libraries [26], with each sample barcoded to allow for multiplex sequencing. The triplicate libraries of empty vector and unc-62 RNAi were separately pooled and sequenced, yielding >14 million mapped reads for each sample (Figure 4B and Table S1).
To identify transcripts with consistently altered expression upon unc-62 RNAi treatment, we developed an approach based upon a Rank Product method used for DNA microarray analysis (Figure S5) [27]. Our analysis yielded 182 transcripts with altered expression upon unc-62 RNAi treatment at a false discovery rate of 10%, with 115 transcripts showing decreased expression and 67 showing increased expression upon unc-62 RNAi (Figure 4C and Dataset S1). Of these 182 unc-62 dependent transcripts, only ten are directly bound by UNC-62 in young adults (factor-specific binding sites with q-value<10−5) indicating that the remainder are indirect targets (Table S3).
We found that three of the eight genes that are bound and activated by UNC-62 are vitellogenin genes. C. elegans contains 6 vitellogenin transcripts (vit-1 through vit-6) that encode the major yolk proteins [28]. These transcripts are specifically expressed in the intestine during the reproductive period, and vitellogenin proteins are exported out of the intestine into oocytes [29]. RNA-seq experiments indicated that expression of all six vitellogenin transcripts is strongly dependent on unc-62 activity, as there is a 4–10-fold decline in expression for vit-2, vit-3, vit-4, vit-5 and vit-6 and a ∼100-fold decline for vit-1 in unc-62 RNAi worms as compared to controls (Figure 5A). unc-62 RNAi also results in a strong decrease in the level of expression of a VIT-2:GFP translational reporter, indicating that protein as well as mRNA levels of vit-2 are dependent upon UNC-62 activity (Figure 5B). RNAi targeting exon 7a of unc-62 results in a decrease of VIT-2:GFP expression, whereas RNAi targeting exon 7b does not, indicating that the unc-62(7a) isoform activates vitellogenin expression (Figure 5B).
In addition to the significant association observed at three of the six vitellogenin loci in the initial ChIP-seq analysis (q-value<10−5)(Table S3), we found that the remaining three vitellogenin loci were bound by UNC-62 in adults with q-value<0.01 (Figure 5A). Interestingly, interaction with vitellogenin promoters is not observed in ChIP-seq of UNC-62 in L3 stage worms, indicating that UNC-62 is bound to the vitellogenin promoters in adults but not young larvae. We validated binding of UNC-62 to the promoter region of the vit-2 and vit-6 genes in adults by ChIP-qPCR, observing 5- to 15-fold enrichment of UNC-62 binding to these regions compared to control promoter regions (Figure S6).
Previous experiments have identified two sequence motifs that are critical sequence elements shared among all 6 vitellogenin promoters: GATA-motif VPE-2 (CTGATAA) and VPE-1 (ATTGACA) [30]. The VPE-1 motif contains the core TGACA sequence bound by Homothorax in Drosophila [31]. We used a de novo motif search among all UNC-62 young adult binding sites and found that a motif with consensus sequence ATTGACA was the most significantly enriched motif, and that the seven base ATTGACA sequence itself was 4.3-fold enriched in UNC-62 binding sites (Figure 5C, Figure S7). The finding that UNC-62 is a direct and necessary regulator of all six vitellogenin genes, combined with the fact that the VPE-1 motif perfectly matches the UNC-62 motif identified from a de novo motif search, suggests that UNC-62 directly binds to the VPE-1 element and activates vitellogenin expression.
The six vitellogenin genes are among the genes that show the greatest decrease in expression with age in the entire C. elegans genome, with a 16-fold decrease in expression between young (day 2) and old (day 11) adulthood (Figure 5D) [10], [32]. Having observed that UNC-62 binds to and activates the vitellogenin genes in young adults, we next asked whether residual levels of vitellogenin expression in old adults were still dependent on UNC-62 activity. We micro-dissected the intestine from day 15 adult worms fed either control or unc-62 RNAi bacteria, and performed qRT-PCR on RNA isolated from the intestines to measure levels of vit-2 and vit-6 expression. Despite expression levels decreasing with age, vit-2 and vit-6 expression declined further when unc-62 was knocked down in old adults (Figure 5E), consistent with residual vitellogenin expression in old age remaining activated by unc-62.
Indirect targets of unc-62 are enriched for genes altered during aging
To further understand the roles of unc-62 during aging, we analyzed the 182 transcripts that show altered expression in unc-62 RNAi-treated worms, consisting of 67 transcripts that increase expression and 115 transcripts that decrease expression. The 115 targets that are (directly or indirectly) activated by wild-type unc-62 show a strong enrichment for genes specifically expressed in the hypodermis (2.6-fold; p<10−5). Notably, 41 of these 115 genes are collagen genes, out of a total of 90 collagen genes quantified in the RNA-seq experiment (29-fold enriched, p<10−15)(Figure 6A). These unc-62-activated genes also include 38 that decrease in expression during normal aging (a 3.9-fold enrichment, p<10−20) [10], including 26 of the 41 unc-62-activated collagen genes (Figure 6A–6B). These findings suggest that changes in gene expression of hypodermal genes upon unc-62 knockdown (particularly, collagen genes) are similar to changes normally observed as adult worms age; i.e., unc-62 RNAi mimics the normal aging process in the hypodermis.
Next, we considered the genes with increased expression upon unc-62 RNAi. At the tissue-level, we found that 32 of these genes show intestine-enriched expression in one or more datasets (2.6-fold enriched, p<10−10). These genes that are repressed by wild-type unc-62 activity (although not bound by UNC-62 directly) also show a 3.2-fold enrichment for genes that decrease in expression with age (p<10−4)(Figure 6C). Thus, unc-62 RNAi indirectly leads to increased expression of genes that normally decline with age and are primarily expressed in the intestine.
A possible explanation for the overall increase in gene expression in the intestine in unc-62 RNAi worms could be that it is an indirect consequence of shutting down vitellogenin gene expression. The vitellogenin transcripts are among the most highly-expressed transcripts in C. elegans; vit-4 and vit-6 are the sole non-ribosomal protein genes among the ten highest expressed genes in young adults, and the six vitellogenin genes comprise ∼3% of all mapped reads in young adults in RNA seq experiments despite being expressed solely in intestinal cells (Table S1). Thus, as an indirect consequence of dramatically reducing vitellogenin gene expression by unc-62 RNAi, transcriptional machinery that was previously allocated to express vitellogenin genes may become available to activate transcription of other genes.
This model posits that the decrease in vitellogenin transcription would lead to a generalized increase of the remaining genes expressed in the intestine. To test this, we identified 1699 genes that are intestine-enriched [23], [33], [34] as well as the subset of 291 that decrease in expression with age [10]. We observed that intestinal genes in general, and particularly those that decrease expression with age, are significantly shifted towards increased expression upon unc-62 RNAi (both p<10−20 by Kolmogorov-Smirnov test)(Figure 6D). In contrast, there was no shift observed either for neuron-enriched or muscle-enriched genes (p>0.01 by Kolmogorov-Smirnov test)(Figure S8). Thus, in addition to the 67 genes that were identified to be significantly altered upon unc-62 RNAi, there is a general increase in expression of intestinal genes in unc-62 RNAi worms that may provide an indirect benefit contributing to longer life. As intestinal gene expression generally declines with age, this broad effect of gene activation by unc-62 RNAi opposes the normal aging process in the intestine.
The above findings suggest that unc-62 knockdown may delay the age-related decay of the intestine. In order to assay the morphology of the intestine, we examined intestinal expression of an opt-2:GFP reporter. opt-2 encodes an intestinal low affinity/high capacity oligopeptide transporter that localizes to the apical membrane of the intestine, and is visible as a single tubular structure tracking with the inner intestinal membrane in young adults (Figure 6E) [35]. During normal aging, this pattern rapidly degrades by day three of age (long before worms begin to die from old age). In unc-62 RNAi worms, the age-related deterioration of intestinal morphology (measured by presence of apical expression of opt-2:GFP) is delayed (Figure 6E). In summary, knockdown of unc-62 increases lifespan, increases the period of healthy intestine morphology, and ameliorates intestinal expression changes that occur during wild-type aging.
Expression of a new class of UNC-62 targets in the intestine in old age
As part of our analysis, we discovered a set of genes bound by UNC-62 in young adults that, during development, are enriched for expression in neuronal cells as opposed to the intestine. We next examined whether expression in the intestine of these targets changes during normal aging, and whether their expression is dependent on unc-62 activity.
To profile expression changes with age in the intestine, we micro-dissected young and old intestines, isolated RNA, and performed qPCR to assay expression of these UNC-62 targets. In these samples, control transcripts known to be specifically expressed in the intestine showed a greater than 10-fold enrichment compared to non-intestinal transcripts (Figure S9). We randomly selected eight UNC-62 young adult targets (that we refer to as Class N) that have neuronal-enriched expression patterns but had significant UNC-62 binding in adults by ChIP-seq, allowing for the possibility that their expression is under UNC-62 control during normal adult aging. As negative controls, we used ten genes with neuronal-enriched expression that were bound by UNC-62 in early development, but no longer identified as significantly bound in adults.
We observed that expression of the eight UNC-62 Class N young adult targets in dissected intestines showed a significant shift towards increased expression in old age (p = 0.013 by Wilcoxon rank-sum test)(Figure 7A). Specifically, five of the eight adult targets increased more than two-fold in expression with age (including one gene (ins-7) which has previous been shown to increase with age in the intestine [36]), whereas none of the ten control targets did so. To determine whether expression of class N genes was dependent on UNC-62 activity, we obtained RNA from micro-dissected intestines from old adults grown on unc-62 RNAi. We found that four (tsp-1, max-1, ins-7, and T05B11.1) of the five class N genes that increase expression by two-fold in intestines with age showed decreased expression upon unc-62 RNAi (Figure 7B). In summary, class N genes are an intriguing set of unc-62-dependent genes that are expressed predominantly in neuronal tissues in development, but then become expressed in the intestine in old adults.
Late-life expression of at least two class N genes (max-1 and ins-7) appears to be detrimental, as RNAi treatment or loss-of-function mutations in these two genes have been shown to extend lifespan. The first, max-1, encodes a PH/MyTH4/FERM domain-containing protein that is required for proper motor axon projections [37], and max-1 RNAi treatment increases lifespan by ∼5% [38]. The second, ins-7, encodes an insulin/IGF-1-like peptide that acts as an agonist for the insulin/IGF-1 receptor DAF-2 [39]. The insulin/IGF-1 signaling pathway ultimately represses the activity of the FOXO-family transcription factor DAF-16. When the activity of the insulin/IGF-1 signaling pathway is low or off, DAF-16/FOXO localizes to the nucleus and activates expression of beneficial genes that extend lifespan [12]. Knockdown of insulin ins-7 activity throughout the worm results in activation of DAF-16/FOXO and extends lifespan, and intestine-specific over-expression of ins-7 in young adults is toxic [36]. We verified that knockdown of ins-7 expression specifically in the intestine can extend lifespan by showing that intestine-specific RNAi of ins-7 results in an 8.3% extension of mean lifespan (p = 7.8×10−8 by log-rank test)(Figure S4B).
One possible mechanism contributing to the lifespan extension seen in unc-62 RNAi worms is that ins-7 expression does not increase, but instead remains at a low level in the intestine in old age. As a result, DAF-16/FOXO activity would not be repressed by ins-7, and could instead remain high in unc-62 RNAi worms as opposed to being repressed in normal aged worms. According to this possibility, unc-62 RNAi should result in activation of DAF-16/FOXO activity and the longevity phenotype should be at least partially dependent on daf-16 activity. First, we found that unc-62 RNAi results in increased expression of a DAF-16 reporter (sod-3:GFP [40]) in aged animals, suggesting an increase of DAF-16/FOXO activity (Figure 7C). As has been previously observed, increased DAF-16/FOXO activity upon unc-62 knockdown was not observed in young adults (when ins-7 is not yet induced in the intestine) [17]. Second, we found that a daf-16 mutation suppresses the longevity caused by unc-62 RNAi, indicating that the benefit conferred by unc-62 knockdown requires daf-16 activity (Figure 7D). These results indicate that extended longevity conferred by unc-62 RNAi involves preventing ins-7 expression in old age, which results in activation of the insulin signaling pathway and repression of DAF-16/FOXO.
unc-62 RNAi extends the lifespan of glp-1 mutants
It has previously been shown that germline signaling can affect aging, as glp-1(e2141) mutants that lack the germline show extended lifespan [41]. We next asked whether the beneficial effects of knockdown of unc-62 were independent of those from removing a germline by comparing the lifespan of an unc-62 RNAi;glp-1(e2141) double mutant to the lifespan of glp-1(e2141) single mutants. We found that unc-62 RNAi increased the mean lifespan of glp-1(e2141) mutants by 28% (p<10−5), indicating that unc-62 activity is still detrimental for lifespan in worms lacking a germline (Figure S4C). Previous studies have shown that vitellogenins continue to be expressed in animals lacking a germline [42], [43]. Thus, the lifespan extension by unc-62 RNAi in a glp-1(e2141) background may reflect the prevention of vitellogenin accumulation. It may also reflect other downstream effects of unc-62, including increased availability of transcriptional resources for intestinal expression and decreased expression of class N genes in old worms.
Discussion
There are a large number of datasets generated by the modENCODE project that profile targets of transcription factors by ChIP-seq. To identify putative aging regulators, we developed an approach in which we query all transcription factors for those with directly-bound targets that are enriched for altered expression during normal aging. In addition to transcription factors ELT-3 and SKN-1 that have previously been characterized as regulators of aging that can modulate lifespan [10], [21], this approach identified factors that had not previously been linked to gene expression changes that occur during normal aging. The identification of HOX co-factor UNC-62 (Homothorax/hth/Meis) was particularly interesting, as knockdown of unc-62 in adults significantly extends lifespan [17] but unc-62 activity is required for embryonic and larval development [13], [14], [15], [16].
We found that unc-62 acts in the intestine to limit lifespan and that intestinal unc-62(7a) expression decreases with age. These results suggest that the age-dependent decrease of unc-62 in the intestine is likely beneficial rather than detrimental for lifespan. With regards to unc-62 expression, the young state is thus not a model for health as compared to the old state, as this would predict that RNAi of unc-62 in young worms would shorten their lifespan. In order to better understand why RNAi of unc-62 results in longer lifespan, we applied a variety of genomics analyses to explore downstream targets of unc-62 in young and old adults. We identified three downstream effects of unc-62 activity that can be linked to the aging process: 1) expression of yolk proteins that accumulate with age, 2) generalized effects on transcription of a large number of intestinal genes, and 3) activation of non-intestinal genes in the intestine in old age. The detrimental effects of unc-62 in adult worms may reflect a combination of these downstream effects, each partially contributing to the overall limitation of lifespan.
Intestinal UNC-62 and target yolk proteins decline with age
The first mechanism that we identify linking unc-62 to aging is through direct activation of vitellogenin genes in the intestine. We find that UNC-62 binds to and activates all six vitellogenin genes, and that the decrease in expression of unc-62(7a) with age parallels the decline in intestinal expression of vitellogenin genes. These results suggest that the decline with age of UNC-62(7a) contributes to decreased expression of vitellogenins with age. In addition to unc-62/Homothorax, we note that it is possible that the age-dependent decline of vitellogenin gene expression may also be due to decreased activity of the GATA transcription factor ELT-2. ELT-2 is a master regulator of intestinal gene expression and also regulates vitellogenin expression through a distinct VPE-2 promoter element [23].
The six vitellogenin genes encode the C. elegans yolk proteins, which are utilized for growth and development of progeny. The C. elegans vitellogenin proteins appear to be toxic in old age, as knockdown of vit-2 and vit-5 by RNAi results in increased lifespan [39]. During the self-fertile reproductive phase of hermaphrodites (through day five of adulthood [44]), the vitellogenin proteins are removed from the mother through egg-laying. After self-fertile reproduction ends at day five, vitellogenins remain necessary for reproduction by cross-fertilization, which occurs up to day 13 [45]. In the absence of cross-fertilization, vitellogenins accumulate in the body cavity, are subjected to oxidative damage and ultimately become one of the most prevalent proteins in aggregates late in life [5], [46], [47]. Similar to C. elegans, vitellogenin proteins in Drosophila form aggregates and accumulate oxidative damage in old age [48]. Thus, decreased vitellogenin protein accumulation in the body cavity presents one downstream effect of unc-62 knockdown that may contribute to extended longevity.
A second mechanism that may link unc-62 RNAi to the aging process is through allocation of intestinal transcriptional resources. We observe a general increase in expression of intestinal genes upon unc-62 RNAi. Previous microarray studies have indicated that intestinal genes show a general decrease in expression with age [10], [49], and unc-62 RNAi has a particularly strong effect on these age-regulated intestinal genes. Thus, while unc-62(+) is necessary for intestinal expression of vitellogenins that are produced to provide food for developing embryos, it leads to repression of other intestinal genes. This is reminiscent of the disposable soma theory of aging, in which resource allocation between the soma and the germ-line provides a balance between maternal health and the health and number of progeny [50].
In contrast to the vitellogenins that are directly bound and activated by UNC-62, most of the unc-62-repressed intestinal genes are not directly bound by UNC-62. One possible mechanism for how unc-62 could indirectly repress these genes involves transcriptional resources in the intestine. As vitellogenins are among the most highly-expressed transcripts in the worm genome, silencing of vitellogenins upon unc-62 RNAi would increase the availability of transcriptional machinery for expression of other genes in the intestine. In this model, unc-62 would provide a molecular mechanism through which the hermaphrodite intestine allocates transcriptional resources for vitellogenin expression to provide for progeny at the expense of general expression of other intestinal genes that help to maintain somatic health of the mother. The general increase in intestinal expression would not be expected upon RNAi of the vitellogenins themselves, as RNAi would not block vitellogenin gene transcription but rather silence protein production through degradation of vitellogenin mRNA. We do not yet know whether the general increase in intestinal expression that we observe upon unc-62 RNAi involves only the increased availability of general transcriptional resources, or whether it also requires the activity of specific intestinal regulators (e.g. other transcription factors such as ELT-2).
In addition to vitellogenins, there are many other genes that are bound by UNC-62 in young adults that decrease expression with age in the intestine. However, it is not clear whether age-related decline of these other targets is caused by decreased levels of unc-62(7a) or by altered activity of other transcriptional regulators. Using RNA-seq, we found that few of these direct target genes showed changes in expression in unc-62 RNAi worms, suggesting either that UNC-62 does not regulate their expression or that our RNA-seq method was not sensitive enough to detect changes in their expression. This apparent incongruity between transcription factor binding and either activation or repression is not unique to UNC-62, as many other transcription factors in various species have shown a low overlap between genes that are directly bound and whose expression is responsive to over-expression or knockdown of that transcription factor [51], [52], [53]. Further studies will be required to explore whether UNC-62 binding plays a role in age-regulation of these other genes.
Although our analysis largely focused on the intestine, genetic experiments indicate that activity of unc-62 in the hypodermis may be important for specifying lifespan. Reduction of unc-62 activity specifically in the hypodermis extends lifespan, but the targets that are directly activated by UNC-62 in the hypodermis are not known. However, the indirect targets of UNC-62 generally decrease expression with age and include a substantial fraction of hypodermal collagen genes. It remains unclear whether altered collagen expression with age is due to altered activity of unc-62 itself or additional co-factors.
An insulin-like gene increases expression and is activated by unc-62 in the aged intestine
A third way that unc-62 can modulate longevity is via its role in the aberrant expression of genes in the intestine in old animals. Our analysis identified a set of five class N UNC-62 target genes that were expressed primarily in neuronal tissues during development, but that increased expression in the intestine as animals grow old. One such gene that is decreased in expression in aged intestines upon unc-62 RNAi treatment is known to contribute to extended longevity, as insulin INS-7 is linked to aging through its role as a signaling molecule in the insulin/IGF signaling pathway [39]. ins-7 is predominantly expressed in neuronal tissues during development, but increases expression in the intestine with age [36]. In the insulin signaling pathway, insulin binding to the DAF-2 insulin/IGF-1 receptor turns on a signaling cascade that represses the activity of the DAF-16/FOXO transcription factor. ins-7 has been shown to repress the activity of the insulin signaling pathway, as knockdown of ins-7 increases lifespan as well as nuclear localization and activity of DAF-16 [36]. We show that ins-7 is directly activated by UNC-62 in old adult intestines, that intestinal ins-7 limits lifespan, and that unc-62 RNAi represses ins-7 and leads to increased DAF-16 activity in old age. These results indicate that knockdown of unc-62 allows DAF-16 to remain active and to benefit aged animals. The results showing that unc-62 is required for expression of ins-7 in old worms are particularly intriguing in light of the previous finding that intestinal ins-7 not only represses DAF-16, but that DAF-16 also represses ins-7 [36]. Thus, the increase of ins-7 with age appears to trigger a double negative feedback loop that would be prevented by the knockdown of ins-7 upon unc-62 RNAi.
If activation of DAF-16 by unc-62 RNAi in old worms is a major cause for extended longevity, then mutations in daf-16 should either partially or completely suppress the longevity phenotype of unc-62 RNAi. Our experiments indicated that daf-16(m26) fully suppresses the longevity phenotype of unc-62 RNAi animals (Figure 7D). However, a previous study by Curran, et al. (2007) using a mutation (eri-1(mg366)) that increases sensitivity to RNAi and observed a weaker suppression; in this study, an eri-1(mg366);daf-16(mgDf47);unc-62 RNAi triple mutant lived 31% longer than the eri-1(mg366);daf-16(mgDf47) double mutant, but 42% shorter than the eri-1(mg366);unc-62 RNAi double mutant [17]. It is possible that this discrepancy is due to the use of the eri-1(mg366) mutation, as this mutation enables RNAi to affect neuronal tissues that are otherwise insensitive. If so, this could suggest an additional role for unc-62 mediating lifespan in neuronal tissues.
unc-62 activates four of the five class N genes that increase expression in the intestine in old age. However, it is unclear whether the increase in expression of these genes with age is caused by changes in UNC-62 activity per se. Neither of the unc-62 splicing isoforms appear to increase expression with age in the intestine; unc-62(7a) decreases with age, and we did not observe expression of unc-62(7b) (or any other unc-62 isoforms) in the intestine. Thus, it appears that expression of these genes requires unc-62 activity, but that their increase in expression with age may be driven by an additional regulator.
The eight class N genes that we tested in this work are a small amount of the total such genes in the genome. Therefore, it is likely that there are many more class N targets with similar behavior that could also contribute to modulation of lifespan by unc-62. These results suggest that the old intestine is not simply altered by accumulation of damage or decreased expression of intestinal genes. Instead, the old intestine also appears to differ from the young intestine by an increased expression of class N genes. Because their expression appears to be detrimental, it seems unlikely that age-dependent expression of class N genes is evolutionarily selected. Rather, the onset of class N gene expression in old worms may occur at an age that is beyond the force of natural selection, when there may not be sufficient selective pressure in aged animals to suppress deleterious changes in gene networks.
Links between development and aging
Recent evidence indicates that transcriptional changes during aging are not only a consequence of damage accumulation, but also reflect altered activity of developmental regulators. Analysis of aging of the human pre-frontal cortex suggests that the majority of age-dependent changes in expression are not unique to aging, but rather mirror both reversals and extensions of developmental expression patterns [54]. In mice, an age-dependent increase in Wnt signaling is associated with a shift from myogenic to fibrogenic lineage phenotypes among mouse muscle stem cells, which is abrogated by the addition of Wnt antagonists [55]. In C. elegans, a regulatory circuit in the hypodermal tissue involving three GATA transcription factor genes, elt-3, elt-5, and elt-6, appears to not only modulate lifespan but is also linked to altered expression of target genes with age [10]. Expression of the elt-5 and elt-6 GATA transcription factors increase with age, leading to a decrease in expression of the elt-3 GATA factor. Lifespan is increased in mutants in which elt-3 expression remains high in old age, suggesting that the age-dependent decrease of elt-3 expression is detrimental to longevity. As the elt-5, elt-6, and elt-3 circuit also plays a key role in specifying the hypodermal fate during development [56], these results suggest that mis-regulation of a key developmental pathway in the hypodermis during aging underlies intrinsic aging of the organism [10].
UNC-62 Homothorax/Meis provides a further example that developmental and reproductive activities of certain transcription factors are not stable in adults, and are not optimized for longer lifespan. To ameliorate or even reverse the aging process, it is not enough to simply prevent damage accumulation. Rather, it may also be necessary to correct developmental regulators responsible for intrinsic changes in the aging transcriptome that limit lifespan.
Materials and Methods
ChIP–seq data and analysis
Transcription factor (TF) target identification from chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) data was obtained from the modENCODE consortium (http://submit.modencode.org/submit/public/list or http://data.modencode.org/) [18], [20], using datasets available as of 04/30/12. TF Binding sites as well as read density files were mapped to WS220 coordinates using scripts available from WormBase (ftp://ftp.sanger.ac.uk/pub2/wormbase/software/Remap-between-versions/) [57].
All 98 ChIP-seq experiments (profiling targets of 57 transcription factors in one or multiple developmental stages) that contained at least 100 binding sites meeting a q-value<10−5 significance cutoff were used. ‘Factor-specific’ binding sites of a transcription factor were defined as those that were significantly enriched in the given factor, and had no position within the binding site that was enriched in 9 or more other transcription factors (out of the 57 total) (Figure S1A–S1B). Binding sites were associated with all transcripts (obtained from Wormbase release WS220) if the point of maximal read density within the peak occurred within the gene body, or was less than 3 kb upstream of the annotated transcription start site (for cases where the upstream intergenic distance was less than 3 kb, only the region up to the neighboring gene was used). Genes significantly altered in expression with age (p-value<10−4) were obtained from [10], with updated mapping of probes to transcripts (http://nemates.org/MA/primers/new_primer_names.html).
Significance of overlap between ChIP-seq targets and various gene lists was determined by a Fisher's Exact test on the 2×2 contingency table using the R statistics program (version 2.11.1), approximated with a Chi-Square test for datasets with expected and observed overlaps greater than 5. P-values for Chi-Square tests were obtained using the Perl Statistics::Distributions module.
ChIP–qPCR validation
To validate UNC-62 association with vit-2 and vit-6 promoters in adults, ChIP-qPCR was performed on independently generated adult lysates from strain OP600. ChIP was performed using the protocol and GFP antibodies generated by the modENCODE project [18], except for the replacement of sepharose beads with Invitrogen Protein G Dynabeads. qPCR primers were generated that flanked the peak of read density from UNC-62 YA ChIP-seq for the vit-2 and vit-6 promoter regions, with three regions bound by hypodermal transcription factor ELT-3 (B0222.8) and pharyngeal and intestinal TF PHA-4 (F32A5.4 and F35G2.1) serving as negative controls. For each, enrichment was determined relative to a non-immunoprecipitated input sample.
To identify sequence motifs in UNC-62 binding sites, the 100 base window surrounding the position of maximum ChIP-seq read density in each UNC-62 binding site was identified as the core binding region. To compare these regions against promoter sequences associated with genes with similar expression and tissue-specificities, 200 bp windows on either side of this core region served as the negative control. A de novo motif search for enriched 7-mers was performed using the RSAT program [58], with motif logos generated using WebLogo [59]. To validate the motif, frequencies of the ATTGACA 7-mer sequence were calculated for core and flanking regions, and enrichment calculated by Fisher's Exact test. As an additional control, this analysis was repeated using Fisher-Yates-shuffled core regions as the negative control, yielding similar enrichments.
RNA–seq profiling of unc-62 knockdown
RNA-seq quantification upon unc-62 RNAi was performed using the sterile TJ1060 (fem-15(b26); spe-9(hc88)) strain that has previously been used for transcriptome studies in adult worms in order to avoid isolating mRNA from developing embryos [10], [60]. Three independent biological replicate batches of approximately ∼10,000 sterile hermaphrodites were grown at the non-permissive (25°C) temperature until the first day of adulthood, at which point half were placed on plates seeded with unc-62 RNAi and half on empty vector RNAi. The worms were then grown for 3 days at 20°C, after which RNA was isolated by Trizol extraction followed by phenol-chloroform purification. Poly-A-purified RNA was then isolated using the Qiagen Oligotex Mini kit, and 3′ end enriched RNA-seq libraries were prepared in accordance with the 3′SEQ method. 400 ng of mRNA was heat-sheared for 7.5 minutes at 85°C to obtain 200–500 bp fragments, and barcoded oligo-dT primers were used to generate cDNA for mRNA fragments located at the 3′ end of transcripts. After second strand synthesis, adaptors were ligated to the 5′ end, and PCR primers (including the barcode sequence) were used in order to enable multiplexed sequencing. The three replicate libraries from worms grown on unc-62 RNAi were then pooled and sequenced in a single flow-cell lane on the Illumina single-end 36 bp HiSeq 2000 platform, with the three libraries from worms grown on control RNAi pooled and sequenced in an additional lane.
A total of ∼71 million and ∼73 million post-filtering reads were obtained for control and unc-62 RNAi respectively, with roughly equal proportions coming from the three biological replicates (Table S1). Reads with 10 or more consecutive A's or T's were removed as likely amplification artifacts, and reads were then mapped to the Wormbase reference WS215 genome and transcriptome using Bowtie version 0.12.7 [61], using the options “-a –best –strata -v 2” to obtain all genomic position(s) to which the read mapped with the least number mismatches (up to 2 allowed). Greater than 14 million reads for each independent replicate were obtained that mapped uniquely to the sense strand of WS215 annotated transcripts (with more than 85% of post-filter reads mapping to C. elegans in each experiment)(Table S1). To determine expression for each transcript in the WS215 Wormbase annotation, we counted the number of reads that uniquely mapped to the sense strand (i.e., they did not map equally well elsewhere in the genome). For genes with multiple annotated isoforms, only the transcript with the highest average expression across the six experiments was used for further analysis.
RNA-seq datasets have been submitted to the Gene Expression Omnibus under series GSE39574 (samples GSM971870-GSM971875).
Significance analysis by rank products
Identification of transcripts with consistently altered expression in all three biological replicates was patterned on the Rank Product method previously described for microarray analysis [27]. To avoid technical noise among transcripts with low expression, the subset of 7333 genes for which 10 or more sequencing reads were observed either in all three empty vector libraries, or all three unc-62 RNAi libraries, was obtained, and a pseudocount of one read was added to each gene in each condition. For each of the three biological replicates, genes were then ranked based on the observed expression fold-change between unc-62 RNAi and control. The rank product for each gene was then calculated as the product of that gene's ranks in each of the three replicate experiments. Repeating this approach on10,000 random permutations of the ranks for each of the three experiments converted these rank product values to E-values and Percent False Positive (PFP) values as previously described [27]. At a PFP cutoff of 10%, 116 genes with decreased expression and 69 genes with increased expression upon unc-62 RNAi were obtained. Three genes that met this false-positive criteria did not change consistently across the three biological replicates (i.e., were not in the top third of fold-changes in each replicate) and were discarded, yielding 115 genes decreased and 67 genes increased upon unc-62 RNAi.
Analysis of genes altered in expression upon unc-62 RNAi
Significantly increased and decreased genes upon unc-62 RNAi were compared against various gene lists as described. For all analyses, genes were discarded if they were not present in the 7333 genes with reliable expression in RNA-seq experiments, or were not assayed in the published dataset (as they were either absent on the microarray or were not yet annotated). For analysis of genes that decrease upon unc-62 RNAi, a dataset of genes with enriched expression in hypodermis (during L3-L4 larval stages) was used [34], and collagen genes were annotated as those belonging to NCBI KOG3544. For analysis of genes that increase upon unc-62 RNAi, intestine-enriched genes were defined as those with significantly enriched expression in intestine relative to whole-worm controls from experiments profiling intestinal gene expression in late embryos [34], L2 larvae [34], L4 larvae [33], or adults [23]. Genes that decrease expression with age were obtained from Budovskaya, et al. [10]. Kolmogorov-Smirnov tests were performed in MATLAB (version R2010b).
Intestine dissection
Intestine dissections were performed as described by McGhee, et al. (2007) [23] with slight modifications. Adult worms were placed in ∼10 uL of buffer containing 100 µl of PBS–EDTA–ATA (125 mM NaCl, 16.6 mM Na2HPO4, 8.4 mM NaH2PO4, 0.1 mM EDTA, 1 mM aurin tricarboxylic acid), supplemented with 0.1 uL of Superasin. 27G1/2 needles were used to cut below the pharyngeal bulb, and intestines were allowed to extrude for ∼2–3 minutes. Single intestines were then isolated with a pulled capillary tube, washed 3 times with PBS-EDTA-ATA-Superasin, and placed in Buffer RLT Plus + Beta-mercaptoethanol (Qiagen RNeasy Plus Micro Kit). RNA was isolated from batches of approximately 100 intestines using the RNeasy Plus Micro Kit. RNA quality for one young and one old sample was validated using the RNA 6000 Pico assay on the Bioanalyzer 2100 (Agilent), obtaining RNA integrity numbers of 10 and 9.6 respectively.
To identify class N (neuronal-enriched) young adult UNC-62 targets, we obtained twenty lists of genes with specific or enriched expression in either specific neuronal subtypes or in all neurons [34], [62], [63], [64], [65]. To determine whether UNC-62 binding in young adults correlated with altered expression in intestines with age, a negative control set of genes with identical expression criteria were identified that were associated with UNC-62 binding sites in L2 or L3 but not YA. The initial experiment compared a pair of samples composed of ∼100 intestines from young (day 1) or old (day 15) worms. Transcripts that were undetectable or did not show reliable signal across qPCR technical replicates in young intestines were discarded, as it cannot be determined whether these represented failed primer pairs or if these transcripts are simply never expressed in intestines at levels reliably detected by qPCR. Screening of 9 class N and 31 control genes yielded 7 class N and 6 control genes that were reliably detected in both young and old intestines, with an additional 1 class N and 4 control genes that (presumably due to age down-regulation) were detected in young but not old intestines. For these 5, a Ct value of 40 was used as an upper bound estimate of old adult expression. For each, expression was calculated relative to previously described control transcript tbb-2 (identical results were observed using alternative control let-70) [66]. Normalization was performed as previously described [67].
Further experiments to test the effect of unc-62 RNAi were performed on two biological replicates of ∼100 intestines for young control (day 1), old control (day 15 grown on empty vector RNAi), and old unc-62 knockdown (day 15 grown on unc-62 RNAi) respectively. For old samples, worms were grown on normal (OP50) E. coli until day 1 of adulthood, and then shifted to RNAi bacteria (as described for lifespan experiments). For analysis using these samples, transcript abundance was calculated using a 256-fold dilution curve performed on cDNA prepared from whole-worms. Technical triplicates of each qPCR reaction were performed.
Strains
To generate the unc-62:GFP translational reporter, fosmid WRM061dC01 was obtained from Geneservice, and an eGFP tag was inserted at the C-terminus of unc-62 by recombineering (Figure 1C) [68]. This reporter was integrated by biolistic bombardment, yielding strain OP600 (unc-119(ed3);Is[unc-62:GFP;unc-119(+)]. In order to better visualize intestinal expression, this strain was crossed with the glo-4(ok623) strain that lacks gut autofluorescence, and Is[unc-62:GFP;unc-119(+)];glo-4(ok623) worms were isolated in the F2 generation (strain SD1897).
To generate isoform-specific reporters for unc-62(7a) and unc-62(7b), a single base in either exon 7a or 7b was mutated to generate an in-frame stop codon (Figure 2A). This approach required the insertion of a Kanamycin cassette (on the reverse strand) in the constitutively spliced intron between exons 8 and 9 to aid selection for the mutated variant. Biolistic bombardment was used to integrate these fosmids, yielding Is[unc-62(7a):GFP;unc-119(+)];unc-119(ed3) (strain OP601) and Is[unc-62(7b):GFP;unc-119(+)];unc-119(ed3) (strain OP602) worms. These strains were similarly crossed to glo-4(ok623) worms to generate Is[unc-62(7a):GFP;unc-119(+)];glo-4(ok623) (strain SD1890) and Is[unc-62(7b):GFP;glo-4(ok623) (strain SD1898) which were used for fluorescence measurements.
To determine whether the full-length translational fusion was functional, unc-119(ed3); Is[unc-62:GFP;unc-119(+)] was crossed into unc-62(s472) animals, and unc-62(s472); Is[unc-62:GFP;unc-119(+)] worms were isolated in the F2 generation using a PCR screening approach, yielding healthy strain SD1880. To further validate wild-type function, the full unc-62:GFP as well as the isoform-specific unc-62(7b):GFP fusions were crossed into unc-62(e644) worms containing a loss of function mutation in exon 7b, yielding (unc-62:GFP;unc-62(e644)(strain SD1887) and unc-62(7b):GFP;unc-62(e644)(strain SD1879) that were isolated by PCR screening of the F2 generation (using allele-specific oligonucleotides targeting the e644 mutation). Both of these rescued strains have viable phenotypes.
Five tissue-specific RNAi strains from previous publications were used to test the tissue of action for unc-62: intestine-specific RNAi strain OLB11, hypodermal-specific RNAi strain NR222, muscle-specific RNAi strain NR350, neuronal- and vulval precursor cell-specific RNAi strain NK742, and uterine-specific RNAi strain NK640. To validate tissue-specific knockdown for intestine-specific and hypodermal-specific strains, these strains were crossed with the unc-62:GFP reporter, yielding SD1954 (Is[unc-62:GFP; unc-119(+)]; rde-1(ne219); Is[lin-26p:rde-1; lin-26p:nls:GFP; rol-6(su1006)]) and SD1955 (Is[unc-62:GFP; unc-119(+)]; rde-1(ne219); Is[elt-2p:rde-1; rol-6(su1006)]). These worms were then placed on control or GFP RNAi as L1 larvae, and imaged as L4 larvae. To assay the effect of unc-62 RNAi on daf-16 activity, we used a previously-described reporter that expresses cytoplasmic GFP from the sod-3 promoter [69].
Lifespan analysis
General C. elegans techniques as well as lifespan experiments were performed as previously described [10]. Unless otherwise noted, lifespan experiments were performed by placing adult worms (with visible eggs) on NGM plates containing 30 uM 5-fluoro-2′-deoxyuridine (FUDR). Deaths before day 7 of adulthood were censured. Log-rank tests were performed using OASIS [70].
RNAi experiments
RNAi clones used were obtained from the Ahringer RNAi library [71] and sequenced to verify proper insertions. To generate RNAi against exons 7a and 7b of unc-62, these exons (as well as ∼50 bp of flanking introns) were PCR amplified and inserted into the L4440 vector at the EcoRV restriction site, and transformed into HT115(de3) E. coli. RNAi knockdown experiments in adults (including lifespan experiments) were performed on NGM plates supplemented with 30 uM FUDR, 100 ug/mL Ampicillin, and 2 mM IPTG to induce dsRNA expression.
Supporting Information
Zdroje
1. HerndonLA, SchmeissnerPJ, DudaronekJM, BrownPA, ListnerKM, et al. (2002) Stochastic and genetic factors influence tissue-specific decline in ageing C. elegans. Nature 419: 808–814.
2. HuangC, XiongC, KornfeldK (2004) Measurements of age-related changes of physiological processes that predict lifespan of Caenorhabditis elegans. Proc Natl Acad Sci U S A 101: 8084–8089.
3. TankEM, RodgersKE, KenyonC (2011) Spontaneous age-related neurite branching in Caenorhabditis elegans. J Neurosci 31: 9279–9288.
4. TothML, MelentijevicI, ShahL, BhatiaA, LuK, et al. (2012) Neurite Sprouting and Synapse Deterioration in the Aging Caenorhabditis elegans Nervous System. J Neurosci 32: 8778–8790.
5. McGeeMD, WeberD, DayN, VitelliC, CrippenD, et al. (2011) Loss of intestinal nuclei and intestinal integrity in aging C. elegans. Aging Cell 10: 699–710.
6. ZahnJM, KimSK (2007) Systems biology of aging in four species. Curr Opin Biotechnol 18: 355–359.
7. GoldenTR, HubbardA, DandoC, HerrenMA, MelovS (2008) Age-related behaviors have distinct transcriptional profiles in Caenorhabditis elegans. Aging Cell 7: 850–865.
8. GoldenTR, MelovS (2004) Microarray analysis of gene expression with age in individual nematodes. Aging Cell 3: 111–124.
9. LundJ, TedescoP, DukeK, WangJ, KimSK, et al. (2002) Transcriptional profile of aging in C. elegans. Curr Biol 12: 1566–1573.
10. BudovskayaYV, WuK, SouthworthLK, JiangM, TedescoP, et al. (2008) An elt-3/elt-5/elt-6 GATA transcription circuit guides aging in C. elegans. Cell 134: 291–303.
11. Sanchez-BlancoA, KimSK (2011) Variable pathogenicity determines individual lifespan in Caenorhabditis elegans. PLoS Genet 7: e1002047 doi:10.1371/journal.pgen.1002047
12. KenyonCJ (2010) The genetics of ageing. Nature 464: 504–512.
13. Van AukenK, WeaverD, RobertsonB, SundaramM, SaldiT, et al. (2002) Roles of the Homothorax/Meis/Prep homolog UNC-62 and the Exd/Pbx homologs CEH-20 and CEH-40 in C. elegans embryogenesis. Development 129: 5255–5268.
14. YangL, SymM, KenyonC (2005) The roles of two C. elegans HOX co-factor orthologs in cell migration and vulva development. Development 132: 1413–1428.
15. JiangY, ShiH, LiuJ (2009) Two Hox cofactors, the Meis/Hth homolog UNC-62 and the Pbx/Exd homolog CEH-20, function together during C. elegans postembryonic mesodermal development. Dev Biol 334: 535–546.
16. PottsMB, WangDP, CameronS (2009) Trithorax, Hox, and TALE-class homeodomain proteins ensure cell survival through repression of the BH3-only gene egl-1. Dev Biol 329: 374–385.
17. CurranSP, RuvkunG (2007) Lifespan regulation by evolutionarily conserved genes essential for viability. PLoS Genet 3: e56 doi:10.1371/journal.pgen.0030056
18. NiuW, LuZJ, ZhongM, SarovM, MurrayJI, et al. (2011) Diverse transcription factor binding features revealed by genome-wide ChIP-seq in C. elegans. Genome Res 21: 245–254.
19. RozowskyJ, EuskirchenG, AuerbachRK, ZhangZD, GibsonT, et al. (2009) PeakSeq enables systematic scoring of ChIP-seq experiments relative to controls. Nat Biotechnol 27: 66–75.
20. GersteinMB, LuZJ, Van NostrandEL, ChengC, ArshinoffBI, et al. (2010) Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project. Science 330: 1775–1787.
21. TulletJM, HertweckM, AnJH, BakerJ, HwangJY, et al. (2008) Direct inhibition of the longevity-promoting factor SKN-1 by insulin-like signaling in C. elegans. Cell 132: 1025–1038.
22. RyooHD, MartyT, CasaresF, AffolterM, MannRS (1999) Regulation of Hox target genes by a DNA bound Homothorax/Hox/Extradenticle complex. Development 126: 5137–5148.
23. McGheeJD, SleumerMC, BilenkyM, WongK, McKaySJ, et al. (2007) The ELT-2 GATA-factor and the global regulation of transcription in the C. elegans intestine. Dev Biol 302: 627–645.
24. QadotaH, InoueM, HikitaT, KoppenM, HardinJD, et al. (2007) Establishment of a tissue-specific RNAi system in C. elegans. Gene 400: 166–173.
25. McGheeJD, FukushigeT, KrauseMW, MinnemaSE, GoszczynskiB, et al. (2009) ELT-2 is the predominant transcription factor controlling differentiation and function of the C. elegans intestine, from embryo to adult. Dev Biol 327: 551–565.
26. BeckAH, WengZ, WittenDM, ZhuS, FoleyJW, et al. (2010) 3′-end sequencing for expression quantification (3SEQ) from archival tumor samples. PLoS ONE 5: e8768 doi:10.1371/journal.pone.0008768
27. BreitlingR, ArmengaudP, AmtmannA, HerzykP (2004) Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett 573: 83–92.
28. SpiethJ, BlumenthalT (1985) The Caenorhabditis elegans vitellogenin gene family includes a gene encoding a distantly related protein. Mol Cell Biol 5: 2495–2501.
29. SpiethJ, MacMorrisM, BrovermanS, GreenspoonS, BlumenthalT (1988) Regulated expression of a vitellogenin fusion gene in transgenic nematodes. Dev Biol 130: 285–293.
30. MacMorrisM, BrovermanS, GreenspoonS, LeaK, MadejC, et al. (1992) Regulation of vitellogenin gene expression in transgenic Caenorhabditis elegans: short sequences required for activation of the vit-2 promoter. Mol Cell Biol 12: 1652–1662.
31. NoyesMB, ChristensenRG, WakabayashiA, StormoGD, BrodskyMH, et al. (2008) Analysis of homeodomain specificities allows the family-wide prediction of preferred recognition sites. Cell 133: 1277–1289.
32. FabianTJ, JohnsonTE (1995) Identification genes that are differentially expressed during aging in Caenorhabditis elegans. J Gerontol A Biol Sci Med Sci 50: B245–253.
33. PauliF, LiuY, KimYA, ChenPJ, KimSK (2006) Chromosomal clustering and GATA transcriptional regulation of intestine-expressed genes in C. elegans. Development 133: 287–295.
34. SpencerWC, ZellerG, WatsonJD, HenzSR, WatkinsKL, et al. (2011) A spatial and temporal map of C. elegans gene expression. Genome Res 21: 325–341.
35. NehrkeK (2003) A reduction in intestinal cell pHi due to loss of the Caenorhabditis elegans Na+/H+ exchanger NHX-2 increases life span. J Biol Chem 278: 44657–44666.
36. MurphyCT, LeeSJ, KenyonC (2007) Tissue entrainment by feedback regulation of insulin gene expression in the endoderm of Caenorhabditis elegans. Proc Natl Acad Sci U S A 104: 19046–19050.
37. HuangX, ChengHJ, Tessier-LavigneM, JinY (2002) MAX-1, a novel PH/MyTH4/FERM domain cytoplasmic protein implicated in netrin-mediated axon repulsion. Neuron 34: 563–576.
38. HamiltonB, DongY, ShindoM, LiuW, OdellI, et al. (2005) A systematic RNAi screen for longevity genes in C. elegans. Genes Dev 19: 1544–1555.
39. MurphyCT, McCarrollSA, BargmannCI, FraserA, KamathRS, et al. (2003) Genes that act downstream of DAF-16 to influence the lifespan of Caenorhabditis elegans. Nature 424: 277–283.
40. OhSW, MukhopadhyayA, DixitBL, RahaT, GreenMR, et al. (2006) Identification of direct DAF-16 targets controlling longevity, metabolism and diapause by chromatin immunoprecipitation. Nat Genet 38: 251–257.
41. Arantes-OliveiraN, ApfeldJ, DillinA, KenyonC (2002) Regulation of life-span by germ-line stem cells in Caenorhabditis elegans. Science 295: 502–505.
42. KimbleJ, SharrockWJ (1983) Tissue-specific synthesis of yolk proteins in Caenorhabditis elegans. Dev Biol 96: 189–196.
43. BantscheffM, RingelB, MadiA, SchnabelR, GlockerMO, et al. (2004) Differential proteome analysis and mass spectrometric characterization of germ line development-related proteins of Caenorhabditis elegans. Proteomics 4: 2283–2295.
44. GrantB, HirshD (1999) Receptor-mediated endocytosis in the Caenorhabditis elegans oocyte. Mol Biol Cell 10: 4311–4326.
45. MendenhallAR, WuD, ParkSK, CypserJR, TedescoPM, et al. (2011) Genetic dissection of late-life fertility in Caenorhabditis elegans. J Gerontol A Biol Sci Med Sci 66: 842–854.
46. NakamuraA, YasudaK, AdachiH, SakuraiY, IshiiN, et al. (1999) Vitellogenin-6 is a major carbonylated protein in aged nematode, Caenorhabditis elegans. Biochem Biophys Res Commun 264: 580–583.
47. DavidDC, OllikainenN, TrinidadJC, CaryMP, BurlingameAL, et al. (2010) Widespread protein aggregation as an inherent part of aging in C. elegans. PLoS Biol 8: e1000450 doi:10.1371/journal.pbio.1000450
48. Fredriksson A, Krogh EJ, Hernebring M, Pettersson E, Javadi A, et al.. (2012) Effects of aging and reproduction on protein quality control in soma and gametes of Drosophila melanogaster. Aging Cell.
49. YoungmanMJ, RogersZN, KimDH (2011) A decline in p38 MAPK signaling underlies immunosenescence in Caenorhabditis elegans. PLoS Genet 7: e1002082 doi:10.1371/journal.pgen.1002082
50. KirkwoodTB (1977) Evolution of ageing. Nature 270: 301–304.
51. LaiF, ChangJS, WuWS (2010) Identifying a Transcription Factor's Regulatory Targets from its Binding Targets. Gene Regul Syst Bio 4: 125–133.
52. MacIsaacKD, LoKA, GordonW, MotolaS, MazorT, et al. (2010) A quantitative model of transcriptional regulation reveals the influence of binding location on expression. PLoS Comput Biol 6: e1000773 doi:10.1371/journal.pcbi.1000773
53. HuZ, KillionPJ, IyerVR (2007) Genetic reconstruction of a functional transcriptional regulatory network. Nat Genet 39: 683–687.
54. SomelM, GuoS, FuN, YanZ, HuHY, et al. (2010) MicroRNA, mRNA, and protein expression link development and aging in human and macaque brain. Genome Res 20: 1207–1218.
55. BrackAS, ConboyMJ, RoyS, LeeM, KuoCJ, et al. (2007) Increased Wnt signaling during aging alters muscle stem cell fate and increases fibrosis. Science 317: 807–810.
56. KohK, RothmanJH (2001) ELT-5 and ELT-6 are required continuously to regulate epidermal seam cell differentiation and cell fusion in C. elegans. Development 128: 2867–2880.
57. HarrisTW, AntoshechkinI, BieriT, BlasiarD, ChanJ, et al. (2010) WormBase: a comprehensive resource for nematode research. Nucleic Acids Res 38: D463–467.
58. Thomas-ChollierM, HerrmannC, DefranceM, SandO, ThieffryD, et al. (2012) RSAT peak-motifs: motif analysis in full-size ChIP-seq datasets. Nucleic Acids Res 40: e31.
59. CrooksGE, HonG, ChandoniaJM, BrennerSE (2004) WebLogo: a sequence logo generator. Genome Res 14: 1188–1190.
60. FabianTJ, JohnsonTE (1994) Production of age-synchronous mass cultures of Caenorhabditis elegans. J Gerontol 49: B145–156.
61. LangmeadB, TrapnellC, PopM, SalzbergSL (2009) Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25.
62. Von StetinaSE, WatsonJD, FoxRM, OlszewskiKL, SpencerWC, et al. (2007) Cell-specific microarray profiling experiments reveal a comprehensive picture of gene expression in the C. elegans nervous system. Genome Biol 8: R135.
63. FoxRM, Von StetinaSE, BarlowSJ, ShafferC, OlszewskiKL, et al. (2005) A gene expression fingerprint of C. elegans embryonic motor neurons. BMC Genomics 6: 42.
64. ColosimoME, BrownA, MukhopadhyayS, GabelC, LanjuinAE, et al. (2004) Identification of thermosensory and olfactory neuron-specific genes via expression profiling of single neuron types. Curr Biol 14: 2245–2251.
65. ZhangY, MaC, DeloheryT, NasipakB, FoatBC, et al. (2002) Identification of genes expressed in C. elegans touch receptor neurons. Nature 418: 331–335.
66. BrockTJ, BrowseJ, WattsJL (2006) Genetic regulation of unsaturated fatty acid composition in C. elegans. PLoS Genet 2: e108 doi:10.1371/journal.pgen.0020108
67. KarlenY, McNairA, PerseguersS, MazzaC, MermodN (2007) Statistical significance of quantitative PCR. BMC Bioinformatics 8: 131.
68. SarovM, SchneiderS, PozniakovskiA, RoguevA, ErnstS, et al. (2006) A recombineering pipeline for functional genomics applied to Caenorhabditis elegans. Nat Methods 3: 839–844.
69. LibinaN, BermanJR, KenyonC (2003) Tissue-specific activities of C. elegans DAF-16 in the regulation of lifespan. Cell 115: 489–502.
70. YangJS, NamHJ, SeoM, HanSK, ChoiY, et al. (2011) OASIS: online application for the survival analysis of lifespan assays performed in aging research. PLoS ONE 6: e23525 doi:10.1371/journal.pone.0023525
71. KamathRS, FraserAG, DongY, PoulinG, DurbinR, et al. (2003) Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature 421: 231–237.
72. RoyPJ, StuartJM, LundJ, KimSK (2002) Chromosomal clustering of muscle-expressed genes in Caenorhabditis elegans. Nature 418: 975–979.
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2013 Číslo 2
Nejčtenější v tomto čísle
- Complex Inheritance of Melanoma and Pigmentation of Coat and Skin in Grey Horses
- Coordination of Chromatid Separation and Spindle Elongation by Antagonistic Activities of Mitotic and S-Phase CDKs
- Autophagy Induction Is a Tor- and Tp53-Independent Cell Survival Response in a Zebrafish Model of Disrupted Ribosome Biogenesis
- Assembly of the Auditory Circuitry by a Genetic Network in the Mouse Brainstem