A 2-million-year-old ecosystem in Greenland uncovered by environmental DNA.

Late Pliocene and Early Pleistocene epochs 3.6 to 0.8 million years ago1 had climates resembling those forecasted under future warming2. Palaeoclimatic records show strong polar amplification with mean annual temperatures of 11-19 °C above contemporary values3,4. The biological communities inhabiting the Arctic during this time remain poorly known because fossils are rare5. Here we report an ancient environmental DNA6 (eDNA) record describing the rich plant and animal assemblages of the Kap København Formation in North Greenland, dated to around two million years ago. The record shows an open boreal forest ecosystem with mixed vegetation of poplar, birch and thuja trees, as well as a variety of Arctic and boreal shrubs and herbs, many of which had not previously been detected at the site from macrofossil and pollen records. The DNA record confirms the presence of hare and mitochondrial DNA from animals including mastodons, reindeer, rodents and geese, all ancestral to their present-day and late Pleistocene relatives. The presence of marine species including horseshoe crab and green algae support a warmer climate than today. The reconstructed ecosystem has no modern analogue. The survival of such ancient eDNA probably relates to its binding to mineral surfaces. Our findings open new areas of genetic research, demonstrating that it is possible to track the ecology and evolution of biological communities from two million years ago using ancient eDNA.

Emergence of methicillin resistance predates the clinical use of antibiotics.

The discovery of antibiotics more than 80 years ago has led to considerable improvements in human and animal health. Although antibiotic resistance in environmental bacteria is ancient, resistance in human pathogens is thought to be a modern phenomenon that is driven by the clinical use of antibiotics1. Here we show that particular lineages of methicillin-resistant Staphylococcus aureus-a notorious human pathogen-appeared in European hedgehogs in the pre-antibiotic era. Subsequently, these lineages spread within the local hedgehog populations and between hedgehogs and secondary hosts, including livestock and humans. We also demonstrate that the hedgehog dermatophyte Trichophyton erinacei produces two ß-lactam antibiotics that provide a natural selective environment in which methicillin-resistant S. aureus isolates have an advantage over susceptible isolates. Together, these results suggest that methicillin resistance emerged in the pre-antibiotic era as a co-evolutionary adaptation of S. aureus to the colonization of dermatophyte-infected hedgehogs. The evolution of clinically relevant antibiotic-resistance genes in wild animals and the connectivity of natural, agricultural and human ecosystems demonstrate that the use of a One Health approach is critical for our understanding and management of antibiotic resistance, which is one of the biggest threats to global health, food security and development.

A single introduction of wild rabbits triggered the biological invasion of Australia.

Biological invasions are a major cause of environmental and economic disruption. While ecological factors are key determinants of their success, the role of genetics has been more challenging to demonstrate. The colonization of Australia by the European rabbit is one of the most iconic and devastating biological invasions in recorded history. Here, we show that despite numerous introductions over a 70-y period, this invasion was triggered by a single release of a few animals that spread thousands of kilometers across the continent. We found genetic support for historical accounts that these were English rabbits imported in 1859 by a settler named Thomas Austin and traced the origin of the invasive population back to his birthplace in England. We also find evidence of additional introductions that established local populations but have not spread geographically. Combining genomic and historical data we show that, contrary to the earlier introductions, which consisted mostly of domestic animals, the invasive rabbits had wild ancestry. In New Zealand and Tasmania, rabbits also became a pest several decades after being introduced. We argue that the common denominator of these invasions was the arrival of a new genotype that was better adapted to the natural environment. These findings demonstrate how the genetic composition of invasive individuals can determine the success of an introduction and provide a mechanism by which multiple introductions can be required for a biological invasion.

Life and Death of Selfish Genes: Comparative Genomics Reveals the Dynamic Evolution of Cytoplasmic Incompatibility.

