Interplay of dFOXO and two ETS-family transcription factors determines lifespan in Drosophila melanogaster.

Forkhead box O (FoxO) transcription factors (TFs) are key drivers of complex transcriptional programmes that determine animal lifespan. FoxOs regulate a number of other TFs, but how these TFs in turn might mediate the anti-ageing programmes orchestrated by FoxOs in vivo is unclear. Here, we identify an E-twenty six (ETS)-family transcriptional repressor, Anterior open (Aop), as regulated by the single Drosophila melanogaster FoxO (dFOXO) in the adult gut. AOP, the functional orthologue of the human Etv6/Tel protein, binds numerous genomic sites also occupied by dFOXO and counteracts the activity of an ETS activator, Pointed (Pnt), to prevent the lifespan-shortening effects of co-activation of dFOXO and PNT. This detrimental synergistic effect of dFOXO and PNT appears to stem from a mis-regulation of lipid metabolism. At the same time, AOP activity in another fly organ, the fat body, has further beneficial roles, regulating genes in common with dfoxo, such as the secreted, non-sensory, odorant binding protein (Obp99b), and robustly extending lifespan. Our study reveals a complex interplay between evolutionarily conserved ETS factors and dFOXO, the functional significance of which may extend well beyond animal lifespan.

Introduction

Forkhead Box O (FoxO) transcription factors (TFs) play a key, evolutionarily conserved role in ageing. Drosophila melanogaster has a single FoxO orthologue (dfoxo) and increasing its activity in certain tissues is sufficient to extend fly lifespan [1]–[4] Furthermore, both dfoxo and the Caenorhabditis elegans othologue, daf-16, are strictly required for lifespan extension upon reduction in insulin/IGF-like signalling (IIS) [5], [6]. This evolutionary conservation appears to extend all the way to yeast on one side, where forkhead-like factors can extend lifespan [7], and to humans on the other, where certain variants of the FoxO3A locus are robustly correlated with longevity [8]–[12].

FoxOs control a plethora of traits at both organismal and cellular levels, including control of cell cycle, cell death, growth and metabolism. In all cases, FoxOs can be viewed as acting to preserve homeostasis [13]. Indeed, numerous processes are remodelled by activation of FoxOs, through regulation of a large number of direct and indirect targets, all acting in concert to preserve homeostasis in old age and extend animal lifespan [14]–[19].

Several studies have examined the targets of FoxOs. A striking finding of these studies is that FoxOs control a range of other cellular regulators. These include secreted endocrine factors, components of intracellular signalling pathways and several TFs [14], [16]–[20]. Transcriptional feedback within the signalling pathway plays a role [21], but in most cases the functions of these other regulators remain unknown, both in isolated cells and, more importantly, in vivo.

The putative roles of TFs regulated by FoxOs are particularly intriguing. Numerous studies have shown that FoxOs interact with a number of unrelated TFs, in a number of ways, with important consequences for the output of both interacting partners [13], [22]. These TFs include Myc, p53, Smads, ß-catenin, and numerous nuclear hormone receptors [22]–[27]. Hence, there is a potential for the TFs regulated by FoxO to profoundly alter FoxO's functional output through interactions with FoxO itself. However, it remains unclear what the role of these interactions is in the whole animal, in vivo and, specifically, what is their role in lifespan?

In this study we set out to elucidate the role played in lifespan by a TF directly regulated by dFOXO. We identify an E-twenty six (ETS) - family transcriptional repressor, Anterior open (Aop), as regulated by dFOXO in the adult Drosophila gut. Aop is the functional orthologue of the human Etv6 gene and, in Drosophila, it is known to counteract the activity of an ETS activator, Pointed (Pnt). We show that Aop acts to prevent the detrimental effects of co-activation of dFOXO and PNT in adult Drosophila gut, and we present evidence that this interaction is mediated by binding to the same genomic locations as dFOXO. AOP activation on its own in the adult fat body can also robustly extend lifespan. Our study reveals a complex interplay between evolutionarily conserved ETS-family TFs and dFOXO in longevity. The significance of this interplay may extend to other physiological processes.

Results

dFOXO regulates distinct genes but similar functions in the adult gut and fat body

dFOXO, like its mammalian orthologues, controls gene expression in a tissue-specific manner [19], [28]–[30]. Hence, to investigate the functional interplay between dFOXO and one of its target TFs, we turned our attention to a tissue-specific, adult-inducible, lifespan-relevant system. Over-expression of dfoxo using the RU486-inducible, S Geneswitch driver [31], robustly extends lifespan [1], [4], [32]–[34]. S restricts dfoxo induction to two specific adult fly organs: the midgut and abdominal fat body (subsequently referred to as gut and fat body; ) [31], the latter functionally equivalent to mammalian white adipose tissue and liver. Both have an evolutionarily conserved role in aging [35], [36], and it is currently unclear whether activation of dfoxo in either organ alone is sufficient to extend lifespan. For these reasons, we chose to identify the TFs regulated by dFOXO in both of these organs.

We micro-dissected mid-guts or carcass-associated thoracic/abdominal fat body of S females (+/− RU486) and determined their mRNA profiles using Affymetrix gene expression arrays (ArrayExpress accession number: E-MTAB-1020). In each case, we controlled for the changes associated with induction of the driver alone (S +/− RU486). 447 genes were differentially expressed in the gut (p value cut-off of 0.00285 corresponding to FDR of 5%, ). We detected fewer significant changes in the fat body, 87 differentially regulated genes (p value cut-off 0.0022, FDR 20%, ), most-likely due to the difficulty of dissecting this loosely-associated tissue. The full list of genes regulated by dFOXO, as well as all other lists mentioned in the paper, are given in . The list included some well-known targets of dFOXO, such as initiation factor 4E binding protein (4ebp) [37], [38] and the Drosophila insulin receptor (dInR) [37], both activated in the gut, and the insulin-regulated kinase Akt [17], [39], induced in both the gut and fat body.

Figure 1

dFOXO targets in the adult gut and fat body.

A Proportional Venn diagram showing the sets of genes that were differentially regulated by dfoxo induction in the gut or the fat body. The number of genes that were bound by GFP-dFOXO within each differentially expressed-gene set are given in black. p values for significant set overlaps are indicated. B Biological process GO categories differentially regulated (p<10−10) in the fat body or gut upon induction of dfoxo as determined by Catmap analysis. Any redundant categories (overlap by more than 75%) were removed, retaining the most specific category. The full list is given in . The intensity of red shows the log10-transformed p-value associated with differential regulation for each category.

The overlap between genes regulated in the same direction in gut and fat body was significant (p = 0.002 for up-regulated, p<10−15 for down-regulated genes), but by no means complete ( ), indicating that dFOXO regulates both related and unrelated sets of genes in the two tissues. This did not appear to be caused by weaker signal obtained from the fat body, since changes unique to the fat body were also detected (e.g. Obp99b; see later). To further explore this, we examined the functions regulated by dFOXO in the two tissues using Catmap, an approach that can detect co-ordinated, subtle changes over many genes, rather than depending on an arbitrarily chosen, differential expression p-value cut-off [40]. Although some differences occurred, the biological process Gene Ontology (GO) categories differentially expressed in the two tissues were similar ( ), indicating that dFOXO regulates similar functions in the gut and in the fat body.

Notably, dFOXO strongly repressed respiratory electron transport chain components in both the gut and fat body (p = 4×10−25 and p = 1×10−21, respectively) and, in particular, components of complex I, which transfers electrons from NADH to ubiquinone (p = 2×10−19 and p = 2×10−17). Indeed, part of the effect of dFOXO on lifespan could be mediated by its repression of the components of complex I, because reducing the electron flow through this complex, by bypassing it, can extend fly lifespan [41], [42].

5 TFs are directly regulated by dFOXO in the adult gut

