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.     

The genomic footprint of climate adaptation in Chironomus riparius

The gradual heterogeneity of climatic factors poses varying selection pressures across geographic distances that leave signatures of clinal variation in the genome. Separating signatures of clinal adaptation from signatures of other evolutionary forces, such as demographic processes, genetic drift and adaptation, to nonclinal conditions of the immediate local environment is a major challenge. Here, we examine climate adaptation in five natural populations of the harlequin fly Chironomus riparius sampled along a climatic gradient across Europe. Our study integrates experimental data, individual genome resequencing, Pool‐Seq data and population genetic modelling. Common‐garden experiments revealed significantly different population growth rates at test temperatures corresponding to the population origin along the climate gradient, suggesting thermal adaptation on the phenotypic level. Based on a population genomic analysis, we derived empirical estimates of historical demography and migration. We used an F_ST outlier approach to infer positive selection across the climate gradient, in combination with an environmental association analysis. In total, we identified 162 candidate genes as genomic basis of climate adaptation. Enriched functions among these candidate genes involved the apoptotic process and molecular response to heat, as well as functions identified in studies of climate adaptation in other insects. Our results show that local climate conditions impose strong selection pressures and lead to genomic adaptation despite strong gene flow. Moreover, these results imply that selection to different climatic conditions seems to converge on a functional level, at least between different insect species.     

Chironomus riparius (Diptera) genome sequencing reveals the impact of minisatellite transposable elements on population divergence

Active transposable elements (TEs) may result in divergent genomic insertion and abundance patterns among conspecific populations. Upon secondary contact, such divergent genetic backgrounds can theoretically give rise to classical Dobzhansky–Muller incompatibilities (DMI), thus contributing to the evolution of endogenous genetic barriers and eventually causing population divergence. We investigated differential TE abundance among conspecific populations of the nonbiting midge Chironomus riparius and evaluated their potential role in causing endogenous genetic incompatibilities between these populations. We focussed on a Chironomus‐specific TE, the minisatellite‐like Cla‐element, whose activity is associated with speciation in the genus. Using a newly generated and annotated draft genome for a genomic study with five natural C. riparius populations, we found highly population‐specific TE insertion patterns with many private insertions. A significant correlation of the pairwise F_ST estimated from genomewide single‐nucleotide polymorphisms (SNPs) and the F_ST estimated from TEs is consistent with drift as the major force driving TE population differentiation. However, the significantly higher Cla‐element F_ST level due to a high proportion of differentially fixed Cla‐element insertions also indicates selection against segregating (i.e. heterozygous) insertions. With reciprocal crossing experiments and fluorescent in situ hybridization of Cla‐elements to polytene chromosomes, we documented phenotypic effects on female fertility and chromosomal mispairings. We propose that the inferred negative selection on heterozygous Cla‐element insertions may cause endogenous genetic barriers and therefore acts as DMI among C. riparius populations. The intrinsic genomic turnover exerted by TEs may thus have a direct impact on population divergence that is operationally different from drift and local adaptation.     

Mechanisms of thermal adaptation and evolutionary potential of conspecific populations to changing environments

Heterogeneous and ever‐changing thermal environments drive the evolution of populations and species, especially when extreme conditions increase selection pressure for traits influencing fitness. However, projections of biological diversity under scenarios of climate change rarely consider evolutionary adaptive potential of natural species. In this study, we tested for mechanistic evidence of evolutionary thermal adaptation among ecologically divergent redband trout populations (Oncorhynchus mykiss gairdneri) in cardiorespiratory function, cellular response and genomic variation. In a common garden environment, fish from an extreme desert climate had significantly higher critical thermal maximum (p < .05) and broader optimum thermal window for aerobic scope (>3°C) than fish from cooler montane climate. In addition, the desert population had the highest maximum heart rate during warming (20% greater than montane populations), indicating improved capacity to deliver oxygen to internal tissues. In response to acute heat stress, distinct sets of cardiac genes were induced among ecotypes, which helps to explain the differences in cardiorespiratory function. Candidate genomic markers and genes underlying these physiological adaptations were also pinpointed, such as genes involved in stress response and metabolic activity (hsp40, ldh‐b and camkk2). These markers were developed into a multivariate model that not only accurately predicted critical thermal maxima, but also evolutionary limit of thermal adaptation in these specific redband trout populations relative to the expected limit for the species. This study demonstrates mechanisms and limitations of an aquatic species to evolve under changing environments that can be incorporated into advanced models to predict ecological consequences of climate change for natural organisms.     

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new‐type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new‐type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population‐based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.     

