The Population Consequences of Sex and Conflict

Natural selection generally operates at the level of the individual, or more specifically at the level of the gene. As a result, individual selection does not always favour traits which benefit the population or species as a whole. The spread of an individual gene may even act to the detriment of the organism in which it finds. Thus selection at the level of the individual can affect processes at the level of the organism, group or even at the level of the species. As most behaviours ultimately affect births, deaths and the distribution of individuals, it seems inevitable that behavioural decisions will have an impact on population dynamics and population densities. Behavioural decisions can often involve costs through allocation of energy into behavioural strategies, such as the investment into armaments involved in fighting over resources or increased mortality due to injury or increased predation risk. Similarly, behaviour may act o to benefit the population, in terms of higher survival and increased fecundity. Examples include increased investment through parental care, choosing a mate based on the nuptial gifts they may supply and choosing territories in the face of competition. Investigating the impact of behaviour on population ecology may seem like a trivial task, but it is likely to have important consequences at different levels. For example, antagonistic behaviour may occasionally become so extreme that it increases the risk of extinction, and such extinction risk may have important implications for conservation. As a corollary, any such behaviour may also act as a macroevolutionary force, weeding out populations with traits which, whilst beneficial to the individuals in the short term, ultimately result in population extinction. In this thesis, I examine how behaviours, specifically conflict and competition over a resource and aspects of behaviour involved in sexual selection, can affect population densities, and what the implications are for the evolution and ecology of the populations in question. It is found that both behaviours related to individual conflict and mating strategies can have an effect at the level of the population, but that various factors, such as a feedback between selection and population densities or macroevolution caused by species extinctions, may act to limit the intensity of conflicts that we observe in nature.

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.

Climate forcing on avian life history

In the 21st century, human-induced global climate change has been highlighted as one of the most serious threats to ecosystems worldwide. According to global climate scenarios, the mean temperature in Finland is expected to increase by 1.8 4.0°C by the end of the century. The regional and seasonal change in temperature has predicted to be spatially and temporally asymmetric, where the High-Arctic and Antarctic areas and winter and spring seasons have been projected to face the highest temperature increase. To understand how species respond to the ongoing climate change, we need to study how climate affects species in different phases of their life cycle. The impact of climate on breeding and migration of eight large-sized bird species was studied in this thesis, taking food availability into account. The findings show that climatic variables have considerable impact on the life-history traits of large-sized birds in northern Europe. The magnitude of climatic effects on migration and breeding was comparable with that of food supply, conventionally regarded as the main factor affecting these life-history traits. Based on the results of this thesis and the current climate scenarios, the following not mutually exclusive responses are possible in the near future. Firstly, asymmetric climate change may result in a mistiming of breeding because mild winters and early spring may lead to earlier breeding, whereas offspring are hatching into colder conditions which elevate mortality. Secondly, climate induced responses can differ between species with different breeding tactics (income vs. capital breeding), so that especially capital breeders can gain advantage on global warming as they can sustain higher energy resources. Thirdly, increasing precipitation has the potential to reduce the breeding success of many species by exposing nestlings to more severe post-hatching conditions and hampering the hunting conditions of parents. Fourthly, decreasing ice cover and earlier ice-break in the Baltic Sea will allow earlier spring migration in waterfowl. In eiders, this can potentially lead to more productive breeding. Fifthly, warming temperatures can favour parents preparing for breeding and increase nestling survival. Lastly, the climate-induced phenological changes in life history events will likely continue. Furthermore, interactions between climate and food resources can be complex and interact with each other. Eiders provide an illustrative example of this complexity, being caught in the crossfire between more benign ice conditions and lower salinity negatively affecting their prime food resource. The general conclusion is that climate is controlling not only the phenology of the species but also their reproductive output, thus affecting the entire population dynamics.

Phylogeography and Population Genetics of the Moor frog (Rana arvalis Nilsson) in Northern Europe

