Effects of environmental variation on ecological and evolutionary dynamics

Environmental variation is a fact of life for all the species on earth: for any population of any particular species, the local environmental conditions are liable to vary in both time and space. In today's world, anthropogenic activity is causing habitat loss and fragmentation for many species, which may profoundly alter the characteristics of environmental variation in remaining habitat. Previous research indicates that, as habitat is lost, the spatial configuration of remaining habitat will increasingly affect the dynamics by which populations are governed. Through the use of mathematical models, this thesis asks how environmental variation interacts with species properties to influence population dynamics, local adaptation, and dispersal evolution. More specifically, we couple continuous-time continuous-space stochastic population dynamic models to landscape models. We manipulate environmental variation via parameters such as mean patch size, patch density, and patch longevity. Among other findings, we show that a mixture of high and low quality habitat is commonly better for a population than uniformly mediocre habitat. This conclusion is justified by purely ecological arguments, yet the positive effects of landscape heterogeneity may be enhanced further by local adaptation, and by the evolution of short-ranged dispersal. The predicted evolutionary responses to environmental variation are complex, however, since they involve numerous conflicting factors. We discuss why the species that have high levels of local adaptation within their ranges may not be the same species that benefit from local adaptation during range expansion. We show how habitat loss can lead to either increased or decreased selection for dispersal depending on the type of habitat and the manner in which it is lost. To study the models, we develop a recent analytical method, Perturbation expansion, to enable the incorporation of environmental variation. Within this context, we use two methods to address evolutionary dynamics: Adaptive dynamics, which assumes mutations occur infrequently so that the ecological and evolutionary timescales can be separated, and via Genotype distributions, which assume mutations are more frequent. The two approaches generally lead to similar predictions yet, exceptionally, we show how the evolutionary response of dispersal behaviour to habitat turnover may qualitatively depend on the mutation rate.

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.

Molecular genetics of Meckel syndrome : Ciliary genes are defective in MKS

Meckel syndrome (MKS, MIM 249000) is a severe developmental disorder that leads to death already in utero or shortly after birth. MKS diagnosis can be established by a careful ultrasound examination already at 11-14 weeks of gestation. The main features of MKS are occipital meningoencephalocele, cystic kidney dysplasia and fibrotic changes of the liver. In addition, polydactyly is frequently reported in the cases. The aim of the study was to characterize the molecular and functional defects in MKS. In this study we were able to identify two major MKS mutations in Finnish population, which cover over 90% of the cases. The first mutation is a 29 bp intronic deletion in the MKS1 gene (c.1483-7_35del) that is found in 70% of the families and the second is a C>T substitution in the coding region of CC2D2A (c.1762C>T), that is found in 20% of the MKS families. Both of these mutations result in abnormal splicing. The discovery of the disease genes has revealed that MKS is caused by primary cilia dysfunction. MKS1 gene has a conserved B9 domain, and it is found in the predicted ciliary proteome. CC2D2A protein is also found in the predicted ciliary proteome and it has a Ca2+ binding domain. The number of genes behind MKS has increased rapidly in the past years and to date, mutations have been identified in five genes (MKS1, TMEM67/MKS3, CEP290/MKS4, RPGRIP1L/MKS5 and CC2D2A/MKS6). Identification of the disease genes mutations has also revealed that MKS is an allelic disorder with other syndromes with overlapping phenotypes. Disorders that are caused by primary cilia dysfunction are collectively known as ciliopathies. Sequence analysis of all the known MKS genes in Finnish and non-Finnish families available to us, where the mutation was still unknown, revealed mutations in 14 out of the 30 families included in the study. When we collected all the reported mutations in MKS genes in different syndromes we could see that there was clearly a genotype-syndrome correlation between the mutations and the syndromes, since the same pair of mutations has never been reported in different syndromes. The basic molecular events behind MKS will not only give us information of this syndrome, but also significant novel information on early fetal development in general.

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.

Communities of wood-inhabiting fungi : Ecological requirements and responses to forest management and fragmentation