Rapid responses to a strong experimental selection for heat hardening in the invasive whitefly Bemisia tabaci MEAM 1

Adaptation to new environments is an important issue for invasive species as colonization depends on evolvability in their new distribution range. Here, we considered the case of the whitefly Bemisia tabaci MEAM 1 (Gennadius) (Hemiptera: Aleyrodidae), a pest that has recently invaded Colombia and where thermal adaptation has been proposed to explain its colonizing ability. An experimental evolution study was conducted to assess the evolutionary potential of B. tabaci in relation to its upper thermal limits, to explain its rapid adaptation during post‐invasion periods. Selection for hardening capacity was conducted in four whitefly populations. We measured thermal responses in relation to fitness components (survival, fecundity, and viability) for 5–7 generations under a strong selection regime. Heat hardening responded rapidly in both sexes. This was expressed as an increase in survival, but not in fecundity or viability. These results suggest that thermal responses for heat hardening are not correlated and evolve independently. Increased survival after few generations of selection points to high adaptive potential in this insect, which leads to rapid post‐invasion adaptation. Our study can help to predict population responses to environmental change and explain the colonizing ability of this pest.     

Neutral and adaptive genomic signatures of rapid poleward range expansion

Many species are expanding their range polewards, and this has been associated with rapid phenotypic change. Yet, it is unclear to what extent this reflects rapid genetic adaptation or neutral processes associated with range expansion, or selection linked to the new thermal conditions encountered. To disentangle these alternatives, we studied the genomic signature of range expansion in the damselfly Coenagrion scitulum using 4950 newly developed genomic SNPs and linked this to the rapidly evolved phenotypic differences between core and (newly established) edge populations. Most edge populations were genetically clearly differentiated from the core populations and all were differentiated from each other indicating independent range expansion events. In addition, evidence for genetic drift in the edge populations, and strong evidence for adaptive genetic variation in association with the range expansion was detected. We identified one SNP under consistent selection in four of the five edge populations and showed that the allele increasing in frequency is associated with increased flight performance. This indicates collateral, non‐neutral evolutionary changes in independent edge populations driven by the range expansion process. We also detected a genomic signature of adaptation to the newly encountered thermal regimes, reflecting a pattern of countergradient variation. The latter signature was identified at a single SNP as well as in a set of covarying SNPs using a polygenic multilocus approach to detect selection. Overall, this study highlights how a strategic geographic sampling design and the integration of genomic, phenotypic and environmental data can identify and disentangle the neutral and adaptive processes that are simultaneously operating during range expansions.     

Dynamics of genomic change during evolutionary rescue in the seed beetle Callosobruchus maculatus

Rapid adaptation can prevent extinction when populations are exposed to extremely marginal or stressful environments. Factors that affect the likelihood of evolutionary rescue from extinction have been identified, but much less is known about the evolutionary dynamics (e.g., rates and patterns of allele frequency change) and genomic basis of successful rescue, particularly in multicellular organisms. We conducted an evolve‐and‐resequence experiment to investigate the dynamics of evolutionary rescue at the genetic level in the cowpea seed beetle, Callosobruchus maculatus, when it is experimentally shifted to a stressful host plant, lentil. Low survival (~1%) at the onset of the experiment caused population decline. But adaptive evolution quickly rescued the population, with survival rates climbing to 69% by the F5 generation and 90% by the F10 generation. Population genomic data showed that rescue likely was caused by rapid evolutionary change at multiple loci, with many alleles fixing or nearly fixing within five generations of selection on lentil. Selection on these loci was only moderately consistent in time, but parallel evolutionary changes were evident in sublines formed after the lentil line had passed through a bottleneck. By comparing estimates of selection and genomic change on lentil across five independent C. maculatus lines (the new lentil‐adapted line, three long‐established lines and one case of failed evolutionary rescue), we found that adaptation on lentil occurred via somewhat idiosyncratic evolutionary changes. Overall, our results suggest that evolutionary rescue in this system can be caused by very strong selection on multiple loci driving rapid and pronounced genomic change.     

