Laboratory selection experiments using Drosophila: what do they really tell us?

An analysis of direct and indirect effects in Drosophila melanogaster undergoing a few cycles of experimental evolution for stress-related traits

The physiological mechanisms underpinning adaptations to starvation and cold stresses have been extensively studied in Drosophila, yet the understanding of correlated changes in stress-related and life-history traits, as well as the energetics of stress tolerance, still remains elusive. To answer the questions empirically in this context, we allowed D. melanogaster to evolve for either increased starvation or cold tolerance (24-generations / regime) in an experimental evolution system, and examined whether selection of either trait affects un-selected stress trait, as well as the impacts potential changes in life-history and mating success-related traits. Our results revealed remarkable changes in starvation/cold tolerance (up to 1.5-fold) as a direct effect of selection, while cold tolerance had been dramatically reduced (1.26-fold) in the starvation tolerant (ST) lines compared to control counterparts, although no such changes were evident in cold-tolerant (CT) lines. ST lines exhibited a higher level of body lipids and a reduced level of trehalose content, while CT lines accumulated a greater levels of body lipid and trehalose contents. Noticeably, we found that selection for starvation or cold tolerance positively correlates with larval development time, longevity, and copulation duration, indicating that these traits are among the most common targets of selection trajectories shaping stress tolerance. Altogether, this study highlights the complexity of mechanisms evolved in ST lines that contribute to enhanced starvation tolerance, but also negatively impact cold tolerance. Nevertheless, mechanisms foraging enhanced cold tolerance in CT lines appear not to target starvation tolerance. Moreover, the parallel changes in life history/mating success traits across stress regimes could indicate some generic pathways evolved in stressful environments, targeting life-history and mating success characteristics to optimize fitness.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

Experimental test of the contributions of initial variation and new mutations to adaptive evolution in a novel environment

Experimental evolution is an approach that allows researchers to study organisms as they evolve in controlled environments. Despite the growing popularity of this approach, there are conceptual gaps among projects that use different experimental designs. One such gap concerns the contributions to adaptation of genetic variation present at the start of an experiment and that of new mutations that arise during an experiment. The primary source of genetic variation has historically depended largely on the study organisms. In the long-term evolution experiment (LTEE) using Escherichia coli, for example, each population started from a single haploid cell, and therefore, adaptation depended entirely on new mutations. Most other microbial evolution experiments have followed the same strategy. By contrast, evolution experiments using multicellular, sexually reproducing organisms typically start with preexisting variation that fuels the response to selection. New mutations may also come into play in later generations of these experiments, but it is generally difficult to quantify their contribution in these studies. Here, we performed an experiment using E. coli to compare the contributions of initial genetic variation and new mutations to adaptation in a new environment. Our experiment had four treatments that varied in their starting diversity, with 18 populations in each treatment. One treatment depended entirely on new mutations, while the other three began with mixtures of clones, whole-population samples, or mixtures of whole-population samples from the LTEE. We tracked a genetic marker associated with different founders in two treatments. These data revealed significant variation in fitness among the founders, and that variation impacted evolution in the early generations of our experiment. However, there were no differences in fitness among the treatments after 500 or 2,000 generations in the new environment, despite the variation in fitness among the founders. These results indicate that new mutations quickly dominated, and eventually they contributed more to adaptation than did the initial variation. Our study thus shows that preexisting genetic variation can have a strong impact on early evolution in a new environment, but new beneficial mutations may contribute more to later evolution and can even drive some initially beneficial variants to extinction.

Urban evolution of thermal physiology in a range-expanding, mycophagous fruit fly, Drosophila tripunctata

In Drosophila spp., their often high number of annual generations, large population sizes and large amounts of standing genetic variation should predispose them to undergo contemporary adaptation to climatic warming. However, a number of laboratory experimental evolution studies in this group of organisms suggest strong limits on the rate and magnitude of contemporary thermal adaptation. Here, we explore this discrepancy by examining the potential for rapid evolutionary divergence between wild populations of Drosophila tripunctata Loew, 1862 from rural and urban sites. We performed a multi-generation common garden study and found evidence for the evolution of higher heat tolerance (critical thermal maximum) in flies from urban populations. We also detected evolutionary divergence in cold resistance (chill coma recovery time), with diminished cold resistance in flies from urban populations, although the effect was weaker than the shift in heat tolerance. Our study provides evidence of contemporary urban thermal adaptation, although the magnitude of phenotypic change lagged the magnitude of environmental temperature change across the urbanization gradient, suggesting potential limits on the evolution of urban thermal physiology.

Disentangling the effects of jasmonate and tissue loss on the sex allocation of an annual plant

Selection through pollinators plays a major role in the evolution of reproductive traits. However, herbivory can also induce changes in plant sexual expression and sexual systems, potentially influencing conditions governing transitions between sexual systems. Previous work has shown that herbivory has a strong effect on sex allocation in the wind-pollinated annual plant Mercurialis annua, likely via responses to resource loss. It is also known that many plants respond to herbivory by inducing signaling, and endogenous responses to it, via the plant hormone jasmonate. Here, we attempt to uncouple the effects of herbivory on sex allocation in M. annua through resource limitation (tissue loss) versus plant responses to jasmonate hormone signaling. We used a two-factorial experiment with four treatment combinations: control, herbivory (25% chronic tissue loss), jasmonate, and combined herbivory and jasmonate. We estimated the effects of tissue loss and defense-inducing hormones on reproductive allocation, male reproductive effort, and sex allocation. Tissue loss caused plants to reduce their male reproductive effort, resulting in changes in total sex allocation. However, application of jasmonate after herbivory reversed its effect on male investment. Our results show that herbivory has consequences on plant sex expression and sex allocation, and that defense-related hormones such as jasmonate can buffer the impacts. We discuss the physiological mechanisms that might underpin the effects of herbivory on sex allocation, and their potential implications for the evolution of plant sexual systems.

