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

The genetic basis of dispersal in a vertebrate metapopulation

Dispersal affects evolutionary processes by changing population size and genetic composition, influencing the viability and persistence of populations. Investigating which mechanisms underlie variation in dispersal phenotypes and whether populations harbour adaptive potential for dispersal is crucial to understanding the eco-evolutionary dynamics of this important trait. Here, we investigate the genetic architecture of dispersal among successfully recruited individuals in an insular metapopulation of house sparrows. We use an extensive long-term individual-based ecological data set and high-density single-nucleotide polymorphism (SNP) genotypes for over 2500 individuals. We conducted a genome-wide association study (GWAS), and found a relationship between dispersal probability and a SNP located near genes known to regulate circadian rhythm, glycogenesis and exercise performance, among other functions. However, this SNP only explained 3.8% of variance, suggesting that dispersal is a polygenic trait. We then used an animal model to estimate heritable genetic variation (σA 2 ), which composes 10% of the total variation in dispersal probability. Finally, we investigated differences in σA 2 across populations occupying ecologically relevant habitat types (farm vs. non-farm) using a genetic groups animal model. We found different adaptive potentials across habitats, with higher mean breeding value, σA 2 , and heritability for the habitat presenting lower dispersal rates, suggesting also different roles of environmental variation. Our results suggest a complex genetic architecture of dispersal and demonstrate that adaptive potential may be environment dependent in key eco-evolutionary traits. The eco-evolutionary implications of such environment dependence and consequent spatial variation are likely to become ever more important with the increased fragmentation and loss of suitable habitats for many natural populations.

Lack of intergenerational reproductive conflict, rather than lack of inclusive fitness benefits, explains absence of post-reproductive lifespan in long-finned pilot whales

Life-history theory suggests that individuals should reproduce until death, yet females of a small number of mammals live for a significant period after ceasing reproduction, a phenomenon known as post-reproductive lifespan. It is thought that the evolution of this trait is facilitated by increasing local relatedness throughout a female's lifetime. This allows older females to gain inclusive fitness through helping their offspring (known as a mother effect) and/or grandoffspring (known as a grandmother effect), rather than gaining direct fitness through reproducing. However, older females may only benefit from stopping reproducing when their direct offspring compete with those of their daughters. Here, we investigate whether a lack of post-reproductive lifespan in long-finned pilot whales (Globicephala melas) results from minimal benefits incurred from the presence of older females, or from a lack of costs resulting from mother-daughter co-reproduction. Using microsatellite data, we conducted parentage analysis on individuals from 25 pods and find that younger females were more likely to have offspring if their mother was present in their pod, indicating that mothers may assist inexperienced daughters to reproduce. However, we found no evidence of reproductive conflict between co-reproducing mothers and daughters, indicating that females may be able to reproduce into old age while simultaneously aiding their daughters in reproduction. This highlights the importance of reproductive conflict in the evolution of a post-reproductive lifespan and demonstrates that mother and grandmother effects alone do not result in the evolution of a post-reproductive lifespan.

Long-Term Artificial Selection for Increased Larval Body Weight of Hermetia illucens in Industrial Settings

Black soldier fly (Hermetia illucens) farming has exponentially increased in recent years due to the ability of its larvae to efficiently convert low-grade organic materials into high-value food, feed, and technical products. There is a need to further improve the efficiency of production, to meet the rising demands for proteins in the feed and food industries under limited resources. One means of improvement is artificial selection, which has been widely applied in plants and in other livestock species. In 2019, a genetic improvement program was started with the aim to increase larval body weight in black soldier fly larvae. In this paper, we present the outcomes of this breeding program after 10, 13, and 16 generations of selection. The performance of the selected body weight line was compared to the base population line over six experimental rounds under different environmental conditions. Under automated production settings, an average increase of +39% in larval weight, +34% in wet crate yield, +26% in dry matter crate yield, +32% in crude protein per crate, and +21% crude fat per crate was achieved in the selected line compared to the base population line. This research demonstrates the potential contribution of artificial selection to improve efficiency when farming black soldier flies in industrial settings. Further research is needed to fully unlock that potential.

Cold winters have morph-specific effects on natal dispersal distance in a wild raptor

