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.

The community of medium-sized carnivores: the interactions between species, habitats and rabies

The composition of the carnivore community influences the different forms of inter-specific interactions. Furthermore, inter-specific interactions of carnivores have important implications for intra-guild competition, epidemiology and strategies of species-specific population management. Zoonooses, such as rabies, are diseases that can be transmitted from wildlife to people. Knowing the ecological characteristics of the species helps us to choose the right preventive actions and to time them accurately. In this thesis, I have studied how raccoon dogs Nyctereutes procyonoides, European badgers Meles meles, red foxes Vulpes vulpes and domestic cats Felis silvestris catus act as members of carnivore community, and how these interactions relate to the transmission risk of rabies. In the study area, these species form a community of medium-sized and rather generalist predators. They live in the same areas, in spatially and temporally overlapping home ranges and use the same habitats and dens and even have similar diets. However, there is no direct evidence of competition. Shared dens point to good tolerance of other species. Numerous observations of animals moving in each other’s proximity give similar clues. However, overlapping home ranges and similar habitat preferences lead to frequent inter-specific contacts, which increase the risk of possible rabies transmission. Also, the new insight of habitat use gained by this study illustrates the similar favouring of deciduous forests and fields by these sympatric medium-sized carnivores, creating a basis for contact zones, i.e. risky habitats for rabies transmission and spread. This study is so far the only simultaneous radio tracking study of raccoon dogs, badgers, foxes and cats. These results give new insight of the interactions in the carnivore community, as well as of the behaviour of each individual species. Also, these results have significant implications for the planning of rabies control. In order to reach viable management decisions, not only one or two species should be taken into consideration, but the whole community. In particular, this changes the perspective to inter-specific contacts, animal densities, densities of individuals susceptible to diseases and the magnitude of preventive actions. Rabies should be considered as a multi-vector disease, at least in Finland and the Baltic states. It is of interest for disease management to be able to model an epizootic with local parameters to reflect the real situation and also to suite best the local management needs.

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.

Archaea in the Mycorrhizosphere of Boreal Forest Trees

Archaea were long thought to be a group of ancient bacteria, which mainly lived in extreme environments. Due to the development of DNA sequencing methods and molecular phylogenetic analyses, it was shown that the living organisms are in fact divided into three domains; the Archaea, Bacteria and the Eucarya. Since the beginning of the previous decade, it was shown that archaea generally inhabit moderate environments and that these non-extremophilic archaea are more ubiquitous than the extremophiles. Group 1 of non-extreme archaea affiliate with the phylum Crenarchaeota. The most commonly found soil archaea belong to the subgroup 1.1b. However, the Crenarchaeota found in the Fennoscandian boreal forest soil belong to the subgroup 1.1c. The organic top layer of the boreal forest soil, the humus, is dominated by ectomycorrhizal fungal hyphae. These colonise virtually all tree fine root tips in the humus layer and have been shown to harbour distinct bacterial populations different from those in the humus. The archaea have also been shown to colonise both boreal forest humus and the rhizospheres of plants. In this work, studies on the archaeal communities in the ectomycorrhizospheres of boreal forest trees were conducted in microcosms. Archaea belonging to the group 1.1c Crenarchaeota and Euryarchaeota of the genera Halobacterium and Methanolobus were detected. The archaea generally colonised fungal habitats, such as ectomycorrhizas and external mycelia, rather than the non-mycorrhizal fine roots of trees. The species of ectomycorrhizal fungus had a great impact on the archaeal community composition. A stable euryarchaeotal community was detected especially in the mycorrhizas, of most of the tested Scots pine colonising ectomycorrhizal fungi. The Crenarchaeota appeared more sporadically in these habitats, but had a greater diversity than the Euryarchaeota. P. involutus mycorrhizas had a higher diversity of 1.1c Crenarchaeota than the other ectomycorrhizal fungi. The detection level of archaea in the roots of boreal trees was generally low although archaea have been shown to associate with roots of different plants. However, alder showed a high diversity of 1.1c Crenarchaeota, exceeding that of any of the tested mycorrhizas. The archaeal 16S rRNA genes detected from the non-mycorrhizal roots were different from those of the P. involutus mycorrhizas. In the phylogenetic analyses, the archaeal 16S rRNA gene sequences obtained from non-mycorrhizal fine roots fell in a separate cluster within the group 1.1c Crenarchaeota than those from the mycorrhizas. When the roots of the differrent tree species were colonised by P. involutus, the diversity and frequency of the archaeal populations of the different tree species were more similar to each other. Both Cren- and Euryarchaeota were enriched in cultures to which C-1 substrates were added. The 1.1c Crenarchaeota grew anaerobically in mineral medium with CH4 and CO2 as the only available C sources, and in yeast extract media with CO2 and CH4 or H2. The crenarchaeotal diversity was higher in aerobic cultures on mineral medium with CH4 or CH3OH than in the anaerobic cultures. Ecological functions of the mycorrhizal 1.1c Crenarchaeota in both anaerobic and aerobic cycling of C-1 compounds were indicated. The phylogenetic analyses did not divide the detected Crenarchaeota into anaerobic and aerobic groups. This may suggest that the mycorrhizospheric crenarchaeotal communities consist of closely related groups of anaerobic and aerobic 1.1c Crenarchaeota, or the 1.1c Crenarchaeota may be facultatively anaerobic. Halobacteria were enriched in non-saline anaerobic yeast extract medium cultures in which CH4 was either added or produced, but were not detected in the aerobic cultures. They may potentially be involved in anaerobic CH4 cycling in ectomycorrhizas. The CH4 production of the mycorrhizal samples was over 10 times higher than for humus devoid of mycorrhizal hyphae, indicating a high CH4 production potential of the mycorrhizal metanogenic community. Autofluorescent methanogenic archaea were detected by microscopy and 16S rRNA gene sequences of the genus Methanolobus were obtained. The archaeal community depended on both tree species and the type of ectomycorrhizal fungus colonising the roots and the Cren- and Euryarchaeota may have different ecological functions in the different parts of the boreal forest tree rhizosphere and mycorrhizosphere. By employing the results of this study, it may be possible to isolate both 1.1c Crenarchaeota as well as non-halophilic halobacteria and aerotolerant methanogens from mycorrhizospheres. These archaea may be used as indicators for change in the boreal forest soil ecosystem due to different factors, such as exploitations of forests and the rise in global temperature. More information about the microbial populations with apparently low cell numbers but significant ecological impacts, such as the boreal forest soil methanogens, may be of crucial importance to counteract human impacts on such globally important ecosystems as the boreal forests.

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.

