NATURAL SELECTION IN POPULATIONS SUBJECT TO A MIGRATION LOAD

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.

Evolution repeats itself in replicate long-term studies in the wild

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Climatic and soil characteristics account for the genetic structure of the invasive cactus moth Cactoblastis cactorum, in its native range in Argentina

Background

Knowledge of the physical and environmental conditions that may limit the migration of invasive species is crucial to assess the potential for expansion outside their native ranges. The cactus moth, Cactoblastis cactorum, is native to South America (Argentina, Paraguay, Uruguay and Brazil) and has been introduced and invaded the Caribbean and southern United States, among other regions. In North America there is an ongoing process of range expansion threatening cacti biodiversity of the genus Opuntia and the commercial profits of domesticated Opuntia ficus-indica.

Methods

To further understand what influences the distribution and genetic structure of this otherwise important threat to native and managed ecosystems, in the present study we combined ecological niche modeling and population genetic analyses to identify potential environmental barriers in the native region of Argentina. Samples were collected on the host with the wider distribution range, O. ficus-indica.

Results

Significant genetic structure was detected using 10 nuclear microsatellites and 24 sampling sites. At least six genetic groups delimited by mountain ranges, salt flats and wetlands were mainly located to the west of the Dry Chaco ecoregion. Niche modeling supports that this region has high environmental suitability where the upper soil temperature and humidity, soil carbon content and precipitation were the main environmental factors that explain the presence of the moth. Environmental filters such as the upper soil layer may be critical for pupal survival and consequently for the establishment of populations in new habitats, whereas the presence of available hosts is a necessary conditions for insect survival, upper soil and climatic characteristics will determine the opportunities for a successful establishment.

Evolutionary dynamics of dispersal and local adaptation in multi-resource landscapes

Dispersal can enable access to resources in new locations. Consequently, traits that govern dispersal probability and dispersal distance may impact an individual's ability to acquire resources. However, spatial variation in the quality or quantity of resources may mediate potential adaptive benefits of novel dispersal traits. Ecological traits (i.e., those that determine how an individual processes resources) will also, by definition, affect how an individual interacts with the resource landscape. In a spatially heterogeneous environment, this creates potential for evolutionary feedbacks between dispersal-related traits and ecological traits. For example, dispersal may introduce individuals to novel resources, at which point there may be selection for local adaptation of ecological traits. Conversely, an individual's ability to utilize different resource types may determine how dispersal impacts fitness. Here, we develop an individual-based model to investigate co-evolution of dispersal and ecological traits in a landscape where multiple resources vary independently across space. We find that: (1) resource specialists can emerge and tend to evolve dispersal strategies suited to the structure of their preferred resource type and (2) generalists, when they emerge, tend to possess intermediate dispersal strategies. Lastly, we note that the effect of dispersal on the evolution of the ecological trait is weaker than vice versa and, as a result, appreciable heterogeneity in the abundance of resources across a landscape will likely obscure a signal of co-evolution.

Eco-evolutionary dynamics: An experimental demonstration in nature

Ecological change drives rapid evolution, which then should feed back to influence ecological change. A new study uses experiments with Timema stick insects to demonstrate such feedbacks in nature, revealing that they can be very rapid, strong, and stabilizing.

A stabilizing eco-evolutionary feedback loop in the wild

There is increasing evidence that evolutionary and ecological processes can operate on the same timescale1,2 (i.e., contemporary time). As such, evolution can be sufficiently rapid to affect ecological processes such as predation or competition. Thus, evolution can influence population, community, and ecosystem-level dynamics. Indeed, studies have now shown that evolutionary dynamics can alter community structure3,4,5,6 and ecosystem function.7,8,9,10 In turn, shifts in ecological dynamics driven by evolution might feed back to affect the evolutionary trajectory of individual species.11 This feedback loop, where evolutionary and ecological changes reciprocally affect one another, is a central tenet of eco-evolutionary dynamics.1,12 However, most work on such dynamics in natural populations has focused on one-way causal associations between ecology and evolution.13 Hence, direct empirical evidence for eco-evolutionary feedback is rare and limited to laboratory or mesocosm experiments.13,14,15,16 Here, we show in the wild that eco-evolutionary dynamics in a plant-feeding arthropod community involve a negative feedback loop. Specifically, adaptation in cryptic coloration in a stick-insect species mediates bird predation, with local maladaptation increasing predation. In turn, the abundance of arthropods is reduced by predation. Here, we experimentally manipulate arthropod abundance to show that these changes at the community level feed back to affect the stick-insect evolution. Specifically, low-arthropod abundance increases the strength of selection on crypsis, increasing local adaptation of stick insects in a negative feedback loop. Our results suggest that eco-evolutionary feedbacks are able to stabilize complex systems by preventing consistent directional change and therefore increasing resilience.

Recruitment and migration patterns reveal a key role for seed banks in the meta-population dynamics of an aquatic plant

Progressive habitat fragmentation threatens plant species with narrow habitat requirements. While local environmental conditions define population growth rates and recruitment success at the patch level, dispersal is critical for population viability at the landscape scale. Identifying the dynamics of plant meta-populations is often confounded by the uncertainty about soil-stored population compartments. We combined a landscape-scale assessment of an amphibious plant's population structure with measurements of dispersal complexity in time to track dispersal and putative shifts in functional connectivity. Using 13 microsatellite markers, we analyzed the genetic structure of extant Oenanthe aquatica populations and their soil seed banks in a kettle hole system to uncover hidden connectivity among populations in time and space. Considerable spatial genetic structure and isolation-by-distance suggest limited gene flow between sites. Spatial isolation and patch size showed minor effects on genetic diversity. Genetic similarity found among extant populations and their seed banks suggests increased local recruitment, despite some evidence of migration and recent colonization. Results indicate stepping-stone dispersal across adjacent populations. Among permanent and ephemeral demes the resulting meta-population demography could be determined by source-sink dynamics. Overall, these spatiotemporal connectivity patterns support mainland-island dynamics in our system, highlighting the importance of persistent seed banks as enduring sources of genetic diversity.

Genomic conservation of crop wild relatives: A case study of citrus

