Living on the edge : Population genetics of Finno-Ugric-speaking humans in North Eurasia

The major aim of this thesis was to examine the origins and distribution of uniparental and autosomal genetic variation among the Finno-Ugric-speaking human populations living in Boreal and Arctic regions of North Eurasia. In more detail, I aimed to disentangle the underlying molecular and population genetic factors which have produced the patterns of uniparental and autosomal genetic diversity in these populations. Among Finno-Ugrics the genetic amalgamation and clinal distribution of West and East Eurasian gene pools were observed within uniparental markers. This admixture indicates that North Eurasia was colonized through Central Asia/ South Siberia by human groups already carrying both West and East Eurasian lineages. The complex combination of founder effects, gene flow and genetic drift underlying the genetic diversity of the Finno-Ugric- speaking populations were emphasized by low haplotype diversity within and among uniparental and biparental markers. A high prevalence of lactase persistence allele among the North Eurasian Finno- Ugric agriculturalist populations was also shown indicating a local adaptation to subsistence change with lactose rich diet. Moreover, the haplotype background of lactase persistence allele among the Finno- Ugric-speakers strongly suggested that the lactase persistence T-13910 mutation was introduced independently more than once to the North Eurasian gene pool. A significant difference in genetic diversity, haplotype structure and LD distribution within the cytochrome P450 CYP2C and CYP2D regions revealed the unique gene pool of the Finno-Ugric Saami created mainly by population genetic processes compared to other Europeans and sub-Saharan Mandenka population. From all studied populations the Saami showed also significantly the highest allele frequency of a CYP2C19 gene mutation causing variable drug reactions. The diversity patterns observed within CYP2C and CYP2D regions emphasize the strong effect of demographic history shaping genetic diversity and LD especially among such small and constant size populations as the Finno-Ugric-speaking Saami. Moreover, the increased LD in Saami due to genetic drift and/or admixture was shown to offer an advantage for further attempts to identify alleles associated to common complex pharmacogenetic traits.

Pharmacogenetics of SLCO1B1: Population genetics and effect on statins