Linking xylem structure and function: the comparative method in from the cold

This article is a Commentary on Savage et al. (2022), 235: 953–964.

Molecular variation of the cytochrome b DNA and protein sequences in Phytoseiulus macropilis and P. persimilis (Acari: Phytoseiidae) reflect population differentiation

Several phytoseiid mite species are important natural enemies used in biological control strategies. In the present study, Cytb mtDNA sequences of various populations of two species, Phytoseiulus macropolis and P. persimilis, were compared to determine whether the specimens collected in Brazil could belong to P. persimilis as this latter species is reported in South America but not in Brazil. The Cytb marker was used because of its high evolution rate, assumed to capture intraspecific variation. No overlap between intra- and interspecific distances was observed but the distances were quite low for interspecific variation. This can be due to the particular biology of Phytoseiulus species and this shows the difficulty to apply a universal threshold in genetic distances to conclude about the existence of one or several species. Cytb mtDNA sequences were also considered to assess intraspecific variation. The DNA sequences of P. persimilis populations were very similar, probably because they all originated from the West Palearctic region or because of a prevalence of commercialized specimens in natura. For P. macropilis, higher genetic distances were observed and differentiation was noted according to geographic location and, to a smaller extent, pyrethroid resistance. To determine how DNA variation might impact the protein function (CytB fragment considered), the amino acid compositions of the populations studied were compared. No diagnostic mutation was observed between pyrethroid resistant and susceptible populations, whereas four mutations were identified between populations of P. macropilis separated by 1300 km (different climatic conditions). The impact of such mutations is discussed but knowledge is scarce, which makes it difficult to root testable hypotheses. The protein analysis clearly opens new perspectives in Phytoseiidae studies.

How well do genetic markers inform about responses to intraspecific admixture? A comparative analysis of microsatellites and RADseq

Background

Fitness consequences of intraspecific genetic admixture can vary from positive to negative depending on the genetic composition of the populations and environmental conditions. Because admixture has potential to influence the success of management and conservation efforts, genetic similarity has been suggested to be used as a proxy to predict the outcome. Studies utilizing microsatellites (a neutral marker) to investigate associations between genetic distance and admixture effects show conflicting results. Marker types that yield information on genome-wide and/or adaptive variation might be more useful for predicting responses to inter-population hybridization. In this study we utilized published data for three populations of pike (Esox lucius) to investigate associations between offspring performance (hatching success) and parental genetic similarity in experimentally purebred and admixed families, based on neutral (microsatellites), genome-wide neutral (RADseq SNPs), and adaptive (SNPs under selection) markers.

Results

Estimated similarity varied among the markers, likely reflecting differences in their inherent properties, but was consistently higher in purebred than admixed families. A significant interaction between marker type and admixture treatment reflected that neutral SNPs yielded higher estimates than adaptive SNPs for admixed families whereas no difference was found for purebred families, which indicates that neutral similarity was not reflective of adaptive similarity. When all samples were pooled, no association between similarity and performance was found for any marker. For microsatellites, similarity was positively correlated with hatching success in purebred families, whereas no association was found in admixed families; however, the direction of the effect differed between the population combinations.

Conclusions

The results strengthen the notion that, as of today, there is no proxy that can reliably predicted the outcome of admixture. This emphasizes the need of further studies to advance knowledge that can shed light on how to safeguard against negative consequences of admixture, and thereby inform management and promote conservation of biological diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00974-3.

Spatial heterogeneity in resources alters selective dynamics in Drosophila melanogaster

Environmental features can alter the behaviors and phenotypes of organisms, influencing the dynamics of natural and sexual selection. Experimental environmental manipulation, particularly when conducted in experiments where the dynamics of the purging of deleterious alleles are compared, has demonstrated both direct and indirect effects on the strength and direction of selection. However, many of these studies are conducted with fairly simplistic environments, where it is not always clear how or why particular forms of spatial heterogeneity influence behavior or selection. Using Drosophila melanogaster, we tested three different spatial environments designed to determine if spatial constraint of critical resources influences the efficiency of natural and sexual selection. We conducted two allele purging experiments to (1) assess effects of these spatial treatments on selective dynamics of six recessive mutations, and (2) determine how these dynamics changed when sexual selection was relaxed and spatial area reduced for two of the mutants. Allele purging dynamics depended on spatial environment, however the patterns of purging rates between the environments differed across distinct deleterious mutations. We also tested two of the mutant alleles, and demonstrate sexual selection increased the purging rate.

Competitive history shapes rapid evolution in a seasonal climate

Eco-evolutionary dynamics will play a critical role in determining species' fates as climatic conditions change. Unfortunately, we have little understanding of how rapid evolutionary responses to climate play out when species are embedded in the competitive communities that they inhabit in nature. We tested the effects of rapid evolution in response to interspecific competition on subsequent ecological and evolutionary trajectories in a seasonally changing climate using a field-based evolution experiment with Drosophila melanogaster Populations of D. melanogaster were either exposed, or not exposed, to interspecific competition with an invasive competitor, Zaprionus indianus, over the summer. We then quantified these populations' ecological trajectories (abundances) and evolutionary trajectories (heritable phenotypic change) when exposed to a cooling fall climate. We found that competition with Z. indianus in the summer affected the subsequent evolutionary trajectory of D. melanogaster populations in the fall, after all interspecific competition had ceased. Specifically, flies with a history of interspecific competition evolved under fall conditions to be larger and have lower cold fecundity and faster development than flies without a history of interspecific competition. Surprisingly, this divergent fall evolutionary trajectory occurred in the absence of any detectible effect of the summer competitive environment on phenotypic evolution over the summer or population dynamics in the fall. This study demonstrates that competitive interactions can leave a legacy that shapes evolutionary responses to climate even after competition has ceased, and more broadly, that evolution in response to one selective pressure can fundamentally alter evolution in response to subsequent agents of selection.