We identified a total of 16 TFs regulated by induction of dfoxo in either the gut or the fat body, including p53 and the nuclear hormone receptor HR96 (see for the full list). To further narrow down the set of interesting candidates, we isolated the TFs encoded by genes directly bound and regulated by dFOXO, thus identifying the immediate second tier of regulators. We determined the genomic regions bound by dFOXO using a GFP-dFOXO fusion protein that is functional in lifespan-extension [33]. We prepared chromatin from RU486-fed S female flies and pulled down the DNA associated with the fusion protein, and hence restricted to gut and fat body, using an anti-GFP antibody. As a control, we performed chromatin immunoprecipitation (ChIP) using the same antibody on chromatin prepared from females over-expressing dfoxo alone (S + RU486).

To confirm that we were detecting tissue-restricted binding, we compared the GFP-dFOXO bound sites with those that we previously identified in whole adults or S2 cells [17]. The dInR locus is transcribed from three promoters under tight spatio-temporal control [43]. In whole flies, dFOXO is detected as bound in the coding region of the gene and is absent from the P1 or P3 regions [17] ( ). The functional significance of this binding in the 3′ region in Drosophila is unclear but the mammalian FoxO proteins are able to act at great distances [20]. On the other hand, in serum-starved S2 cells dFOXO is bound to the P1 promoter but not the coding or P3 regions [17] ( ). In contrast, GFP-dFOXO expressed in the adult gut and fat body was bound to the P3 promoter of dInR ( ), revealing a different pattern of binding in these tissues.

Figure 2

dFOXO binding sites in adult gut and fat body.

A Enrichment of the P3 promoter sequences of the dInR locus in the chromatin samples prepared from RU486-fed S, RU486-fed S (mock for the anti-GFP IP), wild-type 7-day old females or 2-h serum-starved S2 cells, after IP with either anti-GFP antibody or anti-dFOXO antibody, as indicated. Enrichment is expressed relative to U6, as means ± SEM of three biological replicates of chromatin, except for S2 cells where three IPs were performed from the same chromatin sample. ANOVA on log-transformed data detected significant differences (p<10−3), and enrichment in S, after IP with anti-GFP antibody, was greater than all others (t-test, p<10−3). B ChIP-chip traces, showing the enrichment (log2-transformed) of the GFP-dFOXO-immunoprecipitated DNA over total chromatin, are averages of three biological repeats after subtraction of the mock and are shown over a region of chromosome 3R. Red denotes the enrichment associated with peak regions. C Summary of the regulatory relationships between dFOXO and the five TFs it directly induces in the adult gut. Arrows indicate transcriptional activation.

We hybridised three biological repeats of the experimental and control ChIP samples to tiling arrays. We identified ∼1400 genomic regions bound by GFP-dFOXO in the gut and/or the fat body (ArrayExpress accession number: E-MTAB-1021; for examples of peaks identified see ; for a list of all bound locations see ), with highly reproducible ChIP-chip signal across the three biological replicates (). The regions occupied by dFOXO in the gut and fat body were different from those previously identified in whole flies (), further confirming that we were detecting tissue-restricted binding. They were predominantly located in the 5′ end of genes (), indicating promoter-proximal binding, in contrast to enhancer binding observed with mammalian FoxO3 [20]. The regions bound by dFOXO were enriched for forkhead-like binding motifs (), confirming conservation of in vivo binding-sequence preference.

Finally, we identified the genes that are likely to be directly regulated by dFOXO based on their transcriptional responsiveness to induction of dfoxo and proximity (<1 kb) to a GFP-dFOXO bound site ( ). Genes up-regulated in the gut were specifically enriched for the GFP-dFOXO bound genes (p = 5×10−9). This is consistent with the predominant function of dFOXO as a transcriptional activator, conserved in its worm and mammalian orthologues [16], [19], [20]. We also observed GFP-dFOXO binding in the vicinity of the genes regulated by dfoxo in the fat body ( ), but this overlap was not significant. Note that numerous GFP-dFOXO-bound sites could not be associated with specific transcriptional events. This has been observed previously for FoxO factors in flies and other organisms [16], [17], [19], [20], [28], and it is currently unclear whether this is due to technical limitations in associating expression changes to binding events, whether FoxOs detected on these sites are poised for activation under a different set of conditions, or whether some of these sites are not functional.

The set of direct dFOXO targets included five sequence-specific TFs: clockwork orange (cwo), mothers against dpp (Mad), bunched (bun), anterior open (Aop) and CG12054 (summarised in ). All of these were activated by dFOXO in the gut (, for Aop see also ). dFOXO also bound its own locus, and inducing the transgenic and intron-less dfoxo in S females resulted in an increase in unspliced dfoxo (), indicating that dFOXO self-activates ( ).

Figure 3

dFOXO regulates expression of Aop in the adult gut.

A Schematic of the Aop locus with black boxes representing exons, red boxes – regions detected as bound by GFP-dFOXO in the ChIP-chip experiment on induced S females (dFOXO gut/fat body), yellow box - region detected as bound by dFOXO in wild-type females (dFOXO whole fly, data obtained from reference [17]), and black bars – position of amplicons used for qPCR in B. B The enrichment of 5′ or 3′ end of the Aop locus, relative to U6, after anti-GFP IP of chromatin from RU486-induced S females (mock), anti-GFP IP of chromatin from RU486-induced S females (gut/fat body) or anti-dFOXO IP of wild-type female chromatin (whole fly). Means ± SEM of three biological repeats are shown, with enrichment in the mock control set to one. ANOVA on log-transformed data detected significant differences (p = 0.03 per region) and the enrichment of the 5′ region was different in gut/fat body from the mock (one-tailed t-test, p = 6×10−3), while the 3′ region was enriched in the whole fly (one-tailed t-test, p = 5×10−3). C Aop mRNA was quantified relative to Act by qPCR in guts or fat bodies of S, or TIGS>dfoxo flies induced or not with RU486. Boxplots show log-10 derived relative expression with - RU486 values set to zero. Data for S females were analysed with a mixed-effects linear model with dissection batch as a random effect. The effects of RU486, tissue and their interaction was significant (p<0.05) and RU486 caused significant up-regulation of Aop in the gut (one-tailed t-test, n = 3–4, p = 2×10−3) but not the fat body (one-tailed t-test, n = 4, p>0.05). Significant changes were observed in TIGS>dfoxo guts (t-test, n>3, p = 0.02).

dFOXO induces an ETS transcriptional repressor in the adult gut

Among the TFs discovered as directly activated by dFOXO in the gut, we were especially interested in Aop (a.k.a. Yan). This ETS-family transcriptional repressor has a clear human orthologue in the Etv6 (a.k.a. Tel) gene [44]. Aop and Etv6 display clear conservation of function, with Aop involved in tracheal sprouting in flies while Etv6 is involved in the equivalent process of endothelial sprouting in mammals [45]. Interestingly, mammalian FoxOs and Etv6 act in similar physiological processes: both are tumour suppressors required for maintenance of adult haematopoietic stem cells [30], [46]; indicating that the previously-uncharacterised, functional interplay between the two factors may be evolutionarily conserved. Importantly, Drosophila presents a unique opportunity to examine the physiological functions of Aop/Etv6, because there is no known orthologue in yeast or worm.

ChIP-chip revealed that GFP-dFOXO bound in the promoter-proximal, 5′ end of Aop, specifically in the 1st exon and 1st intron ( ). The binding to the 1st exon in the gut/fat body was confirmed with qPCR ( ). Similarly to dInR, we previously observed dFOXO bound to the coding region (3′ end) of this gene in whole flies (see [17] and ). The functional significance of this 3′-end binding is unclear [17].

Induction of dfoxo in the gut and the fat body resulted in significant induction of Aop transcript only in the gut (p = 0.002, ). To confirm the statistical significance of this tissue-restricted effect of dfoxo, we analysed the Aop expression data with a mixed-effects linear model and found a significant difference in the way the two tissues respond to RU486 (p<0.05). To confirm that these effects were cell autonomous, we induced dfoxo solely in the gut using an RU486-inducible, gut-specific driver, TIGS [31]. Feeding RU486 to TIGS>dfoxo females also resulted in induction of Aop transcript in the gut ( ). Hence, dFOXO induces Aop transcription in the gut, most likely through direct binding to the Aop promoter in gut cells.