The ongoing climate change along with increasing levels of pollutants, diseases, habitat loss and fragmentation constitute global threats to the persistence of many populations, species and ecosystems. However, for the long-term persistence of local populations, one of the biggest threats is the intrinsic loss of genetic variation. In order to adapt to changes in the environment, organisms must have a sufficient supply of heritable variation in traits important for their fitness. With a loss of genetic variation, the risk of extinction will increase. For conservational practices, one should therefore understand the processes that shape the genetic population structure and also the broader (historical) phylogenetic patterning of the species in focus. In this thesis, microsatellite markers were applied to study genetic diversity and population differentiation of the protected moor frog (Rana arvalis) in Fennoscandia from both historical (evolutionary) and applied (conservation) perspectives. The results demonstrate that R. arvalis populations are highly structured over rather short geographic distances. Moreover, the results suggest that R. arvalis recolonized Fennoscandia from two directions after the last ice age. This has had implications for the genetic structuring and population differentiation, especially in the northernmost parts where the two lineages have met. Compared to more southern populations, the genetic variation decreases and the interpopulation differentiation increases dramatically towards north. This could be an outcome of serial population bottlenecking along the recolonization route. Also, current isolation and small population sizes increase the effect of drift, thus reinforcing the observed pattern. The same pattern can also be seen in island populations. However, though R. arvalis on the island of Gotland has lost most of its neutral genetic variability, our results indicate that the levels of additive genetic variation have remained high. This conforms to the conjecture that though neutral markers are widely used in conservation purposes, they may be quite uninformative about the levels of genetic variation in ecologically important traits. Finally, the evolutionary impact of the typical amphibian mating behaviour on genetic diversity was investigated. Given the short time available for larval development, it is important that mating takes place as early as possible. The genetic data and earlier capture-recapture data suggest that R. arvalis gather at mating grounds they are familiar with. However, by forming leks in random to relatedness, and having multiple paternities in single clutches, the risk of inbreeding may be minimized in this otherwise highly philopatric species.

Living on voles - plastic life of the Ural owl

Individuals face variable environmental conditions during their life. This may be due to migration, dispersion, environmental changes or, for example, annual variation in weather conditions. Genetic adaptation to a novel environment happens through natural selection. Phenotypic plasticity allows, however, a quick individual response to a new environment. Phenotypic plasticity may also be beneficial for individual if the environment is highly variable. For example, eggs are costly to produce. If the food conditions vary significantly between breeding seasons it is useful to be able to adjust the clutch and egg size according to the food abundance. In this thesis I use Ural owl vole system to study phenotypic plasticity and natural selection using a number of reproduction related traits. The Ural owl (Strix uralensis) is a long-lived and sedentary species. The reproduction and survival of the Ural owl, in fact their whole life, is tied to the dramatically fluctuating vole densities. Ural owls do not cause vole cycles but they have to adjust their behaviour to the rather predictable population fluctuations of these small mammals. Earlier work with this system has shown that Ural owl laying date and clutch size are plastic in relation to vole abundance. Further, individual laying date clutch size reaction norms have been shown to vary in the amount of plasticity. My work extends the knowledge of natural selection and phenotypic plasticity in traits related to reproduction. I show that egg size, timing of the onset of incubation and nest defense aggressiveness are plastic traits with fitness consequences for the Ural owl. Although egg size is in general thought to be a fixed characteristic of an individual, this highly heritable trait in the Ural owl is also remarkably plastic in relation to the changes in vole numbers, Ural owls are laying the largest eggs when their prey is most abundant. Timing of the onset of incubation is an individual-specific property and plastic in relation to clutch size. Timing of incubation is an important underlying cause for asynchronous hatching in birds. Asynchronous hatching is beneficial to offspring survival in Ural owl. Hence, timing of the onset of incubation may also be under natural selection. Ural owl females also adjust their nest defense aggressiveness according to the vole dynamics, being most aggressive in years when they produce the largest broods. Individual females show different levels of nest defense aggressiveness. Aggressiveness is positively correlated with the phenotypic plasticity of aggressiveness. As elevated nest defense aggressiveness is selected for, it may promote the plasticity of aggressive nest defense behaviour. All the studied traits are repeatable or heritable on individual level, and their expression is either directly or indirectly sensitive to changes in vole numbers. My work considers a number of important fitness-related traits showing phenotypic plasticity in all of them. Further, in two chapters I show that there is individual variation in the amount of plasticity exhibited. These findings on plasticity in reproduction related traits suggest that variable environments indeed promote plasticity.

