A historical stepping-stone path for an island-colonizing cactus across a submerged "bridge" archipelago

Here we use population genomic data (ddRAD-Seq) and ecological niche modeling to test biogeographic hypotheses for the divergence of the island-endemic cactus species Cereus insularis Hemsl. (Cereeae; Cactaceae) from its sister species C. fernambucensis Lem. The Cereus insularis grows in the Fernando de Noronha Islands (FNI), a Neotropical archipelago located 350 km off the Brazilian Atlantic Forest (BAF) coast. Phylogeographic reconstructions support a northward expansion by the common ancestor of C. insularis and C. fernambucensis along the mainland BAF coast, with C. insularis diverging from the widespread mainland taxon C. fernambucensis after colonizing FNI in the late Pleistocene. The morphologically distinct C. insularis is monophyletic and nested within C. fernambucensis, as expected from a progenitor-derivative speciation model. We tested alternative biogeographic and demographic hypotheses for the colonization of the FNI using Approximate Bayesian Computation. We found the greatest support for a stepping-stone path that emerged during periods of decreased sea level (the "bridge" hypothesis), in congruence with historical ecological niche modeling that shows highly suitable habitats on stepping-stone islands during glacial periods. The outlier analyses reveal signatures of selection in C. insularis, suggesting a putative role of adaptation driving rapid anagenic differentiation of this species in FNI.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.

Anthropogenic Change and the Process of Speciation

Anthropogenic impacts on the environment alter speciation processes by affecting both geographical contexts and selection patterns on a worldwide scale. Here we review evidence of these effects. We find that human activities often generate spatial isolation between populations and thereby promote genetic divergence but also frequently cause sudden secondary contact and hybridization between diverging lineages. Human-caused environmental changes produce new ecological niches, altering selection in diverse ways that can drive diversification; but changes also often remove niches and cause extirpations. Human impacts that alter selection regimes are widespread and strong in magnitude, ranging from local changes in biotic and abiotic conditions to direct harvesting to global climate change. Altered selection, and evolutionary responses to it, impacts early-stage divergence of lineages, but does not necessarily lead toward speciation and persistence of separate species. Altogether, humans both promote and hinder speciation, although new species would form very slowly relative to anthropogenic hybridization, which can be nearly instantaneous. Speculating about the future of speciation, we highlight two key conclusions: (1) Humans will have a large influence on extinction and "despeciation" dynamics in the short term and on early-stage lineage divergence, and thus potentially speciation in the longer term, and (2) long-term monitoring combined with easily dated anthropogenic changes will improve our understanding of the processes of speciation. We can use this knowledge to preserve and restore ecosystems in ways that promote (re-)diversification, increasing future opportunities of speciation and enhancing biodiversity.

Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species

The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-023-09612-5.

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.

Trophic specialization on unique resources despite limited niche divergence in a celebrated example of sympatric speciation

Trophic niche partitioning is observed in many adaptive radiations and is hypothesized to be a central process underlying species divergence. However, patterns of dietary niche partitioning are inconsistent across radiations and there are few studies of niche partitioning in putative examples of sympatric speciation. Here, we conducted the first quantitative study of dietary niche partitioning using stomach contents and stable isotope analyses in one of the most celebrated examples of sympatric speciation: the cichlid radiation from crater lake Barombi Mbo, Cameroon. We found little evidence for trophic niche partitioning among cichlids, including the nine species coexisting in the narrow littoral zone. Stable isotope analyses supported these conclusions of substantial dietary overlap. Our data, however, did reveal that five of eleven species consume rare dietary items, including freshwater sponge, terrestrial ants, and nocturnal foraging on shrimp. Stomach contents of the spongivore (Pungu maclareni) were 20% freshwater sponge, notable considering that only 0.04% of all fishes consume sponges. Overall, we conclude that cichlid species in lake Barombi Mbo overlap considerably in broad dietary niches-in part due to the large proportion of detritus in the stomach contents of all species-but there is evidence for divergence among species in their diet specializations on unique resources. We speculate that these species may utilize these additional specialized resources during periods of low resource abundance in support of Liem's paradox.

Seascape genomics identify adaptive barriers correlated to tidal amplitude in the shore crab Carcinus maenas

Most marine invertebrates disperse during a planktonic larval stage that may drift for weeks with ocean currents. A challenge for larvae of coastal species is to return to coastal nursery habitats. Shore crab (Carcinus maenas L.) larvae are known to show tidal rhythmicity in vertical migration in tidal areas and circadian rhythmicity in microtidal areas, which seems to increase successful coastal settlement. We studied genome-wide differentiation based on 24,000 single nucleotide polymorphisms of 12 native populations of shore crab sampled from a large tidal amplitude gradient from macrotidal (~8 m) to microtidal (~0.2 m). Dispersal and recruitment success of larvae was assessed with a Lagrangian biophysical model, which showed a strong effect of larval behaviour on long-term connectivity, and dispersal barriers that partly coincided with different tidal environments. The genetic population structure showed a subdivision of the samples into three clusters, which represent micro-, meso- and macrotidal areas. The genetic differentiation was mostly driven by 0.5% outlier loci, which showed strong allelic clines located at the limits between the three tidal areas. Demographic modelling suggested that the two genetic barriers have different origins. Differential gene expression of two clock genes (cyc and pdp1) further highlighted phenotypic differences among genetic clusters that are potentially linked to the differences in larval behaviour. Taken together, our seascape genomic study suggests that tidal regime acts as a strong selection force on shore crab population structure, consistent with larval behaviour affecting dispersal and recruitment success.