Kohonneiden kolesterolipitoisuuksien alentamisessa käytettävien statiinien hyödyt sydän- ja verisuonisairauksien estossa on vahvasti osoitettu ja niiden käyttö on niin Suomessa kuin muuallakin maailmassa kasvanut voimakkaasti – Suomessa statiininkäyttäjiä on noin 600 000. Statiinilääkitys on pitkäaikaisessakin käytössä melko hyvin siedetty, mutta yleisimpinä haittavaikutuksina voi ilmetä lihasheikkoutta, -kipua ja -kramppeja, jotka voivat edetä jopa henkeä uhkaavaksi lihasvaurioksi. Lihashaittariski suurenee suhteessa statiiniannokseen ja plasman statiinipitoisuuksiin. Statiinien plasmapitoisuuksissa, tehossa ja haittavaikutusten ilmenemisessä on suuria potilaskohtaisia eroja. SLCO1B1-geenin koodaama OATP1B1-kuljetusproteiini kuljettaa monia elimistön omia aineita ja lääkeaineita verenkierrosta solukalvon läpi maksasoluun, mm. statiineja, joiden kolesterolia alentava vaikutus ja poistuminen elimistöstä tapahtuvat pääosin maksassa. Erään SLCO1B1-geenin nukleotidimuutoksen (c.521T>C) tiedetään heikentävän OATP1B1:n kuljetustehoa. Tässä väitöskirjatyössä selvitettiin SLCO1B1-geenin perinnöllistä muuntelua suomalaisilla ja eri väestöissä maailmanlaajuisesti. Lisäksi selvitettiin SLCO1B1:n muunnosten vaikutusta eri statiinien pitoisuuksiin (farmakokinetiikka) ja vaikutuksiin (farmakodynamiikka) sekä kolesteroliaineenvaihduntaan. Näihin tutkimuksiin valittiin SLCO1B1-genotyypin perusteella terveitä vapaaehtoisia koehenkilöitä, joille annettiin eri päivinä kerta-annos kutakin tutkittavaa statiinia: fluvastatiinia, pravastatiinia, simvastatiinia, rosuvastatiinia ja atorvastatiinia. Verinäytteistä määritettiin plasman statiinien ja niiden aineenvaihduntatuotteiden sekä kolesterolin ja sen muodostumista ja imeytymistä kuvaavien merkkiaineiden pitoisuuksia. Toiminnallisesti merkittävien SLCO1B1-geenimuunnosten esiintyvyydessä todettiin suuria eroja eri väestöjen välillä. Suomalaisilla SLCO1B1 c.521TC-genotyypin (geenimuunnos toisessa vastinkromosomissa) esiintyvyys oli noin 32 % ja SLCO1B1 c.521CC-genotyypin (geenimuunnos molemmissa vastinkromosomeissa) esiintyvyys noin 4 %. Globaalisti geenimuunnosten esiintyvyys korreloi maapallon leveyspiirien kanssa siten, että matalaan transportteriaktiivisuuteen johtavat muunnokset olivat yleisimpiä pohjoisessa ja korkeaan aktiivisuuteen johtavat päiväntasaajan lähellä asuvilla väestöillä. SLCO1B1-genotyypillä oli merkittävä vaikutus statiinien plasmapitoisuksiin lukuun ottamatta fluvastatiinia. Simvastatiinihapon plasmapitoisuudet olivat keskimäärin 220 %, atorvastatiinin 140 %, pravastatiinin 90 % ja rosuvastatiinin 70 % suuremmat c.521CC-genotyypin omaavilla koehenkilöillä verrattuna normaalin c.521TT-genotyypin omaaviin. Genotyypillä ei ollut merkittävää vaikutusta minkään statiinin tehoon tässä kerta-annostutkimuksessa, mutta geenimuunnoksen kantajilla perustason kolesterolisynteesinopeus oli suurempi. Tulokset osoittavat, että SLCO1B1 c.521T>C geenimuunnos on varsin yleinen suomalaisilla ja muilla ei-afrikkalaisilla väestöillä. Tämä geenimuunnos voi altistaa erityisesti simvastatiinin, mutta myös atorvastatiinin, pravastatiinin ja rosuvastatiinin, aiheuttamille lihashaitoille suurentamalla niiden plasmapitoisuuksia. SLCO1B1:n geenimuunnoksen testaamista voidaan tulevaisuudessa käyttää apuna valittaessa sopivaa statiinilääkitystä ja -annosta potilaalle, ja näin parantaa sekä statiinihoidon turvallisuutta että tehoa.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

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.

Ecology, population genetics and conservation of the African violet

The ongoing rapid fragmentation of tropical forests is a major threat to global biodiversity. This is because many of the tropical forests are so-called biodiversity 'hotspots', areas that host exceptional species richness and concentrations of endemic species. Forest fragmentation has negative ecological and genetic consequences for plant survival. Proposed reasons for plant species' loss in forest fragments are, e.g., abiotic edge effects, altered species interactions, increased genetic drift, and inbreeding depression. To be able to conserve plants in forest fragments, the ecological and genetic processes that threaten the species have to be understood. That is possible only after obtaining adequate information on their biology, including taxonomy, life history, reproduction, and spatial and genetic structure of the populations. In this research, I focused on the African violet (genus Saintpaulia), a little-studied conservation flagship from the Eastern Arc Mountains and Coastal Forests hotspot of Tanzania and Kenya. The main objective of the research was to increase understanding of the life history, ecology and population genetics of Saintpaulia that is needed for the design of appropriate conservation measures. A further aim was to provide population-level insights into the difficult taxonomy of Saintpaulia. Ecological field work was conducted in a relatively little fragmented protected forest in the Amani Nature Reserve in the East Usambara Mountains, in northeastern Tanzania, complemented by population genetic laboratory work and ecological experiments in Helsinki, Finland. All components of the research were conducted with Saintpaulia ionantha ssp. grotei, which forms a taxonomically controversial population complex in the study area. My results suggest that Saintpaulia has good reproductive performance in forests with low disturbance levels in the East Usambara Mountains. Another important finding was that seed production depends on sufficient pollinator service. The availability of pollinators should thus be considered in the in situ management of threatened populations. Dynamic population stage structures were observed suggesting that the studied populations are demographically viable. High mortality of seedlings and juveniles was observed during the dry season but this was compensated by ample recruitment of new seedlings after the rainy season. Reduced tree canopy closure and substrate quality are likely to exacerbate seedling and juvenile mortality, and, therefore, forest fragmentation and disturbance are serious threats to the regeneration of Saintpaulia. Restoration of sufficient shade to enhance seedling establishment is an important conservation measure in populations located in disturbed habitats. Long-term demographic monitoring, which enables the forecasting of a population s future, is also recommended in disturbed habitats. High genetic diversities were observed in the populations, which suggest that they possess the variation that is needed for evolutionary responses in a changing environment. Thus, genetic management of the studied populations does not seem necessary as long as the habitats remain favourable for Saintpaulia. The observed high levels of inbreeding in some of the populations, and the reduced fitness of the inbred progeny compared to the outbred progeny, as revealed by the hand-pollination experiment, indicate that inbreeding and inbreeding depression are potential mechanisms contributing to the extinction of Saintpaulia populations. The relatively weak genetic divergence of the three different morphotypes of Saintpaulia ionantha ssp. grotei lend support to the hypothesis that the populations in the Usambara/lowlands region represent a segregating metapopulation (or metapopulations), where subpopulations are adapting to their particular environments. The partial genetic and phenological integrity, and the distinct trailing habit of the morphotype 'grotei' would, however, justify its placement in a taxonomic rank of its own, perhaps in a subspecific rank.