Dispersal is a key process with crucial implications in spatial distribution, density, and genetic structure of species’ populations. Dispersal strategies can vary according to both individual and environmental features, but putative phenotype-by-environment interactions have rarely been accounted for. Melanin-based color polymorphism is a phenotypic trait associated with specific behavioral and physiological profiles and is, therefore, a good candidate trait to study dispersal tactics in different environments. Here, using a 40 years dataset of a population of color polymorphic tawny owls (Strix aluco), we investigated natal dispersal distance of recruiting gray and pheomelanic reddish-brown (hereafter brown) color morphs in relation to post-fledging winter temperature and individual characteristics. Because morphs are differently sensitive to cold winters, we predicted that morphs’ natal dispersal distances vary according to winter conditions. Winter temperature did not affect the proportion of brown (or gray) among recruits. We found that dispersal distances correlate with winter temperature in an opposite manner in the two morphs. Although the gray morph undertakes larger movements in harsher conditions, likely because it copes better with winter severity, the brown morph disperses shorter distances when winters are harsher. We discuss this morph-specific natal dispersal pattern in the context of competition for territories between morphs and in terms of costs and benefits of these alternative strategies. Our results stress the importance of considering the interaction between phenotype and environment to fully disentangle dispersal movement patterns and provide further evidence that climate affects the behavior and local distribution of this species.

Species-specific effects of thermal stress on the expression of genetic variation across a diverse group of plant and animal taxa under experimental conditions

Assessing the genetic adaptive potential of populations and species is essential for better understanding evolutionary processes. However, the expression of genetic variation may depend on environmental conditions, which may speed up or slow down evolutionary responses. Thus, the same selection pressure may lead to different responses. Against this background, we here investigate the effects of thermal stress on genetic variation, mainly under controlled laboratory conditions. We estimated additive genetic variance (V_A), narrow-sense heritability (h²) and the coefficient of genetic variation (CV_A) under both benign control and stressful thermal conditions. We included six species spanning a diverse range of plant and animal taxa, and a total of 25 morphological and life-history traits. Our results show that (1) thermal stress reduced fitness components, (2) the majority of traits showed significant genetic variation and that (3) thermal stress affected the expression of genetic variation (V_A, h² or CV_A) in only one-third of the cases (25 of 75 analyses, mostly in one clonal species). Moreover, the effects were highly species-specific, with genetic variation increasing in 11 and decreasing in 14 cases under stress. Our results hence indicate that thermal stress does not generally affect the expression of genetic variation under laboratory conditions but, nevertheless, increases or decreases genetic variation in specific cases. Consequently, predicting the rate of genetic adaptation might not be generally complicated by environmental variation, but requires a careful case-by-case consideration.

Dispersibility of the Pale Grass Blue Butterfly Zizeeria maha (Lepidoptera: Lycaenidae) Revealed by One-Individual Tracking in the Field: Quantitative Comparisons between Subspecies and between Sexes

The pale grass blue butterfly Zizeeria maha (Lepidoptera: Lycaenidae) has been used as an environmental indicator species for radioactive pollution after the Fukushima nuclear accident. Here, based on the one-individual tracking method in the field, we examined dispersal-associated and other behavioral traits of this butterfly, focusing on two subspecies, Z. maha argia in mainland Japan and Z. maha okinawana in Okinawa. The accumulated distances in the adult lifespan were 18.9 km and 38.2 km in mainland and Okinawa males, respectively, and 15.0 km and 7.8 km in mainland and Okinawa females, respectively. However, the mean distance from the starting point was only 24.2 m and 21.1 m in the mainland and Okinawa males, respectively, and 13.7 m and 7.4 m in the mainland and Okinawa females, respectively. Some quantitative differences in resting and feeding were found between subspecies and between sexes. The ARIMA (autoregressive integrated moving average) model indicated that the dispersal distance was 52.3 m (99% confidence interval value of 706.6 m) from the starting point in mainland males. These results support the idea that despite some behavioral differences, both subspecies of this butterfly are suitable as an environmental indicator because of the small dispersal ranges.

Landscape permeability and individual variation in a dispersal-linked gene jointly determine genetic structure in the Glanville fritillary butterfly

There is now clear evidence that species across a broad range of taxa harbor extensive heritable variation in dispersal. While studies suggest that this variation can facilitate demographic outcomes such as range expansion and invasions, few have considered the consequences of intraspecific variation in dispersal for the maintenance and distribution of genetic variation across fragmented landscapes. Here, we examine how landscape characteristics and individual variation in dispersal combine to predict genetic structure using genomic and spatial data from the Glanville fritillary butterfly. We used linear and latent factor mixed models to identify the landscape features that best predict spatial sorting of alleles in the dispersal-related gene phosphoglucose isomerase (Pgi). We next used structural equation modeling to test if variation in Pgi mediated gene flow as measured by Fst at putatively neutral loci. In a year when the population was recovering following a large decline, individuals with a genotype associated with greater dispersal ability were found at significantly higher frequencies in populations isolated by water and forest, and these populations showed lower levels of genetic differentiation at neutral loci. These relationships disappeared in the next year when metapopulation density was high, suggesting that the effects of individual variation are context dependent. Together our results highlight that (1) more complex aspects of landscape structure beyond just the configuration of habitat can be important for maintaining spatial variation in dispersal traits and (2) that individual variation in dispersal plays a key role in maintaining genetic variation across fragmented landscapes.