Aop prevents the detrimental co-activation of dfoxo and Pointed

To understand the relationship between dfoxo and its target TF, Aop, we next examined the physiological role played by Aop in the context of tissue-restricted dfoxo induction and lifespan. Since Aop is an essential gene [47], [48], we chose to knock it down in the adult gut and fat body using a short-hairpin RNAi construct with no predicted off-targets [49]. Driving this construct with a ubiquitous, constitutive driver (daughtelessGAL4) resulted in the expected embryonic lethality. Inducing this construct in the adult gut and fat body, using the S driver, reduced the levels of Aop mRNA by ∼70% (p = 0.04, ) but had no major effect on lifespan ( ), revealing that, in a wild-type fly, Aop in these tissues is not limiting for survival.

Figure 4

Aop prevents the detrimental effects of dfoxo and Pnt co-activation.

A Survival of female flies expressing dfoxo, Aop or both under the control of S in the presence or absence or RU486. Log-rank test revealed significant effect of RU486 for S (p<10−4; total dead/censored: − RU486 145/0, + RU486 139/6; median/maximum lifespan: − RU486: 77/87, + RU486 80/90) and S (p<10−4; total dead/censored: − RU486 145/2, + RU486 155/2; median/maximum lifespan: − RU486: 74/90, + RU486 65/79) but not S (p = 0.5; total dead/censored: − RU486 142/0, + RU486 147/4; median/maximum lifespan: − RU486: 72/86, + RU486 74/87). CHP analysis showed that the response to RU486 in S was significantly different from the response in S (p<5×10−15) and S (p<10−7) females. B Survival of female flies expressing dfoxo, PntP1 and Aop under the control of S in the presence or absence or RU486. Log-rank test revealed significant effect of RU486 for S (p<10−4; total dead/censored: − RU486 138/6, + RU486 145/2; median/maximum lifespan: − RU486: 77/90, + RU486 63/83) but not S (p = 0.06; total dead/censored: − RU486 147/0, + RU486 143/2; median/maximum lifespan: − RU486: 84/92, + RU486 79/90) or S (p = 0.1; total dead/censored: − RU486 147/3, + RU486 142/3; median/maximum lifespan: − RU486: 77/90, + RU486 77/83). CHP analysis showed that the response to RU486 in S was significantly different from the response in S (p<10−15) but not in S (p = 0.5) females. C Survival of TIGS>dfoxo female flies in the presence or absence of RU486. Log-rank test detected no significant differences (p = 0.3; total dead/censored: − RU486 152/1, + RU486 149/2; median/maximum lifespan: − RU486: 77/87, + RU486 77/86). D Survival of female flies expressing dfoxo, Pnt under the control of TIGS in the presence or absence or RU486. Log-rank test revealed significant effect of RU486 for TIGS>Pnt (p<10−4; total dead/censored: − RU486 140/2, + RU486 136/7; median/maximum lifespan: − RU486: 74/86, + RU486 35/50), TIGS>Pnt (p<10−4; total dead/censored: − RU486 146/3, + RU486 148/1; median/maximum lifespan: − RU486: 74/85, + RU486 32/45) and TIGS>Pnt (p = 0.1; total dead/censored: − RU486 143/1, + RU486 155/1; median/maximum lifespan: − RU486: 72/85, + RU486 35/55). CHP analysis showed that the response to RU486 in TIGS>Pnt was significantly different from the response in TIGS>Pnt (p = 2×10−4) but not in TIGS>Pnt (p = 0.8) females.

To test if Aop is required for the dfoxo-induced longevity, we simultaneously expressed dfoxo and knocked-down Aop. Surprisingly, while dfoxo alone extended lifespan, this combined treatment was detrimental to the fly ( ), indicating that Aop is required to prevent some toxic effect of dfoxo activation [39]. To determine the statistical significance of this synergistic effect, we used Cox Proportional Hazards (CPH) analysis [50]. This type of survival analysis can determine the significance of individual covariates as well as their interaction. We used “genotype” and “RU486” as individual covariates and found that the response to RU486 was significantly different in S genotype from the response in either S or S genotypes (interaction of “genotype” and “RU486”, p = 5×10−15 and p = 1×10−7, respectively). This confirmed the synthetic interaction between the loss of Aop and induction of dfoxo. Importantly, expressing a control RNAi construct targeting GFP with S driver did not prevent dfoxo from extending lifespan ().

We next examined the mechanism underlying this toxic effect of inducing dfoxo in the absence of Aop. Aop's activity is known to counteract that of Pointed (Pnt), an ETS-family transcriptional activator [51], [52]. These two TFs regulate the same genes through binding to the same regulatory elements but with opposing outcomes [51]–[53]. Hence, the synthetic toxicity of Aop loss- and dfoxo gain-of-function indicated that co-activation of Pnt and dfoxo could be highly detrimental. To test this, we induced the constitutively active form of Pnt (Pnt) [51] in the adult gut and fat body. This resulted in a marginal negative effect on lifespan (p = 0.06) and, similarly to the phenotype of Aop loss and dfoxo gain, the co-induction of dfoxo and Pnt was highly detrimental ( ).

It was possible that dFOXO induced the transcription of Aop in order to avert co-activation of PNT and the resulting detrimental synergistic effect. Indeed, the detrimental effect of combined dFOXO and PNT activity was completely rescued by additional activation of AOP, achieved by expression of a constitutively active form of Aop [48], called Aop, in addition to dfoxo and Pnt ( ). The significance of this effect was confirmed by CPH analysis, which showed that the response to RU486 in S genotype was significantly different from that in S genotype (p<10−15), but not from the response observed in S (p = 0.5). Importantly, induction of Aop did not interfere with the induction of dfoxo or Pnt in the S females compared to the S females ().

dFOXO regulates Aop transcription specifically in the gut ( ). For this reason we next examined whether the genetic interactions we observed between dfoxo and Aop/Pnt could be recapitulated with the gut-specific TIGS driver. This driver has been successfully used by others to determine the contribution of the gut to S-triggered phenotypes [54]. The gut is a highly regionalised organ with several cell types [55], [56] and we confirmed that the two drivers activate transgene expression in similar regions (), both driving expression in the enterocytes (). Feeding RU486 to TIGS>dfoxo females had no significant effect on lifespan ( ), even though dfoxo was induced in the gut () and activated the expression of Aop ( ). With the caveat that the levels of expression achieved with TIGS and S in the gut are not identical (, ) and the two drivers may not drive in completely overlapping subsets of gut cells ( and B), our data indicate that activation of dFOXO in the gut alone is insufficient for lifespan extension.

On the other hand, the induction of Pnt solely in the adult gut was highly detrimental for the fly ( ), while RU486 had no effect on the lifespan of the TIGS-alone or UAS-Pnt-alone controls ( and D). The interaction between dfoxo and Pnt revealed with the TIGS driver was similar to that observed with S: induction of dfoxo exacerbated the toxicity of Pnt and this, in turn, could be remedied by further induction of Aop ( ). While the magnitude of this synergistic effect was smaller with TIGS than S, CHP analysis confirmed that the response to RU486 in TIGS>Pnt females was significantly different from that in TIGS>Pnt (p<10−4), but not in TIGS>Pnt females (p = 0.8). Hence, the dFOXO-PNT-AOP interaction can be mapped specifically to the gut. Note, however, that we cannot exclude the possibility of a similar interaction also occurring in the fat body.

Overall, our data are consistent with dFOXO up-regulating Aop transcription in the gut in order to counteract the activity of PNT. In this way, dFOXO prevents the detrimental effect that would result from dFOXO being active at the same time as PNT in that organ.

dFOXO and AOP share binding locations