Spatial and temporal determinants of Finnish farmland bird populations

During the past decades agricultural intensification has caused dramatic population declines in a wide range of taxa related to farmland habitats, including farmland birds. In this thesis, I studied how boreal farmland landscape characteristics and agricultural land use affect the abundance and diversity of farmland birds using extensive field data collected by territory mapping of breeding farmland birds in various parts of Finland. My results show that the area and openness of agricultural areas are key determinants of farmland bird abundance and distribution. A landscape composition with enough open farmland combined with key habitats such as farmyards and wetland is likely to provide essential prerequisites for the occurrence of a rich farmland avifauna. In Finland, the majority of large areas suitable for open habitat specialists are located in southern and western parts of the country. However, the diversity of the species with an unfavourable conservation status in Europe (SPECs) had notable hotspot areas in northern and north-western agricultural areas. I found that in boreal agroecosystems farmland birds favour fields with springtime vegetative cover, especially agricultural grasslands and set-asides. Hence, in the spring cereal dominated Finnish agroecosystems it is the absence of field vegetation that may limit populations of many farmland bird species. It is likely that the decrease of crops providing vegetative cover in the spring, such as permanent grasslands, cultivated grass, and autumn-sown cereals, has greatly contributed to the declines of Finnish farmland birds. Grass crops have persistently declined in Finland as a consequence of specialization in crop production and the large-scale decline in livestock husbandry. Small-scale non-crop habitats, especially ditches and ditch margins, are also important for many bird species in the Finnish agroecosystems, but have dramatically declined during the last decades. A major problem for farmland bird conservation in Finland is the conflict between landscape structure and agricultural management. Areas with mixed and cattle farming are virtually absent from the large agricultural plains of southern and south-western Finland, where the landscape structure is more likely to be favourable for rich farmland bird assemblages. On the other hand, mixed and cattle farming is still rather frequent in northern and central parts of the country, where the landscape structure is not suitable for many farmland specialist birds requiring open landscapes. My results provide useful guidelines for farmland bird conservation, and imply that considerable attention needs to be paid to landscape factors when selecting areas for various conservational management actions, such as agri-environment schemes. Actions promoting the abundance of set-asides, grass crops, and ditches would markedly benefit Finnish farmland bird populations. Organic farming may benefit farmland birds, but it is not clear how general its beneficial effect is in boreal agroecosystems. The most urgent action aiming to preserve farmland biodiversity would be to support re-introducing and sustaining cattle farming by environmental subsidies. This would be especially beneficial in the southern parts of Finland, where the landscape characteristics and abundance of agricultural areas are most suitable for farmland birds and where cattle farming is currently rare.

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.

Ecological impacts of Phlebiopsis gigantea biocontrol treatment against Heterobasidion spp. as revealed by fungal community profiling and population analyses

