Accumulation of Mutational Load at the Edges of a Species Range

A Historical Misstep: Niche Shift to Specialisation Is Pushing Insular Ginger Towards an Evolutionary Dead End

Niche specialisation is a double-edged sword as it aids species in adapting to a particular environment but makes them more susceptible to environmental change, which may result in species extinction. Although it has long been debated whether niche specialisation necessarily falls into an 'evolutionary dead end', empirical evidence from a population genetics perspective remains scant, especially when comparing both ecological generalists and specialists simultaneously. In this study, we scrutinised two Taiwan endemic gingers (Zingiber pleiostachyum and Z. shuanglongense) to evaluate how their contrasting patterns in niche breadth evolution have shaped their evolutionary trajectories. We utilised a genome-wide sequencing approach to investigate the demographic histories of each species, assess their maladaptation to future climate change, and estimate their mutational loads. Our results revealed distinct demographic histories between these two gingers. Z. shuanglongense, as the specialist, despite an initial increase during the Last Glacial Maximum (~22 Kya), has been subjected to a long-term decrease in effective population size (Ne), while Z. pleiostachyum is on the contrary increasing, leading to a significantly larger current Ne. Furthermore, ecological specialists are much more vulnerable to future climate change and exhibit greater drift-associated deleterious mutations compared to generalists, directly affecting species' fitness. This study strongly supports the idea that the transition in niche breadth towards specialisation will push Z. shuanglongense perilously close to extinction and also sheds light on species conservation within limited migratory space.

The Demographic History of Populations and Genomic Imprinting have Shaped the Transposon Patterns in Arabidopsis lyrata

Purifying selection is expected to prevent the accumulation of transposable elements (TEs) within their host, especially when located in and around genes and if affected by epigenetic silencing. However, positive selection may favor the spread of TEs, causing genomic imprinting under parental conflict, as genomic imprinting allows parent-specific influence over resource accumulation to the progeny. Concomitantly, the number and frequency of TE insertions in natural populations are conditioned by demographic events. In this study, we aimed to test how demography and selective forces interact to affect the accumulation of TEs around genes, depending on their epigenetic silencing, with a particular focus on imprinted genes. To this aim, we compared the frequency and distribution of TEs in Arabidopsis lyrata from Europe and North America. Generally, we found that TE insertions showed a lower frequency when they were inserted in or near genes, especially TEs targeted by epigenetic silencing, suggesting purifying selection at work. We also found that many TEs were lost or got fixed in North American populations during the colonization and the postglacial range expansion from refugia of the species in North America, as well as during the transition to selfing, suggesting a potential "TE load." Finally, we found that silenced TEs increased in frequency and even tended to reach fixation when they were linked to imprinted genes. We conclude that in A. lyrata, genomic imprinting has spread in natural populations through demographic events and positive selection acting on silenced TEs, potentially under a parental conflict scenario.

SNP-RFLP Markers for the Study of Arabidopsis lyrata

Arabidopsis lyrata has become a useful system for the study of comparative genomics, hybridization, polyploidization, and evolutionary transitions from outcrossing to selfing. Previous studies of its mating system have used microsatellite loci, but low allelic diversity, particularly in self-compatible populations characterized by low levels of outcrossing, reduces the utility of these markers for more detailed studies. Here, we aimed to develop population-level SNP markers for A. lyrata ssp. lyrata sampled from a self-compatible population at Rondeau Provincial Park, Ontario, Canada. We performed de novo SNP discovery and identified 6808 putative SNPs from genome-wide sequences of 22 individuals originating from a highly selfing population. Further filtering and marker validation enabled the development of 17 SNP marker loci that can be visualized using standard PCR-RFLP protocols. These markers had average minor-allele frequencies of 0.40 in the target population, and four of seven markers were variable in a small sample from nine other A. lyrata populations. These PCR-RFLP markers have the potential to be useful for the analysis of mating patterns within and beyond the inbred self-compatible populations of A. lyrata studied here and enable the continued development of A. lyrata as a model for studying evolutionary transitions from outcrossing to selfing.

Genomic predictions of invasiveness and adaptability of the cotton bollworm in response to climate change

Agricultural pests cause enormous losses in annual agricultural production. Understanding the evolutionary responses and adaptive capacity of agricultural pests under climate change is crucial for establishing sustainable and environmentally friendly agricultural pest management. In this study, we integrate climate modeling and landscape genomics to investigate the distributional dynamics of the cotton bollworm (Helicoverpa armigera) in the adaptation to local environments and resilience to future climate change. Notably, the predicted inhabitable areas with higher suitability for the cotton bollworm could be eight times larger in the coming decades. Climate change is one of the factors driving the dynamics of distribution and population differentiation of the cotton bollworm. Approximately 19,000 years ago, the cotton bollworm expanded from its ancestral African population, followed by gradual occupations of the European, Asian, Oceanian, and American continents. Furthermore, we identify seven subpopulations with high dispersal and adaptability which may have an increased risk of invasion potential. Additionally, a large number of candidate genes and SNPs linked to climatic adaptation were mapped. These findings could inform sustainable pest management strategies in the face of climate change, aiding future pest forecasting and management planning.

Genomic implications of the repeated shift to self-fertilization across a species' geographic distribution

The ability to self-fertilize often varies among closely related hermaphroditic plant species, though, variation can also exist within species. In the North American Arabidopsis lyrata, the shift from self-incompatibility (SI) to selfing established in multiple regions independently, mostly since recent postglacial range expansion. This has made the species an ideal model for the investigation of the genomic basis of the breakdown of SI and its population genetic consequences. By comparing nearby selfing and outcrossing populations across the entire species' geographic distribution, we investigated variation at the self-incompatibility (S-)locus and across the genome. Furthermore, a diallel crossing experiment on one mixed-mating population was performed to gain insight into the inheritance of mating system variation. We confirmed that the breakdown of SI had evolved in several S-locus backgrounds. The diallel suggested the involvement of biparental contributions with dominance relations. Though, the population-level genome-wide association study did not single out clear-cut candidate genes but several regions with one near the S-locus. On the implication side, selfing as compared to outcrossing populations had less than half of the genomic diversity, while the number and length of runs of homozygosity (ROHs) scaled with the degree of inbreeding. Selfing populations with a history of long expansion had the longest ROHs. The results highlight that mating system shift to selfing, its genetic underpinning and the likely negative genomic consequences for evolutionary potential can be strongly interlinked with past range dynamics.

