|Home | About | Journals | Submit | Contact Us | Français|
Frequent extinctions of local populations in metapopulations create opportunities for migrant females to establish new populations. In a metapopulation of the Glanville fritillary butterfly (Melitaea cinxia), more mobile individuals are more likely to establish new populations, especially in habitat patches that are poorly connected to existing populations. Here we show that flight metabolic rate and the frequency of a specific allele of the metabolic enzyme phosphoglucose isomerase (pgi) were both highest in newly established, isolated populations. Furthermore, genotypes with this pgi allele had elevated flight metabolic rates. These results suggest that genetic variation in pgi or a closely linked locus has a direct effect on flight metabolism, dispersal rate, and thereby on metapopulation dynamics in this species. These results also contribute to an emerging understanding of the mechanisms by which population turnover in heterogeneous landscapes may maintain genetic and phenotypic variation across populations.
Dispersal is one of the key traits determining survival of species in fragmented landscapes. There is an extensive literature examining both theoretically and empirically how changes in landscape structure may influence the evolution of dispersal rate (Clobert et al. 2001; Woiwod et al. 2001; Bullock et al. 2002; Clobert et al. 2004). This work has shown that there is often considerable additive genetic variance for dispersal, and hence dispersal-related traits may respond rapidly to selection (reviewed in Roff & Fairbairn 2001).
The structure of a fragmented landscape may itself contribute to the maintenance of variation in dispersal rate (Travis & Dytham 1999; Hovestadt et al. 2001; Hanski et al. 2004), because dispersive individuals are more likely than sedentary ones to establish new populations and less likely to persist in old populations (Peroni 1994; Olivieri et al. 1995; Taylor & Merriam 1995; Cody & Overton 1996; Piquot et al. 1998; Thomas et al. 1998; Hill et al. 2001; Hanski et al. 2002; Hughes et al. 2003; Hanski et al. 2004). It is clear that variation in dispersal rate has a physiological and genetic basis, but the actual mechanisms remain poorly understood (Roff & Fairbairn 2001). This is especially true for species with continuous variation in dispersal capacity. For wing-dimorphic insects, some of the physiological mechanisms are better known (Zera & Denno 1997; Zhao & Zera 2002) and in some cases wing polymorphism appears to be caused by a single gene with two alleles of unknown molecular function (Roff & Fairbairn 1991; Caillaud et al. 2002; Braendle et al. 2005). However, so far no gene of known molecular function has been found to influence dispersal rate.
Here we report on a candidate gene (Fitzpatrick et al. 2005) and related physiological measurements that are strongly correlated with variation in dispersal rate in the Glanville fritillary butterfly (Melitaea cinxia). Specifically, we relate genetic variation in the glycolytic enzyme phosphoglucose isomerase (pgi) and variation in flight metabolic rate to known spatial variation in dispersal rate in a large and well studied metapopulation of the Glanville fritillary in the Åland Islands in Finland (Hanski et al. 2002; Hanski et al. 2004). This metapopulation occurs in a fragmented landscape consisting of ca 4000 discrete habitat patches, scattered across an area of 50×70km (Hanski et al. 1995; Hanski 1999; Nieminen et al. 2004). Approximately 500 of the patches are occupied in any given year. The rate of population turnover is high, with ca 100 local extinctions each year and a roughly equivalent number of colonizations of unoccupied patches (Nieminen et al. 2004). A typical female visits at most a few other patches apart from the natal patch during her lifetime (Hanski et al. 1994; Kuussaari et al. 1996; Hanski et al. 2004).
Using a spatially realistic model of the evolution of dispersal rate in this metapopulation, Hanski et al. (2004) predicted that females that establish new local populations have a higher dispersal rate than the average female in the metapopulation. Furthermore, this difference should increase with patch isolation, because higher dispersal capacity makes it more likely that an individual will reach a patch with low connectivity to existing populations. On the other hand, with time, the average mobility of individuals in isolated populations should decrease, because highly mobile individuals are likely to emigrate and, due to patch isolation, emigration is not balanced by the arrival of other mobile individuals (Hanski et al. 2004). The model thus predicted a contrasting relationship between average mobility and spatial connectivity in old versus new populations (figure 1). Results consistent with these predictions were found in a mark–recapture study that examined dispersal rate in the field (Hanski et al. 2002) and in a physiological study that measured the [ATP]/[ADP] ratio in flight muscles of female butterflies after a fixed period of forced flight (Hanski et al. 2004). The latter result suggested that ability to rephosphorylate ADP correlates with variation in dispersal rate, leading to the hypothesis that flight metabolism has an immediate effect on variation in dispersal rate.
Flying insects have the highest known mass-specific rates of energy consumption, with certain glycolytic enzymes working at rates close to their maximal flux capacity (Suarez 2000). Nonetheless, genetic variation for flight capacity appears to be rampant, as numerous studies have demonstrated heritable variation for flight endurance (Kent & Rankin 2001; Roff & Fairbairn 2001). Maximal flight performance also shows substantial genetic variation that responds to directional selection (Marden et al. 1997) and maps to specific genomic regions (Montooth et al. 2003). It is reasonable to hypothesize that genes responsible for variation in flight metabolic rate may also be responsible for variation in dispersal rate. Thus, in the present study, we link allelic variation at a metabolic enzyme in the Glanville fritillary to variation in flight physiology among individuals that originate from populations of different ages and spatial connectivities.
We identified pgi as a candidate locus for variation in flight metabolic rate and dispersal rate based on previous studies showing that, in Colias butterflies, flight capacity and female fecundity are correlated with genotypic differences in pgi enzyme kinetics and thermal stability (Watt 1977; Watt 1983; Watt 1992; Hughes & Zalucki 1993). The Åland metapopulation of the Glanville fritillary is known to be polymorphic for pgi (Saccheri et al. 1998). We tested whether female butterflies with different pgi genotypes differ in flight metabolic rate, and whether pgi alleles and flight metabolic rate vary among populations of dissimilar ages and spatial connectivities in a manner predicted by the model of Hanski et al. (2004) and consistent with observed phenotypic variation in dispersal rate.
A total of 176 adult butterflies were caught in the field during the early flight season in 2004, between 2 and 12 June. Half of the individuals were caught from populations that were newly established by migrants in the previous summer, and half from populations older than 5 years. The spatial connectivity of populations varied widely in both groups of populations. Connectivity of population i was calculated as
where dij is the distance in km between populations i and j and Nj is the number of larval groups in population j in autumn 2002. Connectivity, the inverse of isolation, is expected to be proportional to the number of immigrants arriving in a population (for further details see Hanski et al. 2004). Captured butterflies were individually marked and released into a large outdoor cage (32×26×3m), where the vast majority of all matings and ovipositions were recorded under practically natural conditions (described in more detail in Hanski et al. in press). It was not possible to determine the age of the individuals transferred to the cage, but they were all caught in the early flight season and thus were at most a few days old. A more detailed description of the long-term field study and methods used to survey the populations is presented by Nieminen et al. (2004).
From 10 to 16 June, randomly selected females were removed from the cage for one day to measure their flight metabolic rate. Flight metabolic rate was measured using standard respirometry techniques (Lighton 1991) that are described in detail in §S2 in the electronic supplementary material. We stimulated butterflies to fly inside a transparent 1l jar through which dry CO2-free air was pumped at a regulated flow rate of 0.95lmin−1. Air temperature within the jar was nearly invariant during the testing of an individual butterfly, averaging 31.3°C across all of the tests (s.d.=0.7). Butterflies flew in a nearly continuous fashion as long as we applied gentle shakes or taps to get them flying again whenever they alighted. Flight was stimulated in this fashion for 15min, after which the jar was shaded and a steady baseline of resting CO2 emission was reestablished. Many individuals fatigued prior to the end of the 15min flight test; these were clearly cases of an inability to continue to fly rather than unwillingness. In these cases we continued stimulation so that the butterfly flew again as soon as it was able. At the conclusion of the experiment, the butterfly was removed from the jar, allowed to imbibe 25% honey solution, and returned to the cage. Respirometry experiments were performed blindly with regard to the source population of the butterflies. From the data recorded, we subtracted the mean pre-flight CO2 emission rate (resting metabolism) and used standard equations for open-flow respirometry (Lighton 1991) to determine the rate of CO2 emission attributable to flight metabolism.
In separate experiments we have established that the respiratory exchange quotient (the ratio of CO2 emitted to O2 consumed) does not differ from 1.0 over the 15min flight test or between strongly and weakly flying individuals. Other experiments showed that flight metabolic rate in this species remains relatively constant with age until it drops precipitously at the oldest ages (K. Niitepõld, unpublished data).
During 21–24 June, we collected all remaining individuals in the cage and preserved them in liquid nitrogen. Individuals that had been found newly dead or dying before 21 June had been preserved earlier. A total of 75 individuals (32 females, 43 males) were thus preserved and later genotyped for pgi. These 75 individuals originated from 18 new and 21 old populations (one to seven individuals per population), scattered across the Åland Islands (figure S1 in the electronic supplementary material). The material is thus well replicated at the population level.
The sample of 75 individuals was genotyped at the phosphoglucose-isomerase locus (pgi, enzyme commission number EC 126.96.36.199). Half of the thoraces were homogenized in 200μl water. Allelic variation was assessed using 6μl homogenate for cellulose acetate electrophoresis (Hebert & Beaton 1993). Gels were run in Tris Glycine buffer at 200V for 2h in a refrigerator and scored directly after staining.
Individuals from new populations were the offspring of the original colonizers rather than the colonizers themselves and thus the genotype of the colonizers remained unknown. Hence, for data analysis, we contrasted allele frequencies rather than genotype frequencies between old and new populations. This contrast is slightly conservative, because the allele frequencies represent equally the colonizing females and their mates. Females typically mate soon after eclosion and hence before dispersal (Kuussaari et al. 1998).
All egg clusters laid by individual females in the outdoor cage were collected and the eggs in each cluster were counted. We tested whether average clutch size of individual females was correlated with pgi genotype and flight metabolic rate. Prior to the analysis, average clutch size was corrected for body mass and all statistics were weighted by the number of egg clutches laid by a given female.
Among the 75 genotyped individuals, a total of seven alleles could be distinguished at the pgi locus. Allele frequencies were very similar to those recorded in the same metapopulation in 1995 (table 1; Saccheri et al. 1998; Saccheri et al. 2004), suggesting stable allele frequencies over a period of 9 years. Note, however, that our sample is somewhat biased towards new populations, which we specifically sought for the purpose of the present experiment, whereas populations of intermediate age were not sampled.
In our sample, the frequencies of the alleles pgi-d and pgi-f varied significantly with population age and connectivity (figure 2, table 2). The frequency of pgi-f showed the same pattern among populations as predicted by a model for phenotypic variation in dispersal propensity (figure 1, Hanski et al. 2004). That is, pgi-f was especially common in new, isolated populations, whereas it was rare in old, isolated populations (figure 2). Hence the interaction between population age and connectivity was significant in determining the frequencies of this allele (table 2). For pgi-d the pattern was inverse (figure 2, table 2). For the rarer alleles sample sizes were too small to test whether their frequencies depended on population age and connectivity.
We measured the flight metabolic rate of 45 females. A few (four from new populations, three from old populations) were unable to accomplish weight-supported flight and tended to show signs of advanced age or injury. Thus their data were discarded prior to any analyses. The remaining 38 flight-capable butterflies came from 25 distinct populations.
Flight metabolic rate and the ability to maintain flight continuously for 15min were highly variable (figure 3). All butterflies attained their peak metabolic rate early in the experiment, with varying degrees of decline thereafter. Among females from old populations, 13 of 20 fatigued and were unable to maintain flight for the entire 15min. Among females from new populations, a smaller proportion fatigued (8 of 18), but the difference was not significant (two-tailed Fisher's exact test, p=0.34).
Both the peak metabolic rate (ml CO2 h−1) and the total volume of CO2 emitted during the 15min flight test showed a positive relationship with body mass (two-tailed p=0.028 and p=0.090, respectively). To account for this body size effect, we used residuals from these regressions for all further statistical analyses involving metabolic rates.
We used a general linear model to test whether flight metabolic rate varied among populations of different age and connectivity in a manner parallel to predicted variation in dispersal rate. This model (r2=0.20) showed that the interaction between population age and connectivity had a significant effect on peak metabolic rate (p=0.040, figure 4). In particular, females from the most isolated (low connectivity) new populations had the highest peak metabolic rates.
The total volume of CO2 emitted during the 15min flight test was strongly correlated with peak metabolic rate (r2=0.54), but showed no significant association with population age, connectivity, or the interaction between population age and connectivity. Total metabolism was more variable among individuals than was peak metabolism, which reduced the ability to detect statistically significant differences.
Twenty-five of the females tested for flight metabolic rate were subsequently genotyped at pgi. Because the distribution of allele frequencies among populations (figure 2) suggested that genotypes with the pgi-f allele were better dispersers, we first tested whether there was a difference in metabolic rate between individuals with the pgi-f allele (either homozygous or heterozygous) as compared to genotypes without this allele. Peak metabolic rate was higher in genotypes with the pgi-f allele than in genotypes without this allele (two-tailed t-test, t=2.37, d.f.=23, p=0.026, figure 5). At average body mass, females with pgi-f had a peak metabolic rate of 3.44ml CO2 h−1 (s.e.=0.19), compared to 2.85ml CO2 h−1 (s.e.=0.15) in females without pgi-f (a 17% increase). Overall, 20% of the variance in peak metabolic rate was explained by the presence or absence of pgi-f. Genotypes with pgi-f also differed from genotypes without pgi-f for total CO2 produced during 15min flight (two-tailed t-test, t=2.40, d.f.=23, p=0.025). For pgi-d as well as for all rarer alleles there were no differences in metabolic rate between individuals with and without the allele.
Looking at individual genotypes, the sample sizes were very small, but it is noteworthy that the ranking of peak metabolic rate for the three most common genotypes was pgi-df>pgi-ff>pgi-dd (figure 5). This ranking appears to explain why metabolic rate depended on the presence of pgi-f but not of pgi-d, as individuals with the latter allele had a high (heterozygous) or low (homozygous) metabolic rate.
To test whether genotypes with pgi-f had a higher metabolic rate than genotypes without pgi-f independently of the age and connectivity of their population, we used population age, connectivity, and their interaction in a general linear model to examine variation in peak metabolic rate. In this model, the difference between the pgi-f and non-pgi-f was non-significant (p=0.09), though means between the two genotypes were substantially different (mean ±s.e.=0.23±0.18 versus −0.15±0.14). The effects of population age, connectivity, and their interaction were all clearly non-significant (p>0.36). This result indicates that the effect of pgi on peak metabolic rate is unlikely to be due to a spurious correlation of some other factor with the age and spatial location of populations.
Individuals with pgi-f genotypes had significantly higher body mass-adjusted average clutch sizes than individuals without this allele (two-tailed t-test, t=2.85, d.f.=20, p=0.010, figure 5). The means for females with and without pgi-f were 223 (s.e.=12.7) and 184 (9.7) eggs per clutch, with 29% of the total variance in clutch size being explained by the presence of pgi-f. For pgi-d and the rarer alleles, average clutch size did not depend on the presence of the allele. Among the three most common genotypes, the ranking was again pgi-df>pgi-ff>pgi-dd (figure 5), with two of the three contrasts being significant (pgi-dd versus pgi-ff: two-tailed t-test, p=0.021; pgi-dd versus pgi-df: p=0.006; pgi-df versus pgi-ff: p=0.14). The difference in residual clutch size between genotypes with and without pgi-f remained significant (p=0.020) when population age, connectivity, and the interaction between population age and connectivity were included in the model.
Average clutch size was positively related to peak metabolic rate (p=0.045) and total CO2 emitted during the 15min flight test (p=0.035), but both effects were non-significant (p>0.4) when the presence or absence of the pgi-f allele was included as an independent variable. Thus our data consistently indicate that butterflies with the pgi-f allele have both elevated flight metabolism and larger clutch sizes.
We found a significant association between pgi allele frequencies and flight metabolic rate in the Glanville fritillary butterfly, and that both of these traits varied among local populations with dissimilar ages and spatial connectivities in a manner consistent with field-measured variation in dispersal rate (Hanski et al. 2002; Hanski et al. 2004). A previous study on the Glanville fritillary, which showed parallel variation among populations in the [ATP]/[ADP] ratio in flight muscles following a fixed period of flight (Hanski et al. 2004), suggested that the observed variation in dispersal rate may at least partly be caused by variation in flight physiology. We have now extended this finding by linking variation in dispersal rate and flight physiology to genetic variation in an enzyme with a central role in glycolysis and potentially a direct effect on flight physiology (Watt & Boggs 1987; Watt & Dean 2000).
In Colias butterflies, different pgi-genotypes have different enzyme-kinetic properties (Watt 1977; Watt 1983), which lead to different glycolytic fluxes and flight performance among the pgi genotypes (Watt et al. 1983; Watt & Dean 2000). The pgi polymorphism documented here for the Glanville fritillary appears very similar to that found in Colias butterflies. First, the number and frequency distribution of electrophoretically distinguishable alleles was similar to those in Colias (Watt 1977; Watt 1983; Watt et al. 2003). Second, as in Colias (Watt et al. 2003), the allele frequencies have been stable in samples spanning almost a decade. Third, the ranking of flight capacity among the three most common genotypes is very similar in these two species, with heterozygotes ranking highest (Watt et al. 1983). Nonetheless, it is unlikely that the pgi alleles in the Glanville fritillary are identical to those in Colias, as the two genera diverged ca 40Myr (N. Wahlberg, personal communication) and structural differences of alleles across species are evident within the genus Colias alone (Watt et al. 1996). Rather, these results suggest that similar forces have shaped the evolution of pgi variation in different butterfly species.
In summary, we have presented correlational evidence suggesting an effect of the pgi locus on dispersal rate in the Glanville fritillary. This correlational evidence is supported by a mechanistic hypothesis, which involves an effect of pgi on flight physiology and an effect of flight physiology on flight capacity and hence on dispersal rate. To prove that pgi has a major effect on variation in dispersal rate one would have to manipulate the pgi genotype of individuals experimentally, which is currently not possible in the Glanville fritillary. But the hypothesis can also be further corroborated by analyzing the mechanistic consequences of allelic variation in pgi in greater detail, as has been done for Colias by analyzing the enzyme kinetic properties of different pgi genotypes. As we have not yet established this mechanistic link in the Glanville fritillary, we cannot exclude the alternative hypothesis that the observed variation is due to an unknown locus in close linkage with pgi. Furthermore, pgi is unlikely to be the only gene involved in determining dispersal rate, as theory and empirical data suggest that variation in most quantitative traits is governed by a few loci of large effects and many loci of smaller effects (e.g. Paterson et al. 1991; Orr 1998). Nonetheless, the present data suffice to strongly indicate that pgi is a candidate for a gene with a major effect on dispersal rate in the Glanville fritillary.
Allelic variation at pgi and flight metabolic rate were correlated with average clutch size. In Colias, the kinetically superior pgi genotypes have a higher fecundity because they can fly and oviposit at lower temperatures and thus over a longer period of time during the day than the kinetically less competent genotypes (Watt 1992). Similarly, M. Saastamoinen (unpublished data) has observed that Glanville fritillary females with the allele pgi-f lay earlier in the day, when temperatures tend to be higher and when all females typically lay larger clutches. Possibly females with the pgi-f allele become physiologically able to lay earlier in the day than the non-pgi-f females, and can thereby take advantage of the high noon temperatures. Regardless of the actual explanation of why metabolically superior and hence more dispersive females lay larger clutches, this association is consistent with earlier results (Hanski et al. in press) showing no direct trade-off between dispersal rate and fecundity in the Glanville fritillary, in contrast to what is commonly found in wing-dimorphic insects (Roff 1977; Zera & Denno 1997).
Previous studies on the Glanville fritillary have shown that phenotypic and physiological variation related to dispersal rate and flight capacity are maintained among local populations of different ages and spatial connectivities in a manner predicted by an evolutionary model of dispersal (Hanski et al. 2002; Hanski et al. 2004). Here we have shown that the same is true for genetic variation at a locus that appears to be causally related to variation in dispersal rate and physiology. Recent studies on bird populations illustrate other ways in which non-random dispersal combined with spatial variation in the expression of genetic variation can lead to unexpected population differentiation (Garant et al. 2005; Postma & Noordwijk 2005). Though the mechanisms are different, these studies provide further evidence that dispersal and gene flow among distinct local populations does not necessarily homogenize allele frequencies (Slatkin 1987). In our study system, individuals with a higher dispersal rate are more likely to reach unoccupied patches, especially if isolated from existing populations; this creates a non-random distribution of genotypes among populations. As populations persist, dispersive genotypes appear to be lost due to emigration and populations tend to become more sedentary. Genetic variation for dispersal in this system can be maintained by the spatially and temporally dynamic processes of extinction and recolonization (Hanski et al. 2004). Furthermore, as the ecological metapopulation dynamics depend critically on dispersal rate, variation in the dispersal rate is likely to have a feed-back effect on metapopulation dynamics.
This study and the previous work on the Glanville fritillary (Hanski et al. 2002; Hanski et al. 2004) provide the beginnings of a mechanistic explanation of non-random dispersal of butterfly phenotypes in a metapopulation. Frequent extinction of small local populations creates opportunities for migrant females to establish new populations. Our results indicate that females with the pgi-f allele are stronger fliers and hence more likely to establish new populations, especially when the new populations are poorly connected to existing populations and hence the colonizers have to fly for a long distance. The exciting prospect raised by these results is the possibility of developing a truly mechanistic understanding of the spatial and temporal dynamics in genes with major phenotypic effects in metapopulations in real fragmented landscapes. Metapopulation dynamics are a pervasive and general feature of species living in fragmented landscapes (Hanski 1999; Ehrlich & Hanski 2004). Therefore, although this paper examines particular traits of a single species, our approach is likely to be applicable across a wide range of species, including perhaps humans, who show genetic variation in metabolic rates and activity (Kubaszek et al. 2003; Payne & Montgomery 2003; Yang et al. 2003).
We thank the staff of Tvärminne Zoological Station and the Åland 2004 team for support and technical assistance. K. Niitepõld helped with the measurements of flight metabolic rate and shared unpublished data. We profited from comments by C. Liautard, C. W. Wheat, and two anonymous reviewers. T. Roslin and T. J. Kawecki provided statistical advice, and E. Meyke helped with figure preparation. The work was supported by the Academy of Finland (grants no. 38604 and 44887 to I.H., Finnish Centre of Excellence Programme, 2000–2005) and grants IBN-0091040 and IBN-0412651 to J.H.M. from the US National Science Foundation.
As this paper exceeds the maximum length normally permitted, the authors have agreed to contribute to production costs.