Eco-evolutionary spatial dynamics in the Glanville fritillary butterfly

Demography and environment modulate the effects of genetic diversity on extinction risk in a butterfly metapopulation

Linking genetic diversity to extinction is a common goal in genomic studies. Recently, a debate has arisen regarding the importance of genetic variation in conservation as some studies have failed to find associations between genome-wide genetic diversity and extinction risk. However, only rarely are genetic diversity and fitness measured together in the wild, and typically demographic history and environment are ignored. It is therefore difficult to infer whether a lack of an association is real or obscured by confounding factors. To address these shortcomings, we analyzed genetic data from 7,501 individuals with extinction data from 279 meadows and mortality of 1,742 larval nests in a butterfly metapopulation. We found a strong negative association between genetic diversity and extinction when considering only heterozygosity in models. However, this association disappeared when accounting for ecological covariates, suggesting a confounding between demography and genetics and a more complex role for heterozygosity in extinction risk. Modeling interactions between heterozygosity and demographic variables revealed that associations between extinction and heterozygosity were context-dependent. For example, extinction declined with increasing heterozygosity in large, but not currently small populations, although negative associations between heterozygosity, extinction, and mortality were detected in small populations with a recent history of decline. We conclude that low genetic diversity is an important predictor of extinction, predicting >25% increase in extinction beyond ecological factors in certain contexts. These results highlight that inferences about the importance of genetic diversity for population viability should not rely on genomic data alone but require investments in obtaining demographic and environmental data from natural populations.

Evolutionary Game Dynamics with Environmental Feedback in a Network with Two Communities

Recent developments of eco-evolutionary models have shown that evolving feedbacks between behavioral strategies and the environment of game interactions, leading to changes in the underlying payoff matrix, can impact the underlying population dynamics in various manners. We propose and analyze an eco-evolutionary game dynamics model on a network with two communities such that players interact with other players in the same community and those in the opposite community at different rates. In our model, we consider two-person matrix games with pairwise interactions occurring on individual edges and assume that the environmental state depends on edges rather than on nodes or being globally shared in the population. We analytically determine the equilibria and their stability under a symmetric population structure assumption, and we also numerically study the replicator dynamics of the general model. The model shows rich dynamical behavior, such as multiple transcritical bifurcations, multistability, and anti-synchronous oscillations. Our work offers insights into understanding how the presence of community structure impacts the eco-evolutionary dynamics within and between niches.

Genetic architecture of dispersal behaviour in the post-harvest pest and model organism Tribolium castaneum

Dispersal behaviour is an important aspect of the life-history of animals. However, the genetic architecture of dispersal-related traits is often obscure or unknown, even in well studied species. Tribolium castaneum is a globally significant post-harvest pest and established model organism, yet studies of its dispersal have shown ambiguous results and the genetic basis of this behaviour remains unresolved. We combine experimental evolution and agent-based modelling to investigate the number of loci underlying dispersal in T. castaneum, and whether the trait is sex-linked. Our findings demonstrate rapid evolution of dispersal behaviour under selection. We find no evidence of sex-biases in the dispersal behaviour of the offspring of crosses, supporting an autosomal genetic basis of the trait. Moreover, simulated data approximates experimental data under simulated scenarios where the dispersal trait is controlled by one or few loci, but not many loci. Levels of dispersal in experimentally inbred lines, compared with simulations, indicate that a single locus model is not well supported. Taken together, these lines of evidence support an oligogenic architecture underlying dispersal in Tribolium castaneum. These results have implications for applied pest management and for our understanding of the evolution of dispersal in the coleoptera, the world's most species-rich order.

Spatial metapopulation dynamics with local and global colonization

We revisit a spatial metapopulation model on continuous space as a stochastic point pattern dynamics. In the model, local patches as points are distributed with a certain spatial configuration and status of each patch changes stochastically between occupied and empty: an occupied patch becomes empty by local extinction and an empty patch becomes occupied both by local and global colonization. We carry out simulation analysis and derive an analytical model in terms of singlet, pair and triplet probabilities that describe the stochastic dynamics. Using a simple closure that approximates triplet probabilities by singlet and pair probabilities, we show that equilibrium singlet and pair probabilities can be analytically derived. The derived equilibrium properties successfully describe simulation results under a certain condition where the range of local colonization and the proportion of global colonization play key roles. Our model is an extension of the classical non-spatial Levins model to a spatially explicit metapopulation model. We appeal the advantage of point pattern approach to study spatial dynamics in general ecology and call for the need to deepen our understanding of mathematical tools to explore point pattern dynamics.

Different mechanisms underpin the decline in growth of anchovies and sardines of the Bay of Biscay

Declines in individuals' growth in exploited fish species are generally attributed to evolutionary consequences of size-selective fishing or to plastic responses due to constraints set by changing environmental conditions dampening individuals' growth. However, other processes such as growth compensation and non-directional selection can occur and their importance on the overall phenotypic response of exploited populations has largely been ignored. Using otolith growth data collected in European anchovy and sardine of the Bay of Biscay (18 cohorts from 2000 to 2018), we parameterized the breeder's equation to determine whether declines in size-at-age in these species were due to an adaptive response (i.e. related to directional or non-directional selection differentials within parental cohorts) or a plastic response (i.e. related to changes in environmental). We found that growth at age-0 in anchovy declined between parents and their offspring when biomass increased and the selective disappearance of large individuals was high in parents. Therefore, an adaptive response probably occurred in years with high fishing effort and the large increase in biomass after the collapse of this stock maintained this adaptive response subsequently. In sardine offspring, higher growth at age-0 was associated with increasing biomass between parents and offspring, suggesting a plastic response to a bottom-up process (i.e. a change in food quantity or quality). Parental cohorts in which selection favoured individuals with high growth compensation produced offspring high catch up growth rates, which may explain the smaller decline in growth in sardine relative to anchovy. Finally, on non-directional selection differentials were not significantly related to the changes in growth at age-0 and growth compensation at age-1 in both species. Although anchovy and sardine have similar ecologies, the mechanisms underlying the declines in their growth are clearly different. The consequences of the exploitation of natural populations could be long lasting if density-dependent processes follow adaptive changes.

