Modeling ecological and evolutionary processes in spatially structured populations

Many species inhabit fragmented landscapes, resulting either from anthropogenic or from natural processes. The ecological and evolutionary dynamics of spatially structured populations are affected by a complex interplay between endogenous and exogenous factors. The metapopulation approach, simplifying the landscape to a discrete set of patches of breeding habitat surrounded by unsuitable matrix, has become a widely applied paradigm for the study of species inhabiting highly fragmented landscapes. In this thesis, I focus on the construction of biologically realistic models and their parameterization with empirical data, with the general objective of understanding how the interactions between individuals and their spatially structured environment affect ecological and evolutionary processes in fragmented landscapes. I study two hierarchically structured model systems, which are the Glanville fritillary butterfly in the Åland Islands, and a system of two interacting aphid species in the Tvärminne archipelago, both being located in South-Western Finland. The interesting and challenging feature of both study systems is that the population dynamics occur over multiple spatial scales that are linked by various processes. My main emphasis is in the development of mathematical and statistical methodologies. For the Glanville fritillary case study, I first build a Bayesian framework for the estimation of death rates and capture probabilities from mark-recapture data, with the novelty of accounting for variation among individuals in capture probabilities and survival. I then characterize the dispersal phase of the butterflies by deriving a mathematical approximation of a diffusion-based movement model applied to a network of patches. I use the movement model as a building block to construct an individual-based evolutionary model for the Glanville fritillary butterfly metapopulation. I parameterize the evolutionary model using a pattern-oriented approach, and use it to study how the landscape structure affects the evolution of dispersal. For the aphid case study, I develop a Bayesian model of hierarchical multi-scale metapopulation dynamics, where the observed extinction and colonization rates are decomposed into intrinsic rates operating specifically at each spatial scale. In summary, I show how analytical approaches, hierarchical Bayesian methods and individual-based simulations can be used individually or in combination to tackle complex problems from many different viewpoints. In particular, hierarchical Bayesian methods provide a useful tool for decomposing ecological complexity into more tractable components.

The range expansion of the European map butterfly in Finland

Climate change will influence the living conditions of all life on Earth. For some species the change in the environmental conditions that has occurred so far has already increased the risk of extinction, and the extinction risk is predicted to increase for large numbers of species in the future. Some species may have time to adapt to the changing environmental conditions, but the rate and magnitude of the change are too great to allow many species to survive via evolutionary changes. Species responses to climate change have been documented for some decades. Some groups of species, like many insects, respond readily to changes in temperature conditions and have shifted their distributions northwards to new climatically suitable regions. Such range shifts have been well documented especially in temperate zones. In this context, butterflies have been studied more than any other group of species, partly for the reason that their past geographical ranges are well documented, which facilitates species-climate modelling and other analyses. The aim of the modelling studies is to examine to what extent shifts in species distributions can be explained by climatic and other factors. Models can also be used to predict the future distributions of species. In this thesis, I have studied the response to climate change of one species of butterfly within one geographically restricted area. The study species, the European map butterfly (Araschnia levana), has expanded rapidly northwards in Finland during the last two decades. I used statistical and dynamic modelling approaches in combination with field studies to analyse the effects of climate warming and landscape structure on the expansion. I studied possible role of molecular variation in phosphoglucose isomerase (PGI), a glycolytic enzyme affecting flight metabolism and thereby flight performance, in the observed expansion of the map butterfly at two separate expansion fronts in Finland. The expansion rate of the map butterfly was shown to be correlated with the frequency of warmer than average summers during the study period. The result is in line with the greater probability of occurrence of the second generation during warm summers and previous results on this species showing greater mobility of the second than first generation individuals. The results of a field study in this thesis indicated low mobility of the first generation butterflies. Climatic variables alone were not sufficient to explain the observed expansion in Finland. There are also problems in transferring the climate model to new regions from the ones from which data were available to construct the model. The climate model predicted a wider distribution in the south-western part of Finland than what has been observed. Dynamic modelling of the expansion in response to landscape structure suggested that habitat and landscape structure influence the rate of expansion. In southern Finland the landscape structure may have slowed down the expansion rate. The results on PGI suggested that allelic variation in this enzyme may influence flight performance and thereby the rate of expansion. Genetic differences of the populations at the two expansion fronts may explain at least partly the observed differences in the rate of expansion. Individuals with the genotype associated with high flight metabolic rate were most frequent in eastern Finland, where the rate of range expansion has been highest.

