Rapid genomic changes in Drosophila melanogaster adapting to desiccation stress in an experimental evolution system

No single rescue recipe: genome complexities modulate insect response to climate change

Declines in insect populations have gained formidable attention. Given their crucial role in the ecosystem, the causes of declining insect populations must be investigated. However, the insect clade has been associated with low extinction and high diversification rates. It is unlikely that insects underwent mass extinctions in the past. However, the pace of current climate change could make insect populations vulnerable to extinction. We propose genome size (GS) and transposable elements (TEs) to be rough estimates to assess extinction risk. Larger GS and/or proliferating TEs have been associated with adaptation in rapid climate change scenarios. We speculate that unstable, stressful environmental conditions are strongly associated with GS and TE expansion, which could be further correlated with adaptations. Alternately, stressful conditions trigger TE bursts that are not purged in smaller populations. GS and TE loads could be indicators of small effective populations in the wild, likely experiencing bottlenecks or drastic climatic perturbations, which calls for an urgent assessment of extinction risk.

Genetic basis and repeatability for desiccation resistance in Drosophila melanogaster (Diptera: Drosophilidae)

Dehydration is a stress factor for organisms inhabiting natural habitats where water is scarce. Thus, it may be expected that species facing arid environments will develop mechanisms that maximize resistance to desiccation. Insects are excellent models for studying the effects of dehydration as well as the mechanisms and processes that prevent water loss since the effect of desiccation is greater due to the higher area/volume ratio than larger animals. Even though physiological and behavioral mechanisms to cope with desiccation are being understood, the genetic basis underlying the mechanisms related to variation in desiccation resistance and the context-dependent effect remain unsolved. Here we analyze the genetic bases of desiccation resistance in Drosophila melanogaster and identify candidate genes that underlie trait variation. Our quantitative genetic analysis of desiccation resistance revealed sexual dimorphism and extensive genetic variation. The phenotype-genotype association analyses (GWAS) identified 71 candidate genes responsible for total phenotypic variation in desiccation resistance. Half of these candidate genes were sex-specific suggesting that the genetic architecture underlying this adaptive trait differs between males and females. Moreover, the public availability of desiccation data analyzed on the same lines but in a different lab allows us to investigate the reliability and repeatability of results obtained in independent screens. Our survey indicates a pervasive micro-environment lab-dependent effect since we did not detect overlap in the sets of genes affecting desiccation resistance identified between labs.

Gene expression differences consistent with water loss reduction underlie desiccation tolerance of natural Drosophila populations

BACKGROUND

Climate change is one of the main factors shaping the distribution and biodiversity of organisms, among others by greatly altering water availability, thus exposing species and ecosystems to harsh desiccation conditions. However, most of the studies so far have focused on the effects of increased temperature. Integrating transcriptomics and physiology is key to advancing our knowledge on how species cope with desiccation stress, and these studies are still best accomplished in model organisms.

RESULTS

Here, we characterized the natural variation of European D. melanogaster populations across climate zones and found that strains from arid regions were similar or more tolerant to desiccation compared with strains from temperate regions. Tolerant and sensitive strains differed not only in their transcriptomic response to stress but also in their basal expression levels. We further showed that gene expression changes in tolerant strains correlated with their physiological response to desiccation stress and with their cuticular hydrocarbon composition, and functionally validated three of the candidate genes identified. Transposable elements, which are known to influence stress response across organisms, were not found to be enriched nearby differentially expressed genes. Finally, we identified several tRNA-derived small RNA fragments that differentially targeted genes in response to desiccation stress.

CONCLUSIONS

Overall, our results showed that basal gene expression differences across individuals should be analyzed if we are to understand the genetic basis of differential stress survival. Moreover, tRNA-derived small RNA fragments appear to be relevant across stress responses and allow for the identification of stress-response genes not detected at the transcriptional level.

Selection of p53 pathway in adaptive evolution and reproductive isolation in incipient sympatric speciation of Drosophila at Evolution Canyon

