Genomic basis of aging and life-history evolution in Drosophila melanogaster

Ecology, the pace-of-life, epistatic selection and the maintenance of genetic variation in life-history genes

Evolutionary genetics has long struggled with understanding how functional genes under selection remain polymorphic in natural populations. Taking as a starting point that natural selection is ultimately a manifestation of ecological processes, we spotlight an underemphasized and potentially ubiquitous ecological effect that may have fundamental effects on the maintenance of genetic variation. Negative frequency dependency is a well-established emergent property of density dependence in ecology, because the relative profitability of different modes of exploiting or utilizing limiting resources tends to be inversely proportional to their frequency in a population. We suggest that this may often generate negative frequency-dependent selection (NFDS) on major effect loci that affect rate-dependent physiological processes, such as metabolic rate, that are phenotypically manifested as polymorphism in pace-of-life syndromes. When such a locus under NFDS shows stable intermediate frequency polymorphism, this should generate epistatic selection potentially involving large numbers of loci with more minor effects on life-history (LH) traits. When alternative alleles at such loci show sign epistasis with a major effect locus, this associative NFDS will promote the maintenance of polygenic variation in LH genes. We provide examples of the kind of major effect loci that could be involved and suggest empirical avenues that may better inform us on the importance and reach of this process.

Neuronal expression in Drosophila of an evolutionarily conserved metallophosphodiesterase reveals pleiotropic roles in longevity and odorant response

Evolutionarily conserved genes often play critical roles in organismal physiology. Here, we describe multiple roles of a previously uncharacterized Class III metallophosphodiesterase in Drosophila, an ortholog of the MPPED1 and MPPED2 proteins expressed in the mammalian brain. dMpped, the product of CG16717, hydrolyzed phosphodiester substrates including cAMP and cGMP in a metal-dependent manner. dMpped is expressed during development and in the adult fly. RNA-seq analysis of dMppedKO flies revealed misregulation of innate immune pathways. dMppedKO flies showed a reduced lifespan, which could be restored in Dredd hypomorphs, indicating that excessive production of antimicrobial peptides contributed to reduced longevity. Elevated levels of cAMP and cGMP in the brain of dMppedKO flies was restored on neuronal expression of dMpped, with a concomitant reduction in levels of antimicrobial peptides and restoration of normal life span. We observed that dMpped is expressed in the antennal lobe in the fly brain. dMppedKO flies showed defective specific attractant perception and desiccation sensitivity, correlated with the overexpression of Obp28 and Obp59 in knock-out flies. Importantly, neuronal expression of mammalian MPPED2 restored lifespan in dMppedKO flies. This is the first description of the pleiotropic roles of an evolutionarily conserved metallophosphodiesterase that may moonlight in diverse signaling pathways in an organism.

Polygenic adaptation to overnutrition reveals a role for cholinergic signaling in longevity

Overnutrition by high-sugar (HS) feeding reduces both the lifespan and healthspan across taxa. Pressuring organisms to adapt to overnutrition can highlight genes and pathways important for the healthspan in stressful environments. We used an experimental evolution approach to adapt four replicate, outbred population pairs of Drosophila melanogaster to a HS or control diet. Sexes were separated and aged on either diet until mid-life, then mated to produce the next generation, allowing enrichment for protective alleles over time. All HS-selected populations increased their lifespan and were therefore used as a platform to compare allele frequencies and gene expression. Pathways functioning in the nervous system were overrepresented in the genomic data and showed evidence for parallel evolution, although very few genes were the same across replicates. Acetylcholine-related genes, including the muscarinic receptor mAChR-A, showed significant changes in allele frequency in multiple selected populations and differential expression on a HS diet. Using genetic and pharmacological approaches, we show that cholinergic signaling affects Drosophila feeding in a sugar-specific fashion. Together, these results suggest that adaptation produces changes in allele frequencies that benefit animals under conditions of overnutrition and that it is repeatable at the pathway level.

A Single Nucleotide Variant in the PPARγ-homolog Eip75B Affects Fecundity in Drosophila

Single nucleotide polymorphisms are the most common type of genetic variation, but how these variants contribute to the adaptation of complex phenotypes is largely unknown. Experimental evolution and genome-wide association studies have demonstrated that variation in the PPARγ-homolog Eip75B has associated with longevity and life-history differences in the fruit fly Drosophila melanogaster. Using RNAi knockdown, we first demonstrate that reduced expression of Eip75B in adult flies affects lifespan, egg-laying rate, and egg volume. We then tested the effects of a naturally occurring SNP within a cis-regulatory domain of Eip75B by applying two complementary approaches: a Mendelian randomization approach using lines of the Drosophila Genetic Reference Panel, and allelic replacement using precise CRISPR/Cas9-induced genome editing. Our experiments reveal that this natural polymorphism has a significant pleiotropic effect on fecundity and egg-to-adult viability, but not on longevity or other life-history traits. Our results provide a rare functional validation at the nucleotide level and identify a natural allelic variant affecting fitness and life-history adaptation.

Dietary stress remodels the genetic architecture of lifespan variation in outbred Drosophila

Evolutionary theory suggests that lifespan-reducing alleles should be purged from the gene pool, and yet decades of genome-wide association and model organism studies have shown that they persist. One potential explanation is that alleles that regulate lifespan do so only in certain environmental contexts. We exposed outbred Drosophila to control and high-sugar diets and genotyped more than 10,000 adult flies to track allele frequency changes over the course of a single adult lifespan. We identified thousands of lifespan-associated alleles associated with early versus late-life trade-offs, late-onset effects and genotype-by-environment interactions. Remarkably, a third of lifespan-associated genetic variation had environmentally dependent effects on lifespan. We find that lifespan-reducing alleles are often recently derived, have stronger effects on a high-sugar diet and show signatures of selection in wild Drosophila populations, consistent with the evolutionary mismatch hypothesis. Our results provide insight into the highly polygenic and context-dependent genetic architecture of lifespan variation and the evolutionary processes that shape this key trait.

Experimental Life History Evolution Results in Sex-specific Evolution of Gene Expression in Seed Beetles

The patterns of reproductive timing and senescence vary within and across species owing to differences in reproductive strategies, but our understanding of the molecular underpinnings of such variation is incomplete. This is perhaps particularly true for sex differences. We investigated the evolution of sex-specific gene expression associated with life history divergence in replicated populations of the seed beetle Acanthoscelides obtectus, experimentally evolving under (E)arly or (L)ate life reproduction for >200 generations which has resulted in strongly divergent life histories. We detected 1,646 genes that were differentially expressed in E and L lines, consistent with a highly polygenic basis of life history evolution. Only 30% of differentially expressed genes were similarly affected in males and females. The evolution of long life was associated with significantly reduced sex differences in expression, especially in non-reproductive tissues. The expression differences were overall more pronounced in females, in accordance with their greater phenotypic divergence in lifespan. Functional enrichment analysis revealed differences between E and L beetles in gene categories previously implicated in aging, such as mitochondrial function and defense response. The results show that divergent life history evolution can be associated with profound changes in gene expression that alter the transcriptome in a sex-specific way, highlighting the importance of understanding the mechanisms of aging in each sex.

