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

Density-Dependent Selection during Range Expansion Affects Expansion Load in Life History Traits

AbstractModels of range expansion have independently explored fitness consequences of life history trait evolution and increased rates of genetic drift-or "allele surfing"-during spatial spread, but no previous model has examined the interactions between these two processes. Here, using spatially explicit simulations, we explore an ecologically complex range expansion scenario that combines density-dependent selection with allele surfing to asses the genetic and fitness consequences of density-dependent selection on the evolution of life history traits. We demonstrate that density-dependent selection on the range edge acts differently depending on the life history trait and can either diminish or enhance allele surfing. Specifically, we show that selection at the range edge is always weaker at sites affecting competitive ability (K-selected traits) than at sites affecting birth rate (r-selected traits). We then link differences in the frequency of deleterious mutations to differences in the efficacy of selection and rate of mutation accumulation across distinct life history traits. Finally, we demonstrate that the observed fitness consequences of allele surfing depend on the population density in which expansion load is measured. Our work highlights the complex relationship between ecology and expressed genetic load, which will be important to consider when interpreting both experimental and field studies of range expansion.

Fitness effects of a demography-dispersal trade-off in expandingSaccharomyces cerevisiaemats

Fungi expand in space and time to form complex multicellular communities. The mechanisms by which they do so can vary dramatically and determine the life-history and dispersal traits of expanding populations. These traits influence deterministic and stochastic components of evolution, resulting in complex eco-evolutionary dynamics during colony expansion. We perform experiments on budding yeast strains genetically engineered to display rough-surface and smooth-surface phenotypes in colony-like structures called 'mats'. Previously, it was shown that the rough-surface strain has a competitive advantage over the smooth-surface strain when grown on semi-solid media. We experimentally observe the emergence and expansion of segments with a distinct smooth-surface phenotype during rough-surface mat development. We propose a trade-off between dispersal and local carrying capacity to explain the relative fitness of these two phenotypes. Using a modified stepping-stone model, we demonstrate that this trade-off gives the high-dispersing, rough-surface phenotype a competitive advantage from standing variation, but that it inhibits this phenotype's ability to invade a resident smooth-surface population via mutation. However, the trade-off improves the ability of the smooth-surface phenotype to invade in rough-surface mats, replicating the frequent emergence of smooth-surface segments in experiments. Together, these computational and experimental findings advance our understanding of the complex eco-evolutionary dynamics of fungal mat expansion.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Past human expansions shaped the spatial pattern of Neanderthal ancestry

The worldwide expansion of modern humans (Homo sapiens) started before the extinction of Neanderthals (Homo neanderthalensis). Both species coexisted and interbred, leading to slightly higher introgression in East Asians than in Europeans. This distinct ancestry level has been argued to result from selection, but range expansions of modern humans could provide an alternative explanation. This hypothesis would lead to spatial introgression gradients, increasing with distance from the expansion source. We investigate the presence of Neanderthal introgression gradients after past human expansions by analyzing Eurasian paleogenomes. We show that the out-of-Africa expansion resulted in spatial gradients of Neanderthal ancestry that persisted through time. While keeping the same gradient orientation, the expansion of early Neolithic farmers contributed decisively to reducing the Neanderthal introgression in European populations compared to Asian populations. This is because Neolithic farmers carried less Neanderthal DNA than preceding Paleolithic hunter-gatherers. This study shows that inferences about past human population dynamics can be made from the spatiotemporal variation in archaic introgression.

Allele surfing causes maladaptation in a Pacific salmon of conservation concern

How various factors, including demography, recombination or genome duplication, may impact the efficacy of natural selection and the burden of deleterious mutations, is a central question in evolutionary biology and genetics. In this study, we show that key evolutionary processes, including variations in i) effective population size (Ne) ii) recombination rates and iii) chromosome inheritance, have influenced the genetic load and efficacy of selection in Coho salmon (Oncorhynchus kisutch), a widely distributed salmonid species on the west coast of North America. Using whole genome resequencing data from 14 populations at different migratory distances from their southern glacial refugium, we found evidence supporting gene surfing, wherein reduced Ne at the postglacial recolonization front, leads to a decrease in the efficacy of selection and a surf of deleterious alleles in the northernmost populations. Furthermore, our results indicate that recombination rates play a prime role in shaping the load along the genome. Additionally, we identified variation in polyploidy as a contributing factor to within-genome variation of the load. Overall, our results align remarkably well with expectations under the nearly neutral theory of molecular evolution. We discuss the fundamental and applied implications of these findings for evolutionary and conservation genomics.

Fine-scale population structure of the northern hard clam (Mercenaria mercenaria) revealed by genome-wide SNP markers

Aquaculture is growing rapidly worldwide, and sustainability is dependent on an understanding of current genetic variation and levels of connectivity among populations. Genetic data are essential to mitigate the genetic and ecological impacts of aquaculture on wild populations and guard against unintended human-induced loss of intraspecific diversity in aquacultured lines. Impacts of disregarding genetics can include loss of diversity within and between populations and disruption of local adaptation patterns, which can lead to a decrease in fitness. The northern hard clam, Mercenaria mercenaria (Linnaeus, 1758), is an economically valuable aquaculture species along the North American Atlantic and Gulf coasts. Hard clams have a pelagic larval phase that allows for dispersal, but the level of genetic connectivity among geographic areas is not well understood. To better inform the establishment of site-appropriate aquaculture brood stocks, this study used DArTseq™ genotyping by sequencing to characterize the genetic stock structure of wild clams sampled along the east coast of North America and document genetic diversity within populations. Samples were collected from 15 locations from Prince Edward Island, Canada, to South Carolina, USA. Stringent data filtering resulted in 4960 single nucleotide polymorphisms from 448 individuals. Five genetic breaks separating six genetically distinct populations were identified: Canada, Maine, Massachusetts, Mid-Atlantic, Chesapeake Bay, and the Carolinas (F ST 0.003-0.046; p < 0.0001). This is the first study to assess population genetic structure of this economically important hard clam along a large portion of its native range with high-resolution genomic markers, enabling identification of previously unrecognized population structure. Results of this study not only broaden insight into the factors shaping the current distribution of M. mercenaria but also reveal the genetic population dynamics of a species with a long pelagic larval dispersal period along the North American Atlantic and Gulf coasts.