Trait-fitness associations do not predict within-species phenotypic evolution over 2 million years

Long-term patterns of phenotypic change are the cumulative results of tens of thousands to millions of years of evolution. Yet, empirical and theoretical studies of phenotypic selection are largely based on contemporary populations. The challenges in studying phenotypic evolution, in particular trait-fitness associations in the deep past, are barriers to linking micro- and macroevolution. Here, we capitalize on the unique opportunity offered by a marine colonial organism commonly preserved in the fossil record to investigate trait-fitness associations over 2 Myr. We use the density of female polymorphs in colonies of Antartothoa tongima as a proxy for fecundity, a fitness component, and investigate multivariate signals of trait-fitness associations in six time intervals on the backdrop of Pleistocene climatic shifts. We detect negative trait-fitness associations for feeding polymorph (autozooid) sizes, positive associations for autozooid shape but no particular relationship between fecundity and brood chamber size. In addition, we demonstrate that long-term trait patterns are explained by palaeoclimate (as approximated by ∂¹⁸O), and to a lesser extent by ecological interactions (i.e. overgrowth competition and substrate crowding). Our analyses show that macroevolutionary outcomes of trait evolution are not a simple scaling-up from the trait-fitness associations.

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.

Critical Thermal Limits Do Not Vary between Wild-caught and Captive-bred Tadpoles of Agalychnis spurrelli (Anura: Hylidae)

Captive-bred organisms are widely used in ecology, evolution and conservation research, especially in scenarios where natural populations are scarce or at risk of extinction. Yet, it is still unclear whether captivity may alter thermal tolerances, crucial traits to predict species resilience to global warming. Here, we study whether captive-bred tadpoles of the gliding treefrog (Agalychnis spurrelli) show different thermal tolerances than wild-caught individuals. Our results show that there are no differences between critical thermal limits (CTmax and CTmin) of captive-bred and wild-caught tadpoles exposed to three-day acclimatization at 20 °C. Therefore, we suggest that the use of captive-bred amphibians is valid and may be appropriate in experimental comparisons to thermal physiological studies of wild 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.

Ecological genomics of adaptation to unpredictability in experimental rotifer populations

Elucidating the genetic basis of phenotypic variation in response to different environments is key to understanding how populations evolve. Facultatively sexual rotifers can develop adaptive responses to fluctuating environments. In a previous evolution experiment, diapause-related traits changed rapidly in response to two selective regimes (predictable vs unpredictable) in laboratory populations of the rotifer Brachionus plicatilis. Here, we investigate the genomic basis of adaptation to environmental unpredictability in these experimental populations. We identified and genotyped genome-wide polymorphisms in 169 clones from both selective regimes after seven cycles of selection using genotyping by sequencing (GBS). Additionally, we used GBS data from the 270 field clones from which the laboratory populations were established. This GBS dataset was used to identify candidate SNPs under selection. A total of 76 SNPs showed divergent selection, three of which are candidates for being under selection in the particular unpredictable fluctuation pattern studied. Most of the remaining SNPs showed strong signals of adaptation to laboratory conditions. Furthermore, a genotype-phenotype association approach revealed five SNPs associated with two key life-history traits in the adaptation to unpredictability. Our results contribute to elucidating the genomic basis for adaptation to unpredictable environments and lay the groundwork for future evolution studies in rotifers.

Major range loss predicted from lack of heat adaptability in an alpine Drosophila species

Climate warming is threatening biodiversity worldwide. Climate specialists such as alpine species are especially likely to be vulnerable. Adaptation by rapid evolution is the only long-term option for survival of many species, but the adaptive evolutionary potential of heat resistance has not been assessed in an alpine invertebrate. Here, we show that the alpine fly Drosophila nigrosparsa cannot readily adapt to heat stress. Heat-exposed flies from a regime with increased ambient temperature and a regime with increased temperature plus artificial selection for heat tolerance were less heat tolerant than the control group. Increased ambient temperature affected negatively both fitness and competitiveness. Ecological niche models predicted the loss of three quarters of the climatically habitable areas of this fly by the end of this century. Our findings suggest that, alongside with other climate specialists, species from mountainous regions are highly vulnerable to climate warming and unlikely to adapt through evolutionary genetic changes.

Changes in lipid classes of Drosophila melanogaster in response to selection for three stress traits

Laboratory selection on environmental stress traits is an evolutionary approach that is informative in the context of understanding stress adaptation. Here we characterize changes in a lipidome of Drosophila melanogaster in lines selected for increased heat (elevated heat knockdown refractoriness), cold (decreased time to recover from chill-coma) and desiccation survival. Selection for desiccation resistance resulted in changes in multiple lipid classes used to characterize a lipidome. This included a decrease in triacylglycerols (TAGs) which is relevant to interpretation of storage lipid levels in previous D. melanogaster desiccation survival selection experiments. Chill-coma recovery rate selection was expected to show extensive changes in lipid classes, but only phosphatidic acids exhibited significant change. Selection for increased heat knockdown resistance resulted in a substantial change in the abundance of a class of lipids (diacylglycerols) which could play a role in mediating the heat shock response or result in an increase in neutral lipid mobilization.