Adaptation in Outbred Sexual Yeast is Repeatable, Polygenic and Favors Rare Haplotypes

We carried out a 200 generation Evolve and Resequence (E&R) experiment initiated from an outbred diploid recombined 18-way synthetic base population. Replicate populations were evolved at large effective population sizes (>105 individuals), exposed to several different chemical challenges over 12 weeks of evolution, and whole-genome resequenced. Weekly forced outcrossing resulted in an average between adjacent-gene per cell division recombination rate of ∼0.0008. Despite attempts to force weekly sex, roughly half of our populations evolved cheaters and appear to be evolving asexually. Focusing on seven chemical stressors and 55 total evolved populations that remained sexual we observed large fitness gains and highly repeatable patterns of genome-wide haplotype change within chemical challenges, with limited levels of repeatability across chemical treatments. Adaptation appears highly polygenic with almost the entire genome showing significant and consistent patterns of haplotype change with little evidence for long-range linkage disequilibrium in a subset of populations for which we sequenced haploid clones. That is, almost the entire genome is under selection or drafting with selected sites. At any given locus adaptation was almost always dominated by one of the 18 founder's alleles, with that allele varying spatially and between treatments, suggesting that selection acts primarily on rare variants private to a founder or haplotype blocks harboring multiple mutations.

Natural alleles at the Doa locus underpin evolutionary changes in Drosophila lifespan and fecundity

'Evolve and resequence' (E&R) studies in Drosophila melanogaster have identified many candidate loci underlying the evolution of ageing and life history, but experiments that validate the effects of such candidates remain rare. In a recent E&R study we have identified several alleles of the LAMMER kinase Darkener of apricot (Doa) as candidates for evolutionary changes in lifespan and fecundity. Here, we use two complementary approaches to confirm a functional role of Doa in life-history evolution. First, we used transgenic RNAi to study the effects of Doa at the whole-gene level. Ubiquitous silencing of expression in adult flies reduced both lifespan and fecundity, indicating pleiotropic effects. Second, to characterize segregating variation at Doa, we examined four candidate single nucleotide polymorphisms (SNPs; Doa-1, -2, -3, -4) using a genetic association approach. Three candidate SNPs had effects that were qualitatively consistent with expectations based on our E&R study: Doa-2 pleiotropically affected both lifespan and late-life fecundity; Doa-1 affected lifespan (but not fecundity); and Doa-4 affected late-life fecundity (but not lifespan). Finally, the last candidate allele (Doa-3) also affected lifespan, but in the opposite direction from predicted.

From microbes to mammals: The experimental evolution of aging and longevity across species

Senescence, the functional deterioration of cells or organisms associated with increased age, is pervasive across the tree of life. Yet our understanding of the genetic and physiological basis underlying age-related declines in health and reproduction remains limited. Experimental evolution allows empirical examination of the question of why aging occurs; imposing selection for age-specific fitness traits shifts patterns of aging in experimental populations, enabling investigations of the variation underlying senescence and the mechanisms governing it. Whole-genome sequencing of experimentally evolved populations may reveal candidate genomic variants underlying particular aging patterns; unfortunately, most study systems suffer from limitations that weaken associations between genotypes and phenotypes. In this review, we provide a survey of experimental evolution studies that have altered population-level patterns of reproductive timing and senescence in a variety of species. We discuss the specific selection conditions that have increased longevity, the phenotypic responses and trade-offs that accompany these increases, and examine genomic data collected from these experiments. Additionally, we consider how selected field studies complement laboratory experiments on life-history evolution. Finally, we address the strengths and weaknesses of existing study systems, and evaluate which model organisms appear most promising for future genomic investigations of the evolutionary biology of aging.

Combining Metabolomics and Experimental Evolution Reveals Key Mechanisms Underlying Longevity Differences in Laboratory Evolved Drosophila melanogaster Populations

Experimental evolution with Drosophila melanogaster has been used extensively for decades to study aging and longevity. In recent years, the addition of DNA and RNA sequencing to this framework has allowed researchers to leverage the statistical power inherent to experimental evolution to study the genetic basis of longevity itself. Here, we incorporated metabolomic data into to this framework to generate even deeper insights into the physiological and genetic mechanisms underlying longevity differences in three groups of experimentally evolved D. melanogaster populations with different aging and longevity patterns. Our metabolomic analysis found that aging alters mitochondrial metabolism through increased consumption of NAD+ and increased usage of the TCA cycle. Combining our genomic and metabolomic data produced a list of biologically relevant candidate genes. Among these candidates, we found significant enrichment for genes and pathways associated with neurological development and function, and carbohydrate metabolism. While we do not explicitly find enrichment for aging canonical genes, neurological dysregulation and carbohydrate metabolism are both known to be associated with accelerated aging and reduced longevity. Taken together, our results provide plausible genetic mechanisms for what might be driving longevity differences in this experimental system. More broadly, our findings demonstrate the value of combining multiple types of omic data with experimental evolution when attempting to dissect mechanisms underlying complex and highly polygenic traits such as aging.

Innate Immunity as an Executor of the Programmed Death of Individual Organisms for the Benefit of the Entire Population

In humans, over-activation of innate immunity in response to viral or bacterial infections often causes severe illness and death. Furthermore, similar mechanisms related to innate immunity can cause pathogenesis and death in sepsis, massive trauma (including surgery and burns), ischemia/reperfusion, some toxic lesions, and viral infections including COVID-19. Based on the reviewed observations, we suggest that such severe outcomes may be manifestations of a controlled suicidal strategy protecting the entire population from the spread of pathogens and from dangerous pathologies rather than an aberrant hyperstimulation of defense responses. We argue that innate immunity may be involved in the implementation of an altruistic programmed death of an organism aimed at increasing the well-being of the whole community. We discuss possible ways to suppress this atavistic program by interfering with innate immunity and suggest that combating this program should be a major goal of future medicine.

Evolution, Chance, and Aging

Aging has provided fruitful challenges for evolutionary theory, and evolutionary theory has deepened our understanding of aging. A great deal of genetic and molecular data now exists concerning mortality regulation and there is a growing body of knowledge concerning the life histories of diverse species. Assimilating all relevant data into a framework for the evolution of aging promises to significantly advance the field. We propose extensions of some key concepts to provide greater precision when applying these concepts to age-structured contexts. Secondary or byproduct effects of mutations are proposed as an important factor affecting survival patterns, including effects that may operate in small populations subject to genetic drift, widening the possibilities for mutation accumulation and pleiotropy. Molecular and genetic studies have indicated a diverse array of mechanisms that can modify aging and mortality rates, while transcriptome data indicate a high level of tissue and species specificity for genes affected by aging. The diversity of mechanisms and gene effects that can contribute to the pattern of aging in different organisms may mirror the complex evolutionary processes behind aging.

