Genome-Wide Allele Frequency Changes Reveal That Dynamic Metapopulations Evolve Differently

Two important characteristics of metapopulations are extinction-(re)colonization dynamics and gene flow between subpopulations. These processes can cause strong shifts in genome-wide allele frequencies that are generally not observed in "classical" (large, stable, and panmictic) populations. Subpopulations founded by one or a few individuals, the so-called propagule model, are initially expected to show intermediate allele frequencies at polymorphic sites until natural selection and genetic drift drive allele frequencies toward a mutation-selection-drift equilibrium characterized by a negative exponential-like distribution of the site frequency spectrum. We followed changes in site frequency spectrum distribution in a natural metapopulation of the cyclically parthenogenetic pond-dwelling microcrustacean Daphnia magna using biannual pool-seq samples collected over a 5-yr period from 118 ponds occupied by subpopulations of known age. As expected under the propagule model, site frequency spectra in newly founded subpopulations trended toward intermediate allele frequencies and shifted toward right-skewed distributions as the populations aged. Immigration and subsequent hybrid vigor altered this dynamic. We show that the analysis of site frequency spectrum dynamics is a powerful approach to understand evolution in metapopulations. It allowed us to disentangle evolutionary processes occurring in a natural metapopulation, where many subpopulations evolve in parallel. Thereby, stochastic processes like founder and immigration events lead to a pattern of subpopulation divergence, while genetic drift leads to converging site frequency spectrum distributions in the persisting subpopulations. The observed processes are well explained by the propagule model and highlight that metapopulations evolve differently from classical populations.

Scaling up the effects of inbreeding depression from individuals to metapopulations

Inbreeding is common in nature, and many laboratory studies have documented that inbreeding depression can reduce the fitness of individuals. Demonstrating the consequences of inbreeding depression on the growth and persistence of populations is more challenging because populations are often regulated by density- or frequency-dependent selection and influenced by demographic and environmental stochasticity. A few empirical studies have shown that inbreeding depression can increase extinction risk of local populations. The importance of inbreeding depression at the metapopulation level has been conjectured based on population-level studies but has not been evaluated. We quantified the impact of inbreeding depression affecting the fitness of individuals on metapopulation persistence in heterogeneous habitat networks of different sizes and habitat configuration in a context of natural butterfly metapopulations. We developed a spatial individual-based simulation model of metapopulations with explicit genetics. We used Approximate Bayesian Computation to fit the model to extensive demographic, genetic and life-history data available for the well-studied Glanville fritillary butterfly (Melitaea cinxia) metapopulations in the Åland islands in SW Finland. We compared 18 semi-independent habitat networks differing in size and fragmentation. The results show that inbreeding is more frequent in small habitat networks, and consequently, inbreeding depression elevates extinction risks in small metapopulations. Metapopulation persistence and neutral genetic diversity maintained in the metapopulations increase with the total habitat amount in and mean patch size of habitat networks. Dispersal and mating behaviour interact with landscape structure to determine how likely it is to encounter kin while looking for mates. Inbreeding depression can decrease the viability of small metapopulations even when they are strongly influenced by stochastic extinction-colonization dynamics and density-dependent selection. The findings from this study support that genetic factors, in addition to demographic factors, can contribute to extinctions of small local populations and also of metapopulations.

A genetic analysis of dragonfly population structure

Dragonflies reside in both aquatic and terrestrial environments, depending on their life stage, necessitating the conservation of drastically different habitats; however, little is understood about how nymph and adult dragonflies function as metapopulations within connected habitat. We used genetic techniques to examine nymphs and adults within a single metapopulation both spatially and temporally to better understand metapopulation structure and the processes that might influence said structure. We sampled 97 nymphs and 149 adult Sympetrum obtrusum from eight locations, four aquatic, and four terrestrial, at the Pierce Cedar Creek Institute in Southwest Michigan over two summers. We performed AFLP genetic analysis and used the Bayesian analysis program STRUCTURE to detect genetic clusters from sampled individuals. STRUCTURE detected k = u4 populations, in which nymphs and adults from the same locations collected in different years did not necessarily fall into the same clusters. We also evaluated grouping using the statistical clustering analyses NMDS and MRPP. The results of these confirmed findings from STRUCTURE and emphasized differences between adults collected in 2012 and all other generations. These results suggest that both dispersal and a temporal cycle of emergence of nymphs from unique clusters every other year could be influential in structuring dragonfly populations, although our methods were not able to fully distinguish the influences of either force. This study provides a better understanding of local dragonfly metapopulation structure and provides a starting point for future studies to investigate the spatial and temporal mechanisms controlling metapopulation structure. The results of the study should prove informative for managers working to preserve genetic diversity in connected dragonfly metapopulations, especially in the face of increasing anthropogenic landscape changes.

Coral reef fish populations can persist without immigration

