Evolution of longevity improves immunity in Drosophila

Abstract Much has been learned about the genetics of aging from studies in model organisms, but still little is known about naturally occurring alleles that contribute to variation in longevity. For example, analysis of mutants and transgenes has identified insulin signaling as a major regulator of longevity, yet whether standing variation in this pathway underlies microevolutionary changes in lifespan and correlated fitness traits remains largely unclear. Here, we have analyzed the genomes of a set of Drosophila melanogaster lines that have been maintained under direct selection for postponed reproduction and indirect selection for longevity, relative to unselected control lines, for over 35 years. We identified many candidate loci shaped by selection for longevity and late‐life fertility, but – contrary to expectation – we did not find overrepresentation of canonical longevity genes. Instead, we found an enrichment of immunity genes, particularly in the Toll pathway, suggesting that evolutionary changes in immune function might underpin – in part – the evolution of late‐life fertility and longevity. To test whether this genomic signature is causative, we performed functional experiments. In contrast to control flies, long‐lived flies tended to downregulate the expression of antimicrobial peptides upon infection with age yet survived fungal, bacterial, and viral infections significantly better, consistent with alleviated immunosenescence. To examine whether genes of the Toll pathway directly affect longevity, we employed conditional knockdown using in vivo RNAi. In adults, RNAi against the Toll receptor extended lifespan, whereas silencing the pathway antagonist cactus‐–causing immune hyperactivation – dramatically shortened lifespan. Together, our results suggest that genetic changes in the age‐dependent regulation of immune homeostasis might contribute to the evolution of longer life.