Despite much research on forest biodiversity in Fennoscandia, the exact mechanisms of species declines in dead-wood dependent fungi are still poorly understood. In particular, there is only limited information on why certain fungal species have responded negatively to habitat loss and fragmentation, while others have not. Understanding the mechanisms behind species declines would be essential for the design and development of ecologically effective and scientifically informed conservation measures, and management practices that would promote biodiversity in production forests. In this thesis I study the ecology of polypores and their responses to forest management, with a particular focus on why some species have declined more than others. The data considered in the thesis comprise altogether 98,318 dead-wood objects, with 43,085 observations of 174 fungal species. Out of these, 1,964 observations represent 58 red-listed species. The data were collected from 496 sites, including woodland key habitats, clear-cuts with retention trees, mature managed forests, and natural or natural-like forests in southern Finland and Russian Karelia. I show that the most relevant way of measuring resource availability can differ to a great extent between species seemingly sharing the same resources. It is thus critical to measure the availability of resources in a way that takes into account the ecological requirements of the species. The results show that connectivity at the local, landscape and regional scales is important especially for the highly specialized species, many of which are also red-listed. Habitat loss and fragmentation affect not only species diversity but also the relative abundances of the species and, consequently, species interactions and fungal successional pathways. Changes in species distributions and abundances are likely to affect the food chains in which wood-inhabiting fungi are involved, and thus the functioning of the whole forest ecosystem. The findings of my thesis highlight the importance of protecting well-connected, large and high-quality forest areas to maintain forest biodiversity. Small habitat patches distributed across the landscape are likely to contribute only marginally to protection of red-listed species, especially if habitat quality is not substantially higher than in ordinary managed forest, as is the case with woodland key habitats. Key habitats might supplement the forest protection network if they were delineated larger and if harvesting of individual trees was prohibited in them. Taking the landscape perspective into account in the design and development of conservation measures is critical while striving to halt the decline of forest biodiversity in an ecologically effective manner.

Evaluation of methods and applications for behavioural profiling of transgenic mice

During the last 10-15 years interest in mouse behavioural analysis has evolved considerably. The driving force is development in molecular biological techniques that allow manipulation of the mouse genome by changing the expression of genes. Therefore, with some limitations it is possible to study how genes participate in regulation of physiological functions and to create models explaining genetic contribution to various pathological conditions. The first aim of our study was to establish a framework for behavioural phenotyping of genetically modified mice. We established comprehensive battery of tests for the initial screening of mutant mice. These included tests for exploratory and locomotor activity, emotional behaviour, sensory functions, and cognitive performance. Our interest was in the behavioural patterns of common background strains used for genetic manipulations in mice. Additionally we studied the behavioural effect of sex differences, test history, and individual housing. Our findings highlight the importance of careful consideration of genetic background for analysis of mutant mice. It was evident that some backgrounds may mask or modify the behavioural phenotype of mutants and thereby lead to false positive or negative findings. Moreover, there is no universal strain that is equally suitable for all tests, and using different backgrounds allows one to address possible phenotype modifying factors. We discovered that previous experience affected performance in several tasks. The most sensitive traits were the exploratory and emotional behaviour, as well as motor and nociceptive functions. Therefore, it may be essential to repeat some of the tests in naïve animals for assuring the phenotype. Social isolation for a long time period had strong effects on exploratory behaviour, but also on learning and memory. All experiments revealed significant interactions between strain and environmental factors (test history or housing condition) indicating genotype-dependent effects of environmental manipulations. Several mutant line analyses utilize this information. For example, we studied mice overexpressing as well as those lacking extracellular matrix protein heparin-binding growth-associated molecule (HB-GAM), and mice lacking N-syndecan (a receptor for HB-GAM). All mutant mice appeared to be fertile and healthy, without any apparent neurological or sensory defects. The lack of HB-GAM and N-syndecan, however, significantly reduced the learning capacity of the mice. On the other hand, overexpression of HB-GAM resulted in facilitated learning. Moreover, HB-GAM knockout mice displayed higher anxiety-like behaviour, whereas anxiety was reduced in HB-GAM overexpressing mice. Changes in hippocampal plasticity accompanied the behavioural phenotypes. We conclude that HB-GAM and N-syndecan are involved in the modulation of synaptic plasticity in hippocampus and play a role in regulation of anxiety- and learning-related behaviour.