Determining the conditions under which populations may persist requires accurate estimates of demographic parameters, including immigration, local reproductive success, and mortality rates. In marine populations, empirical estimates of these parameters are rare, due at least in part to the pelagic dispersal stage common to most marine organisms. Here, we evaluate population persistence and turnover for a population of orange clownfish, Amphiprion percula, at Kimbe Island in Papua New Guinea. All fish in the population were sampled and genotyped on five occasions at 2-year intervals spanning eight years. The genetic data enabled estimates of reproductive success retained in the same population (reproductive success to self-recruitment), reproductive success exported to other subpopulations (reproductive success to local connectivity), and immigration and mortality rates of sub-adults and adults. Approximately 50% of the recruits were assigned to parents from the Kimbe Island population and this was stable through the sampling period. Stability in the proportion of local and immigrant settlers is likely due to: low annual mortality rates and stable egg production rates, and the short larval stages and sensory capacities of reef fish larvae. Biannual mortality rates ranged from 0.09 to 0.55 and varied significantly spatially. We used these data to parametrize a model that estimated the probability of the Kimbe Island population persisting in the absence of immigration. The Kimbe Island population was found to persist without significant immigration. Model results suggest the island population persists because the largest of the subpopulations are maintained due to having low mortality and high self-recruitment rates. Our results enable managers to appropriately target and scale actions to maximize persistence likelihood as disturbance frequencies increase.

Temporal relationship between genetic and warning signal variation in the aposematic wood tiger moth (Parasemia plantaginis)

Many plants and animals advertise unpalatability through warning signals in the form of colour and shape. Variation in warning signals within local populations is not expected because they are subject to directional selection. However, mounting evidence of warning signal variation within local populations suggests that other selective forces may be acting. Moreover, different selective pressures may act on the individual components of a warning signal. At present, we have a limited understanding about how multiple selection processes operate simultaneously on warning signal components, and even less about their temporal and spatial dynamics. Here, we examined temporal variation of several wing warning signal components (colour, UV-reflectance, signal size and pattern) of two co-occurring colour morphs of the aposematic wood tiger moth (Parasemia plantaginis). Sampling was carried out in four geographical regions over three consecutive years. We also evaluated each morph's temporal genetic structure by analysing mitochondrial sequence data and nuclear microsatellite markers. Our results revealed temporal differences between the morphs for most signal components measured. Moreover, variation occurred differently in the fore- and hindwings. We found no differences in the genetic structure between the morphs within years and regions, suggesting single local populations. However, local genetic structure fluctuated temporally. Negative correlations were found between variation produced by neutrally evolving genetic markers and those of the different signal components, indicating a non-neutral evolution for most warning signal components. Taken together, our results suggest that differential selection on warning signal components and fluctuating population structure can be one explanation for the maintenance of warning signal variation in this aposematic species.

Human-facilitated metapopulation dynamics in an emerging pest species, Cimex lectularius

The number and demographic history of colonists can have dramatic consequences for the way in which genetic diversity is distributed and maintained in a metapopulation. The bed bug (Cimex lectularius) is a re-emerging pest species whose close association with humans has led to frequent local extinction and colonization, that is, to metapopulation dynamics. Pest control limits the lifespan of subpopulations, causing frequent local extinctions, and human-facilitated dispersal allows the colonization of empty patches. Founder events often result in drastic reductions in diversity and an increased influence of genetic drift. Coupled with restricted migration, this can lead to rapid population differentiation. We therefore predicted strong population structuring. Here, using 21 newly characterized microsatellite markers and approximate Bayesian computation (ABC), we investigate simplified versions of two classical models of metapopulation dynamics, in a coalescent framework, to estimate the number and genetic composition of founders in the common bed bug. We found very limited diversity within infestations but high degrees of structuring across the city of London, with extreme levels of genetic differentiation between infestations (F_ST = 0.59). ABC results suggest a common origin of all founders of a given subpopulation and that the numbers of colonists were low, implying that even a single mated female is enough to found a new infestation successfully. These patterns of colonization are close to the predictions of the propagule pool model, where all founders originate from the same parental infestation. These results show that aspects of metapopulation dynamics can be captured in simple models and provide insights that are valuable for the future targeted control of bed bug infestations.

Differences in the aerobic capacity of flight muscles between butterfly populations and species with dissimilar flight abilities

Habitat loss and climate change are rapidly converting natural habitats and thereby increasing the significance of dispersal capacity for vulnerable species. Flight is necessary for dispersal in many insects, and differences in dispersal capacity may reflect dissimilarities in flight muscle aerobic capacity. In a large metapopulation of the Glanville fritillary butterfly in the Åland Islands in Finland, adults disperse frequently between small local populations. Individuals found in newly established populations have higher flight metabolic rates and field-measured dispersal distances than butterflies in old populations. To assess possible differences in flight muscle aerobic capacity among Glanville fritillary populations, enzyme activities and tissue concentrations of the mitochondrial protein Cytochrome-c Oxidase (CytOx) were measured and compared with four other species of Nymphalid butterflies. Flight muscle structure and mitochondrial density were also examined in the Glanville fritillary and a long-distance migrant, the red admiral. Glanville fritillaries from new populations had significantly higher aerobic capacities than individuals from old populations. Comparing the different species, strong-flying butterfly species had higher flight muscle CytOx content and enzymatic activity than short-distance fliers, and mitochondria were larger and more numerous in the flight muscle of the red admiral than the Glanville fritillary. These results suggest that superior dispersal capacity of butterflies in new populations of the Glanville fritillary is due in part to greater aerobic capacity, though this species has a low aerobic capacity in general when compared with known strong fliers. Low aerobic capacity may limit dispersal ability of the Glanville fritillary.

