Local adaptation in a changing world: the roles of gene-flow, mutation, and sexual reproduction

Out to sea: ocean currents and patterns of asymmetric gene flow in an intertidal fish species

Passive dispersal via wind or ocean currents can drive asymmetric gene flow, which influences patterns of genetic variation and the capacity of populations to evolve in response to environmental change. The mangrove rivulus fish (Kryptolebias marmoratus), hereafter "rivulus," is an intertidal fish species restricted to the highly fragmented New World mangrove forests of Central America, the Caribbean, the Bahamas, and Florida. Mangrove patches are biological islands with dramatic differences in both abiotic and biotic conditions compared to adjacent habitat. Over 1,000 individual rivulus across 17 populations throughout its range were genotyped at 32 highly polymorphic microsatellites. Range-wide population genetic structure was evaluated with five complementary approaches that found eight distinct population clusters. However, an analysis of molecular variance indicated significant population genetic structure among regions, populations within regions, sampling locations within populations, and individuals within sampling locations, indicating that rivulus has both broad- and fine-scale genetic differentiation. Integrating range-wide genetic data with biophysical modeling based on 10 years of ocean current data showed that ocean currents and the distance between populations over water drive gene flow patterns on broad scales. Directional migration estimates suggested some significant asymmetries in gene flow that also were mediated by ocean currents and distance. Specifically, populations in the center of the range (Florida Keys) were identified as sinks that received migrants (and alleles) from other populations but failed to export individuals. These populations thus harbor genetic variation, perhaps even from extirpated populations across the range, but ocean currents and complex arrangements of landmasses might prevent the distribution of that genetic variation elsewhere. Hence, the inherent asymmetry of ocean currents shown to impact both genetic differentiation and directional migration rates may be responsible for the complex distribution of genetic variation across the range and observed patterns of metapopulation structure.

Effects of Vertically Heterogeneous Soil Salinity on Genetic Polymorphism and Productivity of the Widespread Halophyte Bassia prostrata

Salinity is one of the environmental factors that affects both productivity and genetic diversity in plant species. Within the soil profile, salinity is a dynamic indicator and significantly changes with depth. The present study examined the effects of the vertical heterogeneity of soil salinity chemistry on the plant height, fresh and dry biomass accumulation, water content, level of genetic polymorphism, and observed and expected heterozygosity in seven populations of halophyte Bassia prostrata in natural habitats. Soil salinity ranged from slight (S_salts = 0.11-0.25%) to extreme (S_salts = 1.35-2.57%). The main contributors to salinity were Na⁺, Ca²⁺, and Mg²⁺. Multivariate analysis revealed that biomass accumulation is positively affected by moderate/high salinity in 20-60 cm soil layers, which may be associated with the salt required for the optimal growth of the halophyte B. prostrata. The formation of seed genetic diversity is negatively affected by slight/moderate salinity in the 0-40 cm layers. An increase in divalent ion content can reduce genetic diversity and increase the local adaptation of B. prostrata to magnesium-calcium sulfate salinity. The effect of the in-depth distribution of soil salinity on productivity and genetic diversity may be related to seasonal variables during biomass accumulation (summer) and seed formation (autumn).

Unravelling the causes and consequences of dispersal syndromes in a wild passerine

Evidence accumulates that dispersal is correlated with individual behavioural phenotype (dispersal syndrome). The evolutionary causes and consequences of such covariation depend on the degree of plasticity versus inheritance of the traits, which requires challenging experiments to implement in mobile organisms. Here, we combine a forced dispersal experiment, natural colonization and longitudinal data to establish if dispersal and aggression levels are integrated and to test their adaptive nature in pied flycatchers (Ficedula hypoleuca). We found that (forced) dispersers behaved more aggressively in their first breeding year after dispersal and decreased their aggression in following years. Strength of dispersal syndrome and direction of fecundity selection on aggression in newly colonized areas varied between years. We propose that the net benefits of aggression for dispersers increase under harsh conditions (e.g. low food abundance). This hypothesis now warrants further testing. Overall, this study provides unprecedented experimental evidence that dispersal syndromes can be remodelled via adaptive plasticity depending on the individuals' local breeding experience and/or year-specific ecological conditions. It highlights the importance of individual behavioural variation in population dynamics.

Plant functional traits reflect different dimensions of species invasiveness