Landscape Structure Affects Metapopulation-Scale Tipping Points

AbstractEven when environments deteriorate gradually, ecosystems may shift abruptly from one state to another. Such catastrophic shifts are difficult to predict and sometimes to reverse (so-called hysteresis). While well studied in simplified contexts, we lack a general understanding of how catastrophic shifts spread in realistically spatially structured landscapes. For different types of landscape structures, including typical terrestrial modular and riverine dendritic networks, we here investigate landscape-scale stability in metapopulations whose patches can locally exhibit catastrophic shifts. We find that such metapopulations usually exhibit large-scale catastrophic shifts and hysteresis and that the properties of these shifts depend strongly on the metapopulation spatial structure and on the population dispersal rate: an intermediate dispersal rate, a low average degree, or a riverine spatial structure can largely reduce hysteresis size. Our study suggests that large-scale restoration is easier with spatially clustered restoration efforts and in populations characterized by an intermediate dispersal rate.

Plant disease risk is modified by multiple global change drivers

Plant diseases are strongly influenced by host biodiversity, spatial structure, and abiotic conditions. All of these are undergoing rapid change, as the climate is warming, habitats are being lost, and nitrogen deposition is changing nutrient dynamics of ecosystems with ensuing consequences for biodiversity. Here, I review examples of plant-pathogen associations to demonstrate how our ability to understand, model and predict disease dynamics is becoming increasingly difficult, as both plant and pathogen populations and communities are undergoing extensive change. The extent of this change is influenced via both direct and combined effects of global change drivers, and especially the latter are still poorly understood. Change at one trophic level is expected to drive change also at the other, and hence feedback loops between plants and their pathogens are expected to drive changes in disease risk both through ecological as well as evolutionary mechanisms. Many of the examples discussed here demonstrate an increase in disease risk as a result of ongoing change, suggesting that unless we successfully mitigate global environmental change, plant disease is going to become an increasingly heavy burden on our societies with far-reaching consequences for food security and functioning of ecosystems.

Movements of a Specialist Butterfly in Relation to Mowing Management of Its Habitat Patches

Over the centuries, mowing and grazing have been crucial for sustaining populations of grassland insects and their overall diversity in Europe. While long-term positive effects of mowing have been studied in more detail, little is known about the direct impacts of mowing on adult butterflies. Here, we explore how different habitat management (mown, recovered after mowing and unmown) affects movements and population estimates of the endangered specialist butterfly Phengaris teleius. The results showed higher dispersal probability from mown (22%) and recovered meadows (16%) than from the unmown ones (9%). However, mowing shortened the average dispersal distances (mown = 102 m, recovered = 198 m, unmown = 246 m) and reduced butterfly population size. In contrast, a larger area of the habitat patches promoted long-distance dispersal and sustained larger populations. We hypothesise that mowing caused depletion of resources and triggered dispersal of poorly adapted individuals. This behaviour is maladaptive and could lead to higher dispersal-related mortality; thus, mowing should be avoided before and during the butterfly flight period. This study suggests that the species’ persistence in a fragmented landscape depends on large, unmown and interconnected habitats that support more viable populations, promote long-distance dispersal, and enable (re)colonisation of vacant patches.

Self-segregation in heterogeneous metapopulation landscapes

Complex interactions are at the root of the population dynamics of many natural systems, particularly for being responsible for the allocation of species and individuals across apposite niches of the ecological landscapes. On the other side, the randomness that unavoidably characterises complex systems has increasingly challenged the niche paradigm providing alternative neutral theoretical models. We introduce a network-inspired metapopulation individual-based model (IBM), hereby named self-segregation, where the density of individuals in the hosting patches (local habitats) drives the individuals spatial assembling while still constrained by nodes' saturation. In particular, we prove that the core-periphery structure of the networked landscape triggers the spontaneous emergence of vacant habitat patches, which segregate the population in multistable patterns of isolated (sub)communities separated by empty patches. Furthermore, a quantisation effect in the number of vacant patches is observed once the total system mass varies continuously, emphasising thus a striking feature of the robustness of population stationary distributions. Notably, our model reproduces the patch vacancy found in the fragmented habitat of the Glanville fritillary butterfly Melitaea cinxia, an endemic species of the Åland islands. We argue that such spontaneous breaking of the natural habitat supports the concept of the highly contentious (Grinnellian) niche vacancy and also suggests a new mechanism for the endogeneous habitat fragmentation and consequently the peripatric speciation.

Limited predictability of body length in a fish population

