Influence of learning on range expansion and adaptation to novel habitats

Learning has been postulated to ‘drive’ evolution, but its influence on adaptive evolution in heterogeneous environments has not been formally examined. We used a spatially explicit individual‐based model to study the effect of learning on the expansion and adaptation of a species to a novel habitat. Fitness was mediated by a behavioural trait (resource preference), which in turn was determined by both the genotype and learning. Our findings indicate that learning substantially increases the range of parameters under which the species expands and adapts to the novel habitat, particularly if the two habitats are separated by a sharp ecotone (rather than a gradient). However, for a broad range of parameters, learning reduces the degree of genetically‐based local adaptation following the expansion and facilitates maintenance of genetic variation within local populations. Thus, in heterogeneous environments learning may facilitate evolutionary range expansions and maintenance of the potential of local populations to respond to subsequent environmental changes.     

Long‐distance dispersal maximizes evolutionary potential during rapid geographic range expansion

Conventional wisdom predicts that sequential founder events will cause genetic diversity to erode in species with expanding geographic ranges, limiting evolutionary potential at the range margin. Here, we show that invasive European starlings (Sturnus vulgaris) in South Africa preserve genetic diversity during range expansion, possibly as a result of frequent long‐distance dispersal events. We further show that unfavourable environmental conditions trigger enhanced dispersal, as indicated by signatures of selection detected across the expanding range. This brings genetic variation to the expansion front, counterbalancing the cumulative effects of sequential founding events and optimizing standing genetic diversity and thus evolutionary potential at range margins during spread. Therefore, dispersal strategies should be highlighted as key determinants of the ecological and evolutionary performances of species in novel environments and in response to global environmental change.     

Fitness declines towards range limits and local adaptation to climate affect dispersal evolution during climate‐induced range shifts

Dispersal ability will largely determine whether species track their climatic niches during climate change, a process especially important for populations at contracting (low‐latitude/low‐elevation) range limits that otherwise risk extinction. We investigate whether dispersal evolution at contracting range limits is facilitated by two processes that potentially enable edge populations to experience and adjust to the effects of climate deterioration before they cause extinction: (i) climate‐induced fitness declines towards range limits and (ii) local adaptation to a shifting climate gradient. We simulate a species distributed continuously along a temperature gradient using a spatially explicit, individual‐based model. We compare range‐wide dispersal evolution during climate stability vs. directional climate change, with uniform fitness vs. fitness that declines towards range limits (RLs), and for a single climate genotype vs. multiple genotypes locally adapted to temperature. During climate stability, dispersal decreased towards RLs when fitness was uniform, but increased when fitness declined towards RLs, due to highly dispersive genotypes maintaining sink populations at RLs, increased kin selection in smaller populations, and an emergent fitness asymmetry that favoured dispersal in low‐quality habitat. However, this initial dispersal advantage at low‐fitness RLs did not facilitate climate tracking, as it was outweighed by an increased probability of extinction. Locally adapted genotypes benefited from staying close to their climate optima; this selected against dispersal under stable climates but for increased dispersal throughout shifting ranges, compared to cases without local adaptation. Dispersal increased at expanding RLs in most scenarios, but only increased at the range centre and contracting RLs given local adaptation to climate.     

Quantifying the effects of migration and mutation on adaptation and demography in spatially heterogeneous environments

How do mutation and gene flow influence population persistence, niche expansion and local adaptation in spatially heterogeneous environments? In this article, we analyse a demographic and evolutionary model of adaptation to an environment containing two habitats in equal frequencies, and we bridge the gap between different theoretical frameworks. Qualitatively, our model yields four qualitative types of outcomes: (i) global extinction of the population, (ii) adaptation to one habitat only, but also adaptation to both habitats with, (iii) specialized phenotypes or (iv) with generalized phenotypes, and we determine the conditions under which each equilibrium is reached. We derive new analytical approximations for the local densities and the distributions of traits in each habitat under a migration–selection–mutation balance, compute the equilibrium values of the means, variances and asymmetries of the local distributions of phenotypes, and contrast the effects of migration and mutation on the evolutionary outcome. We then check our analytical results by solving our model numerically, and also assess their robustness in the presence of demographic stochasticity. Although increased migration results in a decrease in local adaptation, mutation in our model does not influence the values of the local mean traits. Yet, both migration and mutation can have dramatic effects on population size and even lead to metapopulation extinction when selection is strong. Niche expansion, the ability for the population to adapt to both habitats, can also be prevented by small migration rates and a reduced evolutionary potential characterized by rare mutation events of small effects; however, niche expansion is otherwise the most likely outcome. Although our results are derived under the assumption of clonal reproduction, we finally show and discuss the links between our model and previous quantitative genetics models.     