Population genetics of the invasive plant species Impatiens glandulifera in Southern Finland

Biological invasions affect biodiversity worldwide, and, consequently, the invaded ecosystems may suffer from significant losses in economic and cultural values. Impatiens glandulifera Royle (Balsaminaceae) is an invasive annual herb, native to the western Himalayas and introduced into Europe in the 19th century as a garden ornamental plant. The massive invasion of I. glandulifera is due to its high reproductive output, rapid growth and its ability to outcompete native species. In Finland, the first observations regarding the presence of I. glandulifera date from the year 1947, and today it is considered a serious problem in riparian habitats. The aim of this master’s thesis research is to reveal the population genetic structure of I. glandulifera in Finland and to find out whether there have been one or multiple invasions in Finland. The study focuses on investigating the origin of I. glandulifera in Southern Finland, by comparing plant samples from the Helsinki region with those from its native region and other regions of invasion. Samples from four populations in Helsinki and from the United Kingdom, Canada, India and Pakistan were collected and genotyped using 11 microsatellite markers. The genetic analyses were evaluated using the programs Arlequin and Structure. The results of the genetic analyses suggested that I. glandulifera has been introduced to Finland more than once. Multiple introductions are supported by the higher level of genetic diversity detected within and among Finnish populations than would be expected for a single introduction. Results of the Bayesian Structure analysis divided the four Finnish populations into four clusters. This geographical structure was further supported by pairwise Fst values among populations. The causes and potential consequences of such multiple introductions of I. glandulifera in Finland and further perspectives are discussed.

Sindbis virus and Pogosta disease in Finland

