Functional genomics of life history variation in a butterfly metapopulation

In fragmented landscapes, small populations frequently go extinct and new ones are established with poorly understood consequences for genetic diversity and evolution of life history traits. Here, we apply functional genomic tools to an ecological model system, the well-studied metapopulation of the Glanville fritillary butterfly. We investigate how dispersal and colonization select upon existing genetic variation affecting life history traits by comparing common-garden reared 2-day adult females from new populations with those from established older populations. New-population females had higher expression of abdomen genes involved in egg provisioning and thorax genes involved in the maintenance of flight muscle proteins. Physiological studies confirmed that new-population butterflies have accelerated egg maturation, apparently regulated by higher juvenile hormone titer and angiotensin converting enzyme mRNA, as well as enhanced flight metabolism. Gene expression varied between allelic forms of two metabolic genes (Pgi and Sdhd), which themselves were associated with differences in flight metabolic rate, population age and population growth rate. These results identify likely molecular mechanisms underpinning life history variation that is maintained by extinction-colonization dynamics in metapopulations.     

A candidate locus for variation in dispersal rate in a butterfly metapopulation

Frequent extinctions of local populations in metapopulations create opportunities for migrant females to establish new populations. In a metapopulation of the Glanville fritillary butterfly (Melitaea cinxia), more mobile individuals are more likely to establish new populations, especially in habitat patches that are poorly connected to existing populations. Here we show that flight metabolic rate and the frequency of a specific allele of the metabolic enzyme phosphoglucose isomerase (pgi) were both highest in newly established, isolated populations. Furthermore, genotypes with this pgi allele had elevated flight metabolic rates. These results suggest that genetic variation in pgi or a closely linked locus has a direct effect on flight metabolism, dispersal rate, and thereby on metapopulation dynamics in this species. These results also contribute to an emerging understanding of the mechanisms by which population turnover in heterogeneous landscapes may maintain genetic and phenotypic variation across populations.     

The genomic and physiological basis of life history variation in a butterfly metapopulation

Unravelling the mechanisms underlying variation in life history traits is of fundamental importance for our understanding of adaptation by natural selection. While progress has been made in mapping fitness-related phenotypes to genotypes, mainly in a handful of model organisms, functional genomic studies of life history adaptations are still in their infancy. In particular, despite a few notable exceptions, the genomic basis of life history variation in natural populations remains poorly understood. This is especially true for the genetic underpinnings of life history phenotypes subject to diversifying selection driven by ecological dynamics in patchy environments--as opposed to adaptations involving strong directional selection owing to major environmental changes, such as latitudinal gradients, extreme climatic events or transitions from salt to freshwater. In this issue of Molecular Ecology,Wheat et al. (2011) now make a significant leap forward by applying the tools of functional genomics to dispersal-related life history variation in a butterfly metapopulation. Using a combination of microarrays, quantitative PCR and physiological measurements, the authors uncover several metabolic and endocrine factors that likely contribute to the observed life history phenotypes. By identifying molecular candidate mechanisms of fitness variation maintained by dispersal dynamics in a heterogeneous environment,they also begin to address fascinating interactions between the levels of physiology, ecology and evolution.     

Dispersal-related life-history trade-offs in a butterfly metapopulation

1. Recent studies on butterflies have documented apparent evolutionary changes in dispersal rate in response to climate change and habitat change. These studies often assume a trade-off between dispersal rate (or flight capacity) and reproduction, which is the rule in wing-dimorphic species but might not occur equally in wing-monomorphic species such as butterflies. 2. To investigate the relationship between dispersal rate and fecundity in the Glanville fritillary butterfly Melitaea cinxia we recorded lifetime individual movements, matings, ovipositions, and maximal life span in a large (32 x 26 m) population cage in the field. Experimental material was obtained from 20 newly established and 20 old local populations within a large metapopulation in the Aland Islands in Finland. 3. Females of the Glanville fritillary from newly established populations are known to be more dispersive in the field, and in the cage they showed significantly greater mobility, mated earlier, and laid more egg clutches than females from old populations. The dispersive females from new populations exhibited no reduced lifetime fecundity in the cage, but they had a shorter maximal life span than old-population females. 4. These results challenge the dispersal-fecundity trade-off for nonmigratory butterflies but instead suggest a physiological trade-off between high metabolic performance and reduced maximal life span. High metabolic performance may explain high rates of dispersal and oviposition in early life. 5. In fragmented landscapes, an ecological trade-off exists between being more dispersive and hence spending more time in the landscape matrix vs. having more time for reproduction in the habitat. We estimate with a dispersal model parameterized for the Glanville fritillary that the lifetime egg production is 4% smaller on average in the more dispersive butterflies in a representative landscape, with much variation depending on landscape structure in the neighbourhood of the natal patch, from--26 to 45% in the landscape analysed in this paper.     