Too big to purge: persistence of deleterious Mutations in Island populations of the European Barn Owl (Tyto alba)

A key aspect of assessing the risk of extinction/extirpation for a particular wild species or population is the status of inbreeding, but the origin of inbreeding and the current mutational load are also two crucial factors to consider when determining survival probability of a population. In this study, we used samples from 502 barn owls from continental and island populations across Europe, with the aim of quantifying and comparing the level of inbreeding between populations with differing demographic histories. In addition to comparing inbreeding status, we determined whether inbreeding is due to non-random mating or high co-ancestry within the population. We show that islands have higher levels of inbreeding than continental populations, and that this is mainly due to small effective population sizes rather than recent consanguineous mating. We assess the probability that a region is autozygous along the genome and show that this probability decreased as the number of genes present in that region increased. Finally, we looked for evidence of reduced selection efficiency and purging in island populations. Among island populations, we found an increase in numbers of both neutral and deleterious minor alleles, possibly as a result of drift and decreased selection efficiency but we found no evidence of purging.

Lineage Differentiation and Genomic Vulnerability in a Relict Tree From Subtropical Forests

The subtropical forests of East Asia are renowned for their high plant diversity, particularly the abundance of ancient relict species. However, both the evolutionary history of these relict species and their capacity for resilience in the face of impending climatic changes remain unclear. Using whole-genome resequencing data, we investigated the lineage differentiation and demographic history of the relict and endangered tree, Bretschneidera sinensis (Akaniaceae). We employed a combination of population genomic and landscape genomic approaches to evaluate variation in mutation load and genomic offset, aiming to predict how different populations may respond to climate change. Our analysis revealed a profound genomic divergence between the East and West lineages, likely as the result of recurrent bottlenecks due to climatic fluctuations during the glacial period. Furthermore, we identified several genes potentially linked to growth characteristics and hypoxia response that had been subjected to positive selection during the lineage differentiation. Our assessment of genomic vulnerability uncovered a significantly higher mutation load and genomic offset in the edge populations of B. sinensis compared to their core counterparts. This implies that the edge populations are likely to experience the most significant impact from the predicted climate conditions. Overall, our research sheds light on the historical lineage differentiation and contemporary genomic vulnerability of B. sinensis. Broadening our understanding of the speciation history and future resilience of relict and endangered species such as B. sinensis, is crucial in developing effective conservation strategies in anticipation of future climatic changes.

The genetic architecture of repeated local adaptation to climate in distantly related plants

Closely related species often use the same genes to adapt to similar environments. However, we know little about why such genes possess increased adaptive potential and whether this is conserved across deeper evolutionary lineages. Adaptation to climate presents a natural laboratory to test these ideas, as even distantly related species must contend with similar stresses. Here, we re-analyse genomic data from thousands of individuals from 25 plant species as diverged as lodgepole pine and Arabidopsis (~300 Myr). We test for genetic repeatability based on within-species associations between allele frequencies in genes and variation in 21 climate variables. Our results demonstrate significant statistical evidence for genetic repeatability across deep time that is not expected under randomness, identifying a suite of 108 gene families (orthogroups) and gene functions that repeatedly drive local adaptation to climate. This set includes many orthogroups with well-known functions in abiotic stress response. Using gene co-expression networks to quantify pleiotropy, we find that orthogroups with stronger evidence for repeatability exhibit greater network centrality and broader expression across tissues (higher pleiotropy), contrary to the 'cost of complexity' theory. These gene families may be important in helping wild and crop species cope with future climate change, representing important candidates for future study.

Forest genomics in the Caucasus through the lens of its dominant tree species - Fagus orientalis

The last glacial period is known to have greatly influenced the demographic history of temperate forest trees, with important range contractions and post-glacial expansions that led to the formation of multiple genetic lineages and secondary contact zones in the Northern Hemisphere. These dynamics have been extensively studied for European and North American species but are still poorly understood in other temperate regions of rich biodiversity such as the Caucasus. Our study helps filling that gap by deciphering the genomic landscapes of F. orientalis across the South Caucasus. The use of genome-wide data confirmed a past demographic history strongly influenced by the Last Glacial Maximum, revealing two disjunct glacial refugia in the Colchis and Hyrcanian regions. The resulting patterns of genetic diversity, load and differentiation are not always concordant across the region, with genetic load pinpointing the location of the glacial refugia more efficiently than genetic diversity alone. The Hyrcanian forests show depleted genetic diversity and substantial isolation, even if long-distance gene flow is still present with the main centre of diversity in the Greater Caucasus. Finally, we characterize a strong heterogeneity of genetic diversity and differentiation along the species chromosomes, with noticeably a first chromosome showing low diversity and weak differentiation.

Limits to species' range: the tension between local and global adaptation

We know that heritable variation is abundant, and that selection causes all but the smallest populations to rapidly shift beyond their original trait distribution. So then, what limits the range of a species? There are physical constraints and also population genetic limits to the effectiveness of selection, ultimately set by population size. Global adaptation, where the same genotype is favoured over the whole range, is most efficient when based on a multitude of weakly selected alleles and is effective even when local demes are small, provided that there is some gene flow. In contrast, local adaptation is sensitive to gene flow and may require alleles with substantial effect. How can populations combine the advantages of large effective size with the ability to specialise into local niches? To what extent does reproductive isolation help resolve this tension? I address these questions using eco-evolutionary models of polygenic adaptation, contrasting discrete demes with continuousspace.

