Population genetic evidence for cold adaptation in European Drosophila melanogaster populations

Polygenic adaptation contributes to the invasive success of the Colorado potato beetle

How invasive species cope with novel selective pressures with limited genetic variation is a fundamental question in molecular ecology. Several mechanisms have been proposed, but they can lack generality. Here, we addressed an alternative solution, polygenic adaptation, wherein traits that arise from multiple combinations of loci may be less sensitive to loss of variation during invasion. We tested the polygenic signal of environmental adaptation of Colorado potato beetle (CPB) introduced in Eurasia. Population genomic analyses showed declining genetic diversity in the eastward expansion of Eurasian populations, and weak population genetic structure (except for the invasion fronts in Asia). Demographic history showed that all populations shared a strong bottleneck about 100 years ago when CPB was introduced to Europe. Genome scans revealed a suite of genes involved in activity regulation functions that are plausibly related to cold stress, including some well-founded functions (e.g., the activity of phosphodiesterase, the G-protein regulator) and discrete functions. Such polygenic architecture supports the hypothesis that polygenic adaptation and potentially genetic redundancy can fuel the adaptation of CPB despite strong genetic depletion, thus representing a promising general mechanism for resolving the genetic paradox of invasion. More broadly, most complex traits based on polygenes may be less sensitive to invasive bottlenecks, thus ensuring the evolutionary success of invasive species in novel environments.

Cardiac performance in heat-stressed flies of heat-susceptible and heat-resistant Drosophila melanogaster

Thermotolerance is a complex trait that can greatly differ between heat-susceptible (HS) and heat-adapted populations of small insects including Drosophila, with short-term effects after a sub-lethal level of heat stress on many physiological functions. Cardiac performance could accordingly be more robust in heat-resistant (HR) than in HS individuals under heat stress. Here, we tested heart performance under heat-stress effects in two recombinant inbred lines (RIL) of Drosophila melanogaster that dramatically differ in heat knockdown resistance. Heart rate did not strongly differ between heat-susceptible and heat-tolerant flies after a sub-lethal heat stress. Instead, heat-susceptible flies showed a much higher arrhythmia incidence, a longer duration of each heartbeat, and a larger amount of bradycardia than heat-tolerant flies. The highly conserved cardiac proteins SERCA, RyR and NCX that participate in the excitation/contraction coupling, did not differ in activity level between HR and HS flies. Available information for both RIL suggests that heart performance under heat stress may be linked, at least partially, to candidate genes of previously identified quantitative trait loci (QTL) for thermotolerance. This study indicates that HR flies can be genetically more robust in their heart performance than HS flies under even sub-lethal levels of heat stress.

The effects of genetic drift and genomic selection on differentiation and local adaptation of the introduced populations of Aedes albopictus in southern Russia

Background

Asian tiger mosquito Aedes albopictus is an arbovirus vector that has spread from its native habitation areal in Southeast Asia throughout North and South Americas, Europe, and Africa. Ae. albopictus was first detected in the Southern Federal District of the Russian Federation in the subtropical town of Sochi in 2011. In subsequent years, this species has been described in the continental areas with more severe climate and lower winter temperatures.

Methods

Genomic analysis of pooled Ae. albopictus samples collected in the mosquito populations in the coastal and continental regions of the Krasnodar Krai was conducted to look for the genetic changes associated with the spread and potential cold adaptation in Ae. albopictus.

Results

The results of the phylogenetic analysis based on mitochondrial genomes corresponded well with the hypothesis that Ae. albopictus haplotype A1a2a1 was introduced into the region from a single source. Population analysis revealed the role of dispersal and genetic drift in the local adaptation of the Asian tiger mosquito. The absence of shared haplotypes between the samples and high fixation indices suggest that gene flow between samples was heavily restricted. Mitochondrial and genomic differentiation together with different distances between dispersal routes, natural and anthropogenic barriers and local effective population size reduction could lead to difficulties in local climatic adaptations due to reduced selection effectiveness. We have found genomic regions with selective sweep patterns which can be considered as having been affected by recent selection events. The genes located in these regions participate in neural protection, lipid conservation, and cuticle formation during diapause. These processes were shown to be important for cold adaptation in the previous transcriptomic and proteomic studies. However, the population history and relatively low coverage obtained in the present article could have negatively affect sweep detection.

Parallel and Population-specific Gene Regulatory Evolution in Cold-Adapted Fly Populations