Sympatric speciation (SS) refers to the origin of new species within a freely breeding population. The ‘Evolution Canyon’ (EC) in Israel is a natural microsite model for SS of species across phylogenies from viruses and bacteria to mammals, adapting to, and speciating in, interslope microclimates. The cosmopolitan Drosophila melanogaster at EC I, Mount Carmel, is undergoing incipient SS in response to sharply divergent interslope microclimate stresses, including solar radiation, temperature, humidity and pathogenicity. We demonstrated here a selective interslope divergence of single nucleotide polymorphism (SNP) distribution in the Drosophila p53 pathway. This involves a total of 71 genes, which are associated with DNA repair, heat response, and fungal and bacterial resistant pathways. This distribution pattern links the previously observed thermotolerance and ageing divergence of D. melanogaster between the opposite canyon slopes: the south-facing slope (SFS, or African slope: tropical, savannoid and dry) and the abutting north-facing slope (NFS, or European slope; temperate, forested, cool and humid). The genes with interslope-significant differential SNPs link the p53 pathway with pathways related to the responses to microclimates through protein-protein interaction. Moreover, for the first time we provide evidence that the p53 pathway is linked to reproductive isolation, and is thus actively participating in incipient SS of D. melanogaster. This is the first demonstration of a link between the p53 pathway and reproductive isolation, thereby contributing to adaptive incipient sympatric speciation.

Dehydration Dynamics in Terrestrial Arthropods: From Water Sensing to Trophic Interactions

Since the transition from water to land, maintaining water balance has been a key challenge for terrestrial arthropods. We explore factors that allow terrestrial arthropods to survive within a variably dry world and how they shape ecological interactions. Detection of water and hydration is critical for maintaining water content. Efficient regulation of internal water content is accomplished by excretory and osmoregulatory systems that balance water intake and loss. Biochemical and physiological responses are necessary as water content declines to prevent and repair the damage that occurs during dehydration. Desiccation avoidance can occur seasonally or daily via a move to more favorable areas. Dehydration and its avoidance have ecological impacts that extend beyond a single species to alter trophic interactions. As climate changes, evolutionary and ecological processes will be critical to species survival during drought.

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.

Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster

Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).

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.

Eco-genetics of desiccation resistance in Drosophila

Climate change globally perturbs water circulation thereby influencing ecosystems including cultivated land. Both harmful and beneficial species of insects are likely to be vulnerable to such changes in climate. As small animals with a disadvantageous surface area to body mass ratio, they face a risk of desiccation. A number of behavioural, physiological and genetic strategies are deployed to solve these problems during adaptation in various Drosophila species. Over 100 desiccation-related genes have been identified in laboratory and wild populations of the cosmopolitan fruit fly Drosophila melanogaster and its sister species in large-scale and single-gene approaches. These genes are involved in water sensing and homeostasis, and barrier formation and function via the production and composition of surface lipids and via pigmentation. Interestingly, the genetic strategy implemented in a given population appears to be unpredictable. In part, this may be due to different experimental approaches in different studies. The observed variability may also reflect a rich standing genetic variation in Drosophila allowing a quasi-random choice of response strategies through soft-sweep events, although further studies are needed to unravel any underlying principles. These findings underline that D. melanogaster is a robust species well adapted to resist climate change-related desiccation. The rich data obtained in Drosophila research provide a framework to address and understand desiccation resistance in other insects. Through the application of powerful genetic tools in the model organism D. melanogaster, the functions of desiccation-related genes revealed by correlative studies can be tested and the underlying molecular mechanisms of desiccation tolerance understood. The combination of the wealth of available data and its genetic accessibility makes Drosophila an ideal bioindicator. Accumulation of data on desiccation resistance in Drosophila may allow us to create a world map of genetic evolution in response to climate change in an insect genome. Ultimately these efforts may provide guidelines for dealing with the effects of climate-related perturbations on insect population dynamics in the future.

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.

Intercontinental genomic parallelism in multiple three-spined stickleback adaptive radiations