Intuitively, the genetic interactions observed between dfoxo and Aop/Pnt were not consistent with a simple, linear cascade of dFOXO acting on AOP/PNT to influence lifespan ( ). To formally examine this, we used Boolean network modelling. We created network configurations with dFOXO, AOP, PNT and “lifespan” as nodes that can take values of 0 or 1 (inactive/short life or active/long life) and described the relationships between them using Boolean logic operators. We then perturbed the network by fixing either dFOXO or PNT or both as “active” and examined the probability that lifespan will take on the value of 1 (long life) after 1000 state transition using Markov chain simulations. This type of modelling formalised what was intuitively evident: the synergistic negative effect of dfoxo and Pnt co-activation could be recapitulated by a network circuit that includes a feed-forward loop between dfoxo and Pnt with a negative effect on lifespan (described by a NAND operator), but not by a linear cascade ( ).

Figure 5

dFOXO and AOPACT share binding sites across the genome.

A Boolean network models of dFOXO-AOP-PNT-lifespan interaction. The linear model is indicated in grey and the feedforward model in red. The logic gate symbols for AND, NOT and NAND (NOT AND) are used, and in the case shown NAND means that only the combined effect of dFOXO and PNT is detrimental to lifespan. For each model the probability of a positive outcome for lifespan (lifespan = 1; y axis) after 1000 state transition was determined using Markov chain simulations after fixing either dFOXO or PNT or both as active ( = 1; x axis). Only the feedforward model can capture the synergistic negative effects of the two TFs. B ChIP-chip traces, showing the enrichment (log2-transformed) of the FLAG-AOPACT-immunoprecipitated DNA over total chromatin, are averages of three biological repeats after subtraction of the mock control and are shown over a region of chromosome 3R. Red denotes the enrichment associated with peak regions. C Proportional Venn diagram showing the overlap between the genes in the vicinity (<1 kb) of GFP-dFOXO or FLAG-AOPACT binding in the gut/fat body (Hypergeometric test, p<10−5). D Proportional Venn diagram showing the overlap in GFP-dFOXO-bound and FLAG-AOPACT-bound genomic sites. Note that the average overlap is shown, since one peak in one set can overlap multiple peaks in the other due to differing peak lengths: 819 GFP-dFOXO-bound and 810 FLAG-AOPACT-bound sites are in the overlap. Bootstrap analysis revealed the overlap as significant (p<10−3).

Feed-forward loops involve co-regulation of the same targets, and are often accompanied by extensive overlaps in genomic sites bound by the TFs [57]. To gather further evidence for the existence of this feed-forward loop, we determined the relationship between the sites bound by the AOP/PNT couple and dFOXO. We performed ChIP-chip on flies expressing FLAG-AOPACT, since both PNT and AOP are known to bind the same sequences, competing for the same sites [51]–[53].

To facilitate the comparisons with dFOXO, we used the S driver and the RU486 inducer, with the control IP performed with an anti-FLAG antibody on chromatin from flies expressing untagged AOPACT (ArrayExpress accession number: E-MTAB-1306). We discovered ∼4000 genomic regions bound by FLAG-AOPACT neighbouring some 3000 genes, including the dfoxo locus, with good correlation between biological repeats (). An example of peaks is given in . The binding locations were associated with genes and tended to occur in the 5′ region (), consistent with a role in promoter-proximal regulation of transcription. Similarly to the recently published AOP binding locations observed in larvae [58], we found that AOP tended to bind long stretches of DNA, longer than dFOXO (). At the same time, over 80% of the regions bound in adult gut/fat body were distinct from those bound in larvae ().

We found a significant and substantial overlap in the genes bound by GFP-dFOXO and those bound by FLAG-AOPACT (p<10−5, ), indicating that AOP/PNT and dFOXO may regulate numerous genes in common. Furthermore, dFOXO did not just bind in the vicinity of the same genes but actually to the same regions of the DNA: 60% of GFP-dFOXO-bound regions directly overlapped regions of FLAG-AOPACT binding ( ). This striking and significant (p<10−3) overlap in sites bound corroborates the existence of a feed-forward loop between AOP/PNT and dFOXO. Furthermore, the genetic interactions we observed between dfoxo and Pnt/Aop are likely to be mediated by functional interactions on the shared regulatory regions.

Pnt and dfoxo synergistically affect lipid metabolism

The strength of the negative synergistic effect observed between dfoxo and Pnt prompted us to seek out its physiological basis. To initiate this investigation, we looked at the GO categories that are over-represented in genes bound by FLAG-AOPACT. The most over-represented functional category was “lipid particle” (p = 3×10−17, ), a GO category that includes lipid droplets, sites of cellular lipid storage. This was also the most over-represented category in the set of genes bound by both FLAG-AOPACT and GFP-dFOXO (p = 3×10−5, ), indicating that PNT/AOP and dFOXO together regulate lipid metabolism.

Figure 6

dfoxo and Pnt act synergistically on lipid metabolism.

A Top three most significantly enriched GO categories within genes bound by FLAG-AOPACT alone or FLAG-AOPACT and dFOXO, as determined by EASE analysis. The intensity of red indicates log-10 derived p value associated with the enrichment. The complete GO analysis is given in . B CG6296 mRNA was quantified relative to Act by qPCR in the females of the indicated genotypes, induced or not with RU486. Boxplots show log-10 derived relative expression with − RU486 values set to zero. Data were analysed with a linear model and the effects of RU486, genotype and their interaction were significant (n = 4, p<10−4), with S + RU486 condition being different to all others (t-test , p<10−4). C The levels of TAG were quantified under the same conditions as in B. Data were analysed with a linear model and the effects of RU486 (p = 2×10−4), genotype (p<10−4) and their interaction (p = 0.01) were significant (n = 8, where 10 measurements were made and the lowest and highest measurement removed from each group). The S + RU486 condition was different from all others (t-test , p<0.05). The genotypes are denoted with the same colours in B and C, with the legend given in B.

Activation of dFOXO in the aging gut represses two gastric lipases, encoded by lipA (a.k.a. magro) and CG6295 genes [39]. These two lipases are thought to facilitate assimilation of ingested lipids, and their repression by dFOXO results in a decrease in lipid stores and a shortened lifespan [39]. We observed the detrimental effects of dFOXO activation in the gut to be dependent on simultaneous activation of PNT ( ), prompting us to examine whether the expression of either lipA or CG6295 was synergistically regulated by dfoxo and Pnt. We used the S driver, where the synergistic effect on lifespan is pronounced, and found that both dfoxo and Pnt resulted in repression of CG6295, 3-fold and 8-fold, respectively ( ). When the two TF were co-expressed, the reduction in CG6295 was even greater, reaching 80-fold reduction after only 5 days of transgene induction ( ). Indeed, a linear model revealed a significant difference in the way the three genotypes responded to RU486 (p<10−4), and the amount of CG6295 transcript was significantly different between S>PntP1 + RU486 or S + RU486 and S>Pnt + RU486 conditions (p<10−4). Hence dfoxo and Pnt synergistically repress the CG6295 lipase. The effect of the two TFs on lipA was similar but smaller, and there was no significant difference between S + RU486 and S>Pnt + RU486 conditions (p>0.05, ). Note that neither AOPACT nor dFOXO bind in the vicinity of CG6295 and hence the effects on the expression of this gene are likely to involve an intermediate factor.

To examine whether the synergistic effect on the expression of the CG6295 lipase had a significant physiological consequence, we determined the levels of triacylglycerol (TAG) stores in flies expressing dfoxo, Pnt or both, under the control of the S driver. We found that the levels of TAG paralleled the levels of CG6295 transcript ( ) and, using a linear model, we confirmed that the effects of RU486 on the TAG levels in the 3 genotypes were significantly different (p = 0.01), with the TAG being significantly more depleted in S flies than in the other two genotypes after only 5 days of RU486 feeding (p<0.02, ).

The effect on CG6295 expression and TAG levels could be symptomatic of a more general loss of gut integrity. The function of the gut as a barrier is a good surrogate for overall gut integrity and is important for survival [59]. The effectiveness of this barrier can be assessed by feeding flies a food containing a dye (fluoresceine) and scoring the number of flies in which the gut is unable to exclude this dye from the whole body (“smurf” phenotype) [59]. We found that driving dfoxo and Pnt for 3 weeks using S did not increase the proportion of smurfs in the population (), while the effect on lipid metabolism was evident after only 5 days ( ). This confirmed that the effect on lipid metabolism was not due to a general loss of gut integrity. In addition, we also observed no significant changes in feeding behaviour of S>Pnt females after RU486 treatment, as assessed with the proboscis-extension assay ().