Conservation of crop wild relatives is critical for plant breeding and food security. The lack of clarity on the genetic factors that lead to endangered status or extinction create difficulties when attempting to develop concrete recommendations for conserving a citrus wild relative: the wild relatives of crops. Here, we evaluate the conservation of wild kumquat (Fortunella hindsii) using genomic, geographical, environmental, and phenotypic data, and forward simulations. Genome resequencing data from 73 accessions from the Fortunella genus were combined to investigate population structure, demography, inbreeding, introgression, and genetic load. Population structure was correlated with reproductive type (i.e., sexual and apomictic) and with a significant differentiation within the sexually reproducing population. The effective population size for one of the sexually reproducing subpopulations has recently declined to ~1,000, resulting in high levels of inbreeding. In particular, we found that 58% of the ecological niche overlapped between wild and cultivated populations and that there was extensive introgression into wild samples from cultivated populations. Interestingly, the introgression pattern and accumulation of genetic load may be influenced by the type of reproduction. In wild apomictic samples, the introgressed regions were primarily heterozygous, and genome-wide deleterious variants were hidden in the heterozygous state. In contrast, wild sexually reproducing samples carried a higher recessive deleterious burden. Furthermore, we also found that sexually reproducing samples were self-incompatible, which prevented the reduction of genetic diversity by selfing. Our population genomic analyses provide specific recommendations for distinct reproductive types and monitoring during conservation. This study highlights the genomic landscape of a wild relative of citrus and provides recommendations for the conservation of crop wild relatives.

Source-sink dynamics within a complex life history

Source-sink patch dynamics occur when movement from sources stabilizes sinks by compensating for low local vital rates. The mechanisms underlying source-sink dynamics may be complicated in species that undergo transitions between discrete life stages, particularly when stages have overlapping habitat requirements and similar movement abilities. In these species, for example, the demographic effects of movement by one stage may augment or offset the effects of movement by another stage. We used a stream salamander system to investigate patch dynamics within this form of complex life history. Specifically, we tested the hypothesis that the salamander Gyrinophilus porphyriticus experiences source-sink dynamics in riffles and pools, the dominant geomorphic patch types in headwater streams. We estimated stage-specific survival probabilities in riffles and pools and stage-specific movement probabilities between the two patch types using 8 years of capture-recapture data on 4491 individuals, including premetamorphic larvae and postmetamorphic adults. We then incorporated survival and movement probabilities into a stage-structured, two-patch model to determine the demographic interactions between riffles and pools. Monthly survival probabilities of both stages were higher in pools than in riffles. Larvae were more likely to move from riffles to pools, but adults were more likely to move from pools to riffles, despite experiencing much lower survival in riffles. In simulations, eliminating interpatch movements by both stages indicated that riffles are sinks that rely on immigration from pools for stability. Allowing only larvae to move stabilized both patch types, but allowing only adults to move destabilized pools due to the demographic cost of adult emigration. These results indicated that larval movement not only stabilizes riffles, but also offsets the destabilizing effects of maladaptive adult movement. Similar patch dynamics may emerge in any structured population in which movement and local vital rates differ by age, size, or stage. Addressing these forms of internal demographic structure in patch dynamics analyses will help to refine and advance general understanding of spatial ecology.

Dispersal plasticity driven by variation in fitness across species and environmental gradients

Dispersal plasticity, when organisms adjust their dispersal decisions depending on their environment, can play a major role in ecological and evolutionary dynamics, but how it relates to fitness remains scarcely explored. Theory predicts that high dispersal plasticity should evolve when environmental gradients have a strong impact on fitness. Using microcosms, we tested in five species of the genus Tetrahymena whether dispersal plasticity relates to differences in fitness sensitivity along three environmental gradients. Dispersal plasticity was species- and environment-dependent. As expected, dispersal plasticity was generally related to fitness sensitivity, with higher dispersal plasticity when fitness is more affected by environmental gradients. Individuals often preferentially disperse out of low fitness environments, but leaving environments that should yield high fitness was also commonly observed. We provide empirical support for a fundamental, but largely untested, assumption in dispersal theory: the extent of dispersal plasticity correlates with fitness sensitivity to the environment.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

Compiling forty years of guppy research to investigate the factors contributing to (non)parallel evolution

Examples of parallel evolution have been crucial for our understanding of adaptation via natural selection. However, strong parallelism is not always observed even in seemingly similar environments where natural selection is expected to favour similar phenotypes. Leveraging this variation in parallelism within well-researched study systems can provide insight into the factors that contribute to variation in adaptive responses. Here we analyse the results of 36 studies reporting 446 average trait values in Trinidadian guppies, Poecilia reticulata, from different predation regimes. We examine how the extent of predator-driven phenotypic parallelism is influenced by six factors: sex, trait type, rearing environment, ecological complexity, evolutionary history, and time since colonization. Analyses show that parallel evolution in guppies is highly variable and weak on average, with only 24.7% of the variation among populations being explained by predation regime. Levels of parallelism appeared to be especially weak for colour traits, and parallelism decreased with increasing complexity of evolutionary history (i.e., when estimates of parallelism from populations within a single drainage were compared to estimates of parallelism from populations pooled between two major drainages). Suggestive - but not significant - trends that warrant further research include interactions between the sexes and different trait categories. Quantifying and accounting for these and other sources of variation among evolutionary 'replicates' can be leveraged to better understand the extent to which seemingly similar environments drive parallel and nonparallel aspects of phenotypic divergence.

Tracking of SARS-CoV-2 Alpha variant (B.1.1.7) in Palestine

