Transcriptome Analysis Reveals Signature of Adaptation to Landscape Fragmentation

We characterize allelic and gene expression variation between populations of the Glanville fritillary butterfly (Melitaea cinxia) from two fragmented and two continuous landscapes in northern Europe. The populations exhibit significant differences in their life history traits, e.g. butterflies from fragmented landscapes have higher flight metabolic rate and dispersal rate in the field, and higher larval growth rate, than butterflies from continuous landscapes. In fragmented landscapes, local populations are small and have a high risk of local extinction, and hence the long-term persistence at the landscape level is based on frequent re-colonization of vacant habitat patches, which is predicted to select for increased dispersal rate. Using RNA-seq data and a common garden experiment, we found that a large number of genes (1,841) were differentially expressed between the landscape types. Hexamerin genes, the expression of which has previously been shown to have high heritability and which correlate strongly with larval development time in the Glanville fritillary, had higher expression in fragmented than continuous landscapes. Genes that were more highly expressed in butterflies from newly-established than old local populations within a fragmented landscape were also more highly expressed, at the landscape level, in fragmented than continuous landscapes. This result suggests that recurrent extinctions and re-colonizations in fragmented landscapes select a for specific expression profile. Genes that were significantly up-regulated following an experimental flight treatment had higher basal expression in fragmented landscapes, indicating that these butterflies are genetically primed for frequent flight. Active flight causes oxidative stress, but butterflies from fragmented landscapes were more tolerant of hypoxia. We conclude that differences in gene expression between the landscape types reflect genomic adaptations to landscape fragmentation.     

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.     

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.     

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.     

Butterfly species differing in mobility show different structures of dispersal‐related syndromes in the same fragmented landscape

Mobility varies strongly between and within species, reflecting different dispersal strategies. Within species, such differences can imply suites of traits associated in syndromes. Different syndrome structures have been found within species among populations differing in the selective pressures they are exposed to. Similarly, we expect species differing in mobility to show different syndrome structures in response to similar selective pressures such as landscape fragmentation. Using butterflies originating from the same fragmented landscape, we investigated the differences in mobility syndrome between four common butterflies (Pyronia tithonus, Pararge aegeria Maniola jurtina, Pieris rapae) known to differ in their mobility. We expected individuals from the less mobile species to display a resident strategy because of high dispersal cost in this fragmented landscape, and individuals from the more mobile species to display a larger range of movement strategies. Moreover, as syndromes can only be detected whenever individuals differ in their dispersal strategies, we expected mobility syndromes to be observable only in populations where dispersal polymorphism is maintained. We thus expected stronger correlations between mobility‐related traits in more mobile species. Using three mobility tests in controlled conditions designed to measure different components of mobility, we showed that mobility‐related traits were indeed correlated only in the most mobile species. The absence of correlation in the less mobile species may be explained by a low variation in movement strategies, dispersal being counter‐selected.     

Long‐term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia): survey methods,data management,and long‐term population trends

Long‐term observational studies conducted at large (regional) spatial scales contribute to better understanding of landscape effects on population and evolutionary dynamics, including the conditions that affect long‐term viability of species, but large‐scale studies are expensive and logistically challenging to keep running for a long time. Here, we describe the long‐term metapopulation study of the Glanville fritillary butterfly (Melitaea cinxia) that has been conducted since 1991 in a large network of 4000 habitat patches (dry meadows) within a study area of 50 by 70 km in the Åland Islands in Finland. We explain how the landscape structure has been described, including definition, delimitation, and mapping of the habitat patches; methods of field survey, including the logistics, cost, and reliability of the survey; and data management using the EarthCape biodiversity platform. We describe the long‐term metapopulation dynamics of the Glanville fritillary based on the survey. There has been no long‐term change in the overall size of the metapopulation, but the level of spatial synchrony and hence the amplitude of fluctuations in year‐to‐year metapopulation dynamics have increased over the years, possibly due to increasing frequency of exceptional weather conditions. We discuss the added value of large‐scale and long‐term population studies, but also emphasize the need to integrate more targeted experimental studies in the context of long‐term observational studies. For instance, in the case of the Glanville fritillary project, the long‐term study has produced an opportunity to sample individuals for experiments from local populations with a known demographic history. These studies have demonstrated striking differences in dispersal rate and other life‐history traits of individuals from newly established local populations (the offspring of colonizers) versus individuals from old, established local populations. The long‐term observational study has stimulated the development of metapopulation models and provided an opportunity to test model predictions. This combination of empirical studies and modeling has facilitated the study of key phenomena in spatial dynamics, such as extinction threshold and extinction debt.     

Strong dispersal in a parasitoid wasp overwhelms habitat fragmentation and host population dynamics

