Deleterious Mutations Can Surf to High Densities on the Wave Front of an Expanding Population

Landscape-Level Patterns of Elevated FS1 Asian Allele Frequencies in Populations of Gypsy Moth (Lepidoptera: Erebidae) at a Northern U.S. Boundary

From a regulatory perspective, Asian gypsy moth is a species complex consisting of three species of Lymantria and two subspecies of Lymantria dispar (L.), differing from the European subspecies, L. dispar dispar (L.), by having consistently flight-capable females. As such, the invasion potential in North America is thought to exceed that of European gypsy moth. USDA-APHIS therefore has a monitoring program to detect Asian gypsy moth at high-risk introduction pathways. Molecular markers are used to improve the diagnosis of Asian gypsy moth. One such marker, which targets the FS1 locus, detects an allele, FS1-A, prevalent in Asian populations but occurring at low frequencies (3-6%) throughout the European gypsy moth's range in North America. However, some locales, such as Minnesota, exhibit elevated FS1-A frequencies. We studied the distribution of the FS1-A allele in northern Minnesota, 2013-2014, assessing spatial patterns in the distribution of the FS1-A allele using Moran's I and using spatial regression techniques to examine if the FS1-A allele was associated with putative movement pathways. We also used time series analysis to discern if temporal patterns in FS1-A or possible introduction events occurred. Our results indicated that FS1-A occurred randomly in space and time. We found no evidence that elevated FS1-A frequencies were associated with movement pathways or possible immigration events into this region over the two years. Elevated frequencies of the FS1-A allele within this region could be due to genetic drift and allelic surfing along the expanding population front, or to selection of physiological or behavioral traits.

Rapid trait evolution drives increased speed and variance in experimental range expansions

Range expansions are central to two ecological issues reshaping patterns of global biodiversity: biological invasions and climate change. Traditional theory considers range expansion as the outcome of the demographic processes of birth, death and dispersal, while ignoring the evolutionary implications of such processes. Recent research suggests evolution could also play a critical role in determining expansion speed but controlled experiments are lacking. Here we use flour beetles (Tribolium castaneum) to show experimentally that mean expansion speed and stochastic variation in speed are both increased by rapid evolution of traits at the expansion edge. We find that higher dispersal ability and lower intrinsic growth rates evolve at the expansion edge compared with spatially nonevolving controls. Furthermore, evolution of these traits is variable, leading to enhanced variance in speed among replicate population expansions. Our results demonstrate that evolutionary processes must be considered alongside demographic ones to better understand and predict range expansions.

Spatial structure provides unique opportunities for evolution during range expansions. Here, the authors show experimentally using the red flour beetle, Tribolium castaneum, that dispersal and growth can evolve through spatial processes, increasing expansion speed and its variance.

Rapid evolution of dispersal ability makes biological invasions faster and more variable

Genetic variation in dispersal ability may result in the spatial sorting of alleles during range expansion. Recent theory suggests that spatial sorting can favour the rapid evolution of life history traits at expanding fronts, and therefore modify the ecological dynamics of range expansion. Here we test this prediction by disrupting spatial sorting in replicated invasions of the bean beetle Callosobruchus maculatus across homogeneous experimental landscapes. We show that spatial sorting promotes rapid evolution of dispersal distance, which increases the speed and variability of replicated invasions: after 10 generations of range expansion, invasions subject to spatial sorting spread 8.9% farther and exhibit 41-fold more variable spread dynamics relative to invasions in which spatial sorting is suppressed. Correspondingly, descendants from spatially evolving invasions exhibit greater mean and variance in dispersal distance. Our results reveal an important role for rapid evolution during invasion, even in the absence of environmental filters, and argue for evolutionarily informed forecasts of invasive spread by exotic species or climate change migration by native species.

Theory suggests that spatial sorting by dispersal ability can generate evolutionarily accelerated range expansions. Using the bean beetle Callosobruchus maculatus, this study shows that evolution not only increases the speed of range expansion, as predicted, but also increases variability.

Finding the Genomic Basis of Local Adaptation: Pitfalls, Practical Solutions, and Future Directions

Uncovering the genetic and evolutionary basis of local adaptation is a major focus of evolutionary biology. The recent development of cost-effective methods for obtaining high-quality genome-scale data makes it possible to identify some of the loci responsible for adaptive differences among populations. Two basic approaches for identifying putatively locally adaptive loci have been developed and are broadly used: one that identifies loci with unusually high genetic differentiation among populations (differentiation outlier methods) and one that searches for correlations between local population allele frequencies and local environments (genetic-environment association methods). Here, we review the promises and challenges of these genome scan methods, including correcting for the confounding influence of a species' demographic history, biases caused by missing aspects of the genome, matching scales of environmental data with population structure, and other statistical considerations. In each case, we make suggestions for best practices for maximizing the accuracy and efficiency of genome scans to detect the underlying genetic basis of local adaptation. With attention to their current limitations, genome scan methods can be an important tool in finding the genetic basis of adaptive evolutionary change.