Parallelism, the evolution of similar traits in populations diversifying in similar conditions, provides strong evidence of adaptation by natural selection. Many studies of parallelism focus on comparisons of different ecotypes or contrasting environments, defined a priori, which could upwardly bias the apparent prevalence of parallelism. Here, we estimated genomic parallelism associated with components of environmental and phenotypic variation at an intercontinental scale across four freshwater adaptive radiations (Alaska, British Columbia, Iceland and Scotland) of the three-spined stickleback (Gasterosteus aculeatus). We combined large-scale biological sampling and phenotyping with restriction site associated DNA sequencing (RAD-Seq) data from 73 freshwater lake populations and four marine ones (1,380 fish) to associate genome-wide allele frequencies with continuous distributions of environmental and phenotypic variation. Our three main findings demonstrate that (1) quantitative variation in phenotypes and environments can predict genomic parallelism; (2) genomic parallelism at the early stages of adaptive radiations, even at large geographic scales, is founded on standing variation; and (3) similar environments are a better predictor of genome-wide parallelism than similar phenotypes. Overall, this study validates the importance and predictive power of major phenotypic and environmental factors likely to influence the emergence of common patterns of genomic divergence, providing a clearer picture than analyses of dichotomous phenotypes and environments.

Resource quality determines the evolution of resistance and its genetic basis

Parasites impose strong selection on their hosts, but the level of any evolved resistance may be constrained by the availability of resources. However, studies identifying the genomic basis of such resource-mediated selection are rare, particularly in nonmodel organisms. Here, we investigated the role of nutrition in the evolution of resistance to a DNA virus (PiGV), and any associated trade-offs in a lepidopteran pest species (Plodia interpunctella). Through selection experiments and whole-genome resequencing, we identify genetic markers of resistance that vary between the nutritional environments during selection. We do not find consistent evolution of resistance in the presence of virus but rather see substantial variation among replicate populations. Resistance in a low-nutrition environment is negatively correlated with growth rate, consistent with an established trade-off between immunity and development, but this relationship is highly context dependent. Whole-genome resequencing of the host shows that resistance mechanisms are likely to be highly polygenic and although the underlying genetic architecture may differ between high and low-nutrition environments, similar mechanisms are commonly used. As a whole, our results emphasize the importance of the resource environment on influencing the evolution of resistance.

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.

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

Genomic divergence and adaptive convergence in Drosophila simulans from Evolution Canyon, Israel

Biodiversity refugia formed by unique features of the Mediterranean arid landscape, such as the dramatic ecological contrast of "Evolution Canyon," provide a natural laboratory in which local adaptations to divergent microclimate conditions can be investigated. Significant insights have been provided by studies of Drosophila melanogaster diversifying along the thermal gradient in Evolution Canyon, but a comparative framework to survey adaptive convergence across sister species at the site has been lacking. To fill this void, we present an analysis of genomic polymorphism and evolutionary divergence of Drosophila simulans, a close relative of Drosophila melanogaster with which it co-occurs on both slopes of the canyon. Our results show even deeper interslope divergence in D. simulans than in D. melanogaster, with extensive signatures of selective sweeps present in flies from both slopes but enhanced in the population from the hotter and drier south-facing slope. Interslope divergence was enriched for genes related to electrochemical balance and transmembrane transport, likely in response to increased selection for dehydration resistance on the hotter slope. Both species shared genomic regions that underwent major selective sweeps, but the overall level of adaptive convergence was low, demonstrating no shortage of alternative genomic solutions to cope with the challenges of the microclimate contrast. Mobile elements were a major source of genetic polymorphism and divergence, affecting all parts of the genome, including coding sequences of mating behavior-related genes.

Dissection of Complex, Fitness-Related Traits in Multiple Drosophila Mapping Populations Offers Insight into the Genetic Control of Stress Resistance