The Impact of Genetic Surfing on Neutral Genomic Diversity

Range expansions have been common in the history of most species. Serial founder effects and subsequent population growth at expansion fronts typically lead to a loss of genomic diversity along the expansion axis. A frequent consequence is the phenomenon of "gene surfing," where variants located near the expanding front can reach high frequencies or even fix in newly colonized territories. Although gene surfing events have been characterized thoroughly for a specific locus, their effects on linked genomic regions and the overall patterns of genomic diversity have been little investigated. In this study, we simulated the evolution of whole genomes during several types of 1D and 2D range expansions differing by the extent of migration, founder events, and recombination rates. We focused on the characterization of local dips of diversity, or "troughs," taken as a proxy for surfing events. We find that, for a given recombination rate, once we consider the amount of diversity lost since the beginning of the expansion, it is possible to predict the initial evolution of trough density and their average width irrespective of the expansion condition. Furthermore, when recombination rates vary across the genome, we find that troughs are over-represented in regions of low recombination. Therefore, range expansions can leave local and global genomic signatures often interpreted as evidence of past selective events. Given the generality of our results, they could be used as a null model for species having gone through recent expansions, and thus be helpful to correctly interpret many evolutionary biology studies.

Population genetic models for the spatial spread of adaptive variants: A review in light of SARS-CoV-2 evolution

Theoretical population genetics has long studied the arrival and geographic spread of adaptive variants through the analysis of mathematical models of dispersal and natural selection. These models take on a renewed interest in the context of the COVID-19 pandemic, especially given the consequences that novel adaptive variants have had on the course of the pandemic as they have spread through global populations. Here, we review theoretical models for the spatial spread of adaptive variants and identify areas to be improved in future work, toward a better understanding of variants of concern in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) evolution and other contemporary applications. As we describe, characteristics of pandemics such as COVID-19-such as the impact of long-distance travel patterns and the overdispersion of lineages due to superspreading events-suggest new directions for improving upon existing population genetic models.

Hydrodynamic flow and concentration gradients in the gut enhance neutral bacterial diversity

The gut microbiota features important genetic diversity, and the specific spatial features of the gut may shape evolution within this environment. We investigate the fixation probability of neutral bacterial mutants within a minimal model of the gut that includes hydrodynamic flow and resulting gradients of food and bacterial concentrations. We find that this fixation probability is substantially increased, compared with an equivalent well-mixed system, in the regime where the profiles of food and bacterial concentration are strongly spatially dependent. Fixation probability then becomes independent of total population size. We show that our results can be rationalized by introducing an active population, which consists of those bacteria that are actively consuming food and dividing. The active population size yields an effective population size for neutral mutant fixation probability in the gut.

SARS-CoV-2 variant evolution in the United States: High accumulation of viral mutations over time likely through serial Founder Events and mutational bursts

Since the first case of COVID-19 in December 2019 in Wuhan, China, SARS-CoV-2 has spread worldwide and within a year and a half has caused 3.56 million deaths globally. With dramatically increasing infection numbers, and the arrival of new variants with increased infectivity, tracking the evolution of its genome is crucial for effectively controlling the pandemic and informing vaccine platform development. Our study explores evolution of SARS-CoV-2 in a representative cohort of sequences covering the entire genome in the United States, through all of 2020 and early 2021. Strikingly, we detected many accumulating Single Nucleotide Variations (SNVs) encoding amino acid changes in the SARS-CoV-2 genome, with a pattern indicative of RNA editing enzymes as major mutators of SARS-CoV-2 genomes. We report three major variants through October of 2020. These revealed 14 key mutations that were found in various combinations among 14 distinct predominant signatures. These signatures likely represent evolutionary lineages of SARS-CoV-2 in the U.S. and reveal clues to its evolution such as a mutational burst in the summer of 2020 likely leading to a homegrown new variant, and a trend towards higher mutational load among viral isolates, but with occasional mutation loss. The last quartile of 2020 revealed a concerning accumulation of mostly novel low frequency replacement mutations in the Spike protein, and a hypermutable glutamine residue near the putative furin cleavage site. Finally, end of the year data and 2021 revealed the gradual increase to prevalence of known variants of concern, particularly B.1.1.7, that have acquired additional Spike mutations. Overall, our results suggest that predominant viral genomes are dynamically evolving over time, with periods of mutational bursts and unabated mutation accumulation. This high level of existing variation, even at low frequencies and especially in the Spike-encoding region may become problematic when super-spreader events, akin to serial Founder Events in evolution, drive these rare mutations to prominence.

Spatial population genetics in heavily managed species: Separating patterns of historical translocation from contemporary gene flow in white-tailed deer