On the role played by the carrying capacity and the ancestral population size during a range expansion

Most species are structured and various population genetics models have been proposed to investigate their history. For mathematical tractability, most of these models make the simplifying assumption of equilibrium. Here we focus on the properties of a nonequilibrium spatial explicit model, range expansions (REs). Despite their abundance, many details of their genetic consequences need yet to be fully investigated. The model we studied is characterized by four main parameters: the effective population size of each deme (N), the migration rate per generation per deme (m), the time of the expansion (T_exp) and the effective size of the deme from which the expansion started (N_anc). By means of extensive coalescent simulations, we focused on two aspects of range expansions for fixed Nm: (1) the separate influence of N and m and (2) the role of N_anc. We compared our results with an equilibrium stepping stone model and found two main features typical of REs: an excess of rare variants for larger N and a complex interaction between N, T_exp and N_anc in shaping the degree of population differentiation (which depends only on Nm in the stepping stone model). Finally, we developed an approximate Bayesian computation approach to jointly estimate N and m and to infer N_anc. When applied to pseudo-observed data sets, we could correctly recover both N and m (but not N_anc), provided a large number of demes were sampled. These findings highlight how it will be possible to estimate the dispersal rate in nonequilibrium metapopulations by using population genetics approaches.

Genetic signatures of a range expansion in natura: when clones play leapfrog

The genetic consequences of range expansions have generally been investigated at wide geographical and temporal scales, long after the colonization event. A unique ecological system enabled us to both monitor the colonization dynamics and decipher the genetic footprints of expansion over a very short time period. Each year an epidemic of the poplar rust (Melampsora larici‐populina) expands clonally and linearly along the Durance River, in the Alps. The colonization dynamics observed in 2004 showed two phases with different genetic outcomes. Upstream, fast colonization maintained high genetic diversity. Downstream, the colonization wave progressively faltered, diversity eroded, and differentiation increased, as expected under recurrent founder events. In line with the high dispersal abilities of rust pathogens, we provide evidence for leapfrog dispersal of clones. Our results thus emphasize the importance of colonization dynamics in shaping spatial genetic structure in the face of high gene flow.

Metabolic adaptations in a range-expanding arthropod

Despite an increasing number of studies documenting life‐history evolution during range expansions or shifts, we lack a mechanistic understanding of the underlying physiological processes. In this explorative study, we used a metabolomics approach to study physiological changes associated with the recent range expansion of the two‐spotted spider mite (Tetranychus urticae). Mite populations were sampled along a latitudinal gradient from range core to edge and reared under benign common garden conditions for two generations. Using gas chromatography–mass spectrometry, we obtained metabolic population profiles, which showed a gradual differentiation along the latitudinal gradient, indicating (epi)genetic changes in the metabolome in association with range expansion. These changes seemed not related with shifts in the mites’ energetic metabolism, but rather with differential use of amino acids. Particularly, more dispersive northern populations showed lowered concentrations of several essential and nonessential amino acids, suggesting a potential downregulation of metabolic pathways associated with protein synthesis.

Inferring the demographic history underlying parallel genomic divergence among pairs of parasitic and nonparasitic lamprey ecotypes

Understanding the evolutionary mechanisms generating parallel genomic divergence patterns among replicate ecotype pairs remains an important challenge in speciation research. We investigated the genomic divergence between the anadromous parasitic river lamprey (Lampetra fluviatilis) and the freshwater-resident nonparasitic brook lamprey (Lampetra planeri) in nine population pairs displaying variable levels of geographic connectivity. We genotyped 338 individuals with RAD sequencing and inferred the demographic divergence history of each population pair using a diffusion approximation method. Divergence patterns in geographically connected population pairs were better explained by introgression after secondary contact, whereas disconnected population pairs have retained a signal of ancient migration. In all ecotype pairs, models accounting for differential introgression among loci outperformed homogeneous migration models. Generating neutral predictions from the inferred divergence scenarios to detect highly differentiated markers identified greater proportions of outliers in disconnected population pairs than in connected pairs. However, increased similarity in the most divergent genomic regions was found among connected ecotype pairs, indicating that gene flow was instrumental in generating parallelism at the molecular level. These results suggest that heterogeneous genomic differentiation and parallelism among replicate ecotype pairs have partly emerged through restricted introgression in genomic islands.

Allele surfing promotes microbial adaptation from standing variation

The coupling of ecology and evolution during range expansions enables mutations to establish at expanding range margins and reach high frequencies. This phenomenon, called allele surfing, is thought to have caused revolutions in the gene pool of many species, most evidently in microbial communities. It has remained unclear, however, under which conditions allele surfing promotes or hinders adaptation. Here, using microbial experiments and simulations, we show that, starting with standing adaptive variation, range expansions generate a larger increase in mean fitness than spatially uniform population expansions. The adaptation gain results from 'soft' selective sweeps emerging from surfing beneficial mutations. The rate of these surfing events is shown to sensitively depend on the strength of genetic drift, which varies among strains and environmental conditions. More generally, allele surfing promotes the rate of adaptation per biomass produced, which could help developing biofilms and other resource-limited populations to cope with environmental challenges.