Impact Summary
Despite much progress in our understanding of the genetic basis of aging, mainly from studying large-effect mutants, little is known about natural variants that contribute to the evolution of lifespan and related fitness traits. To identify the mechanisms by which longevity evolves, we sequenced a set of D. melanogaster populations that have been undergoing selection for late-life reproduction and postponed senescence, relative to unselected controls, for over 35 years. Instead of an enrichment of evolutionary changes in previously identified "canonical" longevity genes, we found an enrichment of genetically diverged immunity genes, suggesting that variation in immune function contributes to the evolution of lifespan and late-life fertility. To test this hypothesis, we employed immunity assays: longlived flies survived infections better and showed altered age-dependent immune gene expression as compared to control flies. Using in vivo RNAi we confirmed that reduced expression of immune genes extends lifespan while immune overactivation is strongly detrimental.
Despite major progress in our understanding of the genetic basis of aging and life history, especially in model organisms such as yeast, C. elegans, Drosophila, and mice (Guarente and Kenyon 2000;Partridge and Gems 2002;Tatar et al. 2003;Guarente et al. 2008;Kenyon 2010;Flatt and Heyland 2011), the identity and effects of naturally segregating polymorphisms that affect variation in lifespan and correlated fitness traits and which might thus underpin the evolution of longevity and life history remain poorly understood to date (De Luca et al. 2003;Pasyukova et al. 2004;Carbone et al. 2006;Flatt and Schmidt 2009;Paaby et al. 2014;Carnes et al. 2015;Flatt and Partridge 2018).
Several major evolutionarily conserved pathways that regulate lifespan and correlated fitness traits, including insulin/insulinlike growth factor 1 signaling (IIS), have been identified using analyses of large-effect mutants and transgenes in the laboratory (Partridge and Gems 2002;Tatar et al. 2003;Kenyon 2010), but to what extent genes in these "canonical" pathways harbor segregating alleles that affect lifespan is mostly unknown (Flatt and Schmidt 2009;Paaby et al. 2014;Carnes et al. 2015;Flatt and Partridge 2018). For instance, only few studies to date have identified functional effects of segregating IIS polymorphisms upon lifespan and correlated life-history traits in populations of Drosophila (Paaby et al. 2010(Paaby et al. , 2014Remolina et al. 2012) or which contribute to longevity in human centenarians (Suh et al. 2008; Willcox et al. 2008;Flachsbart et al. 2017;Joshi et al. 2017).
Here, we take advantage of a >35-year-long laboratory selection experiment for late-life fertility and increased lifespan in Drosophila melanogaster, first published by colleagues in 1984 (Luckinbill et al. 1984; also see Luckinbill and Clare 1985;Arking 1987), to analyze the genomic footprints underlying the evolution of delayed reproduction and postponed aging. In this long-term selection experiment, replicate lines derived from an outbred base population have been selected for late-life fertility and-indirectly-for increased lifespan by breeding only from flies that survived and were fertile at a relatively old age. In contrast, unselected replicate control lines have been propagated across generations by breeding from flies with a random age at reproduction (for details see Supplementary methods). Selected flies in this experiment have evolved late-life fertility and live 40-50% longer than unselected control flies, yet exhibit reduced early fecundity relative to the controls (see Supplementary methods). Thus, these selection lines are subject to a genetic trade-off between late-life performance (long life, late-life fertility) and early fecundity, as is commonly observed in laboratory evolution experiments that directly or indirectly select for changes in Drosophila lifespan (Luckinbill et al. 1984;Rose 1984;Zwaan et al. 1995;Partridge et al. 1999;Stearns et al. 2000;Remolina et al. 2012).
The central finding from our genomic analysis of this selection experiment is that evolutionary changes in innate immunity contribute to the evolution of late-life performance in fruit flies, probably by improving age-dependent immune homeostasis. Although still little is understood about the mechanistic interplay between immunity and aging (Garschall and Flatt 2018), our analyses suggest that immune function is a major longevity assurance mechanism that can be targeted by selection on standing genetic variation.

OF LONGEVITY
To characterize the genomic signature of longevity we used nextgeneration pool-sequencing (Pool-seq) (Schlötterer et al. 2014) to obtain genome-wide allele frequency estimates from four longlived selection lines and two unselected control lines after ࣙ 144 generations of selection (see Supplementary methods for details). We identified candidate SNPs by comparing allele frequency differentiation between the selection and control regimes with a stringent F ST outlier approach (Lewontin and Krakauer 1973;Akey 2009) (Fig. 1A,B). The majority of SNPs (62.2%) showed no or less differentiation between the selection versus control regime as compared to differentiation within these regimes (selection signal-to-noise ratio ࣘ 0; Fig. 1B,C). We defined SNPs as candidates if they showed very strong, consistent verify that this number of candidate SNPs is not due to chance we applied our candidate criteria to all 6435 possible sets of eight pairwise comparisons; out of these combinations only one set is biologically informative in terms of inferring selection, that is the set of all eight pairwise control versus selection comparisons (see Supplementary methods). No combination of eight pairwise comparisons yielded as many candidate SNPs as this "true" set of comparisons (red bar), with a probability that the "true" number of candidate SNPs is due to chance of P 1.6 × 10 −4 . and significant differentiation in all eight pairwise comparisons between the four selection and two control lines (signal-tonoise ratio 0.9; F ST(selection vs. control) > 0.9; Bonferroni-corrected Fisher's exact test: P < 10 −9 ) (Fig. 1A,B,C). Using this approach, we identified 8205 candidate SNPs in 868 genes dis-tributed across the entire genome ( Fig. 1B; Table S1; genes were defined as the sequence between the ends of the 5' and 3' UTRs plus 1 kb up-and downstream; also see Supplementary methods). Candidate loci appeared to cluster non-randomly in specific genomic regions, suggesting pervasive polygenic selection and/or indirect selection due to "hitchhiking" ("genetic draft") ( Fig. 1B; Table S1). To further validate our set of longevity candidate SNPs and to exclude false positives due to randomness, for example because of genetic drift, we used a combinatorial approach (see Supplementary methods). We found that-when applying our stringent candidate criteria -it is highly unlikely (P 1.6 × 10 −4 ) that this large number of candidate SNPs arose by chance (Fig. 1D).

PARALLELISM
While some mechanisms of longevity are evolutionarily conserved ("shared") among species and thus "public," for example insulin/insulin-like growth factor 1 signaling (IIS), most others are likely to be lineage-specific and thus 'private' (Martin et al. 1996;Partridge and Gems 2002;McElwee et al. 2007). Similarly, at the intraspecific level, parallel and convergent evolution in independent populations might result in the repeated use of the same genes underlying a given trait ("gene reuse") (Conte et al. 2012), but to what extent this might be the case for longevity remains unclear. Addressing this question might give insights into the predictability of the evolution of lifespan at the genetic level (Stern and Orgogozo 2008;Conte et al. 2012).
To examine how frequently the same genes are used by different populations during the evolution of late-life fertility and longevity, we compared our list of candidate genes to those from two other "Evolve and Resequence" studies of Drosophila longevity and correlated life-history traits (Remolina et al. 2012;Carnes et al. 2015). The study by Carnes et al. (2015) provides a genomic analysis of an independent long-term selection experiment by Rose (Rose 1984) similar in duration to ours (Luckinbill et al. 1984), with both selection experiments first published back-to-back in 1984. The other study, by Remolina et al. (2012), performed whole-genome sequencing of a shorter, 50-generationlong selection experiment for longevity. Importantly, both Rose (1984) and Remolina et al. (2012) selected for increased lifespan by postponing reproduction, using a design that is qualitatively identical to ours.
We discovered statistically significant sharing of candidate loci across all possible overlaps among the three datasets (Fig. 2, Table S2), indicating genetic parallelism underlying the evolution of late-life performance. Our dataset contained 147 (11.7%) of the candidate genes of Carnes et al. (2015) and 102 (10.9%) of those of Remolina et al. (2012). Twenty candidate genes (ß2%) were shared across all three studies, representing clear cases of gene reuse during the evolution of longevity and late-life fertility (Fig. 2, Table S2). Thus, as might be expected from a highly complex and polygenic trait such as lifespan (McElwee et al. 2007), most candidate loci tend to be population-specific. However, a small but significant proportion of candidate loci is shared among independent populations, perhaps suggesting the existence of  Table S2 for functional annotations of the shared longevity candidate genes; see Table S5 for statistical details.
"preferred" loci of evolutionary change (Stern and Orgogozo 2008) for longevity. Several of these "high confidence" genes represent promising candidate loci for future functional experiments.
Notably, although each study identified several loci that belong to "canonical" longevity pathways (Guarente and Kenyon 2000;Partridge and Gems 2002;Tatar et al. 2003;Guarente et al. 2008;Kenyon 2010), for example the IIS pathway, the candidate lists and overlaps contain few "classical" lifespan genes that have previously been identified in studies of large-effect mutants and transgenes. This might be due to a lack of standing variation at these "canonical" longevity loci: perhaps these conserved-effect loci have been optimized by selection but are now subject to strong purifying selection (see Remolina et al. 2012;Flatt and Partridge 2018). Thus, while segregating IIS polymorphisms with major effects on life-history traits including lifespan have been identified (Geiger- Thornsberry and Mackay 2004;Paaby et al. 2010Paaby et al. , 2014Flachsbart et al. 2017;Joshi et al. 2017), our results are consistent with the hypothesis that loci in these canonical pathways might be under selective constraints (see Remolina et al. 2012;Flatt and Partridge 2018).
Even though "canonical" longevity loci seem to be underrepresented, many of the overlapping candidate genes that we have identified have strong empirical support from functional genetics, GWAS, QTL, or gene expression studies, with known roles in lifespan determination, somatic maintenance (e.g., resistance against starvation or oxidative stress, immunity, metabolism), and age-specific fecundity (see functional annotations in Table S2). The fact that several candidate loci are known to affect age-specific fecundity is consistent with the age-at-reproduction selection regime used by all three studies and possibly also with genetic trade-offs between early fecundity and lifespan (and/or late-life fecundity) seen in these selection experiments.

IMMUNE FUNCTION
We next sought to characterize the functions of our candidate loci with gene ontology (GO) analysis (Kofler and Schlötterer 2012) (Table S3; considering the ontologies "Biological Function," "Molecular Function," and "Cellular Component"). Interestingly, we found an enrichment of candidate genes associated with "antifungal peptides" with a false discovery rate of ß9% (FDR = 0.085), whereas the term "determination of adult lifespan" had no support (FDR = 1) (Table S3). Immunity against fungi (and gram-positive bacteria) is regulated by Toll signaling (Belvin and Anderson 1996;Lemaitre et al. 1996;De Gregorio et al. 2002;Valanne et al. 2011), and among our candidates we identified several prominent members of this pathway, including the Toll ligand spätzle (spz), the receptor Toll (Tl), the Toll inhibitor cactus (cact), the NFκB transcription factors Dorsalrelated immunity factor (Dif) and dorsal (dl), the upstream serine proteases persephone (psh) and sphinx2, and two regulators of cactus, scalloped (sd) and cactin (Fig. 3, Table S4). The other major immune pathway, the Imd pathway (De Gregorio et al. 2002;Kleino and Silverman 2014;Myllymäki et al. 2014), also harbored several but fewer candidates, including peptidoglycan recognition protein LE (PRGP-LE) and the antimicrobial peptide Cecropin A1 (CecA1) (Fig. 3, Table S4).
The enrichment of immunity genes prompted us to hypothesize that genetic changes in immune function might contribute to the evolution of longevity and correlated fitness traits (DeVeale et al. 2004;Finch 2007). Importantly, Remolina et al. (2012) also found enrichment of genes involved in "defense response to fungus," and Carnes et al. (2015) observed divergence in immune gene expression between long-lived selection and control lines, suggesting that the relation between immunity and lifespan might be general (DeVeale et al. 2004;Finch 2007). While we found a larger number of genes in the Toll pathway, Carnes et al. (2015) and Remolina et al. (2012) found more candidates in the Imd pathway. However, several immune genes are shared across the three studies, despite a relatively small overlap at the individual gene level (Table S4). Immunity might thus represent a general mechanism underlying longevity, with immune genes having pleiotropic effects on lifespan and correlated fitness components.
Despite this compelling commonality across independent experiments, still little is known about how immunity proximately affects longevity and correlated fitness traits; similarly, whether genetic changes in immunity might contribute to the evolution of longer life remains unknown (Garsin et al. 2003;DeVeale et al. 2004;Kurz and Tan 2004;Libert et al. 2006;Troemel et al. 2006;Libert et al. 2008;Fernando et al. 2014;Guo et al. 2014;McCormack et al. 2016;Kounatidis et al. 2017;Loch et al. 2017;Yunger et al. 2017). We therefore aimed to test whether the evolved genomic signature of immune gene enrichment observed in our study -and similarly by Carnes et al. (2015) and Remolina et al. (2012)-might represent a physiological mechanism underlying evolutionary changes in lifespan and late-life fertility.

INDUCTION WITH AGE
We first examined whether the selection and control lines differ in the expression of antimicrobial peptides (AMPs), the major effectors of the innate immune response. We used three AMPs as readouts of Toll and Imd signaling activity, Drosomycin (Drs), Attacin A (AttA), and Diptericin (Dpt). Drs and AttA are regulated by both Toll and Imd signaling, whereas Dpt is mainly regulated by the Imd pathway (De Gregorio et al. 2002). Using quantitative realtime PCR, we determined mRNA levels of young (5-6-day-old) and aged (25-26-day-old) female flies, either without pricking, upon aseptic pricking (mock control) or upon prick infection with Erwinia carotovora carotovora 15 (Ecc 15). Systemic infections with this bacterium induce the expression of all three AMPs assayed here (Lemaitre et al. 1997;Basset et al. 2000;De Gregorio et al. 2002).
Without pricking, control flies upregulated AMP baseline expression with age (Fig. 4A) -a pattern that is commonly observed in wild-type flies and attributed to persistent chronic infection and a prolonged immune response at old age (Seroude et al. 2002;DeVeale et al. 2004;Zerofsky et al. 2005;Ren et al. 2007;Ramsden et al. 2008). In marked contrast to control flies, baseline AMP levels remained constant as a function of age in selected flies (Fig. 4A).
AMP expression also differed substantially between control and selected flies upon infection: at young age, the AMP response was slightly stronger in long-lived flies than in control flies, whereas at old age long-lived flies tended to downregulate AMP induction (Fig. 4b). Thus, unlike aged wild-type flies which upregulate AMPs but suffer from immunosenescence and show signs consistent with chronic inflammation (i.e., reduced infection survival, increased bacterial load, more persistent AMP induction upon infection; see Zerofsky et al. 2005;Ren et al. 2007;Ramsden et al. 2008;Myllymäki et al. 2014 Figure 3. Genes of the Toll and Imd pathways represent longevity candidates. Overview of the Toll and Imd pathways, the two major pathways regulating the humoral innate immune response against fungi and gram-positive bacteria (Toll) and gram-negative bacteria (Imd). Among our longevity candidates we found an enrichment of immunity-related genes (enrichment of GO terms associated with "antifungal peptides"). Longevity candidate genes identified in the Toll and Imd pathways are shown in red. For additional immunityrelated candidate genes see Table S4.
Our results therefore suggest that long-lived flies might have evolved improved age-dependent immune homeostasis and alleviated immunosenescence (DeVeale et al. 2004). These evolutionary changes in immune gene induction might also be linked to the late-life fertility of the long-lived lines. Since in our selection experiment lifespan was selected for by postponing reproduction, the observed differences in immune gene induction between the regimes might be a byproduct of selection for increased late-life fertility in the long-lived selection lines. This would be consistent with the observation that infection reduces fecundity: infection-induced synthesis of AMPs incurs a cost of reproduction in wild-type flies but this cost is abolished in Imd pathway mutants (Zerofsky et al. 2005).

UPON INFECTION
To investigate whether selected and control flies differ in realized immune function we measured their survival after infection with four different pathogens (Fig. 5). Long-lived flies survived infections with a fungus (Beauveria bassiana, Bb), with the Gramnegative bacterium Ecc15 and with the Gram-positive bacterium Enterococcus faecalis (Ef) overall markedly better than control flies (Fig. 5A,B,C,E). Improved survival of long-lived flies was observed for both young and aged flies after infection with Bb and Ecc15, whereas for Ef infection only aged long-lived flies showed increased survival relative to controls (Fig. 5A,B,C,E). Because one of our candidate genes, the JAK/STAT activating cytokine un-paired3 (upd3 ; Table S4), is involved in antiviral immunity (Zhu et al. 2013), we also measured the survival of flies upon infection with Drosophila C virus (DCV). This assay was carried out only with young, not aged flies, but we again found that long-lived flies survived infection with DCV much better than control flies (Fig. 5D,E). The evolution of prolonged lifespan might thus be accompanied-or partly be caused -by selection for improved realized immunity.
Next, we examined the ability of selection and control flies to successfully clear bacterial (Ecc15) infections over a 6-day period postinfection. The ability of control flies to clear an infection was  Table S5.
higher than that of long-lived flies at young age but declined at old age; in contrast, clearance was overall lower in long-lived flies yet did not change with age (Fig. 5F). The lower clearance ability of long-lived selected flies, independent of their age, together with their improved survival upon infection, possibly indicates that they have evolved to be more tolerant to infections than unselected control flies (Best et al. 2008;Schneider 2008, 2012;Felix et al. 2012).

OVERACTIVATION IS DETRIMENTAL
Our results above support the idea that improved age-dependent regulation of immunity contributes to longevity and late-life fertility, but how immune genes affect lifespan is not well studied, especially in Drosophila (DeVeale et al. 2004;Libert et al. 2006;Fernando et al. 2014;Guo et al. 2014;Kounatidis et al. 2017;Loch et al. 2017). For example, previous work has shown that constitutive upregulation of the peptidoglycan recognition proteins PGRP-LE and PGRP-LC causes hyperactivation of Imd signaling and reduces lifespan (DeVeale et al. 2004;Libert et al. 2006). Similarly, several mutants of negative regulators of Imd signaling display shortened lifespan (Fernando et al. 2014;Kounatidis et al. 2017). While we also identified PGRP-LE as a lifespan candidate gene, most immunity genes in our analysis belong to the Toll pathway (Fig. 3, Table S4).
To examine whether Toll signaling affects lifespan, we used transgenic RNAi to silence four longevity candidate genes of the Toll pathway: the ligand spz, the receptor Tl, the inhibitor cact, and the NFκB transcription factor Dif. To prevent deleterious side effects of knocking down these developmentally critical genes (Nüsslein-Volhard and Wieschaus 1980;Belvin and Anderson 1996) we used a mifepristone-inducible daughterless (da)-GeneSwitch(GS)-GAL4 driver (Tricoire et al. 2009) to direct expression of UAS-RNAi constructs against  Table S5.
these genes specifically during adulthood and throughout the fly body. Downregulation of the Tl receptor-but not of its ligand spzmildly but significantly extended lifespan (Fig. 6A,B,C,D), while silencing the antagonist cact-leading to Toll pathway hyperactivation (Lemaitre et al. 1996;Aggarwal and Silverman 2008)drastically reduced lifespan (Fig. 6E,F), similar to the effects of overactivation or derepression of Imd signaling (DeVeale et al.  Table S5. 2004; Libert et al. 2006;Guo et al. 2014;Kounatidis et al. 2017). Interestingly, we found opposite lifespan effects of Dif-RNAi for females (Fig. 6G) and males (Fig. 6H). In agreement with our findings for females, two studies have previously found that a loss-of-function mutant of Dif lives longer than wild-type (Le Bourg et al. 2012;Petersen et al. 2013), but why silencing Dif reduces male lifespan remains unclear. Our results thus establish that downregulation of Toll signaling increases lifespan (albeit weakly so), whereas overactivation of this pathway strongly shortens life.

Expression of the different UAS-RNAi responder constructs was driven with a mifepristone-inducible daughterless(da)-GeneSwitch(GS)-
Our findings for the Toll pathway are also consistent with recent studies of IMD signaling showing that lifespan is extended under conditions of reduced lifetime IMD activity (Loch et al. 2017) or when the IMD AMPs AttacinC (AttC) and Diptericin B (DiptB) are downregulated in the fat body (Lin et al. 2018). The evidence available to date therefore suggests that decreased activity of the immune system can promote lifespan (DeVeale 2004), possibly by reducing the costs of immune deployment (McKean and Lazzaro 2011). Moreover, as we show here, longer lifespan can evolve -at least partly -via evolutionary changes in immunity.

Conclusion
Explaining the genetic basis of variation in longevity is a longstanding problem in evolutionary genetics and the biology of aging (Finch 1990;Rose 1991;Zwaan 1999;Partridge and Gems 2006;Flatt and Schmidt 2009;Flatt and Partridge 2018). Here we have performed a whole-genome sequencing analysis of an over 35-year-long selection experiment for postponed aging and late-life fertility in Drosophila (Luckinbill et al. 1984).
Notably, among the longevity candidate genes identified in our genomic screen, we found an enrichment of immune genes, especially in the Toll pathway. By comparing our data to those from two previous genomic studies of longevity selection in Drosophila (Remolina et al. 2012;Carnes et al. 2015) we infer that-while different studies might identify different immune genes as longevity candidates-immune function likely represents a general processlevel mechanism underlying the evolution of longevity assurance and of late-life performance (Martin et al. 1996;Partridge and Gems 2002;McElwee et al. 2007). This is particularly noteworthy in view of the growing evidence that aging, inflammation and immunity are intricately linked at the molecular level (DeVeale et al. 2004 ;Kurz and Tan 2004;Finch 2007;Salminen et al. 2008;Eleftherianos and Castillo 2012). However, how immunity contributes to longevity and correlated fitness traits is largely unclear.
While aged wild-type flies upregulate immune gene expression (Pletcher et al. 2002;Seroude et al. 2002;Landis et al. 2004), they typically have a reduced capacity to fight off and survive infections, suggesting that they suffer from immune overactivation and immunopathology (Zerofsky et al. 2005;Ren et al. 2007;Ramsden et al. 2008). Here, we show that long-lived flies, by contrast, tend to downregulate the induction of immune effector genes (AMPs) with age yet have substantially improved survivorship upon infection. This seems to confirm that elevated immune gene expression at old age might either be ineffective or even detrimental, perhaps as a consequence of senescent dysregulation of gene expression (Zerofsky et al. 2005;Khan et al. 2017). The downregulation of AMPs seen in the long-lived selection lines might also be a byproduct of selection for late-life fertility in these lines since elevated AMP expression upon infection is known to reduce fecundity (Zerofsky et al. 2005).
Since optimal immunity depends on the balance between efficient clearance of pathogens and limiting immunity-induced damage (Cassedevall and Pirofski 1999;Read et al. 2008;Råberg et al. 2009;Medzhitov et al. 2012), we propose that selection for longevity and late-life fertility leads to improved age-dependent immune homeostasis and alleviates the trade-off between immunity and immunopathology. This trade-off can be decoupled to some degree by tolerance mechanisms (Medzhitov et al. 2012), suggesting that the improved immunity of long-lived flies mightat least in part-be due to increased tolerance. In line with the notion of a trade-off between immunity and immunity-induced damage, work in the mealworm beetle shows that deployment of the immune effector phenoloxidase (PO) causes early-life inflammation, faster aging, and immunopathology at old age, whereas RNAi silencing of PO extends lifespan and improves survival after infection (Khan et al. 2017). This is consistent with the fact that hyperactivation or derepression of Imd signaling (DeVeale et al. 2004;Libert et al. 2006;Fernando et al. 2014;Kounatidis et al. 2017) and, as we observe here, of Toll signaling reduces lifespan. Conversely, we find that adult downregulation of Toll signaling mildly promotes lifespan, similar to recent findings for the Imd pathway (Kounatidis et al. 2017;Lin et al. 2018).
Together, our work reveals the existence of a causal-but mechanistically still poorly understood-link between improved age-dependent immunity and the evolution of longevity and latelife fertility (Garschall and Flatt 2018). This relationship clearly warrants further mechanistic and evolutionary study.

Methods
All methods are given in the Supplementary methods file (see Supporting Information section below), including details of selection and control lines, next-generation sequencing, bioinformatic, and statistical analyses, gene expression analyses, immunity assays, transgenic RNAi, and lifespan assays. in the laboratory; the Bloomington Drosophila Stock Center (BDSC) and the Vienna Drosophila RNAi Center (VDRC) for fly stocks; Véronique Monnier for the da-GS-GAL4 strain; and Luis Teixeira for the DCV strain. Our work was supported by grants from the Austrian Science Foundation (FWF P21498-B11 and W1225) and the Swiss National Science Foundation (SNSF PP00P3 133641 and PP00P3 165836) to T.F. G. S.-M. was supported by a NOS Alive -IGC fellowship.

DATA AVAILABILITY
Sequencing data used for genomic analyses are available from the European Nucleotide Archive (ENA) under accession PRJEB28048 / ERP110212. Raw data for experimental assays are available from Dryad under accession https://doi.org/10.5061/dryad.cp38vj4.

CONFLICT OF INTEREST
The authors declare no conflict of interest.

Supporting Information
Additional supporting information may be found online in the Supporting Information section at the end of the article.
Supplementary methods (pdf). Description of all methods, including details of selection and control lines, next-generation sequencing, bioinformatic and statistical analyses, gene expression analyses, immunity assays, transgenic RNAi and lifespan assays. Table S1. (xls). Longevity candidate SNPs and candidate genes. Table S2. (xls). Shared candidate genes across three independent studies. Table S3. (xls). Gene ontology (GO) analysis of longevity candidate genes. Table S4. (xls). Immunity genes implicated in lifespan and aging. Table S5. (xls). Full statistical details of data analyses shown in the main text.