We leverage two complementary Drosophila melanogaster mapping panels to genetically dissect starvation resistance-an important fitness trait. Using >1600 genotypes from the multiparental Drosophila Synthetic Population Resource (DSPR), we map numerous starvation stress QTL that collectively explain a substantial fraction of trait heritability. Mapped QTL effects allowed us to estimate DSPR founder phenotypes, predictions that were correlated with the actual phenotypes of these lines. We observe a modest phenotypic correlation between starvation resistance and triglyceride level, traits that have been linked in previous studies. However, overlap among QTL identified for each trait is low. Since we also show that DSPR strains with extreme starvation phenotypes differ in desiccation resistance and activity level, our data imply multiple physiological mechanisms contribute to starvation variability. We additionally exploited the Drosophila Genetic Reference Panel (DGRP) to identify sequence variants associated with starvation resistance. Consistent with prior work these sites rarely fall within QTL intervals mapped in the DSPR. We were offered a unique opportunity to directly compare association mapping results across laboratories since two other groups previously measured starvation resistance in the DGRP. We found strong phenotypic correlations among studies, but extremely low overlap in the sets of genomewide significant sites. Despite this, our analyses revealed that the most highly associated variants from each study typically showed the same additive effect sign in independent studies, in contrast to otherwise equivalent sets of random variants. This consistency provides evidence for reproducible trait-associated sites in a widely used mapping panel, and highlights the polygenic nature of starvation resistance.

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.

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.

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.

Spatiotemporal dynamics and genome-wide association genome-wide association analysis of desiccation tolerance in Drosophila melanogaster

Water availability is a major environmental challenge to a variety of terrestrial organisms. In insects, desiccation tolerance varies predictably over spatial and temporal scales and is an important physiological determinant of fitness in natural populations. Here, we examine the dynamics of desiccation tolerance in North American populations of Drosophila melanogaster using: (a) natural populations sampled across latitudes and seasons; (b) experimental evolution in field mesocosms over seasonal time; (c) genome-wide associations to identify SNPs/genes associated with variation for desiccation tolerance; and (d) subsequent analysis of patterns of clinal/seasonal enrichment in existing pooled sequencing data of populations sampled in both North America and Australia. A cline in desiccation tolerance was observed, for which tolerance exhibited a positive association with latitude; tolerance also varied predictably with culture temperature, demonstrating a significant degree of thermal plasticity. Desiccation tolerance evolved rapidly in field mesocosms, although only males showed differences in desiccation tolerance between spring and autumn collections from natural populations. Water loss rates did not vary significantly among latitudinal or seasonal populations; however, changes in metabolic rates during prolonged exposure to dry conditions are consistent with increased tolerance in higher latitude populations. Genome-wide associations in a panel of inbred lines identified twenty-five SNPs in twenty-one loci associated with sex-averaged desiccation tolerance, but there is no robust signal of spatially varying selection on genes associated with desiccation tolerance. Together, our results suggest that desiccation tolerance is a complex and important fitness component that evolves rapidly and predictably in natural populations.

Whole genome SNP-associated signatures of local adaptation in honeybees of the Iberian Peninsula

The availability of powerful high-throughput genomic tools, combined with genome scans, has helped identifying genes and genetic changes responsible for environmental adaptation in many organisms, including the honeybee. Here, we resequenced 87 whole genomes of the honeybee native to Iberia and used conceptually different selection methods (Samβada, LFMM, PCAdapt, iHs) together with in sillico protein modelling to search for selection footprints along environmental gradients. We found 670 outlier SNPs, most of which associated with precipitation, longitude and latitude. Over 88.7% SNPs laid outside exons and there was a significant enrichment in regions adjacent to exons and UTRs. Enrichment was also detected in exonic regions. Furthermore, in silico protein modelling suggests that several non-synonymous SNPs are likely direct targets of selection, as they lead to amino acid replacements in functionally important sites of proteins. We identified genomic signatures of local adaptation in 140 genes, many of which are putatively implicated in fitness-related functions such as reproduction, immunity, olfaction, lipid biosynthesis and circadian clock. Our genome scan suggests that local adaptation in the Iberian honeybee involves variations in regions that might alter patterns of gene expression and in protein-coding genes, which are promising candidates to underpin adaptive change in the honeybee.