Wood decay fungi belonging to the species complex Heterobasidion annosum sensu lato are among the most common and economically important species causing root rot and stem decay in conifers of the northern temperate regions. New infections by these pathogens can be suppressed by tree stump treatments using chemical or biological control agents. In Finland, the corticiaceous fungus Phlebiopsis gigantea has been formulated into a commercial biocontrol agent called Rotstop (Verdera Ltd.). This thesis addresses the ecological impacts of Rotstop biocontrol treatment on the mycoflora of conifer stumps. Locally, fungal communities within Rotstop-treated and untreated stumps were analyzed using a novel method based on DGGE profiling of small subunit ribosomal DNA fragments amplified directly from wood samples. Population analyses for P. gigantea and H. annosum s.l. were conducted to evaluate possible risks associated with local and/or global distribution of the Rotstop strain. Based on molecular community profiling by DGGE, we detected a few individual wood-inhabiting fungal species (OTUs) that seemed to have suffered or benefited from the Rotstop biocontrol treatment. The DGGE analyses also revealed fungal diversity not retrieved by cultivation and some fungal sequence types untypical for decomposing conifer wood. However, statistical analysis of DGGE community profiles obtained from Rotstop-treated and untreated conifer stumps revealed that the Rotstop treatment had not caused a statistically significant reduction in the species diversity of wood-inhabiting fungi within our experimental forest plots. Locally, ISSR genotyping of cultured P. gigantea strains showed that the Rotstop biocontrol strain was capable of surviving up to six years within treated Norway spruce stumps, while in Scots pine stumps it was sooner replaced by successor fungal species. In addition, the spread of resident P. gigantea strains into Rotstop-treated forest stands seemed effective in preventing the formation of genetically monomorphic populations in the short run. On a global scale, we detected a considerable level of genetic differentiation between the interfertile European and North American populations of P. gigantea. These results strongly suggest that local biocontrol strains should be used in order to prevent global spread of P. gigantea and hybrid formation between geographically isolated populations. The population analysis for H. annosum s.l. revealed a collection of Chinese fungal strains that showed a high degree of laboratory fertility with three different allopatric H. annosum s.l. taxa. However, based on the molecular markers, the Chinese strains could be clearly affiliated with the H. parviporum taxonomical cluster, which thus appears to have a continuous distribution range from Europe through southern Siberia to northern China. Keywords: Rotstop, wood decay, DGGE, ISSR fingerprinting, ribosomal DNA

Population dynamics of blue mussels in a variable environment at the edge of their range

Spatial and temporal variation in the abundance of species can often be ascribed to spatial and temporal variation in the surrounding environment. Knowledge of how biotic and abiotic factors operate over different spatial and temporal scales in determining distribution, abundance, and structure of populations lies at the heart of ecology. The major part of the current ecological theory stems from studies carried out in central parts of the distributional range of species, whereas knowledge of how marginal populations function is inadequate. Understanding how marginal populations, living at the edge of their range, function is however in a key position to advance ecology and evolutionary biology as scientific disciplines. My thesis focuses on the factors affecting dynamics of marginal populations of blue mussels (Mytilus edulis) living close to their tolerance limits with regard to salinity. The thesis aims to highlight the dynamics at the edge of the range and contrast these with dynamics in more central parts of the range in order to understand the potential interplay between the central and the marginal part in the focal system. The objectives of the thesis are approached by studies on: (1) factors affecting regional patterns of the species, (2) long-term temporal dynamics of the focal species spaced along a regional salinity gradient, (3) selective predation by increasing populations of roach (Rutilus rutilus) when feeding on their main food item, the blue mussel, (4) the primary and secondary effects of local wave exposure gradients and (5) the role of small-scale habitat heterogeneity as determinants of large-scale pattern. The thesis shows that populations of blue mussels are largely determined by large scale changes in sea water salinity, affecting mainly recruitment success and longevity of local populations. In opposite to the traditional view, the thesis strongly indicate that vertebrate predators strongly affect abundance and size structure of blue mussel populations, and that the role of these predators increases towards the margin where populations are increasingly top-down controlled. The thesis also indicates that the positive role of biogenic habitat modifiers increases towards the marginal areas, where populations of blue mussels are largely recruitment limited. Finally, the thesis shows that local blue mussel populations are strongly dependent on high water turbulence, and therefore, dense populations are constrained to offshore habitats. Finally, the thesis suggests that ongoing sedimentation of rocky shores is detrimental for the species, affecting recruitment success and post-recruit survival, pushing stable mussel beds towards offshore areas. Ongoing large scale changes in the Baltic Sea, especially dilution processes with attendant effects, are predicted to substantially contract the distributional range of the mussel, but also affect more central populations. The thesis shows that in order to understand the functioning of marginal populations, research should (1) strive for multi-scale approaches in order to link ecosystem patterns with ecosystem processes, and (2) challenge the prevailing tenets that origin from research carried out in central areas that may not be valid at the edge.