Approximately 100 years ago, unregulated harvest nearly eliminated white-tailed deer (Odocoileus virginianus) from eastern North America, which subsequently served to catalyze wildlife management as a national priority. An extensive stock-replenishment effort soon followed, with deer broadly translocated among states as a means of re-establishment. However, an unintended consequence was that natural patterns of gene flow became obscured and pretranslocation signatures of population structure were replaced. We applied cutting-edge molecular and biogeographic tools to disentangle genetic signatures of historical management from those reflecting spatially heterogeneous dispersal by evaluating 35,099 single nucleotide polymorphisms (SNPs) derived via reduced-representation genomic sequencing from 1143 deer sampled statewide in Arkansas. We then employed Simpson's diversity index to summarize ancestry assignments and visualize spatial genetic transitions. Using sub-sampled transects across these transitions, we tested clinal patterns across loci against theoretical expectations of their response under scenarios of re-colonization and restricted dispersal. Two salient results emerged: (A) Genetic signatures from historic translocations are demonstrably apparent; and (B) Geographic filters (major rivers; urban centers; highways) now act as inflection points for the distribution of this contemporary ancestry. These results yielded a statewide assessment of contemporary population structure in deer as driven by historic translocations as well as ongoing processes. In addition, the analytical framework employed herein to effectively decipher extant/historic drivers of deer distribution in Arkansas is also applicable for other biodiversity elements with similarly complex demographic histories.

Insights into invasive species from whole-genome resequencing

Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.

Spatial organization in microbial range expansion emerges from trophic dependencies and successful lineages

Evidence suggests that bacterial community spatial organization affects their ecological function, yet details of the mechanisms that promote spatial patterns remain difficult to resolve experimentally. In contrast to bacterial communities in liquid cultures, surface-attached range expansion fosters genetic segregation of the growing population with preferential access to nutrients and reduced mechanical restrictions for cells at the expanding periphery. Here we elucidate how localized conditions in cross-feeding bacterial communities shape community spatial organization. We combine experiments with an individual based mathematical model to resolve how trophic dependencies affect localized growth rates and nucleate successful cell lineages. The model tracks individual cell lineages and attributes these with trophic dependencies that promote counterintuitive reproductive advantages and result in lasting influences on the community structure, and potentially, on its functioning. We examine persistence of lucky lineages in structured habitats where expansion is interrupted by physical obstacles to gain insights into patterns in porous domains.

The Spatial Signature of Introgression After a Biological Invasion With Hybridization

The accumulation of genome-wide molecular data has emphasized the important role of hybridization in the evolution of many organisms, which may carry introgressed genomic segments resulting from past admixture events with other taxa. Despite a number of examples of hybridization occurring during biological invasions, the resulting spatial patterns of genomic introgression remain poorly understood. Preliminary simulation studies have suggested a heterogeneous spatial level of introgression for invasive taxa after range expansion. We investigated in detail the robustness of this pattern and its persistence over time for both invasive and local organisms. Using spatially explicit simulations, we explored the spatial distribution of introgression across the area of colonization of an invasive taxon hybridizing with a local taxon. The general pattern for neutral loci supported by our results is an increasing introgression of local genes into the invasive taxon with the increase in the distance from the source of the invasion and a decreasing introgression of invasive genes into the local taxon. However, we also show there is some variation in this general trend depending on the scenario investigated. Spatial heterogeneity of introgression within a given taxon is thus an expected neutral pattern in structured populations after a biological invasion with a low to moderate amount of hybridization. We further show that this pattern is consistent with published empirical observations. Using additional simulations, we argue that the spatial pattern of Neanderthal introgression in modern humans, which has been documented to be higher in Asia than in Europe, can be explained by a model of hybridization with Neanderthals in Eurasia during the range expansion of modern humans from Africa. Our results support the view that weak hybridization during range expansion may explain spatially heterogeneous introgression patterns without the need to invoke selection.

The spatial Muller's ratchet: Surfing of deleterious mutations during range expansion

During a range expansion, deleterious mutations can "surf" on the colonization front. The resultant decrease in fitness is known as expansion load. An Allee effect is known to reduce the loss of genetic diversity of expanding populations, by changing the nature of the expansion from "pulled" to "pushed". We study the impact of an Allee effect on the formation of an expansion load with a new model, in which individuals have the genetic structure of a Muller's ratchet. A key feature of Muller's ratchet is that the population fatally accumulates deleterious mutations due to the stochastic loss of the fittest individuals, an event called a click of the ratchet. We observe fast clicks of the ratchet at the colonization front owing to small population size, followed by a slow fitness recovery due to migration of fit individuals from the bulk of the population, leading to a transient expansion load. For large population size, we are able to derive quantitative features of the expansion wave, such as the wave speed and the frequency of individuals carrying a given number of mutations. Using simulations, we show that the presence of an Allee effect reduces the rate at which clicks occur at the front, and thus reduces the expansion load.

Philopatry at the frontier: A demographically driven scenario for the evolution of multilevel societies in baboons (Papio)

The baboons (Papio sp.) exhibit marked interspecies variation in social behavior. The thesis presented here argues, first, that male philopatry is a crucial factor, arguably the crucial factor, underlying the other distinctive features (one-male units, multilevel society) shared by hamadryas and Guinea baboons, but not other species of Papio. The second suggestion is that male philopatry as a population norm was not an adaptation to a particular habitat or set of ecological circumstances but evolved in the common ancestor of hamadryas and Guinea baboons as a response to natural selection in the demographic context peculiar to the frontier of a rapidly expanding population. Other derived features of social structure (male-male tolerance, some facultative female dispersal) subsequently evolved to accommodate male philopatry. The mitochondrial genetic population structure of extant baboons preserves a footprint of the initial expansion of 'modern' Papio. Immediately after the expansion, male-philopatric, multilevel populations with a general physical and behavioral resemblance to Guinea baboons occupied the whole northern hemisphere range of the genus. Behavioral and physical autapomorphies of hamadryas baboons evolved in a subpopulation of this ancestral northern base, in response to a less productive habitat of the Horn of Africa. Subsequently, ancestral olive baboons 'reinvented' male dispersal. They and yellow baboons, another male-dispersing species, then replaced most of the male-philopatric northern populations, by male-driven introgression and nuclear swamping.

Shift in ecological strategy helps marginal populations of shepherd's purse (Capsella bursa-pastoris) to overcome a high genetic load