Bringing genomics to the field: An integrative approach to seed sourcing for forest restoration

Premise

Global anthropogenic change threatens the health and productivity of forest ecosystems. Assisted migration and reforestation are tools to help mitigate these impacts. However, questions remain about how to approach sourcing seeds to ensure high establishment and future adaptability.

Methods

Using exome-capture sequencing, we demonstrate a computational approach to finding the best n-sets from a candidate list of seed sources that collectively achieve high genetic diversity (GD) and minimal genetic load (GL), while also increasing evolvability in quantitative traits. The benefits of this three-part strategy (diversity-load-evolvability) are to increase near-term establishment success while also boosting evolutionary potential to respond to future stressors. Members of The Nature Conservancy and the Central Appalachian Spruce Restoration Initiative planted 58,000 seedlings across 255 acres. A subset of seedlings was monitored for establishment success and variation in growth.

Results

The results show gains in GD relative to GL and increases in quantitative genetic variation in seedling growth for pooled vs. single-source restoration. No single "super source" was observed across planting sites; rather, monitoring results demonstrate that pooling of multiple sources helps achieve higher GD:GL and evolvability.

Discussion

Our study shows the potential for integrating genomics into local-scale restoration and the importance of building partnerships between academic researchers and applied conservation managers.

Marine fishes experiencing high-velocity range shifts may not be climate change winners

Climate change is driving the global redistribution of species. A common assumption is that rapid range shifts occur in tandem with overall stable or positive abundance trends throughout the range and thus these species may be considered as climate change 'winners'. However, although establishing the link between range shift velocities and population trends is crucial for predicting climate change impacts it has not been empirically tested. Using 2,572 estimates of changes in marine fish abundance spread across the world's oceans, we show that poleward range shifts are not necessarily associated with positive population trends. Species experiencing high-velocity range shifts seem to experience local population declines irrespective of the position throughout the species range. High range shift velocities of 17 km yr-1 are associated with a 50% decrease in population sizes over a period of 10 yr, which is dramatic compared to the overall stable population trends in non-shifting species. This pattern, however, mostly occurs in populations located in the poleward, colder, portion of the species range. The lack of a positive association between poleward range shift velocities and population trends at the coldest portion of the range contrasts with the view that rapid range shifts safeguard against local population declines. Instead, our work suggests that marine fishes experiencing rapid range shifts could be more vulnerable to climatic change and therefore should be carefully assessed for conservation status.

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.

Limited effects of population age on the genetic structure of spatially isolated forest herb populations in temperate Europe

Due to multiple land-cover changes, forest herb populations residing in forest patches embedded in agricultural landscapes display different ages and, thus, experience differences in genetic exchange, mutation accumulation and genetic drift. The extent of divergence in present-day population genetic structure among these populations of different ages remains unclear, considering their diverse breeding systems and associated pollinators. Answering this question is essential to understand these species' persistence, maintenance of evolutionary potential and adaptability to changing environments. We applied a multi-landscape setup to compare the genetic structure of forest herb populations across forest patches of different ages (18-338 years). We studied the impact on three common slow-colonizer herb species with distinct breeding systems and associated pollinators: Polygonatum multiflorum (outcrossing, long-distance pollinators), Anemone nemorosa (outcrossing, short-distance pollinators) and Oxalis acetosella (mixed breeding). We aimed to assess if in general older populations displayed higher genetic diversity and lower differentiation than younger ones. We also anticipated that P. multiflorum would show the smallest while O. acetosella the largest difference, between old and young populations. We found that older populations had a higher observed heterozygosity (H o) but a similar level of allelic richness (A r) and expected heterozygosity (H e) as younger populations, except for A. nemorosa, which exhibited higher A r and H e in younger populations. As populations aged, their pairwise genetic differentiation measured by D PS decreased independent of species identity while the other two genetic differentiation measures showed either comparable levels between old and young populations (G" ST) or inconsistency among three species (cGD). The age difference of the two populations did not explain their genetic differentiation. Synthesis: We found restricted evidence that forest herb populations with different ages differ in their genetic structure, indicating that populations of different ages can reach a similar genetic structure within decades and thus persist in the long term after habitat disturbance. Despite their distinct breeding systems and associated pollinators, the three studied species exhibited partly similar genetic patterns, suggesting that their common characteristics, such as being slow colonizers or their ability to propagate vegetatively, are important in determining their long-term response to land-cover change.

Genome size variation and evolution during invasive range expansion in an introduced plant

Plants demonstrate exceptional variation in genome size across species, and their genome sizes can also vary dramatically across individuals and populations within species. This aspect of genetic variation can have consequences for traits and fitness, but few studies attributed genome size differentiation to ecological and evolutionary processes. Biological invasions present particularly useful natural laboratories to infer selective agents that might drive genome size shifts across environments and population histories. Here, we test hypotheses for the evolutionary causes of genome size variation across 14 invading populations of yellow starthistle, Centaurea solstitialis, in California, United States. We use a survey of genome sizes and trait variation to ask: (1) Is variation in genome size associated with developmental trait variation? (2) Are genome sizes smaller toward the leading edge of the expansion, consistent with selection for "colonizer" traits? Or alternatively, does genome size increase toward the leading edge of the expansion, consistent with predicted consequences of founder effects and drift? (3) Finally, are genome sizes smaller at higher elevations, consistent with selection for shorter development times? We found that 2C DNA content varied 1.21-fold among all samples, and was associated with flowering time variation, such that plants with larger genomes reproduced later, with lower lifetime capitula production. Genome sizes increased toward the leading edge of the invasion, but tended to decrease at higher elevations, consistent with genetic drift during range expansion but potentially strong selection for smaller genomes and faster development time at higher elevations. These results demonstrate how genome size variation can contribute to traits directly tied to reproductive success, and how selection and drift can shape that variation. We highlight the influence of genome size on dynamics underlying a rapid range expansion in a highly problematic invasive plant.

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

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.