Based on these data, we propose that the detrimental synergy between dfoxo and Pnt arises from a mis-regulation of lipid metabolism resulting in a profound drop in TAG stores. Indeed, S flies were starvation sensitive after 5 days of RU486 feeding, and this sensitivity could be reversed by co-induction of Aop ().

Activation of AOP alone can extend lifespan

We next focused on other possible roles of Aop in Drosophila lifespan. Expression of the activated form of Aop in the flies co-expressing dfoxo and Pnt with S driver had a substantially beneficial effect on lifespan ( ). This was due, in part, to its role in preventing dFOXO and PNT co-activation, but we hypothesised that activation of AOP might additionally increase the lifespan of wild type flies. We induced expression of the activated form of AOP (Aop) using the S driver and the RU486 inducer in an otherwise wild-type adult female. This resulted in significant lifespan-extension (p = 2×10−10), increasing the median by 14% and maximal lifespan by 11% ( ), while RU486 feeding had no effect on the lifespans of the driver- or transgene-alone controls ( and B). This effect was robust and observed in six independent experiments, performed in the course of three years, with average median lifespan extension of 12% ( ). Furthermore, independently generated, FLAG-tagged Aop also extended lifespan ().

Figure 7

Aop extends lifespan.

A Survival of S female flies in the presence or absence or RU486. Log-rank test detected significant differences (p = 2×10−10; total dead/censored: − RU486 145/1, + RU486 129/3; median/maximum lifespan: − RU486: 74/87, + RU486 84/97). B Median lifespan extension achieved by RU486 feeding in S (6 experiments) or TIGS>Aop (3 experiments) females. Log-rank test detected significant extension (p<0.05) of lifespan in 6 out of 6 S and in 1 out of 3 TIGS>Aop trials. In one TIGS>Aop trial lifespan was shortened. C Haemolymph glucose and trehalose in S females fed or not RU486, where RU486 had a significant effect in each case (t-test, p = 0.01 and 0.02 respectively; n = 8 where 10 measurements were made and the highest and lowest measurement removed from each group). D Survival of S or dfoxoΔ/Δ S female flies in the presence or absence or RU486. Log-rank test detected significant differences for both: S (p<10−4; total dead/censored: − RU486 139/3, + RU486 136/7; median/maximum lifespan: − RU486: 77/87, + RU486 86/94) and dfoxoΔ/Δ S (p<10−4; total dead/censored: − RU486 138/0, + RU486 136/4; median/maximum lifespan: − RU486: 56/81, + RU486 72/90). CHP analysis revealed significant effects of RU486 (p<10−15) and dfoxo (p<10−15) but no significant difference in the response to RU486 between the two lines (p = 0.2).

To examine whether this lifespan-extending effect could be localised to the same tissues where the interaction between dfoxo and Pnt/Aop occurs, we drove Aop with the gut specific TIGS driver. Induction of Aop using the TIGS driver did not extend lifespan ( ), despite transgene induction in the gut with both TIGS and S drivers (). Although the induction levels were lower with TIGS (), they must still be physiologically relevant because, even with this driver, Aop could partially remedy the toxicity of Pnt and dfoxo co-activation ( ).

Our data are consistent with two roles of Aop in lifespan modulation: (1) Aop can counteract the negative effect of dfoxo and Pnt co-activation in the gut, (2) its additional activation in the fat body extends wild-type lifespan. The lifespan extension by Aop in flies co-expressing dfoxo and Pnt under S control most likely combines the two beneficial effects of Aop activation, since the effect appears greater in that context than in the otherwise wild-type female (median increased by 22% between S + RU486 and S + RU486, , versus 12% average extension in wild-type, ).

We examined other phenotypes triggered in S females by RU486 feeding, focusing on phenotypes often associated with lifespan extension. RU486 feeding did not cause any significant changes in starvation, hydrogen peroxide or DDT resistance; any changes to whole body trehalose, glycogen or lipid content; feeding or fecundity (). However, we did find that induction of Aop in the adult gut and fat body resulted in increased circulating sugars ( ), a diabetic phenotype potentially indicative of slightly reduced IIS. Pnt is known to promote IIS in larvae and its gain-of-function stimulates clearance of circulating sugars [60], consistent with the converse phenotype we observed in S adults.

Since FoxOs are inhibited by IIS [61], Aop could act to reduce IIS and activate dFOXO, thereby increasing the lifespan of otherwise wild-type females. However, Aop extended lifespan in complete absence of dfoxo ( ), and CPH analysis found no evidence for a difference in the response to RU486 in dfoxoΔ/Δ S versus S females (p = 0.2). Hence, Aop extends wild-type lifespan independently of dfoxo.

dFOXO and AOP share targets in the fat body

We further examined the function of Aop in the fat body, where its activity appears important for lifespan but where it is not regulated by dfoxo. We found that AOP transcript and protein levels were strongly induced in the fat body when S females were fed RU486 (), but in only a portion (∼26±6%) of the fat body cells, suggesting heterogeneity of the tissue ( , for quantification see ). This mosaic expression was not due to the driver, because S drove GFP expression in all the cells of the tissue ().

Figure 8

Aop and dfoxo both regulate a humoral factor, Obp99b, in the fat body.

A AOP was visualised by immunofuorescence in fat bodies of S flies induced or not with RU486. In the merged image, AOP is shown in green, DAPI-stained nuclei in blue. AOP accumulates in nuclei of 26±6% of fat body cells. B Genes regulated by RU486 induction in fat bodies of S flies were determined by microarray analysis and compared to changes in the fat body upon induction of S. Mean log2 fold change caused by induction of Aop (y-axis) is plotted against that caused by induction of dfoxo (x-axis). Genes with significant differential expression (at 20% FDR) upon induction of Aop are shown in red, dfoxo in green and both Aop and dfoxo in blue. The location of the Obp99b gene on the graph is indicated. C The levels of Obp99b mRNA in abdominal fat body mRNAs were determined relative to actin by qPCR upon RU486 induction in S (n = 4) or S (n = 8) female flies. Boxplots show log-10 derived ratios scaled to set the −RU486 condition to 0. Data were analysed using a mixed-effects linear model with genotype and RU486 as main effects and dissection batch as a random effect. Both main effects as well as their interaction were significant (p<0.01), and one-tailed t-test indicated the +RU486 condition had more Obp99b transcript for each genotype (p<0.05). D Obp99b-V5 was detected in haemolymph (1 µl) of S female flies, fed or not RU486, or in whole-fly extracts (equivalent of half a fly) using anti-V5 antibody. The secreted IMP-L2 protein was used as a loading control.

To elucidate the function of Aop in the fat body we used microarray analysis on dissected tissue to identify changes in gene expression upon RU486 feeding in S females (ArrayExpress accession number: E-MTAB-927). The induction of AOPACT in only a few cells of the fat body precluded identification of many genes with significant changes upon RU486 feeding. We identified 8 genes (p value cut off of 0.00017, corresponding to 20% FDR, ). Two of these were also regulated by dfoxo ( ), and this modest overlap was statistically significant (p = 2×10−3). The functions regulated by the two factors, as determined by Catmap analysis () were also partially overlapping. For example, both dfoxo and Aop had a negative effect on respiratory electron transport chain (p = 10−21 and p = 7×10−9 respectively), while Aop alone had an effect on carbohydrate metabolism (p = 9×10−18 for Aop, p>0.05 for dfoxo) (). Interestingly, we found that the lifespan benefits from the induction of Aop, dfoxo or both were not significantly different (). This indicated that the effects of the two factors are not additive and, in turn, that the two are likely to affect lifespan via the same or related physiological processes, in an otherwise wild-type background.