Impact of male condition on his spermatophore and consequences for female reproductive performance in the Glanville fritillary butterfly

In butterflies, male reproductive success is highly related to the quality and the size of the spermatophore transferred to the female. The spermatophore is a capsule produced by the male during copulation, which in many species contains sperm in addition to a nuptial gift, and which is digested by the female after copulation. The nuptial gift may contribute to egg production and offspring quality, and in some cases also to female body maintenance. The production of the spermatophore, however, represents a cost for the male and, in polyandrous species, ejaculates are sometimes allocated adaptively across matings. Nonetheless, although the ecological factors affecting the reproductive success of female butterflies have been the topic of numerous studies, little information exists on the factors affecting males' contribution to reproduction, and the indirect impacts on female fecundity and fitness. We used the Glanville fritillary butterfly, Melitaea cinxia (Linnaeus, 1758) (Nymphalidae), in order to assess variation in male allocation to matings. In this species, smaller males produce smaller spermatophores, but variation in spermatophore size is not correlated with female reproductive success. We show that spermatophore size increases with male age at first mating, decreases with mating frequency and adult food-deprivation, and is not influenced by developmental food-limitation. The length of copulation period does not influence the spermatophore size nor influences the polyandrous mating behavior in this species. Male contribution to his spermatophore size is clearly influenced by his condition and adult-resource at the time of mating. Despite this variation, spermatophore size does not seem to have a direct impact on female reproductive output or mating behavior.

Is the behavioural divergence between range-core and range-edge populations of cane toads (Rhinella marina) due to evolutionary change or developmental plasticity?

Individuals at the leading edge of expanding biological invasions often show distinctive phenotypic traits, in ways that enhance their ability to disperse rapidly and to function effectively in novel environments. Cane toads (Rhinella marina) at the invasion front in Australia exhibit shifts in morphology, physiology and behaviour (directionality of dispersal, boldness, risk-taking). We took a common-garden approach, raising toads from range-core and range-edge populations in captivity, to see if the behavioural divergences observed in wild-caught toads are also evident in common-garden offspring. Captive-raised toads from the invasion vanguard population were more exploratory and bolder (more prone to 'risky' behaviours) than toads from the range core, which suggests that these are evolved, genetic traits. Our study highlights the importance of behaviour as being potentially adaptive in invasive populations and adds these behavioural traits to the increasing list of phenotypic traits that have evolved rapidly during the toads' 80-year spread through tropical Australia.

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.

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.

The evolution of the competition-dispersal trade-off affects α- and β-diversity in a heterogeneous metacommunity

Difference in dispersal ability is a key driver of species coexistence in metacommunities. However, the available frameworks for interpreting species diversity patterns in natura often overlook trade-offs and evolutionary constraints associated with dispersal. Here, we build a metacommunity model accounting for dispersal evolution and a competition-dispersal trade-off. Depending on the distribution of carrying capacities among communities, species dispersal values are distributed either around a single strategy (evolutionarily stable strategy, ESS), or around distinct strategies (evolutionary branching, EB). We show that limited dispersal generates spatial aggregation of dispersal traits in ESS and EB scenarios, and that the competition-dispersal trade-off strengthens the pattern in the EB scenario. Importantly, individuals in larger (respectively (resp.) smaller) communities tend to harbour lower (resp. higher) dispersal, especially under the EB scenario. We explore how dispersal evolution affects species diversity patterns by comparing those from our model to the predictions of a neutral metacommunity model. The most marked difference is detected under EB, with distinctive values of both α- and β-diversity (e.g. the dissimilarity in species composition between small and large communities was significantly larger than neutral predictions). We conclude that, from an empirical perspective, jointly assessing community carrying capacity with species dispersal strategies should improve our understanding of diversity patterns in metacommunities.

Rapid evolution of dispersal ability makes biological invasions faster and more variable

