Stochastic models for mainland-island metapopulations in static and dynamic landscapes

This paper has three primary aims: to establish an effective means for modelling mainland-island metapopulations inhabiting a dynamic landscape: to investigate the effect of immigration and dynamic changes in habitat on metapopulation patch occupancy dynamics; and to illustrate the implications of our results for decision-making and population management. We first extend the mainland-island metapopulation model of Alonso and McKane [Bull. Math. Biol. 64:913-958,2002] to incorporate a dynamic landscape. It is shown, for both the static and the dynamic landscape models, that a suitably scaled version of the process converges to a unique deterministic model as the size of the system becomes large. We also establish that. under quite general conditions, the density of occupied patches, and the densities of suitable and occupied patches, for the respective models, have approximate normal distributions. Our results not only provide us with estimates for the means and variances that are valid at all stages in the evolution of the population, but also provide a tool for fitting the models to real metapopulations. We discuss the effect of immigration and habitat dynamics on metapopulations, showing that mainland-like patches heavily influence metapopulation persistence, and we argue for adopting measures to increase connectivity between this large patch and the other island-like patches. We illustrate our results with specific reference to examples of populations of butterfly and the grasshopper Bryodema tuberculata.

Host-parasitoid relationship in different Cotesia melitaearum and Melitaea cinxia populations around the Baltic Sea

The relationship between hosts and parasites is one of the most studied interactions between living organisms, and it is both universal and common in nature. Parasitoids are special type of parasites whose offspring develop attached to or within a single host organism that it ultimately consumes and kills. Hosts are arthropods and most parasitoids belong to the insect order Hymenoptera. For almost two decades metapopulation research on the Glanville fritillary butterfly (Melitaea cinxia) has been conducted in the Åland Islands, Finland. The studies have been concerned with the population dynamics, evolution, genetics, behavior, natural history and life history characteristics of M. cinxia, as well as with species interacting with the butterfly. The parasitoids of M. cinxia have been under long term studies and much has been learned about specific host-parasitoid interactions during the past decade. The research for this Master s thesis was done in the Åland Islands during summer 2010. I conducted a reciprocal transplant style experiment in order to compare the performance of host butterflies (M. cinxia) under attack by different parasitoid wasps (C. melitaearum). I used hosts and parasitoids from five origins around the Baltic Sea: Öland, Uppland, Åland, Saaremaa and Pikku-Tytärsaari. The host-parasitoid relationship was studied in terms of host susceptibility and parasitoid virulence, addressing specifically the possible effects of inbreeding and local adaptation of both parasitoids and their hosts. I compared various factors such as host defence ratio, parasitoid development rate, cocoon production rate etc. I also conducted a small scale C. melitaearum egg development experiment and C. melitaearum external morphology comparison between different parasitoid populations. The results show that host resistance and parasitoid virulence differ between both host and parasitoid populations. For example, Öland hosts were most susceptible to parasitoids and especially vulnerable to Pikku-Tytärsaari wasps. Pikku-Tytärsaari wasps were most successful in terms of parasitoids virulence and efficiency except in Saaremaa hosts, where the wasp did not succeed. Saaremaa hosts were resistant except towards Åland parasitoids. I did not find any simple pattern concerning host resistance and parasitoid virulence between inbred and outbred populations. Also, the effect of local adaptation was not detected, perhaps because metapopulation processes disturb local adaptation of the studied populations. Morphological comparisons showed differences between studied wasp populations and sexual dimorphism was obvious with females being bigger that males. There were also interesting differences among populations in male and female wing shapes. The results raise many further questions. Especially interesting were Pikku-Tytärsaari wasps that did well in terms of efficiency and virulence. Pikku-Tytärsaari is a small, isolated island in the Gulf of Finland and both the host and parasitoids are extremely inbred. For the host and parasitoid to persist in the island, the host has to have some mechanisms to escape the parasitoid. Further research will be done on the subject to discover the mechanisms of Pikku-Tytärsaari host s ability to escape parasitism. Also, genetic analyses will be conducted in the near future to determine the relatedness of used C. melitaearum populations.