Sindbis virus (SINV) (genus Alphavirus, family Togaviridae) is an enveloped virus with a genome of single-stranded, positive-polarity RNA of 11.7 kilobases. SINV is widespread in Eurasia, Africa, and Australia, but clinical infection only occurs in a few geographically restricted areas, mainly in Northern Europe. In Europe, antibodies to SINV were detected from patients with fever, rash, and arthritis for the first time in the early 1980s in Finland. It became evident that the causative agent of this syndrome, named Pogosta disease, was closely related to SINV. The disease is also found in Sweden (Ockelbo disease) and in Russia (Karelian fever). Since 1974, for unknown reason, the disease has occurred as large outbreaks every seven years in Finland. This study is to a large degree based on the material collected during the 2002 Pogosta disease outbreak in Finland. We first developed SINV IgM and IgG enzyme immunoassays (EIA), based on highly purified SINV, to be used in serodiagnostics. The EIAs correlated well with the hemagglutination inhibition (HI) test, and all individuals showed neutralizing antibodies. The sensitivities of the IgM and IgG EIAs were 97.6% and 100%, and specificities 95.2% and 97.6%, respectively. E1 and E2 envelope glycoproteins of SINV were shown to be recognized by IgM and IgG in the immunoblot early in infection. We isolated SINV from five patients with acute Pogosta disease; one virus strain was recovered from whole blood, and four other strains from skin lesions. The etiology of Pogosta disease was confirmed by these first Finnish SINV strains, also representing the first human SINV isolates from Europe. Phylogenetic analysis indicated that the Finnish SINV strains clustered with the strains previously isolated from mosquitoes in Sweden and Russia, and seemed to have a common ancestor with South-African strains. Northern European SINV strains could be maintained locally in disease-endemic regions, but the phylogenetic analysis also suggests that redistribution of SINV tends to occur in a longitudinal direction, possibly with migratory birds. We searched for SINV antibodies in resident grouse (N=621), whose population crashes have previously coincided with human SINV outbreaks, and in migratory birds (N=836). SINV HI antibodies were found for the first time in birds during their spring migration to Northern Europe, from three individuals: red-backed shrike, robin, and song thrush. Of the grouse, 27.4% were seropositive in 2003, one year after a human outbreak, but only 1.4% of the grouse were seropositive in 2004. Thus, grouse might contribute to the human epidemiology of SINV. A total of 86 patients with verified SINV infection were recruited to the study in 2002. SINV RNA detection or virus isolation from blood and/or skin lesions was successful in eight patients. IgM antibodies became detectable within the first eight days of illness, and IgG within 11 days. The acute phase of Pogosta disease was characterized by arthritis, itching rash, fatigue, mild fever, headache, and muscle pain. Half of the patients reported in self-administered questionnaires joint symptoms to last > 12 months. Physical examination in 49 of these patients three years after infection revealed persistent joint manifestations. Arthritis (swelling and tenderness in physical examination) was diagnosed in 4.1% (2/49) of the patients. Tenderness in palpation or in movement of a joint was found in 14.3% of the patients in the rheumatologic examination, and additional 10.2% complained persisting arthralgia at the interview. Thus, 24.5% of the patients had joint manifestations attributable to the infection three years earlier. A positive IgM antibody response persisted in 3/49 of the patients; both two patients with arthritis were in this group. Persistent symptoms of SINV infection might have considerable public health implications in areas with high seroprevalence. The age-standardized seroprevalence of SINV (1999-2003, N=2529) in the human population in Finland was 5.2%. The seroprevalence was high in North Karelia, Kainuu, and Central Ostrobothnia. The incidence was highest in North Karelia. Seroprevalence in men (6.0%) was significantly higher than in women (4.1%), however, the average annualized incidence in the non-epidemic years was higher in women than in men, possibly indicating that infected men are more frequently asymptomatic. The seroprevalence increased with age, reaching 15.4% in persons aged 60-69 years. The incidence was highest in persons aged 50-59 years.

Dengue virus infection : Diagnostics and molecular epidemiology