Expansion load and the evolutionary dynamics of a species range

Expanding populations incur a mutation burden, the so-called expansion load. Using a mixture of individual-based simulations and analytical modeling, we study the expansion load process in models where population growth depends on the population's fitness (i.e., hard selection). We show that expansion load can severely slow down expansions and limit a species' range, even in the absence of environmental variation. We also study the effect of recombination on the dynamics of a species range and on the evolution of mean fitness on the wave front. If recombination is strong, mean fitness on front approaches an equilibrium value at which the effects of fixed mutations cancel each other out. The equilibrium rate at which new demes are colonized is similar to the rate at which beneficial mutations spread through the core. Without recombination, the dynamics is more complex, and beneficial mutations from the core of the range can invade the front of the expansion, which results in irregular and episodic expansion. Although the rate of adaptation is generally higher in recombining organisms, the mean fitness on the front may be larger in the absence of recombination because high-fitness individuals from the core have a higher chance to invade the front. Our findings have important consequences for the evolutionary dynamics of species ranges as well as on the role and the evolution of recombination during range expansions.

The invasive bighead goby Ponticola kessleri displays large-scale genetic similarities and small-scale genetic differentiation in relation to shipping patterns

Colonization events, range expansions and species invasions leave genetic signatures in the genomes of invasive organisms and produce intricate special patterns. Predictions have been made as to how those patterns arise, but only very rarely, genetic processes can be monitored in real time during range expansions. In an attempt to change that, we track a very recently established invasive population of a fish species, the bighead goby Ponticola kessleri, with high temporal and spatial resolution through 2 years to identify patterns over time. We then compare Swiss and German samples of bighead goby along the river Rhine using microsatellites, mitochondrial D-loop sequences and geometric morphometrics to investigate geographic patterns. We detect weak temporal and strong geographic patterns in the data, which are inconsistent with isolation by distance and indicate long range transport. In search of an explanation for our observations, we analyse the vector properties and travel patterns of commercial vessels on the river Rhine. We present evidence that freshwater cargo ships and tankers are plausible vectors for larvae of invasive goby species. We also present indications that cargo ships and tankers act as differential vectors for this species. In summary, we present genetic data at unique temporal resolution from a vertebrate invasion front and substantiate the paramount role of commercial shipping in freshwater fish translocations.

The genetic backburn: using rapid evolution to halt invasions

The impact of an invasive species depends upon the extent of area across which it ultimately spreads. A powerful strategy for limiting impact, then, is to limit spread, and this can most easily be achieved by managing or reinforcing natural barriers to spread. Using a simulation model, we show that rapid evolutionary increases in dispersal can render permeable an otherwise effective barrier. On the other hand, we also show that, once the barrier is reached, and if it holds, resultant evolutionary decreases in dispersal rapidly make the barrier more effective. Finally, we sketch a strategy--the genetic backburn--in which low-dispersal individuals from the range core are translocated to the nearside of the barrier ahead of the oncoming invasion. We find that the genetic backburn--by preventing invasion front genotypes reaching the barrier, and hastening the evolutionary decrease in dispersal--can make barriers substantially more effective. In our simulations, the genetic backburn never reduced barrier strength, however, the improvement to barrier strength was negligible when there was substantial long-distance dispersal, or when there was no genetic variation for dispersal distance. The improvement in barrier strength also depended on the trade-off between dispersal and competitive ability, with a stronger trade-off conferring greater power to the genetic backburn.

The Limits of Natural Selection in a Nonequilibrium World

Evolutionary theory predicts that factors such as a small population size or low recombination rate can limit the action of natural selection. The emerging field of comparative population genomics offers an opportunity to evaluate these hypotheses. However, classical theoretical predictions assume that populations are at demographic equilibrium. This assumption is likely to be violated in the very populations researchers use to evaluate selection's limits: populations that have experienced a recent shift in population size and/or effective recombination rates. Here we highlight theory and data analyses concerning limitations on the action of natural selection in nonequilibrial populations and argue that substantial care is needed to appropriately test whether species and populations show meaningful differences in selection efficacy. A move toward model-based inferences that explicitly incorporate nonequilibrium dynamics provides a promising approach to more accurately contrast selection efficacy across populations and interpret its significance.

Expanding population edges: theories, traits, and trade-offs

Recent patterns of global change have highlighted the importance of understanding the dynamics and mechanisms of species range shifts and expansions. Unique demographic features, spatial processes, and selective pressures can result in the accumulation and evolution of distinctive phenotypic traits at the leading edges of expansions. We review the characteristics of expanding range margins and highlight possible mechanisms for the appearance of phenotypic differences between individuals at the leading edge and core of the range. The development of life history traits that increase dispersal or reproductive ability is predicted by theory and supported with extensive empirical evidence. Many examples of rapid phenotypic change are associated with trade-offs that may influence the persistence of the trait once expansion ends. Accounting for the effects of edge phenotypes and related trade-offs could be critical for predicting the spread of invasive species and population responses to climate change.