Bioavailability and toxicity of chemicals in inland waters : the importance of prevailing water chemistry and implications for risk assessment

Physicochemical characterization of freshwater samples from Finland, Sweden, the Netherlands, and Spain revealed that water hardness and pH decreased and the quantity and quality of humic substances changed considerably in this geographical series from south to north. Since the ambient water chemistry may affect the availability of chemicals, the total aqueous concentration of a chemical may be insufficient to predict the bioconcentration, subsequent biological response, and thus risk. In addition, organisms could be affected directly by water quality characteristics. In this context the main objective of this thesis was to investigate the bioavailability of selected ecotoxicologically relevant chemicals (cadmium, benzo(a)pyrene, and pyrene) in various European surface waters and to show the importance of certain water chemistry characteristics in interpreting the bioavailability and toxicity results. The bioavailability of cadmium to Daphnia magna was examined in very soft humic lake water. Humic substances as natural ligands decreased the free and bioavailable proportion of cadmium in soft lake water. As a consequence the uptake rate and the acute toxicity decreased compared with the humic-free reference. When the hardness of humic lake water was artificially elevated, the acute toxicity of cadmium decreased, although the proportion of free cadmium increased. The decreased bioavailability of cadmium in hard water was a result of effective competition for uptake by the hardness cations, especially calcium ions. The protective role of humic substances and water hardness against cadmium toxicity was also observed in Lumbriculus variegatus, although D. magna was more sensitive to cadmium. The bioavailability of two polycyclic aromatic hydrocarbons (PAHs), pyrene and benzo(a)pyrene, was studied in European surface waters of varying water chemistry. Humic substances acted as complexing ligands with both PAHs, but the bioavailability of the more lipophilic benzo(a)pyrene to D. magna was affected more by humic substances than that of pyrene. In addition, not only the quantity of humic substances, but also their quality affected the bioavailability of benzo(a)pyrene. Nevertheless, the humic substances played a protective role in the photo-enhanced toxicity of pyrene under UV-B radiation. Water hardness had no effect on pyrene toxicity. Results indicate that the typical physicochemical characteristics of boreal freshwaters should be considered carefully in local and regional risk assessment of chemicals concerning the Fennoscandian region.

Harmful algae in the planktonic food web of the Baltic Sea