Trait-based invasiveness studies typically categorize exotic species as invasive or noninvasive, implicitly assuming species form two homogenous groups. However, species can become invasive in different ways (e.g., high abundance, fast spread), likely relying on different functional traits to do so. As such, binary classification may obscure traits associated with invasiveness. We tested whether (1) the way in which invasiveness is quantified influences its correlation with functional traits and (2) different demography-based metrics are related to different sets of traits. Using a case study of 251 herbs exotic to Victoria, Australia, we quantified species' invasiveness using 10 metrics: four continuous, demography-based dimensions of invasiveness (spread rate, local abundance, geographic and environmental range sizes) and six binary classifications of invasiveness (based on alternative sources and invasion criteria). We examined the correlation between species' invasiveness and a set of four traits known to relate to plant demography and invasion. Then, we examined whether different demographic dimensions of invasiveness were better explained by different sets of traits. We found that the way invasiveness was quantified was important: different traits were linked with different invasiveness metrics, and some traits showed opposite effects across metrics. Species with fast spread were either tall with small seeds (i.e., good colonizers), or had heavy, animal-dispersed seeds. Plants with a large environmental range had greater plasticity for some traits. Locally abundant plants had low SLA and heavy seeds (i.e., strong competitors). Animal dispersal was also key to reach a large geographic range. No traits were consistently related to the six binary classifications. Our results indicate that exotic plants are invasive in different ways and rely on different combinations of traits to be so. Some traits (e.g., seed mass) had complex relationships with invasion: they apparently promote, hampered, or had no influence on different dimensions of invasiveness. Our findings are consistent with the notion that plant species use strategies that may be near optimal under some, but not all, ecological conditions. Compared to binary classifications of invasiveness, the use of invasiveness dimensions advances clearer hypothesis testing in invasion science.

A general mathematical method for predicting spatio-temporal correlations emerging from agent-based models

Agent-based models are used to study complex phenomena in many fields of science. While simulating agent-based models is often straightforward, predicting their behaviour mathematically has remained a key challenge. Recently developed mathematical methods allow the prediction of the emerging spatial patterns for a general class of agent-based models, whereas the prediction of spatio-temporal pattern has been thus far achieved only for special cases. We present a general and mathematically rigorous methodology that allows deriving the spatio-temporal correlation structure for a general class of individual-based models. To do so, we define an auxiliary model, in which each agent type of the primary model expands to three types, called the original, the past and the new agents. In this way, the auxiliary model keeps track of both the initial and current state of the primary model, and hence the spatio-temporal correlations of the primary model can be derived from the spatial correlations of the auxiliary model. We illustrate the agreement between analytical predictions and agent-based simulations using two example models from theoretical ecology. In particular, we show that the methodology is able to correctly predict the dynamical behaviour of a host–parasite model that shows spatially localized oscillations.

Signatures of selection in a recent invasion reveal adaptive divergence in a highly vagile invasive species

A detailed understanding of population genetics in invasive populations helps us to identify drivers of successful alien introductions. Here, we investigate putative signals of selection in Australian populations of invasive common starlings, Sturnus vulgaris, and seek to understand how these have been influenced by introduction history. We used reduced representation sequencing to determine population structure, and identify Single Nucleotide Polymorphisms (SNPs) that are putatively under selection. We found that since their introduction into Australia, starling populations have become genetically differentiated despite the potential for high levels of dispersal, and that starlings have responded to selective pressures imposed by a wide range of environmental conditions across their geographic range. Isolation by distance appears to have played a strong role in determining genetic substructure across the starling's Australian range. Analyses of candidate SNPs that are putatively under selection indicated that aridity, precipitation and temperature may be important factors driving adaptive variation across the starling's invasive range in Australia. However, we also noted that the historic introduction regime may leave footprints on sites flagged as being under adaptive selection, and encourage critical interpretation of selection analyses in non-native populations.

Predicting Thermal Adaptation by Looking Into Populations' Genomic Past