Fine Mapping without Phenotyping: Identification of Selection Targets in Secondary Evolve and Resequence Experiments

Evolve and Resequence (E&R) studies investigate the genomic selection response of populations in an Experimental Evolution setup. Despite the popularity of E&R, empirical studies in sexually reproducing organisms typically suffer from an excess of candidate loci due to linkage disequilibrium, and single gene or SNP resolution is the exception rather than the rule. Recently, so-called "secondary E&R" has been suggested as promising experimental follow-up procedure to confirm putatively selected regions from a primary E&R study. Secondary E&R provides also the opportunity to increase mapping resolution by allowing for additional recombination events, which separate the selection target from neutral hitchhikers. Here, we use computer simulations to assess the effect of different crossing schemes, population size, experimental duration, and number of replicates on the power and resolution of secondary E&R. We find that the crossing scheme and population size are crucial factors determining power and resolution of secondary E&R: A simple crossing scheme with few founder lines consistently outcompetes crossing schemes where evolved populations from a primary E&R experiment are mixed with a complex ancestral founder population. Regardless of the experimental design tested, a population size of at least 4,800 individuals, which is roughly five times larger than population sizes in typical E&R studies, is required to achieve a power of at least 75%. Our study provides an important step toward improved experimental designs aiming to characterize causative SNPs in Experimental Evolution studies.

The Genomic Architecture of Adaptation to Larval Malnutrition Points to a Trade-off with Adult Starvation Resistance in Drosophila

Periods of nutrient shortage impose strong selection on animal populations. Experimental studies of genetic adaptation to nutrient shortage largely focus on resistance to acute starvation at adult stage; it is not clear how conclusions drawn from these studies extrapolate to other forms of nutritional stress. We studied the genomic signature of adaptation to chronic juvenile malnutrition in six populations of Drosophila melanogaster evolved for 150 generations on an extremely nutrient-poor larval diet. Comparison with control populations evolved on standard food revealed repeatable genomic differentiation between the two set of population, involving >3,000 candidate SNPs forming >100 independently evolving clusters. The candidate genomic regions were enriched in genes implicated in hormone, carbohydrate, and lipid metabolism, including some with known effects on fitness-related life-history traits. Rather than being close to fixation, a substantial fraction of candidate SNPs segregated at intermediate allele frequencies in all malnutrition-adapted populations. This, together with patterns of among-population variation in allele frequencies and estimates of Tajima's D, suggests that the poor diet results in balancing selection on some genomic regions. Our candidate genes for tolerance to larval malnutrition showed a high overlap with genes previously implicated in acute starvation resistance. However, adaptation to larval malnutrition in our study was associated with reduced tolerance to acute adult starvation. Thus, rather than reflecting synergy, the shared genomic architecture appears to mediate an evolutionary trade-off between tolerances to these two forms of nutritional stress.

Meta-analysis suggests the microbiome responds to Evolve and Resequence experiments in Drosophila melanogaster

BACKGROUND

Experimental evolution has a long history of uncovering fundamental insights into evolutionary processes, but has largely neglected one underappreciated component--the microbiome. As eukaryotic hosts evolve, the microbiome may also respond to selection. However, the microbial contribution to host evolution remains poorly understood. Here, we re-analyzed genomic data to characterize the metagenomes from ten Evolve and Resequence (E&R) experiments in Drosophila melanogaster to determine how the microbiome changed in response to host selection.

RESULTS

Bacterial diversity was significantly different in 5/10 studies, primarily in traits associated with metabolism or immunity. Duration of selection did not significantly influence bacterial diversity, highlighting the importance of associations with specific host traits.

CONCLUSIONS

Our genomic re-analysis suggests the microbiome often responds to host selection; thus, the microbiome may contribute to the response of Drosophila in E&R experiments. We outline important considerations for incorporating the microbiome into E&R experiments. The E&R approach may provide critical insights into host-microbiome interactions and fundamental insight into the genomic basis of adaptation.

Convergent Adaptation to Quantitative Host Resistance in a Major Plant Pathogen

Plant pathogens can adapt to quantitative resistance, eroding its effectiveness. The aim of this work was to reveal the genomic basis of adaptation to such a resistance in populations of the fungus Pseudocercospora fijiensis, a major devastating pathogen of banana, by studying convergent adaptation on different cultivars. Samples from P. fijiensis populations showing a local adaptation pattern on new banana hybrids with quantitative resistance were compared, based on a genome scan approach, with samples from traditional and more susceptible cultivars in Cuba and the Dominican Republic. Whole-genome sequencing of pools of P. fijiensis isolates (pool-seq) sampled from three locations per country was conducted according to a paired population design. The findings of different combined analyses highly supported the existence of convergent adaptation on the study cultivars between locations within but not between countries. Five to six genomic regions involved in this adaptation were detected in each country. An annotation analysis and available biological data supported the hypothesis that some genes within the detected genomic regions may play a role in quantitative pathogenicity, including gene regulation. The results suggested that the genetic basis of fungal adaptation to quantitative plant resistance is at least oligogenic, while highlighting the existence of specific host-pathogen interactions for this kind of resistance.IMPORTANCE Understanding the genetic basis of pathogen adaptation to quantitative resistance in plants has a key role to play in establishing durable strategies for resistance deployment. In this context, a population genomic approach was developed for a major plant pathogen (the fungus Pseudocercospora fijiensis causing black leaf streak disease of banana) whereby samples from new resistant banana hybrids were compared with samples from more susceptible conventional cultivars in two countries. A total of 11 genomic regions for which there was strong evidence of selection by quantitative resistance were detected. An annotation analysis and available biological data supported the hypothesis that some of the genes within these regions may play a role in quantitative pathogenicity. These results suggested a polygenic basis of quantitative pathogenicity in this fungal pathogen and complex molecular plant-pathogen interactions in quantitative disease development involving several genes on both sides.

Transposable Element Landscape in Drosophila Populations Selected for Longevity

Transposable elements (TEs) inflict numerous negative effects on health and fitness as they replicate by integrating into new regions of the host genome. Even though organisms employ powerful mechanisms to demobilize TEs, transposons gradually lose repression during aging. The rising TE activity causes genomic instability and was implicated in age-dependent neurodegenerative diseases, inflammation, and the determination of lifespan. It is therefore conceivable that long-lived individuals have improved TE silencing mechanisms resulting in reduced TE expression relative to their shorter-lived counterparts and fewer genomic insertions. Here, we test this hypothesis by performing the first genome-wide analysis of TE insertions and expression in populations of Drosophila melanogaster selected for longevity through late-life reproduction for 50–170 generations from four independent studies. Contrary to our expectation, TE families were generally more abundant in long-lived populations compared with nonselected controls. Although simulations showed that this was not expected under neutrality, we found little evidence for selection driving TE abundance differences. Additional RNA-seq analysis revealed a tendency for reducing TE expression in selected populations, which might be more important for lifespan than regulating genomic insertions. We further find limited evidence of parallel selection on genes related to TE regulation and transposition. However, telomeric TEs were genomically and transcriptionally more abundant in long-lived flies, suggesting improved telomere maintenance as a promising TE-mediated mechanism for prolonging lifespan. Our results provide a novel viewpoint indicating that reproduction at old age increases the opportunity of TEs to be passed on to the next generation with little impact on longevity.