Limited dispersal and local adaptation promote allopatric speciation in a biodiversity hotspot

Recently diverged or diverging populations can offer unobstructed insights into early barriers to gene flow during the initial stages of speciation. The current study utilised a novel insect system (order Mantophasmatodea) to shed light on the early drivers of speciation. The members of this group have limited dispersal abilities, small allopatric distributions and strong habitat associations in the Cape Floristic Region biodiversity hotspot in South Africa. Sister taxa from the diverse family Austrophasmatidae were chosen as focal species (Karoophasma biedouwense, K. botterkloofense). Population genetics and Generalized Dissimilarity Modelling (GDM) were used to characterise spatial patterns of genetic variation and evaluate the contribution of environmental factors to population divergence and speciation. Extensive sampling confirmed the suspected allopatry of these taxa. However, hybrids were identified in a narrow region occurring between the species' distributions. Strong population structure was found over short geographic distances; particularly in K. biedouwense in which geographic distance accounted for 32% of genetic variation over a scale of 50 km (r = .56, p < .001). GDM explained 42%-78% of the deviance in observed genetic dissimilarities. Geographic distance was consistently indicated to be important for between species and within population differentiation, suggesting that limited dispersal ability may be an important neutral driver of divergence. Temperature, altitude, precipitation and vegetation were also indicated as important factors, suggesting the possible role of adaptation to local environmental conditions for species divergence. The discovery of the hybrid-zone, and the multiple allopatric species pairs in Austrophasmatidae support the idea that this could be a promising group to further our understanding of speciation modes.

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.

Speciation across the Tree of Life

Much of what we know about speciation comes from detailed studies of well-known model systems. Although there have been several important syntheses on speciation, few (if any) have explicitly compared speciation among major groups across the Tree of Life. Here, we synthesize and compare what is known about key aspects of speciation across taxa, including bacteria, protists, fungi, plants, and major animal groups. We focus on three main questions. Is allopatric speciation predominant across groups? How common is ecological divergence of sister species (a requirement for ecological speciation), and on what niche axes do species diverge in each group? What are the reproductive isolating barriers in each group? Our review suggests the following patterns. (i) Based on our survey and projected species numbers, the most frequent speciation process across the Tree of Life may be co-speciation between endosymbiotic bacteria and their insect hosts. (ii) Allopatric speciation appears to be present in all major groups, and may be the most common mode in both animals and plants, based on non-overlapping ranges of sister species. (iii) Full sympatry of sister species is also widespread, and may be more common in fungi than allopatry. (iv) Full sympatry of sister species is more common in some marine animals than in terrestrial and freshwater ones. (v) Ecological divergence of sister species is widespread in all groups, including ~70% of surveyed species pairs of plants and insects. (vi) Major axes of ecological divergence involve species interactions (e.g. host-switching) and habitat divergence. (vii) Prezygotic isolation appears to be generally more widespread and important than postzygotic isolation. (viii) Rates of diversification (and presumably speciation) are strikingly different across groups, with the fastest rates in plants, and successively slower rates in animals, fungi, and protists, with the slowest rates in prokaryotes. Overall, our study represents an initial step towards understanding general patterns in speciation across all organisms.

Genetic and phenotypic variation exhibit both predictable and stochastic patterns across an intertidal fish metapopulation

Interactions among selection, gene flow, and drift affect the trajectory of adaptive evolution. In natural populations, the direction and magnitude of these processes can be variable across different spatial, temporal, or ontogenetic scales. Consequently, variability in evolutionary processes affects the predictability or stochasticity of microevolutionary outcomes. We studied an intertidal fish, Bathygobius cocosensis (Bleeker, 1854), to understand how space, time, and life stage structure genetic and phenotypic variation in a species with potentially extensive dispersal and a complex life cycle (larval dispersal preceding benthic recruitment). We sampled juvenile and adult life stages, at three sites, over three years. Genome-wide SNPs uncovered a pattern of chaotic genetic patchiness, that is, weak-but-significant patchy spatial genetic structure that was variable through time and between life stages. Outlier locus analyses suggested that targets of spatially divergent selection were mostly temporally variable, though a significant number of spatial outlier loci were shared between life stages. Head shape, a putatively ecologically responsive (adaptive) phenotype in B. cocosensis also exhibited high temporal variability within sites. However, consistent spatial relationships between sites indicated that environmental similarities among sites may generate predictable phenotype distributions across space. Our study highlights the complex microevolutionary dynamics of marine systems, where consideration of multiple ecological dimensions can reveal both predictable and stochastic patterns in the distributions of genetic and phenotypic variation. Such considerations probably apply to species that possess short, complex life cycles, have large dispersal potential and fecundities, and that inhabit heterogeneous environments.