Cytoplasmic incompatibility is a selfish reproductive manipulation induced by the endosymbiont Wolbachia in arthropods. In males Wolbachia modifies sperm, leading to embryonic mortality in crosses with Wolbachia-free females. In females, Wolbachia rescues the cross and allows development to proceed normally. This provides a reproductive advantage to infected females, allowing the maternally transmitted symbiont to spread rapidly through host populations. We identified homologs of the genes underlying this phenotype, cifA and cifB, in 52 of 71 new and published Wolbachia genome sequences. They are strongly associated with cytoplasmic incompatibility. There are up to seven copies of the genes in each genome, and phylogenetic analysis shows that Wolbachia frequently acquires new copies due to pervasive horizontal transfer between strains. In many cases, the genes have subsequently acquired loss-of-function mutations to become pseudogenes. As predicted by theory, this tends to occur first in cifB, whose sole function is to modify sperm, and then in cifA, which is required to rescue the cross in females. Although cif genes recombine, recombination is largely restricted to closely related homologs. This is predicted under a model of coevolution between sperm modification and embryonic rescue, where recombination between distantly related pairs of genes would create a self-incompatible strain. Together, these patterns of gene gain, loss, and recombination support evolutionary models of cytoplasmic incompatibility.

Genome Reduction Is Associated with Bacterial Pathogenicity across Different Scales of Temporal and Ecological Divergence.

Emerging bacterial pathogens threaten global health and food security, and so it is important to ask whether these transitions to pathogenicity have any common features. We present a systematic study of the claim that pathogenicity is associated with genome reduction and gene loss. We compare broad-scale patterns across all bacteria, with detailed analyses of Streptococcus suis, an emerging zoonotic pathogen of pigs, which has undergone multiple transitions between disease and carriage forms. We find that pathogenicity is consistently associated with reduced genome size across three scales of divergence (between species within genera, and between and within genetic clusters of S. suis). Although genome reduction is also found in mutualist and commensal bacterial endosymbionts, genome reduction in pathogens cannot be solely attributed to the features of their ecology that they share with these species, that is, host restriction or intracellularity. Moreover, other typical correlates of genome reduction in endosymbionts (reduced metabolic capacity, reduced GC content, and the transient expansion of nonfunctional elements) are not consistently observed in pathogens. Together, our results indicate that genome reduction is a consistent correlate of pathogenicity in bacteria.

Prevalence and polymorphism of a mussel transmissible cancer in Europe.

Transmissible cancers are parasitic malignant cell lineages that have acquired the ability to infect new hosts from the same species, or sometimes related species. First described in dogs and Tasmanian devils, transmissible cancers were later discovered in some marine bivalves affected by a leukaemia-like disease. In Mytilus mussels, two lineages of bivalve transmissible neoplasia (BTN) have been described to date (MtrBTN1 and MtrBTN2), both of which emerged in a Mytilus trossulus founder individual. Here, we performed extensive screening of genetic chimerism, a hallmark of transmissible cancer, by genotyping 106 single nucleotide polymorphisms of 5,907 European Mytilus mussels. Genetic analysis allowed us to simultaneously obtain the genotype of hosts - Mytilus edulis, M. galloprovincialis or hybrids - and the genotype of tumours of heavily infected individuals. In addition, a subset of 222 individuals were systematically genotyped and analysed by histology to screen for possible nontransmissible cancers. We detected MtrBTN2 at low prevalence in M. edulis, and also in M. galloprovincialis and hybrids although at a much lower prevalence. No MtrBTN1 or new BTN were found, but eight individuals with nontransmissible neoplasia were observed at a single polluted site on the same sampling date. We observed a diversity of MtrBTN2 genotypes that appeared more introgressed or more ancestral than MtrBTN1 and reference healthy M. trossulus individuals. The observed polymorphism is probably due to somatic null alleles caused by structural variations or point mutations in primer-binding sites leading to enhanced detection of the host alleles. Despite low prevalence, two sublineages divergent by 10% fixed somatic null alleles and one nonsynonymous mtCOI (mitochondrial cytochrome oxidase I) substitution are cospreading in the same geographical area, suggesting a complex diversification of MtrBTN2 since its emergence and host species shift.