Identifying longevity associated genes by integrating gene expression and curated annotations

Aging is a complex process with poorly understood genetic mechanisms. Recent studies have sought to classify genes as pro-longevity or anti-longevity using a variety of machine learning algorithms. However, it is not clear which types of features are best for optimizing classification performance and which algorithms are best suited to this task. Further, performance assessments based on held-out test data are lacking. We systematically compare five popular classification algorithms using gene ontology and gene expression datasets as features to predict the pro-longevity versus anti-longevity status of genes for two model organisms (C. elegans and S. cerevisiae) using the GenAge database as ground truth. We find that elastic net penalized logistic regression performs particularly well at this task. Using elastic net, we make novel predictions of pro- and anti-longevity genes that are not currently in the GenAge database.

Aging is a complex process with poorly understood genetic mechanisms. Recent studies have sought to classify genes as pro-longevity or anti-longevity using a variety of machine learning algorithms. However, it is not clear which types of features are best for optimizing classification performance and which algorithms are best suited to this task. Further, performance assessments based on held-out test data are lacking. We systematically compare five popular classification algorithms using gene ontology and gene expression datasets as features to predict the pro-longevity versus anti-longevity status of genes for two model organisms (nematode worms and yeast) using the GenAge database as ground truth. We find that elastic net penalized logistic regression performs particularly well at this task. Using elastic net, we make novel predictions of pro- and anti-longevity genes that are not currently in the GenAge database.

Deleterious mutations show increasing negative effects with age in Drosophila melanogaster

BACKGROUND

In order for aging to evolve in response to a declining strength of selection with age, a genetic architecture that allows for mutations with age-specific effects on organismal performance is required. Our understanding of how selective effects of individual mutations are distributed across ages is however poor. Established evolutionary theories assume that mutations causing aging have negative late-life effects, coupled to either positive or neutral effects early in life. New theory now suggests evolution of aging may also result from deleterious mutations with increasing negative effects with age, a possibility that has not yet been empirically explored.

RESULTS

To directly test how the effects of deleterious mutations are distributed across ages, we separately measure age-specific effects on fecundity for each of 20 mutations in Drosophila melanogaster. We find that deleterious mutations in general have a negative effect that increases with age and that the rate of increase depends on how deleterious a mutation is early in life.

CONCLUSIONS

Our findings suggest that aging does not exclusively depend on genetic variants assumed by the established evolutionary theories of aging. Instead, aging can result from deleterious mutations with negative effects that amplify with age. If increasing negative effect with age is a general property of deleterious mutations, the proportion of mutations with the capacity to contribute towards aging may be considerably larger than previously believed.

Detecting selected haplotype blocks in evolve and resequence experiments

Shifting from the analysis of single nucleotide polymorphisms to the reconstruction of selected haplotypes greatly facilitates the interpretation of evolve and resequence (E&R) experiments. Merging highly correlated hitchhiker SNPs into haplotype blocks reduces thousands of candidates to few selected regions. Current methods of haplotype reconstruction from Pool‐seq data need a variety of data‐specific parameters that are typically defined ad hoc and require haplotype sequences for validation. Here, we introduce haplovalidate, a tool which detects selected haplotypes in Pool‐seq time series data without the need for sequenced haplotypes. Haplovalidate makes data‐driven choices of two key parameters for the clustering procedure, the minimum correlation between SNPs constituting a cluster and the window size. Applying haplovalidate to simulated E&R data reliably detects selected haplotype blocks with low false discovery rates. Importantly, our analyses identified a restriction of the haplotype block‐based approach to describe the genomic architecture of adaptation. We detected a substantial fraction of haplotypes containing multiple selection targets. These blocks were considered as one region of selection and therefore led to underestimation of the number of selection targets. We demonstrate that the separate analysis of earlier time points can significantly increase the separation of selection targets into individual haplotype blocks. We conclude that the analysis of selected haplotype blocks has great potential for the characterization of the adaptive architecture with E&R experiments.

Genetic Basis of Increased Lifespan and Postponed Senescence in Drosophila melanogaster

Limited lifespan and senescence are near-universal phenomena. These quantitative traits exhibit variation in natural populations due to the segregation of many interacting loci and from environmental effects. Due to the complexity of the genetic control of lifespan and senescence, our understanding of the genetic basis of variation in these traits is incomplete. Here, we analyzed the pattern of genetic divergence between long-lived (O) Drosophila melanogaster lines selected for postponed reproductive senescence and unselected control (B) lines. We quantified the productivity of the O and B lines and found that reproductive senescence is maternally controlled. We therefore chose 57 candidate genes that are expressed in ovaries, 49 of which have human orthologs, and assessed the effects of RNA interference in ovaries and accessary glands on lifespan and reproduction. All but one candidate gene affected at least one life history trait in one sex or productivity week. In addition, 23 genes had antagonistic pleiotropic effects on lifespan and productivity. Identifying evolutionarily conserved genes affecting increased lifespan and delayed reproductive senescence is the first step toward understanding the evolutionary forces that maintain segregating variation at these loci in nature and may provide potential targets for therapeutic intervention to delay senescence while increasing lifespan.

Context-dependent genetic architecture of Drosophila life span

Understanding the genetic basis of variation in life span is a major challenge that is difficult to address in human populations. Evolutionary theory predicts that alleles affecting natural variation in life span will have properties that enable them to persist in populations at intermediate frequencies, such as late-life-specific deleterious effects, antagonistic pleiotropic effects on early and late-age fitness components, and/or sex- and environment-specific or antagonistic effects. Here, we quantified variation in life span in males and females reared in 3 thermal environments for the sequenced, inbred lines of the Drosophila melanogaster Genetic Reference Panel (DGRP) and an advanced intercross outbred population derived from a subset of DGRP lines. Quantitative genetic analyses of life span and the micro-environmental variance of life span in the DGRP revealed significant genetic variance for both traits within each sex and environment, as well as significant genotype-by-sex interaction (GSI) and genotype-by-environment interaction (GEI). Genome-wide association (GWA) mapping in both populations implicates over 2,000 candidate genes with sex- and environment-specific or antagonistic pleiotropic allelic effects. Over 1,000 of these genes are associated with variation in life span in other D. melanogaster populations. We functionally assessed the effects of 15 candidate genes using RNA interference (RNAi): all affected life span and/or micro-environmental variance of life span in at least one sex and environment and exhibited sex-and environment-specific effects. Our results implicate novel candidate genes affecting life span and suggest that variation for life span may be maintained by variable allelic effects in heterogeneous environments.