How Obstacles Perturb Population Fronts and Alter Their Genetic Structure

As populations spread into new territory, environmental heterogeneities can shape the population front and genetic composition. We focus here on the effects of an important building block of heterogeneous environments, isolated obstacles. With a combination of experiments, theory, and simulation, we show how isolated obstacles both create long-lived distortions of the front shape and amplify the effect of genetic drift. A system of bacteriophage T7 spreading on a spatially heterogeneous Escherichia coli lawn serves as an experimental model system to study population expansions. Using an inkjet printer, we create well-defined replicates of the lawn and quantitatively study the population expansion of phage T7. The transient perturbations of the population front found in the experiments are well described by a model in which the front moves with constant speed. Independent of the precise details of the expansion, we show that obstacles create a kink in the front that persists over large distances and is insensitive to the details of the obstacle's shape. The small deviations between experimental findings and the predictions of the constant speed model can be understood with a more general reaction-diffusion model, which reduces to the constant speed model when the obstacle size is large compared to the front width. Using this framework, we demonstrate that frontier genotypes just grazing the side of an isolated obstacle increase in abundance, a phenomenon we call 'geometry-enhanced genetic drift', complementary to the founder effect associated with spatial bottlenecks. Bacterial range expansions around nutrient-poor barriers and stochastic simulations confirm this prediction. The effect of the obstacle on the genealogy of individuals at the front is characterized by simulations and rationalized using the constant speed model. Lastly, we consider the effect of two obstacles on front shape and genetic composition of the population illuminating the effects expected from complex environments with many obstacles.

Variability of innate immune system genes in Native American populations-relationship with history and epidemiology

OBJECTIVES

The immune system of a host, defending him/her against invading pathogens, has two main subsystems: innate immunity and acquired immunity. There are several evidences showing that Native American populations are immunologically different from non-Native populations. Our aim was to describe the variability of innate immune system genes in Native American populations.

MATERIALS AND METHODS

We investigated heterozygozities and patterns of population differentiation (FST ) of 14 polymorphisms related to the innate immune response in five Native American populations (Aché, Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante) and the results were compared with the three major world population data (YRI, CEU, and CHB) available at the 1,000 genomes database.

RESULTS

Mean heterozygosities ranged between 0.241 ± 0.057 (Aché) and 0.343 ± 0.033 (Kaingang), but no significant differences were observed (Friedman test, P = 0.197). Mean heterozygosities were also not significantly different when Amerindians were pooled and compared with the 1000 genomes populations (Friedman test, P = 0.506). When the Native American populations were grouped as Amerindians, a significantly higher FST value (0.194) was observed between the Amerindian and African populations. The Ewens-Watterson neutrality test showed that these markers are not under strong selective pressure.

DISCUSSION

Native American populations present similar levels of heterozygosity as those of other continents, but are different from Africans in the frequency of polymorphisms of innate immune genes. This higher differentiation is probably due to demographic processes that occurred during the out-of-Africa event.

Contrasting Genomic Diversity in Two Closely Related Postharvest Pathogens: Penicillium digitatum and Penicillium expansum

Penicillium digitatum and Penicillium expansum are two closely related fungal plant pathogens causing green and blue mold in harvested fruit, respectively. The two species differ in their host specificity, being P. digitatum restricted to citrus fruits and P. expansum able to infect a wide range of fruits after harvest. Although host-specific Penicillium species have been found to have a smaller gene content, it is so far unclear whether these different host specificities impact genome variation at the intraspecific level. Here we assessed genome variation across four P. digitatum and seven P. expansum isolates from geographically distant regions. Our results show very high similarity (average 0.06 SNPs [single nucleotide polymorphism] per kb) between globally distributed isolates of P. digitatum pointing to a recent expansion of a single lineage. This low level of genetic variation found in our samples contrasts with the higher genetic variability observed in the similarly distributed P. expansum isolates (2.44 SNPs per kb). Patterns of polymorphism in P. expansum indicate that recombination exists between genetically diverged strains. Consistent with the existence of sexual recombination and heterothallism, which was unknown for this species, we identified the two alternative mating types in different P. expansum isolates. Patterns of polymorphism in P. digitatum indicate a recent clonal population expansion of a single lineage that has reached worldwide distribution. We suggest that the contrasting patterns of genomic variation between the two species reflect underlying differences in population dynamics related with host specificities and related agricultural practices. It should be noted, however, that this results should be confirmed with a larger sampling of strains, as new strains may broaden the diversity so far found in P. digitatum.

Mitochondrial DNA genetic diversity and LCT-13910 and deltaF508 CFTR alleles typing in the medieval sample from Poland