Geospatial environmental data modelling applications using remote sensing, GIS and spatial statistics

This thesis presents novel modelling applications for environmental geospatial data using remote sensing, GIS and statistical modelling techniques. The studied themes can be classified into four main themes: (i) to develop advanced geospatial databases. Paper (I) demonstrates the creation of a geospatial database for the Glanville fritillary butterfly (Melitaea cinxia) in the Åland Islands, south-western Finland; (ii) to analyse species diversity and distribution using GIS techniques. Paper (II) presents a diversity and geographical distribution analysis for Scopulini moths at a world-wide scale; (iii) to study spatiotemporal forest cover change. Paper (III) presents a study of exotic and indigenous tree cover change detection in Taita Hills Kenya using airborne imagery and GIS analysis techniques; (iv) to explore predictive modelling techniques using geospatial data. In Paper (IV) human population occurrence and abundance in the Taita Hills highlands was predicted using the generalized additive modelling (GAM) technique. Paper (V) presents techniques to enhance fire prediction and burned area estimation at a regional scale in East Caprivi Namibia. Paper (VI) compares eight state-of-the-art predictive modelling methods to improve fire prediction, burned area estimation and fire risk mapping in East Caprivi Namibia. The results in Paper (I) showed that geospatial data can be managed effectively using advanced relational database management systems. Metapopulation data for Melitaea cinxia butterfly was successfully combined with GPS-delimited habitat patch information and climatic data. Using the geospatial database, spatial analyses were successfully conducted at habitat patch level or at more coarse analysis scales. Moreover, this study showed it appears evident that at a large-scale spatially correlated weather conditions are one of the primary causes of spatially correlated changes in Melitaea cinxia population sizes. In Paper (II) spatiotemporal characteristics of Socupulini moths description, diversity and distribution were analysed at a world-wide scale and for the first time GIS techniques were used for Scopulini moth geographical distribution analysis. This study revealed that Scopulini moths have a cosmopolitan distribution. The majority of the species have been described from the low latitudes, sub-Saharan Africa being the hot spot of species diversity. However, the taxonomical effort has been uneven among biogeographical regions. Paper III showed that forest cover change can be analysed in great detail using modern airborne imagery techniques and historical aerial photographs. However, when spatiotemporal forest cover change is studied care has to be taken in co-registration and image interpretation when historical black and white aerial photography is used. In Paper (IV) human population distribution and abundance could be modelled with fairly good results using geospatial predictors and non-Gaussian predictive modelling techniques. Moreover, land cover layer is not necessary needed as a predictor because first and second-order image texture measurements derived from satellite imagery had more power to explain the variation in dwelling unit occurrence and abundance. Paper V showed that generalized linear model (GLM) is a suitable technique for fire occurrence prediction and for burned area estimation. GLM based burned area estimations were found to be more superior than the existing MODIS burned area product (MCD45A1). However, spatial autocorrelation of fires has to be taken into account when using the GLM technique for fire occurrence prediction. Paper VI showed that novel statistical predictive modelling techniques can be used to improve fire prediction, burned area estimation and fire risk mapping at a regional scale. However, some noticeable variation between different predictive modelling techniques for fire occurrence prediction and burned area estimation existed.

Behavioural, population, and genetic processes affecting metapopulation dynamics of the Glanville fritillary butterfly

In my thesis I have been studying the effects of population fragmentation and extinction-recolonization dynamics on genetic and evolutionary processes in the Glanville fritillary butterfly (Melitaea cinxia). By conducting crosses within and among newly-colonized populations and using several fitness measures, I found a strong decrease in fitness following colonization by a few related individuals, and a strong negative relationship between parental relatedness and offspring fitness. Thereafter, I was interested in determining the number and relatedness of individuals colonizing new populations, which I did using a set of microsatellites I had previously developed for this species. Additionally, I am interested in the evolution of key life-history traits. By following the lifetime reproductive success of males emerging at different times in a semi-natural setup, I demonstrated that protandry is adaptive in males, and I was able to rule out, for M. cinxia, alternative incidental hypotheses evoked to explain the evolution of protandry in insects. Finally, in work I did together with Prof. Hanna Kokko, I am proposing bet-hedging as a new mechanism that could explain the evolution of polyandry in M. cinxia.