Large-scale genomic analysis of antimicrobial resistance in the zoonotic pathogen Streptococcus suis.

BACKGROUND: Antimicrobial resistance (AMR) is among the gravest threats to human health and food security worldwide. The use of antimicrobials in livestock production can lead to emergence of AMR, which can have direct effects on humans through spread of zoonotic disease. Pigs pose a particular risk as they are a source of zoonotic diseases and receive more antimicrobials than most other livestock. Here we use a large-scale genomic approach to characterise AMR in Streptococcus suis, a commensal found in most pigs, but which can also cause serious disease in both pigs and humans. RESULTS: We obtained replicated measures of Minimum Inhibitory Concentration (MIC) for 16 antibiotics, across a panel of 678 isolates, from the major pig-producing regions of the world. For several drugs, there was no natural separation into 'resistant' and 'susceptible', highlighting the need to treat MIC as a quantitative trait. We found differences in MICs between countries, consistent with their patterns of antimicrobial usage. AMR levels were high even for drugs not used to treat S. suis, with many multidrug-resistant isolates. Similar levels of resistance were found in pigs and humans from regions associated with zoonotic transmission. We next used whole genome sequences for each isolate to identify 43 candidate resistance determinants, 22 of which were novel in S. suis. The presence of these determinants explained most of the variation in MIC. But there were also interesting complications, including epistatic interactions, where known resistance alleles had no effect in some genetic backgrounds. Beta-lactam resistance involved many core genome variants of small effect, appearing in a characteristic order. CONCLUSIONS: We present a large dataset allowing the analysis of the multiple contributing factors to AMR in S. suis. The high levels of AMR in S. suis that we observe are reflected by antibiotic usage patterns but our results confirm the potential for genomic data to aid in the fight against AMR.

How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels.

The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study 'replicated' instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.

The distribution of epistasis on simple fitness landscapes.

Fitness interactions between mutations can influence a population's evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA (small nucleolar RNA). Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.

The Phylogeny of Rickettsia Using Different Evolutionary Signatures: How Tree-Like is Bacterial Evolution?

Rickettsia is a genus of intracellular bacteria whose hosts and transmission strategies are both impressively diverse, and this is reflected in a highly dynamic genome. Some previous studies have described the evolutionary history of Rickettsia as non-tree-like, due to incongruity between phylogenetic reconstructions using different portions of the genome. Here, we reconstruct the Rickettsia phylogeny using whole-genome data, including two new genomes from previously unsampled host groups. We find that a single topology, which is supported by multiple sources of phylogenetic signal, well describes the evolutionary history of the core genome. We do observe extensive incongruence between individual gene trees, but analyses of simulations over a single topology and interspersed partitions of sites show that this is more plausibly attributed to systematic error than to horizontal gene transfer. Some conflicting placements also result from phylogenetic analyses of accessory genome content (i.e., gene presence/absence), but we argue that these are also due to systematic error, stemming from convergent genome reduction, which cannot be accommodated by existing phylogenetic methods. Our results show that, even within a single genus, tests for gene exchange based on phylogenetic incongruence may be susceptible to false positives.

Rates of vaccine evolution show strong effects of latency: implications for varicella zoster virus epidemiology.

Varicella-zoster virus (VZV) causes chickenpox and shingles, and is found in human populations worldwide. The lack of temporal signal in the diversity of VZV makes substitution rate estimates unreliable, which is a barrier to understanding the context of its global spread. Here, we estimate rates of evolution by studying live attenuated vaccines, which evolved in 22 vaccinated patients for known periods of time, sometimes, but not always undergoing latency. We show that the attenuated virus evolves rapidly (∼ 10(-6) substitutions/site/day), but that rates decrease dramatically when the virus undergoes latency. These data are best explained by a model in which viral populations evolve for around 13 days before becoming latent, but then undergo no replication during latency. This implies that rates of viral evolution will depend strongly on transmission patterns. Nevertheless, we show that implausibly long latency periods are required to date the most recent common ancestor of extant VZV to an "out-of-Africa" migration with humans, as has been previously suggested.