Dengue is a mosquito-borne viral disease caused by the four dengue virus serotypes (DENV-1-4) and is currently considered as the most important arthropod-borne viral disease in the world. Nearly half of the human population lives in risk areas, and 50-100 million infections occur yearly according to World Health Organization. The disease can vary from a mild febrile disease to severe haemorrhagic fever and shock. A secondary infection with heterologous serotype increases the risk for severe disease outcome. During the last three decades the impact of dengue has dramatically increased in the endemic areas including the tropics and subtropics of the world. The current situation with massive epidemics of severe disease forms has been associated with socio-ecological changes that have increased the transmission and enabled the co-circulation of different serotypes. Consequently, an increase of dengue has also been observed in travelers visiting these areas. Currently approximately 30 cases are diagnosed yearly in Finnish travelers. In travelers dengue is rarely a life-threatening disease, however in the current study, a fatality was documented in a young Finnish patient who experienced a prolonged primary dengue infection. To improve particularly early laboratory diagnostics, a novel real-time RT-PCR method was developed for the detection of DENV-1-4 RNA based on TaqMan chemistry. The method was shown to be sensitive and specific for detecting DENV RNA and suitable for diagnostic use. The newly developed real-time RT-PCR was compared to other available early diagnostic methods including IgM and NS1 antigen detection using a panel of selected patient samples. The results suggest that the best diagnostic rates are achieved by a combination of IgM with RNA or NS1 detection. The dengue virus strains studied here included the first DENV strains isolated from serum samples of Finnish travelers collected in 2000-2005. The results of sequence analysis demonstrated that the 11 isolates included all four DENV serotypes and presented a global sample of DENV strains from different geographical areas including Asia, Africa and South America. In the present study sequence analysis was also carried out for a collection of 23 novel DENV-2 isolates from Venezuelan patients collected in 1999-2005. The Venezuelan DENV-2 exclusively represented the American-Asian genotype, suggesting that no foreign DENV-2 lineages have recently been introduced to the country. The results also suggest that the DENV-2 viruses detected earlier from Venezuela have been maintained in the area where they have evolved into several lineages. This is in contrast to the pattern observed in some other dengue endemic areas, where introductions of novel virus types and lineages are frequently detected.

Contested Lands: Land Disputes in Semi-Arid Parts of Northern Tanzania : Case Studies of the Loliondo and Sale Divisions

This study Contested Lands: Land disputes in semi-arid parts of northern Tanzania. Case Studies of the Loliondo and Sale Division in the Ngorongoro District concentrates on describing the specific land disputes which took place in the 1990s in the Loliondo and Sale Divisions of the Ngorongoro District in northern Tanzania. The study shows the territorial and historical transformation of territories and property and their relation to the land disputes of the 1990s'. It was assumed that land disputes have been firstly linked to changing spatiality due to the zoning policies of the State territoriality and, secondly, they can be related to the State control of property where the ownership of land property has been redefined through statutory laws. In the analysis of the land disputes issues such as use of territoriality, boundary construction and property claims, in geographical space, are highlighted. Generally, from the 1980s onwards, increases in human population within both Divisions have put pressure on land/resources. This has led to the increased control of land/resource, to the construction of boundaries and finally to formalized land rights on village lands of the Loliondo Division. The land disputes have thus been linked to the use of legal power and to the re-creation of the boundary (informal or formal) either by the Maasai or the Sonjo on the Loliondo and Sale village lands. In Loliondo Division land disputes have been resource-based and related to multiple allocations of land or game resource concessions. Land disputes became clearly political and legal struggles with an ecological reference.Land disputes were stimulated when the common land/resource rights on village lands of the Maasai pastoralists became regulated and insecure. The analysis of past land disputes showed that space-place tensions on village lands can be presented as a platform on which spatial and property issues with complex power relations have been debated. The reduction of future land disputes will succeed only when/if local property rights to land and resources are acknowledged, especially in rural lands of the Tanzanian State.

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

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

Enhanced processing of SPOT multispectral satellite imagery for environmental monitoring and modelling