Effects of management and landscape structure on biodiversity in boreal agricultural farmland

Intensified agricultural practises introduced after the Second World War are identified as a major cause of global biodiversity declines. In several European countries agri-environment support schemes have been introduced to counteract the ongoing biodiversity declines. Farmers participating in agri-environment schemes are financially compensated for decreasing the intensity of farming practises leading to smaller yields and lower income. The Finnish agri-environment support scheme is composed of a set of measures, such as widened field margins along main ditches (obligatory measure), management of features increasing landscape diversity, management of semi-natural grasslands, and organic farming (special agreement measures). The magnitude of the benefits for biodiversity depends on landscape context and the properties of individual schemes. In this thesis I studied whether one agri-environment scheme, organic farming, is beneficial for species diversity and abundance of diurnal lepidopterans, bumblebees, carabid beetles and arable weeds. I found that organic farming did not enhance species richness of selected insect taxa, although bumblebee species richness tended to be higher in organic farms. Abundance of lepidopterans and bumblebees was not enhanced by organic farming, but carabid beetle abundance was higher in mixed farms with both cereal crop production and animal husbandry. Both species richness and abundance of arable weeds were higher in organic farms. My second objective was to study how landscape structure shapes farmland butterfly communities. I found that the percentage of habitat specialists and species with poor dispersal abilities in butterfly assemblages decreased with increasing arable field cover, leading to a dramatic decrease in butterfly beta diversity. In field boundaries local species richness of butterflies was linearly related to landscape species richness in geographic regions with high arable field cover, indicating that butterfly species richness in field boundaries is more limited by landscape factors than local habitat factors. In study landscapes containing semi-natural grasslands the relationship decelerated at high landscape species richness, suggesting that local species richness of butterflies in field boundaries is limited by habitat factors (demanding habitat specialists that occurred in semi-natural grasslands were absent in field margins). My results suggest that management options in field margins will affect mainly generalists, and species with good dispersal abilities, in landscapes with high arable field cover. Habitat specialists and species with poor dispersal abilities may benefit of management options if these are applied in the vicinity of source populations.

Phylogenetics of Myrmica ants and their social parasites

Understanding the overwhelming diversity of life calls for complex organisational schemes. The field of systematics may thus be seen as the cornerstone of evolutionary biology. In the last few decades, systematics has been rejuvenated through the introduction of molecular methods such as DNA barcoding and multi-gene phylogenetic approaches. These methods may shed new light on established taxonomic ideas and problems. For example, the classification of ants has aroused much debate due to reinterpretation of morphological characters or contradictions between molecular data and morphology. Only in the last few years a consensus was reached regarding the phylogeny of ant subfamilies. However, the situation remains deplorable for lower taxonomic ranks such as subfamilies, tribes and genera. This thesis describes the systematics and evolution of the Holarctic ant genus Myrmica and the tribe to which it belongs, Myrmicini. Using barcoding, molecular-phylogenetic data and divergence time estimations, it addresses questions regarding the taxonomy, morphology and biogeography of this group. Furthermore, the interrelationships between socially parasitic Myrmica species and their hosts (other species in the genus) were inferred. The phylogeny suggests that social parasitism evolved several times in Myrmica. Finally, this thesis investigated whether coevolution shaped the phylogeny of socially parasitic Maculinea butterflies that live inside Myrmica colonies. No evidence was found for coevolution.