Molecular evolution offers an insightful theory to interpret the genomic consequences of thermal adaptation to previous events of climate change beyond range shifts. However, disentangling often mixed footprints of selective and demographic processes from those due to lineage sorting, recombination rate variation, and genomic constrains is not trivial. Therefore, here we condense current and historical population genomic tools to study thermal adaptation and outline key developments (genomic prediction, machine learning) that might assist their utilization for improving forecasts of populations' responses to thermal variation. We start by summarizing how recent thermal-driven selective and demographic responses can be inferred by coalescent methods and in turn how quantitative genetic theory offers suitable multi-trait predictions over a few generations via the breeder's equation. We later assume that enough generations have passed as to display genomic signatures of divergent selection to thermal variation and describe how these footprints can be reconstructed using genome-wide association and selection scans or, alternatively, may be used for forward prediction over multiple generations under an infinitesimal genomic prediction model. Finally, we move deeper in time to comprehend the genomic consequences of thermal shifts at an evolutionary time scale by relying on phylogeographic approaches that allow for reticulate evolution and ecological parapatric speciation, and end by envisioning the potential of modern machine learning techniques to better inform long-term predictions. We conclude that foreseeing future thermal adaptive responses requires bridging the multiple spatial scales of historical and predictive environmental change research under modern cohesive approaches such as genomic prediction and machine learning frameworks.

An evolutionary perspective on marine invasions

Species distributions are rapidly changing as human globalization increasingly moves organisms to novel environments. In marine systems, species introductions are the result of a number of anthropogenic mechanisms, notably shipping, aquaculture/mariculture, the pet and bait trades, and the creation of canals. Marine invasions are a global threat to human and non-human populations alike and are often listed as one of the top conservation concerns worldwide, having ecological, evolutionary, and social ramifications. Evolutionary investigations of marine invasions can provide crucial insight into an introduced species' potential impacts in its new range, including: physiological adaptation and behavioral changes to exploit new environments; changes in resident populations, community interactions, and ecosystems; and severe reductions in genetic diversity that may limit evolutionary potential in the introduced range. This special issue focuses on current research advances in the evolutionary biology of marine invasions and can be broadly classified into a few major avenues of research: the evolutionary history of invasive populations, post-invasion reproductive changes, and the role of evolution in parasite introductions. Together, they demonstrate the value of investigating marine invasions from an evolutionary perspective, with benefits to both fundamental and applied evolutionary biology at local and broad scales.

Simulating the effects of local adaptation and life history on the ability of plants to track climate shifts

Many studies have examined the impact of dispersal on local adaptation, but much less attention has been paid to how local adaptation influences range shifts. The aim of this study was to test how local adaptation might affect climate-driven range shifts in plants, and if this might differ between plants with different life histories. Simulated range shift dynamics were compared for hypothetical annual, perennial and tree species, each comprised of either one plastic genotype or six locally adapted genotypes. The landscape consists of shifting climate bands made up of 20 × 20 m patches containing multiple individuals. Effects of seed dispersal, breadth of the plastic species' tolerance, steepness of the climate gradient and rate of the climate shift are also examined. Local adaptation increased the equilibrium range size and aided range shifts by boosting fitness near range edges. However, when the rate of climate change was doubled on a steep gradient, locally adapted trees exhibited a higher percent loss of range during the climate shift. The plastic annual species with short dispersal was unable to recover its range size even after the climate stabilized, while the locally adapted annuals tracked climate change well. The results suggest that in most situations local adaptation and longer dispersal distances will be advantageous, though not necessarily sufficient, for tracking suitable climates. However, local adaptation might put species with long generation times at greater risk when climate shifts are very rapid. If confirmed by empirical tests, these results suggest that identifying variation between species in how fitness varies along climate gradients and in these key demographic rates might aid in prioritizing management actions.

Mitotic Recombination and Rapid Genome Evolution in the Invasive Forest Pathogen Phytophthora ramorum

Alien species are often successful invaders in new environments, despite the introduction of a few isolates with a reduced genetic pool. This is called the genetic paradox of invasion. We found two mechanisms by which the invasive forest pathogen causing sudden oak and sudden larch death can evolve. Extensive mitotic recombination producing runs of homozygosity generates genotypic diversity even in the absence of sexual reproduction, and rapid turnover of genes in the non-core, or nonessential portion of genome not shared by all isolates, allows pathogenicity genes to evolve rapidly or be eliminated while retaining essential genes. Mitotic recombination events occur in genomic hot spots, resulting in similar ROH patterns in different isolates or groups; one ROH, independently generated in two different groups, was enriched in pathogenicity genes and may be a target for selection. This provides important insights into the evolution of invasive alien pathogens and their potential for adaptation and future persistence.