Sexual Dimorphisms in Innate Immunity and Responses to Infection in Drosophila melanogaster

The sexes show profound differences in responses to infection and the development of autoimmunity. Dimorphisms in immune responses are ubiquitous across taxa, from arthropods to vertebrates. Drosophila melanogaster shows strong sex dimorphisms in immune system responses at baseline, upon pathogenic challenge, and over aging. We have performed an exhaustive survey of peer-reviewed literature on Drosophila immunity, and present a database of publications indicating the sex(es) analyzed in each study. While we found a growing interest in the community in adult immunity and in reporting both sexes, the main body of work in this field uses only one sex, or does not stratify by sex. We synthesize evidence for sexually dimorphic responses to bacterial, viral, and fungal infections. Dimorphisms may be mediated by distinct immune compartments, and we review work on sex differences in behavioral, epithelial, cellular, and systemic (fat body-mediated) immunity. Emerging work on sexually dimorphic aging of immune tissues, immune senescence, and inflammation are examined. We consider evolutionary drivers for sex differences in immune investment, highlight the features of Drosophila biology that make it particularly amenable to studies of immune dimorphisms, and discuss areas for future exploration.

Diverse biological processes coordinate the transcriptional response to nutritional changes in a Drosophila melanogaster multiparent population

Background

Environmental variation in the amount of resources available to populations challenge individuals to optimize the allocation of those resources to key fitness functions. This coordination of resource allocation relative to resource availability is commonly attributed to key nutrient sensing gene pathways in laboratory model organisms, chiefly the insulin/TOR signaling pathway. However, the genetic basis of diet-induced variation in gene expression is less clear.

Results

To describe the natural genetic variation underlying nutrient-dependent differences, we used an outbred panel derived from a multiparental population, the Drosophila Synthetic Population Resource. We analyzed RNA sequence data from multiple female tissue samples dissected from flies reared in three nutritional conditions: high sugar (HS), dietary restriction (DR), and control (C) diets. A large proportion of genes in the experiment (19.6% or 2471 genes) were significantly differentially expressed for the effect of diet, and 7.8% (978 genes) for the effect of the interaction between diet and tissue type (LRT, P_adj. < 0.05). Interestingly, we observed similar patterns of gene expression relative to the C diet, in the DR and HS treated flies, a response likely reflecting diet component ratios. Hierarchical clustering identified 21 robust gene modules showing intra-modularly similar patterns of expression across diets, all of which were highly significant for diet or diet-tissue interaction effects (FDR P_adj. < 0.05). Gene set enrichment analysis for different diet-tissue combinations revealed a diverse set of pathways and gene ontology (GO) terms (two-sample t-test, FDR < 0.05). GO analysis on individual co-expressed modules likewise showed a large number of terms encompassing many cellular and nuclear processes (Fisher exact test, P_adj. < 0.01). Although a handful of genes in the IIS/TOR pathway including Ilp5, Rheb, and Sirt2 showed significant elevation in expression, many key genes such as InR, chico, most insulin peptide genes, and the nutrient-sensing pathways were not observed.

Conclusions

Our results suggest that a more diverse network of pathways and gene networks mediate the diet response in our population. These results have important implications for future studies focusing on diet responses in natural populations.

Life-History Evolution and the Genetics of Fitness Components in Drosophila melanogaster

Life-history traits or "fitness components"-such as age and size at maturity, fecundity and fertility, age-specific rates of survival, and life span-are the major phenotypic determinants of Darwinian fitness. Analyzing the evolution and genetics of these phenotypic targets of selection is central to our understanding of adaptation. Due to its simple and rapid life cycle, cosmopolitan distribution, ease of maintenance in the laboratory, well-understood evolutionary genetics, and its versatile genetic toolbox, the "vinegar fly" Drosophila melanogaster is one of the most powerful, experimentally tractable model systems for studying "life-history evolution." Here, I review what has been learned about the evolution and genetics of life-history variation in D. melanogaster by drawing on numerous sources spanning population and quantitative genetics, genomics, experimental evolution, evolutionary ecology, and physiology. This body of work has contributed greatly to our knowledge of several fundamental problems in evolutionary biology, including the amount and maintenance of genetic variation, the evolution of body size, clines and climate adaptation, the evolution of senescence, phenotypic plasticity, the nature of life-history trade-offs, and so forth. While major progress has been made, important facets of these and other questions remain open, and the D. melanogaster system will undoubtedly continue to deliver key insights into central issues of life-history evolution and the genetics of adaptation.

Distinct genomic signals of lifespan and life history evolution in response to postponed reproduction and larval diet in Drosophila

Reproduction and diet are two major factors controlling the physiology of aging and life history, but how they interact to affect the evolution of longevity is unknown. Moreover, although studies of large-effect mutants suggest an important role of nutrient sensing pathways in regulating aging, the genetic basis of evolutionary changes in lifespan remains poorly understood. To address these questions, we analyzed the genomes of experimentally evolved Drosophila melanogaster populations subjected to a factorial combination of two selection regimes: reproductive age (early versus postponed), and diet during the larval stage ("low," "control," "high"), resulting in six treatment combinations with four replicate populations each. Selection on reproductive age consistently affected lifespan, with flies from the postponed reproduction regime having evolved a longer lifespan. In contrast, larval diet affected lifespan only in early-reproducing populations: flies adapted to the "low" diet lived longer than those adapted to control diet. Here, we find genomic evidence for strong independent evolutionary responses to either selection regime, as well as loci that diverged in response to both regimes, thus representing genomic interactions between the two. Overall, we find that the genomic basis of longevity is largely independent of dietary adaptation. Differentiated loci were not enriched for "canonical" longevity genes, suggesting that naturally occurring genic targets of selection for longevity differ qualitatively from variants found in mutant screens. Comparing our candidate loci to those from other "evolve and resequence" studies of longevity demonstrated significant overlap among independent experiments. This suggests that the evolution of longevity, despite its presumed complex and polygenic nature, might be to some extent convergent and predictable.

Multiple genetic loci affect place learning and memory performance in Drosophila melanogaster