The evolution and genetics of virus host shifts.

Emerging viral diseases are often the product of a host shift, where a pathogen jumps from its original host into a novel species. Phylogenetic studies show that host shifts are a frequent event in the evolution of most pathogens, but why pathogens successfully jump between some host species but not others is only just becoming clear. The susceptibility of potential new hosts can vary enormously, with close relatives of the natural host typically being the most susceptible. Often, pathogens must adapt to successfully infect a novel host, for example by evolving to use different cell surface receptors, to escape the immune response, or to ensure they are transmitted by the new host. In viruses there are often limited molecular solutions to achieve this, and the same sequence changes are often seen each time a virus infects a particular host. These changes may come at a cost to other aspects of the pathogen's fitness, and this may sometimes prevent host shifts from occurring. Here we examine how these evolutionary factors affect patterns of host shifts and disease emergence.

Local interspecies introgression is the main cause of extreme levels of intraspecific differentiation in mussels.

Structured populations, and replicated zones of contact between species, are an ideal opportunity to study regions of the genome with unusual levels of differentiation; and these can illuminate the genomic architecture of species isolation, and the spread of adaptive alleles across species ranges. Here, we investigated the effects of gene flow on divergence and adaptation in the Mytilus complex of species, including replicated parental populations in quite distant geographical locations. We used target enrichment sequencing of 1269 contigs of a few kb each, including some genes of known function, to infer gene genealogies at a small chromosomal scale. We show that geography is an important determinant of the genomewide patterns of introgression in Mytilus and that gene flow between different species, with contiguous ranges, explained up to half of the intraspecific outliers. This suggests that local introgression is both widespread and tends to affect larger chromosomal regions than purely intraspecific processes. We argue that this situation might be common, and this implies that genome scans should always consider the possibility of introgression from sister species, unsampled differentiated backgrounds, or even extinct relatives, for example Neanderthals in humans. The hypothesis that reticulate evolution over long periods of time contributes widely to adaptation, and to the spatial and genomic reorganization of genetic backgrounds, needs to be more widely considered to make better sense of genome scans.

The incidence of bacterial endosymbionts in terrestrial arthropods.

Intracellular endosymbiotic bacteria are found in many terrestrial arthropods and have a profound influence on host biology. A basic question about these symbionts is why they infect the hosts that they do, but estimating symbiont incidence (the proportion of potential host species that are actually infected) is complicated by dynamic or low prevalence infections. We develop a maximum-likelihood approach to estimating incidence, and testing hypotheses about its variation. We apply our method to a database of screens for bacterial symbionts, containing more than 3600 distinct arthropod species and more than 150 000 individual arthropods. After accounting for sampling bias, we estimate that 52% (CIs: 48-57) of arthropod species are infected with Wolbachia, 24% (CIs: 20-42) with Rickettsia and 13% (CIs: 13-55) with Cardinium. We then show that these differences stem from the significantly reduced incidence of Rickettsia and Cardinium in most hexapod orders, which might be explained by evolutionary differences in the arthropod immune response. Finally, we test the prediction that symbiont incidence should be higher in speciose host clades. But while some groups do show a trend for more infection in species-rich families, the correlations are generally weak and inconsistent. These results argue against a major role for parasitic symbionts in driving arthropod diversification.

Possible footprints of APOBEC3F and/or other APOBEC3 deaminases, but not APOBEC3G, on HIV-1 from patients with acute/early and chronic infections.