The high- and low-activity forms of human plasma phospholipid transfer protein (PLTP)

Plasma phospholipid transfer protein (PLTP) plays a crucial role in high-density lipoprotein (HDL) metabolism and reverse cholesterol transport (RCT). It mediates the generation of pre-beta-HDL particles, enhances the cholesterol efflux from peripheral cells to pre-beta-HDL, and metabolically maintains the plasma HDL levels by facilitating the transfer of post-lipolytic surface remnants of triglyceride-rich lipoproteins to HDL. In addition to the antiatherogenic properties, recent findings indicate that PLTP has also proatherogenic characteristics, and that these opposite characteristics of PLTP are dependent on the site of PLTP expression and action. In human plasma, PLTP exists in a high-activity (HA-PLTP) and a low-activity form (LA-PLTP), which are associated with macromolecular complexes of different size and composition. The aims of this thesis were to isolate the two PLTP forms from human plasma, to characterize the molecular complexes in which the HA- and LA-PLTP reside, and to study the interactions of the PLTP forms with apolipoproteins (apo) and the ability of apolipoproteins to regulate PLTP activity. In addition, we aimed to study the distribution of the two PLTP forms in a Finnish population sample as well as to find possible regulatory factors for PLTP by investigating the influence of lipid and glucose metabolism on the balance between the HA- and LA-PLTP. For these purposes, an enzyme-linked immunosorbent assay (ELISA) capable of determining the serum total PLTP concentration and quantitating the two PLTP forms separately was developed. In this thesis, it was demonstrated that the HA-PLTP isolated from human plasma copurified with apoE, whereas the LA-PLTP formed a complex with apoA-I. The separation of these two PLTP forms was carried out by a dextran sulfate (DxSO4)-CaCl2 precipitation of plasma samples before the mass determination. A similar immunoreactivity of the two PLTP forms in the ELISA could be reached after a partial sample denaturation by SDS. Among normolipidemic Finnish individuals, the mean PLTP mass was 6.6 +/- 1.5 mg/l and the mean PLTP activity 6.6 +/- 1.7 umol/ml/h. Of the serum PLTP concentration, almost 50% represented HA-PLTP. The results indicate that plasma HDL levels could regulate PLTP concentration, while PLTP activity could be regulated by plasma triglyceride-rich very low-density lipoprotein (VLDL) concentration. Furthermore, new evidence is presented that PLTP could also play a role in glucose metabolism. Finally, both PLTP forms were found to interact with apoA-I, apoA-IV, and apoE. In addition, both apoE and apoA-IV, but not apoA-I, were capable of activating the LA-PLTP. These findings suggest that the distribution of the HA- and LA-PLTP in human plasma is subject to dynamic regulation by apolipoproteins.