Invasive alien species often have reduced genetic diversity and must adapt to new environments. Given the success of many invasions, this is sometimes called the genetic paradox of invasion. Phytophthora ramorum is invasive, limited to asexual reproduction within four lineages, and presumed clonal. It is responsible for sudden oak death in the United States, sudden larch death in Europe, and ramorum blight in North America and Europe. We sequenced the genomes of 107 isolates to determine how this pathogen can overcome the invasion paradox. Mitotic recombination (MR) associated with transposons and low gene density has generated runs of homozygosity (ROH) affecting 2,698 genes, resulting in novel genotypic diversity within the lineages. One ROH enriched in effectors was fixed in the NA1 lineage. An independent ROH affected the same scaffold in the EU1 lineage, suggesting an MR hot spot and a selection target. Differences in host infection between EU1 isolates with and without the ROH suggest that they may differ in aggressiveness. Non-core regions (not shared by all lineages) had signatures of accelerated evolution and were enriched in putative pathogenicity genes and transposons. There was a striking pattern of gene loss, including all effectors, in the non-core EU2 genome. Positive selection was observed in 8.0% of RxLR and 18.8% of Crinkler effector genes compared with 0.9% of the core eukaryotic gene set. We conclude that the P. ramorum lineages are diverging via a rapidly evolving non-core genome and that the invasive asexual lineages are not clonal, but display genotypic diversity caused by MR.

Different Transcriptomic Responses to Thermal Stress in Heat-Tolerant and Heat-Sensitive Pacific Abalones Indicated by Cardiac Performance

The Pacific abalone Haliotis discus hannai is one of the most economically important mollusks in China. Even though it has been farmed in southern China for almost 20 years, summer mortality remains the most challengeable problem for Pacific abalone aquaculture recently. Here, we determined the different heat tolerance ability for five selective lines of H. discus hannai by measuring the cardiac performance and Arrhenius breakpoint temperature (ABT). The Red line (RL) and Yangxia line (YL) were determined as the most heat-sensitive and most heat-tolerant line, respectively. Heart rates for RL were significantly lower than those of the YL at the same temperature (p < 0.05). The differentially expressed genes (DEGs), which were enriched in several pathways including cardiac muscle contraction, glutathione metabolism and oxidative phosphorylation, were identified between RL and YL at control temperature (20°C) and heat stress temperature (28.5°C, the ABT of the RL) by RNA-seq method. In the RL, 3370 DEGs were identified between the control and the heat-stress temperature, while only 1351 DEGs were identified in YL between these two temperature tests. Most of these DEGs were enriched in the pathways such as protein processing in endoplasmic reticulum, nucleotide binding and oligomerization domain (NOD) like receptor signaling, and ubiquitin mediated proteolysis. Notably, the most heat-tolerant line YL used an effective heat-protection strategy based on moderate transcriptional changes and regulation on the expression of key genes.

Comparing the consequences of natural selection, adaptive phenotypic plasticity, and matching habitat choice for phenotype-environment matching, population genetic structure, and reproductive isolation in meta-populations

Organisms commonly experience significant spatiotemporal variation in their environments. In response to such heterogeneity, different mechanisms may act that enhance ecological performance locally. However, depending on the nature of the mechanism involved, the consequences for populations may differ greatly. Building on a previous model that investigated the conditions under which different adaptive mechanisms (co)evolve, this study compares the ecological and evolutionary population consequences of three very different responses to environmental heterogeneity: matching habitat choice (directed gene flow), adaptive plasticity (associated with random gene flow), and divergent natural selection. Using individual‐based simulations, we show that matching habitat choice can have a greater adaptive potential than plasticity or natural selection: it allows for local adaptation while protecting genetic polymorphism despite global mating or strong environmental changes. Our simulations further reveal that increasing environmental fluctuations and unpredictability generally favor the emergence of specialist genotypes but that matching habitat choice is better at preventing local maladaptation by individuals. This confirms that matching habitat choice can speed up the genetic divergence among populations, cause indirect assortative mating via spatial clustering, and hence even facilitate sympatric speciation. This study highlights the potential importance of directed dispersal in local adaptation and speciation, stresses the difficulty of deriving its operation from nonexperimental observational data alone, and helps define a set of ecological conditions which should favor its emergence and subsequent detection in nature.

Population genetics and adaptation to climate along elevation gradients in invasive Solidago canadensis