Genetic variation in dispersal ability may result in the spatial sorting of alleles during range expansion. Recent theory suggests that spatial sorting can favour the rapid evolution of life history traits at expanding fronts, and therefore modify the ecological dynamics of range expansion. Here we test this prediction by disrupting spatial sorting in replicated invasions of the bean beetle Callosobruchus maculatus across homogeneous experimental landscapes. We show that spatial sorting promotes rapid evolution of dispersal distance, which increases the speed and variability of replicated invasions: after 10 generations of range expansion, invasions subject to spatial sorting spread 8.9% farther and exhibit 41-fold more variable spread dynamics relative to invasions in which spatial sorting is suppressed. Correspondingly, descendants from spatially evolving invasions exhibit greater mean and variance in dispersal distance. Our results reveal an important role for rapid evolution during invasion, even in the absence of environmental filters, and argue for evolutionarily informed forecasts of invasive spread by exotic species or climate change migration by native species.

Theory suggests that spatial sorting by dispersal ability can generate evolutionarily accelerated range expansions. Using the bean beetle Callosobruchus maculatus, this study shows that evolution not only increases the speed of range expansion, as predicted, but also increases variability.

Individual variation changes dispersal distance and area requirements of a checkerspot butterfly

Individual variation in movement can have important consequences for spatial population dynamics. For instance, individual variation increases leptokurtosis in dispersal distance, such that more individuals move very short and very long distances relative to a homogeneous population. We quantified individual variation in movement of the Baltimore checkerspot butterfly (Euphydryas phaeton) to investigate its importance for two conservation-related metrics: the expected dispersal distance and the critical minimum patch size, or the smallest area within which a population can persist based on loss due to emigration. All movement parameters showed among-individual variation, with the greatest variation in move lengths and time spent resting. Correlations in among-individual movement parameters indicated that some butterflies were generally more mobile than others. We incorporated empirically estimated movement and demographic parameters into two individual-based models (IBMs), one with homogeneity in movement among individuals, and one with heterogeneity in movement. As expected, individual variation in movement increased the leptokurtosis of lifetime movement distance; the maximum difference in distance moved was substantial (-850 m vs. -5800 m) and is likely of significance for conservation. Individual variation also affected the critical minimum patch size, but the difference (-0.7 ha vs. -0.5 ha) is unlikely to be ecologically significant. Notably, populations with individual variation had higher growth rates in small patches and lower growth rates in large patches, a logical consequence of increased leptokurtosis. Individual variation in movement is fairly straightforward to quantify using mixed effects models and is important for spatial population dynamics, thus we encourage its inclusion in studies of other systems.

Temperature- and sex-related effects of serine protease alleles on larval development in the Glanville fritillary butterfly

The body reserves of adult Lepidoptera are accumulated during larval development. In the Glanville fritillary butterfly, larger body size increases female fecundity, but in males fast larval development and early eclosion, rather than large body size, increase mating success and hence fitness. Larval growth rate is highly heritable, but genetic variation associated with larval development is largely unknown. By comparing the Glanville fritillary population living in the Åland Islands in northern Europe with a population in Nantaizi in China, within the source of the post-glacial range expansion, we identified candidate genes with reduced variation in Åland, potentially affected by selection under cooler climatic conditions than in Nantaizi. We conducted an association study of larval growth traits by genotyping the extremes of phenotypic trait distributions for 23 SNPs in 10 genes. Three genes in clip-domain serine protease family were associated with larval growth rate, development time and pupal weight. Additive effects of two SNPs in the prophenoloxidase-activating proteinase-3 (ProPO3) gene, related to melanization, showed elevated growth rate in high temperature but reduced growth rate in moderate temperature. The allelic effects of the vitellin-degrading protease precursor gene on development time were opposite in the two sexes, one genotype being associated with long development time and heavy larvae in females but short development time in males. Sexually antagonistic selection is here evident in spite of sexual size dimorphism.

The roles of demography and genetics in the early stages of colonization

Colonization success increases with the size of the founding group. Both demographic and genetic factors underlie this relationship, yet because genetic diversity normally increases with numbers of individuals, their relative importance remains unclear. Furthermore, their influence may depend on the environment and may change as colonization progresses from establishment through population growth and then dispersal. We tested the roles of genetics, demography and environment in the founding of Tribolium castaneum populations. Using three genetic backgrounds (inbred to outbred), we released individuals of four founding sizes (2-32) into two environments (natal and novel), and measured establishment success, initial population growth and dispersal. Establishment increased with founding size, whereas population growth was shaped by founding size, genetic background and environment. Population growth was depressed by inbreeding at small founding sizes, but growth rates were similar across genetic backgrounds at large founding size, an interaction indicating that the magnitude of the genetic effects depends upon founding population size. Dispersal rates increased with genetic diversity. These results suggest that numbers of individuals may drive initial establishment, but that subsequent population growth and spread, even in the first generation of colonization, can be driven by genetic processes, including both reduced growth owing to inbreeding depression, and increased dispersal with increased genetic diversity.

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.

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.