We attempted to confirm the resemblance of a local medieval population and to reconstruct their contribution to the formation of the modern Polish population at the DNA level. The HVR I mtDNA sequence and two nuclear alleles, LCT-13910C/T SNP and deltaF508 CFTR, were chosen as markers since the distribution of selected nuclear alleles varies among ethnic groups. A total of 47 specimens were selected from a medieval cemetery in Cedynia (located in the western Polish lowland). Regarding the HVR I profile, the analyzed population differed from the present-day population (P = 0.045, F(st) = 0.0103), in contrast to lactase persistence (LP) based on the LCT-13910T allele, thus indicating the lack of notable frequency changes of this allele during the last millennium (P = 0.141). The sequence of the HVR I mtDNA fragment allowed to identify six major haplogroups including H, U5, T, K, and HV0 within the medieval population of Cedynia which are common in today's central Europe. An analysis of haplogroup frequency and its comparison with modern European populations shows that the studied medieval population is more closely related to Finno-Ugric populations than to the present Polish population. Identification of less common haplogroups, i.e., Z and U2, both atypical of the modern Polish population and of Asian origin, provides evidence for some kind of connections between the studied and foreign populations. Furthermore, a comparison of the available aDNA sequences from medieval Europe suggests that populations differed from one another and a number of data from other locations are required to find out more about the features of the medieval gene pool profile.

Evolution of pathogen virulence across space during an epidemic

We explore pathogen virulence evolution during the spatial expansion of an infectious disease epidemic in the presence of a novel host movement trade-off, using a simple, spatially explicit mathematical model. This work is motivated by empirical observations of the Mycoplasma gallisepticum invasion into North American house finch (Haemorhous mexicanus) populations; however, our results likely have important applications to other emerging infectious diseases in mobile hosts. We assume that infection reduces host movement and survival and that across pathogen strains the severity of these reductions increases with pathogen infectiousness. Assuming these trade-offs between pathogen virulence (host mortality), pathogen transmission, and host movement, we find that pathogen virulence levels near the epidemic front (that maximize wave speed) are lower than those that have a short-term growth rate advantage or that ultimately prevail (i.e., are evolutionarily stable) near the epicenter and where infection becomes endemic (i.e., that maximize the pathogen basic reproductive ratio). We predict that, under these trade-offs, less virulent pathogen strains will dominate the periphery of an epidemic and that more virulent strains will increase in frequency after invasion where disease is endemic. These results have important implications for observing and interpreting spatiotemporal epidemic data and may help explain transient virulence dynamics of emerging infectious diseases.

Survival of mutations arising during invasions

When a neutral mutation arises in an invading population, it quickly either dies out or ‘surfs’, i.e. it comes to occupy almost all the habitat available at its time of origin. Beneficial mutations can also surf, as can deleterious mutations over finite time spans. We develop descriptive statistical models that quantify the relationship between the probability that a mutation will surf and demographic parameters for a cellular automaton model of surfing. We also provide a simple analytic model that performs well at predicting the probability of surfing for neutral and beneficial mutations in one dimension. The results suggest that factors – possibly including even abiotic factors – that promote invasion success may also increase the probability of surfing and associated adaptive genetic change, conditioned on such success.

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation--or clines--have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by-product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.

The genetic signature of rapid range expansions: How dispersal, growth and invasion speed impact heterozygosity and allele surfing

As researchers collect spatiotemporal population and genetic data in tandem, models that connect demography and dispersal to genetics are increasingly relevant. The dominant spatiotemporal model of invasion genetics is the stepping-stone model which represents a gradual range expansion in which individuals jump to uncolonized locations one step at a time. However, many range expansions occur quickly as individuals disperse far from currently colonized regions. For these types of expansion, stepping-stone models are inappropriate. To more accurately reflect wider dispersal in many organisms, we created kernel-based models of invasion genetics based on integrodifference equations. Classic theory relating to integrodifference equations suggests that the speed of range expansions is a function of population growth and dispersal. In our simulations, populations that expanded at the same speed but with spread rates driven by dispersal retained more heterozygosity along axes of expansion than range expansions with rates of spread that were driven primarily by population growth. To investigate surfing we introduced mutant alleles in wave fronts of simulated range expansions. In our models based on random mating, surfing alleles remained at relatively low frequencies and surfed less often compared to previous results based on stepping-stone simulations with asexual reproduction.

Impact of range expansions on current human genomic diversity

The patterns of population genetic diversity depend to a large extent on past demographic history. Most human populations are known to have gone recently through a series of range expansions within and out of Africa, but these spatial expansions are rarely taken into account when interpreting observed genomic diversity, possibly because they are difficult to model. Here we review available evidence in favour of range expansions out of Africa, and we discuss several of their consequences on neutral and selected diversity, including some recent observations on an excess of rare neutral and selected variants in large samples. We further show that in spatially subdivided populations, the sampling strategy can severely impact the resulting genetic diversity and be confounded by past demography. We conclude that ignoring the spatial structure of human population can lead to some misinterpretations of extant genetic diversity.

Population genomics of a symbiont in the early stages of a pest invasion