The outcome of species range expansion depends on the interplay of demographic, environmental and genetic factors. Self-fertilizing species usually show a higher invasive ability than outcrossers but selfing and bottlenecks during colonization also lead to an increased genetic load. The relationship between genomic and phenotypic characteristics of expanding populations has, hitherto, rarely been tested experimentally. We analysed how accessions of the shepherd's purse, Capsella bursa-pastoris, from the colonization front or from the core of the natural range performed under increasing density of competitors. First, accessions from the front showed a lower fitness than those from the core. Second, for all accessions, competitor density impacted negatively both vegetative growth and fruit production. However, despite their higher genetic load and lower absolute performances, accessions from the front were less affected by competition than accessions from the core. This seems to be due to phenotypic trade-offs and a shift in phenology that allow accessions from the front to avoid competition.

Dissection of the mutation accumulation process during bacterial range expansions

BACKGROUND

Recent experimental work has shown that the evolutionary dynamics of bacteria expanding across space can differ dramatically from what we expect under well-mixed conditions. During spatial expansion, deleterious mutations can accumulate due to inefficient selection on the expansion front, potentially interfering with and modifying adaptive evolutionary processes.

RESULTS

We used whole genome sequencing to follow the genomic evolution of 10 mutator Escherichia coli lines during 39 days ( ~ 1650 generations) of a spatial expansion, which allowed us to gain a temporal perspective on the interaction of adaptive and non-adaptive evolutionary processes during range expansions. We used elastic net regression to infer the positive or negative effects of mutations on colony growth. The colony size, measured after three day of growth, decreased at the end of the experiment in all 10 lines, and mutations accumulated at a nearly constant rate over the whole experiment. We find evidence that beneficial mutations accumulate primarily at an early stage of the experiment, leading to a non-linear change of colony size over time. Indeed, the rate of colony size expansion remains almost constant at the beginning of the experiment and then decreases after ~ 12 days of evolution. We also find that beneficial mutations are enriched in genes encoding transport proteins, and genes coding for the membrane structure, whereas deleterious mutations show no enrichment for any biological process.

CONCLUSIONS

Our experiment shows that beneficial mutations target specific biological functions mostly involved in inter or extra membrane processes, whereas deleterious mutations are randomly distributed over the whole genome. It thus appears that the interaction between genetic drift and the availability or depletion of beneficial mutations determines the change in fitness of bacterial populations during range expansion.

Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence

The prevalence and intensity of parasites in wild hosts varies across space and is a key determinant of infection risk in humans, domestic animals and threatened wildlife. Because the immune system serves as the primary barrier to infection, replication and transmission following exposure, we here consider the environmental drivers of immunity. Spatial variation in parasite pressure, abiotic and biotic conditions, and anthropogenic factors can all shape immunity across spatial scales. Identifying the most important spatial drivers of immunity could help pre-empt infectious disease risks, especially in the context of how large-scale factors such as urbanization affect defence by changing environmental conditions. We provide a synthesis of how to apply macroecological approaches to the study of ecoimmunology (i.e. macroimmunology). We first review spatial factors that could generate spatial variation in defence, highlighting the need for large-scale studies that can differentiate competing environmental predictors of immunity and detailing contexts where this approach might be favoured over small-scale experimental studies. We next conduct a systematic review of the literature to assess the frequency of spatial studies and to classify them according to taxa, immune measures, spatial replication and extent, and statistical methods. We review 210 ecoimmunology studies sampling multiple host populations. We show that whereas spatial approaches are relatively common, spatial replication is generally low and unlikely to provide sufficient environmental variation or power to differentiate competing spatial hypotheses. We also highlight statistical biases in macroimmunology, in that few studies characterize and account for spatial dependence statistically, potentially affecting inferences for the relationships between environmental conditions and immune defence. We use these findings to describe tools from geostatistics and spatial modelling that can improve inference about the associations between environmental and immunological variation. In particular, we emphasize exploratory tools that can guide spatial sampling and highlight the need for greater use of mixed-effects models that account for spatial variability while also allowing researchers to account for both individual- and habitat-level covariates. We finally discuss future research priorities for macroimmunology, including focusing on latitudinal gradients, range expansions and urbanization as being especially amenable to large-scale spatial approaches. Methodologically, we highlight critical opportunities posed by assessing spatial variation in host tolerance, using metagenomics to quantify spatial variation in parasite pressure, coupling large-scale field studies with small-scale field experiments and longitudinal approaches, and applying statistical tools from macroecology and meta-analysis to identify generalizable spatial patterns. Such work will facilitate scaling ecoimmunology from individual- to habitat-level insights about the drivers of immune defence and help predict where environmental change may most alter infectious disease risk.

Mutational and Selective Processes Involved in Evolution during Bacterial Range Expansions

Bacterial populations have been shown to accumulate deleterious mutations during spatial expansions that overall decrease their fitness and ability to grow. However, it is unclear if and how they can respond to selection in face of this mutation load. We examine here if artificial selection can counteract the negative effects of range expansions. We examined the molecular evolution of 20 mutator lines selected for fast expansions (SEL) and compared them to 20 other mutator lines freely expanding without artificial selection (CONTROL). We find that the colony size of all 20 SEL lines have increased relative to the ancestral lines, unlike CONTROL lines, showing that enough beneficial mutations are produced during spatial expansions to counteract the negative effect of expansion load. Importantly, SEL and CONTROL lines have similar numbers of mutations indicating that they evolved for the same number of generations and that increased fitness is not due to a purging of deleterious mutations. We find that loss of function mutations better explain the increased colony size of SEL lines than nonsynonymous mutations or a combination of the two. Interestingly, most loss of function mutations are found in simple sequence repeats (SSRs) located in genes involved in gene regulation and gene expression. We postulate that such potentially reversible mutations could play a major role in the rapid adaptation of bacteria to changing environmental conditions by shutting down expensive genes and adjusting gene expression.