Constraints on adaptation of Escherichia coli to mixed‐resource environments increase over time

Can a population evolved in two resources reach the same fitness in both as specialist populations evolved in each of the individual resources? This question is central to theories of ecological specialization, the maintenance of genetic variation, and sympatric speciation, yet relatively few experiments have examined costs of generalism over long‐term adaptation. We tested whether selection in environments containing two resources limits a population's ability to adapt to the individual resources by comparing the fitness of replicate Escherichia coli populations evolved for 6000 generations in the presence of glucose or lactose alone (specialists), or in varying presentations of glucose and lactose together (generalists). We found that all populations had significant fitness increases in both resources, though the magnitude and rate of these increases differed. For the first 4000 generations, most generalist populations increased in fitness as quickly in the individual resources as the corresponding specialist populations. From 5000 generations, however, a widespread cost of adaptation affected all generalists, indicating a growing constraint on their abilities to adapt to two resources simultaneously. Our results indicate that costs of generalism are prevalent, but may influence evolutionary trajectories only after a period of cost‐free adaptation.     

Local adaptation primes cold‐edge populations for range expansion but not warming‐induced range shifts

According to theory, edge populations may be poised to expand species’ ranges if they are locally adapted to extreme conditions, or ill‐suited to colonise beyond‐range habitat if their offspring are genetically and competitively inferior. We tested these contrasting predictions by transplanting low‐, mid‐, and high‐elevation (edge) populations of an annual plant throughout and above its elevational distribution. Seed from poor‐quality edge habitat (one of two transects) had inferior emergence, but edge seeds also had adaptive phenology (both transects). High‐elevation plants flowered earlier, required less heat accumulation to mature seed, and so achieved higher lifetime fitness at and above the range edge. Experimental warming improved fitness above the range, but eliminated the advantage of local cold‐edge populations, supporting recent models in which cold‐adapted edge populations do not facilitate warming‐induced range shifts. The highest above‐range fitness was achieved by a ‘super edge phenotype’ from a neighbouring mountain, suggesting key adaptations exist regionally even if absent from local edge populations.     

Identification of mega‐environments in Europe and effect of allelic variation at maturity E loci on adaptation of European soybean

Soybean cultivation holds great potential for a sustainable agriculture in Europe, but adaptation remains a central issue. In this large mega‐environment (MEV) study, 75 European cultivars from five early maturity groups (MGs 000–II) were evaluated for maturity‐related traits at 22 locations in 10 countries across Europe. Clustering of the locations based on phenotypic similarity revealed six MEVs in latitudinal direction and suggested several more. Analysis of maturity identified several groups of cultivars with phenotypic similarity that are optimally adapted to the different growing regions in Europe. We identified several haplotypes for the allelic variants at the E1, E2, E3 and E4 genes, with each E haplotype comprising cultivars from different MGs. Cultivars with the same E haplotype can exhibit different flowering and maturity characteristics, suggesting that the genetic control of these traits is more complex and that adaptation involves additional genetic pathways, for example temperature requirement. Taken together, our study allowed the first unified assessment of soybean‐growing regions in Europe and illustrates the strong effect of photoperiod on soybean adaptation and MEV classification, as well as the effects of the E maturity loci for soybean adaptation in Europe.     

Influence of Paleolithic range contraction,admixture and long‐distance dispersal on genetic gradients of modern humans in Asia

Cavalli‐Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east‐west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially‐explicit computer simulations of genetic data under the following scenarios: (a) variable levels of admixture between Paleolithic and Neolithic populations, (b) migration through long‐distance dispersal (LDD), (c) Paleolithic range contraction induced by the last glacial maximum (LGM), and (d) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analysed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east‐west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east‐west geographic axis of this continent.     

Life history traits variation in heterogeneous environment: The case of a freshwater snail resistance to pond drying