This study deals with algal species occurring commonly in the Baltic Sea: haptophyte Prymnesium parvum, dinoflagellates Dinophysis acuminata, D. norvegica and D. rotundata, and cyanobacterium Nodularia spumigena. The hypotheses are connected to the toxicity of the species, to the factors determining toxicity, to the consequences of toxicity and to the transfer of toxins in the aquatic food web. Since the Baltic Sea is severely eutrophicated, the fast-growing haptophytes have potential in causing toxic blooms. In our studies, the toxicity (as haemolytic activity) of the haptophyte P. parvum was highest under phosphorus-limited conditions, but the cells were toxic also under nitrogen limitation and under nutrient-balanced growth conditions. The cellular nutrient ratios were tightly related to the toxicity. The stoichiometric flexibility for cellular phosphorus quota was higher than for nitrogen, and nitrogen limitation led to decreased biomass. Negative allelopathic effects on another algae (Rhodomonas salina) could be observed already at low P. parvum cell densities, whereas immediate lysis of R. salina cells occurred at P. parvum cell densities corresponding to natural blooms. Release of dissolved organic carbon from the R. salina cells was measured within 30 minutes, and an increase in bacterial number and biomass was measured within 23 h. Because of the allelopathic effect, formation of a P. parvum bloom may accelerate after a critical cell density is reached and the competing species are eliminated. A P. parvum bloom indirectly stimulates bacterial growth, and alters the functioning of the planktonic food web by increasing the carbon transfer through the microbial loop. Our results were the first reports on DSP toxins in Dinophysis cells in the Gulf of Finland and on PTX-2 in the Baltic Sea. Cellular toxin contents in Dinophysis spp. ranged from 0.2 to 149 pg DTX-1 cell-1 and from 1.6 to 19.9 pg PTX-2 cell-1 in the Gulf of Finland. D. norvegica was found mainly around the thermocline (max. 200 cells L-1), whereas D. acuminata was found in the whole mixed layer (max. 7 280 cells L-1). Toxins in the sediment trap corresponded to 1 % of DTX-1 and 0.01 % PTX-2 of the DSP pool in the suspended matter. This indicates that the majority of the DSP toxins does not enter the benthic community, but is either decomposed in the water column, or transferred to higher trophic levels in the planktonic food chain. We found that nodularin, produced by Nodularia spumigena, was transferred to the copepod Eurytemora affinis through three pathways: by grazing on filaments of small Nodularia, directly from the dissolved pool, and through the microbial food web by copepods grazing on ciliates, dinoflagellates and heterotrophic nanoflagellates. The estimated proportion of the microbial food web in nodularin transfer was 22-45 % and 71-76 % in our two experiments, respectively. This highlights the potential role of the microbial food web in the transfer of toxins in the planktonic food web.

DNA Packaging and Host Cell Lysis: Late Events in Bacteriophage PRD1 Infection

The object of this study is a tailless internal membrane-containing bacteriophage PRD1. It has a dsDNA genome with covalently bound terminal proteins required for replication. The uniqueness of the structure makes this phage a desirable object of research. PRD1 has been studied for some 30 years during which time a lot of information has accumulated on its structure and life-cycle. The two least characterised steps of the PRD1 life-cycle, the genome packaging and virus release are investigated here. PRD1 shares the main principles of virion assembly (DNA packaging in particular) and host cell lysis with other dsDNA bacteriophages. However, this phage has some fascinating individual peculiarities, such as DNA packaging into a membrane vesicle inside the capsid, absence of apparent portal protein, holin inhibitor and procapsid expansion. In the course of this study we have identified the components of the DNA packaging vertex of the capsid, and determined the function of protein P6 in packaging. We managed to purify the procapsids for an in vitro packaging system, optimise the reaction and significantly increase its efficiency. We developed a new method to determine DNA translocation and were able to quantify the efficiency and the rate of packaging. A model for PRD1 DNA packaging was also proposed. Another part of this study covers the lysis of the host cell. As other dsDNA bacteriophages PRD1 has been proposed to utilise a two-component lysis system. The existence of this lysis system in PRD1 has been proven by experiments using recombinant proteins and the multi-step nature of the lysis process has been established.

Effect of long-wave UV radiation on mouse melanoma: an in vitro and in vivo study