Geographic origin and timing of colonization of the Pacific Coast of North America by the rocky shore gastropod Littorina sitkana

The demographic history of a species can have a lasting impact on its contemporary population genetic structure. Northeastern Pacific (NEP) populations of the rocky shore gastropod Littorina sitkana have very little mitochondrial DNA (mtDNA) sequence diversity and show no significant population structure despite lacking dispersive planktonic larvae. A contrasting pattern of high mtDNA diversity in the northwestern Pacific (NWP) suggests that L. sitkana may have recently colonized the NEP from the NWP via stepping-stone colonization through the Aleutian-Commander Archipelago (ACA) following the end of the last glacial 20,000 years ago. Here, we use multi-locus sequence data to test that hypothesis using a combination of descriptive statistics and population divergence modeling aimed at resolving the timing and the geographic origin of NEP populations. Our results show that NEP populations share a common ancestor with a population of L. sitkana on the Kamchatka Peninsula ∼46,900 years ago and that NEP populations diverged from each other ∼21,400 years ago. A more recent population divergence between Kamchatka and NEP populations, than between Kamchatka and other populations in the NWP, suggests that the ACA was the most probable dispersal route. Taking into account the confidence intervals for the estimates, we conservatively estimate that L. sitkana arrived in the NEP between 107,400 and 4,100 years ago, a range of dates that is compatible with post-glacial colonization of the NEP. Unlike other congeners that are relatively abundant in the Pleistocene fossil record of the NEP, only one report of L. sitkana exists from the NEP fossil record. Although broadly consistent with the molecular data, the biogeographic significance of these fossils is difficult to evaluate, as the shells cannot be distinguished from the closely-related congener L. subrotundata.

Expansion history and environmental suitability shape effective population size in a plant invasion

The margins of an expanding range are predicted to be challenging environments for adaptation. Marginal populations should often experience low effective population sizes (N_e ) where genetic drift is high due to demographic expansion and/or census population size is low due to unfavourable environmental conditions. Nevertheless, invasive species demonstrate increasing evidence of rapid evolution and potential adaptation to novel environments encountered during colonization, calling into question whether significant reductions in N_e are realized during range expansions in nature. Here we report one of the first empirical tests of the joint effects of expansion dynamics and environment on effective population size variation during invasive range expansion. We estimate contemporary values of N_e using rates of linkage disequilibrium among genome-wide markers within introduced populations of the highly invasive plant Centaurea solstitialis (yellow starthistle) in North America (California, USA), and within native Eurasian populations. As predicted, we find that N_e within the invaded range is positively correlated with both expansion history (time since founding) and habitat quality (abiotic climate). History and climate had independent additive effects with similar effect sizes, indicating an important role for both factors in this invasion. These results support theoretical expectations for the population genetics of range expansion, though whether these processes can ultimately arrest the spread of an invasive species remains an unanswered question.

Environmental heterogeneity can tip the population genetics of range expansions

The population genetics of most range expansions is thought to be shaped by the competition between Darwinian selection and random genetic drift at the range margins. Here, we show that the evolutionary dynamics during range expansions is highly sensitive to additional fluctuations induced by environmental heterogeneities. Tracking mutant clones with a tunable fitness effect in bacterial colonies grown on randomly patterned surfaces we found that environmental heterogeneity can dramatically reduce the efficacy of selection. Time-lapse microscopy and computer simulations suggest that this effect arises generically from a local 'pinning' of the expansion front, whereby stretches of the front are slowed down on a length scale that depends on the structure of the environmental heterogeneity. This pinning focuses the range expansion into a small number of 'lucky' individuals with access to expansion paths, altering the neutral evolutionary dynamics and increasing the importance of chance relative to selection.

Next-generation phylogeography of the cockle Cerastoderma glaucum: Highly heterogeneous genetic differentiation in a lagoon species

AIM

Coastal lagoons form an intriguing example of fragmented marine habitats. Restricted gene flow among isolated populations of lagoon species may promote their genetic divergence and may thus provide a first step toward speciation. In the present study, the population genetic structure of the lagoon cockle Cerastoderma glaucum has been investigated to clarify the complex phylogeographic pattern found in previous studies, to localize major genetic breaks, and to discuss their origin and maintenance.

LOCATION

The Atlantic and Mediterranean coasts, including the Baltic, North Sea, and Black Sea.

METHODS

A total of 204 C. glaucum individuals from 14 populations were genotyped using restriction site-associated DNA sequencing (RADseq). The genetic diversity, divergence, and structure were analyzed using genome-wide single nucleotide polymorphisms (SNPs). Phylogenetic relationships were inferred under a coalescent model using svdquartets.

RESULTS

The RADseq approach allowed inferring phylogeographic relationships with an unprecedented resolution. Three deeply divergent lineages were identified within C. glaucum that are separated by many genetic barriers: one lineage in the Aegean-Black Sea region, one in the Ionian Sea, and the last one widely distributed from the Western Mediterranean to the Baltic Sea. The nested branching pattern displayed on the species tree largely agrees with the likely scenario of C. glaucum postglacial expansion from the Mediterranean to the Baltic Sea.

MAIN CONCLUSION

The genetic differentiations between geographically separated lagoons proved to be strong, highlighting the evolutionary influence of these naturally fragmented habitats. The postglacial expansion created complex patterns of spatial segregation of genetic diversity with allele frequency gradients in many outlier loci, but also discrepancies between the nuclear and mitochondrial genetic markers that probably arose from genetic surfing of mitochondrial variation.

Why can't we make research grant allocation systems more consistent? A personal opinion

Uncertainty is expected to enter into our grant allocation processes at many points, not limited to those directly involving assessment by peers. The selection of grants for funding is thus prodigiously low in statistical power and will remain so. The replacement of current systems with some form of lottery, as has been proposed, seriously risks weakening the quality of applications. Opportunities exist for agencies to encourage and reward greater clarity and innovation in research outcomes.