Between but Not Within-Species Variation in the Distribution of Fitness Effects

New mutations provide the raw material for evolution and adaptation. The distribution of fitness effects (DFE) describes the spectrum of effects of new mutations that can occur along a genome, and is, therefore, of vital interest in evolutionary biology. Recent work has uncovered striking similarities in the DFE between closely related species, prompting us to ask whether there is variation in the DFE among populations of the same species, or among species with different degrees of divergence, that is whether there is variation in the DFE at different levels of evolution. Using exome capture data from six tree species sampled across Europe we characterized the DFE for multiple species, and for each species, multiple populations, and investigated the factors potentially influencing the DFE, such as demography, population divergence, and genetic background. We find statistical support for the presence of variation in the DFE at the species level, even among relatively closely related species. However, we find very little difference at the population level, suggesting that differences in the DFE are primarily driven by deep features of species biology, and those evolutionarily recent events, such as demographic changes and local adaptation, have little impact.

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.

Purging due to self-fertilization does not prevent accumulation of expansion load

As species expand their geographic ranges, colonizing populations face novel ecological conditions, such as new environments and limited mates, and suffer from evolutionary consequences of demographic change through bottlenecks and mutation load accumulation. Self-fertilization is often observed at species range edges and, in addition to countering the lack of mates, is hypothesized as an evolutionary advantage against load accumulation through increased homozygosity and purging. We study how selfing impacts the accumulation of genetic load during range expansion via purging and/or speed of colonization. Using simulations, we disentangle inbreeding effects due to demography versus due to selfing and find that selfers expand faster, but still accumulate load, regardless of mating system. The severity of variants contributing to this load, however, differs across mating system: higher selfing rates purge large-effect recessive variants leaving a burden of smaller-effect alleles. We compare these predictions to the mixed-mating plant Arabis alpina, using whole-genome sequences from refugial outcrossing populations versus expanded selfing populations. Empirical results indicate accumulation of expansion load along with evidence of purging in selfing populations, concordant with our simulations, suggesting that while purging is a benefit of selfing evolving during range expansions, it is not sufficient to prevent load accumulation due to range expansion.

Drivers of strong isolation and small effective population size at a leading range edge of a widespread plant

Climate change has influenced species distributions worldwide with upward elevational shifts observed in many systems. Leading range edge populations, like those at upper elevation limits, are crucial for climate change responses but can exhibit low genetic diversity due to founder effects, isolation, or limited outbreeding. These factors can hamper local adaptation at range limits. Using the widespread herb, Argentina anserina, we measured ecological attributes (population density on the landscape, area of population occupancy, and plant and flower density) spanning a 1000 m elevation gradient, with high elevation populations at the range limit. We measured vegetative clonal potential in the greenhouse for populations spanning the gradient. We combined these data with a ddRAD-seq dataset to test the hypotheses that high elevation populations would exhibit ecological and genomic signatures of leading range edge populations. We found that population density on the landscape declined towards the high elevation limit, as is expected towards range edges. However, plant density was elevated within edge populations. In the greenhouse, high elevation plants exhibited stronger clonal potential than low elevation plants, likely explaining increased plant density in the field. Phylogeographic analysis supported more recent colonization of high elevation populations which were also more genetically isolated, had more extreme heterozygote excess and had smaller effective population size than low. Results support that colonization of high elevations was likely accompanied by increased asexuality, contributing to a decline in effective population size. Despite high plant density in leading edge populations, their small effective size, isolation and clonality could constrain adaptive potential.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

Intraspecific variation in reproductive barriers between two closely related Arabidopsis sister species

Reproductive isolation (RI) is a critical component of speciation and varies strongly in timing and strength among different sister taxa, depending on, for example the geography of speciation and divergence time. However, these factors may also produce variation in timing and strength among populations within species. Here we tested for variation in the expression of RI among replicate population pairs between the sister taxa Arabidopsis lyrata subsp. lyrata and A. arenicola. While the former is predominantly outcrossing, the latter is predominantly selfing. We focused on intrinsic prezygotic and postzygotic RI as both species occur largely in allopatry. We assessed RI by performing within-population crosses and interspecific between-population crosses, and by raising offspring. RI was generally high between all interspecific population pairs, but it varied in timing and strength depending on population history. Prezygotic isolation was strongest between the closest-related population pair, while early postzygotic isolation was high for all other population pairs. Furthermore, the timing and strength of RI depended strongly on cross direction. Our study provides empirical support that reproductive barriers between species are highly variable among population pairs and asymmetric within population pairs, and this variation seems to follow patterns typically described across species pairs.

Recent speciation associated with range expansion and a shift to self-fertilization in North American Arabidopsis

The main processes classically evoked for promoting reproductive isolation and speciation are geographic separation reducing gene flow among populations, divergent selection, and chance genomic change. In a case study, we present evidence that the additional factors of climate change, range expansion and a shift in mating towards inbreeding can initiate the processes leading to parapatric speciation. At the end of the last Pleistocene glaciation cycle, the North American plant Arabidopsis lyrata expanded its range and concomitantly lost its reproductive mode of outcrossing multiple times. We show that in one of the newly colonized areas, the self-fertilizing recolonization lineage of A. lyrata gave rise to selfing A. arenicola, which expanded its range to subarctic and arctic Canada and Greenland, while the parental species remained restricted to temperate North America. Despite the vast range expansion by the new species, mutational load did not increase, probably because of selfing and quasi-clonal selection. We conclude that such peripheral parapatric speciation combined with range expansion and inbreeding may be an important but so far overlooked mode of speciation.

Evaluating the persistence and utility of five wild Vitis species in the context of climate change