Adaptive zones shape the magnitude of premating reproductive isolation in Timema stick insects

Simpson's fossil-record inspired model of 'adaptive zones' proposes that evolution is dominated by small fluctuations within adaptive zones, occasionally punctuated by larger shifts between zones. This model can help explain why the process of population divergence often results in weak or moderate reproductive isolation (RI), rather than strong RI and distinct species. Applied to the speciation process, the adaptive zones hypothesis makes two inter-related predictions: (i) large shifts between zones are relatively rare, (ii) when large shifts do occur they generate stronger RI than shifts within zones. Here, we use ecological, phylogenetic and behavioural data to test these predictions in Timema stick insects. We show that host use in Timema is dominated by moderate shifts within the systematic divisions of flowering plants and conifers, with only a few extreme shifts between these divisions. However, when extreme shifts occur, they generate greater RI than do more moderate shifts. Our results support the adaptive zones model, and suggest that the net contribution of ecological shifts to diversification is dependent on both their magnitude and frequency. We discuss the generality of our findings in the light of emerging evidence from diverse taxa that the evolution of RI is not always the only factor determining the origin of species diversity. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

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.

The Contemporary Evolution of Fitness

The rate of evolution of population mean fitness informs how selection acting in contemporary populations can counteract environmental change and genetic degradation (mutation, gene flow, drift, recombination). This rate influences population increases (e.g., range expansion), population stability (e.g., cryptic eco-evolutionary dynamics), and population recovery (i.e., evolutionary rescue). We review approaches for estimating such rates, especially in wild populations. We then review empirical estimates derived from two approaches: mutation accumulation (MA) and additive genetic variance in fitness (IAw). MA studies inform how selection counters genetic degradation arising from deleterious mutations, typically generating estimates of <1% per generation. IAw studies provide an integrated prediction of proportional change per generation, nearly always generating estimates of <20% and, more typically, <10%. Overall, considerable, but not unlimited, evolutionary potential exists in populations facing detrimental environmental or genetic change. However, further studies with diverse methods and species are required for more robust and general insights.

Divergence in cryptic leaf colour provides local camouflage in an alpine plant

The efficacy of camouflage through background matching is highly environment-dependent, often resulting in intraspecific colour divergence in animals to optimize crypsis in different visual environments. This phenomenon is largely unexplored in plants, although several lines of evidence suggest they do use crypsis to avoid damage by herbivores. Using Corydalis hemidicentra, an alpine plant with cryptic leaf colour, we quantified background matching between leaves and surrounding rocks in five populations based on an approximate model of their butterfly enemy's colour perception. We also investigated the pigment basis of leaf colour variation and the association between feeding risk and camouflage efficacy. We show that plants exhibit remarkable colour divergence between populations, consistent with differences in rock appearances. Leaf colour varies because of a different quantitative combination of two basic pigments-chlorophyll and anthocyanin-plus different air spaces. As expected, leaf colours are better matched against their native backgrounds than against foreign ones in the eyes of the butterfly. Furthermore, improved crypsis tends to be associated with a higher level of feeding risk. These results suggest that divergent cryptic leaf colour may have evolved to optimize local camouflage in various visual environments, extending our understanding of colour evolution and intraspecific phenotype diversity in plants.

Gene flow and metacommunity arrangement affects coevolutionary dynamics at the mutualism-antagonism interface

Interspecific interactions are affected by community context and, as a consequence, show spatial variation in magnitude and sign. The selective forces imposed by interactions at the mutualism-antagonism interface are a consequence of the traits involved and their matching between species. If mutualistic and antagonistic communities are linked by gene flow, coevolution between a pair of interacting species is influenced by how selection varies in space. Here we investigate the effects of metacommunity arrangement, i.e. patterns of connection between communities and the number of communities, on the coevolutionary dynamics between two species for which the sign and magnitude of the interaction varies across the landscape. We quantify coevolutionary outcome as an index that can be decomposed into the contribution of intraspecific genetic diversity and interspecific interaction. We show that polymorphisms and mismatches are an expected outcome, which is influenced by spatial structure, interaction strength and the degree of gene flow. The index describes how variation is distributed within and between species, and provides information on the directionality of the mismatches and polymorphisms. Finally, we argue that depending on metacommunity arrangement, some communities have disproportionate roles in maintaining genetic diversity, with implications for the coevolution of interacting species in a fragmented landscape.

Observational evidence that maladaptive gene flow reduces patch occupancy in a wild insect metapopulation