Gene exchange between two divergent species of the fungal human pathogen, Coccidioides

The fungal genus Coccidioides is composed of two species, Coccidioides immitis and Coccidioides posadasii. These two species are the causal agents of coccidioidomycosis, a pulmonary disease also known as valley fever. The two species are thought to have shared genetic material due to gene exchange in spite of their long divergence. To quantify the magnitude of shared ancestry between them, we analyzed the genomes of a population sample from each species. Next, we inferred what is the expected size of shared haplotypes that might be inherited from the last common ancestor of the two species and find a cutoff to find what haplotypes have conclusively been exchanged between species. Finally, we precisely identified the breakpoints of the haplotypes that have crossed the species boundary and measure the allele frequency of each introgression in this sample. We find that introgressions are not uniformly distributed across the genome. Most, but not all, of the introgressions segregate at low frequency. Our results show that divergent species can share alleles, that species boundaries can be porous, and highlight the need for a systematic exploration of gene exchange in fungal species.

Causes and Consequences of Spatial Within-Host Viral Spread

The spread of viral pathogens both between and within hosts is inherently a spatial process. While the spatial aspects of viral spread at the epidemiological level have been increasingly well characterized, the spatial aspects of viral spread within infected hosts are still understudied. Here, with a focus on influenza A viruses (IAVs), we first review experimental studies that have shed light on the mechanisms and spatial dynamics of viral spread within hosts. These studies provide strong empirical evidence for highly localized IAV spread within hosts. Since mathematical and computational within-host models have been increasingly used to gain a quantitative understanding of observed viral dynamic patterns, we then review the (relatively few) computational modeling studies that have shed light on possible factors that structure the dynamics of spatial within-host IAV spread. These factors include the dispersal distance of virions, the localization of the immune response, and heterogeneity in host cell phenotypes across the respiratory tract. While informative, we find in these studies a striking absence of theoretical expectations of how spatial dynamics may impact the dynamics of viral populations. To mitigate this, we turn to the extensive ecological and evolutionary literature on range expansions to provide informed theoretical expectations. We find that factors such as the type of density dependence, the frequency of long-distance dispersal, specific life history characteristics, and the extent of spatial heterogeneity are critical factors affecting the speed of population spread and the genetic composition of spatially expanding populations. For each factor that we identified in the theoretical literature, we draw parallels to its analog in viral populations. We end by discussing current knowledge gaps related to the spatial component of within-host IAV spread and the potential for within-host spatial considerations to inform the development of disease control strategies.

Evolution at two time frames: Polymorphisms from an ancient singular divergence event fuel contemporary parallel evolution

When environments change, populations may adapt surprisingly fast, repeatedly and even at microgeographic scales. There is increasing evidence that such cases of rapid parallel evolution are fueled by standing genetic variation, but the source of this genetic variation remains poorly understood. In the saltmarsh beetle Pogonus chalceus, short-winged 'tidal' and long-winged 'seasonal' ecotypes have diverged in response to contrasting hydrological regimes and can be repeatedly found along the Atlantic European coast. By analyzing genomic variation across the beetles' distribution, we reveal that alleles selected in the tidal ecotype are spread across the genome and evolved during a singular and, likely, geographically isolated divergence event, within the last 190 Kya. Due to subsequent admixture, the ancient and differentially selected alleles are currently polymorphic in most populations across its range, which could potentially allow for the fast evolution of one ecotype from a small number of random individuals, as low as 5 to 15, from a population of the other ecotype. Our results suggest that cases of fast parallel ecological divergence can be the result of evolution at two different time frames: divergence in the past, followed by repeated selection on the same divergently evolved alleles after admixture. These findings highlight the importance of an ancient and, likely, allopatric divergence event for driving the rate and direction of contemporary fast evolution under gene flow. This mechanism is potentially driven by periods of geographic isolation imposed by large-scale environmental changes such as glacial cycles.

Mutation load dynamics during environmentally-driven range shifts

The fitness of spatially expanding species has been shown to decrease over time and space, but specialist species tracking their changing environment and shifting their range accordingly have been little studied. We use individual-based simulations and analytical modeling to compare the impact of range expansions and range shifts on genetic diversity and fitness loss, as well as the ability to recover fitness after either a shift or expansion. We find that the speed of a shift has a strong impact on fitness evolution. Fastest shifts show the strongest fitness loss per generation, but intermediate shift speeds lead to the strongest fitness loss per geographic distance. Range shifting species lose fitness more slowly through time than expanding species, however, their fitness measured at equal geographic distances from the source of expansion can be considerably lower. These counter-intuitive results arise from the combination of time over which selection acts and mutations enter the system. Range shifts also exhibit reduced fitness recovery after a geographic shift and may result in extinction, whereas range expansions can persist from the core of the species range. The complexity of range expansions and range shifts highlights the potential for severe consequences of environmental change on species survival.

Decadal stability in genetic variation and structure in the intertidal seaweed Fucus serratus (Heterokontophyta: Fucaceae)

BACKGROUND

The spatial distribution of genetic diversity and structure has important implications for conservation as it reveals a species' strong and weak points with regard to stability and evolutionary capacity. Temporal genetic stability is rarely tested in marine species other than commercially important fishes, but is crucial for the utility of temporal snapshots in conservation management. High and stable diversity can help to mitigate the predicted northward range shift of seaweeds under the impact of climate change. Given the key ecological role of fucoid seaweeds along rocky shores, the positive effect of genetic diversity may reach beyond the species level to stabilize the entire intertidal ecosystem along the temperate North Atlantic. In this study, we estimated the effective population size, as well as temporal changes in genetic structure and diversity of the seaweed F. serratus using 22 microsatellite markers. Samples were taken across latitudes and a range of temperature regimes at seven locations with decadal sampling (2000 and 2010).