Crop wild relatives (CWRs) have the capacity to contribute novel traits to agriculture. Given climate change, these contributions may be especially vital for the persistence of perennial crops, because perennials are often clonally propagated and consequently do not evolve rapidly. By studying the landscape genomics of samples from five Vitis CWRs (V. arizonica, V. mustangensis, V. riparia, V. berlandieri and V. girdiana) in the context of projected climate change, we addressed two goals. The first was to assess the relative potential of different CWR accessions to persist in the face of climate change. By integrating species distribution models with adaptive genetic variation, additional genetic features such as genomic load and a phenotype (resistance to Pierce's Disease), we predicted that accessions from one species (V. mustangensis) are particularly well-suited to persist in future climates. The second goal was to identify which CWR accessions may contribute to bioclimatic adaptation for grapevine (V. vinifera) cultivation. To do so, we evaluated whether CWR accessions have the allelic capacity to persist if moved to locations where grapevines are cultivated in the United States. We identified six candidates from V. mustangensis and hypothesized that they may prove useful for contributing alleles that can mitigate climate impacts on viticulture. By identifying candidate germplasm, this study takes a conceptual step toward assessing the genomic and bioclimatic characteristics of CWRs.

Reduced pollinator service in small populations of Arabidopsis lyrata at its southern range limit

Even though a high fraction of angiosperm plants depends on animal pollinators for sexual reproduction, little is known how pollinator service changes across the ranges of plant species and whether it may contribute to range limits. Here, we tested for variation in pollinator service in the North American Arabidopsis lyrata from its southern to northern range edge and evaluated the driving mechanisms. We monitored insect pollinators using time-lapse cameras in 13 populations over two years and spotted 67 pollinating insect taxa, indicating the generalist nature of this plant-pollinator system. Pollinator service was highest at intermediate local flower densities and higher in large compared to small plant populations. Southern populations had generally smaller population sizes, and visitation rate and pollination ratio decreased with latitude. We also found that pollinator visitation was positively correlated with the richness of other flowering plants. This study indicates that plant populations at southern range edges receive only marginal pollinator service if they are small, and the effect of lower pollination is also detectable within populations across the range when the local flower density is low. Results, therefore, suggest the potential for an Allee effect in pollination that manifests itself across spatial scales.

The evolution of recombination in self-fertilizing organisms

Cytological data from flowering plants suggest that the evolution of recombination rates is affected by the mating system of organisms, as higher chiasma frequencies are often observed in self-fertilizing species compared with their outcrossing relatives. Understanding the evolutionary cause of this effect is of particular interest, as it may shed light on the selective forces favoring recombination in natural populations. While previous models showed that inbreeding may have important effects on selection for recombination, existing analytical treatments are restricted to the case of loosely linked loci and weak selfing rates, and ignore the stochastic effect of genetic interference (Hill-Robertson effect), known to be an important component of selection for recombination in randomly mating populations. In this article, we derive general expressions quantifying the stochastic and deterministic components of selection acting on a mutation affecting the genetic map length of a whole chromosome along which deleterious mutations occur, valid for arbitrary selfing rates. The results show that selfing generally increases selection for recombination caused by interference among mutations as long as selection against deleterious alleles is sufficiently weak. While interference is often the main driver of selection for recombination under tight linkage or high selfing rates, deterministic effects can play a stronger role under intermediate selfing rates and high recombination, selecting against recombination in the absence of epistasis, but favoring recombination when epistasis is negative. Individual-based simulation results indicate that our analytical model often provides accurate predictions for the strength of selection on recombination under partial selfing.

The relative role of plasticity and demographic history in Capsella bursa-pastoris: a common garden experiment in Asia and Europe

The colonization success of a species depends on the interplay between its phenotypic plasticity, adaptive potential and demographic history. Assessing their relative contributions during the different phases of a species range expansion is challenging, and requires large-scale experiments. Here, we investigated the relative contributions of plasticity, performance and demographic history to the worldwide expansion of the shepherd's purse, Capsella bursa-pastoris. We installed two large common gardens of the shepherd's purse, a young, self-fertilizing, allopolyploid weed with a worldwide distribution. One common garden was located in Europe, the other in Asia. We used accessions from three distinct genetic clusters (Middle East, Europe and Asia) that reflect the demographic history of the species. Several life-history traits were measured. To explain the phenotypic variation between and within genetic clusters, we analysed the effects of (i) the genetic clusters, (ii) the phenotypic plasticity and its association to fitness and (iii) the distance in terms of bioclimatic variables between the sampling site of an accession and the common garden, i.e. the environmental distance. Our experiment showed that (i) the performance of C. bursa-pastoris is closely related to its high phenotypic plasticity; (ii) within a common garden, genetic cluster was a main determinant of phenotypic differences; and (iii) at the scale of the experiment, the effect of environmental distance to the common garden could not be distinguished from that of genetic clusters. Phenotypic plasticity and demographic history both play important role at different stages of range expansion. The success of the worldwide expansion of C. bursa-pastoris was undoubtedly influenced by its strong phenotypic plasticity.

Environment dependence of the expression of mutational load and species' range limits

Theoretical and empirical research on the causes of species' range limits suggest the contribution of several intrinsic and extrinsic factors, with potentially complex interactions among them. An intrinsic factor proposed by recent theory is mutational load increasing towards range edges because of genetic drift. Furthermore, environmental quality may decline towards range edges and enhance the expression of load. Here, we tested whether the expression of mutational load associated with range limits in the North American plant Arabidopsis lyrata was enhanced under stressful environmental conditions by comparing the performance of within- versus between-population crosses at common garden sites across the species' distribution and beyond. Heterosis, reflecting the expression of load, increased with heightened estimates of genomic load and with environmental stress caused by warming, but the interaction was not significant. We conclude that range-edge populations suffer from a twofold genetic Allee effect caused by increased mutational load and stress-dependent load linked to general heterozygote deficiency, but there is no synergistic effect between them.

A review on trade-offs at the warm and cold ends of geographical distributions

Species' range limits are ubiquitous. This suggests that the evolution of the ecological niche is constrained in general and at the edges of distributions in particular. While there may be many ecological and genetic reasons for this phenomenon, here we focus on the potential role of trade-offs. We performed a literature search on evidence for trade-offs associated with geographical or elevational range limits. The majority of trade-offs were reported as relevant at either the cold end of species' distribution (n = 19), the warm or dry end (n = 19) or both together (n = 14). One common type of trade-off involved accelerating growth or development (27%), often at the cost of small size. Another common type involved resistance to or tolerance of climatic extremes that occur at certain periods of the year (64%), often at the cost of small size or reduced growth. Trade-offs overlapped with some of the classic trade-offs reported in life-history evolution or thermal adaptation. The results highlight several general insights about species' niches and ranges, and we outline how future research should better integrate the ecological context and test for the presence of microevolutionary trade-offs. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Gene surfing of underdominant alleles promotes formation of hybrid zones

The distribution of genetic diversity over geographical space has long been investigated in population genetics and serves as a useful tool to understand evolution and history of populations. Within some species or across regions of contact between two species, there are instances where there is no apparent ecological determinant of sharp changes in allele frequencies or divergence. To further understand these patterns of spatial genetic structure and potential species divergence, we model the establishment of clines that occur due to the surfing of underdominant alleles during range expansions. We provide analytical approximations for the fixation probability of underdominant alleles at expansion fronts and demonstrate that gene surfing can lead to clines in one-dimensional range expansions. We extend these results to multiple loci via a mixture of analytical theory and individual-based simulations. We study the interaction between the strength of selection against heterozygotes, migration rates, and local recombination rates on the formation of stable hybrid zones. Clines created by surfing at different loci can attract each other and align after expansion, if they are sufficiently close in space and in terms of recombination distance. Our findings suggest that range expansions can set the stage for parapatric speciation due to the alignment of multiple selective clines, even in the absence of ecologically divergent selection. This article is part of the theme issue 'Species' ranges in the face of changing environments (part I)'.

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

The evolutionary genomics of species' responses to climate change

Climate change is a threat to biodiversity. One way that this threat manifests is through pronounced shifts in the geographical range of species over time. To predict these shifts, researchers have primarily used species distribution models. However, these models are based on assumptions of niche conservatism and do not consider evolutionary processes, potentially limiting their accuracy and value. To incorporate evolution into the prediction of species' responses to climate change, researchers have turned to landscape genomic data and examined information about local genetic adaptation using climate models. Although this is an important advancement, this approach currently does not include other evolutionary processes-such as gene flow, population dispersal and genomic load-that are critical for predicting the fate of species across the landscape. Here, we briefly review the current practices for the use of species distribution models and for incorporating local adaptation. We next discuss the rationale and theory for considering additional processes, reviewing how they can be incorporated into studies of species' responses to climate change. We summarize with a conceptual framework of how manifold layers of information can be combined to predict the potential response of specific populations to climate change. We illustrate all of the topics using an exemplar dataset and provide the source code as potential tutorials. This Perspective is intended to be a step towards a more comprehensive integration of population genomics with climate change science.

Early evolution of reproductive isolation: A case of weak inbreeder/strong outbreeder leads to an intraspecific hybridization barrier in Arabidopsis lyrata

Reproductive strategies play a major role in plant speciation. Notably, transitions from outcrossing to selfing may lead to relaxed sexual selection and parental conflict. Shifts in mating systems can affect maternal and paternal interests, and thus parent-specific influence on endosperm development, leading to reproductive isolation: if selfing and outcrossing species hybridize, the resulting seeds may not be viable due to endosperm failure. Nevertheless, it remains unclear how the switch in mating systems can impact reproductive isolation between recently diverged lineages, that is, during the process of speciation. We investigated this question using Arabidopsis lyrata, which recently transitioned to selfing (10,000 years ago) in certain North American populations, where European populations remain outcrossing. We performed reciprocal crosses between selfers and outcrossers, and measured seed viability and endosperm development. We show that parental genomes in the hybrid seed negatively interact, as predicted by parental conflict. This leads to extensive hybrid seed lethality associated with endosperm cellularization disturbance. Our results suggest that this is primarily driven by divergent evolution of the paternal genome between selfers and outcrossers. In addition, we observed other hybrid seed defects, suggesting that sex-specific interests are not the only processes contributing to postzygotic reproductive isolation.

Adaptation across geographic ranges is consistent with strong selection in marginal climates and legacies of range expansion

Every species experiences limits to its geographic distribution. Some evolutionary models predict that populations at range edges are less well adapted to their local environments due to drift, expansion load, or swamping gene flow from the range interior. Alternatively, populations near range edges might be uniquely adapted to marginal environments. In this study, we use a database of transplant studies that quantify performance at broad geographic scales to test how local adaptation, site quality, and population quality change from spatial and climatic range centers toward edges. We find that populations from poleward edges perform relatively poorly, both on average across all sites (15% lower population quality) and when compared to other populations at home (31% relative fitness disadvantage), consistent with these populations harboring high genetic load. Populations from equatorial edges also perform poorly on average (18% lower population quality) but, in contrast, outperform foreign populations (16% relative fitness advantage), suggesting that populations from equatorial edges have strongly adapted to unique environments. Finally, we find that populations from sites that are thermally extreme relative to the species' niche demonstrate strong local adaptation, regardless of their geographic position. Our findings indicate that both nonadaptive processes and adaptive evolution contribute to variation in adaptation across species' ranges.

Maintenance of Adaptive Dynamics and No Detectable Load in a Range-Edge Outcrossing Plant Population

During range expansion, edge populations are expected to face increased genetic drift, which in turn can alter and potentially compromise adaptive dynamics, preventing the removal of deleterious mutations and slowing down adaptation. Here, we contrast populations of the European subspecies Arabidopsis lyrata ssp. petraea, which expanded its Northern range after the last glaciation. We document a sharp decline in effective population size in the range-edge population and observe that nonsynonymous variants segregate at higher frequencies. We detect a 4.9% excess of derived nonsynonymous variants per individual in the range-edge population, suggesting an increase of the genomic burden of deleterious mutations. Inference of the fitness effects of mutations and modeling of allele frequencies under the explicit demographic history of each population predicts a depletion of rare deleterious variants in the range-edge population, but an enrichment for fixed ones, consistent with the bottleneck effect. However, the demographic history of the range-edge population predicts a small net decrease in per-individual fitness. Consistent with this prediction, the range-edge population is not impaired in its growth and survival measured in a common garden experiment. We further observe that the allelic diversity at the self-incompatibility locus, which ensures strict outcrossing and evolves under negative frequency-dependent selection, has remained unchanged. Genomic footprints indicative of selective sweeps are broader in the Northern population but not less frequent. We conclude that the outcrossing species A. lyrata ssp. petraea shows a strong resilience to the effect of range expansion.

Drivers of linkage disequilibrium across a species' geographic range

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species’ range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Nearby genomic variants are often co-inherited because of limited recombination. The extent of non-random association of alleles at different loci is called linkage disequilibrium (LD) and is commonly used in genomic analyses, for example to detect regions under selection or to determine effective population size. Here we reversed testing and addressed how demographic history may affect LD within a species. Using genomic data from more than a thousand individuals of North American Arabidopsis lyrata from across the entire species’ range, we quantified the effect of postglacial range expansion and a shift in mating system from outcrossing to selfing on LD. We show that both factors lead to increased LD, and that the maximal effect of range expansion is comparable with a shift in mating system to selfing. Heightened LD involves deleterious mutations, and therefore, LD can also serve as an indicator of mutation accumulation. Furthermore, we provide evidence that some genes experienced stronger increases in LD possibly due to selection associated with the two demographic changes. Our results provide a novel and broad view on the evolutionary factors shaping LD that may also apply to the very many species that underwent postglacial range expansion.

Demography, genetic diversity and expansion load in the colonizing species Leontodon longirostris (Asteraceae) throughout its native range

Unravelling the evolutionary processes underlying range expansions is fundamental to understand the distribution of organisms, as well as to predict their future responses to environmental change. Predictions for range expansions include a loss of genetic diversity and an accumulation of deleterious alleles along the expansion axis, which can decrease fitness at the range-front (expansion load). In plants, empirical studies supporting expansion load are scarce, and its effects remain to be tested outside a few model species. Leontodon longirostris is a colonizing Asteraceae with a widespread distribution in the Western Mediterranean, providing a particularly interesting system to gain insight into the factors that can enhance or mitigate expansion load. In this study, we produced a first genome draft for the species, covering 418 Mbp (~53% of the genome). Although incomplete, this draft was suitable to design a targeted sequencing of ~1.5 Mbp in 238 L. longirostris plants from 21 populations distributed along putative colonization routes in the Iberian Peninsula. Inferred demographic history supports a range expansion from southern Iberia around 40,000 years ago, reaching northern Iberia around 25,000 years ago. The expansion was accompanied by a loss of genetic diversity and a significant increase in the proportion of putatively deleterious mutations. However, levels of expansion load in L. longirostris were smaller than those found in other plant species, which can be explained, at least partially, by its high dispersal ability, the self-incompatible mating system, and the fact that the expansion occurred along a strong environmental cline.

Polygenic adaptation of rosette growth in Arabidopsis thaliana

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its evolution in natural populations. Here, we investigate evolutionary and environmental factors shaping variation in the growth rate of Arabidopsis thaliana. We used plant diameter as a proxy to monitor plant growth over time in environments that mimicked latitudinal differences in the intensity of natural light radiation, across a set of 278 genotypes sampled within four broad regions, including an outgroup set of genotypes from China. A field experiment conducted under natural conditions confirmed the ecological relevance of the observed variation. All genotypes markedly expanded their rosette diameter when the light supply was decreased, demonstrating that environmental plasticity is a predominant source of variation to adapt plant size to prevailing light conditions. Yet, we detected significant levels of genetic variation both in growth rate and growth plasticity. Genome-wide association studies revealed that only 2 single nucleotide polymorphisms associate with genetic variation for growth above Bonferroni confidence levels. However, marginally associated variants were significantly enriched among genes with an annotated role in growth and stress reactions. Polygenic scores computed from marginally associated variants confirmed the polygenic basis of growth variation. For both light regimes, phenotypic divergence between the most distantly related population (China) and the various regions in Europe is smaller than the variation observed within Europe, indicating that the evolution of growth rate is likely to be constrained by stabilizing selection. We observed that Spanish genotypes, however, reach a significantly larger size than Northern European genotypes. Tests of adaptive divergence and analysis of the individual burden of deleterious mutations reveal that adaptive processes have played a more important role in shaping regional differences in rosette growth than maladaptive evolution.

The rate at which plants grow is a major functional trait in plant ecology. However, little is known about its genetic variation in natural populations. Here, we investigate genetic and environmental factors shaping variation in the growth rate of Arabidopsis thaliana and ask whether genetic variation in plant growth contributes to adaptation to local environmental conditions. We grew plants under two light regimes that mimic latitudinal differences in the intensity of natural light radiation, and measured plant diameter as it grew over time. When the light supply was decreased, plant diameter grew more slowly but reached a markedly larger final size, confirming that plants can adjust their growth to prevailing light conditions. Yet, we also detected significant levels of genetic variation both in growth rate and in how the growth dynamics is adjusted to the light conditions. We show that this variation is encoded by many loci of small effect that are hard to locate in the genome but overall significantly enriched among genes associated with growth and stress reactions. We further observe that Spanish genotypes tended to reach, on average, a significantly larger rosette size than Northern European genotypes. Tests of adaptive divergence indicate that these differences may reflect adaptation to local environmental conditions.

Shape of population interfaces as an indicator of mutational instability in coexisting cell populations

Cellular populations such as avascular tumors and microbial biofilms may 'invade' or grow into surrounding populations. The invading population is often comprised of a heterogeneous mixture of cells with varying growth rates. The population may also exhibit mutational instabilities, such as a heavy deleterious mutation load in a cancerous growth. We study the dynamics of a heterogeneous, mutating population competing with a surrounding homogeneous population, as one might find in a cancerous invasion of healthy tissue. We find that the shape of the population interface serves as an indicator for the evolutionary dynamics within the heterogeneous population. In particular, invasion front undulations become enhanced when the invading population is near a mutational meltdown transition or when the surrounding 'bystander' population is barely able to reinvade the mutating population. We characterize these interface undulations and the effective fitness of the heterogeneous population in one- and two-dimensional systems.

Demographic Processes Linked to Genetic Diversity and Positive Selection across a Species' Range

Demography determines the strength of genetic drift, which generally reduces genetic variation and the efficacy of selection. Here, we disentangled the importance of demographic processes at a local scale (census size and mating system) and at a species-range scale (old split between population clusters, recolonization after the last glaciation cycle, and admixture) in determining within-population genomic diversity and genomic signatures of positive selection. Analyses were based on re-sequence data from 52 populations of North American Arabidopsis lyrata collected across its entire distribution. The mating system and range dynamics since the last glaciation cycle explained around 60% of the variation in genomic diversity among populations and 52% of the variation in the signature of positive selection. Diversity was lowest in selfing compared with outcrossing populations and in areas further away from glacial refugia. In parallel, reduced positive selection was found in selfing populations and in populations with a longer route of postglacial range expansion. The signature of positive selection was also reduced in populations without admixture. We conclude that recent range expansion can have a profound influence on diversity in coding and non-coding DNA, similar in magnitude to the shift toward selfing. Distribution limits may in fact be caused by reduced effective population size and compromised positive selection in recently colonized parts of the range.

What Do We Really Know About Adaptation at Range Edges?

Recent theory and empirical evidence have provided new insights regarding how evolutionary forces interact to shape adaptation at stable and transient range margins. Predictions regarding trait divergence at leading edges are frequently supported. However, declines in fitness at and beyond edges show that trait divergence has sometimes been insufficient to maintain high fitness, so identifying constraints to adaptation at range edges remains a key challenge. Indirect evidence suggests that range expansion may be limited by adaptive genetic variation, but direct estimates of genetic constraints at and beyond range edges are still scarce. Sequence data suggest increased genetic load in edge populations in several systems, but its causes and fitness consequences are usually poorly understood. The balance between maladaptive and positive effects of gene flow on fitness at range edges deserves further study. It is becoming increasingly clear that characterizations about degree of adaptation based solely on geographical peripherality are unsupported.

Westward range expansion from middle latitudes explains the Mississippi River discontinuity in a forest herb of eastern North America

It is often expected that temperate plants have expanded their geographical ranges northward from primarily southern refugia. Evidence for this hypothesis is mixed in eastern North American species, and there is increasing support for colonization from middle latitudes. We studied genome-wide patterns of variation in RADseq loci to test hypotheses concerning range expansion in a North American forest herb (Campanula americana). First, spatial patterns of genetic differentiation were determined. Then phylogenetic relationships and divergence times were estimated. Spatial signatures of genetic drift were also studied to identify the directionality of recent range expansion and its geographical origins. Finally, spatially explicit scenarios for the spread of plants across the landscape were compared, using variation in the population mutation parameter and Tajima's D. We found strong longitudinal subdivision, with populations clustering into groups west and east of the Mississippi River. While the southeastern region was probably part of a diverse Pleistocene refugium, there is little evidence that range expansion involved founders from these southern locales. Instead, declines in genetic diversity and the loss of rare alleles support a westward colonization wave from a middle latitude refugium near the southern Appalachian Mountains, with subsequent expansion from a Pleistocene staging ground in the Mississippi River Valley (0.51-1.27 million years ago). These analyses implicate stepping stone colonization from middle latitudes as an important mechanism of species range expansion in eastern North America. This study further demonstrates the utility of population genetics as a tool to infer the routes travelled by organisms during geographical range expansion.

Eco-evolutionary dynamics of range expansion

Understanding the movement of species' ranges is a classic ecological problem that takes on urgency in this era of global change. Historically treated as a purely ecological process, range expansion is now understood to involve eco-evolutionary feedbacks due to spatial genetic structure that emerges as populations spread. We synthesize empirical and theoretical work on the eco-evolutionary dynamics of range expansion, with emphasis on bridging directional, deterministic processes that favor evolved increases in dispersal and demographic traits with stochastic processes that lead to the random fixation of alleles and traits. We develop a framework for understanding the joint influence of these processes in changing the mean and variance of expansion speed and its underlying traits. Our synthesis of recent laboratory experiments supports the consistent role of evolution in accelerating expansion speed on average, and highlights unexpected diversity in how evolution can influence variability in speed: results not well predicted by current theory. We discuss and evaluate support for three classes of modifiers of eco-evolutionary range dynamics (landscape context, trait genetics, and biotic interactions), identify emerging themes, and suggest new directions for future work in a field that stands to increase in relevance as populations move in response to global change.

Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

Reconstruction of a species’ past demographic history from genetic data can highlight historical factors that have shaped the distribution of genetic diversity along its genome and its geographic range. Here, we combine genotyping-by-sequencing with demographic modelling to address these issues in the Coho salmon, a Pacific salmon of conservation concern in some parts of its range, notably in the south. Our demographic reconstructions reveal a linear decrease in genetic diversity toward the north of the species range, supporting the hypothesis of a northern route of postglacial recolonization from a single major southern refugium. As predicted by theory, we also observed a higher proportion of deleterious mutations in the most distant populations from this refugium. Beyond this general pattern, among-site variation in the proportion of deleterious mutations is consistent with different local trends in effective population sizes. Our results highlight the potential importance of understanding historical factors that have shaped geographic patterns of the distribution of deleterious mutations in order to implement effective management programs for the conservation of wild populations. Such fundamental knowledge of human historical demography is now having major impacts on health sciences, and we argue it is time to integrate such approaches in conservation science as well.

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.