Invasive species often depend on microbial symbionts, but few studies have examined the evolutionary dynamics of symbionts during the early stages of an invasion. The insect Megacopta cribraria and its bacterial nutritional symbiont Candidatus Ishikawaella capsulata invaded the southeastern US in 2009. While M. cribraria was initially discovered on wild kudzu plants, it was found as a pest on soybeans within 1 year of infestation. Because prior research suggests Ishikawaella confers the pest status--that is, the ability to thrive on soybeans--in some Megacopta species, we performed a genomic study on Ishikawaella from US. Megacopta cribraria populations to understand the role of the symbiont in driving host plant preferences. We included Ishikawaella samples collected in the first days of the invasion in 2009 and from 23 locations across the insect's 2011 US range. The 0.75 Mb symbiont genome revealed only 47 fixed differences from the pest-conferring Ishikawaella in Japan, with only one amino acid change in a nutrition-provisioning gene. This similarity, along with a lack of fixed substitutions in the US symbiont population, indicates that Ishikawella likely arrived in the US capable of being a soybean pest. Analyses of allele frequency changes between 2009 and 2011 uncover signatures of both positive and negative selection and suggest that symbionts on soybeans and kudzu experience differential selection for genes related to nutrient provisioning. Our data reveal the evolutionary trajectory of an important insect-bacteria symbiosis in the early stages of an invasion, highlighting the role microbial symbionts may play in the spread of invasive species.

The deleterious mutation load is insensitive to recent population history

Human populations have undergone dramatic changes in population size in the past 100,000 years, including recent rapid growth. How these demographic events have affected the burden of deleterious mutations in individuals and the frequencies of disease mutations in populations remains unclear. We use population genetic models to show that recent human demography has likely had little impact on the average burden of deleterious mutations. This prediction is supported by two exome sequence datasets showing that individuals of west African and European ancestry carry very similar burdens of damaging mutations. We further show that for many diseases, rare alleles are unlikely to contribute a large fraction of the heritable variation, and therefore the impact of recent growth is likely to be modest. However, for those diseases that have a direct impact on fitness, strongly deleterious rare mutations likely do play an important role, and recent growth will have increased their impact.

Human pharmacogenomic variation of antihypertensive drugs: from population genetics to personalized medicine

Aim: To investigate the human pharmacogenetic variation related to antihypertensive drugs, providing a survey of functional interpopulation differences in hypertension pharmacogenes. Materials & methods: The study was divided into two stages. In the first stage, we analyzed 1249 variants located in 57 hypertension pharmacogenes. This first-stage analysis confirmed that geographic origin strongly affects hypertension pharmacogenomic variation and that 31 pharmacogenes are geographically differentiated. In the second stage, we focused our attention on the ethnic-differentiated pharmacogenes, investigating 55,521 genetic variants. In silico analyses were performed to predict the effect of genetic variation. Results: Our analyses indicated functional interpopulation differences, suggesting insight into the mechanisms of antihypertensive drug response. Moreover, our data suggested that rare variants mainly determine the functionality of genes related to antihypertensive drugs. Conclusion: Our study provided important knowledge about the genetics of the antihypertensive drug response, suggesting that next-generation sequencing technologies may develop reliable pharmacogenetic tests for antihypertensive drugs.

Original submitted 19 September 2013; Revision submitted 14 November 2013

On the accumulation of deleterious mutations during range expansions

We investigate the effect of spatial range expansions on the evolution of fitness when beneficial and deleterious mutations cosegregate. We perform individual-based simulations of 1D and 2D range expansions and complement them with analytical approximations for the evolution of mean fitness at the edge of the expansion. We find that deleterious mutations accumulate steadily on the wave front during range expansions, thus creating an expansion load. Reduced fitness due to the expansion load is not restricted to the wave front, but occurs over a large proportion of newly colonized habitats. The expansion load can persist and represent a major fraction of the total mutation load for thousands of generations after the expansion. The phenomenon of expansion load may explain growing evidence that populations that have recently expanded, including humans, show an excess of deleterious mutations. To test the predictions of our model, we analyse functional genetic diversity in humans and find patterns that are consistent with our model.

Effects of local adaptation and interspecific competition on species' responses to climate change

Local adaptation and species interactions have been shown to affect geographic ranges; therefore, we need models of climate impact that include both factors. To identify possible dynamics of species when including these factors, we ran simulations of two competing species using an individual-based, coupled map-lattice model using a linear climatic gradient that varies across latitude and is warmed over time. Reproductive success is governed by an individual's adaptation to local climate as well as its location relative to global constraints. In exploratory experiments varying the strength of adaptation and competition, competition reduces genetic diversity and slows range change, although the two species can coexist in the absence of climate change and shift in the absence of competitors. We also found that one species can drive the other to extinction, sometimes long after climate change ends. Weak selection on local adaptation and poor dispersal ability also caused surfing of cooler-adapted phenotypes from the expanding margin backwards, causing loss of warmer-adapted phenotypes. Finally, geographic ranges can become disjointed, losing centrally-adapted genotypes. These initial results suggest that the interplay between local adaptation and interspecific competition can significantly influence species' responses to climate change, in a way that demands future research.

Rapid shifts in dispersal behavior on an expanding range edge