The population dynamics of a parasite depend on species traits, host dynamics and the environment. Those dynamics are reflected in the genetic structure of the population. Habitat fragmentation has a greater impact on parasites than on their hosts because resource distribution is increasingly fragmented for species at higher trophic levels. This could lead to either more or less genetic structure than the host, depending on the relative dispersal rates of species. We examined the spatial genetic structure of the parasitoid wasp Hyposoter horticola, and how it was influenced by dispersal, host population dynamics and habitat fragmentation. The host, the Glanville fritillary butterfly, lives as a metapopulation in a fragmented landscape in the Åland Islands, Finland. We collected wasps throughout the 50 by 70 km archipelago and determined the genetic diversity, spatial population structure and genetic differentiation using 14 neutral DNA microsatellite loci. We compared the genetic structure of the wasp with that of the host butterfly using published genetic data collected over the shared landscape. Using maternity assignment, we also identified full‐siblings among the sampled parasitoids to estimate the dispersal range of individual females. We found that because the parasitoid is dispersive, it has low genetic structure, is not very sensitive to habitat fragmentation and has less spatial genetic structure than its butterfly host. The wasp is sensitive to regional rather than local host dynamics, and there is a geographic mosaic landscape for antagonistic co‐evolution of host resistance and parasite virulence.     

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.     

Thermal biology of flight in a butterfly: genotype,flight metabolism,and environmental conditions

Knowledge of the effects of thermal conditions on animal movement and dispersal is necessary for a mechanistic understanding of the consequences of climate change and habitat fragmentation. In particular, the flight of ectothermic insects such as small butterflies is greatly influenced by ambient temperature. Here, variation in body temperature during flight is investigated in an ecological model species, the Glanville fritillary butterfly (Melitaea cinxia). Attention is paid on the effects of flight metabolism, genotypes at candidate loci, and environmental conditions. Measurements were made under a natural range of conditions using infrared thermal imaging. Heating of flight muscles by flight metabolism has been presumed to be negligible in small butterflies. However, the results demonstrate that Glanville fritillary males with high flight metabolic rate maintain elevated body temperature better during flight than males with a low rate of flight metabolism. This effect is likely to have a significant influence on the dispersal performance and fitness of butterflies and demonstrates the possible importance of intraspecific physiological variation on dispersal in other similar ectothermic insects. The results also suggest that individuals having an advantage in low ambient temperatures can be susceptible to overheating at high temperatures. Further, tolerance of high temperatures may be important for flight performance, as indicated by an association of heat‐shock protein (Hsp70) genotype with flight metabolic rate and body temperature at takeoff. The dynamics of body temperature at flight and factors affecting it also differed significantly between female and male butterflies, indicating that thermal dynamics are governed by different mechanisms in the two sexes. This study contributes to knowledge about factors affecting intraspecific variation in dispersal‐related thermal performance in butterflies and other insects. Such information is needed for predictive models of the evolution of dispersal in the face of habitat fragmentation and climate change.     

Heritability of flight and resting metabolic rates in the Glanville fritillary butterfly

Dispersal capacity is a key life‐history trait especially in species inhabiting fragmented landscapes. Evolutionary models predict that, given sufficient heritable variation, dispersal rate responds to natural selection imposed by habitat loss and fragmentation. Here, we estimate phenotypic variance components and heritability of flight and resting metabolic rates (RMRs) in an ecological model species, the Glanville fritillary butterfly, in which flight metabolic rate (FMR) is known to correlate strongly with dispersal rate. We modelled a two‐generation pedigree with the animal model to distinguish additive genetic variance from maternal and common environmental effects. The results show that FMR is significantly heritable, with additive genetic variance accounting for about 40% of total phenotypic variance; thus, FMR has the potential to respond to selection on dispersal capacity. Maternal influences on flight metabolism were negligible. Heritability of flight metabolism was context dependent, as in stressful thermal conditions, environmentally induced variation dominated over additive genetic effects. There was no heritability in RMR, which was instead strongly influenced by maternal effects. This study contributes to a mechanistic understanding of the evolution of dispersal‐related traits, a pressing question in view of the challenges posed to many species by changing climate and fragmentation of natural habitats.     

Variation in migration propensity among individuals maintained by landscape structure

Metapopulation dynamics lead to predictable patterns of habitat occupancy, population density and trophic structure in relation to landscape features such as habitat patch size and isolation. Comparable patterns may occur in behavioural, physiological and life‐history traits but remain little studied. In the Glanville fritillary butterfly, females in newly established populations were more mobile than females in old populations. Among females from new populations, mobility decreased with increasing connectivity (decreasing isolation), but in females from old populations mobility increased with connectivity. The [ATP]/[ADP] ratio of flight muscles following controlled activity showed the same pattern as mobility in relation to population age and connectivity, suggesting that physiological differences in flight metabolic performance contribute to the observed variation in mobility. We demonstrate with an evolutionary metapopulation model parameterised for the Glanville fritillary that increasing spatial variation in landscape structure increases variance in mobility among individuals in a metapopulation, supporting the general notion that complex landscape structure maintains life‐history variation.     