Gene flow between populations may either support local adaptation by supplying genetic variation on which selection may act, or counteract it if maladapted alleles arrive faster than can be purged by selection. Although both such effects have been documented within plant species’ native ranges, how the balance of these forces influences local adaptation in invasive plant populations is less clear, in part because introduced species often have lower genetic variation initially but also tend to have good dispersal abilities. To evaluate the extent of gene flow and adaptation to local climate in invasive populations of Solidago canadensis, and the implications of this for range expansion, we compared population differentiation at microsatellite and chloroplast loci for populations across Switzerland and assessed the effect of environmental transfer distance using common gardens. We found that while patterns of differentiation at neutral genetic markers suggested that populations are connected through extensive pollen and seed movement, common-garden plants nonetheless exhibited modest adaptation to local climate conditions. Growth rate and flower production declined with climatic distance from a plant's home site, with clones from colder home sites performing better at or above the range limit. Such adaptation in invasive species is likely to promote further spread, particularly under climate change, as the genotypes positioned near the range edge may be best able to take advantage of lengthening growing seasons to expand the range.

Environmental and anthropogenic drivers of connectivity patterns: A basis for prioritizing conservation efforts for threatened populations

Ecosystem fragmentation and habitat loss have been the focus of landscape management due to restrictions on contemporary connectivity and dispersal of populations. Here, we used an individual approach to determine the drivers of genetic differentiation in caribou of the Canadian Rockies. We modelled the effects of isolation by distance, landscape resistance and predation risk and evaluated the consequences of individual migratory behaviour (seasonally migratory vs. sedentary) on gene flow in this threatened species. We applied distance-based and reciprocal causal modelling approaches, testing alternative hypotheses on the effects of geographic, topographic, environmental and local population-specific variables on genetic differentiation and relatedness among individuals. Overall, gene flow was restricted to neighbouring local populations, with spatial coordinates, local population size, groups and elevation explaining connectivity among individuals. Landscape resistance, geographic distances and predation risk were correlated with genetic distances, with correlations threefold higher for sedentary than for migratory caribou. As local caribou populations are increasingly isolated, our results indicate the need to address genetic connectivity, especially for populations with individuals displaying different migratory behaviours, whilst maintaining quality habitat both within and across the ranges of threatened populations.

What drivers phenotypic divergence in Leymus chinensis (Poaceae) on large-scale gradient, climate or genetic differentiation?

Elucidating the driving factors among-population divergence is an important task in evolutionary biology, however the relative contribution from natural selection and neutral genetic differentiation has been less debated. A manipulation experiment was conducted to examine whether the phenotypic divergence of Leymus chinensis depended on climate variations or genetic differentiations at 18 wild sites along a longitudinal gradient from 114 to 124°E in northeast China and at common garden condition of transplantation. Demographical, morphological and physiological phenotypes of 18 L. chinensis populations exhibited significant divergence along the gradient, but these divergent variations narrowed significantly at the transplantation. Moreover, most of the phenotypes were significantly correlated with mean annual precipitation and temperature in wild sites, suggesting that climatic variables played vital roles in phenotypic divergence of the species. Relative greater heterozygosity (H_E), genotype evenness (E) and Shannon-Wiener diversity (I) in western group of populations suggested that genetic differentiation also drove phenotypic divergence of the species. However, neutral genetic differentiation (F_ST = 0.041) was greatly lower than quantitative differentiation (Q_ST = 0.199), indicating that divergent selection/climate variable was the main factor in determining the phenotypic divergence of the species along the large-scale gradient.

Expression Differentiation Is Constrained to Low-Expression Proteins over Ecological Timescales

Protein expression level is one of the strongest predictors of protein sequence evolutionary rate, with high-expression protein sequences evolving at slower rates than low-expression protein sequences largely because of constraints on protein folding and function. Expression evolutionary rates also have been shown to be negatively correlated with expression level across human and mouse orthologs over relatively long divergence times (i.e., ∼100 million years). Long-term evolutionary patterns, however, often cannot be extrapolated to microevolutionary processes (and vice versa), and whether this relationship holds for traits evolving under directional selection within a single species over ecological timescales (i.e., <5000 years) is unknown and not necessarily expected. Expression is a metabolically costly process, and the expression level of a particular protein is predicted to be a tradeoff between the benefit of its function and the costs of its expression. Selection should drive the expression level of all proteins close to values that maximize fitness, particularly for high-expression proteins because of the increased energetic cost of production. Therefore, stabilizing selection may reduce the amount of standing expression variation for high-expression proteins, and in combination with physiological constraints that may place an upper bound on the range of beneficial expression variation, these constraints could severely limit the availability of beneficial expression variants. To determine whether rapid-expression evolution was restricted to low-expression proteins owing to these constraints on highly expressed proteins over ecological timescales, we compared venom protein expression levels across mainland and island populations for three species of pit vipers. We detected significant differentiation in protein expression levels in two of the three species and found that rapid-expression differentiation was restricted to low-expression proteins. Our results suggest that various constraints on high-expression proteins reduce the availability of beneficial expression variants relative to low-expression proteins, enabling low-expression proteins to evolve and potentially lead to more rapid adaptation.