Dispersal and related life history traits in the Glanville fritillary butterfly

Most studies of life history evolution are based on the assumption that species exist at equilibrium and spatially distinct separated populations. In reality, this is rarely the case, as populations are often spatially structured with ephemeral local populations. Therefore, the characteristics of metapopulations should be considered while studying factors affecting life history evolution. Theoretical studies have examined spatial processes shaping the evolution of life history traits to some extent, but there is little empirical data and evidence to investigate model predictions. In my thesis I have tried to bridge the gap between theoretical and empirical studies by using the well-known Glanville fritillary (Melitaea cinxia) metapopulation as a model system. The long-term persistence of classic metapopulations requires sufficient dispersal to establish new local populations to compensate for local extinctions. Previous studies on the Glanville fritillary have shown that females establishing new populations are not a random sample from the metapopulation, but they are in fact more dispersive than females in old populations. Many other life-history traits, such as body size, fecundity and lifespan, may be related to dispersal rate. Therefore, I examined a range of correlated traits for their evolutionary and ecological consequences. I was particularly interested in how the traits vary under natural environmental conditions, hence all studies were conducted in a large (32 x 26 m) outdoor population cage built upon a natural habitat patch. Individuals for the experiments were sampled from newly-established and old populations within a large metapopulation. Results show that females originating from newly-established populations had higher within-habitat patch mobility than females from old populations. I showed that dispersal rate is heritable and that flight activity is related to variation in a gene encoding the glycolytic enzyme phosphoglucose isomerase. Both among-individual and among-population variation in dispersal are correlated with the reproductive performance of females, though I found no evidence for a trade-off between dispersal and fecundity in terms of lifetime egg production or clutch size. Instead, the results suggest that highly dispersive females from newly-established populations have a shorter lifespan than females from old populations, and that dispersive females may pay a cost in terms of reduced lifetime reproductive success due to increased time spent outside habitat patches. In summary, the results of this thesis show that genotype-dependent dispersal rate correlates with other life history traits in the Glanville fritillary, and that the rapid turnover of local populations (extinctions and re-colonisations) is likely to be the mechanism that maintains phenotypic variation in many life history traits at the metapopulation level.

Flight metabolic rate in the Glanville fritillary butterfly