The Taita Hills in southeastern Kenya form the northernmost part of Africa’s Eastern Arc Mountains, which have been identified by Conservation International as one of the top ten biodiversity hotspots on Earth. As with many areas of the developing world, over recent decades the Taita Hills have experienced significant population growth leading to associated major changes in land use and land cover (LULC), as well as escalating land degradation, particularly soil erosion. Multi-temporal medium resolution multispectral optical satellite data, such as imagery from the SPOT HRV, HRVIR, and HRG sensors, provides a valuable source of information for environmental monitoring and modelling at a landscape level at local and regional scales. However, utilization of multi-temporal SPOT data in quantitative remote sensing studies requires the removal of atmospheric effects and the derivation of surface reflectance factor. Furthermore, for areas of rugged terrain, such as the Taita Hills, topographic correction is necessary to derive comparable reflectance throughout a SPOT scene. Reliable monitoring of LULC change over time and modelling of land degradation and human population distribution and abundance are of crucial importance to sustainable development, natural resource management, biodiversity conservation, and understanding and mitigating climate change and its impacts. The main purpose of this thesis was to develop and validate enhanced processing of SPOT satellite imagery for use in environmental monitoring and modelling at a landscape level, in regions of the developing world with limited ancillary data availability. The Taita Hills formed the application study site, whilst the Helsinki metropolitan region was used as a control site for validation and assessment of the applied atmospheric correction techniques, where multiangular reflectance field measurements were taken and where horizontal visibility meteorological data concurrent with image acquisition were available. The proposed historical empirical line method (HELM) for absolute atmospheric correction was found to be the only applied technique that could derive surface reflectance factor within an RMSE of < 0.02 ps in the SPOT visible and near-infrared bands; an accuracy level identified as a benchmark for successful atmospheric correction. A multi-scale segmentation/object relationship modelling (MSS/ORM) approach was applied to map LULC in the Taita Hills from the multi-temporal SPOT imagery. This object-based procedure was shown to derive significant improvements over a uni-scale maximum-likelihood technique. The derived LULC data was used in combination with low cost GIS geospatial layers describing elevation, rainfall and soil type, to model degradation in the Taita Hills in the form of potential soil loss, utilizing the simple universal soil loss equation (USLE). Furthermore, human population distribution and abundance were modelled with satisfactory results using only SPOT and GIS derived data and non-Gaussian predictive modelling techniques. The SPOT derived LULC data was found to be unnecessary as a predictor because the first and second order image texture measurements had greater power to explain variation in dwelling unit occurrence and abundance. The ability of the procedures to be implemented locally in the developing world using low-cost or freely available data and software was considered. The techniques discussed in this thesis are considered equally applicable to other medium- and high-resolution optical satellite imagery, as well the utilized SPOT data.

Ecological consequences of genetic modifications - an invasion analysis approach

Biological invasions are considered as one of the greatest threats to biodiversity, as they may lead to disruption and homogenization of natural communities, and in the worst case, to native species extinctions. The introduction of gene modified organisms (GMOs) to agricultural, fisheries and forestry practices brings them into contact with natural populations. GMOs may appear as new invasive species if they are able to (1) invade into natural habitats or (2) hybridize with their wild relatives. The benefits of GMOs, such as increased yield or decreased use of insecticides or herbicides in cultivation, may thus be reduced due the potential risks they may cause. A careful ecological risk analysis therefore has to precede any responsible GMO introduction. In this thesis I study ecological invasion in relation to GMOs, and what kind of consequences invasion may have in natural populations. A set of theoretical models that combine life-history evolution, population dynamics, and population genetics were developed for the hazard identification part of ecological risks assessment of GMOs. In addition, the potential benefits of GMOs in management of an invasive pest were analyzed. In the first study I showed that a population that is fluctuating due to scramble-type density dependence (due to, e.g., nutrient competition in plants) may be invaded by a population that is relatively more limited by a resource (e.g., light in plants) that is a cause of contest-type density dependence. This result emphasises the higher risk of invasion in unstable environments. The next two studies focused on escape of a growth hormone (GH) transgenic fish into a natural population. The results showed that previous models may have given too pessimistic a view of the so called Trojan gene -effect, where the invading genotype is harmful for the population as a whole. The previously suggested population extinctions did not occur in my studies, since the changes in mating preferences caused by the GH-fish were be ameliorated by decreased level of competition. The GH-invaders may also have to exceed a threshold density before invasion can be successful. I also showed that the prevalence of mature parr (aka. sneaker) strategy among GH-fish may have clear effect on invasion outcome. The fourth study assessed the risks and developed methods against the invasion of the Colorado Potato Beetle (CPB, Leptinotarsa decemlineata). I showed that the eradication of CPB is most important for the prevention of their establishment, but the cultivation of transgenic Bt-potato could also be effective. In general, my results emphasise that invasion of transgenic species or genotypes to be possible under certain realistic conditions and resulting in competitive exclusion, population decline through outbreeding depression and genotypic displacement of native species. Ecological risk assessment should regard the decline and displacement of the wild genotype by an introduced one as a consequence that is as serious as the population extinction. It will also be crucial to take into account different kinds of behavioural differences among species when assessing the possible hazards that GMOs may cause if escaped. The benefits found of GMO crops effectiveness in pest management may also be too optimistic since CPB may evolve resistance to Bt-toxin. The models in this thesis could be further applied in case specific risk assessment of GMOs by supplementing them with detailed data of the species biology, the effect of the transgene introduced to the species, and also the characteristics of the populations or the environments in the risk of being invaded.