Correlated extinctions, colonizations and population fluctuations in a highly connected ringlet butterfly metapopulation

 The persistence of metapopulations is likely to be highly dependent on whether population dynamics are correlated among habitat patches as a result of migration between patches and spatially-correlated environmental stochasticity (weather effects). We examined whether population dynamics of the ringlet butterfly, Aphantopus hyperantus, were synchronous in an area of approximately 0.5 km², with respect to extinction, colonization and population fluctuations. Monks Wood Butterfly Monitoring Scheme transect count data from 1973 to 1995, revealed (A) a major environmental perturbation, the drought of 1976, which caused synchronized extinctions of A. hyperantus in subsequent years, (B) synchronized recolonization in years following the large number of apparent extinctions, and (C) population changes by A. hyperantus were highly correlated in many of the 14 sections of the transect, presumably reflecting similar responses to environmental stochasticity, and the exchange of individuals among sections. However, extinction and population synchrony depended on habitat type. Following the 1976 drought, A. hyperantus apparently became extinct from the most open and most shady habitats it occupied, with some persistence in habitats of intermediate shading, thus showing retraction to core populations in central parts of an environmental gradient, albeit with an average shift to relatively open habitat. Populations at extreme ends of the environmental gradient occupied by A. hyperantus fluctuated least synchronously, suggesting a potential buffering effect of habitat heterogeneity, but this was not crucial to survival after the 1976 drought. Thus, not all habitats are equally important to persistence. Correlated temporal dynamics, variation in habitat quality and the interaction between habitat quality and temporal environmental stochasticity are important determinants of metapopulation persistence and should be incorporated in metapopulation models. Received: 26 April 1996 / Accepted: 17 July 1996     

Short-term studies underestimate 30-generation changes in a butterfly metapopulation

Most studies of rare and endangered species are based on work carried out within one generation, or over one to a few generations of the study organism. We report the results of a study that spans 30 generations (years) of the entire natural range of a butterfly race that is endemic to 35 km(2) of north Wales, UK. Short-term studies (surveys in single years and dynamics over 4 years) of this system led to the prediction that the regional distribution would be quite stable, and that colonization and extinction dynamics would be relatively unimportant. However, a longer-term study revealed unexpectedly high levels of population turnover (local extinction and colonization), affecting 18 out of the 20 patches that were occupied at any time during the period. Modelling the system (using the 'incidence function model' (IFM) for metapopulations) also showed higher levels of colonization and extinction with increasing duration of the study. The longer-term dynamics observed in this system can be compared, at a metapopulation level, with the increased levels of variation observed with increasing time that have been observed in single populations. Long-term changes may arise from local changes in the environment that make individual patches more or less suitable for the butterfly, or from unusual colonization or extinction events that take metapopulations into alternative states. One implication is that metapopulation and population viability analyses based on studies that cover only a few animal or plant generations may underestimate extinction threats.     

Spatial variation in population growth rate and community structure affects local and regional dynamics