Analysis of population genetic structure and gene flow in an annual plant before and after a rapid evolutionary response to drought

The impact of environmental change on population structure is not well understood. This study aimed to examine the effect of a climate change event on gene flow over space and time in two populations of Brassica rapa that evolved more synchronous flowering times over 5 years of drought in southern California. Using plants grown from seeds collected before and after the drought, we estimated genetic parameters within and between populations and across generations. We expected that with greater temporal opportunity to cross-pollinate, due to reduced phenological isolation, these populations would exhibit an increase in gene flow following the drought. We found low but significant FST, but no change in FST or Nm across the drought, in contrast to predictions. Bayesian analysis of these data indicates minor differentiation between the two populations but no noticeable change in structure before and after the shift in flowering times. However, we found high and significant levels of FIS, indicating that inbreeding likely occurred in these populations despite self-incompatibility in B. rapa. In this system, we did not find an impact of climate change on gene flow or population structuring. The contribution of gene flow to adaptive evolution may vary by system, however, and is thus an important parameter to consider in further studies of natural responses to environmental change.

Host ecotype generates evolutionary and epidemiological divergence across a pathogen metapopulation

The extent and speed at which pathogens adapt to host resistance varies considerably. This presents a challenge for predicting when--and where--pathogen evolution may occur. While gene flow and spatially heterogeneous environments are recognized to be critical for the evolutionary potential of pathogen populations, we lack an understanding of how the two jointly shape coevolutionary trajectories between hosts and pathogens. The rust pathogen Melampsora lini infects two ecotypes of its host plant Linum marginale that occur in close proximity yet in distinct populations and habitats. In this study, we found that within-population epidemics were different between the two habitats. We then tested for pathogen local adaptation at host population and ecotype level in a reciprocal inoculation study. Even after controlling for the effect of spatial structure on infection outcome, we found strong evidence of pathogen adaptation at the host ecotype level. Moreover, sequence analysis of two pathogen infectivity loci revealed strong genetic differentiation by host ecotype but not by distance. Hence, environmental variation can be a key determinant of pathogen population genetic structure and coevolutionary dynamics and can generate strong asymmetry in infection risks through space.

ALADYN - a spatially explicit, allelic model for simulating adaptive dynamics

ALADYN is a freely available cross-platform C++ modeling framework for stochastic simulation of joint allelic and demographic dynamics of spatially-structured populations. Juvenile survival is linked to the degree of match between an individual's phenotype and the local phenotypic optimum. There is considerable flexibility provided for the demography of the considered species and the genetic architecture of the traits under selection. ALADYN facilitates the investigation of adaptive processes to spatially and/or temporally changing conditions and the resulting niche and range dynamics. To our knowledge ALADYN is so far the only model that allows a continuous resolution of individuals' locations in a spatially explicit landscape together with the associated patterns of selection.

Genomic variability as a driver of plant-pathogen coevolution?

Pathogens apply one of the strongest selective pressures in plant populations. Understanding plant-pathogen coevolution has therefore been a major research focus for at least sixty years [1]. Recent comparative genomic studies have revealed that the genes involved in plant defense and pathogen virulence are among the most polymorphic in the respective genomes. Which fraction of this diversity influences the host-pathogen interaction? Do coevolutionary dynamics maintain variation? Here we review recent literature on the evolutionary and molecular processes that shape this variation, focusing primarily on gene-for-gene interactions. In summarizing theoretical and empirical studies of the processes that shape this variation in natural plant and pathogen populations, we find a disconnect between the complexity of ecological interactions involving hosts and their myriad microbes, and the models that describe them.