Dispersal is a highly important life history trait. In fragmented landscapes the long-term persistence of populations depends on dispersal. Evolution of dispersal is affected by costs and benefits and these may differ between different landscapes. This results in differences in the strength and direction of natural selection on dispersal in fragmented landscapes. Dispersal has been shown to be a nonrandom process that is associated with traits such as flight ability in insects. This thesis examines genetic and physiological traits affecting dispersal in the Glanville fritillary butterfly (Melitaea cinxia). Flight metabolic rate is a repeatable trait representing flight ability. Unlike in many vertebrates, resting metabolic rate cannot be used as a surrogate of maximum metabolic rate as no strong correlation between the two was found in the Glanville fritillary. Resting and flight metabolic rate are affected by environmental variables, most notably temperature. However, only flight metabolic rate has a strong genetic component. Molecular variation in the much-studied candidate locus phosphoglucose isomerase (Pgi), which encodes the glycolytic enzyme PGI, has an effect on carbohydrate metabolism in flight. This effect is temperature dependent: in low to moderate temperatures individuals with the heterozygous genotype at the single nucleotide polymorphism (SNP) AA111 have higher flight metabolic rate than the common homozygous genotype. At high temperatures the situation is reversed. This finding suggests that variation in enzyme properties is indeed translated to organismal performance. High-resolution data on individual female Glanville fritillaries moving freely in the field were recorded using harmonic radar. There was a strong positive correlation between flight metabolic rate and dispersal rate. Flight metabolic rate explained one third of the observed variation in the one-hour movement distance. A fine-scaled analysis of mobility showed that mobility peaked at intermediate ambient temperatures but the two common Pgi genotypes differed in their reaction norms to temperature. As with flight metabolic rate, heterozygotes at SNP AA111 were the most active genotype in low to moderate temperatures. The results show that molecular variation is associated with variation in dispersal rate through the link of flight physiology under the influence of environmental conditions. The evolutionary pressures for dispersal differ between males and females. The effect of flight metabolic rate on dispersal was examined in both sexes in field and laboratory conditions. The relationship between flight metabolic rate and dispersal rate in the field and flight duration in the laboratory were found to differ between the two sexes. In females the relationship was positive, but in males the longest distances and flight durations were recorded for individuals with low flight metabolic rate. These findings may reflect male investment in mate locating. Instead of dispersing, males with high flight metabolic rate may establish territories and follow a perching strategy when locating females and hence move less on the landscape level. Males with low metabolic rate may be forced to disperse due to low competitive success or may show adaptations to an alternative strategy: patrolling. In the light of life history trade-offs and the rate of living theory having high metabolic rate may carry a cost in the form of shortened lifespan. Experiments relating flight metabolic rate to longevity showed a clear correlation in the opposite direction: high flight metabolic rate was associated with long lifespan. This suggests that individuals with high metabolic rate do not pay an extra physiological cost for their high flight capacity, rather there are positive correlations between different measures of fitness. These results highlight the importance of condition.

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 nonreproducing 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 genotype!environment interactions, and these trade-offs are likely to contribute to the maintenance of Pgi polymorphism in the natural populations.

Association between molecular variation in the phosphoglucose isomerase (Pgi) locus and migration propensity in the Glanville fritillary butterfly

Habitat fragmentation produces patches of suitable habitat surrounded by unfavourable matrix habitat. A species may persist in such a fragmented landscape in an equilibrium between the extinctions and recolonizations of local populations, thus forming a metapopulation. Migration between local populations is necessary for the long-term persistence of a metapopulation. The Glanville fritillary butterfly (Melitaea cinxia) forms a metapopulation in the Åland islands in Finland. There is migration between the populations, the extent of which is affected by several environmental factors and variation in the phenotype of individual butterflies. Different allelic forms of the glycolytic enzyme phosphoglucose isomerase (Pgi) has been identified as a possible genetic factor influencing flight performance and migration rate in this species. The frequency of a certain Pgi allele, Pgi-f, follows the same pattern in relation to population age and connectivity as migration propensity. Furthermore, variation in flight metabolic performance, which is likely to affect migration propensity, has been linked to genetic variation in Pgi or a closely linked locus. The aim of this study was to investigate the association between Pgi genotype and the migration propensity in the Glanville fritillary both at the individual and population levels using a statistical modelling approach. A mark-release-recapture (MRR) study was conducted in a habitat patch network of M. cinxia in Åland to collect data on the movements of individual butterflies. Larval samples from the study area were also collected for population level examinations. Each butterfly and larva was genotyped at the Pgi locus. The MRR data was parameterised with two mathematical models of migration: the Virtual Migration Model (VM) and the spatially explicit diffusion model. VM model predicted and observed numbers of emigrants from populations with high and low frequencies of Pgi-f were compared. Posterior predictive data sets were simulated based on the parameters of the diffusion model. Lack-of-fit of observed values to the model predicted values of several descriptors of movements were detected, and the effect of Pgi genotype on the deviations was assessed by randomizations including the genotype information. This study revealed a possible difference in the effect of Pgi genotype on migration propensity between the two sexes in the Glanville fritillary. The females with and males without the Pgi-f allele moved more between habitat patches, which is probably related to differences in the function of flight in the two sexes. Females may use their high flight capacity to migrate between habitat patches to find suitable oviposition sites, whereas males may use it to acquire mates by keeping a territory and fighting off other intruding males, possibly causing them to emigrate. The results were consistent across different movement descriptors and at the individual and population levels. The effect of Pgi is likely to be dependent on the structure of the landscape and the prevailing environmental conditions.