Ecologists and population geneticists have long suspected that the diversity of living organisms was connected to the structure of their environment. In heterogeneous environments, diversifying selection combined to restricted gene flow may indeed lead to locally adapted populations. The freshwater snail, Galba truncatula, is a good model to address this question because it is present in a heterogeneous environment composed of temporary and permanent waters. In order to test the selective importance of those environments, we proposed here to measure survival of lineages from both habitats during drought episodes. To this purpose, we experimentally submitted adults and juveniles individuals from both habitats to drought. We found a difference in desiccation resistance between temporary and permanents waters only for adults. Adults from temporary habitats were found more resistant to drought. This divergence in desiccation resistance seems to explain the unexpected life history traits differences between habitats observed.     

Breach of the northern Rocky Mountain geoclimatic barrier: initiation of range expansion by the mountain pine beetle

Aim Our aim is to examine the historical breach of the geoclimatic barrier of the Rocky Mountains by the mountain pine beetle (Dendroctonus ponderosae Hopkins). This recent range expansion from west of the North American continental divide into the eastern boreal forest threatens to provide a conduit to naïve pine hosts in eastern North America. We examine the initial expansion events and determine potential mechanism(s) of spread by comparing spread patterns in consecutive years to various dispersal hypotheses such as: (1) meso‐scale atmospheric dispersal of insects from source populations south‐west of the Rocky Mountains in British Columbia (i.e. their historical range), (2) anthropogenic transport of infested plant material, and (3) spread of insect populations across adjacent stands via corridors of suitable habitat. Location British Columbia, Canada. Methods We explore potential mechanism(s) of invasion of the mountain pine beetle using spatial point process models for the initial 3 years of landscape‐level data collection, 2004–2006. Specifically, we examine observed patterns of infestation relative to covariates reflecting various dispersal hypotheses. We select the most parsimonious models for each of the initial 3 years of invasion using information criteria statistics. Results The initial range expansion and invasion of the beetle was characterized by aerial deposition along a strong north‐west to south‐east gradient, with additional aerial deposition and localized dispersal from persisting populations in following years. Main conclusions Following deposition of a wave front of mountain pine beetles parallel to the Rocky Mountains via meso‐scale atmospheric dispersal, the areas of highest intensity of infestations advanced up to 25 km north‐east towards jack pine (Pinus banksiana) habitat in a single year. There appeared to be no association between putative anthropogenic movement of infested materials and initial range expansion of the mountain pine beetle across the continental divide.     

The magnitude of local adaptation under genotype‐dependent dispersal

Dispersal moves individuals from patches where their immediate ancestors were successful to sites where their genotypes are untested. As a result, dispersal generally reduces fitness, a phenomenon known as “migration load.” The strength of migration load depends on the pattern of dispersal and can be dramatically lessened or reversed when individuals move preferentially toward patches conferring higher fitness. Evolutionary ecologists have long modeled nonrandom dispersal, focusing primarily on its effects on population density over space, the maintenance of genetic variation, and reproductive isolation. Here, we build upon previous work by calculating how the extent of local adaptation and the migration load are affected when individuals differ in their dispersal rate in a genotype‐dependent manner that alters their match to their environment. Examining a one‐locus, two‐patch model, we show that local adaptation occurs through a combination of natural selection and adaptive dispersal. For a substantial portion of parameter space, adaptive dispersal can be the predominant force generating local adaptation. Furthermore, genetic load may be largely averted with adaptive dispersal whenever individuals move before selective deaths occur. Thus, to understand the mechanisms driving local adaptation, biologists must account for the extent and nature of nonrandom, genotype‐dependent dispersal, and the potential for adaptation via spatial sorting of genotypes.     

Estimating range expansion of wildlife in heterogeneous landscapes: A spatially explicit state‐space matrix model coupled with an improved numerical integration technique

Dispersal as well as population growth is a key demographic process that determines population dynamics. However, determining the effects of environmental covariates on dispersal from spatial‐temporal abundance proxy data is challenging owing to the complexity of model specification for directional dispersal permeability and the extremely high computational loads for numerical integration. In this paper, we present a case study estimating how environmental covariates affect the dispersal of Japanese sika deer by developing a spatially explicit state‐space matrix model coupled with an improved numerical integration technique (Markov chain Monte Carlo with particle filters). In particular, we explored the environmental drivers of inhomogeneous range expansion, characteristic of animals with short dispersal. Our model framework successfully reproduced the complex population dynamics of sika deer, including rapid changes in densely populated areas and distribution fronts within a decade. Furthermore, our results revealed that the inhomogeneous range expansion of sika deer seemed to be primarily caused by the dispersal process (i.e., movement barriers in fragmented forests) rather than population growth. Our state‐space matrix model enables the inference of population dynamics for a broad range of organisms, even those with low dispersal ability, in heterogeneous landscapes, and could address many pressing issues in conservation biology and ecosystem management.     