UNLABELLED: Members of the apolipoprotein B mRNA-editing enzyme-catalytic polypeptide-like-3 (APOBEC3) innate cellular cytidine deaminase family, particularly APOBEC3F and APOBEC3G, can cause extensive and lethal G-to-A mutations in HIV-1 plus-strand DNA (termed hypermutation). It is unclear if APOBEC3-induced mutations in vivo are always lethal or can occur at sublethal levels that increase HIV-1 diversification and viral adaptation to the host. The viral accessory protein Vif counteracts APOBEC3 activity by binding to APOBEC3 and promoting proteasome degradation; however, the efficiency of this interaction varies, since a range of hypermutation frequencies are observed in HIV-1 patient DNA. Therefore, we examined "footprints" of APOBEC3G and APOBEC3F activity in longitudinal HIV-1 RNA pol sequences from approximately 3,000 chronically infected patients by determining whether G-to-A mutations occurred in motifs that were favored or disfavored by these deaminases. G-to-A mutations were more frequent in APOBEC3G-disfavored than in APOBEC3G-favored contexts. In contrast, mutations in APOBEC3F-disfavored contexts were relatively rare, whereas mutations in contexts favoring APOBEC3F (and possibly other deaminases) occurred 16% more often than average G-to-A mutations. These results were supported by analyses of >500 HIV-1 env sequences from acute/early infection. IMPORTANCE: Collectively, our results suggest that APOBEC3G-induced mutagenesis is lethal to HIV-1, whereas mutagenesis caused by APOBEC3F and/or other deaminases may result in sublethal mutations that might facilitate viral diversification. Therefore, Vif-specific cytotoxic T lymphocyte (CTL) responses and drugs that manipulate the interplay between Vif and APOBEC3 may have beneficial or detrimental clinical effects depending on how they affect the binding of Vif to various members of the APOBEC3 family.

Wolbachia in butterflies and moths: geographic structure in infection frequency.

INTRODUCTION: Butterflies and moths (Lepidoptera) constitute one of the most diverse insect orders, and play an important role in ecosystem function. However, little is known in terms of their bacterial communities. Wolbachia, perhaps the most common and widespread intracellular bacterium on Earth, can manipulate the physiology and reproduction of its hosts, and is transmitted vertically from mother to offspring, or sometimes horizontally between species. While its role in some hosts has been studied extensively, its incidence across Lepidoptera is poorly understood. A recent analysis using a beta-binomial model to infer the between-species distribution of prevalence estimated that approximately 40 % of arthropod species are infected with Wolbachia, but particular taxonomic groups and ecological niches seem to display substantially higher or lower incidences. In this study, we took an initial step and applied a similar, maximum likelihood approach to 300 species of Lepidoptera (7604 individuals from 660 populations) belonging to 17 families and 10 superfamilies, and sampled from 36 countries, representing all continents excluding Antarctica. RESULTS: Approximately a quarter to a third of individuals appear to be infected with Wolbachia, and around 80 % of Lepidoptera species are infected at a non-negligible frequency. This incidence estimate is very high compared to arthropods in general. Wolbachia infection in Lepidoptera is shown to vary between families, but there is no evidence for closely related groups to show similar infection levels. True butterflies (Papilionoidea) are overrepresented in our data, however, our estimates show this group can be taken as a representative for the other major lepidopteran superfamilies. We also show substantial variation in infection level according to geography - closer locations tend to show similar infection levels. We further show that variation in geography is due to a latitudinal gradient in Wolbachia infection, with lower frequencies towards higher latitudes. CONCLUSIONS: Our comprehensive survey of Wolbachia infection in Lepidoptera suggests that infection incidence is very high, and provides evidence that climate and geography are strong predictors of infection frequency.

Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster.

Wolbachia are maternally inherited symbiotic bacteria, commonly found in arthropods, which are able to manipulate the reproduction of their host in order to maximise their transmission. The evolutionary history of endosymbionts like Wolbachia can be revealed by integrating information on infection status in natural populations with patterns of sequence variation in Wolbachia and host mitochondrial genomes. Here we use whole-genome resequencing data from 290 lines of Drosophila melanogaster from North America, Europe, and Africa to predict Wolbachia infection status, estimate relative cytoplasmic genome copy number, and reconstruct Wolbachia and mitochondrial genome sequences. Overall, 63% of Drosophila strains were predicted to be infected with Wolbachia by our in silico analysis pipeline, which shows 99% concordance with infection status determined by diagnostic PCR. Complete Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with strict vertical transmission through the maternal cytoplasm and imperfect transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most recent common ancestor of all Wolbachia and mitochondrial genomes in D. melanogaster dates to around 8,000 years ago. We find evidence for a recent global replacement of ancestral Wolbachia and mtDNA lineages, but our data suggest that the derived wMel lineage arose several thousand years ago, not in the 20th century as previously proposed. Our data also provide evidence that this global replacement event is incomplete and is likely to be one of several similar incomplete replacement events that have occurred since the out-of-Africa migration that allowed D. melanogaster to colonize worldwide habitats. This study provides a complete genomic analysis of the evolutionary mode and temporal dynamics of the D. melanogaster-Wolbachia symbiosis, as well as important resources for further analyses of the impact of Wolbachia on host biology.