As surges of the COVID-19 pandemic continue globally, including in Palestine, several new SARS-CoV-2 variants have been introduced. This expansion has impacted transmission, disease severity, virulence, diagnosis, therapy, and natural and vaccine-induced immunity. Here, 183 whole genome sequences (WGS) were analyzed, of which 129 were from Palestinian cases, 62 of which were collected in 11 Palestinian districts between October 2020 and April 2021 and sequenced completely. A dramatic shift from the wild type to the Alpha variant (B 1.1.7) was observed within a short period of time. Cluster mapping revealed statistically significant clades in two main Palestinian cities, Al-Khalil (Monte Carlo hypothesis test-Poisson model, P = 0.00000000012) and Nablus (Monte Carlo hypothesis test-Poisson model, P = 0.014 and 0.015). The phylogenetic tree showed three main clusters of SARS-CoV-2 with high bootstrap values (>90). However, population genetics analysis showed a genetically homogenous population supported by low Wright's F-statistic values (Fst <0.25), high gene flow (Nm > 3), and statistically insignificant Tajima's D values (Tajima's test, neutrality model prediction, P = 0.02). The Alpha variant, rapidly replaced the wild type, causing a major surge that peaked in April 2021, with an increased COVID-19 mortality rate, especially, in the Al-Khalil and Nablus districts. The source of introduction remains uncertain, despite the minimal genetic variation. The study substantiates the use of WGS for SARS-CoV-2 surveillance as an early warning system to track down new variants requiring effective control.

Target enrichment and extensive population sampling help untangle the recent, rapid radiation of Oenothera sect. Calylophus

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking non-coding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic and morphological diversity in the group. We found that the combination of exons and flanking non-coding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation.

Genetic load: genomic estimates and applications in non-model animals

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Gene flow influences the genomic architecture of local adaptation in six riverine fish species

Understanding how gene flow influences adaptive divergence is important for predicting adaptive responses. Theoretical studies suggest that when gene flow is high, clustering of adaptive genes in fewer genomic regions would protect adaptive alleles from recombination and thus be selected for, but few studies have tested it with empirical data. Here, we used restriction site-associated sequencing to generate genomic data for six fish species with contrasting life histories from six reaches of the Upper Mississippi River System, USA. We used four differentiation-based outlier tests and three genotype-environment association analyses to define neutral single nucleotide polymorphisms (SNPs) and outlier SNPs that were putatively under selection. We then examined the distribution of outlier SNPs along the genome and investigated whether these SNPs were found in genomic islands of differentiation and inversions. We found that gene flow varied among species, and outlier SNPs were clustered more tightly in species with higher gene flow. The two species with the highest overall F_ST (0.0303-0.0720) and therefore lowest gene flow showed little evidence of clusters of outlier SNPs, with outlier SNPs in these species spreading uniformly across the genome. In contrast, nearly all outlier SNPs in the species with the lowest F_ST (0.0003) were found in a single large putative inversion. Two other species with intermediate gene flow (F_ST  ~ 0.0025-0.0050) also showed clustered genomic architectures, with most islands of differentiation clustered on a few chromosomes. Our results provide important empirical evidence to support the hypothesis that increasingly clustered architecture of local adaptation is associated with high gene flow.

Evolutionary and phylogenetic analyses of the barcoding region suggest geographical relationships among Blastocystis sp., ST3 in humans

Blastocystis sp., has 21 distinct subtypes of which ST3 thought to be the most prevalent subtype. This study aims to analyze the global variations of ST3. In total, 496 sequences with more than 400 nucleotides from Asia, Europe, Africa, and America were included in this study. Results show that allele 34 was the most prevalent allele in all continents. The lowest and highest allele diversity were observed in Europe and Africa, respectively. The nucleotide diversity ranged from 0.0077 in Europe to 0.02 in Africa, and haplotype diversity ranged from 0.461 in America to 0.6 in Africa. The haplotype network and Bayesian structure showed at least two major clusters including Asia and Europe-Africa-America. Tajima's D values for all continents were negative and statistically significant, indicating an excess of rare nucleotide variants. Similarly, the Fu's FS test showed negative values for all regions, indicating an excess of rare haplotypes. Pairwise F_ST exhibited a high genetic differentiation between Asia and other continents. Mismatch analysis for all populations showed a unimodal distribution. Our findings indicate that there are two probable major clusters of Blastocystis sp. ST3, a cluster which is shared between Europe, Africa, and America, and a cluster which is restricted to Asia.

Parallel shape divergence between ecotypes of the limpet Nacella concinna along the Antarctic Peninsula: a new model species for parallel evolution?

Parallel phenotypic divergence is the independent differentiation between phenotypes of the same lineage or species occupying ecologically similar environments in different populations. We tested in the Antarctic limpet Nacella concinna the extent of parallel morphological divergence in littoral and sublittoral ecotypes throughout its distribution range. These ecotypes differ in morphological, behavioural and physiological characteristics. We studied the lateral and dorsal outlines of shells and the genetic variation of the mitochondrial gene Cytochrome Oxidase subunit I from both ecotypes in 17 sample sites along more than 2,000 km. The genetic data indicate that both ecotypes belong to a single evolutionary lineage. The magnitude and direction of phenotypic variation differ between ecotypes across sample sites; completely parallel ecotype-pairs (i.e., they diverge in the same magnitude and in the same direction) were detected in 84.85% of lateral and 65.15% in dorsal view comparisons. Besides, specific traits (relative shell height, position of shell apex, and elliptical/pear-shape outline variation) showed high parallelism. We observed weak morphological covariation between the two shape shell views, indicating that distinct evolutionary forces and environmental pressures could be acting on this limpet shell shape. Our results demonstrate there is a strong parallel morphological divergence pattern in N. concinna along its distribution, making this Antarctic species a suitable model for the study of different evolutionary forces shaping the shell evolution of this limpet.

The maintenance of standing genetic variation: Gene flow vs. selective neutrality in Atlantic stickleback fish

Adaptation to derived habitats often occurs from standing genetic variation. The maintenance within ancestral populations of genetic variants favourable in derived habitats is commonly ascribed to long‐term antagonism between purifying selection and gene flow resulting from hybridization across habitats. A largely unexplored alternative idea based on quantitative genetic models of polygenic adaptation is that variants favoured in derived habitats are neutral in ancestral populations when their frequency is relatively low. To explore the latter, we first identify genetic variants important to the adaptation of threespine stickleback fish (Gasterosteus aculeatus) to a rare derived habitat—nutrient‐depleted acidic lakes—based on whole‐genome sequence data. Sequencing marine stickleback from six locations across the Atlantic Ocean then allows us to infer that the frequency of these derived variants in the ancestral habitat is unrelated to the likely opportunity for gene flow of these variants from acidic‐adapted populations. This result is consistent with the selective neutrality of derived variants within the ancestor. Our study thus supports an underappreciated explanation for the maintenance of standing genetic variation, and calls for a better understanding of the fitness consequences of adaptive variation across habitats and genomic backgrounds.

Biodiversity, resilience and the stability of evolutionary systems

Various macro-evolutionary phenomena, such as long-term stability punctuated by bursts of evolution, are difficult to explain via the micro-evolutionary process of weak selection acting steadily on individual mutations. In contrast, bursts of change are expected if evolutionary systems are complex and balanced, with occasional disruption of balance. Such disruption represents the collapse of resilience, akin to the snapping of an elastic band. It can be driven by external factors, or by self-propagating feedback loops internal to a system. Thus, evolutionary resilience could help explain how evolution generates broader patterns of biodiversity. We outline evidence and tests for this hypothesis, which emphasizes the processes balancing evolution, as urged fifty years ago in ecological genetics and via modern results in a range of systems.

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.

Rapid microgeographic evolution in response to climate change

Environmental change is predicted to accelerate into the future and will exert strong selection pressure on biota. Although many species may be fated to extinction, others may survive through their capacity to evolve rapidly at highly localized (i.e., microgeographic) scales. Yet, even as new examples have been discovered, the limits to such evolutionary responses have not often been evaluated. One of the first examples of microgeographic variation involved pond populations of wood frogs (Rana sylvatica). Although separated by just tens to hundreds of meters, these populations exhibited countergradient variation in intrinsic embryonic development rates when reared in a common garden. We repeated this experiment 17 years (approximately six to nine generations) later and found that microgeographic variation persists in contemporary populations. Furthermore, we found that contemporary embryos have evolved to develop 14-19% faster than those in 2001. Structural equation models indicate that the predominant cause for this response is likely due to changes in climate over the intervening 17 years. Despite potential for rapid and fine-scale evolution, demographic declines in populations experiencing the greatest changes in climate and habitat imply a limit to the species' ability to mitigate extreme environmental change.

Are immigrants outbred and unrelated? Testing standard assumptions in a wild metapopulation

Immigration into small recipient populations is expected to alleviate inbreeding and increase genetic variation, and hence facilitate population persistence through genetic and/or evolutionary rescue. Such expectations depend on three standard assumptions: that immigrants are outbred, unrelated to existing natives at arrival, and unrelated to each other. These assumptions are rarely explicitly verified, including in key field systems in evolutionary ecology. Yet, they could be violated due to non-random or repeated immigration from adjacent small populations. We combined molecular genetic marker data for 150-160 microsatellite loci with comprehensive pedigree data to test the three assumptions for a song sparrow (Melospiza melodia) population that is a model system for quantifying effects of inbreeding and immigration in the wild. Immigrants were less homozygous than existing natives on average, with mean homozygosity that closely resembled outbred natives. Immigrants can therefore be considered outbred on the focal population scale. Comparisons of homozygosity of real or hypothetical offspring of immigrant-native, native-native and immigrant-immigrant pairings implied that immigrants were typically unrelated to existing natives and to each other. Indeed, immigrants' offspring would be even less homozygous than outbred individuals on the focal population scale. The three standard assumptions of population genetic and evolutionary theory were consequently largely validated. Yet, our analyses revealed some deviations that should be accounted for in future analyses of heterosis and inbreeding depression, implying that the three assumptions should be verified in other systems to probe patterns of non-random or repeated dispersal and facilitate precise and unbiased estimation of key evolutionary parameters.

Adaptation in a heterogeneous environment I: persistence versus extinction

Understanding how a diversity of plants in agroecosystems affects the adaptation of pathogens is a key issue in agroecology. We analyze PDE systems describing the dynamics of adaptation of two phenotypically structured populations, under the effects of mutation, selection and migration in a two-patch environment, each patch being associated with a different phenotypic optimum. We consider two types of growth functions that depend on the n-dimensional phenotypic trait: either local and linear or nonlocal nonlinear. In both cases, we obtain existence and uniqueness results as well as a characterization of the large-time behaviour of the solution (persistence or extinction) based on the sign of a principal eigenvalue. We show that migration between the two environments decreases the chances of persistence, with in some cases a 'lethal migration threshold' above which persistence is not possible. Comparison with stochastic individual-based simulations shows that the PDE approach accurately captures this threshold. Our results illustrate the importance of cultivar mixtures for disease prevention and control.

Evolutionary mismatch along salinity gradients in a Neotropical water strider

The evolution of local adaptation is crucial for the in situ persistence of populations in changing environments. However, selection along broad environmental gradients could render local adaptation difficult, and might even result in maladaptation. We address this issue by quantifying fitness trade-offs (via common garden experiments) along a salinity gradient in two populations of the Neotropical water strider Telmatometra withei-a species found in both fresh (FW) and brackish (BW) water environments across Panama. We found evidence for local adaptation in the FW population in its home FW environment. However, the BW population showed only partial adaptation to the BW environment, with a high magnitude of maladaptation along naturally occurring salinity gradients. Indeed, its overall fitness was ~60% lower than that of the ancestral FW population in its home environment, highlighting the role of phenotypic plasticity, rather than local adaptation, in high salinity environments. This suggests that populations seemingly persisting in high salinity environments might in fact be maladapted, following drastic changes in salinity. Thus, variable selection imposed by salinization could result in evolutionary mismatch, where the fitness of a population is displaced from its optimal environment. Understanding the fitness consequences of persisting in fluctuating salinity environments is crucial to predict the persistence of populations facing increasing salinization. It will also help develop evolutionarily informed management strategies in the context of global change.

Standing genetic variation in laboratory populations of insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis (Diptera: Psychodidae: Phlebotominae) for the evolution of resistance

Insecticides can exert strong selection on insect pest species, including those that vector diseases, and have led to rapid evolution of resistance. Despite such rapid evolution, relatively little is known about standing genetic variation for resistance in insecticide-susceptible populations of many species. To help fill this knowledge gap, we generated genotyping-by-sequencing data from insecticide-susceptible Phlebotomus papatasi and Lutzomyia longipalpis sand flies that survived or died from a sub-diagnostic exposure to either permethrin or malathion using a modified version of the Centers for Disease Control and Prevention bottle bioassay. Multi-locus genome-wide association mapping methods were used to quantify standing genetic variation for insecticide resistance in these populations and to identify specific alleles associated with insecticide survival. For each insecticide treatment, we estimated the proportion of the variation in survival explained by the genetic data (i.e., "chip" heritability) and the number and contribution of individual loci with measurable effects. For all treatments, survival to an insecticide exposure was heritable with a polygenic architecture. Both P. papatasi and L. longipalpis had alleles for survival that resided within many genes throughout their genomes. The implications for resistance conferred by many alleles, as well as inferences made about the utility of laboratory insecticide resistance association studies compared to field observations, are discussed.

Suppressing evolution in genetically engineered systems through repeated supplementation

Genetically engineered organisms are prone to evolve in response to the engineering. This evolution is often undesirable and can negatively affect the purpose of the engineering. Methods that maintain the stability of engineered genomes are therefore critical to the successful design and use of genetically engineered organisms. One potential method to limit unwanted evolution is by taking advantage of the ability of gene flow to counter local adaption, a process of supplementation. Here, we investigate the feasibility of supplementation as a mechanism to offset the evolutionary degradation of a transgene in three model systems: a bioreactor, a gene drive, and a transmissible vaccine. In each model, continual introduction from a stock is used to balance mutation and selection against the transgene. Each system has its unique features. The bioreactor system is especially tractable and has a simple answer: The level of supplementation required to maintain the transgene at a frequency p ^ is approximatelyp ^ s , where s is the selective disadvantage of the transgene. Supplementation is also feasible in the transmissible vaccine case but is probably not practical to prevent the evolution of resistance against a gene drive. We note, however, that the continual replacement of even a small fraction of a large population can be challenging, limiting the usefulness of supplementation as a means of controlling unwanted evolution.

Immigration counter-acts local micro-evolution of a major fitness component: Migration-selection balance in free-living song sparrows

Ongoing adaptive evolution, and resulting “evolutionary rescue” of declining populations, requires additive genetic variation in fitness. Such variation can be increased by gene flow resulting from immigration, potentially facilitating evolution. But, gene flow could in fact constrain rather than facilitate local adaptive evolution if immigrants have low additive genetic values for local fitness. Local migration‐selection balance and micro‐evolutionary stasis could then result. However, key quantitative genetic effects of natural immigration, comprising the degrees to which gene flow increases the total local additive genetic variance yet counteracts local adaptive evolutionary change, have not been explicitly quantified in wild populations. Key implications of gene flow for population and evolutionary dynamics consequently remain unclear. Our quantitative genetic analyses of long‐term data from free‐living song sparrows (Melospiza melodia) show that mean breeding value for local juvenile survival to adulthood, a major component of fitness, increased across cohorts more than expected solely due to drift. Such micro‐evolutionary change should be expected given nonzero additive genetic variance and consistent directional selection. However, this evolutionary increase was counteracted by negative additive genetic effects of recent immigrants, which increased total additive genetic variance but prevented a net directional evolutionary increase in total additive genetic value. These analyses imply an approximate quantitative genetic migration‐selection balance in a major fitness component, and hence demonstrate a key mechanism by which substantial additive genetic variation can be maintained yet decoupled from local adaptive evolutionary change.

Limited divergent adaptation despite a substantial environmental cline in wild pea

Isolation by environment (IBE) is a widespread phenomenon in nature. It is commonly expected that the degree of difference among environments is proportional to the level of divergence between populations in their respective environments. It is therefore assumed that a species' genetic diversity displays a pattern of IBE in the presence of a strong environmental cline if gene flow does not mitigate isolation. We tested this common assumption by analysing the genetic diversity and demographic history of Pisum fulvum, which inhabits contrasting habitats in the southern Levant and is expected to display only minor migration rates between populations, making it an ideal test case. Ecogeographical and subpopulation structure were analysed and compared. The correlation of genetic with environmental distances was calculated to test the effect of isolation by distance and IBE and detect the main drivers of these effects. Historical effective population size was estimated using stairway plot. Limited overlap of ecogeographical and genetic clustering was observed, and correlation between genetic and environmental distances was statistically significant but small. We detected a sharp decline of effective population size during the last glacial period. The low degree of IBE may be the result of genetic drift due to a past bottleneck. Our findings contradict the expectation that strong environmental clines cause IBE in the absence of extensive gene flow.

Gene flow relates to evolutionary divergence among populations at the range margin

Background

Morphological differentiation between populations resulting from local adaptations to environmental conditions is likely to be more pronounced in populations with increasing genetic isolation. In a previous study a positive clinal variation in body size was observed in isolated Roesel’s bush-cricket, Metrioptera roeselii, populations, but were absent from populations within a continuous distribution at the same latitudinal range. This observational study inferred that there was a phenotypic effect of gene flow on climate-induced selection in this species.

Methods

To disentangle genetic versus environmental drivers of population differences in morphology, we measured the size of four different body traits in wild-caught individuals from the two most distinct latitudinally-matched pairs of populations occurring at about 60°N latitude in northern Europe, characterised by either restricted or continuous gene flow, and corresponding individuals raised under laboratory conditions.

Results

Individuals that originated from the genetically isolated populations were always bigger (femur, pronotum and genital appendages) when compared to individuals from latitudinally-matched areas characterised by continuous gene flow between populations. The magnitude of this effect was similar for wild-caught and laboratory-reared individuals. We found that previously observed size cline variation in both male and female crickets was likely to be the result of local genetic adaptation rather than phenotypic plasticity.

Conclusions

This strongly suggests that restricted gene flow is of major importance for frequencies of alleles that participate in climate-induced selection acting to favour larger phenotypes in isolated populations towards colder latitudes.

Mixed support for gene flow as a constraint to local adaptation and contributor to the limited geographic range of an endemic salamander

Understanding mechanisms that underlie species range limits is at the core of evolutionary ecology. Asymmetric gene flow between larger core populations and smaller edge populations can swamp local adaptation at the range edge and inhibit further range expansion. However, empirical tests of this theory are exceedingly rare. We tested the hypothesis that asymmetric gene flow can constrain local adaptation and thereby species' range limits in an endemic US salamander (Ambystoma barbouri) by determining if gene flow is asymmetric between the core and peripheries of the species' geographic distribution and testing whether local adaptation is swamped at range edges with a reciprocal transplant experiment. Using putatively neutral loci from populations across three core-to-edge transects that covered nearly the entire species' geographic range, we found evidence for asymmetric, core-to-edge gene flow along western and northern transects, but not along a southern transect. Subsequently, the reciprocal transplant experiment suggested that northern and western edge populations are locally adapted despite experiencing asymmetric gene flow, yet have lower fitness in their respective home regions than those of centre population. Conversely, southern populations exhibit low deme quality, experiencing high mortality regardless of where they were reared, probably due to harsher edge habitat conditions. Consequently, we provide rare species-wide evidence that local adaptation can occur despite asymmetric gene flow, though migration from the core may prohibit range expansion by reducing fitness in edge populations. Further, our multitransect study shows that multiple, nonmutually exclusive mechanisms can lead to range limits within a single species.

Exploring context dependency in eco-evolutionary patterns with the stick insect Timema cristinae

Rapid evolution can influence the ecology of populations, communities, and ecosystems, but the importance of evolution for ecological dynamics remains unclear, largely because the contexts in which evolution is powerful are poorly resolved. Here, we carry out a large observational study to test hypotheses about context dependency of eco-evolutionary patterns previously identified on the stick insect Timema cristinae. Experiments and observations conducted in 2011 and 2012 documented predator-mediated negative effects of camouflage maladaptation (i.e., evolutionary dynamics) on: (a) T. cristinae abundance and, (b) species richness and abundance of other arthropods. Here we show that camouflage maladaptation does not correlate with T. cristinae abundance and, instead, is associated with increased abundance and species richness of cohabitating arthropods. We furthermore find that plants with high levels of Timema maladaptation tend to have higher foliar nitrogen, that is, higher nutritional value, and more positive mass-abundance slopes in the coexisting arthropod communities. We propose explanations for the observed contrasting results, such as negative density- and frequency-dependent selection, feedbacks between herbivore abundance and plant nutritional quality, and common effects of predation pressure on selection and prey abundance. Our results demonstrate the utility of observational studies to assess the context dependency of eco-evolutionary dynamics patterns and provide testable hypotheses for future work.

Variation in breeding practices and geographic isolation drive subpopulation differentiation, contributing to the loss of genetic diversity within dog breed lineages

Background

Discrete breed ideals are not restricted to delimiting dog breeds from another, but also are key drivers of subpopulation differentiation. As genetic differentiation due to population fragmentation results in increased rates of inbreeding and loss of genetic diversity, detecting and alleviating the reasons of population fragmentation can provide effective tools for the maintenance of healthy dog breeds.

Results

Using a genome-wide SNP array, we detected genetic differentiation to subpopulations in six breeds, Belgian Shepherd, English Greyhound, Finnish Lapphund, Italian Greyhound, Labrador Retriever and Shetland Sheepdog, either due to geographical isolation or as a result of differential breeding strategies. The subpopulation differentiation was strongest in show dog lineages.

Conclusions

Besides geographical differentiation caused by founder effect and lack of gene flow, selection on champion looks or restricted pedigrees is a strong driver of population fragmentation. Artificial barriers for gene flow between the different subpopulations should be recognized, their necessity evaluated critically and perhaps abolished in order to maintain genetic diversity within a breed. Subpopulation differentiation might also result in false positive signals in genome-wide association studies of different traits.

Lay summary

Purebred dogs are, by definition, reproductively isolated from other breeds. However, similar isolation can also occur within a breed due to conflicting breeder ideals and geographic distances between the dog populations. We show here that both of these examples can contribute to breed division, with subsequent loss of genetic variation in the resulting breed lineages. Breeders should avoid creating unnecessary boundaries between breed lineages and facilitate the exchange of dogs between countries.

Springtail phylogeography highlights biosecurity risks of repeated invasions and intraregional transfers among remote islands

Human-mediated transport of species outside their natural range is a rapidly growing threat to biodiversity, particularly for island ecosystems that have evolved in isolation. The genetic structure underpinning island populations will largely determine their response to increased transport and thus help to inform biosecurity management. However, this information is severely lacking for some groups, such as the soil fauna. We therefore analysed the phylogeographic structure of an indigenous and an invasive springtail species (Collembola: Poduromorpha), each distributed across multiple remote sub-Antarctic islands, where human activity is currently intensifying. For both species, we generated a genome-wide SNP data set and additionally analysed all available COI barcodes. Genetic differentiation in the indigenous springtail Tullbergia bisetosa is substantial among (and, to a lesser degree, within) islands, reflecting low dispersal and historic population fragmentation, while COI patterns reveal ancestral signatures of postglacial recolonization. This pronounced geographic structure demonstrates the key role of allopatric divergence in shaping the region's diversity and highlights the vulnerability of indigenous populations to genetic homogenization via human transport. For the invasive species Hypogastrura viatica, nuclear genetic structure is much less apparent, particularly for islands linked by regular shipping, while diverged COI haplotypes indicate multiple independent introductions to each island. Thus, human transport has likely facilitated this species' persistence since its initial colonization, through the ongoing introduction and inter-island spread of genetic variation. These findings highlight the different evolutionary consequences of human transport for indigenous and invasive soil species. Crucially, both outcomes demonstrate the need for improved intraregional biosecurity among remote island systems, where the policy focus to date has been on external introductions.

The role of common ancestry and gene flow in the evolution of human-directed play behaviour in dogs

Among-population variance of phenotypic traits is of high relevance for understanding evolutionary mechanisms that operate in relatively short timescales, but various sources of nonindependence, such as common ancestry and gene flow, can hamper the interpretations. In this comparative analysis of 138 dog breeds, we demonstrate how such confounders can independently shape the evolution of a behavioural trait (human-directed play behaviour from the Dog Mentality Assessment project). We combined information on genetic relatedness and haplotype sharing to reflect common ancestry and gene flow, respectively, and entered these into a phylogenetic mixed model to partition the among-breed variance of human-directed play behaviour while also accounting for within-breed variance. We found that 75% of the among-breed variance was explained by overall genetic relatedness among breeds, whereas 15% could be attributed to haplotype sharing that arises from gene flow. Therefore, most of the differences in human-directed play behaviour among breeds have likely been caused by constraints of common ancestry as a likely consequence of past selection regimes. On the other hand, gene flow caused by crosses among breeds has played a minor, but not negligible role. Our study serves as an example of an analytical approach that can be applied to comparative situations where the effects of shared origin and gene flow require quantification and appropriate statistical control in a within-species/among-population framework. Altogether, our results suggest that the evolutionary history of dog breeds has left remarkable signatures on the among-breed variation of a behavioural phenotype.

Understanding Maladaptation by Uniting Ecological and Evolutionary Perspectives

Evolutionary biologists have long trained their sights on adaptation, focusing on the power of natural selection to produce relative fitness advantages while often ignoring changes in absolute fitness. Ecologists generally have taken a different tack, focusing on changes in abundance and ranges that reflect absolute fitness while often ignoring relative fitness. Uniting these perspectives, we articulate various causes of relative and absolute maladaptation and review numerous examples of their occurrence. This review indicates that maladaptation is reasonably common from both perspectives, yet often in contrasting ways. That is, maladaptation can appear strong from a relative fitness perspective, yet populations can be growing in abundance. Conversely, resident individuals can appear locally adapted (relative to nonresident individuals) yet be declining in abundance. Understanding and interpreting these disconnects between relative and absolute maladaptation, as well as the cases of agreement, is increasingly critical in the face of accelerating human-mediated environmental change. We therefore present a framework for studying maladaptation, focusing in particular on the relationship between absolute and relative fitness, thereby drawing together evolutionary and ecological perspectives. The unification of these ecological and evolutionary perspectives has the potential to bring together previously disjunct research areas while addressing key conceptual issues and specific practical problems.

Quantifying maladaptation during the evolution of sexual dimorphism

Females and males have distinct trait optima, resulting in selection for sexual dimorphism. However, most traits have strong cross-sex genetic correlations, which constrain evolutionary divergence between the sexes and lead to protracted periods of maladaptation during the evolution of sexual dimorphism. While such constraints are thought to be costly in terms of individual and population fitness, it remains unclear how severe such costs are likely to be. Building upon classical models for the 'cost of selection' in changing environments (sensu Haldane), we derived a theoretical expression for the analogous cost of evolving sexual dimorphism; this cost is a simple function of genetic (co)variances of female and male traits and sex differences in trait optima. We then conducted a comprehensive literature search, compiled quantitative genetic data from a diverse set of traits and populations, and used them to quantify costs of sexual dimorphism in the light of our model. For roughly 90% of traits, costs of sexual dimorphism appear to be modest, and comparable to the costs of fixing one or a few beneficial substitutions. For the remaining traits (approx. 10%), sexual dimorphism appears to carry a substantial cost-potentially orders of magnitude greater than costs of selection during adaptation to environmental changes.

Implications of existing local (mal)adaptations for ecological forecasting under environmental change

Standing genetic variation represents a genetic load on population fitness but can also support a rapid response to short-term environmental change, and the greatest potential source of such standing genetic variation typically exists among locally adapted populations living along an environmental gradient. Here, we develop a spatially explicit simulation model to quantify the contribution of existing genetic variation arising from migration-mutation-selection-drift balance to time to extinction under environmental change. Simulations reveal that local adaptation across a species range associated with an underlying environmental gradient could extend time to extinction by nearly threefold irrespective of the rate of environmental change. The potential for preadapted alleles to increase the rate of adaptation changes the relative importance of established extinction risk factors; in particular, it reduced the importance of the breadth of environmental tolerance and it increased the relative importance of fecundity. Although migration of preadapted alleles generally increased persistence time, it decreased it at rates of environmental change close to the critical rate of change by creating a population bottleneck, which ultimately limited the rate at which de novo mutations could arise. An analysis of the extinction dynamics further revealed that one consequence of gene flow is the potential to maximize population growth rate in at least part of the species range, which is likely to have consequences for forecasting the consequences of ecological interactions. Our study shows that predictions of persistence time change fundamentally when existing local adaptations are explicitly taken into account, underscoring the need to preserve and manage genetic diversity.

Phenotype-environment mismatch in metapopulations-Implications for the maintenance of maladaptation at the regional scale

Maladaptation is widespread in natural populations. However, maladaptation has most often been associated with absolute population decline in local habitats rather than on a spectrum of relative fitness variation that can assist natural populations in their persistence at larger regional scales. We report results from a field experiment that tested for relative maladaptation between-pond habitats with spatial heterogeneity and (a)symmetric selection in pH. In the experiment, we quantified relative maladaptation in a copepod metapopulation as a mismatch between the mean population phenotype and the optimal trait value that would maximize mean population fitness under either stable or fluctuating pH environmental conditions. To complement the field experiment, we constructed a metapopulation model that addressed both relative (distance from the optimum) and absolute (negative population growth) maladaptation, with the aim of forecasting maladaptation to pH at the regional scale in relation to spatial structure (environmental heterogeneity and connectivity) and temporal environmental fluctuations. The results from our experiment indicated that maladaptation to pH at the regional scale depended on the asymmetry of the fitness surface at the local level. The results from our metapopulation model revealed how dispersal and (a)symmetric selection can operate on the fitness surface to maintain maladaptive phenotype-environment mismatch at local and regional scales in a metapopulation. Environmental stochasticity resulted in the maintenance of maladaptation that was robust to dispersal, but also revealed an interaction between the asymmetry in selection and environmental correlation. Our findings emphasize the importance of maladaptation for planning conservation strategies that can support adaptive potential in fragmented and changing landscapes.

Causes of maladaptation

Evolutionary biologists tend to approach the study of the natural world within a framework of adaptation, inspired perhaps by the power of natural selection to produce fitness advantages that drive population persistence and biological diversity. In contrast, evolution has rarely been studied through the lens of adaptation's complement, maladaptation. This contrast is surprising because maladaptation is a prevalent feature of evolution: population trait values are rarely distributed optimally; local populations often have lower fitness than imported ones; populations decline; and local and global extinctions are common. Yet we lack a general framework for understanding maladaptation; for instance in terms of distribution, severity, and dynamics. Similar uncertainties apply to the causes of maladaptation. We suggest that incorporating maladaptation-based perspectives into evolutionary biology would facilitate better understanding of the natural world. Approaches within a maladaptation framework might be especially profitable in applied evolution contexts - where reductions in fitness are common. Toward advancing a more balanced study of evolution, here we present a conceptual framework describing causes of maladaptation. As the introductory article for a Special Feature on maladaptation, we also summarize the studies in this Issue, highlighting the causes of maladaptation in each study. We hope that our framework and the papers in this Special Issue will help catalyze the study of maladaptation in applied evolution, supporting greater understanding of evolutionary dynamics in our rapidly changing world.

Transient local adaptation and source-sink dynamics in experimental populations experiencing spatially heterogeneous environments

Local adaptation is determined by the strength of selection and the level of gene flow within heterogeneous landscapes. The presence of benign habitat can act as an evolutionary stepping stone for local adaptation to challenging environments by providing the necessary genetic variation. At the same time, migration load from benign habitats will hinder adaptation. In a community context, interspecific competition is expected to select against maladapted migrants, hence reducing migration load and facilitating adaptation. As the interplay between competition and spatial heterogeneity on the joint ecological and evolutionary dynamics of populations is poorly understood, we performed an evolutionary experiment using the herbivore spider mite Tetranychus urticae as a model. We studied the species's demography and local adaptation in a challenging environment that consisted of an initial sink (pepper plants) and/or a more benign environment (cucumber plants). Half of the experimental populations were exposed to a competitor, the congeneric T. ludeni. We show that while spider mites only adapted to the challenging pepper environment when it was spatially interspersed with benign cucumber habitat, this adaptation was only temporary and disappeared when the populations in the benign cucumber environment were expanding and spilling-over to the challenging pepper environment. Although the focal species outcompeted the competitor after about two months, a negative effect of competition on the focal species's performance persisted in the benign environment. Adaptation to challenging habitat in heterogeneous landscapes thus highly depends on demography and source-sink dynamics, but also on competitive interactions with other species, even if they are only present for a short time span.

Population structure promotes the evolution of costly sex in artificial gene networks

We build on previous observations that Hill-Robertson interference generates an advantage of sex that, in structured populations, can be large enough to explain the evolutionary maintenance of costly sex. We employed a gene network model that explicitly incorporates interactions between genes. Mutations in the gene networks have variable effects that depend on the genetic background in which they appear. Consequently, our simulations include two costs of sex-recombination and migration loads-that were missing from previous studies of the evolution of costly sex. Our results suggest a critical role for population structure that lies in its ability to align the long- and short-term advantages of sex. We show that the addition of population structure favored the evolution of sex by disproportionately decreasing the equilibrium mean fitness of asexual populations, primarily by increasing the strength of Muller's Ratchet. Population structure also increased the ability of the short-term advantage of sex to counter the primary limit to the evolution of sex in the gene network model-recombination load. On the other hand, highly structured populations experienced migration load in the form of Dobzhansky-Muller incompatibilities, decreasing the effective rate of migration between demes and, consequently, accelerating the accumulation of drift load in the sexual populations.

Evidence of local adaptation in a waterfall-climbing Hawaiian goby fish derived from coupled biophysical modeling of larval dispersal and post-settlement selection

Background

Local adaptation of marine and diadromous species is thought to be a product of larval dispersal, settlement mortality, and differential reproductive success, particularly in heterogeneous post-settlement habitats. We evaluated this premise with an oceanographic passive larval dispersal model coupled with individual-based models of post-settlement selection and reproduction to infer conditions that underlie local adaptation in Sicyopterus stimpsoni, an amphidromous Hawaiian goby known for its ability to climb waterfalls.

Results

Our model results demonstrated that larval dispersal is spatio-temporally asymmetric, with more larvae dispersed from the southeast (the Big Island) to northwest (Kaua‘i) along the archipelago, reflecting prevailing conditions such as El Niño/La Niña oscillations. Yet connectivity is nonetheless sufficient to result in homogenous populations across the archipelago. We also found, however, that ontogenetic shifts in habitat can give rise to adaptive morphological divergence when the strength of predation-driven post-settlement selection crosses a critical threshold. Notably, our simulations showed that larval dispersal is not the only factor determining the likelihood of morphological divergence. We found adaptive potential and evolutionary trajectories of S. stimpsoni were greater on islands with stronger environmental gradients and greater variance in larval cohort morphology due to fluctuating immigration.

Conclusions

Contrary to expectation, these findings indicate that immigration can act in concert with selection to favor local adaptation and divergence in species with marine larval dispersal. Further development of model simulations, parameterized to reflect additional empirical estimates of abiotic and biotic factors, will help advance our understanding of the proximate and ultimate mechanisms driving adaptive evolution, population resilience, and speciation in marine-associated species.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1413-4) contains supplementary material, which is available to authorized users.