Changes in gene regulation at multiple levels may comprise an important share of the molecular changes underlying adaptive evolution in nature. However, few studies have assayed within- and between-population variation in gene regulatory traits at a transcriptomic scale, and therefore inferences about the characteristics of adaptive regulatory changes have been elusive. Here, we assess quantitative trait differentiation in gene expression levels and alternative splicing (intron usage) between three closely related pairs of natural populations of Drosophila melanogaster from contrasting thermal environments that reflect three separate instances of cold tolerance evolution. The cold-adapted populations were known to show population genetic evidence for parallel evolution at the SNP level, and here we find evidence for parallel expression evolution between them, with stronger parallelism at larval and adult stages than for pupae. We also implement a flexible method to estimate cis- vs trans-encoded contributions to expression or splicing differences at the adult stage. The apparent contributions of cis- vs trans-regulation to adaptive evolution vary substantially among population pairs. While two of three population pairs show a greater enrichment of cis-regulatory differences among adaptation candidates, trans-regulatory differences are more likely to be implicated in parallel expression changes between population pairs. Genes with significant cis-effects are enriched for signals of elevated genetic differentiation between cold- and warm-adapted populations, suggesting that they are potential targets of local adaptation. These findings expand our knowledge of adaptive gene regulatory evolution and our ability to make inferences about this important and widespread process.

False and true positives in arthropod thermal adaptation candidate gene lists

Genome-wide studies are prone to false positives due to inherently low priors and statistical power. One approach to ameliorate this problem is to seek validation of reported candidate genes across independent studies: genes with repeatedly discovered effects are less likely to be false positives. Inversely, genes reported only as many times as expected by chance alone, while possibly representing novel discoveries, are also more likely to be false positives. We show that, across over 30 genome-wide studies that reported Drosophila and Daphnia genes with possible roles in thermal adaptation, the combined lists of candidate genes and orthologous groups are rapidly approaching the total number of genes and orthologous groups in the respective genomes. This is consistent with the expectation of high frequency of false positives. The majority of these spurious candidates have been identified by one or a few studies, as expected by chance alone. In contrast, a noticeable minority of genes have been identified by numerous studies with the probabilities of such discoveries occurring by chance alone being exceedingly small. For this subset of genes, different studies are in agreement with each other despite differences in the ecological settings, genomic tools and methodology, and reporting thresholds. We provide a reference set of presumed true positives among Drosophila candidate genes and orthologous groups involved in response to changes in temperature, suitable for cross-validation purposes. Despite this approach being prone to false negatives, this list of presumed true positives includes several hundred genes, consistent with the "omnigenic" concept of genetic architecture of complex traits.

Temperature, rainfall and wind variables underlie environmental adaptation in natural populations of Drosophila melanogaster

While several studies in a diverse set of species have shed light on the genes underlying adaptation, our knowledge on the selective pressures that explain the observed patterns lags behind. Drosophila melanogaster is a valuable organism to study environmental adaptation because this species originated in Southern Africa and has recently expanded worldwide, and also because it has a functionally well-annotated genome. In this study, we aimed to decipher which environmental variables are relevant for adaptation of D. melanogaster natural populations in Europe and North America. We analysed 36 whole-genome pool-seq samples of D. melanogaster natural populations collected in 20 European and 11 North American locations. We used the BayPass software to identify single nucleotide polymorphisms (SNPs) and transposable elements (TEs) showing signature of adaptive differentiation across populations, as well as significant associations with 59 environmental variables related to temperature, rainfall, evaporation, solar radiation, wind, daylight hours, and soil type. We found that in addition to temperature and rainfall, wind related variables are also relevant for D. melanogaster environmental adaptation. Interestingly, 23%-51% of the genes that showed significant associations with environmental variables were not found overly differentiated across populations. In addition to SNPs, we also identified 10 reference transposable element insertions associated with environmental variables. Our results showed that genome-environment association analysis can identify adaptive genetic variants that are undetected by population differentiation analysis while also allowing the identification of candidate environmental drivers of adaptation.

Genomic Analysis of European Drosophila melanogaster Populations Reveals Longitudinal Structure, Continent-Wide Selection, and Previously Unknown DNA Viruses

Genetic variation is the fuel of evolution, with standing genetic variation especially important for short-term evolution and local adaptation. To date, studies of spatiotemporal patterns of genetic variation in natural populations have been challenging, as comprehensive sampling is logistically difficult, and sequencing of entire populations costly. Here, we address these issues using a collaborative approach, sequencing 48 pooled population samples from 32 locations, and perform the first continent-wide genomic analysis of genetic variation in European Drosophila melanogaster. Our analyses uncover longitudinal population structure, provide evidence for continent-wide selective sweeps, identify candidate genes for local climate adaptation, and document clines in chromosomal inversion and transposable element frequencies. We also characterize variation among populations in the composition of the fly microbiome, and identify five new DNA viruses in our samples.

