Eco-evolutionary Dynamics Research Articles

Sharp habitat shifts, evolutionary tipping points and rescue: Quantifying the perilous path of a specialist species towards a refugium in a changing environment

Specialist species thriving under specific environmental conditions in narrow geographic ranges are widely recognized as heavily threatened by climate deregulation. Many might rely on both their potential to adapt and to disperse towards a refugium to avoid extinction. It is thus crucial to understand the influence of environmental conditions on the unfolding process of adaptation. Here, I study the eco-evolutionary dynamics of a sexually reproducing specialist species in a two-patch quantitative genetic model with moving optima. Thanks to a separation of ecological and evolutionary time scales and the phase-line study of the selection gradient, I derive the critical environmental speed for persistence, which reflects how the existence of a refugium impacts extinction patterns and how it relates to the cost of dispersal. Moreover, the analysis provides key insights about the dynamics that arise on the path towards this refugium. I show that after an initial increase of population size, there exists a critical environmental speed above which the species crosses a tipping point, resulting into an abrupt habitat switch. In addition, when selection for local adaptation is strong, this habitat switch passes through an evolutionary “death valley”, leading to a phenomenon related to evolutionary rescue, which can promote extinction for lower environmental speeds than the critical one.

The eco-evolutionary dynamics of two strategic species: From the predator-prey to the innocent-spreader rumor model

Species frequently engage in both competitive and cooperative interactions, delicately balancing these dynamics to optimize their chances of survival and reproduction. While competition drives individuals to compete for limited resources, cooperation can emerge as a strategic response, mitigating risk and enhancing collective payoff. To bridge theoretical game approaches such as payoff, cooperation, and defections in ecological systems, we propose a two-species predator–prey model inspired by the principles and variations of the prisoner’s dilemma game. We comprehensively address and analytically verify all stable strategic states, exploring the role of payoff parameters both individually and collectively. Additionally, we investigate the effect of free space. Beyond ecological contexts, we present a model of rumor propagation within a social system to establish connections with the prisoner’s dilemma game. In both systems, our primary focus is to discuss strategies and enhance the cooperative factor within the system, given its crucial importance across diverse environments.

Eco-evolutionary dynamics for finite populations and the noise-induced reversal of selection

Eco-evolutionary dynamics for finite populations and the noise-induced reversal of selection

A pangenomic atlas reveals eco-evolutionary dynamics that shape type VI secretion systems in plant-pathogenic Ralstonia.

Soilborne Ralstonia solanacearum species complex (RSSC) pathogens disrupt microbial communities as they invade roots and fatally wilt plants. RSSC pathogens secrete antimicrobial toxins using a type VI secretion system (T6SS). To investigate how evolution and ecology have shaped the T6SS of these bacterial pathogens, we analyzed the T6SS gene content and architecture across the RSSC and their evolutionary relatives. Our analysis reveals that two ecologically similar Burkholderiaceae taxa, xylem-pathogenic RSSC and Paracidovorax, have convergently evolved to wield large arsenals of T6SS toxins. To understand the mechanisms underlying genomic enrichment of T6SS toxins, we compiled an atlas of 1,066 auxiliary T6SS toxin clusters ("aux" clusters) across 99 high-quality RSSC genomes. We classified 25 types of aux clusters with toxins that predominantly target lipids, nucleic acids, or unknown cellular substrates. The aux clusters were located in diverse genetic neighborhoods and had complex phylogenetic distributions, suggesting frequent horizontal gene flow. Phages and other mobile genetic elements account for most of the aux cluster acquisition on the chromosome but very little on the megaplasmid. Nevertheless, RSSC genomes were more enriched in aux clusters on the megaplasmid. Although the single, ancestral T6SS was broadly conserved in the RSSC, the T6SS has been convergently lost in atypical, non-soilborne lineages. Overall, our data suggest dynamic interplay between the lifestyle of RSSC lineages and the evolution of T6SSes with robust arsenals of toxins. This pangenomic atlas poises the RSSC as an emerging, tractable model to understand the role of the T6SS in shaping pathogen populations.IMPORTANCEWe explored the eco-evolutionary dynamics that shape the inter-microbial warfare mechanisms of a globally significant plant pathogen, the Ralstonia solanacearum species complex. We discovered that most Ralstonia wilt pathogens have evolved extensive and diverse repertoires of type VI secretion system-associated antimicrobial toxins. These expansive toxin arsenals potentially enhance the ability of Ralstonia pathogens to invade plant microbiomes, enabling them to rapidly colonize and kill their host plants. We devised a classification system to categorize the Ralstonia toxins. Interestingly, many of the toxin gene clusters are encoded on mobile genetic elements, including prophages, which may be mutualistic symbionts that enhance the inter-microbial competitiveness of Ralstonia wilt pathogens. Moreover, our findings suggest that the convergent loss of this multi-gene trait contributes to genome reduction in two vector-transmitted lineages of Ralstonia pathogens. Our findings demonstrate that the interplay between microbial ecology and pathogen lifestyle shapes the evolution of a genetically complex antimicrobial weapon.