APOBEC3G-induced hypermutation of human immunodeficiency virus type-1 is typically a discrete "all or nothing" phenomenon.

The rapid evolution of Human Immunodeficiency Virus (HIV-1) allows studies of ongoing host-pathogen interactions. One key selective host factor is APOBEC3G (hA3G) that can cause extensive and inactivating Guanosine-to-Adenosine (G-to-A) mutation on HIV plus-strand DNA (termed hypermutation). HIV can inhibit this innate anti-viral defense through binding of the viral protein Vif to hA3G, but binding efficiency varies and hypermutation frequencies fluctuate in patients. A pivotal question is whether hA3G-induced G-to-A mutation is always lethal to the virus or if it may occur at sub-lethal frequencies that could increase viral diversification. We show in vitro that limiting-levels of hA3G-activity (i.e. when only a single hA3G-unit is likely to act on HIV) produce hypermutation frequencies similar to those in patients and demonstrate in silico that potentially non-lethal G-to-A mutation rates are ∼10-fold lower than the lowest observed hypermutation levels in vitro and in vivo. Our results suggest that even a single incorporated hA3G-unit is likely to cause extensive and inactivating levels of HIV hypermutation and that hypermutation therefore is typically a discrete "all or nothing" phenomenon. Thus, therapeutic measures that inhibit the interaction between Vif and hA3G will likely not increase virus diversification but expand the fraction of hypermutated proviruses within the infected host.

Inferring host range dynamics from comparative data: the protozoan parasites of new world monkeys.

Uncovering the ecological determinants of parasite host range is a central goal of comparative parasitology and infectious disease ecology. But while parasites are often distributed nonrandomly across the host phylogeny, such patterns are difficult to interpret without a genealogy for the parasite samples and without knowing what sorts of ecological dynamics might lead to what sorts of nonrandomness. We investigated inferences from comparative data, using presence/absence records from protozoan parasites of the New World monkeys. We first demonstrate several distinct types of phylogenetic signal in these data, showing, for example, that parasite species are clustered on the host tree and that closely related host species harbor similar numbers of parasite species. We then show that all of these patterns can be generated by a single, simple dynamical model, in which parasite host range changes more rapidly than host speciation/extinction and parasites preferentially colonize uninfected host species that are closely related to their existing hosts. Fitting this model to data, we then estimate its parameters. Finally, we caution that quite different ecological processes can lead to similar signatures but show how phylogenetic variation in host susceptibility can be distinguished from a tendency for parasites to colonize closely related hosts. Our new process-based analyses, which estimate meaningful parameters, should be useful for inferring the determinants of parasite host range and transmission success.

Standing and flowing: the complex origins of adaptive variation.

A population faced with a new selection pressure can only adapt if appropriate genetic variation is available. This genetic variation might come from new mutations or from gene exchange with other populations or species, or it might already segregate in the population as standing genetic variation (which might itself have arisen from either mutation or gene flow). Understanding the relative importance of these sources of adaptive variation is a fundamental issue in evolutionary genetics (Orr & Betancourt ; Barrett & Schluter ; Gladyshev et al. ) and has practical implications for conservation, plant and animal breeding, biological control and infectious disease prevention (e.g. Robertson ; Soulé & Wilcox ; Prentis et al. ; Pennings ). In this issue of Molecular Ecology, Roesti et al. () make an important contribution to this longstanding debate.