RAPID LABORATORY EVOLUTION OF ADULT LIFE-HISTORY TRAITS IN DROSOPHILA MELANOGASTER IN RESPONSE TO TEMPERATURE

Three replicate lines of Drosophila melanogaster were cultured at each of two temperatures (16.5°C and 25°C) in population cages for 4 yr. The lifespans of both sexes and the fecundity and fertility of the females were then measured at both experimental temperatures. The characters showed evidence of adaptation; flies of both sexes from each selection regime showed higher longevity, and females showed higher fecundity and fertility, than flies from the other selection regime when they were tested at the experimental temperature at which they had evolved. Calculation of intrinsic rates of increase under different assumptions about the rate of population increase showed that the difference between the lines from the two selection regimes became less the higher the rate of population increase, because the lines were more similar in early adulthood than they were later. Despite the increased adaptation of the low-temperature lines to the low temperature, like the high temperature lines they produced progeny at a higher rate at the higher temperature. The lines may have independently evolved adaptations to their respective thermal regimes during the experiment, or there may have been a trade-off between adaptation to the two temperatures, or mutation pressure may have lowered adaptation to the temperature that the flies no longer encountered.     

Rapid,nonparallel genomic evolution of Brassica rapa (field mustard) under experimental drought

While we know that climate change can potentially cause rapid phenotypic evolution, our understanding of the genetic basis and degree of genetic parallelism of rapid evolutionary responses to climate change is limited. In this study, we combined the resurrection approach with an evolve-and-resequence design to examine genome-wide evolutionary changes following drought. We exposed genetically similar replicate populations of the annual plant Brassica rapa derived from a field population in southern California to four generations of experimental drought or watered conditions in a greenhouse. Genome-wide sequencing of ancestral and descendant population pools identified hundreds of SNPs that showed evidence of rapidly evolving in response to drought. Several of these were in stress response genes, and two were identified in a prior study of drought response in this species. However, almost all genetic changes were unique among experimental populations, indicating that the evolutionary changes were largely nonparallel, despite the fact that genetically similar replicates of the same founder population had experienced controlled and consistent selection regimes. This nonparallelism of evolution at the genetic level is potentially because of polygenetic adaptation allowing for multiple different genetic routes to similar phenotypic outcomes. Our findings help to elucidate the relationship between rapid phenotypic and genomic evolution and shed light on the degree of parallelism and predictability of genomic evolution to environmental change.     

EVOLUTIONARY ADAPTATION TO TEMPERATURE. V. ADAPTIVE MECHANISMS AND CORRELATED RESPONSES IN EXPERIMENTAL LINES OF ESCHERICHIA COLI

We previously demonstrated temperature-specific genetic adaptation in experimental lines of Escherichia coli. Six initially identical populations were propagated for 2000 generations under each of five regimes: constant 20°C, 32°C, 37°C, and 42°C, and a daily switch between 32°C and 42°C. Glucose was the sole carbon source in all cases. Here, we examine the physiological bases of adaptation to determine whether the same mechanisms evolved among the replicate lines within each thermal regime and across different regimes. Specifically, we investigate whether changes in glucose transport may account for the temperature-specific adaptation. We compared each line's direct response of fitness to glucose with its correlated response to maltose; glucose and maltose enter the cell by different pathways, but their catabolism is identical. Except for lines maintained at the ancestral temperature (37°C), almost all derived lines had significantly different fitnesses (relative to their common ancestor) in glucose and maltose, supporting the hypothesis that adaptation involved changes in glucose transport. An alternative explanation, that maltose transport decayed by genetic drift, appears unlikely for reasons that are discussed. Although most lines showed evidence of temperature-specific adaptation in glucose transport, several different mechanisms may underlie these improvements, as indicated by heterogeneity in correlated responses (across temperatures and substrates) among replicate lines adapted to the same regime. This heterogeneity provides a latent pool of genetic variation for responding to environmental change.     