Recent theoretical studies have identified chaotic dynamics in eco-evolutionary models. Yet, empirical evidence for eco-evolutionary chaos in natural ecosystems is lacking. In this study, we combine analyses of empirical data and an eco-evolutionary model to uncover chaotic dynamics of body length in a fish population (northeast Arctic cod: Gadus morhua). Consistent with chaotic attractors, the largest Lyapunov exponent (LE) of empirical data is positive, and approximately matches the LE of the model calculation, thus suggesting the potential for chaotic dynamics in this fish population. We also find that the autocorrelation function (ACF) of both empirical data and eco-evolutionary model shows a similar lag of approximately 7 years. Our combined analyses of natural time series and mathematical models suggest that chaotic dynamics of a phenotypic trait may be driven by trait evolution. This finding supports a growing theory that eco-evolutionary feedbacks can produce chaotic dynamics.

The importance of eco-evolutionary dynamics for predicting and managing insect range shifts

Evolutionary change impacts the rate at which insect pests, pollinators, or disease vectors expand or contract their geographic ranges. Although evolutionary changes, and their ecological feedbacks, strongly affect these risks and associated ecological and economic consequences, they are often underappreciated in management efforts. Greater rigor and scope in study design, coupled with innovative technologies and approaches, facilitates our understanding of the causes and consequences of eco-evolutionary dynamics in insect range shifts. Future efforts need to ensure that forecasts allow for demographic and evolutionary change and that management strategies will maximize (or minimize) the adaptive potential of range-shifting insects, with benefits for biodiversity and ecosystem services.

Multi-surveyor capture-mark-recapture as a powerful tool for butterfly population monitoring in the pre-imaginal stage

For many elusive insect species, which are difficult to cover by standard monitoring schemes, innovative survey methods are needed to gain robust data on abundance and population trends. We suggest a monitoring of overwintering larvae for the endangered nymphalid butterfly Limenitis reducta. We tested different removal and capture-mark-recapture (CMR) approaches in a field study in the "Alb-Donau" region, Germany. Classical removal and CMR studies require movement of the organisms under study, but in our approach, we replaced movement of the study organisms by random movement of multiple different surveyors. We tested the validity of the approach by comparing detection frequencies from our field data with simulated detections. Our results indicate that multi-surveyor removal/CMR techniques are suitable for estimating abundance of overwintering L. reducta larvae. Depending on surveyor experience, the average detection probability ranged between 16% for novices and 35% for experts. The uncertainty of population estimates increased with a decrease in personnel expenditure. Estimated larval densities on a spruce clear-cut varied between one and three individuals per 100 m2, probably related to habitat conditions. We suggest a CMR approach with three to four trained surveyors for the monitoring of L. reducta populations in the overwintering stage. Compared with previous sampling methods, our approach is a powerful tool with clear advantages: long survey period, estimates of the absolute population size accompanied by uncertainty measures, and estimates of overwinter mortality. The proposed method can be adapted and used for several different butterfly species, other insect taxa with specific immobile life stages, and some sessile organisms, for example, elusive plants, fungi, or corals.

Climate-driven variation in biotic interactions provides a narrow and variable window of opportunity for an insect herbivore at its ecological margin

Climate-driven geographic range shifts have been associated with transitions between dietary specialism and generalism at range margins. The mechanisms underpinning these often transient niche breadth modifications are poorly known, but utilization of novel resources likely depends on phenological synchrony between the consumer and resource. We use a climate-driven range and host shift by the butterfly Aricia agestis to test how climate-driven changes in host phenology and condition affect phenological synchrony, and consider implications for host use. Our data suggest that the perennial plant that was the primary host before range expansion is a more reliable resource than the annual Geraniaceae upon which the butterfly has become specialized in newly colonized parts of its range. In particular, climate-driven phenological variation in the novel host Geranium dissectum generates a narrow and variable 'window of opportunity' for larval productivity in summer. Therefore, although climatic change may allow species to shift hosts and colonise novel environments, specialization on phenologically limited hosts may not persist at ecological margins as climate change continues. We highlight the potential role for phenological (a)synchrony in determining lability of consumer-resource associations at range margins and the importance of considering causes of synchrony in biotic interactions when predicting range shifts. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

The genome sequence of the Glanville fritillary, Melitaea cinxia (Linnaeus, 1758)

We present a genome assembly from an individual male Melitaea cinxia (the Glanville fritillary; Arthropoda; Insecta; Lepidoptera; Nymphalidae). The genome sequence is 499 megabases in span. The complete assembly is scaffolded into 31 chromosomal pseudomolecules, with the Z sex chromosome assembled. Gene annotation of this assembly on Ensembl has identified 13,666 protein coding genes.

Rescaling Biology: Increasing Integration Across Biological Scales and Subdisciplines to Enhance Understanding and Prediction

The organization of the living world covers a vast range of spatiotemporal scales, from molecules to the biosphere, seconds to centuries. Biologists working within specialized subdisciplines tend to focus on different ranges of scales. Therefore, developing frameworks that enable testing questions and predictions of scaling require sufficient understanding of complex processes across biological subdisciplines and spatiotemporal scales. Frameworks that enable scaling across subdisciplines would ideally allow us to test hypotheses about the degree to which explicit integration across spatiotemporal scales is needed for predicting the outcome of biological processes. For instance, how does genomic variation within populations allow us to explain community structure? How do the dynamics of cellular metabolism translate to our understanding of whole-ecosystem metabolism? Do patterns and processes operate seamlessly across biological scales, or are there fundamental laws of biological scaling that limit our ability to make predictions from one scale to another? Similarly, can sub-organismal structures and processes be sufficiently understood in isolation of potential feedbacks from the population, community, or ecosystem levels? And can we infer the sub-organismal processes from data on the population, community, or ecosystem scale? Concerted efforts to develop more cross-disciplinary frameworks will open doors to a more fully integrated field of biology. In this paper we discuss how we might integrate across scales, specifically by 1. Identifying scales and boundaries, 2. Determining analogous units and processes across scales, 3. Developing frameworks to unite multiple scales, and 4. Extending frameworks to new empirical systems.