Learning and memory are critical functions for all animals, giving individuals the ability to respond to changes in their environment. Within populations, individuals vary, however the mechanisms underlying this variation in performance are largely unknown. Thus, it remains to be determined what genetic factors cause an individual to have high learning ability and what factors determine how well an individual will remember what they have learned. To genetically dissect learning and memory performance, we used the Drosophila synthetic population resource (DSPR), a multiparent mapping resource in the model system Drosophila melanogaster, consisting of a large set of recombinant inbred lines (RILs) that naturally vary in these and other traits. Fruit flies can be trained in a "heat box" to learn to remain on one side of a chamber (place learning) and can remember this (place memory) over short timescales. Using this paradigm, we measured place learning and memory for ~49 000 individual flies from over 700 DSPR RILs. We identified 16 different loci across the genome that significantly affect place learning and/or memory performance, with 5 of these loci affecting both traits. To identify transcriptomic differences associated with performance, we performed RNA-Seq on pooled samples of seven high performing and seven low performing RILs for both learning and memory and identified hundreds of genes with differences in expression in the two sets. Integrating our transcriptomic results with the mapping results allowed us to identify nine promising candidate genes, advancing our understanding of the genetic basis underlying natural variation in learning and memory performance.

A clinal polymorphism in the insulin signaling transcription factor foxo contributes to life-history adaptation in Drosophila

A fundamental aim of adaptation genomics is to identify polymorphisms that underpin variation in fitness traits. In Drosophila melanogaster, latitudinal life-history clines exist on multiple continents and make an excellent system for dissecting the genetics of adaptation. We have previously identified numerous clinal single-nucleotide polymorphism in insulin/insulin-like growth factor signaling (IIS), a pathway known from mutant studies to affect life history. However, the effects of natural variants in this pathway remain poorly understood. Here we investigate how two clinal alternative alleles at foxo, a transcriptional effector of IIS, affect fitness components (viability, size, starvation resistance, fat content). We assessed this polymorphism from the North American cline by reconstituting outbred populations, fixed for either the low- or high-latitude allele, from inbred DGRP lines. Because diet and temperature modulate IIS, we phenotyped alleles across two temperatures (18°C, 25°C) and two diets differing in sugar source and content. Consistent with clinal expectations, the high-latitude allele conferred larger body size and reduced wing loading. Alleles also differed in starvation resistance and expression of insulin-like receptor, a transcriptional target of FOXO. Allelic reaction norms were mostly parallel, with few GxE interactions. Together, our results suggest that variation in IIS makes a major contribution to clinal life-history adaptation.

Slow life-history strategies are associated with negligible actuarial senescence in western Palaearctic salamanders

Actuarial senescence has been viewed for a long time as an inevitable and uniform process. However, the work on senescence has mainly focused on endotherms with deterministic growth and low regeneration capacity during the adult stage, leading to a strong taxonomic bias in the study of ageing. Recent studies have highlighted that senescence could indeed display highly variable trajectories that correlate with species life-history traits. Slow life histories and indeterminate growth seem to be associated with weak and late senescence. Furthermore, high regenerative abilities could lead to negligible senescence in ectotherms. However, demographic data for species that would allow testing of these hypotheses are scarce. Here, we investigated senescence patterns in 'true salamanders' from the western Palaearctic. Our results showed that salamanders have slow life histories and that they experience negligible senescence. This pattern was consistent at both intra- and interspecific levels, suggesting that the absence of senescence may be a phylogenetically conserved trait. The regenerative capacities of salamanders, in combination with other physiological and developmental features such as an indeterminate growth and a low metabolic rate, probably explain why these small ectotherms have lifespans similar to that of large endotherms and, in contrast with most amniotes, undergo negligible senescence. Our study seriously challenges the idea that senescence is a ubiquitous phenomenon in the tree of life.

Drosophila transcriptomics with and without ageing

The genomic basis of ageing still remains unknown despite being a topic of study for many years. Here, we present data from 20 experimentally evolved laboratory populations of Drosophila melanogaster that have undergone two different life-history selection regimes. One set of ten populations demonstrates early ageing whereas the other set of ten populations shows postponed ageing. Additionally, both types of populations consist of five long standing populations and five recently derived populations. Our primary goal was to determine which genes exhibit changes in expression levels by comparing the female transcriptome of the two population sets at two different time points. Using three different sets of increasingly restrictive criteria, we found that 2.1-15.7% (82-629 genes) of the expressed genes are associated with differential ageing between population sets. Conversely, a comparison of recently derived populations to long-standing populations reveals little to no transcriptome differentiation, suggesting that the recent selection regime has had a larger impact on the transcriptome than its more distant evolutionary history. In addition, we found very little evidence for significant enrichment for functional attributes regardless of the set of criteria used. Relative to previous ageing studies, we find little overlap with other lists of aging related genes. The disparity between our results and previously published results is likely due to the high replication used in this study coupled with our use of highly differentiated populations. Our results reinforce the notion that the use of genomic, transcriptomic, and phenotypic data to uncover the genetic basis of a complex trait like ageing can benefit from experimental designs that use highly replicated, experimentally-evolved populations.

Plastic Senescence in the Honey Bee and the Disposable Soma Theory

The demonstration of life span plasticity in natural populations would provide a powerful test of evolutionary theories of senescence. Plastic senescence is not easily explained by mutation accumulation or antagonistic pleiotropy but is a corollary of the disposable soma theory. The life span differences among castes of the eusocial Hymenoptera are potentially some of the most striking and extreme examples of life span plasticity. Although these differences are often assumed to be plastic, this has never been demonstrated conclusively because differences in life span may be caused by the proximate effects of different levels of environmental hazard experienced by castes. Here age-dependent and age-independent components of instantaneous mortality rates of the honey bee (Apis mellifera) were estimated from published life tables for natural and seminatural populations to determine whether differences in life span between queens and workers and between different types of workers are indeed plastic. These differences in life span were found to be due to differences in the rate of actuarial senescence, which correlate positively with the rate of extrinsic mortality, in accordance with the central prediction of evolutionary theories of senescence. Although all three evolutionary theories of senescence could in principle explain such plastic senescence, given differential gene expression between castes or life stages, only the disposable soma theory adequately explains the adaptive regulation of somatic maintenance in response to different environmental conditions that appears to underlie life span plasticity.

Natural and Artificial Selection for Parasitoid Resistance in Drosophila melanogaster Leave Different Genetic Signatures

Adaptation of complex traits depends on standing genetic variation at multiple loci. The allelic variants that have positive fitness effects, however, can differ depending on the genetic background and the selective pressure. Previously, we interrogated the Drosophila melanogaster genome at the population level for polymorphic positions and identified 215 single nucleotide polymorphisms (SNPs) that had significantly changed in frequency after experimental evolution for increased parasitoid resistance. In the current study, we follow up on 11 of these SNPs as putative targets of the experimental selection process (Jalvingh et al., 2014). We study the patterns of genetic variation for these SNPs in several European field populations. Furthermore, we associate the genetic variation of these SNPs to variation in resistance against the parasitoid Asobara tabida, by determining the individual phenotype and SNP genotype for 144 individuals from four Selection lines and four non-selected Control lines and for 400 individuals from 12 Field lines that differ in parasitoid resistance. For the Selection lines we additionally monitored the changes in allele frequencies throughout the five generations of experimental selection. For three genes, mbl (Zn-finger protein), mthl4 (G-protein coupled receptor) and CG17287 (protein-cysteine S-palmitoyltransferase) individual SNP genotypes were significantly associated with resistance level in the Selection and Control lines. Additionally, the minor allele in mbl and mthl4 were consistently and gradually favored throughout the five generations of experimental evolution. However, none of these alleles did appear to be associated to high resistance in the Field lines. We suggest that, within field populations, selection for parasitoid resistance is a gradual process that involves co-adapted gene complexes. Fast artificial selection, however, enforces the sudden cumulating of particular alleles that confer high resistance (genetic sweep). We discuss our findings in the context of local adaptation.