Renal neuroendocrine control of desiccation and cold tolerance by Drosophila suzukii

BACKGROUND

Neuropeptides are central to the regulation of physiological and behavioural processes in insects, directly impacting cold and desiccation survival. However, little is known about the control mechanisms governing these responses in Drosophila suzukii. The close phylogenetic relationship of D. suzukii with Drosophila melanogaster allows, through genomic and functional studies, an insight into the mechanisms directing stress tolerance in D. suzukii.

RESULTS

Capability (Capa), leucokinin (LK), diuretic hormone 44 (DH44 ) and DH31 neuropeptides demonstrated a high level of conservation between D. suzukii and D. melanogaster with respect to peptide sequences, neuronal expression, receptor localisation, and diuretic function in the Malpighian tubules. Despite D. suzukii's ability to populate cold environments, it proved sensitive to both cold and desiccation. Furthermore, in D. suzukii, Capa acts as a desiccation- and cold stress-responsive gene, while DH44 gene expression is increased only after desiccation exposure, and the LK gene after nonlethal cold stress recovery.

CONCLUSION

This study provides a comparative investigation into stress tolerance mediation by neuroendocrine signalling in two Drosophila species, providing evidence that similar signalling pathways control fluid secretion in the Malpighian tubules. Identifying processes governing specific environmental stresses affecting D. suzukii could lead to the development of targeted integrated management strategies to control insect pest populations. © 2017 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Directional selection to improve the sterile insect technique: Survival and sexual performance of desiccation resistant Anastrepha ludens strains

The sterile insect technique (SIT) is an effective, environmentally friendly method for insect control whose success depends on the sexual performance and survival of sterile males. These two parameters are influenced by environmental conditions of target areas, and releasing insects with a higher tolerance to stressful environments can improve SIT efficiency. Directional selection can be used to produce insect strains with higher tolerance to extreme environmental conditions, such as low humidity, for extended periods. We evaluated, under field cage conditions, the sexual competitiveness, sexual compatibility, and survival of strains of Anastrepha ludens (Loew) selected for desiccation resistance to determine the value of directional selection as a possible approach to enhance SIT efficiency. Fly strains (selected and unselected and those mass-reared) were exposed to stressful conditions of low humidity and food and water deprivation for 24 hr before test. As a control, mild conditions without the stressors were used. No differences in sexual competitiveness and sexual compatibility between selected, nonselected, and mass-reared strains were observed when previously exposed to mild conditions. Thus, selection for desiccation resistance does not modified negatively the sexual performance. However, when insects were exposed to stressful conditions, males of selected strains sexually outperform mass-reared males. Additionally, selected strains presented higher survival than mass-reared flies. The approach to integrate directional selection with other technologies in the SIT as well as the implications of using a desiccation-selected strain in the current pest management program is discussed.

Identifying consistent allele frequency differences in studies of stratified populations

With increasing application of pooled‐sequencing approaches to population genomics robust methods are needed to accurately quantify allele frequency differences between populations. Identifying consistent differences across stratified populations can allow us to detect genomic regions under selection and that differ between populations with different histories or attributes. Current popular statistical tests are easily implemented in widely available software tools which make them simple for researchers to apply. However, there are potential problems with the way such tests are used, which means that underlying assumptions about the data are frequently violated.

These problems are highlighted by simulation of simple but realistic population genetic models of neutral evolution and the performance of different tests are assessed. We present alternative tests (including Generalised Linear Models [GLMs] with quasibinomial error structure) with attractive properties for the analysis of allele frequency differences and re‐analyse a published dataset.

The simulations show that common statistical tests for consistent allele frequency differences perform poorly, with high false positive rates. Applying tests that do not confound heterogeneity and main effects significantly improves inference. Variation in sequencing coverage likely produces many false positives and re‐scaling allele frequencies to counts out of a common value or an effective sample size reduces this effect.