Theory predicts that dispersal throughout metapopulations has a variety of consequences for the abundance and distribution of species. Immigration is predicted to increase abundance and habitat patch occupancy, but gene flow can have both positive and negative demographic consequences. Here, we address the eco-evolutionary effects of dispersal in a wild metapopulation of the stick insect Timema cristinae, which exhibits variable degrees of local adaptation throughout a heterogeneous habitat patch network of two host-plant species. To disentangle the ecological and evolutionary contributions of dispersal to habitat patch occupancy and abundance, we contrasted the effects of connectivity to populations inhabiting conspecific host plants and those inhabiting the alternate host plant. Both types of connectivity should increase patch occupancy and abundance through increased immigration and sharing of beneficial alleles through gene flow. However, connectivity to populations inhabiting the alternate host-plant species may uniquely cause maladaptive gene flow that counters the positive demographic effects of immigration. Supporting these predictions, we find the relationship between patch occupancy and alternate-host connectivity to be significantly smaller in slope than the relationship between patch occupancy and conspecific-host connectivity. Our findings illustrate the ecological and evolutionary roles of dispersal in driving the distribution and abundance of species.

Spatial pattern of invasion and the evolutionary responses of native plant species

Invasive plant species can have a strong negative impact on the resident native species, likely imposing new selective pressures on them. Altered selective pressures may result in evolutionary changes in some native species, reducing competitive exclusion and allowing for coexistence with the invader. Native genotypes that are able to coexist with strong invaders may represent a valuable resource for management efforts. A better understanding of the conditions under which native species are more, or less, likely to adapt to an invader is necessary to incorporate these eco-evolutionary dynamics into management strategies. We propose that the spatial structure of invasion, in particular the size and isolation of invaded patches, is one factor which can influence the evolutionary responses of native species through modifying gene flow and the strength of selection. We present a conceptual model in which large, dense, and well-connected patches result in a greater likelihood of native species adaptation. We also identify characteristics of the interacting species that may influence the evolutionary response of native species to invasion and outline potential management implications. Identifying areas of rapid evolutionary change may offer one additional tool to managers in their effort to conserve biodiversity in the face of invasion.

Navigating the pitfalls and promise of landscape genetics

The field of landscape genetics has been evolving rapidly since its emergence in the early 2000s. New applications, techniques and criticisms of techniques appear like clockwork with each new journal issue. The developments are an encouraging, and at times bewildering, sign of progress in an exciting new field of study. However, we suggest that the rapid expansion of landscape genetics has belied important flaws in the development of the field, and we add an air of caution to this breakneck pace of expansion. Specifically, landscape genetic studies often lose sight of the fundamental principles and complex consequences of gene flow, instead favouring simplistic interpretations and broad inferences not necessarily warranted by the data. Here, we describe common pitfalls that characterize such studies, and provide practical guidance to improve landscape genetic investigation, with careful consideration of inferential limits, scale, replication, and the ecological and evolutionary context of spatial genetic patterns. Ultimately, the utility of landscape genetics will depend on translating the relationship between gene flow and landscape features into an understanding of long-term population outcomes. We hope the perspective presented here will steer landscape genetics down a more scientifically sound and productive path, garnering a field that is as informative in the future as it is popular now.

Genetic and morphological consequences of Quaternary glaciations: A relic barbel lineage (Luciobarbus pallaryi, Cyprinidae) of Guir Basin (Algeria)

Climatic variations during the Quaternary period had a considerable impact on landscapes and habitat fragmentation (rivers) in North Africa. These historical events can have significant consequences on the genetic structure of the populations. Indeed, geographically separated and genetically isolated populations tend to differentiate themselves through time, eventually becoming distinct lineages, allowing new species to emerge in later generations. The aim of the present study is to use genetic and morphological techniques to evaluate the major role of the Saalian glaciation (Middle Quaternary) in the establishment of the geographic space and in the evolution of the intraspecific genetic diversity, by tracing the demographic history of barbels belonging to the Luciobarbus pallaryi (Cyprinidae) species in the Guir Basin (Algeria). In this context, two populations, from two distinct and isolated sites, were studied. Analysis of the cytochrome b (cyt b) mitochondrial markers and of the "D-loop" control region has shown that the "upstream" and "downstream" Guir populations are genetically differentiated. The molecular analyses suggest that the upstream population was disconnected from this hydrographic system during the Saalian glaciation period of the Quaternary. Subsequently, it was isolated in the foggaras underground waters in the Great Western Erg, at approximately 320 000 years BP, creating, through a bottleneck effect, a new allopatric lineage referred to as "Adrar". Conversely, the high genetic diversity in the upstream Guir (Bechar) population suggests that the stock is globally in expansion. These barbels (n=52) were also examined with meristic, morphometric, osteological, and biological features. These data also reveal a complete discrimination between the two populations, with a remarkable and distinctive behavioural adaptation for the Adrar specimens: neoteny.