Interestingly, the most strongly up-regulated gene in the fat body upon Aop induction was also up-regulated by dFOXO in the same tissue, albeit less strongly ( ). These observations were confirmed by qPCR ( ). This gene encodes a signal peptide followed by an odorant binding protein (Obp) domain - Obp99b. Non-sensory Obps are thought to bind and carry lipophilic compounds, fulfilling similar functions to the mammalian lipocalin family proteins, rather than acting in odour perception [62], [63]. Obp99b could represent a humoral factor, induced by both dfoxo and Aop, mediating inter-tissue communication. Such inter-tissue communication is known to be important in the physiological effects of tissue-restricted dfoxo induction [4]. To establish that this Obp was actually secreted in vivo we expressed a transgene encoding a V5-tagged Obp99b with the S driver. Feeding RU486 to these flies confirmed the predicted enrichment of Obp99b-V5 in the haemolymph ( ) and revealed that its transcriptional up-regulation is enough to release it in circulation. Hence, in the fat body, Aop and dfoxo regulate a common humoral factor, whose function warrants further investigation.

Discussion

Our study clearly demonstrates the in-vivo complexity of interactions that occur between FoxOs and the TFs they regulate. It emphasises the need to untangle the tissue-specific transcriptional networks within which FoxOs act in order to understand the role of FoxO factors in lifespan.

We show that dFOXO activates the transcription of an ETS repressor, Aop, in adult Drosophila gut to antagonise an ETS activator, Pnt, and avert the detrimental effects of the co-activation of dFOXO and PNT. This interaction most likely takes place at the gene promoters/enhances where AOP may directly displace PNT from a dFOXO bound region. Whether this interaction is facilitated by a direct protein-protein interaction between dFOXO and AOP, or by juxtaposition of dFOXO and AOP binding sites, remains to be determined.

Others have shown that induction of dfoxo solely in the gut is detrimental for lifespan [39]. In our outbred, wild-caught background, this detrimental effect is conditional on co-activation of Pnt. Despite these differences, the deleterious physiological outcome is likely to be based on mis-regulation of lipid metabolism genes and a drop in TAG stores in both cases. In parallel, the activation of AOP in the gut/fat body of otherwise wild-type females has only two phenotypes: increase in lifespan and increase in circulating sugars. The increase in circulating sugars is a metabolic phenotype, generally not detrimental to fly lifespan as it can be observed in a number of long-lived IIS mutants [64], [65]. Together with the effects observed upon dfoxo/Pnt co-activation, it indicates that the Aop/Pnt couple have metabolic functions in the adult fly.

FoxOs are regulated by AKT [61], while ETS factors, such as AOP and PNT, are regulated by ERK [51]. In turn, both ERK and AKT are activated in response to activation of receptor tyrosine kinases (RTK) [66]. Viewed in this context, the interaction between dFOXO, inhibited by AKT, and AOP, inhibited by ERK, insures that the two branches down-stream of RTKs are coordinated. Our findings reveal the mis-coordination of the two branches to be highly detrimental, with the cross-talk between dfoxo and Aop set up to prevent it. This mechanism for coordination of the activities of the two branches, via interactions between FoxO and ETS factors, is potentially relevant in numerous contexts beyond Drosophila lifespan.

Dissecting the relationship between dFOXO and AOP led us to the discovery of a second beneficial role for AOP and the identification of Aop as a longevity determinant in its own right. Aop has an important and well-established role in Drosophila development. The gene was initially isolated as encoding a negative regulator of neural development, in the context of the Drosphila eye [47]. Subsequent studies revealed it to be a general inhibitor of differentiation [48], [67]. Its role in the adult has not been investigated and its effect on lifespan was completely unsuspected.

Several contingent observations lead us to speculate that the role of Aop as a fly longevity gene, as well as its interaction with dfoxo, will be conserved in its mammal orthologue Etv6: Aop and Etv6 display clear conservation of function in other physiological processes [45]; the common roles of dfoxo and Aop in fly lifespan are paralleled by shared roles of FoxO factors and Etv6 in mammals [30], [46]; the neuroendocrine axes controlling growth are important in mammalian lifespan [68] and Etv6 has been identified as associated with human height in genome-wide association studies [69], [70].

Materials and Methods

Fly lines, husbandry, lifespan and physiological assays

S [1], [31], TIGS [31] UAS-dfoxo [1], UAS-GFP-dfoxo [32], UAS-Aop [48], UAS-Obp99b-V5.5 [71], UAS-Pnt [51], dfoxoΔ [6] and UAS-GFPRNAi [72] were backcrossed at least 6 times into the wild-type, outbred, Dahomey population carrying the w mutation, which had been cured of Wolbachia infection over five years ago, and frequently outcrossed back into the same wild-type population. The UAS-Aop construct was obtained from the TRiP collection (HMS01256) [49], and was tracked by PCR during backcrosses as above. The Dahomey stock was collected in 1970 in Dahomey (now Benin) and has been kept in population cages maintaining its lifespan and fecundity at levels similar to freshly caught stocks. Combinations of transgenes/mutants were created using standard fly genetic techniques while avoiding population bottlenecks. To create UAS-FLAG-Aop, the Aop open reading frame was amplified from genomic DNA of flies carrying UAS-Aop (primers are given in ) and cloned into pENTR-D-TOPO vector, confirmed by sequencing and transferred to pTFW P-element-based vector, integrated into the fly genome and the transgene backcrossed as above.

The lines were maintained, and all experiments performed, at 25°C with 60% humidity and 12 h∶12 h light∶dark cycle, on sugar-yeast-agar (1SYA) food [73]. Experimental flies developed at standardised densities and once-mated females we sorted on day two of adulthood onto food containing 200 µM RU486 (Sigma) or control food as required. Lifespans were performed as described [1]. Heamolymph extraction and measurements of circulating glucose and trahalose, and other phenotyping, were performed as described [4], [74], [75]. Examination of the “smurf” phenotype was carried out using fluorescein in 1SYA +/− RU486 as described [59], except the flies were kept on the food for 18 h.

Chromatin immunoprecipitation

Biological triplicates were done for all fly chromatin preparations. For each experiment all the batching was done so that the treatments to be compared were carried out in parallel. Chromatin was prepared from 7-day old S or S (experimental) or S or S (mock control) females fed RU486-containing food from day two of adulthood, as follows: 1000 females were crushed to a fine powder under liquid nitrogen and re-suspended in 6 ml of PBS supplemented with Protease Inhibitors Cocktail (10 µl per ml, Sigma). Cross-linking was performed with 0.5% formaldehyde for 10 min and quenched with addition of 1.5 ml of 2.5 M glycine. The cross-linked chromatin was recovered by centrifugation and washed twice with FA/SDS buffer (50 mM Hepes-KOH pH 7.5, 150 mM NaCl, 1 mM EDTA, 0.1% Na deoxycholate, 0.1% SDS, 1% Triton –X-100 and 1 mM PMSF) re-suspended in the same and incubated for 1 h at 4°C. Chromatin was again recovered by centrifugation and sheared to an average size of 400 bp by sonication, giving on average 6 ml of chromatin in FA/SDS.

For anti-GFP IPs, anti-GFP antibody (Abcam, 1 µl per IP) was bound to Protein-G Dynabeads (Invitrogen) and incubated for 2 h at room temperature with 450 µl of chromatin. Beads were washed once with FA/SDS, 3 times with FA/SDS containing 500 mM NaCl, once with 10 mM Tris-HCl pH 8, 250 mM LiCl, 1 mM EDTA, 1% NP40 and 0.5% Na deoxycholate, and once with TE. For array hybridisation, the entire IP after volume reduction, or 50 ng of total chromatin DNA, were amplified two times (Whole Genome Amplification kit, Sigma) as per manufacturer's instructions. The material from the IP was hybridised against the input material. The labelling and hybridisations were carried out by Nimblegen Systems, using custom Drosophila whole-genome tiling arrays with probes spaced approximately every 300 bp, as described [76]. ChIP-chip data were normalized using the LIMMA package [77] in Bioconductor [78], applying loess normalization within each array and quantile normalization between arrays. Replicate information was pooled by taking the median probe value for each set of arrays and was smoothed along each chromosome using a running median within a window of three probes. Experimental signal was adjusted by mock control (pre-immune serum) data by direct subtraction of median probe intensity values. Peaks were called using the Ringo package [79] in R, using a y0 threshold of 0.97 and a distance cut-off of 600 bp.