RESULTS

Across latitudes, genetic structure and diversity remained stable over 5-10 generations. Stable small-scale structure enhanced regional diversity throughout the species' range. In accordance with its biogeographic history, effective population size and diversity peaked in the species' mid-range in Brittany (France), and declined towards its leading and trailing edge to the north and south. At the species' southern edge, multi-locus-heterozygosity displayed a strong decline from 1999 to 2010.

CONCLUSION

Temporally stable genetic structure over small spatial scales is a potential driver for local adaptation and species radiation in the genus Fucus. Survival and adaptation of the low-diversity leading edge of F. serratus may be enhanced by regional gene flow and 'surfing' of favorable mutations or impaired by the accumulation of deleterious mutations. Our results have clear implications for the conservation of F. serratus at its genetically unique southern edge in Northwest Iberia, where increasing temperatures are likely the major cause for the decline not only of F. serratus, but also other intertidal and subtidal macroalgae. We expect that F. serratus will disappear from Northwest Iberia by 2100 if genetic rescue is not induced by the influx of genetic variation from Brittany.

Behavioural divergence during biological invasions: a study of cane toads (Rhinella marina) from contrasting environments in Hawai'i

Invasive species must deal with novel challenges, both from the alien environment and from pressures arising from range expansion per se (e.g. spatial sorting). Those conditions can create geographical variation in behaviour across the invaded range, as has been documented across regions of Australia invaded by cane toads; range-edge toads are more exploratory and willing to take risks than are conspecifics from the range-core. That behavioural divergence might be a response to range expansion and invasion per se, or to the different environments encountered. Climate differs across the cane toads' invasion range from the wet tropics of Queensland to the seasonally dry climates of northwestern Western Australia. The different thermal and hydric regimes may affect behavioural traits via phenotypic plasticity or through natural selection. We cannot tease apart the effects of range expansion versus climate in an expanding population but can do so in a site where the colonizing species was simultaneously released in all suitable areas, thus removing any subsequent phase of range expansion. Cane toads were introduced to Hawai'i in 1932; and thence to Australia in 1935. Toads were released in all major sugarcane-growing areas in Hawai'i within a 12-month period. Hence, Hawai'ian cane toads provide an opportunity to examine geographical divergence in behavioural traits in a climatically diverse region (each island has both wet and dry sides) in the absence of range expansion subsequent to release. We conducted laboratory-based behavioural trials testing exploration, risk-taking and response to novelty using field-caught toads from the wet and dry sides of two Hawai'ian islands (Oahu and Hawai'i). Toads from the dry side of Oahu had a higher propensity to take risks than did toads from the dry side of Hawai'i. Toads from Oahu were also more exploratory than were conspecifics from the island of Hawai'i. However, toads from wet versus dry climates were similar in all behaviours that we scored, suggesting that founder effects, genetic drift, or developmentally plastic responses to ecological factors other than climate may have driven behavioural divergence between islands.

When increasing population density can promote the evolution of metabolic cooperation

Microbial cooperation drives ecological and epidemiological processes and is affected by the ecology and demography of populations. Population density influences the selection for cooperation, with spatial structure and the type of social dilemma, namely public-goods production or self-restraint, shaping the outcome. While existing theories predict that in spatially structured environments increasing population density can select either for or against cooperation, experimental studies with both public-goods production and self-restraint systems have only ever shown that increasing population density favours cheats. We suggest that the disparity between theory and empirical studies results from experimental procedures not capturing environmental conditions predicted by existing theories to influence the outcome. Our study resolves this issue and provides the first experimental evidence that high population density can favour cooperation in spatially structured environments for both self-restraint and public-goods production systems. Moreover, using a multi-trait mathematical model supported by laboratory experiments we extend this result to systems where the self-restraint and public-goods social dilemmas interact. We thus provide a systematic understanding of how the strength of interaction between the two social dilemmas and the degree of spatial structure within an environment affect selection for cooperation. These findings help to close the current gap between theory and experiments.

Comparative Phylogenomic Assessment of Mitochondrial Introgression among Several Species of Chipmunks (Tamias)

Many species are not completely reproductively isolated, resulting in hybridization and genetic introgression. Organellar genomes, such as those derived from mitochondria (mtDNA) and chloroplasts, introgress frequently in natural systems; however, the forces shaping patterns of introgression are not always clear. Here, we investigate extensive mtDNA introgression in western chipmunks, focusing on species in the Tamias quadrivittatus group from the central and southern Rocky Mountains. Specifically, we investigate the role of selection in driving patterns of introgression. We sequenced 51 mtDNA genomes from six species and combine these sequences with other published genomic data to yield annotated mitochondrial reference genomes for nine species of chipmunks. Genomic characterization was performed using a series of molecular evolutionary and phylogenetic analyses to test protein-coding genes for positive selection. We fit a series of maximum likelihood models using a model-averaging approach, assessed deviations from neutral expectations, and performed additional tests to search for codons under the influence of selection. We found no evidence for positive selection among these genomes, suggesting that selection has not been the driving force of introgression in these species. Thus, extensive mtDNA introgression among several species of chipmunks likely reflects genetic drift of introgressed alleles in historically fluctuating populations.

Range Expansion and the Origin of USA300 North American Epidemic Methicillin-Resistant Staphylococcus aureus