Long-term spatiotemporal genetic structure of an accidental parasitoid introduction, and local changes in prevalence of its associated Wolbachia symbiont

Population bottlenecks associated with founder events strongly impact the establishment and genetic makeup of populations. In addition to their genotype, founding individuals also bring along parasites, as well as symbionts that can manipulate the phenotype of their host, affecting the host population establishment, dynamics and evolution. Thus, to understand introduction, invasion, and spread, we should identify the roles played by accompanying symbionts. In 1991, the parasitoid wasp, Hyposoter horticola, and its associated hyperparasitoid were accidentally introduced from the main Åland islands, Finland, to an isolated island in the archipelago, along with their host, the Glanville fritillary butterfly. Though the receiving island was unoccupied, the butterfly was present on some of the small islands in the vicinity. The three introduced species have persisted locally ever since. A strain of the endosymbiotic bacterium Wolbachia has an intermediate prevalence in the parasitoid H. horticola across the main Åland population. The infection increases its susceptibility of to hyperparasitism. We investigated the establishment and spread of the parasitoid, along with patterns of prevalence of its symbiont using 323 specimens collected between 1992 and 2013, from five localities across Åland, including the source and introduced populations. Using 14 microsatellites and one mitochondrial marker, we suggest that the relatively diverse founding population and occasional migration between islands might have facilitated the persistence of all isolated populations, despite multiple local population crashes. We also show that where the hyperparasitoid is absent, and thus selection against infected wasp genotypes is relaxed, there is near-fixation of Wolbachia.

Metapopulation Patterns of Iberian Butterflies Revealed by Fuzzy Logic

Metapopulation theory considers that the populations of many species are fragmented into patches connected by the migration of individuals through an interterritorial matrix. We applied fuzzy set theory and environmental favorability (F) functions to reveal the metapopulational structure of the 222 butterfly species in the Iberian Peninsula. We used the sets of contiguous grid cells with high favorability (F ≥ 0.8), to identify the favorable patches for each species. We superimposed the known occurrence data to reveal the occupied and empty favorable patches, as unoccupied patches are functional in a metapopulation dynamics analysis. We analyzed the connectivity between patches of each metapopulation by focusing on the territory of intermediate and low favorability for the species (F < 0.8). The friction that each cell opposes to the passage of individuals was computed as 1-F. We used the r.cost function of QGIS to calculate the cost of reaching each cell from a favorable patch. The inverse of the cost was computed as connectivity. Only 126 species can be considered to have a metapopulation structure. These metapopulation structures are part of the dark biodiversity of butterflies because their identification is not evident from the observation of the occurrence data but was revealed using favorability functions.

Multidimensional plasticity in the Glanville fritillary butterfly: larval performance is temperature, host and family specific

Variation in environmental conditions during development can lead to changes in life-history traits with long-lasting effects. Here, we study how variation in temperature and host plant (i.e. the consequences of potential maternal oviposition choices) affects a suite of life-history traits in pre-diapause larvae of the Glanville fritillary butterfly. We focus on offspring survival, larval growth rates and relative fat reserves, and pay specific attention to intraspecific variation in the responses (G × E × E). Globally, thermal performance and survival curves varied between diets of two host plants, suggesting that host modifies the temperature impact, or vice versa. Additionally, we show that the relative fat content has a host-dependent, discontinuous response to developmental temperature. This implies that a potential switch in resource allocation, from more investment in growth at lower temperatures to storage at higher temperatures, is dependent on the larval diet. Interestingly, a large proportion of the variance in larval performance is explained by differences among families, or interactions with this variable. Finally, we demonstrate that these family-specific responses to the host plant remain largely consistent across thermal environments. Together, the results of our study underscore the importance of paying attention to intraspecific trait variation in the field of evolutionary ecology.

Alleles in metabolic and oxygen-sensing genes are associated with antagonistic pleiotropic effects on life history traits and population fitness in an ecological model insect

Genes with opposing effects on fitness at different life stages are the mechanistic basis for evolutionary theories of aging and life history. Examples come from studies of mutations in model organisms, but there is little knowledge of genetic bases of life history tradeoffs in natural populations. Here, we test the hypothesis that alleles affecting oxygen sensing in Glanville fritillary butterflies have opposing effects on larval versus adult fitness-related traits. Intermediate-frequency alleles in Succinate dehydrogenase d, and to a lesser extent Hypoxia inducible factor 1α, are associated in larvae with variation in metabolic rate and activation of the hypoxia inducible factor (HIF) pathway, which affects tracheal development and delivery of oxygen to adult flight muscles. A dominant Sdhd allele is likely to cause antagonistic pleiotropy for fitness through its opposing effects on larval metabolic and growth rate versus adult flight and dispersal, and may have additional effects arising from sensitivity to low-iron host plants. Prior results in Glanville fritillaries indicate that fitness of alleles in Sdhd and another antagonistically pleiotropic metabolic gene, Phosphoglucose isomerase, depend strongly on the size and distribution of host plant patches. Hence, these intermediate-frequency alleles are involved in ecoevolutionary dynamics involving life history tradeoffs.