The skin cancer incidence has increased substantially over the past decades and the role of ultraviolet (UV) radiation in the etiology of skin cancer is well established. Ultraviolet B radiation (280-320 nm) is commonly considered as the more harmful part of the UV-spectrum due to its DNA-damaging potential and well-known carcinogenic effects. Ultraviolet A radiation (320-400 nm) is still regarded as a relatively low health hazard. However, UVA radiation is the predominant component in sunlight, constituting more than 90% of the environmentally relevant solar ultraviolet radiation. In the light of the recent scientific evidence, UVA has been shown to have genotoxic and immunologic effects, and it has been proposed that UVA plays a significant role in the development of skin cancer. Due to the popularity of skin tanning lamps, which emit high intensity UVA radiation and because of the prolonged sun tanning periods with the help of effective UVB blockers, the potential deleterious effects of UVA has emerged as a source of concern for public health. The possibility that UV radiation may affect melanoma metastasis has not been addressed before. UVA radiation can modulate various cellular processes, some of which might affect the metastatic potential of melanoma cells. The aim of the present study was to investigate the possible role of UVA irradiation on the metastatic capacity of mouse melanoma both in vitro and in vivo. The in vitro part of the study dealt with the enhancement of the intercellular interactions occurring either between tumor cells or between tumor cells and endothelial cells after UVA irradiation. The use of the mouse melanoma/endothelium in vitro model showed that a single-dose of UVA to melanoma cells causes an increase in melanoma cell adhesiveness to non-irradiated endothelium after 24-h irradiation. Multiple-dose irradiation of melanoma cells already increased adhesion at a 1-h time-point, which suggests the possible cumulative effect of multiple doses of UVA irradiation. This enhancement of adhesiveness might lead to an increase in binding tumor cells to the endothelial lining of vasculature in various internal organs if occurring also in vivo. A further novel observation is that UVA induced both decline in the expression of E-cadherin adhesion molecule and increase in the expression of the N-cadherin adhesion molecule. In addition, a significant decline in homotypic melanoma-melanoma adhesion (clustering) was observed, which might result in the reduction of E-cadherin expression. The aim of the in vivo animal study was to confirm the physiological significance of previously obtained in vitro results and to determine whether UVA radiation might increase melanoma metastasis in vivo. The use of C57BL/6 mice and syngeneic melanoma cell lines B16-F1 and B16-F10 showed that mice, which were i.v. injected with B16-F1 melanoma cells and thereafter exposed to UVA developed significantly more lung metastases when compared with the non-UVA-exposed group. To study the mechanism behind this phenomenon, the direct effect of UVA-induced lung colonization capacity was examined by the in vitro exposure of B16-F1 cells. Alternatively, the UVA-induced immunosuppression, which might be involved in increased melanoma metastasis, was measured by standard contact hypersensitivity assay (CHS). It appears that the UVA-induced increase of metastasis in vivo might be caused by a combination of UVA-induced systemic immunosuppression, and to the lesser extent, it might be caused by the increased adhesiveness of UVA irradiated melanoma cells. Finally, the UVA effect on gene expression in mouse melanoma was determined by a cDNA array, which revealed UVA-induced changes in the 9 differentially expressed genes that are involved in angiogenesis, cell cycle, stress-response, and cell motility. These results suggest that observed genes might be involved in cellular response to UVA and a physiologically relevant UVA dose have previously unknown cellular implications. The novel results presented in this thesis offer evidence that UVA exposure might increase the metastatic potential of the melanoma cells present in blood circulation. Considering the wellknown UVA-induced deleterious effects on cellular level, this study further supports the notion that UVA radiation might have more potential impact on health than previously suggested. The possibility of the pro-metastatic effects of UVA exposure might not be of very high significance for daily exposures. However, UVA effects might gain physiological significance following extensive sunbathing or solaria tanning periods. Whether similar UVA-induced pro-metastatic effects occur in people sunbathing or using solaria remains to be determined. In the light of the results presented in this thesis, the avoidance of solaria use could be well justified.

Missing-In-Metastasis (MIM) regulates cell morphology by promoting plasma membrane and actin cytoskeleton dynamics