Dispersal biology at an invasion front differs from that of populations within the range core, because novel evolutionary and ecological processes come into play in the nonequilibrium conditions at expanding range edges. In a world where species' range limits are changing rapidly, we need to understand how individuals disperse at an invasion front. We analyzed an extensive dataset from radio-tracking invasive cane toads (Rhinella marina) over the first 8 y since they arrived at a site in tropical Australia. Movement patterns of toads in the invasion vanguard differed from those of individuals in the same area postcolonization. Our model discriminated encamped versus dispersive phases within each toad's movements and demonstrated that pioneer toads spent longer periods in dispersive mode and displayed longer, more directed movements while they were in dispersive mode. These analyses predict that overall displacement per year is more than twice as far for toads at the invasion front compared with those tracked a few years later at the same site. Studies on established populations (or even those a few years postestablishment) thus may massively underestimate dispersal rates at the leading edge of an expanding population. This, in turn, will cause us to underpredict the rates at which invasive organisms move into new territory and at which native taxa can expand into newly available habitat under climate change.

Human genetic variation of CYP450 superfamily: analysis of functional diversity in worldwide populations

AIM

The present study aimed to investigate the human genetic diversity of the CYP450 superfamily in order to identify functional interethnic differences and analyze the role of CYP450 enzymes in human adaptation.

MATERIALS & METHODS

A computational analysis of genetic and functional differences of the 57 CYP450 genes was performed using the Human Genome Diversity Project and HapMap data; comprising approximately 1694 individuals belonging to 62 human populations.

RESULTS

Twenty-six CYP450 SNPs with F-statistics significantly different than the general distribution were identified. Some showed high differentiation among human populations, suggesting that functional interethnic differences may be present. Indeed, some of these are significantly associated with drug response or disease risk. Furthermore, our data highlighted that TBXAS1 and genes in CYP3A cluster may have a role in some processes of human adaptation.

CONCLUSION

Our study provided an analysis of genetic diversity of CYP450 superfamily, identifying functional differences among ethnic groups and their related clinical phenotypes.

Adaptive evolution during an ongoing range expansion: the invasive bank vole (Myodes glareolus) in Ireland

Range expansions are extremely common, but have only recently begun to attract attention in terms of their genetic consequences. As populations expand, demes at the wave front experience strong genetic drift, which is expected to reduce genetic diversity and potentially cause 'allele surfing', where alleles may become fixed over a wide geographical area even if their effects are deleterious. Previous simulation models show that range expansions can generate very strong selective gradients on dispersal, reproduction, competition and immunity. To investigate the effects of range expansion on genetic diversity and adaptation, we studied the population genomics of the bank vole (Myodes glareolus) in Ireland. The bank vole was likely introduced in the late 1920s and is expanding its range at a rate of ~2.5 km/year. Using genotyping-by-sequencing, we genotyped 281 bank voles at 5979 SNP loci. Fourteen sample sites were arranged in three transects running from the introduction site to the wave front of the expansion. We found significant declines in genetic diversity along all three transects. However, there was no evidence that sites at the wave front had accumulated more deleterious mutations. We looked for outlier loci with strong correlations between allele frequency and distance from the introduction site, where the direction of correlation was the same in all three transects. Amongst these outliers, we found significant enrichment for genic SNPs, suggesting the action of selection. Candidates for selection included several genes with immunological functions and several genes that could influence behaviour.

The fate of cooperation during range expansions

Species expand their geographical ranges following an environmental change, long range dispersal, or a new adaptation. Range expansions not only bring an ecological change, but also affect the evolution of the expanding species. Although the dynamics of deleterious, neutral, and beneficial mutations have been extensively studied in expanding populations, the fate of alleles under frequency-dependent selection remains largely unexplored. The dynamics of cooperative alleles are particularly interesting because selection can be both frequency and density dependent, resulting in a coupling between population and evolutionary dynamics. This coupling leads to an increase in the frequency of cooperators at the expansion front, and, under certain conditions, the entire front can be taken over by cooperators. Thus, a mixed population wave can split into an expansion wave of only cooperators followed by an invasion wave of defectors. After the splitting, cooperators increase in abundance by expanding into new territories faster than they are invaded by defectors. Our results not only provide an explanation for the maintenance of cooperation but also elucidate the effect of eco-evolutionary feedback on the maintenance of genetic diversity during range expansions. When cooperators do not split away, we find that defectors can spread much faster with cooperators than they would be able to on their own or by invading cooperators. This enhanced rate of expansion in mixed waves could counterbalance the loss of genetic diversity due to the founder effect for mutations under frequency-dependent selection. Although we focus on cooperator-defector interactions, our analysis could also be relevant for other systems described by reaction-diffusion equations.