Genomic signatures of experimental adaptive radiation in Drosophila

Abiotic environmental factors play a fundamental role in determining the distribution, abundance and adaptive diversification of species. Empowered by new technologies enabling rapid and increasingly accurate examination of genomic variation in populations, researchers may gain new insights into the genomic background of adaptive radiation and stress resistance. We investigated genomic variation across generations of large-scale experimental selection regimes originating from a single founder population of Drosophila melanogaster, diverging in response to ecologically relevant environmental stressors: heat shock, heat knock down, cold shock, desiccation and starvation. When compared to the founder population, and to parallel unselected controls, there were more than 100,000 single nucleotide polymorphisms (SNPs) displaying consistent allelic changes in response to selective pressures across generations. These SNPs were found in both coding and noncoding sequences, with the highest density in promoter regions, and involved a broad range of functionalities, including molecular chaperoning by heat-shock proteins. The SNP patterns were highly stressor-specific despite considerable variation among line replicates within each selection regime, as reflected by a principal component analysis, and co-occurred with selective sweep regions. Only ~15% of SNPs with putatively adaptive changes were shared by at least two selective regimes, while less than 1% of SNPs diverged in opposite directions. Divergent stressors driving evolution in the experimental system of adaptive radiation left distinct genomic signatures, most pronounced in starvation and heat-shock selection regimes.

Gene duplicates resolving sexual conflict rapidly evolved essential gametogenesis functions

Males and females have different fitness optima but share the vast majority of their genomes, causing an inherent genetic conflict between the two sexes that must be resolved to achieve maximal population fitness. We show that two tandem duplicate genes found specifically in Drosophila melanogaster are sexually antagonistic, but rapidly evolved sex-specific functions and expression patterns that mitigate their antagonistic effects. We use copy-specific knockouts and rescue experiments to show that Apollo (Apl) is essential for male fertility but detrimental to female fertility, in addition to its important role in development, while Artemis (Arts) is essential for female fertility but detrimental to male fertility. Further analyses show that Apl and Arts have essential roles in spermatogenesis and oogenesis. These duplicates formed ~200,000 years ago, underwent a strong selective sweep and lost most expression in the antagonized sex. These data provide direct evidence that gene duplication allowed rapid mitigation of sexual conflict by allowing Apl and Arts to evolve essential sex-specific reproductive functions and complementary expression in male and female gonads.

Laboratory maintenance does not alter ecological and physiological patterns among species: a Drosophila case study

Large comparative studies in animal ecology, physiology and evolution often use animals reared in the laboratory for many generations; however, the relevance of these studies hinges on the assumption that laboratory populations are still representative for their wild living conspecifics. In this study, we investigate whether laboratory-maintained and freshly collected animal populations are fundamentally different and whether data from laboratory-maintained animals are valid to use in large comparative investigations of ecological and physiological patterns. Here, we obtained nine species of Drosophila with paired populations of laboratory-maintained and freshly collected flies. These species, representing a range of ecotypes, were assayed for four stress-tolerance, two body-size traits and six life-history traits. For all of these traits, we observed small differences in species-specific comparisons between field and laboratory populations; however, these differences were unsystematic and laboratory maintenance did not eclipse fundamental species characteristics. To investigate whether laboratory maintenance influence the general patterns in comparative studies, we correlated stress tolerance and life-history traits with environmental traits for the laboratory-maintained and freshly collected populations. Based on this analysis, we found that the comparative physiological and ecological trait correlations are similar irrespective of provenience. This finding is important for comparative biology in general because it validates comparative meta-analyses based on laboratory-maintained populations.

An attempt to select non-genetic variation in resistance to starvation and reduced chill coma recovery time in Drosophila melanogaster

Phenotypic variance is attributed to genetic and non-genetic factors, and only the former are supposed to be inherited and thus suitable for the action of selection. Although increasing amounts of data suggest that non-genetic variability may be inherited, we have limited empirical data in animals. Here, we performed an artificial selection experiment using Drosophila melanogaster inbred lines. We quantified the response to selection for a decrease in chill coma recovery time and an increase in starvation resistance. We observed a weak response to selection in the inbred and outbred lines, with variability across lines. At the end of the selection process, differential expression was detected for some genes associated with epigenetics, the piRNA pathway and canalization functions. As the selection process can disturb the canalization process and increase the phenotypic variance of developmental traits, we also investigated possible effects of the selection process on the number of scutellar bristles, fluctuating asymmetry levels, and fitness estimates. These results suggest that, contrary to what was shown in plants, selection of non-genetic variability is not straightforward in Drosophila and appears to be strongly genotype-dependent.

Plasticity and cross-tolerance to heterogeneous environments: divergent stress responses co-evolved in an African fruit fly