Genomic signatures of introgression between commercial and native bumblebees, Bombus terrestris, in western Iberian Peninsula-Implications for conservation and trade regulation

Human-mediated introductions of species may have profound impacts on native ecosystems. One potential impact with largely unforeseen consequences is the potential admixture of introduced with autochthonous species through hybridization. Throughout the world, bumblebees have been deliberately introduced for crop pollination with known negative impacts on native pollinators. Given the likely allochthonous origin of commercial bumblebees used in Portugal (subspecies Bombus terrestris terrestris and B. t. dalmatinus), our aim was to assess their putative introgression with the native Iberian subspecies B. terrestris lusitanicus. We analysed one mitochondrial gene, cytochrome c oxidase subunit I (COX1) and genomic data involving thousands of genome-wide restriction-site-associated DNA markers (RAD-seq). In the mitochondrial COX1 analyses, we detected one relatively common haplotype in commercial bumblebees, also present in wild samples collected nearby the greenhouses where the commercial hives are used. In the RAD-seq analysis, we found a clear genetic differentiation between native and commercial lineages. Furthermore, we detected candidate hybrids in the wild, as well as putatively escaped commercial bumblebees, some of which being potentially fertile males. Although we cannot assess directly the fitness effects of introgressed alleles, there is a risk of maladaptive allele introgression to the local bumblebee subspecies, which can negatively impact autochthon populations. One immediate recommendation to farmers is for the proper disposal of hive boxes, after their use in greenhouses, so as to minimize the risk of escapees contaminating native populations. On the other hand, the feasibility of using local subspecies B. t. lusitanicus, preferably with local production, should be evaluated.