The cells of multicellular organisms have differentiated to carry out specific functions that are often accompanied by distinct cell morphology. The actin cytoskeleton is one of the key regulators of cell shape subsequently controlling multiple cellular events including cell migration, cell division, endo- and exocytosis. A large set of actin regulating proteins has evolved to achieve and tightly coordinate this wide range of functions. Some actin regulator proteins have so-called house keeping roles and are essential for all eukaryotic cells, but some have evolved to meet the requirements of more specialized cell-types found in higher organisms enabling complex functions of differentiated organs, such as liver, kidney and brain. Often processes mediated by the actin cytoskeleton, like formation of cellular protrusions during cell migration, are intimately linked to plasma membrane remodeling. Thus, a close cooperation between these two cellular compartments is necessary, yet not much is known about the underlying molecular mechanisms. This study focused on a vertebrate-specific protein called missing-in-metastasis (MIM), which was originally characterized as a metastasis suppressor of bladder cancer. We demonstrated that MIM regulates the dynamics of actin cytoskeleton via its WH2 domain, and is expressed in a cell-type specific manner. Interestingly, further examination showed that the IM-domain of MIM displays a novel membrane tubulation activity, which induces formation of filopodia in cells. Following studies demonstrated that this membrane deformation activity is crucial for cell protrusions driven by MIM. In mammals, there are five members of IM-domain protein family. Functions and expression patterns of these family members have remained poorly characterized. To understand the physiological functions of MIM, we generated MIM knockout mice. MIM-deficient mice display no apparent developmental defects, but instead suffer from progressive renal disease and increased susceptibility to tumors. This indicates that MIM plays a role in the maintenance of specific physiological functions associated with distinct cell morphologies. Taken together, these studies implicate MIM both in the regulation of the actin cytoskeleton and the plasma membrane. Our results thus suggest that members of MIM/IRSp53 protein family coordinate the actin cytoskeleton:plasma membrane interface to control cell and tissue morphogenesis in multicellular organisms.

The Human UDP-Glucuronosyltransferases: Studies on Substrate Binding and Catalytic Mechanism

PDF

Where and how to conserve: Extending the scope of spatial reserve network design

Ongoing habitat loss and fragmentation threaten much of the biodiversity that we know today. As such, conservation efforts are required if we want to protect biodiversity. Conservation budgets are typically tight, making the cost-effective selection of protected areas difficult. Therefore, reserve design methods have been developed to identify sets of sites, that together represent the species of conservation interest in a cost-effective manner. To be able to select reserve networks, data on species distributions is needed. Such data is often incomplete, but species habitat distribution models (SHDMs) can be used to link the occurrence of the species at the surveyed sites to the environmental conditions at these locations (e.g. climatic, vegetation and soil conditions). The probability of the species occurring at unvisited location is next predicted by the model, based on the environmental conditions of those sites. The spatial configuration of reserve networks is important, because habitat loss around reserves can influence the persistence of species inside the network. Since species differ in their requirements for network configuration, the spatial cohesion of networks needs to be species-specific. A way to account for species-specific requirements is to use spatial variables in SHDMs. Spatial SHDMs allow the evaluation of the effect of reserve network configuration on the probability of occurrence of the species inside the network. Even though reserves are important for conservation, they are not the only option available to conservation planners. To enhance or maintain habitat quality, restoration or maintenance measures are sometimes required. As a result, the number of conservation options per site increases. Currently available reserve selection tools do however not offer the ability to handle multiple, alternative options per site. This thesis extends the existing methodology for reserve design, by offering methods to identify cost-effective conservation planning solutions when multiple, alternative conservation options are available per site. Although restoration and maintenance measures are beneficial to certain species, they can be harmful to other species with different requirements. This introduces trade-offs between species when identifying which conservation action is best applied to which site. The thesis describes how the strength of such trade-offs can be identified, which is useful for assessing consequences of conservation decisions regarding species priorities and budget. Furthermore, the results of the thesis indicate that spatial SHDMs can be successfully used to account for species-specific requirements for spatial cohesion - in the reserve selection (single-option) context as well as in the multi-option context. Accounting for the spatial requirements of multiple species and allowing for several conservation options is however complicated, due to trade-offs in species requirements. It is also shown that spatial SHDMs can be successfully used for gaining information on factors that drive a species spatial distribution. Such information is valuable to conservation planning, as better knowledge on species requirements facilitates the design of networks for species persistence. This methods and results described in this thesis aim to improve species probabilities of persistence, by taking better account of species habitat and spatial requirements. Many real-world conservation planning problems are characterised by a variety of conservation options related to protection, restoration and maintenance of habitat. Planning tools therefore need to be able to incorporate multiple conservation options per site, in order to continue the search for cost-effective conservation planning solutions. Simultaneously, the spatial requirements of species need to be considered. The methods described in this thesis offer a starting point for combining these two relevant aspects of conservation planning.

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

The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.