Host plant defense produces species-specific alterations to flight muscle protein structure and flight-related fitness traits of two armyworms

Insects manifest phenotypic plasticity in their development and behavior in response to plant defenses, via molecular mechanisms that produce tissue-specific changes. Phenotypic changes might vary between species that differ in their preferred hosts and these effects could extend beyond larval stages. To test this, we manipulated the diet of southern armyworm (SAW; Spodoptera eridania) and fall armyworm (FAW; Spodoptera frugiperda) using a tomato mutant for jasmonic acid plant defense pathway (def1), and wild-type plants, and then quantified gene expression of Troponin t (Tnt) and flight muscle metabolism of the adult insects. Differences in Tnt spliceform ratios in insect flight muscles correlate with changes to flight muscle metabolism and flight muscle output. We found that SAW adults reared on induced def1 plants had a higher relative abundance (RA) of the A isoform of Troponin t (Tnt A) in their flight muscles; in contrast, FAW adults reared on induced def1 plants had a lower RA of Tnt A in their flight muscles compared with adults reared on def1 and controls. Although mass-adjusted flight metabolic rate showed no independent host plant effects in either species, higher flight metabolic rates in SAW correlated with increased RA of Tnt A Flight muscle metabolism also showed an interaction of host plants with Tnt A in both species, suggesting that host plants might be influencing flight muscle metabolic output by altering Tnt This study illustrates how insects respond to variation in host plant chemical defense by phenotypic modifications to their flight muscle proteins, with possible implications for dispersal.

The gut bacterial community affects immunity but not metabolism in a specialist herbivorous butterfly

Plant tissues often lack essential nutritive elements and may contain a range of secondary toxic compounds. As nutritional imbalance in food intake may affect the performances of herbivores, the latter have evolved a variety of physiological mechanisms to cope with the challenges of digesting their plant-based diet. Some of these strategies involve living in association with symbiotic microbes that promote the digestion and detoxification of plant compounds or supply their host with essential nutrients missing from the plant diet. In Lepidoptera, a growing body of evidence has, however, recently challenged the idea that herbivores are nutritionally dependent on their gut microbial community. It is suggested that many of the herbivorous Lepidopteran species may not host a resident microbial community, but rather a transient one, acquired from their environment and diet. Studies directly testing these hypotheses are however scarce and come from an even more limited number of species.By coupling comparative metabarcoding, immune gene expression, and metabolomics analyses with experimental manipulation of the gut microbial community of prediapause larvae of the Glanville fritillary butterfly (Melitaea cinxia, L.), we tested whether the gut microbial community supports early larval growth and survival, or modulates metabolism or immunity during early stages of development.We successfully altered this microbiota through antibiotic treatments and consecutively restored it through fecal transplants from conspecifics. Our study suggests that although the microbiota is involved in the up-regulation of an antimicrobial peptide, it did not affect the life history traits or the metabolism of early instars larvae.This study confirms the poor impact of the microbiota on diverse life history traits of yet another Lepidoptera species. However, it also suggests that potential eco-evolutionary host-symbiont strategies that take place in the gut of herbivorous butterfly hosts might have been disregarded, particularly how the microbiota may affect the host immune system homeostasis.

Evolutionary origins for ecological patterns in space

Historically, many biologists assumed that evolution and ecology acted independently because evolution occurred over distances too great to influence most ecological patterns. Today, evidence indicates that evolution can operate over a range of spatial scales, including fine spatial scales. Thus, evolutionary divergence across space might frequently interact with the mechanisms that also determine spatial ecological patterns. Here, we synthesize insights from 500 eco-evolutionary studies and develop a predictive framework that seeks to understand whether and when evolution amplifies, dampens, or creates ecological patterns. We demonstrate that local adaptation can alter everything from spatial variation in population abundances to ecosystem properties. We uncover 14 mechanisms that can mediate the outcome of evolution on spatial ecological patterns. Sometimes, evolution amplifies environmental variation, especially when selection enhances resource uptake or patch selection. The local evolution of foundation or keystone species can create ecological patterns where none existed originally. However, most often, we find that evolution dampens existing environmental gradients, because local adaptation evens out fitness across environments and thus counteracts the variation in associated ecological patterns. Consequently, evolution generally smooths out the underlying heterogeneity in nature, making the world appear less ragged than it would be in the absence of evolution. We end by highlighting the future research needed to inform a fully integrated and predictive biology that accounts for eco-evolutionary interactions in both space and time.

Steering eco-evolutionary game dynamics with manifold control

Feedback loops between population dynamics of individuals and their ecological environment are ubiquitously found in nature and have shown profound effects on the resulting eco-evolutionary dynamics. By incorporating linear environmental feedback law into the replicator dynamics of two-player games, recent theoretical studies have shed light on understanding the oscillating dynamics of the social dilemma. However, the detailed effects of more general nonlinear feedback loops in multi-player games, which are more common especially in microbial systems, remain unclear. Here, we focus on ecological public goods games with environmental feedbacks driven by a nonlinear selection gradient. Unlike previous models, multiple segments of stable and unstable equilibrium manifolds can emerge from the population dynamical systems. We find that a larger relative asymmetrical feedback speed for group interactions centred on cooperators not only accelerates the convergence of stable manifolds but also increases the attraction basin of these stable manifolds. Furthermore, our work offers an innovative manifold control approach: by designing appropriate switching control laws, we are able to steer the eco-evolutionary dynamics to any desired population state. Our mathematical framework is an important generalization and complement to coevolutionary game dynamics, and also fills the theoretical gap in guiding the widespread problem of population state control in microbial experiments.