The road to opportunities: landscape change promotes body-size divergence in a highly mobile species

Landscape change provides a suitable framework for investigating population-level responses to novel ecological pressures. However, relatively little attention has been paid to examine the potential influence of landscape change on the geographic scale of population differentiation. Here, we tested for morphological differentiation of red-necked nightjars Caprimulgus ruficollis breeding in a managed property and a natural reserve situated less than 10 km apart. At both sites, we also estimated site fidelity over 5 years and quantified the potential foraging opportunities for nightjars. Breeding birds in the managed habitat were significantly larger in size—as indexed by keel length—than those in the natural one. However, there were no significant differences in wing or tail length. Immigration from neighboring areas was almost negligible and, furthermore, no individual (out of 1130 captures overall) exchanged habitats between years, indicating strong site fidelity. Food supply for nightjars was equally abundant in both habitats, but the availability of foraging sites was remarkably higher in the managed property. As a result, nightjars—particularly fledglings—in the latter habitat benefited from increased foraging opportunities in relation to those in the natural site. It seems likely that the fine-scale variation in nightjar morphology reflects a phenotypic response to unequal local conditions, since non-random dispersal or differential mortality had been determined not to be influential. High site fidelity appears to contribute to the maintenance of body-size differences between the two habitats. Results from this nightjar population highlight the potential of human-induced landscape change to promote population-level responses at exceedingly small geographic scales.

Population genomics of divergence among extreme and intermediate color forms in a polymorphic insect

Geographic variation in insect coloration is among the most intriguing examples of rapid phenotypic evolution and provides opportunities to study mechanisms of phenotypic change and diversification in closely related lineages. The bumble bee Bombus bifarius comprises two geographically disparate color groups characterized by red‐banded and black‐banded abdominal pigmentation, but with a range of spatially and phenotypically intermediate populations across western North America. Microsatellite analyses have revealed that B. bifarius in the USA are structured into two major groups concordant with geography and color pattern, but also suggest ongoing gene flow among regional populations. In this study, we better resolve the relationships among major color groups to better understand evolutionary mechanisms promoting and maintaining such polymorphism. We analyze >90,000 and >25,000 single‐nucleotide polymorphisms derived from transcriptome (RNAseq) and double digest restriction site associated DNA sequencing (ddRAD), respectively, in representative samples from spatial and color pattern extremes in B. bifarius as well as phenotypic and geographic intermediates. Both ddRAD and RNAseq data illustrate substantial genome‐wide differentiation of the red‐banded (eastern) color form from both black‐banded (western) and intermediate (central) phenotypes and negligible differentiation among the latter populations, with no obvious admixture among bees from the two major lineages. Results thus indicate much stronger background differentiation among B. bifarius lineages than expected, highlighting potential challenges for revealing loci underlying color polymorphism from population genetic data alone. These findings will have significance for resolving taxonomic confusion in this species and in future efforts to investigate color‐pattern evolution in B. bifarius and other polymorphic bumble bee species.

Evolutionary processes driving spatial patterns of intraspecific genetic diversity in river ecosystems

Describing, understanding and predicting the spatial distribution of genetic diversity is a central issue in biological sciences. In river landscapes, it is generally predicted that neutral genetic diversity should increase downstream, but there have been few attempts to test and validate this assumption across taxonomic groups. Moreover, it is still unclear what are the evolutionary processes that may generate this apparent spatial pattern of diversity. Here, we quantitatively synthesized published results from diverse taxa living in river ecosystems, and we performed a meta-analysis to show that a downstream increase in intraspecific genetic diversity (DIGD) actually constitutes a general spatial pattern of biodiversity that is repeatable across taxa. We further demonstrated that DIGD was stronger for strictly waterborne dispersing than for overland dispersing species. However, for a restricted data set focusing on fishes, there was no evidence that DIGD was related to particular species traits. We then searched for general processes underlying DIGD by simulating genetic data in dendritic-like river systems. Simulations revealed that the three processes we considered (downstream-biased dispersal, increase in habitat availability downstream and upstream-directed colonization) might generate DIGD. Using random forest models, we identified from simulations a set of highly informative summary statistics allowing discriminating among the processes causing DIGD. Finally, combining these discriminant statistics and approximate Bayesian computations on a set of twelve empirical case studies, we hypothesized that DIGD were most likely due to the interaction of two of these three processes and that contrary to expectation, they were not solely caused by downstream-biased dispersal.

Signatures of selection in the three-spined stickleback along a small-scale brackish water - freshwater transition zone