1. Theory predicting that populations with high maximum rates of increase (r(max)) will be less stable, and that metapopulations with high average r(max) will be less synchronous, was tested using a small protist, Bodo, that inhabits pitcher plant leaves (Sarracenia purpurea L.). The effects of predators and resources on these relationships were also determined. 2. Abundance data collected for a total of 60 populations of Bodo, over a period of 3 months, at six sites in three bogs in eastern Canada, were used to test these predictions. Mosquitoes were manipulated in half the leaves partway through the season to increase the range of predation rates. 3. Dynamics differed greatly among leaves and sites, but most populations exhibited one or more episodes of rapid increase followed by a population crash. Estimates of r(max) obtained using a linear mixed-effects model, ranged from 1 x 5 to 2 x 7 per day. Resource levels (captured insect) and midge abundances affected r(max). 4. Higher r(max) was associated with greater temporal variability and lower synchrony as predicted. However, in contrast to expectations, populations with higher r(max) also had lower mean abundance and were more suppressed by predators. 5. This study demonstrates that the link between r(max) and temporal variability is key to understanding the dynamics of populations that spend little time near equilibrium, and to predicting and interpreting the effects of community structure on the dynamics of such populations.     

Size and genetic composition of the colonizing propagules in a butterfly metapopulation

The amount and spatial distribution of genetic variation that is maintained in a metapopulation depends critically on the colonization process. Here, we use molecular markers to determine the number and genetic relatedness of individuals establishing new local populations in a large metapopulation of the Glanville fritillary butterfly Melitaea cinxia. The empirical results are compared with the predictions of a dispersal model based on a diffusion approximation of correlated random walk, which serves as a base‐line hypothesis about the rate and pattern of colonization. The results show that half of the new local populations consisted of a single larval group of full sibs and hence necessarily of the offspring of a single female. If the colonization involved two or more larval groups, these were usually oviposited by two different females that were unrelated to each other. The pattern of colonizations is thus intermediate between the propagule pool and the migrant pool models. These results elucidate the generation of genetic stochasticity, which may influence the dynamics of small populations. The dispersal model predicted well the pattern of habitat occupancy and the pattern of colonizations in relation to landscape structure, though which particular habitat patches became colonized was influenced also by measures of habitat quality not included in the model.     

Microsatellite variation and population structure of a recovering Tree frog (Hyla arborea L.) metapopulation

Numbers and sizes of populations of the European tree frog in The Netherlands have dramatically decreased in the second half of the last century due to extensive habitat destruction and fragmentation. We have studied the genetic structure of a slowly recovering meta-population. Strong genetic differentiation, estimated at eight microsatellite loci, was found between clusters of populations (F _st-values above 0.2). Within clusters, consisting of ponds within a radius of about 5 km, European tree frog populations were less differentiated (F _st<0.08) and exact tests showed that most of the ponds within clusters were not significantly differentiated. Although local population sizes have been increasing since 1989, and some new ponds have been colonised in the direct vicinity of ponds that have been populated continuously, little evidence for gene flow between clusters of ponds was found (only one exception). Furthermore, levels of genetic diversity were low compared to populations in comparable areas elsewhere in Europe. Therefore, a continuous conservation effort is needed to prevent any further loss of genetic diversity. The alleviation of important barriers to dispersal between the clusters should be given a high priority for the restoration of the meta-population as a whole.     

Long distance dispersal and landscape occupancy in a metapopulation of the cranberry fritillary butterfly

Movements between habitat patches in a patchy population of the butterfly Boloria aquilonaris were monitored using capture-mark-recapture methods during three successive generations. For each data set, the inverse cumulative proportion of individuals moving 100 m distance classes was fitted to the negative exponential function and the inverse power function. In each case, the negative exponential function provided a better fit than the inverse power function. Two dispersal kernels were generated using both negative exponential and inverse power functions. These dispersal kernels were used to predict movements between 14 suitable sites in a landscape of 220 km². The negative exponential function generated a dispersal kernel predicting extremely low probabilities for movements exceeding 1 km. The inverse power function generated probabilities predicting that between site movements were possible, according to metapopulation size. CMR studies in the landscape revealed that long distance movements occurred at each generation, corresponding to predictions of the inverse power function dispersal kernel. A total of 26 movements between sites (up to 13 km) were detected, together with recolonisation of empty sites. The spatial scale of the metapopulation dynamics is larger than ever reported on butterflies and long distance movements clearly matter to the persistence of this species in a highly fragmented landscape.     