Influence of developmental conditions on immune function and dispersal-related traits in the Glanville fritillary (Melitaea cinxia) butterfly

Organisms in the wild are constantly faced with a wide range of environmental variability, such as fluctuation in food availability. Poor nutritional conditions influence life-histories via individual resource allocation patterns, and trade-offs between competing traits. In this study, we assessed the influence of food restriction during development on the energetically expensive traits flight metabolic rate (proxy of dispersal ability), encapsulation rate (proxy of immune defence), and lifespan using the Glanville fritillary butterfly, Melitaea cinxia, as a model organism. Additionally, we examined the direct costs of flight on individual immune function, and whether those costs increase under restricted environmental conditions. We found that nutritional restriction during development enhanced adult encapsulations rate, but reduced both resting and flight metabolic rates. However, at the individual level metabolic rates were not associated with encapsulation rate. Interestingly, individuals that were forced to fly prior to the immune assays had higher encapsulation rates than individuals that had not flown, suggesting that flying itself enhances immune response. Finally, in the control group encapsulation rate correlated positively with lifespan, whereas in the nutritional restriction group there was no relationship between these traits, suggesting that the association between encapsulation rate on adult lifespan was condition-dependent. Thus stressful events during both larval development (food limitation) and adulthood (forced flight) induce increased immune response in the adult butterflies, which may allow individuals to cope with stressful events later on in life.

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.

Plastic larval development in a butterfly has complex environmental and genetic causes and consequences for population dynamics

1. In insects, the length of larval development time typically influences adult body size and individual fitness, and hence development time can be expected to respond in an adaptive manner to variation in environmental conditions. In the wild, larval growth may be influenced by individual condition, which can be affected by population-level parameters such as population density and abundance and quality of resources. 2. We sampled larvae of the Glanville fritillary butterfly (Melitaea cinxia) from 514 local populations across a large metapopulation before the winter diapause and reared the larvae in common garden conditions after diapause. Here, we report that small post-diapause larvae prolonged their development via an extra larval instar, apparently to compensate for their 'bad start' after diapause. The number of instars was additionally a plastic response to environmental conditions, as the frequency of the extra instar increased under cooler thermal conditions. 3. The benefit of the extra instar is clear, as it allows individuals to develop into larger adults, but the cost is delayed adult eclosion, which is likely to select against the extra instar especially in males, in which early eclosion is critical for mating success. In support of this, the frequency of the extra instar was significantly lower in males (7%) than in females (42%). 4. Polymorphisms in three genes, serpin-1, vitellin-degrading protease precursor and phosphoglucose isomerase, which are known to influence development in insects, were associated with the occurrence of the extra instar. 5. At the level of local populations, the frequency of the extra instar was higher in newly established populations than that in old local ones, possibly reflecting maternal effects, as new populations are often established by females with heavy investment in dispersal. The frequency of the extra instar in turn correlated with the change in population size over 1 year and the risk of local extinction in the natural metapopulation of the Glanville fritillary. 6. Our results highlight the importance of the physiological condition of individuals in shaping subsequent life-history events and even population dynamics.

Environmentally induced dispersal-related life-history syndrome in the tropical butterfly, Bicyclus anynana

Dispersal is a key process for understanding the persistence of populations as well as the capacity of organisms to respond to environmental change. Therefore, understanding factors that may facilitate or constrain the evolution of dispersal is of crucial interest. Assessments of phenotypic variation in various behavioural, physiological and morphological traits related to insect dispersal and flight performance are common, yet very little is known about the genetic associations among these traits. We have used experiments on the butterfly Bicyclus anynana to estimate genetic variation and covariation in seven behavioural, physiological and morphological traits related to flight potential and hence dispersal. Our goal was to characterize the heritabilities and genetic correlations among these traits and thus to understand more about the evolution of dispersal-related life-history syndromes in butterflies. Using a version of the animal model, we showed that all of the traits varied between the sexes, and most were either positively or negatively (phenotypically and/or genetically) correlated with body size. Heritable variation was present in most traits, with the highest heritabilities estimated for body mass and thorax ratio. The variance in flight activity among multiple measurements for the same individual was high even after controlling for the prevailing environmental conditions, indicating the importance of behavioural switching and/or inherent randomness associated with this type of movement. A number of dispersal-related traits showed phenotypic correlations among one another, but only a few of these were associated with significant genetic correlations indicating that covariances between these traits in Bicyclus anynana are mainly environmentally induced.

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.