Population size changes and extinction risk of populations driven by mutant interactors

Spontaneous random mutations are an important source of variation in populations. Many evolutionary models consider mutants with a fixed fitness, chosen from a fitness distribution without considering microscopic interactions among the residents and mutants. Here, we go beyond this and consider "mutant interactors," which lead to new interactions between the residents and invading mutants that can affect the population size and the extinction risk of populations. We model microscopic interactions between individuals by using a dynamic interaction matrix, the dimension of which increases with the emergence of a new mutant and decreases with extinction. The new interaction parameters of the mutant follow a probability distribution around the payoff entries of its ancestor. These new interactions can drive the population away from the previous equilibrium and lead to changes in the population size. Thus, the population size is an evolving property rather than an externally controlled variable. We calculate the average population size of our stochastic system over time and quantify the extinction risk of the population by the mean time to extinction.

Dispersal syndromes can impact ecosystem functioning in spatially structured freshwater populations

Dispersal can strongly influence ecological and evolutionary dynamics. Besides the direct contribution of dispersal to population dynamics, dispersers often differ in their phenotypic attributes from non-dispersers, which leads to dispersal syndromes. The consequences of such dispersal syndromes have been widely explored at the population and community level; however, to date, ecosystem-level effects remain unclear. Here, we examine whether dispersing and resident individuals of two different aquatic keystone invertebrate species have different contributions to detrital processing, a key function in freshwater ecosystems. Using experimental two-patch systems, we found no difference in leaf consumption rates with dispersal status of the common native species Gammarus fossarum. In Dikerogammarus villosus, however, a Ponto-Caspian species now expanding throughout Europe, dispersers consumed leaf litter at roughly three times the rate of non-dispersers. Furthermore, this put the contribution of dispersing D. villosus to leaf litter processing on par with native G. fossarum, after adjusting for differences in organismal size. Given that leaf litter decomposition is a key function in aquatic ecosystems, and the rapid species turnover in freshwater habitats with range expansions of non-native species, this finding suggests that dispersal syndromes may have important consequences for ecosystem functioning.

Ecological feedback on diffusion dynamics

Spatial patterns are ubiquitous across different scales of organization in ecological systems. Animal coat pattern, spatial organization of insect colonies and vegetation in arid areas are prominent examples from such diverse ecologies. Typically, pattern formation has been described by reaction-diffusion equations, which consider individuals dispersing between subpopulations of a global pool. This framework applied to public goods game nicely showed the endurance of populations via diffusion and generation of spatial patterns. However, how the spatial characteristics, such as diffusion, are related to the eco-evolutionary process as well as the nature of the feedback from evolution to ecology and vice versa, has been so far neglected. We present a thorough analysis of the ecologically driven evolutionary dynamics in a spatially extended version of ecological public goods games. Furthermore, we show how these evolutionary dynamics feed back into shaping the ecology, thus together determining the fate of the system.

Frontiers in Metapopulation Biology: The Legacy of Ilkka Hanski

This review of metapopulation biology has a special focus on Professor Ilkka Hanski's (1953–2016) research. Hanski made seminal contributions to both empirical and theoretical metapopulation biology throughout his scientific career. Hanski's early research focused on ecological aspects of metapopulation biology, in particular how the spatial structure of a landscape influences extinction thresholds and how habitat loss and fragmentation can result in extinction debt. Hanski then used the Glanville fritillary system as a natural laboratory within which he studied genetic and evolutionary processes, such as the influence of inbreeding on extinction risk and variation in selection for dispersal traits generated by landscape variation. During the last years of his career, Hanski's work was in the forefront of the rapidly developing field of eco-evolutionary dynamics. Hanski was a pioneer in showing how molecular-level underpinnings of trait variation can explain why evolutionary change can occur rapidly in natural populations and how these changes can subsequently influence ecological dynamics.

Evolutionarily stable communities: a framework for understanding the role of trait evolution in the maintenance of diversity

Biological diversity depends on the interplay between evolutionary diversification and ecological mechanisms allowing species to coexist. Current research increasingly integrates ecology and evolution over a range of timescales, but our common conceptual framework for understanding species coexistence requires better incorporation of evolutionary processes. Here, we focus on the idea of evolutionarily stable communities (ESCs), which are theoretical endpoints of evolution in a community context. We use ESCs as a unifying framework to highlight some important but under-appreciated theoretical results, and we review empirical research relevant to these theoretical predictions. We explain how, in addition to generating diversity, evolution can also limit diversity by reducing the effectiveness of coexistence mechanisms. The coevolving traits of competing species may either diverge or converge, depending on whether the number of species in the community is low (undersaturated) or high (oversaturated) relative to the ESC. Competition in oversaturated communities can lead to extinction or neutrally coexisting, ecologically equivalent species. It is critical to consider trait evolution when investigating fundamental ecological questions like the strength of different coexistence mechanisms, the feasibility of ecologically equivalent species, and the interpretation of different patterns of trait dispersion.

Synchronization and entrainment of metapopulations: A trade-off among time-induced heterogeneity, dispersal, and seasonal force