Plastic adjustments of physiological tolerance to a particular stressor can result in fitness benefits for resistance that might manifest not only in that same environment but also be advantageous when faced with alternative environmental stressors, a phenomenon termed 'cross-tolerance'. The nature and magnitude of cross-tolerance responses can provide important insights into the underlying genetic architecture, potential constraints on or versatility of an organism's stress responses. In this study, we tested for cross-tolerance to a suite of abiotic factors that likely contribute to setting insect population dynamics and geographic range limits: heat, cold, desiccation and starvation resistance in adult Ceratitis rosa following acclimation to all these isolated individual conditions prior to stress assays. Traits of stress resistance scored included critical thermal (activity) limits, chill coma recovery time (CCRT), heat knockdown time (HKDT), desiccation and starvation resistance. In agreement with other studies, we found that acclimation to one stress typically increased resistance for that same stress experienced later in life. A more novel outcome, however, is that here we also found substantial evidence for cross-tolerance. For example, we found an improvement in heat tolerance (critical thermal maxima, CT_max ) following starvation or desiccation hardening and improved desiccation resistance following cold acclimation, indicating pronounced cross-tolerance to these environmental stressors for the traits examined. We also found that two different traits of the same stress resistance differed in their responsiveness to the same stress conditions (e.g. HKDT was less cross-resistant than CT_max ). The results of this study have two major implications that are of broader importance: (i) that these traits likely co-evolved to cope with diverse or simultaneous stressors, and (ii) that a set of common underlying physiological mechanisms might exist between apparently divergent stress responses in this species. This species may prove to be a valuable model for future work on the evolutionary and mechanistic basis of cross-tolerance.

The role of genetic diversity and past-history selection pressures in the susceptibility of Chironomus riparius populations to environmental stress

Natural populations experiencing intense selection and genetic drift may exhibit limited potential to adapt to environmental change. The present study addresses the following aspects of the "genetic erosion" hypothesis in the midge Chironomus riparius: does long-term mercury (Hg) contamination affect the Hg tolerance of midge populations inhabiting such impacted areas? If so, is there any fitness cost under changing environmental conditions? And does genetic impoverishment influence the susceptibility of C. riparius to cope with environmental stressful conditions? For this end, we tested the acute and chronic tolerance to Hg and salinity in four C. riparius populations differing in their levels of genetic diversity (assessed through microsatellite markers) and past-histories of Hg exposure. Results showed that the midge population collected from a heavily Hg-contaminated site had higher Hg tolerance compared to the population collected from a closely-located reference site suggesting directional selection for Hg-tolerant traits in its native environment despite no genetic erosion in the field. No increased susceptibility under changing environmental conditions of salinity stress was observed. Moreover, results also showed that populations with higher genetic diversity performed better in the partial life-cycle assays providing evidence on the key role that genetic diversity plays as mediator of populations' susceptibility to environmental stress. Our findings are discussed in terms of the suitability of C. riparius as a model organism in evolutionary toxicology studies as well as the validity of ecotoxicological assessments using genetically eroded laboratory populations.

Nucleotide diversity inflation as a genome-wide response to experimental lifespan extension in Drosophila melanogaster

Background

Evolutionary theory predicts that antagonistically selected alleles, such as those with divergent pleiotropic effects in early and late life, may often reach intermediate population frequencies due to balancing selection, an elusive process when sought out empirically. Alternatively, genetic diversity may increase as a result of positive frequency-dependent selection and genetic purging in bottlenecked populations.

Results

While experimental evolution systems with directional phenotypic selection typically result in at least local heterozygosity loss, we report that selection for increased lifespan in Drosophila melanogaster leads to an extensive genome-wide increase of nucleotide diversity in the selected lines compared to replicate control lines, pronounced in regions with no or low recombination, such as chromosome 4 and centromere neighborhoods. These changes, particularly in coding sequences, are most consistent with the operation of balancing selection and the antagonistic pleiotropy theory of aging and life history traits that tend to be intercorrelated. Genes involved in antioxidant defenses, along with multiple lncRNAs, were among those most affected by balancing selection. Despite the overwhelming genetic diversification and the paucity of selective sweep regions, two genes with functions important for central nervous system and memory, Ptp10D and Ank2, evolved under positive selection in the longevity lines.

Conclusions

Overall, the ‘evolve-and-resequence’ experimental approach proves successful in providing unique insights into the complex evolutionary dynamics of genomic regions responsible for longevity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3485-0) contains supplementary material, which is available to authorized users.

Genomic Trajectories to Desiccation Resistance: Convergence and Divergence Among Replicate Selected Drosophila Lines

Adaptation to environmental stress is critical for long-term species persistence. With climate change and other anthropogenic stressors compounding natural selective pressures, understanding the nature of adaptation is as important as ever in evolutionary biology. In particular, the number of alternative molecular trajectories available for an organism to reach the same adaptive phenotype remains poorly understood. Here, we investigate this issue in a set of replicated Drosophila melanogaster lines selected for increased desiccation resistance-a classical physiological trait that has been closely linked to Drosophila species distributions. We used pooled whole-genome sequencing (Pool-Seq) to compare the genetic basis of their selection responses, using a matching set of replicated control lines for characterizing laboratory (lab-)adaptation, as well as the original base population. The ratio of effective population size to census size was high over the 21 generations of the experiment at 0.52-0.88 for all selected and control lines. While selected SNPs in replicates of the same treatment (desiccation-selection or lab-adaptation) tended to change frequency in the same direction, suggesting some commonality in the selection response, candidate SNP and gene lists often differed among replicates. Three of the five desiccation-selection replicates showed significant overlap at the gene and network level. All five replicates showed enrichment for ovary-expressed genes, suggesting maternal effects on the selected trait. Divergence between pairs of replicate lines for desiccation-candidate SNPs was greater than between pairs of control lines. This difference also far exceeded the divergence between pairs of replicate lines for neutral SNPs. Overall, while there was overlap in the direction of allele frequency changes and the network and functional categories affected by desiccation selection, replicates showed unique responses at all levels, likely reflecting hitchhiking effects, and highlighting the challenges in identifying candidate genes from these types of experiments when traits are likely to be polygenic.