Parallel ecological and evolutionary responses to selection in a natural bacterial community

Evolution can occur over ecological timescales, suggesting a potentially important role for rapid evolution in shaping community trait distributions. However, evidence of concordant eco-evolutionary dynamics often comes from in vitro studies of highly simplified communities, and measures of ecological and evolutionary dynamics are rarely directly comparable. Here, we quantified how ecological species sorting and rapid evolution simultaneously shape community trait distributions by tracking within- and between-species changes in a key trait in a complex bacterial community. We focused on the production of siderophores; bacteria use these costly secreted metabolites to scavenge poorly soluble iron and to detoxify environments polluted with toxic nonferrous metals. We found that responses to copper-imposed selection within and between species were ultimately the same-intermediate siderophore levels were favored-and occurred over similar timescales. Despite being a social trait, this level of siderophore production was selected regardless of whether species evolved in isolation or in a community context. Our study suggests that evolutionary selection can play a pivotal role in shaping community trait distributions within natural, highly complex, bacterial communities. Furthermore, trait evolution may not always be qualitatively affected by interactions with other community members.

Improved high quality sand fly assemblies enabled by ultra low input long read sequencing

Phlebotomine sand flies are the vectors of leishmaniasis, a neglected tropical disease. High-quality reference genomes are an important tool for understanding the biology and eco-evolutionary dynamics underpinning disease epidemiology. Previous leishmaniasis vector reference sequences were limited by sequencing technologies available at the time and inadequate for high-resolution genomic inquiry. Here, we present updated reference assemblies of two sand flies, Phlebotomus papatasi and Lutzomyia longipalpis. These chromosome-level assemblies were generated using an ultra-low input library protocol, PacBio HiFi long reads, and Hi-C technology. The new P. papatasi reference has a final assembly span of 351.6 Mb and contig and scaffold N50s of 926 kb and 111.8 Mb, respectively. The new Lu. longipalpis reference has a final assembly span of 147.8 Mb and contig and scaffold N50s of 1.09 Mb and 40.6 Mb, respectively. Benchmarking Universal Single-Copy Orthologue (BUSCO) assessments indicated 94.5% and 95.6% complete single copy insecta orthologs for P. papatasi and Lu. longipalpis. These improved assemblies will serve as an invaluable resource for future genomic work on phlebotomine sandflies.

Eco-evolutionary dynamics of structured populations in periodically fluctuating environments: a G function approach.

Understanding the ecological and evolutionary dynamics of populations is critical for both basic and applied purposes in a variety of biological contexts. Although several modeling frameworks have been developed to simulate eco-evolutionary dynamics, many fewer address how to model structured populations. In a prior paper, we put forth the first modeling approach to simulate eco-evolutionary dynamics in structured populations under the G function modeling framework. However, this approach does not allow for accurate simulation under fluctuating environmental conditions. To address this limitation, we draw on the study of periodic differential equations to propose a modified approach that uses a different definition of fitness more suitable for fluctuating environments. We illustrate this method with a simple toy model of life history trade-offs. The generality of this approach allows it to be used in a variety of biological contexts.

Neutral genetic structuring of pathogen populations during rapid adaptation.