Three ways of assessing metapopulation structure in the butterfly Plebejus argus

1. Three independent methods were used to investigate population structure in the butterfly Plebejus argus . First, migration and dispersal ability were measured by mark–release–recapture in seven adjacent habitat patches, and by release of butterflies in unoccupied habitat. Secondly, colonization of newly created habitat was observed over 7 years. Finally, genetic differentiation of local populations within a metapopulation was investigated. Sampled local populations included parts of the mark–release–recapture study area.
2. Plebejus argus is relatively sedentary: the maximum movement detected was 395 m, and only 2% of individuals moved further than 100 m between recaptures on different days. None the less, adjacent local populations in the mark–release–recapture study area were linked by occasional migration, with ≈ 1.4% of individuals moving between patches separated by 13–200 m.
3. Despite low mobility, observed colonizations occurred rapidly over distances of 1 km. Because P. argus occurs at high population densities, 1.4% migration can generate enough migrants to colonize newly suitable habitat quickly at this spatial scale.
4. Mark–release–recapture data were used to predict that there would be limited genetic differentiation through drift between local populations at this spatial scale. The prediction was supported by allele frequency data for the same local populations.
5. Genetic differentiation often indicates higher levels of migration than are revealed by the movements of marked individuals. This study shows that when experimental releases and extensive marking are undertaken in areas that are large relative to most movements, indirect measures of gene flow and direct measures of dispersal can concur.
6. Evidence from the three different approaches was complementary, indicating that P. argus occurs as metapopulations within the study area.     

Bayesian state-space modeling of metapopulation dynamics in the Glanville fritillary butterfly

The complexity of mathematical models of ecological dynamics varies greatly, and it is often difficult to judge what would be the optimal level of complexity in a particular case. Here we compare the parameter estimates, model fits, and predictive abilities of two models of metapopulation dynamics: a detailed individual-based model (IBM) and a population-based stochastic patch occupancy model (SPOM) derived from the IBM. The two models were fitted to a 17-year time series of data for the Glanville fritillary butterfly (Melitaea cinxia) inhabiting a network of 72 small meadows. The data consisted of biannual counts of larval groups (IBM) and the annual presence or absence of local populations (SPOM). The models were fitted using a Bayesian state-space approach with a hierarchical random effect structure to account for observational, demographic, and environmental stochasticities. The detection probability of larval groups (IBM) and the probability of false zeros of local populations (SPOM) in the observation models were simultaneously estimated from the time-series data and independent control data. Prior distributions for dispersal parameters were obtained from a separate analysis of mark–recapture data. Both models fitted the data about equally, but the results were more precise for the IBM than for the SPOM. The two models yielded similar estimates for a random effect parameter describing habitat quality in each patch, which were correlated with independent empirical measures of habitat quality. The modeling results showed that variation in habitat quality influenced patch occupancy more through the effects on movement behavior at patch edges than on carrying capacity, whereas the latter influenced the mean population size in occupied patches. The IBM and the SPOM explained 63% and 45%, respectively, of the observed variation in the fraction of occupied habitat area among 75 independent patch networks not used in parameter estimation. We conclude that, while carefully constructed, detailed models can have better predictive ability than simple models, this advantage comes with the cost of greatly increased data requirements and computational challenges. Our results illustrate how complex models can be helpful in facilitating the construction of effective simpler models.     

Potential metapopulation structure and the effects of habitat quality on population size of the endangered False Ringlet butterfly

The False Ringlet (Coenonympha oedippus) is a European butterfly species, endangered due to the severe loss and fragmentation of its habitat. In Hungary, two remaining populations of the butterfly occur in lowland Purple Moorgrass meadows. We studied a metapopulation occupying twelve habitat patches in Central Hungary. Our aim was to reveal what measures of habitat quality affect population size and density of this metapopulation, estimate dispersal parameters and describe phenology of subpopulations. Local population sizes and dispersal parameters were estimated from an extensive mark–release–recapture dataset, while habitat quality was characterized by groundwater level, cover of grass tussocks, bush cover, height of vegetation and grass litter at each habitat patch. The estimated size of the metapopulation was more than 3,000 individuals. We estimated a low dispersal capacity, especially for females, indicating a very low probability of (re)colonization. Butterfly abundance and density in local populations increased with higher grass litter, lower groundwater level and larger area covered by tussocks. We suppose that these environmental factors affect butterfly abundance by determining the microclimatic conditions for both larvae and adult butterflies. Our results suggest that the long-term preservation of the studied metapopulation needs the maintenance of high quality habitat patches by appropriate mowing regime and water regulation. Management also should facilitate dispersal to strengthen metapopulation structure with creating stepping-stones or gradually increase habitat quality in present matrix.     