Heritability of dispersal rate and other life history traits in the Glanville fritillary butterfly

Knowing the variances and heritabilities (h(2)) of life history traits in populations living under natural conditions is necessary for a mechanistic understanding of respective evolutionary processes. I estimated heritabilities of several life history traits, including dispersal rate, body mass, age at first reproduction, egg mass, clutch size and lifetime reproductive success, in the Glanville fritillary butterfly (Melitaea cinxia) using parent-offspring regression. Experiments were conducted under field conditions in a large population cage (32 x 26 m). Heritability estimates ranged from zero to almost one and several were significantly different from zero. Body size for both sexes, female age at first reproduction and egg weight were all moderately to highly heritable, whereas heritabilities were low or non-existent in clutch size and lifetime egg production. Heritability estimates for dispersal rate varied between the sexes, so that dispersal was heritable from mother to her female offspring only. This finding is consistent with previous results showing that the F1 female but not male offspring of females that naturally established new populations in the field are significantly more dispersive than butterflies in old populations.     

Metapopulation-level adaptation of insect host plant preference and extinction-colonization dynamics in heterogeneous landscapes

Species living in highly fragmented landscapes typically occur as metapopulations with frequent turnover of local populations. The turnover rate depends on population sizes and connectivities, but it may also depend on the phenotypic and genotypic composition of populations. The Glanville fritillary butterfly (Melitaea cinxia) in Finland uses two host plant species, which show variation in their relative abundances at two spatial scales: locally among individual habitat patches and regionally among networks of patches. Female butterflies in turn exhibit spatial variation in genetically determined host plant preference within and among patch networks. Emigration, immigration and establishment of new populations have all been shown to be strongly influenced by the match between the host plant composition of otherwise suitable habitat patches and the host plant preference of migrating butterflies. The evolutionary consequences of such biased migration and colonization with respect to butterfly phenotypes might differ depending on spatial configuration and plant species composition of the patches in heterogeneous patch networks. Using a spatially realistic individual-based model we show that the model-predicted evolution of host plant preference due to biased migration explains a significant amount of the observed variation in host plant use among metapopulations living in dissimilar networks. This example illustrates how the ecological extinction-colonization dynamics may be linked with the evolutionary dynamics of life history traits in metapopulations.     

Landscape structure shapes habitat finding ability in a butterfly

Land-use intensification and habitat fragmentation is predicted to impact on the search strategies animals use to find habitat. We compared the habitat finding ability between populations of the speckled wood butterfly (Pararge aegeria L.) from landscapes that differ in degree of habitat fragmentation. Naïve butterflies reared under standardized laboratory conditions but originating from either fragmented agricultural landscapes or more continuous forested landscapes were released in the field, at fixed distances from a target habitat patch, and their flight paths were recorded. Butterflies originating from fragmented agricultural landscapes were better able to find a woodlot habitat from a distance compared to conspecifics from continuous forested landscapes. To manipulate the access to olfactory information, a subset of individuals from both landscape types were included in an antennae removal experiment. This confirmed the longer perceptual range for butterflies from agricultural landscapes and indicated the significance of both visual and olfactory information for orientation towards habitat. Our results are consistent with selection for increased perceptual range in fragmented landscapes to reduce dispersal costs. An increased perceptual range will alter the functional connectivity and thereby the chances for population persistence for the same level of structural connectivity in a fragmented landscape.     

Mobility and lifetime fecundity in new versus old populations of the Glanville fritillary butterfly

Life history theory often assumes a trade-off between dispersal and reproduction, and such a trade-off is commonly observed in wing-dimorphic insects. The results are less consistent for wing-monomorphic species, for which it is more difficult to assess dispersal capacity and rate. Three replicate experiments were carried out in consecutive years on the Glanville fritillary butterfly in a large outdoor population cage to study the relationship between lifetime egg production and mobility. The experimental material included females originating from newly-established and old populations, as previous studies have shown dispersal capacity to depend on population age. There was a consistent and significant interaction between mobility and population age, such that in newly-established populations mobile females had higher fecundity than less mobile females, while in old populations there was no such relationship. As selection favours individuals with the highest fecundity, selection pressure on mobility is likely to be different between the two population types, which may contribute to maintenance of variation in dispersal rate in the metapopulation as a whole. Several other female traits also affected lifetime fecundity, including lifespan, number of matings and date of eclosion, although these effects were not consistent across the years. These results highlight the importance of conducting experiments in more than one year before generalizing about patterns in life history variation.     