Pathogen species are experiencing strong joint demographic and selective events, especially when they adapt to a new host, for example through overcoming plant resistance. Stochasticity in the founding event and the associated demographic variations hinder our understanding of the expected evolutionary trajectories and the genetic structure emerging at both neutral and selected loci. What would be the typical genetic signatures of such a rapid adaptation event is not elucidated. Here, we build a demogenetic model to monitor pathogen population dynamics and genetic evolution on two host compartments (susceptible and resistant). We design our model to fit two plant pathogen life cycles, 'with' and 'without' host alternation. Our aim is to draw a typology of eco-evolutionary dynamics. Using time-series clustering, we identify three main scenarios: 1) small variations in the pathogen population size and small changes in genetic structure, 2) a strong founder event on the resistant host that in turn leads to the emergence of genetic structure on the susceptible host, and 3) evolutionary rescue that results in a strong founder event on the resistant host, preceded by a bot- tleneck on the susceptible host. We pinpoint differences between life cycles with notably more evolutionary rescue 'with' host alternation. Beyond the selective event itself, the demographic trajectory imposes specific changes in the genetic structure of the pathogen population. Most of these genetic changes are transient, with a signature of resistance overcoming that vanishes within a few years only. Considering time-series is therefore of utmost importance to accurately decipher pathogen evolution.

Regional diversity and leaf microbiome interactions of the fungal maize pathogen Exserohilum turcicum in Switzerland: A metagenomic analysis.

The spread and adaptation of fungal plant pathogens in agroecosystems are facilitated by environmental homogeneity. Metagenomic sequencing of infected tissues allowed us to monitor eco-evolutionary dynamics and interactions between host, pathogen and plant microbiome. Exserohilum turcicum, the causal agent of northern corn leaf blight (NCLB) in maize, is distributed in multiple clonal lineages throughout Europe. To characterize regional pathogen diversity, we conducted metagenomic DNA sequencing on 241 infected leaf samples from the highly susceptible Swiss maize landrace Rheintaler Ribelmais, collected over 3 years (2016-2018) from an average of 14 agricultural farms within the Swiss Rhine Valley. All major European clonal lineages of E. turcicum were identified. Lineages differ by their mating types which indicates potential for sexual recombination and rapid evolution of new pathogen strains, although we found no evidence of recent recombination. The associated eukaryotic and prokaryotic leaf microbiome exhibited variation in taxonomic diversity between years and locations and is likely influenced by local weather conditions. A network analysis revealed distinct clusters of eukaryotic and prokaryotic taxa that correlates with the frequency of E. turcicum sequencing reads, suggesting causal interactions. Notably, the yeast genus Metschnikowia exhibited a strongly negative association with E. turcicum, supporting its known potential as biological control agent against fungal pathogens. Our findings show that metagenomic sequencing is a useful tool for analysing the role of environmental factors and potential pathogen-microbiome interactions in shaping pathogen dynamics and evolution, suggesting their potential for effective pathogen management strategies.

Integrating eDNA metabarcoding and citizen science enhances avian ecological research.

Research Highlight: Bourbour etal., 2024. Feeding en route: Prey availability and traits influence prey selection by an avian predator on migration. Journal of Animal Ecology. Diet selection plays a key role in the eco-evolutionary dynamics of animals, exhibiting substantial variability across species, environments and seasons. The complex interplay between movement capability, hunting strategies, habitat use, prey traits and availability shapes the foraging patterns of avian predators. However, detailed information on how these birds exploit their extensive territories remains limited. In this study, Bourbour etal. utilised a novel integration of eDNA metabarcoding and citizen science to explore predator-prey interactions between migrating sharp-shinned hawks (Accipiter striatus) and an ephemeral avian prey community along North America's Pacific flyway. The research identified 1396 detections from the diet (65 species) of 588 migrating sharp-shinned hawks. Hawks' diet composition correlated with prey abundance indices sourced from the eBird database throughout the migration season, highlighting the significant impact of prey availability-shaped by migration tendency, flocking behaviour, and habitat-on raptor-songbird interactions. Notably, the study also found significant differences in prey size between male and female hawks, indicating that sexual dimorphism has led to diverse foraging strategies during migration. These findings underscore the potential of combining eDNA metabarcoding with citizen science to deepen our understanding of the foraging ecology of highly mobile and wide-ranging birds, as well as to monitor complex and vast ecosystems.