A Meta-Analysis Evaluating the Relationship between Aquatic Contaminants and Chironomid Larval Deformities in Laboratory Studies

Chironomid larval deformities have been widely investigated as an aquatic pollution toxicity end point. Field chironomid surveys often show a spatial association between contaminants and deformities, suggesting contaminants cause deformities. However, over 40 years of laboratory assays have not been able to confirm this causality. We therefore conducted a review of the literature and meta-analysis, in order to (A) assess whether trends across assays indicated dose-response effects, (B) characterize the consistency of results, and (C) investigate whether experimental issues and publication bias were contributing to inconsistency and/or reducing confidence in results. The experimental issues we investigated were extraneous nonchemical laboratory stressors (which may mask or interact with chemical effects), and mortality (which can confound deformity results). Our meta-analysis of the most commonly tested chemicals suggested dose-response effects for copper, but not lead or zinc. However, we also found substantial inconsistency across studies. Both mortality and extraneous stressors were potentially contributing to this inconsistency, reducing confidence in most published data. We observed no evidence of publication bias. We conclude that any causal link between contaminants and deformities remains uncertain, and suggest improved experimental and data reporting procedures to better assess this relationship.

Ancestral population reconstitution from isofemale lines as a tool for experimental evolution

Experimental evolution is a powerful tool to study adaptation under controlled conditions. Laboratory natural selection experiments mimic adaptation in the wild with better‐adapted genotypes having more offspring. Because the selected traits are frequently not known, adaptation is typically measured as fitness increase by comparing evolved populations against an unselected reference population maintained in a laboratory environment. With adaptation to the laboratory conditions and genetic drift, however, it is not clear to what extent such comparisons provide unbiased estimates of adaptation. Alternatively, ancestral variation could be preserved in isofemale lines that can be combined to reconstitute the ancestral population. Here, we assess the impact of selection on alleles segregating in newly established Drosophila isofemale lines. We reconstituted two populations from isofemale lines and compared them to two original ancestral populations (AP) founded from the same lines shortly after collection. No significant allele frequency changes could be detected between both AP and simulations showed that drift had a low impact compared to Pool‐Seq‐associated sampling effects. We conclude that laboratory selection on segregating variation in isofemale lines is too weak to have detectable effects, which validates ancestral population reconstitution from isofemale lines as an unbiased approach for measuring adaptation in evolved populations.

Artificial selection on chill-coma recovery time in Drosophila melanogaster: Direct and correlated responses to selection

Artificial selection can be used to create populations with extreme phenotypic responses to environmental stressors. When artificial selection is applied to a single component of a stress response, this selection may result in correlated responses in other stress responses, a phenomenon called cross-tolerance, which is ultimately controlled by the genetic correlations among traits. We selected for extreme responses to cold tolerance by selecting for chill-coma recovery time from a single temperate population of Drosophila melanogaster. Chill-coma recovery time is a common metric of low, but non-lethal, cold temperature tolerance. Replicated divergent artificial selection was applied to a genetically variable base population for 31 generations, resulting in two cold resistant, two cold susceptible, and two unselected control lines. To quantify the relationship between selection on chill-coma recovery and other metrics of thermal performance, we also measured survivorship after acute cold exposure, survivorship after chronic cold exposure, survivorship after cold exposure following a pre-treatment period (rapid cold hardening), starvation tolerance, and heat tolerance. We find that chill-coma recovery time is heritable within this population and that there is an asymmetric response to increased and decreased chill-coma recovery time. Surprisingly, we found no cross-tolerances between selection on chill-coma recovery time and the other environmental stress response traits. These results suggest that although artificial selection has dramatically altered chill-coma recovery time, the correlated response to selection on other stress response phenotypes has been negligible. The lack of a correlated response suggests that chill-coma recovery time in these selection lines is likely genetically independent from measures of cold survivorship tested here.

The Influence of Insecticide Resistance, Age, Sex, and Blood Feeding Frequency on Thermal Tolerance of Wild and Laboratory Phenotypes of Anopheles funestus (Diptera: Culicidae)

Resistance to insecticides is a global phenomenon and is increasing at an unprecedented rate. How resistant and susceptible strains of malaria vectors might differ in terms of life history and basic biology is often overlooked, despite the potential importance of such information in light of changing climates. Here, we investigated the upper thermal limits (ULT50) of wild and laboratory strains of Anopheles funestus Giles mosquitoes, including resistance status, sex, age, and blood feeding status as potential factors influencing ULT50. No significant differences in ULT50 were observed between strains displaying different resistance patterns, nor was there a significant difference between wild and laboratory strains. In some instances, strains showed a senescence response, displaying decreased ULT50 with an increase in age, and differences between males and females (females displaying higher ULT50 than males). Blood feeding did not seem to influence ULT50 in any way. For An. funestus, it seems evident that there is no cost to resistance despite what is displayed in other anopheline species. This could have significant impacts for vector control, with resistant populations of An. funestus performing just as well, if not better, than susceptible strains, especially under changing environmental conditions such as those expected to occur with climate change.

Parallel trait adaptation across opposing thermal environments in experimental Drosophila melanogaster populations