Functional Analysis of a Putative Target of Spatially Varying Selection in the Menin1 Gene of Drosophila melanogaster

While significant effort has been devoted to investigating the potential influence of spatially varying selection on genomic variation, relatively little effort has been devoted to experimental analysis of putative variants or genes experiencing such selection. Previous population genetic work identified an amino acid polymorphism in the Mnn1 gene as one of the most strongly latitudinally differentiated SNPs in the genome of Drosophila melanogaster in the United States and Australia. Here we report the results of our transgenic analysis of this amino acid polymorphism. Genotypes carrying alternative Mnn1 alleles differed in multiple phenotypes in a direction generally consistent with phenotypic differences previously observed along latitudinal clines. These results support inferences from earlier population genomic work that this variant influences fitness, and support the idea that the alleles exhibiting clines may be likely to have pleiotropic effects that are correlated along the axes favored by natural selection.

Selection for reproduction under short photoperiods changes diapause-associated traits and induces widespread genomic divergence

The incidence of reproductive diapause is a critical aspect of life history in overwintering insects from temperate regions. Much has been learned about the timing, physiology and genetics of diapause in a range of insects, but how the multiple changes involved in this and other photoperiodically regulated traits are interrelated is not well understood. We performed quasinatural selection on reproduction under short photoperiods in a northern fly species, Drosophila montana, to trace the effects of photoperiodic selection on traits regulated by the photoperiodic timer and / or by a circadian clock system. Selection changed several traits associated with reproductive diapause, including the critical day length for diapause (CDL), the frequency of diapausing females under photoperiods that deviate from daily 24 h cycles and cold tolerance, towards the phenotypes typical of lower latitudes. However, selection had no effect on the period of free-running locomotor activity rhythm regulated by the circadian clock in fly brain. At a genomic level, selection induced extensive divergence between the selection and control line replicates in 16 gene clusters involved in signal transduction, membrane properties, immunologlobulins and development. These changes resembled ones detected between latitudinally divergent D. montana populations in the wild and involved SNP divergence associated with several genes linked with diapause induction. Overall, our study shows that photoperiodic selection for reproduction under short photoperiods affects diapause-associated traits without disrupting the central clock network generating circadian rhythms in fly locomor activity.

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.

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.

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.

Geographical analysis of diapause inducibility in European Drosophila melanogaster populations

Seasonal overwintering in insects represents an adaptation to stressful environments and in European Drosophila melanogaster females, low temperatures and short photoperiods can induce an ovarian diapause. Diapause may represent a recent (<15Ky) adaptation to the colonisation of temperate Europe by D. melanogaster from tropical sub-Saharan Africa, because African D. melanogaster and the sibling species D. simulans, have been reported to fail to undergo diapause. Over the past few centuries, D. melanogaster have also invaded North America and Australia, and eastern populations on both continents show a predictable latitudinal cline in diapause induction. In Europe however, a new diapause-enhancing timeless allele, ls-tim, is observed at high levels in southern Italy (∼80%), where it appears to have arisen and has spread throughout the continent with a frequency of ∼20% in Scandinavia. Given the phenotype of ls-tim and its geographical distribution, we might predict that it would work against any latitudinal cline in diapause induction within Europe. Indeed we reveal that any latitudinal cline for diapause in Europe is very weak, as predicted by ls-tim frequencies. In contrast, we determine ls-tim frequencies in North America and observe that they would be expected to strengthen the latitudinal pattern of diapause. Our results reveal how a newly arisen mutation, can, via the stochastic nature of where it initially arose, blur an otherwise adaptive geographical pattern.

α-amanitin resistance in Drosophila melanogaster: A genome-wide association approach

We investigated the mechanisms of mushroom toxin resistance in the Drosophila Genetic Reference Panel (DGRP) fly lines, using genome-wide association studies (GWAS). While Drosophila melanogaster avoids mushrooms in nature, some lines are surprisingly resistant to α-amanitin—a toxin found solely in mushrooms. This resistance may represent a pre-adaptation, which might enable this species to invade the mushroom niche in the future. Although our previous microarray study had strongly suggested that pesticide-metabolizing detoxification genes confer α-amanitin resistance in a Taiwanese D. melanogaster line Ama-KTT, none of the traditional detoxification genes were among the top candidate genes resulting from the GWAS in the current study. Instead, we identified Megalin, Tequila, and widerborst as candidate genes underlying the α-amanitin resistance phenotype in the North American DGRP lines, all three of which are connected to the Target of Rapamycin (TOR) pathway. Both widerborst and Tequila are upstream regulators of TOR, and TOR is a key regulator of autophagy and Megalin-mediated endocytosis. We suggest that endocytosis and autophagy of α-amanitin, followed by lysosomal degradation of the toxin, is one of the mechanisms that confer α-amanitin resistance in the DGRP lines.