Cooperation is beneficial for the species as a whole, but, at the level of an individual, defection pays off. Natural selection is then expected to favor defectors and eliminate cooperation. This prediction is in stark contrast with the abundance of cooperation at all levels of biological systems: from bacterial biofilms to ecosystems and human societies. Several explanations have been proposed to resolve this paradox, including direct reciprocity and group selection. Our work, however, builds upon an observation that natural selection on cooperators might depend both on their relative frequency in the population and on the population density. We find that this feedback between the population and evolutionary dynamics can substantially increase the frequency of cooperators at the front of an expanding population, and can even lead to a splitting of cooperators from defectors. After splitting, only cooperators colonize new territories, while defectors slowly invade them from behind. Since range expansions are very common in nature, our work provides a new explanation of the maintenance of cooperation. More generally, the phenomena we describe could be of interest in other situations when coexisting entities spread in space, be it species in ecology or diffusing and reacting molecules in chemical kinetics.

The coalescent in boundary-limited range expansions

Habitat ranges of most species shift over time, for instance due to climate change, human intervention, or adaptation. These demographic changes often have drastic population genetic effects, such as a stochastic resampling of the gene pool through the "surfing" phenomenon. Most models assume that the speed of range expansions is only limited by the dispersal ability of the colonizing species and its reproductive potential. While such models of "phenotype-limited" expansions apply to species invasions, it is clear that many range expansions are limited rather by the slow motion of habitat boundaries, as driven for instance by global warming. Here, we develop a coalescent model to study the genetic impact of such "boundary-limited" range expansions. Our simulations and analysis show that the resulting loss of genetic diversity is markedly lower than in species invasions if large carrying capacities can be maintained up to the habitat frontier. Counterintuitively, we find that the total loss of diversity does not depend on the speed of the range expansion: Slower expansions have a smaller rate of loss, but also last longer. Boundary-limited range expansions exhibit a characteristic genetic footprint and should therefore be distinguished from range expansions limited only by intrinsic characteristics of the species.

Front acceleration by dynamic selection in Fisher population waves

We introduce a minimal model of population range expansion in which the phenotypes of individuals present no selective advantage and differ only in their diffusion rate. We show that such neutral phenotypic variability (i.e., that does not modify the growth rate) alone can yield phenotype segregation at the front edge, even in absence of genetic noise, and significantly impact the dynamical properties of the expansion wave. We present an exact asymptotic traveling wave solution and show analytically that phenotype segregation accelerates the front propagation. The results are compatible with field observations such as invasions of cane toads in Australia or bush crickets in Britain.

Genomic data reveal a complex making of humans

In the last few years, two paradigms underlying human evolution have crumbled. Modern humans have not totally replaced previous hominins without any admixture, and the expected signatures of adaptations to new environments are surprisingly lacking at the genomic level. Here we review current evidence about archaic admixture and lack of strong selective sweeps in humans. We underline the need to properly model differential admixture in various populations to correctly reconstruct past demography. We also stress the importance of taking into account the spatial dimension of human evolution, which proceeded by a series of range expansions that could have promoted both the introgression of archaic genes and background selection.

Human genomic disease variants: a neutral evolutionary explanation

Many perspectives on the role of evolution in human health include nonempirical assumptions concerning the adaptive evolutionary origins of human diseases. Evolutionary analyses of the increasing wealth of clinical and population genomic data have begun to challenge these presumptions. In order to systematically evaluate such claims, the time has come to build a common framework for an empirical and intellectual unification of evolution and modern medicine. We review the emerging evidence and provide a supporting conceptual framework that establishes the classical neutral theory of molecular evolution (NTME) as the basis for evaluating disease- associated genomic variations in health and medicine. For over a decade, the NTME has already explained the origins and distribution of variants implicated in diseases and has illuminated the power of evolutionary thinking in genomic medicine. We suggest that a majority of disease variants in modern populations will have neutral evolutionary origins (previously neutral), with a relatively smaller fraction exhibiting adaptive evolutionary origins (previously adaptive). This pattern is expected to hold true for common as well as rare disease variants. Ultimately, a neutral evolutionary perspective will provide medicine with an informative and actionable framework that enables objective clinical assessment beyond convenient tendencies to invoke past adaptive events in human history as a root cause of human disease.

Allee effect promotes diversity in traveling waves of colonization

Most mathematical studies on expanding populations have focused on the rate of range expansion of a population. However, the genetic consequences of population expansion remain an understudied body of theory. Describing an expanding population as a traveling wave solution derived from a classical reaction-diffusion model, we analyze the spatio-temporal evolution of its genetic structure. We show that the presence of an Allee effect (i.e., a lower per capita growth rate at low densities) drastically modifies genetic diversity, both in the colonization front and behind it. With an Allee effect (i.e., pushed colonization waves), all of the genetic diversity of a population is conserved in the colonization front. In the absence of an Allee effect (i.e., pulled waves), only the furthest forward members of the initial population persist in the colonization front, indicating a strong erosion of the diversity in this population. These results counteract commonly held notions that the Allee effect generally has adverse consequences. Our study contributes new knowledge to the surfing phenomenon in continuous models without random genetic drift. It also provides insight into the dynamics of traveling wave solutions and leads to a new interpretation of the mathematical notions of pulled and pushed waves.