Applications of gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other teleost fish species

Evolutionary genetics incorporates traditional population genetics and studies of the origins of genetic variation by mutation and recombination, and the molecular evolution of genomes. Among the primary forces that have potential to affect the genetic variation within and among populations, including those that may lead to adaptation and speciation, are genetic drift, gene flow, mutations and natural selection. The main challenges in knowing the genetic basis of evolutionary changes is to distinguish the adaptive selection forces that cause existent DNA sequence variants and also to identify the nucleotide differences responsible for the observed phenotypic variation. To understand the effects of various forces, interpretation of gene sequence variation has been the principal basis of many evolutionary genetic studies. The main aim of this thesis was to assess different forms of teleost gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other species. Firstly, the level of Darwinian adaptive evolution affected coding regions of the growth hormone (GH) gene during the teleost evolution was investigated based on the sequence data existing in public databases. Secondly, a target gene approach was used to identify within population variation in the growth hormone 1 (GH1) gene in salmon. Then, a new strategy for single nucleotide polymorphisms (SNPs) discovery in salmonid fishes was introduced, and, finally, the usefulness of a limited number of SNP markers as molecular tools in several applications of population genetics in Atlantic salmon was assessed. This thesis showed that the gene sequences in databases can be utilized to perform comparative studies of molecular evolution, and some putative evidence of the existence of Darwinian selection during the teleost GH evolution was presented. In addition, existent sequence data was exploited to investigate GH1 gene variation within Atlantic salmon populations throughout its range. Purifying selection is suggested to be the predominant evolutionary force controlling the genetic variation of this gene in salmon, and some support for gene flow between continents was also observed. The novel approach to SNP discovery in species with duplicated genome fragments introduced here proved to be an effective method, and this may have several applications in evolutionary genetics with different species - e.g. when developing gene-targeted markers to investigate quantitative genetic variation. The thesis also demonstrated that only a few SNPs performed highly similar signals in some of the population genetic analyses when compared with the microsatellite markers. This may have useful applications when estimating genetic diversity in genes having a potential role in ecological and conservation issues, or when using hard biological samples in genetic studies as SNPs can be applied with relatively highly degraded DNA.

The influence of eutrophication on sexual selection in sticklebacks

Jordens ekologiska system undergår för tillfället stora förändringar pga. människans aktiviteter. Ett växande antal studier visar att dessa förändringar påverkar naturliga och sexuella urvalet och därmed evolutiva processer. Målet med detta arbete var att undersöka effekterna av omgivningsförändringar på sexuella urvalet genom att använda den ökade övergödningen inom storpiggen Gasterosteus aculeatus lekområden som modell system. Sexuella urvalet är en viktig evolutiv kraft med följder på populations- och artnivå (Kapitel 1). Avhandlingens olika delar fokuserar på övergödningens effekter på upptäckandet av partners, användningen av visuella- och doftsignaler i partnersval, och fördelningen av parningsframgången mellan bobyggande hanar. I Kapitel II och III simuleras hur grumlighet orsakad av fytoplankton påverkar hastigheten med vilken potentiella partners påträffas, genom effekter på synligheten. Resultaten visar att normala algblomningar i Östersjön har en måttlig effekt på finnandet av potentiella partners. Detta tyder på att algblomningarna troligen inte kommer att minska på selektiva parningen pga. ökade sökkostnader. I Kapitel IV visas att storspiggen ändrar relativa användningen av olika signaler när vattnets grumlighet ökar; visuella signaler minskar i betydelse medan doftsignaler ökar i betydelse. Samtidigt underlättas användandet av doftsignaler av ändringar i vattnets kemiska sammansättning då fotosyntesen intensifieras (Kapitel V). Lek i övergödda vatten kan ändå vara kostsamt både på individ- och populationsnivån, då parasiterade hanar, som troligen är dåligt genetiskt anpassade till sin miljö, lyckas få mer ägg i sina bon än friskare hanar som troligen är av högre genetisk kvalitet (Kapitel VI). Övergödningen påverkar således partnersval och konkurrensen om partners genom att påverka upptäckandet av potentiella partners, evalueringen av partners och fördelningen av partners inom lekområdena. De följder detta kan ha för evolutionen av sexuellt selekterad egenskaper och för populationers dynamik och livskraft är dock oklara. Avhandlingen visar på svårigheten att förutse följderna av omgivningsförändringar för sexuella urvalet och effekterna på individ och populationsnivå.