Anti-FLAG IPs were performed essentially as above, except anti-FLAG antibody (Sigma, 3 µl per IP) was used on 1350 µl of chromatin. Subsequent processing was the same except that 100 µM dUTP was spiked in during amplification, and chromatin and IP were separately hybridized to Affymetrix Drosophila Tiling 2.0R arrays. The IP to chromatin ratios were computationally determined after normalisation using the Starr package [80], and peaks identified as above, except that smoothing was over 5 probes.

Genes were associated to peaks if any gene feature occurred within 1 kb upstream or downstream of the outermost peak probes. The overlap between regions bound by GFP-dFOXO and FLAG-AOPACT was compared to the random distribution by bootstrap analysis (10000 iterations looking for marginal basepair overlap in block bootstrap) [81] and the Z-score and p value calculated.

Chromatin was also prepared from wild-type female flies and serum-starved S2 cells and IP performed with the affinity-purified anti-dFOXO antibody [32] as described [17]. Primers used for qPCR are given in .

Microarray and mRNA expression analysis

For microarray analysis of dfoxo regulated genes at least four and up to five biological repeats of 7-day old S or S females that had been fed RU486-containing or control food, were dissected removing either the mid-gut (50 per sample), or the abdominal/thorasic fat body (25 per sample) as associated with the cuticle, and the RNA extracted with RNeasy (Qiagen). The fat body dissections were performed by isolating and opening the fly thorax and/or abdomen, removing the internal organs with minimal disturbance of the fat body and keeping the fat body associated with the cuticle. We estimate that such dissection of the fat body results in ∼80% purity. For microarray analysis of Aop regulated genes four biological repeats of 7-day old S females that had been fed RU486-containing or control food, were dissected removing the abdominal fat body (25 per sample) as associated with the cuticle and the RNA extracted with RNeasy (Qiagen). All dissections and subsequent processing were done in batches so that the samples to be compared were always processed in parallel. The RNA was further processed into cRNA using standard Affymetrix protocols and hybridized to the Affymetrix Drosophila Genome 2.0 Genechip at SHWFGF/Glasgow Polyomics (University of Glasgow). The data were analysed in R. They were summarised and normalised using RMA as implemented in the LIMMA package [82]–[84]. Differential expression was assessed using linear models and the empirical Bayes moderated t-statistic implemented in LIMMA. For dfoxo-regulated gene analysis the effects in the linear model were: “dissection batch”, “S+RU486”, “S” and “S+RU486”; for Aop-regulated genes: “dissection batch” and “S+RU486”. Present/Absent calls were performed with Mas5 and transcripts present on at least four arrays for fat body, or at least five for gut were kept for further analysis. FDR was controlled using the described procedure [85]. For Gene Ontology enrichment, Catmap [40], or DAVID EASE [86] analyses were performed.

For qPCR analysis RNA was isolated from five to 10 midguts, five fat bodies or five guts and fat bodies with TRIZOL, converted into cDNA and qPCR performed as described [17]. The primers are given in .

Boolean network modeling

The network was modelled using BoolNet in R [87].

Immunofluorescence and western blots

Abdominal fat bodies as associated with the cuticle were dissected and immunofluorescence was performed using a monoclonal mouse anti-AOP antibody [48] at 1∶100 dilution using the described protocol [74]. For GFP detection, dissected body parts were fixed for 10 min before visualisation. Images were captured on Zeiss LSM 700, and quantified using ImageJ. For western blots, proteins were extracted from whole flies, or from 10 fat bodies, with TCA, or obtained by denaturing adult haemplymph in Lemmini buffer [74], separated by SDS-PAGE, transferred to nitrocellulose and probed with anti-V5 antibody (1∶500, Sigma), anti-IMP-L2 (1∶2500) [74], anti-ACTIN (1∶1000, Abcam) or anti-AOP (1∶500).

Statistical analysis