Flight metabolic rate has contrasting effects on dispersal in the two sexes of the Glanville fritillary butterfly

Evolution of dispersal is affected by context-specific costs and benefits. One example is sex-biased dispersal in mammals and birds. While many such patterns have been described, the underlying mechanisms are poorly understood. Here, we study genetic and phenotypic traits that affect butterfly flight capacity and examine how these traits are related to dispersal in male and female Glanville fritillary butterflies (Melitaea cinxia). We performed two mark-recapture experiments to examine the associations of individuals' peak flight metabolic rate (MR(peak)) and Pgi genotype with their dispersal in the field. In a third experiment, we studied tethered flight in the laboratory. MR(peak) was negatively correlated with dispersal distance in males but the trend was positive in females, and the interaction between MR(peak) and sex was significant for long-distance dispersal. A similar but nonsignificant trend was found in relation to molecular variation at Pgi, which encodes a glycolytic enzyme: the genotype associated with high MR(peak) tended to be less dispersive in males but more dispersive in females. The same pattern was repeated in the tethered flight experiment: the relationship between MR(peak) and flight duration was positive in females but negative in males. These results suggest that females with high flight capacity are superior in among-population dispersal, which facilitates the spatial spreading of their reproductive effort. In contrast, males with high flight capacity may express territorial behaviour, and thereby increase the number of matings, whereas inferior males may be forced to disperse. Thus, flight capacity has opposite associations with dispersal rate in the two sexes.

Behavioural responses to human-induced environmental change

The initial response of individuals to human-induced environmental change is often behavioural. This can improve the performance of individuals under sudden, large-scale perturbations and maintain viable populations. The response can also give additional time for genetic changes to arise and, hence, facilitate adaptation to new conditions. On the other hand, maladaptive responses, which reduce individual fitness, may occur when individuals encounter conditions that the population has not experienced during its evolutionary history, which can decrease population viability. A growing number of studies find human disturbances to induce behavioural responses, both directly and by altering factors that influence fitness. Common causes of behavioural responses are changes in the transmission of information, the concentration of endocrine disrupters, the availability of resources, the possibility of dispersal, and the abundance of interacting species. Frequent responses are alterations in habitat choice, movements, foraging, social behaviour and reproductive behaviour. Behavioural responses depend on the genetically determined reaction norm of the individuals, which evolves over generations. Populations first respond with individual behavioural plasticity, whereafter changes may arise through innovations and the social transmission of behavioural patterns within and across generations, and, finally, by evolution of the behavioural response over generations. Only a restricted number of species show behavioural adaptations that make them thrive in severely disturbed environments. Hence, rapid human-induced disturbances often decrease the diversity of native species, while facilitating the spread of invasive species with highly plastic behaviours. Consequently, behavioural responses to human-induced environmental change can have profound effects on the distribution, adaptation, speciation and extinction of populations and, hence, on biodiversity. A better understanding of the mechanisms of behavioural responses and their causes and consequences could improve our ability to predict the effects of human-induced environmental change on individual species and on biodiversity.

Locomotion in response to shifting climate zones: not so fast

Although a species' locomotor capacity is suggestive of its ability to escape global climate change, such a suggestion is not necessarily straightforward. Species vary substantially in locomotor capacity, both ontogenetically and within/among populations, and much of this variation has a genetic basis. Accordingly, locomotor capacity can and does evolve rapidly, as selection experiments demonstrate. Importantly, even though this evolution of locomotor capacity may be rapid enough to escape changing climate, genetic correlations among traits (often due to pleiotropy) are such that successful or rapid dispersers are often limited in colonization or reproductive ability, which may be viewed as a trade-off. The nuanced assessment of this variation and evolution is reviewed for well-studied models: salmon, flying versus flightless insects, rodents undergoing experimental evolution, and metapopulations of butterflies. This work reveals how integration of physiology with population biology and functional genomics can be especially informative.

A meta-analysis of dispersal in butterflies