Local adaptation is often obvious when gene flow is impeded, such as observed at large spatial scales and across strong ecological contrasts. However, it becomes less certain at small scales such as between adjacent populations or across weak ecological contrasts, when gene flow is strong. While studies on genomic adaptation tend to focus on the former, less is known about the genomic targets of natural selection in the latter situation. In this study, we investigate genomic adaptation in populations of the three-spined stickleback Gasterosteus aculeatus L. across a small-scale ecological transition with salinities ranging from brackish to fresh. Adaptation to salinity has been repeatedly demonstrated in this species. A genome scan based on 87 microsatellite markers revealed only few signatures of selection, likely owing to the constraints that homogenizing gene flow puts on adaptive divergence. However, the detected loci appear repeatedly as targets of selection in similar studies of genomic adaptation in the three-spined stickleback. We conclude that the signature of genomic selection in the face of strong gene flow is weak, yet detectable. We argue that the range of studies of genomic divergence should be extended to include more systems characterized by limited geographical and ecological isolation, which is often a realistic setting in nature.

How maladaptation can structure biodiversity: eco-evolutionary island biogeography

Current research on eco-evolutionary dynamics is mainly concerned with understanding the role of rapid (or 'contemporary') evolution in structuring ecological patterns. We argue that the current eco-evolutionary research program, which focuses largely on natural selection, should be expanded to more explicitly consider other evolutionary processes such as gene flow. Because multiple evolutionary processes interact to generate quantitative variation in the degree of local maladaptation, we focus on how studying the ecological effects of maladaptation will lead to a more comprehensive view of how evolution can influence ecology. We explore how maladaptation can influence ecology through the lens of island biogeography theory, which yields some novel predictions, such as patch isolation increasing species richness.

Locally adapted traits maintained in the face of high gene flow

Gene flow between phenotypically divergent populations can disrupt local adaptation or, alternatively, may stimulate adaptive evolution by increasing genetic variation. We capitalised on historical Trinidadian guppy transplant experiments to test the phenotypic effects of increased gene flow caused by replicated introductions of adaptively divergent guppies, which were translocated from high- to low-predation environments. We sampled two native populations prior to the onset of gene flow, six historic introduction sites, introduction sources and multiple downstream points in each basin. Extensive gene flow from introductions occurred in all streams, yet adaptive phenotypic divergence across a gradient in predation level was maintained. Descendants of guppies from a high-predation source site showed high phenotypic similarity with native low-predation guppies in as few as ~12 generations after gene flow, likely through a combination of adaptive evolution and phenotypic plasticity. Our results demonstrate that locally adapted phenotypes can be maintained despite extensive gene flow from divergent populations.

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.

Morphological divergence driven by predation environment within and between species of Brachyrhaphis fishes

Natural selection often results in profound differences in body shape among populations from divergent selective environments. Predation is a well-studied driver of divergence, with predators having a strong effect on the evolution of prey body shape, especially for traits related to escape behavior. Comparative studies, both at the population level and between species, show that the presence or absence of predators can alter prey morphology. Although this pattern is well documented in various species or population pairs, few studies have tested for similar patterns of body shape evolution at multiple stages of divergence within a taxonomic group. Here, we examine morphological divergence associated with predation environment in the livebearing fish genus Brachyrhaphis. We compare differences in body shape between populations of B. rhabdophora from different predation environments to differences in body shape between B. roseni and B. terrabensis (sister species) from predator and predator free habitats, respectively. We found that in each lineage, shape differed between predation environments, consistent with the hypothesis that locomotor function is optimized for either steady swimming (predator free) or escape behavior (predator). Although differences in body shape were greatest between B. roseni and B. terrabensis, we found that much of the total morphological diversification between these species had already been achieved within B. rhabdophora (29% in females and 47% in males). Interestingly, at both levels of divergence we found that early in ontogenetic development, females differed in shape between predation environments; however, as females matured, their body shapes converged on a similar phenotype, likely due to the constraints of pregnancy. Finally, we found that body shape varies with body size in a similar way, regardless of predation environment, in each lineage. Our findings are important because they provide evidence that the same source of selection can drive similar phenotypic divergence independently at multiple divergence levels.

Evolution of camouflage drives rapid ecological change in an insect community

BACKGROUND

Evolutionary change in individual species has been hypothesized to have far-reaching consequences for entire ecological communities, and such coupling of ecological and evolutionary dynamics ("eco-evolutionary dynamics") has been demonstrated for a variety systems. However, the general importance of evolutionary dynamics for ecological dynamics remains unclear. Here, we investigate how spatial patterns of local adaptation in the stick insect Timema cristinae, driven by the interaction between multiple evolutionary processes, structure metapopulations, communities, and multitrophic interactions.

RESULTS

Observations of a wild T. cristinae metapopulation show that locally imperfect camouflage reduces population size and that the effect of such maladaptation is comparable to the effects of more traditional ecological factors, including habitat patch size and host-plant species identity. Field manipulations of local adaptation and bird predation support the hypothesis that maladaptation reduces population size through an increase in bird predation. Furthermore, these field experiments show that maladaptation in T. cristinae and consequent increase in bird predation reduce the pooled abundance and species richness of the co-occurring arthropod community, and ultimately cascade to decrease herbivory on host plants. An eco-evolutionary model of the observational data demonstrates that the demographic cost of maladaptation decreases habitat patch occupancy by T. cristinae but enhances metapopulation-level adaptation.