Interplay of Virulence Factors in Complement Resistance of Streptococcus pneumoniae

Streptococcus pneumoniae (pneumococcus) is a normal inhabitant of the human nasopharynx. Symptoms occur in only a small proportion of those who become carriers, but the ubiquity of the organism in the human population results in a large burden of disease. S. pneumoniae is the leading bacterial cause of pneumonia, sepsis, and meningitis worldwide, causing the death of a million children each year. Middle-ear infection is the most common clinical manifestation of mucosal pneumococcal infections. In invasive disease, S. pneumoniae gains access to the bloodstream and spreads to normally sterile parts of the body. The progression from asymptomatic colonization to disease depends on factors characteristic of specific pneumococcal strains as well as the status of host defenses. The polysaccharide capsule surrounding the bacterium is considered to be the most important factor affecting the virulence of pneumococci. It protects pneumococci from phagocytosis and also may determine its affinity to the respiratory epithelium. S. pneumoniae as a species comprises more than 90 different capsular serotypes, but not all of them are equally prevalent in human diseases. Invasive serotypes are rarely isolated from healthy carriers, but relatively often cause invasive disease. Serotypes that are carried asymptomatically for a long time behave like opportunistic pathogens, causing disease in patients who have impaired immune defenses. The complement system is a collection of blood and cell surface proteins that act as a major primary defense against invading microbes. Phagocytic cells with receptors for complement proteins can engulf and destroy pneumococcal cells opsonized with these proteins. S. pneumoniae has evolved a number of ways to subvert mechanisms of innate immunity, and this is likely to contribute to its pathogenicity. The capsular serotype, proteins essential for virulence, as well the genotype, may all influence the ability of pneumococcus to resist complement and its potential to cause disease. Immunization with conjugate vaccines produces opsonic antibodies, which enhance complement deposition and clearance of the bacteria. The pneumococcal vaccine included in the Finnish national immunization program in 2010 contains the most common serotypes causing invasive disease. Clinical data suggest that protection from middle-ear infection and possibly also from invasive disease depends largely on the capsular serotype, for reasons hitherto unknown. The general aim of this thesis is to assess the relative roles of the pneumococcal capsule and virulence proteins in complement evasion and subsequent opsonophagocytic killing. The main question is whether differences between serotypes to resist complement explain the different abilities of serotypes to cause disease. The importance of particular virulence factors to the complement resistance of a strain may vary depending on its genotype. Prior studies have evaluated the effect of the capsule and virulence proteins on complement resistance of S. pneumoniae by comparing only a few strains. In this thesis, the role of pneumococcal virulence factors in the complement resistance of the bacterium was studied in several genotypically different strains. The ability of pneumococci to inhibit deposition of the complement protein C3 on the bacterial surface was found to depend on the capsular serotype as well as on other features of the bacteria. The results suggest that pneumococcal histidine triad (Pht) proteins may play a role in complement inhibition, but their contribution depends on the bacterial genotype. The capsular serotype was found to influence complement resistance more than the bacterial genotype. A higher concentration of anticapsular antibodies was required for the opsonophagocytic killing of serotypes resistant to C3 deposition. The invasive serotypes were more resistant to C3 deposition than the opportunistic serotypes, suggesting that the former are better adapted to resist immune mechanisms controlling the development of invasive disease. The different susceptibilities of serotypes to complement deposition, opsonophagocytosis, and resultant antibody-mediated protection should be taken into account when guidelines for serological correlates for vaccine efficacy evaluations are made. The results of this thesis suggest that antibodies in higher quantity or quality are needed for efficient protection against the invasive serotypes.