On the road to genomically defining bacterial intra-species units.

Over almost three decades, average nucleotide identity (ANI) analysis has been instrumental in operationally defining species in bacteria. However, barely any attention has been paid to soundly defining intra-species units employing ANI analyses until recently. Notably, some very recent publications are good steps forward in that direction. The level of granularity provided by these intra-species units will be relevant to understanding the eco-evolutionary dynamics and transmission of bacterial lineages and mobile genetic elements, antibiotic resistance, and virulence genes. These intra-species units will undoubtedly advance the genomic epidemiology of many bacterial pathogens. In the coming years, we anticipate that many studies will implement ANI-based definitions of different intra-species units, such as strains or sequence types, for many different bacterial species.

Rising above the noise: the influence of population dynamics on the evolution of acoustic signaling

Acoustic signaling is employed by many sexually reproducing species to select for mates and enhance fitness. However, signaling in dense populations can create an auditory background, or chorus, which may interfere with a signal receiver’s phonotactic selectivity, or the ability to distinguish individual signals. Feedback between the strength of an individual’s signal, phonotactic selectivity, and population size, may interact in complex ways to impact the evolution of signaling within a population, potentially leading to the emergence of silence. Here we formulate a general model that captures the dynamic feedback between individual acoustic signalers, phonotactic selectivity, and the population-level chorus to explore the eco-evolutionary dynamics of an acoustic trait within a population. We find that population dynamics have a significant influence on the evolutionary dynamics of the signaling trait, and that very sharp transitions separate conspicuous from silent populations. Our framework also reveals that increased phonotactic selectivity promotes the stability of signaling populations, and that transitions from signaling to silence are prone to hysteresis. We suggest that understanding the relationship between factors influencing population size, such as environmental productivity, as well as factors influencing phonotactic selectivity, such as anthropogenic noise, are central to understanding the complex mosaic of acoustically signaling and silent populations.

Attracting pollinators vs escaping herbivores: eco-evolutionary dynamics of plants confronted with an ecological trade-off

Attracting pollinators vs escaping herbivores: eco-evolutionary dynamics of plants confronted with an ecological trade-off

Genetic architecture of heritable leaf microbes.

Host-associated microbiomes are shaped by both their environment and host genetics, and often impact host performance. The scale of host genetic variation important to microbes is largely unknown yet fundamental to the community assembly of host-associated microbiomes, with implications for the eco-evolutionary dynamics of microbes and hosts. Using Ipomoea hederacea, ivyleaf morning glory, we generated matrilines differing in quantitative genetic variation and leaf shape, which is controlled by a single Mendelian locus. We then investigated the relative roles of Mendelian and quantitative genetic variation in structuring the leaf microbiome and how these two sources of genetic variation contributed to microbe heritability. We found that despite large effects of the environment, both Mendelian and quantitative genetic host variation contribute to microbe heritability and that the cumulative small effect genomic differences due to matriline explained as much or more microbial variation than a single large effect Mendelian locus. Furthermore, our results are the first to suggest that leaf shape itself contributes to variation in the abundances of some phyllosphere microbes.IMPORTANCEWe investigated how host genetic variation affects the assembly of Ipomoea hederacea's natural microbiome. We found that the genetic architecture of leaf-associated microbiomes involves both quantitative genetic variation and Mendelian traits, with similar contributions to microbe heritability. The existence of Mendelian and quantitative genetic variation for host-associated microbes means that plant evolution at the leaf shape locus or other quantitative genetic loci has the potential to shape microbial abundance and community composition.

Pushed waves, trailing edges, and extreme events: Eco-evolutionary dynamics of a geographic range shift in the owl limpet, Lottia gigantea.