Life history of the Glanville fritillary butterfly in fragmented versus continuous landscapes

Habitat loss and fragmentation threaten the long-term viability of innumerable species of plants and animals. At the same time, habitat fragmentation may impose strong natural selection and lead to evolution of life histories with possible consequences for demographic dynamics. The Baltic populations of the Glanville fritillary butterfly (Melitaea cinxia) inhabit regions with highly fragmented habitat (networks of small dry meadows) as well as regions with extensive continuous habitat (calcareous alvar grasslands). Here, we report the results of common garden studies on butterflies originating from two highly fragmented landscapes (FL) in Finland and Sweden and from two continuous landscapes (CL) in Sweden and Estonia, conducted in a large outdoor cage (32 by 26 m) and in the laboratory. We investigated a comprehensive set of 51 life-history traits, including measures of larval growth and development, flight performance, and adult reproductive behavior. Seventeen of the 51 traits showed a significant difference between fragmented versus CL. Most notably, the growth rate of postdiapause larvae and several measures of flight capacity, including flight metabolic rate, were higher in butterflies from fragmented than CL. Females from CL had shorter intervals between consecutive egg clutches and somewhat higher life-time egg production, but shorter longevity, than females from FL. These results are likely to reflect the constant opportunities for oviposition in females living in continuous habitats, while the more dispersive females from FL allocate more resources to dispersal capacity at the cost of egg maturation rate. This study supports theoretical predictions about small population sizes and high rate of population turnover in fragmented habitats selecting for increased rate of dispersal, but the results also indicate that many other life-history traits apart from dispersal are affected by the degree of habitat fragmentation.

High genetic load in an old isolated butterfly population

We investigated inbreeding depression and genetic load in a small (N(e) ∼ 100) population of the Glanville fritillary butterfly (Melitaea cinxia), which has been completely isolated on a small island [Pikku Tytärsaari (PT)] in the Baltic Sea for at least 75 y. As a reference, we studied conspecific populations from the well-studied metapopulation in the Åland Islands (ÅL), 400 km away. A large population in Saaremaa, Estonia, was used as a reference for estimating genetic diversity and N(e). We investigated 58 traits related to behavior, development, morphology, reproductive performance, and metabolism. The PT population exhibited high genetic load (L = 1 - W(PT)/W(ÅL)) in a range of fitness-related traits including adult weight (L = 0.12), flight metabolic rate (L = 0.53), egg viability (L = 0.37), and lifetime production of eggs in an outdoor population cage (L = 0.70). These results imply extensive fixation of deleterious recessive mutations, supported by greatly reduced diversity in microsatellite markers and immediate recovery (heterosis) of egg viability and flight metabolic rate in crosses with other populations. There was no significant inbreeding depression in most traits due to one generation of full-sib mating. Resting metabolic rate was significantly elevated in PT males, which may be related to their short lifespan (L = 0.25). The demographic history and the effective size of the PT population place it in the part of the parameter space in which models predict mutation accumulation. This population exemplifies the increasingly common situation in fragmented landscapes, in which small and completely isolated populations are vulnerable to extinction due to high genetic load.

Eco-evolutionary spatial dynamics in the Glanville fritillary butterfly

Demographic population dynamics, gene flow, and local adaptation may influence each other and lead to coupling of ecological and evolutionary dynamics, especially in species inhabiting fragmented heterogeneous environments. Here, I review long-term research on eco-evolutionary spatial dynamics in the Glanville fritillary butterfly inhabiting a large network of approximately 4,000 meadows in Finland. The metapopulation persists in a balance between frequent local extinctions and recolonizations. The genetic spatial structure as defined by neutral markers is much more coarse-grained than the demographic spatial structure determined by the fragmented habitat, yet small-scale spatial structure has important consequences for the dynamics. I discuss three examples of eco-evolutionary spatial dynamics. (i) Extinction-colonization metapopulation dynamics influence allele frequency changes in the phosphoglucose isomerase (Pgi) gene, which leads to strong associations between genetic variation in Pgi and dispersal, recolonization, and local population dynamics. (ii) Inbreeding in local populations increases their risk for extinction, whereas reciprocal effects between inbreeding, population size, and emigration represent likely eco-evolutionary feedbacks. (iii) Genetically determined female oviposition preference for two host plant species exhibits a cline paralleling a gradient in host plant relative abundances, and host plant preference of dispersing females in relation to the host plant composition of habitat patches influences immigration (gene flow) and recolonization (founder events). Eco-evolutionary spatial dynamics in heterogeneous environments may not lead to directional evolutionary changes unless the environment itself changes, but eco-evolutionary dynamics may contribute to the maintenance of genetic variation attributable to fluctuating selection in space and time.