Evolutionary determinants of population differences in population growth rate × habitat temperature interactions in Chironomus riparius

Little is known about intraspecific variation in fitness performance in response to thermal stress among natural populations and how this relates to evolutionary aspects of species ecology. In this study, population growth rate (PGR; a composite fitness measure) varied among five natural Chironomus riparius populations sampled across a climatic gradient when subjected to three temperature treatments reflecting the typical range of summer habitat temperatures (20, 24 and 28 °C). The variation could be explained by a complex model including effects of genetic drift, genetic diversity and adaptation to average temperature during the warmest month, in addition to experimental temperature. All populations suffered a decrease in PGR from 20 to 28 °C and ΔPGR was significantly correlated with the respective average habitat temperature in the warmest month—populations from warmer areas showing lower ΔPGR. This implies that long-term exposure to higher temperatures in the warmest month (the key reproductive period for C. riparius) is likely to be a key selective force influencing fitness at higher temperatures. A comparison of phenotypic divergence and neutral genetic differentiation revealed that one phenotypic trait—the number of fertile egg masses per female—appeared to be under positive selection in some populations. Our findings support a role for response to temperature selection along a climatic gradient and suggest population history is a key determinant of intraspecific fitness variation. We stress the importance of integrating different types of data (climatic, experimental, genetic) in order to understand the effects of global climate change on biodiversity.     

Adaptive divergence in embryonic thermal plasticity among Atlantic salmon populations

In the context of global changes, the long‐term viability of populations of endangered ectotherms may depend on their adaptive potential and ability to cope with temperature variations. We measured responses of Atlantic salmon embryos from four populations to temperature variations and used a Q_ST–F_ST approach to study the adaptive divergence among these populations. Embryos were reared under two experimental conditions: a low temperature regime at 4 °C until eyed‐stage and 10 °C until the end of embryonic development and a high temperature regime with a constant temperature of 10 °C throughout embryonic development. Significant variations among populations and population × temperature interactions were observed for embryo survival, incubation time and length. Q_ST was higher than F_ST in all but one comparison suggesting an important effect of divergent selection. Q_ST was also higher under the high‐temperature treatment than at low temperature for length and survival due to a higher variance among populations under the stressful warmer treatment. Interestingly, heritability was lower for survival under high temperature in relation to a lower additive genetic variance under that treatment. Overall, these results reveal an adaptive divergence in thermal plasticity in embryonic life stages of Atlantic salmon suggesting that salmon populations may differentially respond to temperature variations induced by climate change. These results also suggest that changes in temperature may alter not only the adaptive potential of natural populations but also the selection regimes among them.     

Thermal adaptation and phenotypic plasticity in a warming world: Insights from common garden experiments on Alaskan sockeye salmon

An important unresolved question is how populations of coldwater‐dependent fishes will respond to rapidly warming water temperatures. For example, the culturally and economically important group, Pacific salmon (Oncorhynchus spp.), experience site‐specific thermal regimes during early development that could be disrupted by warming. To test for thermal local adaptation and heritable phenotypic plasticity in Pacific salmon embryos, we measured the developmental rate, survival, and body size at hatching in two populations of sockeye salmon (Oncorhynchus nerka) that overlap in timing of spawning but incubate in contrasting natural thermal regimes. Using a split half‐sibling design, we exposed embryos of 10 families from each of two populations to variable and constant thermal regimes. These represented both experienced temperatures by each population, and predicted temperatures under plausible future conditions based on a warming scenario from the downscaled global climate model (MIROC A1B scenario). We did not find evidence of thermal local adaptation during the embryonic stage for developmental rate or survival. Within treatments, populations hatched within 1 day of each other, on average, and among treatments, did not differ in survival in response to temperature. We did detect plasticity to temperature; embryos developed 2.5 times longer (189 days) in the coolest regime compared to the warmest regime (74 days). We also detected variation in developmental rates among families within and among temperature regimes, indicating heritable plasticity. Families exhibited a strong positive relationship between thermal variability and phenotypic variability in developmental rate but body length and mass at hatching were largely insensitive to temperature. Overall, our results indicated a lack of thermal local adaptation, but a presence of plasticity in populations experiencing contrasting conditions, as well as family‐specific heritable plasticity that could facilitate adaptive change.     