As climatic variation re-shapes global biodiversity, understanding eco-evolutionary feedbacks during species range shifts is of increasing importance. Theory on range expansions distinguishes between two different forms: "pulled" and "pushed" waves. Pulled waves occur when the source of the expansion comes from low-density peripheral populations, while pushed waves occur when recruitment to the expanding edge is supplied by high-density populations closer to the species' core. How extreme events shape pushed/pulled wave expansion events, as well as trailing-edge declines/contractions, remains largely unexplored. We examined eco-evolutionary responses of a marine invertebrate (the owl limpet, Lottia gigantea) that increased in abundance during the 2014-2016 marine heatwaves near the poleward edge of its geographic range in the northeastern Pacific. We used whole-genome sequencing from 19 populations across >11 degrees of latitude to characterize genomic variation, gene flow, and demographic histories across the species' range. We estimated present-day dispersal potential and past climatic stability to identify how contemporary and historical seascape features shape genomic characteristics. Consistent with expectations of a pushed wave, we found little genomic differentiation between core and leading-edge populations, and higher genomic diversity at range edges. A large and well-mixed population in the northern edge of the species' range is likely a result of ocean current anomalies increasing larval settlement and high-dispersal potential across biogeographic boundaries. Trailing-edge populations have higher differentiation from core populations, possibly driven by local selection and limited gene flow, as well as high genomic diversity likely as a result of climatic stability during the Last Glacial Maximum. Our findings suggest that extreme events can drive poleward range expansions that carry the adaptive potential of core populations, while also cautioning that trailing-edge extirpations may threaten unique evolutionary variation. This work highlights the importance of understanding how both trailing and leading edges respond to global change and extreme events.

Competitive interactions affect introgression and population viability amidst maladaptive hybridization.

The deliberate release of captive-bred individuals, the accidental escape of domesticated strains, or the invasion of closely related conspecifics into wild populations can all lead to introgressive hybridization, which poses a challenge for conservation and wildlife management. Rates of introgression and the magnitude of associated demographic impacts vary widely across ecological contexts. However, the reasons for this variation remain poorly understood. One rarely considered phenomenon in this context is soft selection, wherein relative trait values determine success in intraspecific competition for a limiting resource. Here we develop an eco-genetic model explicitly focussed on understanding the influence of such competitive interactions on the eco-evolutionary dynamics of wild populations experiencing an influx of foreign/domesticated individuals. The model is applicable to any taxon that experiences natural or human-mediated inputs of locally maladapted genotypes ('intrusion'), in addition to phenotype-dependent competition for a limiting resource (e.g. breeding sites, feeding territories). The effects of both acute and chronic intrusion depended strongly on the relative competitiveness of intruders versus locals. When intruders were competitively inferior, density-dependent regulation limited their reproductive success (ability to compete for limited spawning sites), which prevented strong introgression or population declines from occurring. In contrast, when intruders were competitively superior, this amplified introgression and led to increased maladaptation of the admixed population. This had negative consequences for population size and population viability. The results were sensitive to the intrusion level, the magnitude of reproductive excess, trait heritability and the extent to which intruders were maladapted relative to locals. Our findings draw attention to under-appreciated interactions between phenotype-dependent competitive interactions and maladaptive hybridization, which may be critical to determining the impact captive breeding programmes and domesticated escapees can have on otherwise self-sustaining wild populations.

Pseudogenization of NK3 homeobox 2 (Nkx3.2) in monotremes provides insight into unique gastric anatomy and physiology.

The enzymatic breakdown and regulation of food passage through the vertebrate antral stomach and pyloric sphincter (antropyloric region) is a trait conserved over 450 million years. Development of the structures involved is underpinned by a highly conserved signalling pathway involving the hedgehog, bone morphogenetic protein and Wingless/Int-1 (Wnt) protein families. Monotremes are one of the few vertebrate lineages where acid-based digestion has been lost, and this is consistent with the lack of genes for hydrochloric acid secretion and gastric enzymes in the genomes of the platypus (Ornithorhynchus anatinus) and short-beaked echidna (Tachyglossus aculeatus) . Furthermore, these species feature unique gastric phenotypes, both with truncated and aglandular antral stomachs and the platypus with no pylorus. Here, we explore the genetic underpinning of monotreme gastric phenotypes, investigating genes important in antropyloric development using the newest monotreme genomes (mOrnAna1.pri.v4 and mTacAcu1) together with RNA-seq data. We found that the pathway constituents are generally conserved, but surprisingly, NK3 homeobox 2 (Nkx3.2) was pseudogenized in both platypus and echidna. We speculate that the unique sequence evolution of Grem1 and Bmp4 sequences in the echidna lineage may correlate with their pyloric-like restriction and that the convergent loss of gastric acid and stomach size genotypes and phenotypes in teleost and monotreme lineages may be a result of eco-evolutionary dynamics. These findings reflect the effects of gene loss on phenotypic evolution and further elucidate the genetic control of monotreme stomach anatomy and physiology.