Dispersal has recently gained much attention because of its crucial role in the conservation and evolution of species facing major environmental changes such as habitat loss and fragmentation, climate change, and their interactions. Butterflies have long been recognized as ideal model systems for the study of dispersal and a huge amount of data on their ability to disperse has been collected under various conditions. However, no single 'best' method seems to exist leading to the co-occurrence of various approaches to study butterfly mobility, and therefore a high heterogeneity among data on dispersal across this group. Accordingly, we here reviewed the knowledge accumulated on dispersal and mobility in butterflies, to detect general patterns. This meta-analysis specifically addressed two questions. Firstly, do the various methods provide a congruent picture of how dispersal ability is distributed across species? Secondly, is dispersal species-specific? Five sources of data were analysed: multisite mark-recapture experiments, genetic studies, experimental assessments, expert opinions, and transect surveys. We accounted for potential biases due to variation in genetic markers, sample sizes, spatial scales or the level of habitat fragmentation. We showed that the various dispersal estimates generally converged, and that the relative dispersal ability of species could reliably be predicted from their relative vagrancy (records of butterflies outside their normal habitat). Expert opinions gave much less reliable estimates of realized dispersal but instead reflected migration propensity of butterflies. Within-species comparisons showed that genetic estimates were relatively invariable, while other dispersal estimates were highly variable. This latter point questions dispersal as a species-specific, invariant trait.

Heritability of and strong single gene (Pgi) effects on life-history traits in the Glanville fritillary butterfly

We estimated broad-sense heritabilities (H(2)) of 13 female and seven male life-history traits of the Glanville fritillary butterfly (Melitaea cinxia) under semi-natural conditions in a large outdoor population cage. The analysis was based on full-sib families collected as young larvae in the field and reared under common garden conditions. We found significant genetic variance in female lifespan, fecundity, number of matings and host-plant preference as well as in male body mass and mobility. Apart from host-plant preference, female traits that were more strongly correlated with lifetime reproductive success (LRS; measured as total number of eggs laid) had higher H(2). LRS itself exhibited significant heritability. Host-plant preference had very high H(2), consistent with a previously reported genetically determined geographical cline in host-plant preference in the study area. Lifespan and egg hatching rate were significantly associated with a SNP in the coding region of the Pgi gene, for which there is previous evidence for balancing selection. Selection on Pgi, which furthermore shows spatial and temporal variation, may maintain genetic variance in fitness-related life-history traits. In contrast, we found no strong evidence for life-history trade-offs.

Modelling single nucleotide effects in phosphoglucose isomerase on dispersal in the Glanville fritillary butterfly: coupling of ecological and evolutionary dynamics

Dispersal comprises a complex life-history syndrome that influences the demographic dynamics of especially those species that live in fragmented landscapes, the structure of which may in turn be expected to impose selection on dispersal. We have constructed an individual-based evolutionary sexual model of dispersal for species occurring as metapopulations in habitat patch networks. The model assumes correlated random walk dispersal with edge-mediated behaviour (habitat selection) and spatially correlated stochastic local dynamics. The model is parametrized with extensive data for the Glanville fritillary butterfly. Based on empirical results for a single nucleotide polymorphism (SNP) in the phosphoglucose isomerase (Pgi) gene, we assume that dispersal rate in the landscape matrix, fecundity and survival are affected by a locus with two alleles, A and C, individuals with the C allele being more mobile. The model was successfully tested with two independent empirical datasets on spatial variation in Pgi allele frequency. First, at the level of local populations, the frequency of the C allele is the highest in newly established isolated populations and the lowest in old isolated populations. Second, at the level of sub-networks with dissimilar numbers and connectivities of patches, the frequency of C increases with decreasing network size and hence with decreasing average metapopulation size. The frequency of C is the highest in landscapes where local extinction risk is high and where there are abundant opportunities to establish new populations. Our results indicate that the strength of the coupling of the ecological and evolutionary dynamics depends on the spatial scale and is asymmetric, demographic dynamics having a greater immediate impact on genetic dynamics than vice versa.

Weight and nutrition affect pre-mRNA splicing of a muscle gene associated with performance, energetics and life history

A fundamental feature of gene expression in multicellular organisms is the production of distinct transcripts from single genes by alternative splicing (AS), which amplifies protein and functional diversity. In spite of the likely consequences for organismal biology, little is known about how AS varies among individuals or responds to body condition, environmental variation or extracellular signals in general. Here we show that evolutionarily conserved AS of troponin-t in flight muscle of adult moths responds in a quantitative fashion to experimental manipulation of larval nutrition and adult body weight. Troponin-t (Tnt) isoform composition is known to affect muscle force and power output in other animals, and is shown here to be associated with the thorax mass-specific rate of energy consumption during flight. Loading of adults with external weights for 5 days caused an AS response nearly identical to equal increases in actual body weight. In addition, there were effects of larval feeding history on adult Tnt isoform composition that were independent of body weight, with moths from poorer larval feeding regimes producing isoform profiles associated with reduced muscle performance and energy consumption rate. Thus, Tnt isoform composition in striated muscle is responsive to both weight-sensing and nutrition-sensing mechanisms, with consequent effects on function. In free-living butterflies, Tnt isoform composition was also associated with activity level and very strongly with the rate of egg production. Overall, these results show that AS of a muscle gene responds in a quantitative fashion to whole-organism variables, which apparently serves to coordinate muscle strength and energy expenditure with body condition and life history.