The importance of pre-and postcopulatory sexual selection promoting adaptation to increasing temperatures

Global temperatures are increasing rapidly affecting species globally. Understanding if and how different species can adapt fast enough to keep up with increasing temperatures is of vital importance. One mechanism that can accelerate adaptation and promote evolutionary rescue is sexual selection. Two different mechanisms by which sexual selection can facilitate adaptation are pre- and postcopulatory sexual selection. However, the relative effects of these different forms of sexual selection in promoting adaptation are unknown. Here, we present the results from an experimental study in which we exposed fruit flies Drosophila melanogaster to either no mate choice or 1 of 2 different sexual selection regimes (pre- and postcopulatory sexual selection) for 6 generations, under different thermal regimes. Populations showed evidence of thermal adaptation under precopulatory sexual selection, but this effect was not detected in the postcopulatory sexual selection and the no choice mating regime. We further demonstrate that sexual dimorphism decreased when flies evolved under increasing temperatures, consistent with recent theory predicting more sexually concordant selection under environmental stress. Our results suggest an important role for precopulatory sexual selection in promoting thermal adaptation and evolutionary rescue.     

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.     

Experimental evolution for generalists and specialists reveals multivariate genetic constraints on thermal reaction norms

Theory predicts the emergence of generalists in variable environments and antagonistic pleiotropy to favour specialists in constant environments, but empirical data seldom support such generalist–specialist trade‐offs. We selected for generalists and specialists in the dung fly Sepsis punctum (Diptera: Sepsidae) under conditions that we predicted would reveal antagonistic pleiotropy and multivariate trade‐offs underlying thermal reaction norms for juvenile development. We performed replicated laboratory evolution using four treatments: adaptation at a hot (31 °C) or a cold (15 °C) temperature, or under regimes fluctuating between these temperatures, either within or between generations. After 20 generations, we assessed parental effects and genetic responses of thermal reaction norms for three correlated life‐history traits: size at maturity, juvenile growth rate and juvenile survival. We find evidence for antagonistic pleiotropy for performance at hot and cold temperatures, and a temperature‐mediated trade‐off between juvenile survival and size at maturity, suggesting that trade‐offs associated with environmental tolerance can arise via intensified evolutionary compromises between genetically correlated traits. However, despite this antagonistic pleiotropy, we found no support for the evolution of increased thermal tolerance breadth at the expense of reduced maximal performance, suggesting low genetic variance in the generalist–specialist dimension.     

Increases in the evolutionary potential of upper thermal limits under warmer temperatures in two rainforest Drosophila species

Tropical and subtropical species represent the majority of biodiversity. These species are predicted to lack the capacity to evolve higher thermal limits in response to selection imposed by climatic change. However, these assessments have relied on indirect estimates of adaptive capacity, using conditions that do not reflect environmental changes projected under climate change. Using a paternal half‐sib full‐sib breeding design, we estimated the additive genetic variance and narrow‐sense heritability for adult upper thermal limits in two rainforest‐restricted species of Drosophila reared under two thermal regimes, reflecting increases in seasonal temperature projected for the Wet Tropics of Australia and under standard laboratory conditions (constant 25°C). Estimates of additive genetic variation and narrow‐sense heritability for adult heat tolerance were significantly different from zero in both species under projected summer, but not winter or constant, thermal regimes. In contrast, significant broad‐sense genetic variation was apparent in all thermal regimes for egg‐to‐adult viability. Environment‐dependent changes in the expression of genetic variation for adult upper thermal limits suggest that predicting adaptive responses to climate change will be difficult. Estimating adaptive capacity under conditions that do not reflect future environmental conditions may provide limited insight into evolutionary responses to climate change.     

Signatures of local adaptation along environmental gradients in a range‐expanding damselfly (Ischnura elegans)

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range‐expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site‐associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V‐ATPase) and visual processes (long‐wavelength‐sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate‐mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.