Parallel Evolution of Cold Tolerance within Drosophila melanogaster

Drosophila melanogaster originated in tropical Africa before expanding into strikingly different temperate climates in Eurasia and beyond. Here, we find elevated cold tolerance in three distinct geographic regions: beyond the well-studied non-African case, we show that populations from the highlands of Ethiopia and South Africa have significantly increased cold tolerance as well. We observe greater cold tolerance in outbred versus inbred flies, but only in populations with higher inversion frequencies. Each cold-adapted population shows lower inversion frequencies than a closely-related warm-adapted population, suggesting that inversion frequencies may decrease with altitude in addition to latitude. Using the FST-based "Population Branch Excess" statistic (PBE), we found only limited evidence for parallel genetic differentiation at the scale of ∼4 kb windows, specifically between Ethiopian and South African cold-adapted populations. And yet, when we looked for single nucleotide polymorphisms (SNPs) with codirectional frequency change in two or three cold-adapted populations, strong genomic enrichments were observed from all comparisons. These findings could reflect an important role for selection on standing genetic variation leading to "soft sweeps". One SNP showed sufficient codirectional frequency change in all cold-adapted populations to achieve experiment-wide significance: an intronic variant in the synaptic gene Prosap. Another codirectional outlier SNP, at senseless-2, had a strong association with our cold trait measurements, but in the opposite direction as predicted. More generally, proteins involved in neurotransmission were enriched as potential targets of parallel adaptation. The ability to study cold tolerance evolution in a parallel framework will enhance this classic study system for climate adaptation.

A Variable Genetic Architecture of Melanic Evolution in Drosophila melanogaster

Unraveling the genetic architecture of adaptive phenotypic divergence is a fundamental quest in evolutionary biology. In Drosophila melanogaster, high-altitude melanism has evolved in separate mountain ranges in sub-Saharan Africa, potentially as an adaptation to UV intensity. We investigated the genetic basis of this melanism in three populations using a new bulk segregant analysis mapping method. We identified 19 distinct QTL regions from nine mapping crosses, with several QTL peaks overlapping between two or all populations, and yet different crosses involving the same melanic population commonly yielded distinct QTL. The strongest QTL often overlapped well-known pigmentation genes, but we typically did not find wide signals of genetic differentiation (F_ST) between lightly and darkly pigmented populations at these genes. Instead, we found small numbers of highly differentiated SNPs at the probable causative genes. A simulation analysis showed that these patterns of polymorphism were consistent with selection on standing genetic variation. Overall, our results suggest that, even for potentially simpler traits like pigmentation, the complexity of adaptive trait evolution poses important challenges for QTL mapping and population genetic analysis.

Genomics of clinal variation in Drosophila: disentangling the interactions of selection and demography

Clines in phenotypes and genotype frequencies across environmental gradients are commonly taken as evidence for spatially varying selection. Classical examples include the latitudinal clines in various species of Drosophila, which often occur in parallel fashion on multiple continents. Today, genomewide analysis of such clinal systems provides a fantastic opportunity for unravelling the genetics of adaptation, yet major challenges remain. A well-known but often neglected problem is that demographic processes can also generate clinality, independent of or coincident with selection. A closely related issue is how to identify true genic targets of clinal selection. In this issue of Molecular Ecology, three studies illustrate these challenges and how they might be met. Bergland et al. report evidence suggesting that the well-known parallel latitudinal clines in North American and Australian D. melanogaster are confounded by admixture from Africa and Europe, highlighting the importance of distinguishing demographic from adaptive clines. In a companion study, Machado et al. provide the first genomic comparison of latitudinal differentiation in D. melanogaster and its sister species D. simulans. While D. simulans is less clinal than D. melanogaster, a significant fraction of clinal genes is shared between both species, suggesting the existence of convergent adaptation to clinaly varying selection pressures. Finally, by drawing on several independent sources of evidence, Božičević et al. identify a functional network of eight clinal genes that are likely involved in cold adaptation. Together, these studies remind us that clinality does not necessarily imply selection and that separating adaptive signal from demographic noise requires great effort and care.