Significant effects of Pgi genotype and body reserves on lifespan in the Glanville fritillary butterfly

Individuals with a particular variant of the gene phosphoglucose isomerase (Pgi) have been shown to have superior dispersal capacity and fecundity in the Glanville fritillary butterfly (Melitaea cinxia), raising questions about the mechanisms that maintain polymorphism in this gene in the field. Here, we investigate how variation in the Pgi genotype affects female and male life history under controlled conditions. The most striking effect is the longer lifespan of genotypes with high dispersal capacity, especially in non-reproducing females. Butterflies use body reserves for somatic maintenance and reproduction, but different resources (in thorax versus abdomen) are used under dissimilar conditions, with some interactions with the Pgi genotype. These results indicate life-history trade-offs that involve resource allocation and genotypexenvironment interactions, and these trade-offs are likely to contribute to the maintenance of Pgi polymorphism in the natural populations.

Tracking butterfly movements with harmonic radar reveals an effect of population age on movement distance

We used harmonic radar to track freely flying Glanville fritillary butterfly (Melitaea cinxia) females within an area of 30 ha. Butterflies originated from large and continuous populations in China and Estonia, and from newly established or old (> 5 years) small local populations in a highly fragmented landscape in Finland. Caterpillars were raised under common garden conditions and unmated females were tested soon after eclosion. The reconstructed flight paths for 66 individuals comprised a total distance of 51 km with high spatial resolution. Butterflies originating from large continuous populations and from old local populations in Finland exhibited similar movement behaviors, whereas butterflies originating from newly established local populations in the fragmented landscape in Finland moved significantly more than the others. There was no difference in the lengths of individual flight bouts, but the new-population females flew more frequently, resulting in longer daily movement tracks. The flight activity of all individuals was affected by environmental conditions, peaking at 19-23 degrees C (depending on population type), in the early afternoon, and during calm weather. Butterflies from all population types showed a strong tendency to follow habitat edges between the open study area and the neighboring woodlands.

Rapidly developing functional genomics in ecological model systems via 454 transcriptome sequencing

Next generation sequencing technology affords new opportunities in ecological genetics. This paper addresses how an ecological genetics research program focused on a phenotype of interest can quickly move from no genetic resources to having various functional genomic tools. 454 sequencing and its error rates are discussed, followed by a review of de novo transcriptome assemblies focused on the first successful de novo assembly which happens to be in an ecological model system (the Glanville fritillary butterfly). The potential future developments in 454 sequencing are also covered. Particular attention is paid to the difficulties ecological geneticists are likely to encounter through reviewing relevant studies in both model and non-model systems. Various post-sequencing issues and applications of 454 generated data are presented (e.g. database management, microarray construction, molecular marker and candidate gene development). How to use species with genomic resources to inform study of those without is also discussed. In closing, some of the drawbacks of 454 sequencing are presented along with future prospects of this technology.

Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing

We present a de novo assembly of a eukaryote transcriptome using 454 pyrosequencing data. The Glanville fritillary butterfly (Melitaea cinxia; Lepidoptera: Nymphalidae) is a prominent species in population biology but had no previous genomic data. Sequencing runs using two normalized complementary DNA collections from a genetically diverse pool of larvae, pupae, and adults yielded 608,053 expressed sequence tags (mean length = 110 nucleotides), which assembled into 48,354 contigs (sets of overlapping DNA segments) and 59,943 singletons. BLAST comparisons confirmed the accuracy of the sequencing and assembly, and indicated the presence of c. 9000 unique genes, along with > 6000 additional microarray-confirmed unannotated contigs. Average depth of coverage was 6.5-fold for the longest 4800 contigs (348-2849 bp in length), sufficient for detecting large numbers of single nucleotide polymorphisms. Oligonucleotide microarray probes designed from the assembled sequences showed highly repeatable hybridization intensity and revealed biological differences among individuals. We conclude that 454 sequencing, when performed to provide sufficient coverage depth, allows de novo transcriptome assembly and a fast, cost-effective, and reliable method for development of functional genomic tools for nonmodel species. This development narrows the gap between approaches based on model organisms with rich genetic resources vs. species that are most tractable for ecological and evolutionary studies.