CONCLUSIONS

The results demonstrate a pervasive effect of ongoing evolution in a spatial context on population and community dynamics. The eco-evolutionary model makes testable predictions about the influence of the spatial configuration of the patch network on metapopulation size and the spatial scale of adaptation.

Sexual dimorphism dominates divergent host plant use in stick insect trophic morphology

BACKGROUND

Clear examples of ecological speciation exist, often involving divergence in trophic morphology. However, substantial variation also exists in how far the ecological speciation process proceeds, potentially linked to the number of ecological axes, traits, or genes subject to divergent selection. In addition, recent studies highlight how differentiation might occur between the sexes, rather than between populations. We examine variation in trophic morphology in two host-plant ecotypes of walking-stick insects (Timema cristinae), known to have diverged in morphological traits related to crypsis and predator avoidance, and to have reached an intermediate point in the ecological speciation process. Here we test how host plant use, sex, and rearing environment affect variation in trophic morphology in this species using traditional multivariate, novel kernel density based and Bayesian morphometric analyses.

RESULTS

Contrary to expectations, we find limited host-associated divergence in mandible shape. Instead, the main predictor of shape variation is sex, with secondary roles of population of origin and rearing environment.

CONCLUSION

Our results show that trophic morphology does not strongly contribute to host-adapted ecotype divergence in T. cristinae and that traits can respond to complex selection regimes by diverging along different intraspecific lines, thereby impeding progress toward speciation.

Adaptive divergence in seed color camouflage in contrasting soil environments

Although adaptive plant population divergence across contrasting soil conditions is often driven by abiotic soil factors, natural enemies may also contribute. Cryptic matching to the native soil color is a form of defensive camouflage that seeds can use to avoid detection by seed predators. The legume Acmispon wrangelianus occurs across a variety of gray-green serpentine soils and brown nonserpentine soils. Quantitative digital image analysis of seed and soil colors was used to test whether genetically based seed color is a closer match to the color of the native soil than to the color of other nearby soils. Lineages bear seeds that more closely match the color of their native serpentine or nonserpentine soil type than the opposing soil type. Further, even within a soil type, lineages bear seeds with a closer color match to the soil at their native site than to other sites. The striking concordance between seed and native soil color suggests that natural selection for locally camouflaged seed color morphs, probably driven by seed predators, may maintain adaptive divergence in pigmentation, despite the opportunity for migration between soil environments.

Geographically multifarious phenotypic divergence during speciation

Speciation is an important evolutionary process that occurs when barriers to gene flow evolve between previously panmictic populations. Although individual barriers to gene flow have been studied extensively, we know relatively little regarding the number of barriers that isolate species or whether these barriers are polymorphic within species. Herein, we use a series of field and lab experiments to quantify phenotypic divergence and identify possible barriers to gene flow between the butterfly species Lycaeides idas and Lycaeides melissa. We found evidence that L. idas and L. melissa have diverged along multiple phenotypic axes. Specifically, we identified major phenotypic differences in female oviposition preference and diapause initiation, and more moderate divergence in mate preference. Multiple phenotypic differences might operate as barriers to gene flow, as shown by correlations between genetic distance and phenotypic divergence and patterns of phenotypic variation in admixed Lycaeides populations. Although some of these traits differed primarily between species (e.g., diapause initiation), several traits also varied among conspecific populations (e.g., male mate preference and oviposition preference).

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

Background

Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated (‘allopatric’) populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations.

Results

As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation.

Conclusions

Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of ‘genomic islands of speciation’.

Genomic consequences of multiple speciation processes in a stick insect

Diverse geographical modes and mechanisms of speciation are known, and individual speciation genes have now been identified. Despite this progress, genome-wide outcomes of different evolutionary processes during speciation are less understood. Here, we integrate ecological and spatial information, mating trials, transplantation data and analysis of 86 130 single nucleotide polymorphisms (SNPs) in eight populations (28 pairwise comparisons) of Timema cristinae stick insects to test the effects of different factors on genomic divergence in a system undergoing ecological speciation. We find patterns consistent with effects of numerous factors, including geographical distance, gene flow, divergence in host plant use and climate, and selection against maladaptive hybridization (i.e. reinforcement). For example, the number of highly differentiated ‘outlier loci’, allele-frequency clines and the overall distribution of genomic differentiation were recognizably affected by these factors. Although host use has strong effects on phenotypic divergence and reproductive isolation, its effects on genomic divergence were subtler and other factors had pronounced effects. The results demonstrate how genomic data can provide new insights into speciation and how genomic divergence can be complex, yet predictable. Future work could adopt experimental, mapping and functional approaches to directly test which genetic regions are affected by selection and determine their physical location in the genome.

Evidence of adaptation from ancestral variation in young populations of beach mice