Statistical analyses were performed in JMP (version 9) software (SAS Institute), R or Excel (Microsoft). qPCR data were log-transformed to fit a normal distribution and analysed with ANOVA, linear models or mixed effects linear models, followed by selected pair-wise comparisons using t-test. Survival analysis was performed with Log-rank, CPH or mixed-effects CPH methods in JMP, or in R using survival and coxme packages (Terry Therneau, http://CRAN.R-project.org/package=survival, http://CRAN.R-project.org/package=coxme). Analysis of smurf data was performed with a generalised linear model in R and quasibinomial distribution. Analysis of feeding was performed in JMP, using a generalised linear model, binomial distribution and over-dispersion adjustment. Significance of set overlaps was determined with Hypergeometric distribution. Details of tests performed are given in figure legends.

Array data

The unprocessed array data are available from ArrayExpress under accession numbers: E-MTAB-1020, E-MTAB-1021, E-MTAB-927, and E-MTAB-1306.

A dfoxo mRNA was quantified relative to Act by qPCR in the guts or fat bodies of S or TIGS>dfoxo females, induced or not with RU486. Boxplots show log-10 derived relative expression with − RU486 values set to zero. Gut: data were analysed with a linear model and the effect of genotype was not significant (p = 0.08), while the effect of RU486 (p<10−4) and the interaction (p = 0.02) were significant. Each −RU486 was significantly different to its respective +RU486 condition, and the two +RU486 conditions were also significantly different (t-test, n = 3–6, p<0.05). Fat body: The −RU486 and +RU486 conditions were significantly different (t-test, n = 3–5, p = 8×10−4) B Correlation of log-transformed, scaled, raw intensity values (ratio of ChIP to input) for the three biological repeats of anti-GFP ChIP-chip performed on chromatin from S + RU486 (ChIP) or S + RU486 (mock) females for the probes within GFP-dFOXO bound peaks. C The sites bound by GFP-dFOXO when induced by the S driver in the gut and fat body were compared to the locations bound by endogenous dFOXO in whole flies [17]. D The distribution across genomic features of the GFP-dFOXO binding in the gut and fat body was compared to the random distribution determined from simulation of 1000 random peak sets, of identical size, length and chromosomal distribution and the Z-scores calculated. The frequency of occurrence of all the shown features was significantly different from random (p<10−3). Upstream and downstream refer to 1 kb from an annotated gene. “Genes” refers to regions containing annotated genes, as opposed to gene-free regions. E Motifs identified by MEME from the sequences bound by dFOXO in the gut and fat body (above: E value = 1.6×10−264, below: E value = 9.7×10−7) that were similar to other forkhead motifs (e.g above: Foxl1 secondary p = 2×10−3, below: Foxj1 primary p = 4.5×10−7). F The relative expression levels of the 5 TFs directly regulated by dfoxo induction in the gut. The expression levels are from the microarray experiment profiling gut mRNAs in S females with and without RU486 and are scaled to the − RU486 control for each gene. G The intron-containing dfoxo mRNA was quantified relative to Actin mRNA in whole S flies in presence or absence of RU486. t-test indicated significant difference between the two conditions (n = 4, p = 0.01). Note the mRNA expressed from the transgene is intron-less.

(PDF)

Click here for additional data file.

A Aop mRNA was quantified relative to Actin mRNA in guts and fat bodies (combined) of S flies in presence or absence of RU486. t-test indicated significant difference between the two conditions (n = 4, p = 0.04). B Survival of female flies expressing RNAi construct targeting GFP in combination with dfoxo. Log-Rank test detected significant differences in survival with or without RU486 for S but not S females (S: p>0.05; total dead/censored: − RU486 117/1, + RU486 114/1; median/maximum lifespan: − RU486: 73/85, + RU486 76/87; S: p<3×10−3; total dead/censored: − RU486 114/1, + RU486 136/5; median/maximum lifespan: − RU486: 78/90, + RU486 87/97). CPH analysis revealed no significant effect of RU486, significant effect of dfoxo (p = 7×10−4) and significant interaction (p = 5×10−5). C Aop, dfoxo and Pnt mRNA was quantified relative to Actin mRNA in S and S females in presence or absence of RU486. Boxplots show log-10 derived relative expression with − RU486 values set to zero. Data (n = 4) for each transcript were analysed with a linear model. Aop: There was a significant effect of genotype (p = 0.003) and significant effect of RU486 by genotype interaction (p = 0.003), where the levels in S + RU486 were significantly different to all others (p<0.05). dfoxo: Only the effect of RU484 was significant (p = 0.003). Pnt: Only the effect of RU486 was significant (p<10−4). Hence, the induction of dfoxo and Pnt is not different between the genotypes.

(PDF)

Click here for additional data file.

A GFP expression was visualised in midguts of S and TIGS>GFP females. Red squares denote the anterior and the mid regions of the midgut were both drivers activate transgene expression. The settings used to capture images for the two drivers are the same and the GFP intensities between the two drivers can be compared. B Parts of the anterior midgut showing GFP expression in individual cells observed in S and TIGS>GFP females. Both drivers express in the enterocytes, which can be recognised by the intense DAPI staining and large nuclei. Note the settings used to capture the images for S and TIGS>GFP females were not the same so that the intensity of GFP cannot be compared between drivers. C Survival of female flies carrying TIGS alone in the presence or absence of RU486. Log-rank detected no significant differences (p>0.05, n≈130). D Survival of female flies carrying UAS-Pnt alone in the presence or absence of RU486. Log-rank detected no significant differences (p>0.05, n≈140).

(PDF)

Click here for additional data file.

A Correlation of log-transformed ratios (ChIP to input) of scaled, raw intensity values for the three biological repeats of anti-FLAG ChIP-chip performed on chromatin from S + RU486 (ChIP) or S + RU486 (mock) females for the probes within FLAG-AOPACT bound peaks. B The distribution across genomic features of the FLAG-AOPACT binding in the gut and fat body was compared to the random distribution by bootstrap analysis [81] and the Z-scores calculated. The frequency of occurrence of all the shown features was significantly different from random (p<10−3). Upstream and downstream refer to 1 kb from an annotated gene. “Genes” refers to regions containing annotated genes, as opposed to gene-free regions. C The frequency of the lengths of regions bound by GFP-dFOXO (green) or FLAG-AOPACT (red) was plotted as smoothed density. D Percentage of regions bound by GFP-dFOXO or FLAG-AOPACT in adult gut/fat body that overlap the regions bound by AOP in larvae [58]. “6×” and “stage 11” refer to top 3% of the peaks from two different experiments described by Webber and colleagues [58].

(PDF)

Click here for additional data file.

A lipA mRNA was quantified relative to Act by qPCR in the females of the indicated genotypes, induced or not with RU486. Boxplots show log-10 derived relative expression with − RU486 values set to 0. Data were analysed with a linear model and the effects of RU486, genotype and their interaction were significant (n = 4, p<10−3), however, the S + RU486 condition was not significantly different from the S + RU486 condition (t-test , p = 0.2). B The proportion of “smurfs” (flies with impaired gut barrier function, n = 80–95) in the noted genotypes after 3 weeks of RU486 feeding, or the − RU486 condition. Data were analysed with a generalised linear model and quasibinomial distribution, and no significant effect of RU486, genotype or their interaction was detected (p>0.25). C Proportion of S females feeding after 5 days of feeding on RU486 or control food. Data were analysed with a generalised linear model, binomial distribution adjusted for over-dispersion, and no significant effect of RU486 was detected (p = 0.1). D Starvation survival of S and S females after 5 days of feeding on RU486 or the control food. RU486 significantly reduced the survival of S females (Log-rank, n∼100, p = 4×10−4) but not of the S females (Log-rank, n∼100, p = 0.4). CHP analysis detected a significant effect of RU486 (p = 9×10−3), but marginal effect of genotype (p = 0.06) and marginal RU486 by genotype interaction (p = 0.08).

(PDF)

Click here for additional data file.

A Survival of S-alone control females in the absence or presence of RU486. Log-rank test detected no significant differences (p>0.05, n≈150). B Survival of UAS-Aop-alone control females in the absence or presence of RU486. Log-rank test detected no significant differences (p>0.05, n≈150). C Lifespan of S female flies in presence or absence of RU486. The two conditions were different by Log-rank test (p = 3×10−5; total dead/censored: − RU486 143/4, + RU486 147/0; median/maximum lifespan: − RU486: 71/83, + RU486 76/87). D Aop mRNA was quantified relative to Act by qPCR in the guts of S or TIGS>Aop females, induced or not with RU486. Boxplots show log-10 derived relative expression with −RU486 values set to zero. Data (n = 2–4) were analysed with a linear model and the effects of genotype, RU486 and their interaction were all significant(p<10−4). Each + RU486 was significantly different to its respective − RU486 condition, and the two + RU486 conditions were also significantly different (t-test, p<0.05).

(PDF)

Click here for additional data file.

S female flies were placed on appropriate food on day two and either frozen on day 7 for metabolic assays or kept for stress and fecundity assays. For stress assays, 7-day old females were placed on food containing H2O2 (no yeast) or DDT, or were starved on food containing agar alone, and the number of dead flies scored over time. There were no significant differences in survival by Long-rank test (n≈100, p>0.05) in any of the conditions. Levels of trahalose were determined as number of moles of trehalose per fly weight, of lipids as weight of TAG per fly weight, of glycogen as weight of glycogen per fly weight and all are shown as percentage of no-RU486 control. There were no significant differences by t-test (n = 10, p>0.05). Eggs laid per female per day (averaged per vial of 10 females) were determined once per week for the first 4 weeks and summed to give an estimate of lifetime fecundity. There were no significant differences (t-test, n = 10, p>0.05). Feeding was assessed in 8-day old females and no significant differences were detected.

(PDF)

Click here for additional data file.

A Aop mRNA (left) or AOP protein (right) levels were determined against Act/ACT control by qPCR or western blots in the fat bodies of S females fed RU486 or control food. qPCR data were log-10 derived and significant difference found between −RU486 and + RU486 conditions (t-test, p = 8×10−4). B Quantification of the results presented in Figure 8A. AOP was visualised by immunofuorescence in fat bodies of S flies induced or not with RU486. Intensity of nuclear AOP staining was quantified from confocal images using Image J (intensity is presented on an arbitrary scale). Average intensity from at least 3 cells from each animal was calculated for each biological repeat, and location of the nucleus was determined from DAPI staining. For + RU486 samples, the nuclei were binned into “high” or “low” AOP stained ones by visual inspection and quantified separately. Means ± SEM are shown. n = 3 animals − RU486, n = 4 for + RU486; t-test detected a significant difference (p<10−4) between the “high AOP” nuclei on + RU486 and nuclei on − RU486. The proportion of highly stained nuclei was 26%+/−6% in 4 animals fed RU486 examined. C S drives expression in all the cells of the fat body. S flies were fed or not RU486 from day two until day 7 when GFP expression in the fat body was determined by confocal microscopy. GFP was detected uniformly in all the cells of the abdominal fat body.

(PDF)

Click here for additional data file.

A Biological process GO categories differentially regulated (p<10−10) in the fat body upon induction of dfoxo or Aop as determined by Catmap analysis. Any redundant categories (overlap by more than 75%) were removed, retaining the most specific category. The full list is given in . The intensity of red shows the log10-transformed p-value associated with differential regulation for each category. B Median lifespan extension caused by induction or dfoxo, Aop or both with RU486 and the S driver in two independent experimental trials. In each case RU486 had a significant, positive effect on survival (Log-rank, p<10−3). The survival data (n = 1664 deaths/23 censors) were analysed with a mixed effects CPH model, with experimental trial as random effect, and the effect of genotype (p<0.05) and RU486 (p = 10−12) were significant but their interaction was not (p>0.1), revealing that the two factors do not have additive effects.

(PDF)

Click here for additional data file.

Lists of genes and genomic locations and full Ease and Catmap results mentioned in the manuscript.

(ZIP)

Click here for additional data file.

Sequences of primers used in this study.

(PDF)

Click here for additional data file.