Many researchers are interested in identifying allele frequencies that vary consistently across replicates to identify loci underlying phenotypic responses to selection or natural variation in phenotypes. Popular methods that have been suggested for this task perform poorly in simulations. Overall, quasibinomial GLMs perform better and also have the attractive feature of allowing correction for multiple testing by standard procedures and are easily extended to other designs.

The cuticular nature of corneal lenses in Drosophila melanogaster

The dioptric visual system relies on precisely focusing lenses that project light onto a neural retina. While the proteins that constitute the lenses of many vertebrates are relatively well characterized, less is known about the proteins that constitute invertebrate lenses, especially the lens facets in insect compound eyes. To address this question, we used mass spectrophotometry to define the major proteins that comprise the corneal lenses from the adult Drosophila melanogaster compound eye. This led to the identification of four cuticular proteins: two previously identified lens proteins, drosocrystallin and retinin, and two newly identified proteins, Cpr66D and Cpr72Ec. To determine which ommatidial cells contribute each of these proteins to the lens, we conducted in situ hybridization at 50% pupal development, a key age for lens secretion. Our results confirm previous reports that drosocrystallin and retinin are expressed in the two primary corneagenous cells-cone cells and primary pigment cells. Cpr72Ec and Cpr66D, on the other hand, are more highly expressed in higher order interommatidial pigment cells. These data suggest that the complementary expression of cuticular proteins give rise to the center vs periphery of the corneal lens facet, possibly facilitating a refractive gradient that is known to reduce spherical aberration. Moreover, these studies provide a framework for future studies aimed at understanding the cuticular basis of corneal lens function in holometabolous insect eyes.

Revisiting Adaptive Potential, Population Size, and Conservation

Additive genetic variance (V_A) reflects the potential for evolutionary shifts and can be low for some traits or populations. High V_A is critical for the conservation of threatened species under selection to facilitate adaptation. Theory predicts tight associations between population size and V_A, but data from some experimental models, and managed and natural populations do not always support this prediction. However, V_A comparisons often have low statistical power, are undertaken in highly controlled environments distinct from natural habitats, and focus on traits with limited ecological relevance. Moreover, investigations of V_A typically fail to consider rare alleles, genetic load, or linkage disequilibrium, resulting in deleterious effects associated with favored alleles in small populations. Large population size remains essential for ensuring adaptation.

Genomics of Compensatory Adaptation in Experimental Populations of Aspergillus nidulans

Knowledge of the number and nature of genetic changes responsible for adaptation is essential for understanding and predicting evolutionary trajectories. Here, we study the genomic basis of compensatory adaptation to the fitness cost of fungicide resistance in experimentally evolved strains of the filamentous fungus Aspergillus nidulans The original selection experiment tracked the fitness recovery of lines founded by an ancestral strain that was resistant to fludioxonil, but paid a fitness cost in the absence of the fungicide. We obtained whole-genome sequence data for the ancestral A. nidulans strain and eight experimentally evolved strains. We find that fludioxonil resistance in the ancestor was likely conferred by a mutation in histidine kinase nikA, part of the two-component signal transduction system of the high-osmolarity glycerol (HOG) stress response pathway. To compensate for the pleiotropic negative effects of the resistance mutation, the subsequent fitness gains observed in the evolved lines were likely caused by secondary modification of HOG pathway activity. Candidate genes for the compensatory fitness increases were significantly overrepresented by stress response functions, and some were specifically associated with the HOG pathway itself. Parallel evolution at the gene level was rare among evolved lines. There was a positive relationship between the predicted number of adaptive steps, estimated from fitness data, and the number of genomic mutations, determined by whole-genome sequencing. However, the number of genomic mutations was, on average, 8.45 times greater than the number of adaptive steps inferred from fitness data. This research expands our understanding of the genetic basis of adaptation in multicellular eukaryotes and lays out a framework for future work on the genomics of compensatory adaptation in A. nidulans.

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.

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.