Thermal stress is a pervasive selective agent in natural populations that impacts organismal growth, survival, and reproduction. Drosophila melanogaster exhibits a variety of putatively adaptive phenotypic responses to thermal stress in natural and experimental settings; however, accompanying assessments of fitness are typically lacking. Here, we quantify changes in fitness and known thermal tolerance traits in replicated experimental D. melanogaster populations following more than 40 generations of evolution to either cyclic cold or hot temperatures. By evaluating fitness for both evolved populations alongside a reconstituted starting population, we show that the evolved populations were the best adapted within their respective thermal environments. More strikingly, the evolved populations exhibited increased fitness in both environments and improved resistance to both acute heat and cold stress. This unexpected parallel response appeared to be an adaptation to the rapid temperature changes that drove the cycling thermal regimes, as parallel fitness changes were not observed when tested in a constant thermal environment. Our results add to a small, but growing group of studies that demonstrate the importance of fluctuating temperature changes for thermal adaptation and highlight the need for additional work in this area.

Strong Costs and Benefits of Winter Acclimatization in Drosophila melanogaster

Studies on thermal acclimation in insects are often performed on animals acclimated in the laboratory under conditions that are not ecologically relevant. Costs and benefits of acclimation responses under such conditions may not reflect costs and benefits in natural populations subjected to daily and seasonal temperature fluctuations. Here we estimated costs and benefits in thermal tolerance limits in relation to winter acclimatization of Drosophila melanogaster. We sampled flies from a natural habitat during winter in Denmark (field flies) and compared heat and cold tolerance of these to that of flies collected from the same natural population, but acclimated to 25 °C or 13 °C in the laboratory (laboratory flies). We further obtained thermal performance curves for egg-to-adult viability of field and laboratory (25 °C) flies, to estimate possible cross-generational effects of acclimation. We found much higher cold tolerance and a lowered heat tolerance in field flies compared to laboratory flies reared at 25 °C. Flies reared in the laboratory at 13 °C exhibited the same thermal cost-benefit relations as the winter acclimatized flies. We also found a cost of winter acclimatization in terms of decreased egg-to-adult viability at high temperatures of eggs laid by winter acclimatized flies. Based on our findings we suggest that winter acclimatization in nature can induce strong benefits in terms of increased cold tolerance. These benefits can be reproduced in the laboratory under ecologically relevant rearing and testing conditions, and should be incorporated in species distribution modelling. Winter acclimatization also leads to decreased heat tolerance. This may create a mismatch between acclimation responses and the thermal environment, e.g. if temperatures suddenly increase during spring, under current and expected more variable future climatic conditions.

Whole-Genome Resequencing of Experimental Populations Reveals Polygenic Basis of Egg-Size Variation in Drosophila melanogaster

Complete genome resequencing of populations holds great promise in deconstructing complex polygenic traits to elucidate molecular and developmental mechanisms of adaptation. Egg size is a classic adaptive trait in insects, birds, and other taxa, but its highly polygenic architecture has prevented high-resolution genetic analysis. We used replicated experimental evolution in Drosophila melanogaster and whole-genome sequencing to identify consistent signatures of polygenic egg-size adaptation. A generalized linear-mixed model revealed reproducible allele frequency differences between replicated experimental populations selected for large and small egg volumes at approximately 4,000 single nucleotide polymorphisms (SNPs). Several hundred distinct genomic regions contain clusters of these SNPs and have lower heterozygosity than the genomic background, consistent with selection acting on polymorphisms in these regions. These SNPs are also enriched among genes expressed in Drosophila ovaries and many of these genes have well-defined functions in Drosophila oogenesis. Additional genes regulating egg development, growth, and cell size show evidence of directional selection as genes regulating these biological processes are enriched for highly differentiated SNPs. Genetic crosses performed with a subset of candidate genes demonstrated that these genes influence egg size, at least in the large genetic background. These findings confirm the highly polygenic architecture of this adaptive trait, and suggest the involvement of many novel candidate genes in regulating egg size.

Physiological basis of starvation resistance in Drosophila leontia: analysis of sexual dimorphism

Geographically varying starvation stress has often been considered as a natural selector that constrains between-population differences for starvation resistance (SR) in Drosophila species. On the Indian subcontinent, a dozen Drosophila species have shown clinal variations in SR across latitude, but the evolved physiological basis of such contrasting adaptations is largely unknown. In the present study, I untangled the physiological basis of sex-specific as well as between-population divergence for SR in D. leontia, collected across a latitudinal transect of the Indian subcontinent (11°45'-31°19'N). Secondly, I tested the assumptions that hardening to starvation stress facilitates an increased survival under subsequent lethal levels of starvation, and such plastic effects differ between the sexes. I observed several interesting results. In contrast to a steeper cline of starvation-related traits with latitude in females, a shallower gradient was observed for males. Females stored higher (~1.3-fold) dry-mass-specific levels of body lipids and glycogen contents, and utilized these both of these energy resources under starvation stress, whereas the starved males metabolized only body lipids as a source of energy. Conversely, the rate of body lipid utilization and threshold need were considerably higher in females as compared with males. Between-population differences were significant for storage levels of energy reserves only, but not for other avenues (rate of metabolite utilization and threshold need) of SR for both sexes. These findings indicate that multiple pathways shape the physiological basis of sexual dimorphism for SR in D. leontia. Further, single or multiple bouts of starvation hardening conferred an increased longevity (~4-9 h; P<0.001) under subsequent lethal levels of starvation stress for females only, and such plastic responses were consistent with a decrease in rate of metabolite utilization. Nevertheless, between-population effects were non-significant for absolute hardening capacity (AHC=KSR-C). Altogether, these findings suggest that similar evolutionary constraints have resulted in divergent genetic as well as plastic responses to evolve adaptations under starvation stress, and account for the observed sexual dimorphism for basal SR in D. leontia.

Artificial selection on larval growth curves in Tribolium: correlated responses and constraints