Long-Distance Benefits of Marine Reserves: Myth or Reality?

Long-distance (>40-km) dispersal from marine reserves is poorly documented; yet, it can provide essential benefits such as seeding fished areas or connecting marine reserves into networks. From a meta-analysis, we suggest that the spatial scale of marine connectivity is underestimated due to the limited geographic extent of sampling designs. We also found that the largest marine reserves (>1000km²) are the most isolated. These findings have important implications for the assessment of evolutionary, ecological, and socio-economic long-distance benefits of marine reserves. We conclude that existing methods to infer dispersal should consider the up-to-date genomic advances and also expand the spatial scale of sampling designs. Incorporating long-distance connectivity in conservation planning will contribute to increase the benefits of marine reserve networks.

Local adaptation in dispersal in multi-resource landscapes

The distribution of resources in space has important consequences for the evolution of dispersal-related traits. Dispersal moderates patterns of gene flow and, consequently, the potential for local adaptation to spatially differentiated resource types. We lack both models and experiments that evaluate how dispersal evolves in landscapes with multiple resources. Here, we investigate the evolution of dispersal in landscapes that contain two resource types that differ in their spatial autocorrelations. Individuals may possess ecological traits that give them a fitness advantage on one or the other resource. We find that resources differing in their spatial autocorrelation select for different optimal dispersal strategies and, further, that some multi-resource landscapes can support the stable coexistence of distinct dispersal strategies. Whether divergence in dispersal strategies between resource specialists occurs depends on the underlying structure of the resources and the degree of linkage between dispersal strategies and ecological specialization. This work indicates that the spatial autocorrelation of resources is an important factor in determining when evolutionary branching is likely to occur, and sheds light on when secondary isolating mechanisms should arise between locally adapted specialists.