The evolutionary potential of diet-dependent effects on lifespan and fecundity in a multi-parental population of Drosophila melanogaster

The nutritional conditions experienced by a population have a major role in shaping trait evolution in many taxa. Constraints exerted by nutrient limitation or nutrient imbalance can influence the maximal value that fitness components such as reproduction and lifespan attains, and organisms may shift how resources are allocated to different structures and functions in response to changes in nutrition. Whether the phenotypic changes associated with changes in nutrition represent an adaptive response is largely unknown. Further, it is unclear whether the response of fitness components to diet even has the potential to evolve in most systems. In this study, we use an admixed multi-parental population of Drosophila melanogaster reared in three different diet conditions to estimate quantitative genetic parameters for lifespan and fecundity. We find significant genetic variation for both traits in our population and show that lifespan has moderate to high heritabilities within diets. Genetic correlations for lifespan between diets were significantly less than one, demonstrating a strong genotype by diet interaction. These findings demonstrate substantial standing genetic variation in our population that is comparable to natural populations and highlights the potential for adaptation to changing nutritional environments.

Mate choice and gene expression signatures associated with nutritional adaptation in the medfly (Ceratitis capitata)

Evolutionary responses to nutrition are key to understanding host shifts and the resulting potential for reproductive isolation. Experimental evolution has previously been used to describe the responses of the medfly (Ceratitis capitata) to larval diets with different nutritional properties. Within 30 generations this led to divergence in larval development time, egg to adult survival and adaptation in adult body size. Here we used mRNA-seq to identify differences in gene expression patterns in these same populations, using males from the 60^th generation of nutritional selection. We validated differential expression by using qRT-PCR and found that genes linked to metabolism, oxidative phosphorylation and proteolysis were significantly over-represented among the differentially expressed genes. The results provide the first genome-wide survey of the putative mechanisms underpinning evolved responses to nutritional adaptation. In addition, we tested the hypothesis that nutritional adaptation can alter mating patterns. We found evidence for assortative mating by diet at generation 60, but not 90. Hence, the pattern was variable across generations and there was no evidence overall for any isolating mating divergence between the lines. Overall, the results provide insight into the mechanisms underpinning dietary adaptation and extend our knowledge of which traits represent core responses to nutritional selection.

Pervasive Linked Selection and Intermediate-Frequency Alleles Are Implicated in an Evolve-and-Resequencing Experiment of Drosophila simulans

We develop analytical and simulation tools for evolve-and-resequencing experiments and apply them to a new study of rapid evolution in Drosophila simulans Likelihood test statistics applied to pooled population sequencing data suggest parallel evolution of 138 SNPs across the genome. This number is reduced by orders of magnitude from previous studies (thousands or tens of thousands), owing to differences in both experimental design and statistical analysis. Whole genome simulations calibrated from Drosophila genetic data sets indicate that major features of the genome-wide response could be explained by as few as 30 loci under strong directional selection with a corresponding hitchhiking effect. Smaller effect loci are likely also responding, but are below the detection limit of the experiment. Finally, SNPs showing strong parallel evolution in the experiment are intermediate in frequency in the natural population (usually 30-70%) indicative of balancing selection in nature. These loci also exhibit elevated differentiation among natural populations of D. simulans, suggesting environmental heterogeneity as a potential balancing mechanism.

Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis

Background

Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes.

Results

Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes.

Conclusions

Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5310-9) contains supplementary material, which is available to authorized users.

Evolution of longevity improves immunity in Drosophila

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.

Horizons in the evolution of aging

Between the 1930s and 50s, evolutionary biologists developed a successful theory of why organisms age, firmly rooted in population genetic principles. By the 1980s the evolution of aging had a secure experimental basis. Since the force of selection declines with age, aging evolves due to mutation accumulation or a benefit to fitness early in life. Here we review major insights and challenges that have emerged over the last 35 years: selection does not always necessarily decline with age; higher extrinsic (i.e., environmentally caused) mortality does not always accelerate aging; conserved pathways control aging rate; senescence patterns are more diverse than previously thought; aging is not universal; trade-offs involving lifespan can be 'broken'; aging might be 'druggable'; and human life expectancy continues to rise but compressing late-life morbidity remains a pressing challenge.

MimicrEE2: Genome-wide forward simulations of Evolve and Resequencing studies

Evolve and Resequencing (E&R) studies allow us to monitor adaptation at the genomic level. By sequencing evolving populations at regular time intervals, E&R studies promise to shed light on some of the major open questions in evolutionary biology such as the repeatability of evolution and the molecular basis of adaptation. However, data interpretation, statistical analysis and the experimental design of E&R studies increasingly require simulations of evolving populations, a task that is difficult to accomplish with existing tools, which may i) be too slow, ii) require substantial reformatting of data, iii) not support an adaptive scenario of interest or iv) not sufficiently capture the biology of the used model organism. Therefore we developed MimicrEE2, a multi-threaded Java program for genome-wide forward simulations of evolving populations. MimicrEE2 enables the convenient usage of available genomic resources, supports biological particulars of model organism frequently used in E&R studies and offers a wide range of different adaptive models (selective sweeps, polygenic adaptation, epistasis). Due to its user-friendly and efficient design MimicrEE2 will facilitate simulations of E&R studies even for small labs with limited bioinformatics expertise or computational resources. Additionally, the scripts provided for executing MimicrEE2 on a computer cluster permit the coverage even of a large parameter space. MimicrEE2 runs on any computer with Java installed. It is distributed under the GPLv3 license at https://sourceforge.net/projects/mimicree2/.

The interplay between immunity and aging in Drosophila

Here, we provide a brief review of the mechanistic connections between immunity and aging—a fundamental biological relationship that remains poorly understood—by considering two intertwined questions: how does aging affect immunity, and how does immunity affect aging? On the one hand, aging contributes to the deterioration of immune function and predisposes the organism to infections (“immuno-senescence”). On the other hand, excessive activation of the immune system can accelerate degenerative processes, cause inflammation and immunopathology, and thus promote aging (“inflammaging”). Interestingly, several recent lines of evidence support the hypothesis that restrained or curbed immune activity at old age (that is, optimized age-dependent immune homeostasis) might actually improve realized immune function and thereby promote longevity. We focus mainly on insights from Drosophila, a powerful genetic model system in which both immunity and aging have been extensively studied, and conclude by outlining several unresolved questions in the field.