Small-scale patterns in snowmelt timing affect gene flow and the distribution of genetic diversity in the alpine dwarf shrub Salix herbacea

Current threats to biodiversity, such as climate change, are thought to alter the within-species genetic diversity among microhabitats in highly heterogeneous alpine environments. Assessing the spatial organization and dynamics of genetic diversity within species can help to predict the responses of organisms to environmental change. In this study, we evaluated whether small-scale heterogeneity in snowmelt timing restricts gene flow between microhabitats in the common long-lived dwarf shrub Salix herbacea L. We surveyed 273 genets across 12 early- and late-snowmelt sites (that is, ridges and snowbeds) in the Swiss Alps for phenological variation over 2 years and for genetic variation using seven SSR markers. Phenological differentiation triggered by differences in snowmelt timing did not correlate with genetic differentiation between microhabitats. On the contrary, extensive gene flow appeared to occur between microhabitats and slightly less extensively among adjacent mountains. However, ridges exhibited significantly lower levels of genetic diversity than snowbeds, and patterns of effective population size (Ne) and migration (Nem) between microhabitats were strongly asymmetric, with ridges acting as sources and snowbeds as sinks. As no recent genetic bottlenecks were detected in the studied sites, this asymmetry is likely to reflect current meta-population dynamics of the species dominated by gene flow via seeds rather than ancient re-colonization after the last glacial period. Overall, our results suggest that seed dispersal prevents snowmelt-driven genetic isolation, and snowbeds act as sinks of genetic diversity. We discuss the consequences of such small-scale variation in gene flow and diversity levels for population responses to climate change.

Between migration load and evolutionary rescue: dispersal, adaptation and the response of spatially structured populations to environmental change

The evolutionary potential of populations is mainly determined by population size and available genetic variance. However, the adaptability of spatially structured populations may also be affected by dispersal: positively by spreading beneficial mutations across sub-populations, but negatively by moving locally adapted alleles between demes. We develop an individual-based, two-patch, allelic model to investigate the balance between these opposing effects on a population's evolutionary response to rapid climate change. Individual fitness is controlled by two polygenic traits coding for local adaptation either to the environment or to climate. Under conditions of selection that favour the evolution of a generalist phenotype (i.e. weak divergent selection between patches) dispersal has an overall positive effect on the persistence of the population. However, when selection favours locally adapted specialists, the beneficial effects of dispersal outweigh the associated increase in maladaptation for a narrow range of parameter space only (intermediate selection strength and low linkage among loci), where the spread of beneficial climate alleles is not strongly hampered by selection against non-specialists. Given that local selection across heterogeneous and fragmented landscapes is common, the complex effect of dispersal that we describe will play an important role in determining the evolutionary dynamics of many species under rapidly changing climate.

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high- and low-altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high- and low-altitude sites using three common temperature treatments. Local adaptation was assessed using Q(ST)-F(ST) analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high-altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.

Invasion dynamics of a fish-free landscape by brown trout (Salmo trutta)

Metapopulation dynamics over the course of an invasion are usually difficult to grasp because they require large and reliable data collection, often unavailable. The invasion of the fish-free freshwater ecosystems of the remote sub-Antarctic Kerguelen Islands following man-made introductions of brown trout (Salmo trutta) in the 1950's is an exception to this rule. Benefiting from a full long term environmental research monitoring of the invasion, we built a Bayesian dynamic metapopulation model to analyze the invasion dynamics of 85 river systems over 51 years. The model accounted for patch size (river length and connections to lakes), alternative dispersal pathways between rivers, temporal trends in dynamics, and uncertainty in colonization date. The results show that the model correctly represents the observed pattern of invasion, especially if we assume a coastal dispersal pathway between patches. Landscape attributes such as patch size influenced the colonization function, but had no effect on propagule pressure. Independently from patch size and distance between patches, propagule pressure and colonization function were not constant through time. Propagule pressure increased over the course of colonization, whereas the colonization function decreased, conditional on propagule pressure. The resulting pattern of this antagonistic interplay is an initial rapid invasion phase followed by a strong decrease in the invasion rate. These temporal trends may be due to either adaptive processes or environmental gradients encountered along the colonization front. It was not possible to distinguish these two hypotheses. Because invasibility of Kerguelen Is. freshwater ecosystems is very high due to the lack of a pre-existing fish fauna and minimal human interference, our estimates of invasion dynamics represent a blueprint for the potential of brown trout invasiveness in pristine environments. Our conclusions shed light on the future of polar regions where, because of climate change, fish-free ecosystems become increasingly accessible to invasion by fish species.