To understand how organisms adapt to novel habitats, which involves both demographic and selective events, we require knowledge of the evolutionary history of populations and also selected alleles. There are still few cases in which the precise mutations (and hence, defined alleles) that contribute to adaptive change have been identified in nature; one exception is the genetic basis of camouflaging pigmentation of oldfield mice (Peromyscus polionotus) that have colonized the sandy dunes of Florida's Gulf Coast. To quantify the genomic impact of colonization as well as the signature of selection, we resequenced 5000 1.5-kb noncoding loci as well as a 160-kb genomic region surrounding the melanocortin-1 receptor (Mc1r), a gene that contributes to pigmentation differences, in beach and mainland populations. Using a genome-wide phylogenetic approach, we recovered a single monophyletic group comprised of beach mice, consistent with a single colonization event of the Gulf Coast. We also found evidence of a severe founder event, estimated to have occurred less than 3000 years ago. In this demographic context, we show that all beach subspecies share a single derived light Mc1r allele, which was likely selected from standing genetic variation that originated in the mainland. Surprisingly, we were unable to identify a clear signature of selection in the Mc1r region, despite independent evidence that this locus contributes to adaptive coloration. Nonetheless, these data allow us to reconstruct and compare the evolutionary history of populations and alleles to better understand how adaptive evolution, following the colonization of a novel habitat, proceeds in nature.

The evolution of sensory divergence in the context of limited gene flow in the bumblebee bat

The sensory drive theory of speciation predicts that populations of the same species inhabiting different environments can differ in sensory traits, and that this sensory difference can ultimately drive speciation. However, even in the best-known examples of sensory ecology driven speciation, it is uncertain whether the variation in sensory traits is the cause or the consequence of a reduction in levels of gene flow. Here we show strong genetic differentiation, no gene flow and large echolocation differences between the allopatric Myanmar and Thai populations of the world's smallest mammal, Craseonycteris thonglongyai, and suggest that geographic isolation most likely preceded sensory divergence. Within the geographically continuous Thai population, we show that geographic distance has a primary role in limiting gene flow rather than echolocation divergence. In line with sensory-driven speciation models, we suggest that in C. thonglongyai, limited gene flow creates the suitable conditions that favour the evolution of sensory divergence via local adaptation.

Populations of the same species living in different habitats can differ in sensory traits driving speciation, but it is not known if this variation limits gene flow. Here, a genetic and acoustic study of the bumblebee bat suggests that geographic distance, instead of echolocation divergence, limits gene flow.

EVOLUTION OF PREY POLYMORPHISM INDUCED BY LEARNING PREDATORS

A prey species using crypsis to avoid predators has the opportunity to evolve polymorphic crypsis when it is being exposed to two (or more) habitats with different backgrounds. Here, we investigate when this phenomenon can occur, in a simulation study with a sexually reproducing prey and a predator that can learn to find hiding prey, represented by an artificial neural network. Initially, the prey is well adapted to one habitat, but tries to expand its range by invading another, different, habitat. This can cause the prey to evolve toward an intermediate phenotype, equally cryptic in both habitats. The prey can also fail in adapting to its new environment, and stay the same. Alternatively, it can evolve polymorphic crypsis. We find that the evolutionary outcome depends on the amount of dispersal between the habitats, with polymorphic crypsis evolving for low dispersal rates, an intermediate phenotype will evolve for intermediate dispersal rates and no adaptation to the new habitat will occur for high dispersal rates. The distribution of phenotypes of the prey will also vary for different dispersal rates, with narrow distributions for low and high dispersal rate and a wide distribution for intermediate dispersal rates.

Fusion-fission experiments in Aphidius: evolutionary split without isolation in response to environmental bimodality

Studying host-based divergence naturally maintained by a balance between selection and gene flow can provide valuable insights into genetic underpinnings of host adaptation and ecological speciation in parasites. Selection-gene flow balance is often postulated in sympatric host races, but direct experimental evidence is scarce. In this study, we present such evidence obtained in host races of Aphidius ervi, an important hymenopteran agent of biological control of aphids in agriculture, using a novel fusion-fission method of gene flow perturbation. In our study, between-race genetic divergence was obliterated by means of advanced hybridisation, followed by a multi-generation exposure of the resulting genetically uniform hybrid swarm to a two-host environment. This fusion-fission procedure was implemented under two contrasting regimes of between-host gene flow in two replicated experiments involving different racial pairs. Host-based genetic fission in response to environmental bimodality occurred in both experiments in as little as six generations of divergent adaptation despite continuous gene flow. We demonstrate that fission recovery of host-based divergence evolved faster and hybridisation-induced linkage disequilibrium decayed slower under restricted (6.7%) compared with unrestricted gene flow, directly pointing at a balance between gene flow and divergent selection. We also show, in four separate tests, that random drift had no or little role in the observed genetic split. Rates and patterns of fission divergence differed between racial pairs. Comparative linkage analysis of these differences is currently under way to test for the role of genomic architecture of adaptation in ecology-driven divergent evolution.