Body size is often constrained from evolving. Although artificial selection on body size in insects frequently results in a sizable response, these responses usually bear fitness costs. Further, these experiments tend to select only on size at one landmark age, rather than selecting for patterns of growth over the whole larval life stage. To address whether constraints may be caused by larval growth patterns rather than final size, we implemented a function-valued (FV) trait method of selection, in which entire larval growth curves from Tribolium were artificially selected. The selection gradient function used was previously predicted to give the maximal response and was implemented using a novel selection index in the FV framework. Results indicated a significant response after one generation of selection, but no response in subsequent generations. Correlated responses included increased mortality, increased critical weight, and decreased development time (DT). The lack of response in size and development time after the first generation was likely caused by increased mortality suffered in selected lines; we demonstrated that the selection criterion caused both increased body size and increased mortality. We conclude that artificial selection on continuous traits using FV methods is very efficient and that the constraint of body size evolution is likely caused by a suite of trade-offs with other traits.

Fast evolutionary genetic differentiation during experimental colonizations

Founder effects during colonization of a novel environment are expected to change the genetic composition of populations, leading to differentiation between the colonizer population and its source population. Another expected outcome is differentiation among populations derived from repeated independent colonizations starting from the same source. We have previously detected significant founder effects affecting rate of laboratory adaptation among Drosophila subobscura laboratory populations derived from the wild. We also showed that during the first generations in the laboratory, considerable genetic differentiation occurs between foundations. The present study deepens that analysis, taking into account the natural sampling hierarchy of six foundations, derived from different locations, different years and from two samples in one of the years. We show that striking stochastic effects occur in the first two generations of laboratory culture, effects that produce immediate differentiation between foundations, independent of the source of origin and despite similarity among all founders. This divergence is probably due to powerful genetic sampling effects during the first few generations of culture in the novel laboratory environment, as a result of a significant drop in Ne. Changes in demography as well as high variance in reproductive success in the novel environment may contribute to the low values of Ne. This study shows that estimates of genetic differentiation between natural populations may be accurate when based on the initial samples collected in the wild, though considerable genetic differentiation may occur in the very first generations of evolution in a new, confined environment. Rapid and significant evolutionary changes can thus occur during the early generations of a founding event, both in the wild and under domestication, effects of interest for both scientific and conservation purposes.

Effective population size and evolutionary dynamics in outbred laboratory populations of Drosophila

Census population size, sex-ratio and female reproductive success were monitored in 10 laboratory populations of Drosophila melanogaster selected for different ages of reproduction. With this demographic information, we estimated eigenvalue, variance and probability of allele loss effective population sizes. We conclude that estimates of effective size based on genefrequency change at a few loci are biased downwards. We analysed the relative roles of selection and genetic drift in maintaining genetic variation in laboratory populations of Drosophila. We suggest that rare, favourable genetic variants in our laboratory populations have a high chance of being lost if their fitness effect is weak, e.g. 1% or less. However, if the fitness effect of this variation is 10% or greater, these rare variants are likely to increase to high frequency. The demographic information developed in this study suggests that some of our laboratory populations harbour more genetic variation than expected. One explanation for this finding is that part of the genetic variation in these outbred laboratory Drosophila populations may be maintained by some form of balancing selection. We suggest that, unlike bacteria, medium-term adaptation of laboratory populations of fruit flies is not primarily driven by new mutations, but rather by changes in the frequency of preexisting alleles.

A Drosophila laboratory evolution experiment points to low evolutionary potential under increased temperatures likely to be experienced in the future

The ability to respond evolutionarily to increasing temperatures is important for survival of ectotherms in a changing climate. Recent studies suggest that upper thermal limits may be evolutionary constrained. We address this hypothesis in a laboratory evolution experiment, encompassing ecologically relevant thermal regimes. To examine the potential for species to respond to climate change, we exposed replicate populations of Drosophila melanogaster to increasing temperatures (0.3 °C every generation) for 20 generations, whereas corresponding replicate control populations were held at benign thermal conditions throughout the experiment. We hypothesized that replicate populations exposed to increasing temperatures would show increased resistance to warm and dry environments compared with replicate control populations. Contrasting replicate populations held at the two thermal regimes showed (i) an increase in desiccation resistance and a decline in heat knock-down resistance in replicate populations exposed to increasing temperatures, (ii) similar egg-to-adult viability and fecundity in replicate populations from the two thermal regimes, when assessed at high stressful temperatures and (iii) no difference in nucleotide diversity between thermal regimes. The limited scope for adaptive evolutionary responses shown in this study highlights the challenges faced by ectotherms under climate change.

Confirming candidate genes for longevity by RT-qPCR using two different genetic backgrounds and selection methods

Elucidating genes that affect life span or that can be used as biomarkers for ageing has received attention in diverse studies in recent years. Using model organisms and various approaches several genes have been linked to the longevity phenotype. For Drosophila melanogaster those studies have usually focussed on one sex and on flies originating from one genetic background, and results from different studies often do not overlap. Using D. melanogaster selected for increased longevity we aimed to find robust longevity related genes by examining gene expression in both sexes of flies originating from different genetic backgrounds. Further, we compared expression changes across three ages, when flies were young, middle aged or old, to examine how candidate gene expression changes with the onset of ageing. We selected 10 genes based on their expression differences in prior microarray studies. For about 50% of these we confirmed their potential as a candidate longevity gene. We found one robust candidate gene for longevity, CG32638. Three other genes, CG8934, mRpS10 and Spn43Ad, showed a tendency to be involved in life span determination in both backgrounds tested.