Genome-Wide Analysis of Starvation-Selected Drosophila melanogaster-A Genetic Model of Obesity

Experimental evolution affords the opportunity to investigate adaptation to stressful environments. Studies combining experimental evolution with whole-genome resequencing have provided insight into the dynamics of adaptation and a new tool to uncover genes associated with polygenic traits. Here, we selected for starvation resistance in populations of Drosophila melanogaster for over 80 generations. In response, the starvation-selected lines developed an obese condition, storing nearly twice the level of total lipids than their unselected controls. Although these fats provide a ∼3-fold increase in starvation resistance, the imbalance in lipid homeostasis incurs evolutionary cost. Some of these tradeoffs resemble obesity-associated pathologies in mammals including metabolic depression, low activity levels, dilated cardiomyopathy, and disrupted sleeping patterns. To determine the genetic basis of these traits, we resequenced genomic DNA from the selected lines and their controls. We found 1,046,373 polymorphic sites, many of which diverged between selection treatments. In addition, we found a wide range of genetic heterogeneity between the replicates of the selected lines, suggesting multiple mechanisms of adaptation. Genome-wide heterozygosity was low in the selected populations, with many large blocks of SNPs nearing fixation. We found candidate loci under selection by using an algorithm to control for the effects of genetic drift. These loci were mapped to a set of 382 genes, which associated with many processes including nutrient response, catabolic metabolism, and lipid droplet function. The results of our study speak to the evolutionary origins of obesity and provide new targets to understand the polygenic nature of obesity in a unique model system.

Genetic Dissection of Nutrition-Induced Plasticity in Insulin/Insulin-Like Growth Factor Signaling and Median Life Span in a Drosophila Multiparent Population

The nutritional environments that organisms experience are inherently variable, requiring tight coordination of how resources are allocated to different functions relative to the total amount of resources available. A growing body of evidence supports the hypothesis that key endocrine pathways play a fundamental role in this coordination. In particular, the insulin/insulin-like growth factor signaling (IIS) and target of rapamycin (TOR) pathways have been implicated in nutrition-dependent changes in metabolism and nutrient allocation. However, little is known about the genetic basis of standing variation in IIS/TOR or how diet-dependent changes in expression in this pathway influence phenotypes related to resource allocation. To characterize natural genetic variation in the IIS/TOR pathway, we used >250 recombinant inbred lines (RILs) derived from a multiparental mapping population, the Drosophila Synthetic Population Resource, to map transcript-level QTL of genes encoding 52 core IIS/TOR components in three different nutritional environments [dietary restriction (DR), control (C), and high sugar (HS)]. Nearly all genes, 87%, were significantly differentially expressed between diets, though not always in ways predicted by loss-of-function mutants. We identified cis (i.e., local) expression QTL (eQTL) for six genes, all of which are significant in multiple nutrient environments. Further, we identified trans (i.e., distant) eQTL for two genes, specific to a single nutrient environment. Our results are consistent with many small changes in the IIS/TOR pathways. A discriminant function analysis for the C and DR treatments identified a pattern of gene expression associated with the diet treatment. Mapping the composite discriminant function scores revealed a significant global eQTL within the DR diet. A correlation between the discriminant function scores and the median life span (r = 0.46) provides evidence that gene expression changes in response to diet are associated with longevity in these RILs.

Drosophila simulans: A Species with Improved Resolution in Evolve and Resequence Studies

The combination of experimental evolution with high-throughput sequencing of pooled individuals—i.e., evolve and resequence (E&R)—is a powerful approach to study adaptation from standing genetic variation under controlled, replicated conditions. Nevertheless, E&R studies in Drosophila melanogaster have frequently resulted in inordinate numbers of candidate SNPs, particularly for complex traits. Here, we contrast the genomic signature of adaptation following ∼60 generations in a novel hot environment for D. melanogaster and D. simulans. For D. simulans, the regions carrying putatively selected loci were far more distinct, and thus harbored fewer false positives, than those in D. melanogaster. We propose that species without segregating inversions and higher recombination rates, such as D. simulans, are better suited for E&R studies that aim to characterize the genetic variants underlying the adaptive response.

How to get the most bang for your buck: the evolution and physiology of nutrition-dependent resource allocation strategies

All organisms use resources to grow, survive and reproduce. The supply of these resources varies widely across landscapes and time, imposing ultimate constraints on the maximal trait values for allocation-related traits. In this review, we address three key questions fundamental to our understanding of the evolution of allocation strategies and their underlying mechanisms. First, we ask: how diverse are flexible resource allocation strategies among different organisms? We find there are many, varied, examples of flexible strategies that depend on nutrition. However, this diversity is often ignored in some of the best-known cases of resource allocation shifts, such as the commonly observed pattern of lifespan extension under nutrient limitation. A greater appreciation of the wide variety of flexible allocation strategies leads directly to our second major question: what conditions select for different plastic allocation strategies? Here, we highlight the need for additional models that explicitly consider the evolution of phenotypically plastic allocation strategies and empirical tests of the predictions of those models in natural populations. Finally, we consider the question: what are the underlying mechanisms determining resource allocation strategies? Although evolutionary biologists assume differential allocation of resources is a major factor limiting trait evolution, few proximate mechanisms are known that specifically support the model. We argue that an integrated framework can reconcile evolutionary models with proximate mechanisms that appear at first glance to be in conflict with these models. Overall, we encourage future studies to: (i) mimic ecological conditions in which those patterns evolve, and (ii) take advantage of the 'omic' opportunities to produce multi-level data and analytical models that effectively integrate across physiological and evolutionary theory.

Clear: Composition of Likelihoods for Evolve and Resequence Experiments

The advent of next generation sequencing technologies has made whole-genome and whole-population sampling possible, even for eukaryotes with large genomes. With this development, experimental evolution studies can be designed to observe molecular evolution "in action" via evolve-and-resequence (E&R) experiments. Among other applications, E&R studies can be used to locate the genes and variants responsible for genetic adaptation. Most existing literature on time-series data analysis often assumes large population size, accurate allele frequency estimates, or wide time spans. These assumptions do not hold in many E&R studies. In this article, we propose a method-composition of likelihoods for evolve-and-resequence experiments (Clear)-to identify signatures of selection in small population E&R experiments. Clear takes whole-genome sequences of pools of individuals as input, and properly addresses heterogeneous ascertainment bias resulting from uneven coverage. Clear also provides unbiased estimates of model parameters, including population size, selection strength, and dominance, while being computationally efficient. Extensive simulations show that Clear achieves higher power in detecting and localizing selection over a wide range of parameters, and is robust to variation of coverage. We applied the Clear statistic to multiple E&R experiments, including data from a study of adaptation of Drosophila melanogaster to alternating temperatures and a study of outcrossing yeast populations, and identified multiple regions under selection with genome-wide significance.