Demographic and environmental heterogeneities are prevalent across many natural systems. Earlier studies on metapopulation models have mostly considered heterogeneities either in the demographic parameters or in the interaction strength and topology between the spatially separated patches. In contrast, here we study the dynamics of a metapopulation model where each of the uncoupled patches has different periods of oscillations (period mismatch). We show different synchronization dynamics governed by both period mismatch and dispersal in neighboring patches. Indeed, we find both appearance and disappearance of phase synchronization, quasiperiodic oscillations, and period doubling of limit cycle. We also quantify the effect of seasonal variation (entrainment) and dispersal on species synchrony using phase-response curve and a synchrony measure, which thereof identify the influence of stochasticity on species persistence through trade-off mechanisms. Our results show that trade-offs among period mismatch, dispersal, and external force can drive entrained oscillations as well as asynchronous population dynamics that structure ecological communities.

A general method for estimating seed dormancy and colonisation in annual plants from the observation of existing flora

In plant ecology, characterising colonisation and extinction in plant metapopulations is challenging due to the non-detectable seed bank that allows plants to emerge after several years of absence. In this study, we used a Hidden Markov Model to characterise seed dormancy, colonisation and germination solely from the presence-absence of standing flora. Applying the model to data from a long-term survey of 38 annual weeds across France, we identified three homogeneous functional groups: (1) species persisting preferentially through spatial colonisation, (2) species persisting preferentially through seed dormancy and (3) a mix of both strategies. These groups are consistent with existing ecological knowledge, demonstrating that ecologically meaningful parameters can be estimated from simple presence-absence observations. These results indicate that such studies could contribute to the design of weed management strategies. They also open the possibility of testing life-history theories such as the dormancy/colonisation trade-off in natura.

Enzyme polymorphism, oxygen and injury: a lipidomic analysis of flight-induced oxidative damage in a succinate dehydrogenase d (Sdhd)-polymorphic insect

When active tissues receive insufficient oxygen to meet metabolic demand, succinate accumulates and has two fundamental effects: it causes ischemia-reperfusion injury while also activating the hypoxia-inducible factor pathway (HIF). The Glanville fritillary butterfly (Melitaea cinxia) possesses a balanced polymorphism in Sdhd, shown previously to affect HIF pathway activation and tracheal morphology and used here to experimentally test the hypothesis that variation in succinate dehydrogenase affects oxidative injury. We stimulated butterflies to fly continuously in a respirometer (3 min duration), which typically caused episodes of exhaustion and recovery, suggesting a potential for cellular injury from hypoxia and reoxygenation in flight muscles. Indeed, flight muscle from butterflies flown on consecutive days had lipidome profiles similar to those of rested paraquat-injected butterflies, but distinct from those of rested untreated butterflies. Many butterflies showed a decline in flight metabolic rate (FMR) on day 2, and there was a strong inverse relationship between the ratio of day 2 to day 1 FMR and the abundance of sodiated adducts of phosphatidylcholines and co-enzyme Q (CoQ). This result is consistent with elevation of sodiated lipids caused by disrupted intracellular ion homeostasis in mammalian tissues after hypoxia-reperfusion. Butterflies carrying the Sdhd M allele had a higher abundance of lipid markers of cellular damage, but the association was reversed in field-collected butterflies, where focal individuals typically flew for seconds at a time rather than continuously. These results indicate that Glanville fritillary flight muscles can be injured by episodes of high exertion, but injury severity appears to be determined by an interaction between SDH genotype and behavior (prolonged versus intermittent flight).

Human-Mediated Gene Flow Contributes to Metapopulation Genetic Structure of the Pathogenic Fungus Alternaria alternata from Potato

Metapopulation structure generated by recurrent extinctions and recolonizations plays an important role in the evolution of species but is rarely considered in agricultural systems. In this study, generation and mechanism of metapopulation structure were investigated by microsatellite assaying 725 isolates of Alternaria alternata sampled from potato hosts at 16 locations across China. We found a single major cluster, no isolate-geography associations and no bottlenecks in the A. alternata isolates, suggesting a metapopulation genetic structure of the pathogen. We also found weak isolation-by-distance, lower among than within cropping region population differentiation, concordant moving directions of potato products and net gene flow and the highest gene diversity in the region with the most potato imports. These results indicate that in addition to natural dispersal, human-mediated gene flow also contributes to the generation and dynamics of the metapopulation genetic structure of A. alternata in China. Metapopulation structure increases the adaptive capacity of the plant pathogen as a result of enhanced genetic variation and reduced population fragmentation. Consequently, rigid quarantine regulations may be required to reduce population connectivity and the evolutionary potential of A. alternata and other pathogens with a similar population dynamics for a sustainable plant disease management.

Extinction dynamics from metastable coexistences in an evolutionary game

Deterministic evolutionary game dynamics can lead to stable coexistences of different types. Stochasticity, however, drives the loss of such coexistences. This extinction is usually accompanied by population size fluctuations. We investigate the most probable extinction trajectory under such fluctuations by mapping a stochastic evolutionary model to a problem of classical mechanics using the Wentzel-Kramers-Brillouin (WKB) approximation. Our results show that more abundant types in a coexistence may be more likely to go extinct first, in good agreement with previous results. The distance between the coexistence and extinction points is not a good predictor of extinction either. Instead, the WKB method correctly predicts the type going extinct first.

Coevolution of patch-type dependent emigration and patch-type dependent immigration