Evolutionary branching in complex landscapes

Divergent adaptation to different environments can promote speciation, and it is thus important to consider spatial structure in models of speciation. Earlier theoretical work, however, has been limited to particularly simple types of spatial structure (linear environmental gradients and spatially discrete metapopulations), leaving unaddressed the effects of more realistic patterns of landscape heterogeneity, such as nonlinear gradients and spatially continuous patchiness. To elucidate the consequences of such complex landscapes, we adapt an established spatially explicit individual-based model of evolutionary branching. We show that branching is most probable at intermediate levels of various types of heterogeneity and that different types of heterogeneity have, to some extent, additive effects in promoting branching. In contrast to such additivity, we find a novel refugium effect in which refugia in hostile environments provide opportunities for colonization, thus increasing the probability of branching in patchy landscapes. Effects of patchiness depend on the scale of patches relative to dispersal. Providing a needed connection to empirical research on biodiversity and conservation policy, we introduce empirically accessible spatial environmental metrics that quantitatively predict a landscape's branching propensity.

Limited evolutionary rescue of locally adapted populations facing climate change

Dispersal is a key determinant of a population's evolutionary potential. It facilitates the propagation of beneficial alleles throughout the distributional range of spatially outspread populations and increases the speed of adaptation. However, when habitat is heterogeneous and individuals are locally adapted, dispersal may, at the same time, reduce fitness through increasing maladaptation. Here, we use a spatially explicit, allelic simulation model to quantify how these equivocal effects of dispersal affect a population's evolutionary response to changing climate. Individuals carry a diploid set of chromosomes, with alleles coding for adaptation to non-climatic environmental conditions and climatic conditions, respectively. Our model results demonstrate that the interplay between gene flow and habitat heterogeneity may decrease effective dispersal and population size to such an extent that substantially reduces the likelihood of evolutionary rescue. Importantly, even when evolutionary rescue saves a population from extinction, its spatial range following climate change may be strongly narrowed, that is, the rescue is only partial. These findings emphasize that neglecting the impact of non-climatic, local adaptation might lead to a considerable overestimation of a population's evolvability under rapid environmental change.

Contemporary divergence in early life history in grayling (Thymallus thymallus)

Background

Following colonization of new habitats and subsequent selection, adaptation to environmental conditions might be expected to be rapid. In a mountain lake in Norway, Lesjaskogsvatnet, more than 20 distinct spawning demes of grayling have been established since the lake was colonized, some 20-25 generations ago. The demes spawn in tributaries consistently exhibiting either colder or warmer temperature conditions during spawning in spring and subsequent early development during early summer. In order to explore the degree of temperature-related divergence in early development, a multi-temperature common-garden experiment was performed on embryos from four different demes experiencing different spring temperatures.

Results

Early developmental characters were measured to test if individuals from the four demes respond differently to the treatment temperatures. There was clear evidence of among-deme differences (genotype - environment interactions) in larval growth and yolk-to-body-size conversion efficiency. Under the cold treatment regime, larval growth rates were highest for individuals belonging to cold streams. Individuals from warm streams had the highest yolk-consumption rate under cold conditions. As a consequence, yolk-to-body-mass conversion efficiency was highest for cold-deme individuals under cold conditions. As we observed response parallelism between individuals from demes belonging to similar thermal groups for these traits, some of the differentiation seems likely to result from local adaptation

Conclusion

The observed differences in length at age during early larval development most likely have a genetic component, even though both directional and random processes are likely to have influenced evolutionary change in the demes under study.

Applications and implications of neutral versus non-neutral markers in molecular ecology

The field of molecular ecology has expanded enormously in the past two decades, largely because of the growing ease with which neutral molecular genetic data can be obtained from virtually any taxonomic group. However, there is also a growing awareness that neutral molecular data can provide only partial insight into parameters such as genetic diversity, local adaptation, evolutionary potential, effective population size, and taxonomic designations. Here we review some of the applications of neutral versus adaptive markers in molecular ecology, discuss some of the advantages that can be obtained by supplementing studies of molecular ecology with data from non-neutral molecular markers, and summarize new methods that are enabling researchers to generate data from genes that are under selection.