The USA300 North American epidemic (USA300-NAE) clone of methicillin-resistant Staphylococcus aureus has caused a wave of severe skin and soft tissue infections in the United States since it emerged in the early 2000s, but its geographic origin is obscure. Here we use the population genomic signatures expected from the serial founder effects of a geographic range expansion to infer the origin of USA300-NAE and identify polymorphisms associated with its spread. Genome sequences from 357 isolates from 22 U.S. states and territories and seven other countries are compared. We observe two significant signatures of range expansion, including decreases in genetic diversity and increases in derived allele frequency with geographic distance from the Pennsylvania region. These signatures account for approximately half of the core nucleotide variation of this clone, occur genome wide, and are robust to heterogeneity in temporal sampling of isolates, human population density, and recombination detection methods. The potential for positive selection of a gyrA fluoroquinolone resistance allele and several intergenic regions, along with a 2.4 times higher recombination rate in a resistant subclade, is noted. These results are the first to show a pattern of genetic variation that is consistent with a range expansion of an epidemic bacterial clone, and they highlight a rarely considered but potentially common mechanism by which genetic drift may profoundly influence bacterial genetic variation.

The process of geographic spread of an origin population by a series of smaller populations can result in distinctive patterns of genetic variation. We detect these patterns for the first time with an epidemic bacterial clone and use them to uncover the clone’s geographic origin and variants associated with its spread. We study the USA300 clone of methicillin-resistant Staphylococcus aureus, which was first noticed in the early 2000s and subsequently became the leading cause of skin and soft tissue infections in the United States. The eastern United States is the most likely origin of epidemic USA300. Relatively few variants, which include an antibiotic resistance mutation, have persisted during this clone’s spread. Our study suggests that an early chapter in the genetic history of this epidemic bacterial clone was greatly influenced by random subsampling of isolates during the clone’s geographic spread.

Is the behavioural divergence between range-core and range-edge populations of cane toads (Rhinella marina) due to evolutionary change or developmental plasticity?

Individuals at the leading edge of expanding biological invasions often show distinctive phenotypic traits, in ways that enhance their ability to disperse rapidly and to function effectively in novel environments. Cane toads (Rhinella marina) at the invasion front in Australia exhibit shifts in morphology, physiology and behaviour (directionality of dispersal, boldness, risk-taking). We took a common-garden approach, raising toads from range-core and range-edge populations in captivity, to see if the behavioural divergences observed in wild-caught toads are also evident in common-garden offspring. Captive-raised toads from the invasion vanguard population were more exploratory and bolder (more prone to 'risky' behaviours) than toads from the range core, which suggests that these are evolved, genetic traits. Our study highlights the importance of behaviour as being potentially adaptive in invasive populations and adds these behavioural traits to the increasing list of phenotypic traits that have evolved rapidly during the toads' 80-year spread through tropical Australia.

Phase and antigenic variation govern competition dynamics through positioning in bacterial colonies

Cellular positioning towards the surface of bacterial colonies and biofilms can enhance dispersal, provide a selective advantage due to increased nutrient and space availability, or shield interior cells from external stresses. Little is known about the molecular mechanisms that govern bacterial positioning. Using the type IV pilus (T4P) of Neisseria gonorrhoeae, we tested the hypothesis that the processes of phase and antigenic variation govern positioning and thus enhance bacterial fitness in expanding gonococcal colonies. By independently tuning growth rate and T4P-mediated interaction forces, we show that the loss of T4P and the subsequent segregation to the front confers a strong selective advantage. Sequencing of the major pilin gene of the spatially segregated sub-populations and an investigation of the spatio-temporal population dynamics was carried out. Our findings indicate that pilin phase and antigenic variation generate a standing variation of pilin sequences within the inoculation zone, while variants associated with a non-piliated phenotype segregate to the front of the growing colony. We conclude that tuning of attractive forces by phase and antigenic variation is a powerful mechanism for governing the dynamics of bacterial colonies.

The colonization and divergence patterns of Brandt's vole (Lasiopodomys brandtii) populations reveal evidence of genetic surfing

Background

The colonial habit of Brandt’s vole (Lasiopodomys brandtii) differs from that of most other species of the genus Microtus. The demographic history of this species and the patterns shaping its current genetic structure remain unknown. Here, we explored patterns of genetic differentiation and infered the demographic history of Brandt’s vole populations through analyses of nuclear microsatellite and D-loop sequences.

Results

Phylogenetic analyses divided the sampled populations into three main clusters, which represent the southeastern, northeastern and western parts of the total range in Mongolia and China. Molecular data revealed an ancestral area located in the southeast of the extant range, in the Xilinguole District, Inner Mongolia, China, from where Brandt’s vole populations began expanding. A gene flow analysis suggested that the most likely colonization route was from the ancestral area and was followed by subsequent northeastward and westward range expansions. We identified decreases in genetic diversity with increasing distance from the founder population within the newly occupied regions (northeastern and western regions), clinal patterns in the allele frequencies, alleles that were rare in the original area that have become common in the newly occupied regions, and higher genetic differentiation in the expanded range compared with the original one.

Conclusion

Our results indicate that L. brandtii most likely originated from the southeastern part of its current geographic range, and subsequently colonized into the northeastern and western parts by expansion. The genetic patterns among the derived populations and with respect to the original population are consistent with that expected under genetic surfing models, which indicated that genetic drift, rather than gene flow, is the predominant factor underlying the genetic structure of expanding Brandt’s vole populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0995-y) contains supplementary material, which is available to authorized users.

Role of dispersal in resistance evolution and spread

Gene flow via immigration affects rate of evolution of resistance to a pest management tactic, while emigration from a resistant population can spread resistance alleles spatially. Whether resistance detected across the landscape reflects ongoing de novo evolution in different hotspots or spread from a single focal population can determine the most effective mitigation strategy. Pest dispersal dynamics determine the spatio-temporal scale at which mitigation tactics must be applied to contain or reverse resistance in an area. Independent evolution of resistance in different populations appears common but not universal. Conversely, spatial spread appears to be almost inevitable. However, rate and scale of spread depends largely on dispersal dynamics and interplay with factors such as fitness costs, spatially variable selection pressure and whether resistance alleles are spreading through an established population or being carried by populations colonizing new territory.