Genetic spatial structure in a butterfly metapopulation correlates better with past than present demographic structure

The Glanville fritillary butterfly ( Melitaea cinxia ) has been studied in the Åland Islands in Finland since 1991, where it occurs as a classic metapopulation in a large network of 4000 dry meadows. Much ecological work has been conducted on this species, but population genetic studies have been hampered by paucity of suitable genetic markers. Here, using single nucleotide polymorphisms and microsatellites developed for the Glanville fritillary, we examine the correspondence between the demographic and genetic spatial structures. Given the dynamic nature of the metapopulation, the current genetic spatial structure may bear a signal of past changes in population sizes and past patterns of gene flow rather than reflect the current demographic structure or landscape structure. We analyse this question with demographic data for 10 years, using the Rand index to assess the similarity between the genetic, demographic, and landscape spatial structures. Our results show that the current genetic spatial structure is better explained by the past rather than by the current demographic spatial structure or by the spatial configuration of the habitat in the landscape. Furthermore, current genetic diversity is significantly explained by past metapopulation sizes. The time lag between major demographic events and change in the genetic spatial structure and diversity has implications for the study of spatial dynamics.     

Dynamic populations in a dynamic landscape: the metapopulation structure of the marsh fritillary butterfly

The marsh fritillary butterfly Euphydryas aurinia is an endangered species in most of northern Europe. We describe the metapopulation structure of E. aurinia in Finland, where the species has declined drastically in the past decades. We found two types of habitat patches suitable for the species: semi-permanent meadows and transient clearcuts in the forest. Patch area was the most significant variable predicting the occurrence of E. aurinia in a habitat patch. The species tended to be found in young rather than old clearcuts, apparently because the vegetation became too high in the latter. We used the incidence function model to simulate the metapopulation dynamics of E. aurinia in its dynamic landscape and discovered that the continued presence of the semi-permanent meadows is essential for the survival of the species in the study area in southeast Finland.     

Spatial and temporal variability in a butterfly population

Summary The dynamics of a butterfly (Plebejus argus) population were analysed at two levels, (i) the population as a whole and (ii) sections within the population. Some sections of the population fluctuated out of synchrony with others, such that the variability [SD Log(Density+1)] shown by the population as a whole was less than the variability shown by each part of the population — overall temporal variability was dampened by spatial asynchrony. Since observed population variability depends on the spatial scale that is sampled, comparisons of population variability among taxa should be carried out only with caution. Implications for island biogeography and conservation biology are discussed.     

Non-random dispersal in the butterfly Maniola jurtina: implications for metapopulation models

The dispersal patterns of animals are important in metapopulation ecology because they affect the dynamics and survival of populations. Theoretical models assume random dispersal but little is known in practice about the dispersal behaviour of individual animals or the strategy by which dispersers locate distant habitat patches. In the present study, we released individual meadow brown butterflies (Maniola jurtina) in a non-habitat and investigated their ability to return to a suitable habitat. The results provided three reasons for supposing that meadow brown butterflies do not seek habitat by means of random flight. First, when released within the range of their normal dispersal distances, the butterflies orientated towards suitable habitat at a higher rate than expected at random. Second, when released at larger distances from their habitat, they used a non-random, systematic, search strategy in which they flew in loops around the release point and returned periodically to it. Third, butterflies returned to a familiar habitat patch rather than a non-familiar one when given a choice. If dispersers actively orientate towards or search systematically for distant habitat, this may be problematic for existing metapopulation models, including models of the evolution of dispersal rates in metapopulations.