The effect of frequency‐dependent selection on resistance and tolerance to herbivory

Negative frequency‐dependent selection (FDS), where rare genotypes are favoured by selection, is commonly invoked as a mechanism explaining the maintenance of genetic variation in plant defences. However, empirical tests of FDS in plant–herbivore interactions are lacking. We evaluated whether the oviposition preference of the specialist herbivore Lema daturaphila is a mechanism through which this herbivore can exert FDS on its host plant Datura stramonium. The frequency of contrasting resistance–tolerance strategies was manipulated within experimental plots, and the plants were exposed to a similar initial density of their natural herbivore. Herbivore oviposition preference and final density, as well as plant damage and seed production, were estimated. Overall, we found that the high‐resistant–low‐tolerant genotypes produced four times more seeds when common than when rare, whereas the high‐tolerant–low‐resistant genotypes achieved twice its fitness when rare than when common. This pattern was the result of differential oviposition preferences. In addition, when the high‐resistant–low‐tolerant genotypes were common, there was a three‐fold decreased in herbivore final density which led to a decrease in damage level by 10%. Thus, in our experiment positive FDS seems to favour resistance over tolerance. We discuss how this result would change if the extent of herbivore local adaptation and damage modify the pattern of positive FDS acting on resistance and the optimal allocation to tolerance.     

Genetic diversity and structuring across the range of a widely distributed ladybird: focus on rear‐edge populations phenotypically divergent

Population genetics and phenotypic structures are often predicted to vary along the geographic range of a species. This phenomenon would be accentuated for species with large range areas, with discontinuities and marginal populations. We herein compare the genetic patterns of central populations of Coccinella septempunctata L. with those of two phenotypically differentiated populations considered as rear‐edge populations and subspecies based on phenotype (Algeria and Japan). According to the central‐marginal model and expected characteristics of rear‐edge populations, we hypothesize that these rear‐edge populations have (1) a reduced genetic diversity, resulting from their relative isolation over long periods of time, (2) a higher population genetic differentiation, explained by low contemporary gene flow levels, and (3) a relationship between genetic diversity characteristics and phenotypes, due to historical isolation and/or local adaptation. Based on genotyping of 28 populations for 18 microsatellite markers, several levels of regional genetic diversity and differentiation are observed between and within populations, according to their localization: low within‐population genetic diversity and higher genetic differentiation of rear‐edge populations. The genetic structuring clearly dissociates the Algerian and Eastern Asia populations from the others. Geographical patterns of genetic diversity and differentiation support the hypothesis of the central‐marginal model. The pattern observed is in agreement with the phenotypic structure across species range. A clear genetic break between populations of Algeria, the Eastern Asia, and the remaining populations is a dominant feature of the data. Differential local adaptations, absence of gene flow between marginal and central populations, and/or incapacity to mate after colonization, have contributed to their distinct genotypic and phenotypic characteristics.     

Inferring the geographic origin of a range expansion: Latitudinal and longitudinal coordinates inferred from genomic data in an ABC framework with the program x‐origin

Climatic or environmental change is not only driving distributional shifts in species today, but it has also caused distributions to expand and contract in the past. Inferences about the geographic locations of past populations especially regions that served as refugia (i.e., source populations) and migratory routes are a challenging endeavour. Refugial areas may be evidenced from fossil records or regions of temporal stability inferred from ecological niche models. Genomic data offer an alternative and broadly applicable source of information about the locality of refugial areas, especially relative to fossil data, which are either unavailable or incomplete for most species. Here, we present a pipeline we developed (called x ‐origin ) for statistically inferring the geographic origin of range expansion using a spatially explicit coalescent model and an approximate Bayesian computation testing framework. In addition to assessing the probability of specific latitudinal and longitudinal coordinates of refugial or source populations, such inferences can also be made accounting for the effects of temporal and spatial environmental heterogeneity, which may impact migration routes. We demonstrate x ‐origin with an analysis of genomic data collected in the Collared pika that underwent postglacial expansion across Alaska, as well as present an assessment of its accuracy under a known model of expansion to validate the approach.