Destabilized host-parasite dynamics in newly founded populations.

When species disperse into previously unoccupied habitats, new populations encounter unfamiliar species interactions such as altered parasite loads. Theory predicts that newly founded populations should exhibit destabilized eco-evolutionary fluctuations in infection rates and immune traits. However, to understand founder effects biologists typically rely on retrospective studies of range expansions, missing early-generation infection dynamics. To remedy this, we experimentally founded whole-lake populations of threespine stickleback. Infection rates were temporally stable in native source lakes. In contrast, newly founded populations exhibit destabilized host-parasite dynamics: high starting infection rates led to increases in a heritable immune trait (peritoneal fibrosis), suppressing infection rates. The resulting temporal auto-correlation between infection and immunity suggest that newly founded populations can exhibit rapid host-parasite eco-evolutionary dynamics.

Species interactions and eco-evolutionary dynamics of dispersal: the diversity dependence of dispersal.

Dispersal plays a pivotal role in the eco-evolutionary dynamics of spatially structured populations, communities and ecosystems. As an individual-based trait, dispersal is subject to both plasticity and evolution. Its dependence on conditions and context is well understood within single-species metapopulations. However, species do not exist in isolation; they interact locally through various horizontal and vertical interactions. While the significance of species interactions is recognized for species coexistence and food web functioning, our understanding of their influence on regional dynamics, such as their impact on spatial dynamics in metacommunities and meta-food webs, remains limited. Building upon insights from behavioural and community ecology, we aim to elucidate biodiversity as both a driver and an outcome of connectivity. By synthesizing conceptual, theoretical and empirical contributions from global experts in the field, we seek to explore how a more mechanistic understanding of diversity-dispersal relationships influences the distribution of species in spatially and temporally changing environments. Our findings highlight the importance of explicitly considering interspecific interactions as drivers of dispersal, thus reshaping our understanding of fundamental dynamics including species coexistence and the emergent dynamics of metacommunities and meta-ecosystems. We envision that this initiative will pave the way for advanced forecasting approaches to understanding biodiversity dynamics under the pressures of global change. This article is part of the theme issue 'Diversity-dependence of dispersal: interspecific interactions determine spatial dynamics'.

Eco-evolutionary dynamics of cooperative antimicrobial resistance in a population of fluctuating volume and size

Antimicrobial resistance to drugs (AMR), a global threat to human and animal health, is often regarded as resulting from cooperative behaviour. Moreover, microbes generally evolve in volatile environments that, together with demographic fluctuations (birth and death events), drastically alter population size and strain survival. Motivated by the need to better understand the evolution of AMR, we study a population of time-varying size consisting of two competing strains, one drug-resistant and one drug-sensitive, subject to demographic and environmental variability. This is modelled by a binary carrying capacity randomly switching between mild and harsh environmental conditions, and driving the fluctuating volume (total amount of nutrients and antimicrobials at fixed concentration), and thus the size of the community (number of resistant and sensitive cells). We assume that AMR is a shared public good when the concentration of resistant cells exceeds a fixed concentration cooperation threshold, above which the sensitive strain has a growth advantage, whereas resistant cells dominate below it. Using computational means, and devising an analytical treatment (built on suitable quenched and annealed averaging procedures), we fully characterise the influence of fluctuations on the eco-evolutionary dynamics of AMR, and notably obtain specific strain fixation and long-lasting coexistence probabilities as a function of the environmental variation rate and cooperation threshold. We find that microbial strains tend to coexistence, but demographic fluctuations eventually lead to the extinction of resistant or sensitive cells for small or large values of the concentration cooperation threshold, respectively. This also holds for dynamic environments, whose specific properties determine the extinction timescale.