Management of semi-natural grasslands for butterfly and moth communities

This work focuses on the factors affecting species richness, abundance and species composition of butterflies and moths in Finnish semi-natural grasslands, with a special interest in the effects of grazing management. In addition, an aim was set at evaluating the effectiveness of the support for livestock grazing in semi-natural grasslands, which is included in the Finnish agri-environment scheme. In the first field study, butterfly and moth communities in resumed semi-natural pastures were com-pared to old, annually grazed and abandoned previous pastures. Butterfly and moth species compo-sition in restored pastures resembled the compositions observed in old pastures after circa five years of resumed cattle grazing, but diversity of butterflies and moths in resumed pastures remained at a lower level compared with old pastures. None of the butterfly and moth species typical of old pas-tures had become more abundant in restored pastures compared with abandoned pastures. There-fore, it appears that restoration of butterfly and moth communities inhabiting semi-natural grass-lands requires a longer time that was available for monitoring in this study. In the second study, it was shown that local habitat quality has the largest impact on the occurrence and abundance of butterflies and moths compared to the effects of grassland patch area and connec-tivity of the regional grassland network. This emphasizes the importance of current and historical management of semi-natural grasslands on butterfly and moth communities. A positive effect of habitat connectivity was observed on total abundance of the declining butterflies and moths, sug-gesting that these species have strongest populations in well-connected habitat networks. Highest species richness and peak abundance of most individual species of butterflies and moths were generally observed in taller grassland vegetation compared with vascular plants, suggesting a preference towards less intensive management in insects. These differences between plants and their insect herbivores may be understood in the light of both (1) the higher structural diversity of tall vegetation and (2) weaker tolerance of disturbances by herbivorous insects due to their higher trophic level compared to plants. The ecological requirements of all species and species groups inhabiting semi-natural grasslands are probably never met at single restricted sites. Therefore, regional implementation of management to create differently managed areas is imperative for the conservation of different species and species groups dependent on semi-natural grasslands. With limited resources it might be reasonable to focus much of the management efforts in the densest networks of suitable habitat to minimise the risk of extinction of the declining species.

Phylogeny of a taxonomically difficult group and evolution of host location mechanism

The parasitic wasps are one of the largest insect groups and their life histories are remarkably variable. Common to all parasitic wasps is that they kill their hosts, which are usually beetles, butterflies and sometimes spiders. Hosts are often at a larval or pupal stage and live in concealed conditions, such as in plant tissue. Parasitic wasps have two main ways of finding their host. 1) They can detect chemical compounds emitted by damaged plant material or released by larvae living in plant tissue, and 2) detect the larvae by sound vibrations. Even though pupae are immobile and silent, and therefore do not cause vibration, parasitoids have, however, adapted to find passive developmental stages by producing vibration themselves by knocking the substrate with their antennae, and then detecting the echoes with their legs. This echolocation allows a parasitoid to locate its potential hosts that are deeply buried in wood. This study focuses on the relationships of the subfamily Cryptinae (Hymenoptera: Ichneumonidae) and related taxa, and the evolution of host location mechanism. There are no earlier studies of the phylogeny of the Cryptinae, and the position of related taxa are unclear. According to the earlier classification, which is entirely intuitional, the Cryptinae is divided into three tribes: Cryptini, Hemigasterini and Phygadeuontini. Further, these tribes are subdiveded into numerous subtribes. This work, based on molecular characters, shows that the cryptine tribes Cryptini, Phygadeuon¬tini and Hemigasterini come out largely as monophyletic groups, thus agreeing with the earlier classification. The earlier subtribal classification had no support. In addition, it is shown that modified antennal structures are associated with host usage of wood-boring coleopteran hosts. The cryptines have a clear modification series on their antennal tips from a simply tip to a hammer-like structure. The species with strongly modified antennae belong mostly to the tribe Cryptini and they utilise wood-boring beetles as hosts. Also, field observations on insect behaviour support this result.

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.

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.

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.