The three phases of dispersal - emigration, transfer and immigration - are affecting each other and the former and latter decisions may depend on patch types. Despite the inevitable fact of the complexity of the dispersal process, patch-type dependencies of dispersal decisions modelled as emigration and immigration are usually missing in theoretical dispersal models. Here, I investigate the coevolution of patch-type dependent emigration and patch-type dependent immigration in an extended Hamilton-May model. The dispersing population inhabits a landscape structured into many patches of two types and disperses during a continuous-time season. The trait under consideration is a four dimensional vector consisting of two values for emigration probability from the patches and two values for immigration probability into the patches of each type. Using the adaptive dynamics approach I show that four qualitatively different dispersal strategies may evolve in different parameter regions, including a counterintuitive strategy, where patches of one type are fully dispersed from (emigration probability is one) but individuals nevertheless always immigrate into them during the dispersal season (immigration probability is one). I present examples of evolutionary branching in a wide parameter range, when the patches with high local death rate during the dispersal season guarantee a high expected disperser output. I find that two dispersal strategies can coexist after evolutionary branching: a strategy with full immigration only into the patches with high expected disperser output coexists with a strategy that immigrates into any patch. Stochastic simulations agree with the numerical predictions. Since evolutionary branching is also found when immigration evolves alone, the present study is adding coevolutionary constraints on the emigration traits and hence finds that the coevolution of a higher dimensional trait sometimes hinders evolutionary diversification.

Genetics of dispersal

Dispersal is a process of central importance for the ecological and evolutionary dynamics of populations and communities, because of its diverse consequences for gene flow and demography. It is subject to evolutionary change, which begs the question, what is the genetic basis of this potentially complex trait? To address this question, we (i) review the empirical literature on the genetic basis of dispersal, (ii) explore how theoretical investigations of the evolution of dispersal have represented the genetics of dispersal, and (iii) discuss how the genetic basis of dispersal influences theoretical predictions of the evolution of dispersal and potential consequences.

Dispersal has a detectable genetic basis in many organisms, from bacteria to plants and animals. Generally, there is evidence for significant genetic variation for dispersal or dispersal‐related phenotypes or evidence for the micro‐evolution of dispersal in natural populations. Dispersal is typically the outcome of several interacting traits, and this complexity is reflected in its genetic architecture: while some genes of moderate to large effect can influence certain aspects of dispersal, dispersal traits are typically polygenic. Correlations among dispersal traits as well as between dispersal traits and other traits under selection are common, and the genetic basis of dispersal can be highly environment‐dependent.

By contrast, models have historically considered a highly simplified genetic architecture of dispersal. It is only recently that models have started to consider multiple loci influencing dispersal, as well as non‐additive effects such as dominance and epistasis, showing that the genetic basis of dispersal can influence evolutionary rates and outcomes, especially under non‐equilibrium conditions. For example, the number of loci controlling dispersal can influence projected rates of dispersal evolution during range shifts and corresponding demographic impacts. Incorporating more realism in the genetic architecture of dispersal is thus necessary to enable models to move beyond the purely theoretical towards making more useful predictions of evolutionary and ecological dynamics under current and future environmental conditions. To inform these advances, empirical studies need to answer outstanding questions concerning whether specific genes underlie dispersal variation, the genetic architecture of context‐dependent dispersal phenotypes and behaviours, and correlations among dispersal and other traits.

Spatial and temporal genetic structure at the fourth trophic level in a fragmented landscape

A fragmented habitat becomes increasingly fragmented for species at higher trophic levels, such as parasitoids. To persist, these species are expected to possess life-history traits, such as high dispersal, that facilitate their ability to use resources that become scarce in fragmented landscapes. If a specialized parasitoid disperses widely to take advantage of a sparse host, then the parasitoid population should have lower genetic structure than the host. We investigated the temporal and spatial genetic structure of a hyperparasitoid (fourth trophic level) in a fragmented landscape over 50 × 70 km, using microsatellite markers, and compared it with the known structures of its host parasitoid, and the butterfly host which lives as a classic metapopulation. We found that population genetic structure decreases with increasing trophic level. The hyperparasitoid has fewer genetic clusters (K = 4), than its host parasitoid (K = 15), which in turn is less structured than the host butterfly (K = 27). The genetic structure of the hyperparasitoid also shows temporal variation, with genetic differentiation increasing due to reduction of the population size, which reduces the effective population size. Overall, our study confirms the idea that specialized species must be dispersive to use a fragmented host resource, but that this adaptation has limits.

Genetic effects on life-history traits in the Glanville fritillary butterfly

BACKGROUND

Adaptation to local habitat conditions may lead to the natural divergence of populations in life-history traits such as body size, time of reproduction, mate signaling or dispersal capacity. Given enough time and strong enough selection pressures, populations may experience local genetic differentiation. The genetic basis of many life-history traits, and their evolution according to different environmental conditions remain however poorly understood.

METHODS

We conducted an association study on the Glanville fritillary butterfly, using material from five populations along a latitudinal gradient within the Baltic Sea region, which show different degrees of habitat fragmentation. We investigated variation in 10 principal components, cofounding in total 21 life-history traits, according to two environmental types, and 33 genetic SNP markers from 15 candidate genes.

RESULTS

We found that nine SNPs from five genes showed strong trend for trait associations (p-values under 0.001 before correction). These associations, yet non-significant after multiple test corrections, with a total number of 1,086 tests, were consistent across the study populations. Additionally, these nine genes also showed an allele frequency difference between the populations from the northern fragmented versus the southern continuous landscape.

DISCUSSION

Our study provides further support for previously described trait associations within the Glanville fritillary butterfly species across different spatial scales. Although our results alone are inconclusive, they are concordant with previous studies that identified these associations to be related to climatic changes or habitat fragmentation within the Åland population.