Eco-evolutionary dynamics of dispersal in spatially heterogeneous environments

Ecology Letters (2011) 14: 1025-1034 ABSTRACT: Evolutionary changes in natural populations are often so fast that the evolutionary dynamics may influence ecological population dynamics and vice versa. Here we construct an eco-evolutionary model for dispersal by combining a stochastic patch occupancy metapopulation model with a model for changes in the frequency of fast-dispersing individuals in local populations. We test the model using data on allelic variation in the gene phosphoglucose isomerase (Pgi), which is strongly associated with dispersal rate in the Glanville fritillary butterfly. Population-specific measures of immigration and extinction rates and the frequency of fast-dispersing individuals among the immigrants explained 40% of spatial variation in Pgi allele frequency among 97 local populations. The model clarifies the roles of founder events and gene flow in dispersal evolution and resolves a controversy in the literature about the consequences of habitat loss and fragmentation on the evolution of dispersal.     

From individual behavior to metapopulation dynamics: unifying the patchy population and classic metapopulation models

Spatially structured populations in patchy habitats show much variation in migration rate, from patchy populations in which individuals move repeatedly among habitat patches to classic metapopulations with infrequent migration among discrete populations. To establish a common framework for population dynamics in patchy habitats, we describe an individual-based model (IBM) involving a diffusion approximation of correlated random walk of individual movements. As an example, we apply the model to the Glanville fritillary butterfly (Melitaea cinxia) inhabiting a highly fragmented landscape. We derive stochastic patch occupancy model (SPOM) approximations for the IBMs assuming pure demographic stochasticity, uncorrelated environmental stochasticity, or completely correlated environmental stochasticity in local dynamics. Using realistic parameter values for the Glanville fritillary, we show that the SPOMs mimic the behavior of the IBMs well. The SPOMs derived from IBMs have parameters that relate directly to the life history and behavior of individuals, which is an advantage for model interpretation and parameter estimation. The modeling approach that we describe here provides a unified framework for patchy populations with much movements among habitat patches and classic metapopulations with infrequent movements.     

Fragmentation reduces regional‐scale spatial genetic structure in a wind‐pollinated tree because genetic barriers are removed

Gene flow strongly influences the regional genetic structuring of plant populations. Seed and pollen dispersal patterns can respond differently to the increased isolation resulting from habitat fragmentation, with unpredictable consequences for gene flow and population structuring. In a recently fragmented landscape we compared the pre‐ and post‐fragmentation genetic structure of populations of a tree species where pollen and seed dispersal respond differentially to forest fragmentation generated by flooding. Castanopsis sclerophylla is wind‐pollinated, with seeds that are dispersed by gravity and rodents. Using microsatellites, we found no significant difference in genetic diversity between pre‐ and post‐fragmentation cohorts. Significant genetic structure was observed in pre‐fragmentation cohorts, due to an unknown genetic barrier that had isolated one small population. Among post‐fragmentation cohorts this genetic barrier had disappeared and genetic structure was significantly weakened. The strengths of genetic structuring were at a similar level in both cohorts, suggesting that overall gene flow of C. sclerophylla has been unchanged by fragmentation at the regional scale. Fragmentation has blocked seed dispersal among habitats, but this appears to have been compensated for by enhanced pollen dispersal, as indicated by the disappearance of a genetic barrier, probably as a result of increased wind speeds and easier pollen movement over water. Extensive pollen flow can counteract some negative effects of fragmentation and assist the long‐term persistence of small remnant populations.     

Effect of habitat fragmentation on dispersal in the butterfly Proclossiana eunomia

Comparison of dispersal rates of the bog fritillary butterfly between continuous and fragmented landscapes indicates that between patch dispersal is significantly lower in the fragmented landscape, while population densities are of the same order of magnitude. Analyses of the dynamics of the suitable habitat for the butterfly in the fragmented landscape reveal a severe, non linear increase in spatial isolation of patches over a time period of 30 years (i.e. 30 butterfly generations), but simulations of the butterfly metapopulation dynamics using a structured population model show that the lower dispersal rates in the fragmented landscape are far above the critical threshold leading to metapopulation extinction. These results indicate that changes in individual behaviour leading to the decrease of dispersal rates in the fragmented landscape were rapidly selected for when patch spatial isolation increased. The evidence of such an adaptive answer to habitat fragmentation suggests that dispersal mortality is a key factor for metapopulation persistence in fragmented landscapes. We emphasise that landscape spatial configuration and patch isolation have to be taken into account in the debate about large-scale conservation strategies.