Conceptual and statistical modelling of environmental effects in population dynamics

Population dynamics are generally viewed as the result of intrinsic (purely density dependent) and extrinsic (environmental) processes. Both components, and potential interactions between those two, have to be modelled in order to understand and predict dynamics of natural populations; a topic that is of great importance in population management and conservation. This thesis focuses on modelling environmental effects in population dynamics and how effects of potentially relevant environmental variables can be statistically identified and quantified from time series data. Chapter I presents some useful models of multiplicative environmental effects for unstructured density dependent populations. The presented models can be written as standard multiple regression models that are easy to fit to data. Chapters II IV constitute empirical studies that statistically model environmental effects on population dynamics of several migratory bird species with different life history characteristics and migration strategies. In Chapter II, spruce cone crops are found to have a strong positive effect on the population growth of the great spotted woodpecker (Dendrocopos major), while cone crops of pine another important food resource for the species do not effectively explain population growth. The study compares rate- and ratio-dependent effects of cone availability, using state-space models that distinguish between process and observation error in the time series data. Chapter III shows how drought, in combination with settling behaviour during migration, produces asymmetric spatially synchronous patterns of population dynamics in North American ducks (genus Anas). Chapter IV investigates the dynamics of a Finnish population of skylark (Alauda arvensis), and point out effects of rainfall and habitat quality on population growth. Because the skylark time series and some of the environmental variables included show strong positive autocorrelation, the statistical significances are calculated using a Monte Carlo method, where random autocorrelated time series are generated. Chapter V is a simulation-based study, showing that ignoring observation error in analyses of population time series data can bias the estimated effects and measures of uncertainty, if the environmental variables are autocorrelated. It is concluded that the use of state-space models is an effective way to reach more accurate results. In summary, there are several biological assumptions and methodological issues that can affect the inferential outcome when estimating environmental effects from time series data, and that therefore need special attention. The functional form of the environmental effects and potential interactions between environment and population density are important to deal with. Other issues that should be considered are assumptions about density dependent regulation, modelling potential observation error, and when needed, accounting for spatial and/or temporal autocorrelation.

Human genetic variation in the Baltic Sea region: Features of population history and natural selection

In this thesis, the genetic variation of human populations from the Baltic Sea region was studied in order to elucidate population history as well as evolutionary adaptation in this region. The study provided novel understanding of how the complex population level processes of migration, genetic drift, and natural selection have shaped genetic variation in North European populations. Results from genome-wide, mitochondrial DNA and Y-chromosomal analyses suggested that the genetic background of the populations of the Baltic Sea region lies predominantly in Continental Europe, which is consistent with earlier studies and archaeological evidence. The late settlement of Fennoscandia after the Ice Age and the subsequent small population size have led to pronounced genetic drift, especially in Finland and Karelia but also in Sweden, evident especially in genome-wide and Y-chromosomal analyses. Consequently, these populations show striking genetic differentiation, as opposed to much more homogeneous pattern of variation in Central European populations. Additionally, the eastern side of the Baltic Sea was observed to have experienced eastern influence in the genome-wide data as well as in mitochondrial DNA and Y-chromosomal variation – consistent with linguistic connections. However, Slavic influence in the Baltic Sea populations appears minor on genetic level. While the genetic diversity of the Finnish population overall was low, genome-wide and Y-chromosomal results showed pronounced regional differences. The genetic distance between Western and Eastern Finland was larger than for many geographically distant population pairs, and provinces also showed genetic differences. This is probably mainly due to the late settlement of Eastern Finland and local isolation, although differences in ancestral migration waves may contribute to this, too. In contrast, mitochondrial DNA and Y-chromosomal analyses of the contemporary Swedish population revealed a much less pronounced population structure and a fusion of the traces of ancient admixture, genetic drift, and recent immigration. Genome-wide datasets also provide a resource for studying the adaptive evolution of human populations. This study revealed tens of loci with strong signs of recent positive selection in Northern Europe. These results provide interesting targets for future research on evolutionary adaptation, and may be important for understanding the background of disease-causing variants in human populations.