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.

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.

Sedimentary zooplankton remains as indicators of lake ecological quality and trophic structure

To protect and restore lake ecosystems under threats posed by the increasing human population, information on their ecological quality is needed. Lake sediments provide a data rich archive that allows identification of various biological components present prior to anthropogenic alterations as well as a constant record of changes. By providing a longer dimension of time than any ongoing monitoring programme, palaeolimnological methods can help in understanding natural variability and long-term ecological changes in lakes. As zooplankton have a central role in the lake food web, their remains can potentially provide versatile information on past trophic structure. However, various taphonomic processes operating in the lakes still raise questions concerning how subfossil assemblages reflect living communities. This thesis work aimed at improving the use of sedimentary zooplankton remains in the reconstruction of past zooplankton communities and the trophic structure in lakes. To quantify interspecific differences in the accumulation of remains, the subfossils of nine pelagic zooplankton taxa in annually laminated sediments were compared with monitoring results for live zooplankton in Lake Vesijärvi. This lake has a known history of eutrophication and recovery, which resulted from reduced external loading and effective fishing of plankti-benthivorous fish. The response of zooplankton assemblages to these known changes was resolved using annually laminated sediments. The generality of the responses observed in Lake Vesijärvi were further tested with a set of 31 lakes in Southern Finland, relating subfossils in surface sediments to contemporary water quality and fish density, as well as to lake morphometry. The results demonstrated differential preservation and retention of cladoceran species in the sediment. Daphnia, Diaphanosoma and Ceriodaphnia were clearly underrepresented in the sediment samples in comparison to well-preserved Bosmina species, Chydorus, Limnosida and Leptodora. For well-preserved species, the annual net accumulation rate was similar to or above the expected values, reflecting effective sediment focusing and accumulation in the deepest part of the lake. The decreased fish density and improved water quality led to subtle changes in zooplankton community composition. The abundance of Diaphanosoma and Limnosida increased after the reduction in fish density, while Ceriodaphnia and rotifers decreased. The most sensitive indicator of fish density was the mean size of Daphnia ephippia and Bosmina (E.) crassicornis ephippia and carapaces. The concentration of plant-associated species increased, reflecting expanding littoral vegetation along with increasing transparency. Several of the patterns observed in Lake Vesijärvi could also be found within the set of 31 lakes. According to this thesis work, the most useful cladoceran-based indices for nutrient status and planktivorous fish density in Finnish lakes were the relative abundances of certain pelagic taxa, and the mean size of Bosmina spp. carapaces, especially those of Bosmina (E.) cf. coregoni. The abundance of plant-associated species reflected the potential area for aquatic plants. Lake morphometry and sediment organic content, however, explained a relatively high proportion of the variance in the species data, and more studies are needed to quantify lake-specific differences in the accumulation and preservation of remains. Commonly occurring multicollinearity between environmental variables obstructs the cladoceran-based reconstruction of single environmental variables. As taphonomic factors and several direct and indirect structuring forces in lake ecosystems simultaneously affect zooplankton, the subfossil assemblages should be studied in a holistic way before making final conclusions about the trophic structure and the change in lake ecological quality.

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.

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.

Pharmacogenetic Variation at CYP2D6, CYP2C9, and CYP2C19: Population Genetic and Forensic Aspects

Pharmacogenetics deals with genetically determined variation in drug response. In this context, three phase I drug-metabolizing enzymes, CYP2D6, CYP2C9, and CYP2C19, have a central role, affecting the metabolism of about 20-30% of clinically used drugs. Since genes coding for these enzymes in human populations exhibit high genetic polymorphism, they are of major pharmacogenetic importance. The aims of this study were to develop new genotyping methods for CYP2D6, CYP2C9, and CYP2C19 that would cover the most important genetic variants altering the enzyme activity, and, for the first time, to describe the distribution of genetic variation at these loci on global and microgeographic scales. In addition, pharmacogenetics was applied to a postmortem forensic setting to elucidate the role of genetic variation in drug intoxications, focusing mainly on cases related to tricyclic antidepressants, which are commonly involved in fatal drug poisonings in Finland. Genetic variability data were obtained by genotyping new population samples by the methods developed based on PCR and multiplex single-nucleotide primer extension reaction, as well as by collecting data from the literature. Data consisted of 138, 129, and 146 population samples for CYP2D6, CYP2C9, and CYP2C19, respectively. In addition, over 200 postmortem forensic cases were examined with respect to drug and metabolite concentrations and genotypic variation at CYP2D6 and CYP2C19. The distribution of genetic variation within and among human populations was analyzed by descriptive statistics and variance analysis and by correlating the genetic and geographic distances using Mantel tests and spatial autocorrelation. The correlation between phenotypic and genotypic variation in drug metabolism observed in postmortem cases was also analyzed statistically. The genotyping methods developed proved to be informative, technically feasible, and cost-effective. Detailed molecular analysis of CYP2D6 genetic variation in a global survey of human populations revealed that the pattern of variation was similar to those of neutral genomic markers. Most of the CYP2D6 diversity was observed within populations, and the spatial pattern of variation was best described as clinal. On the other hand, genetic variants of CYP2D6, CYP2C9, and CYP2C19 associated with altered enzymatic activity could reach extremely high frequencies in certain geographic regions. Pharmacogenetic variation may also be significantly affected by population-specific demographic histories, as seen within the Finnish population. When pharmacogenetics was applied to a postmortem forensic setting, a correlation between amitriptyline metabolic ratios and genetic variation at CYP2D6 and CYP2C19 was observed in the sample material, even in the presence of confounding factors typical for these cases. In addition, a case of doxepin-related fatal poisoning was shown to be associated with a genetic defect at CYP2D6. Each of the genes studied showed a distinct variation pattern in human populations and high frequencies of altered activity variants, which may reflect the neutral evolution and/or selective pressures caused by dietary or environmental exposure. The results are relevant also from the clinical point of view since the genetic variation at CYP2D6, CYP2C9, and CYP2C19 already has a range of clinical applications, e.g. in cancer treatment and oral anticoagulation therapy. This study revealed that pharmacogenetics may also contribute valuable information to the medicolegal investigation of sudden, unexpected deaths.

Statistical analysis of the associations of hereditary factors with the risk of Type 1 diabetes and Diabetic Nephropathy based on the familial data collected through ascertaiment from population-based registers

In genetic epidemiology, population-based disease registries are commonly used to collect genotype or other risk factor information concerning affected subjects and their relatives. This work presents two new approaches for the statistical inference of ascertained data: a conditional and full likelihood approaches for the disease with variable age at onset phenotype using familial data obtained from population-based registry of incident cases. The aim is to obtain statistically reliable estimates of the general population parameters. The statistical analysis of familial data with variable age at onset becomes more complicated when some of the study subjects are non-susceptible, that is to say these subjects never get the disease. A statistical model for a variable age at onset with long-term survivors is proposed for studies of familial aggregation, using latent variable approach, as well as for prospective studies of genetic association studies with candidate genes. In addition, we explore the possibility of a genetic explanation of the observed increase in the incidence of Type 1 diabetes (T1D) in Finland in recent decades and the hypothesis of non-Mendelian transmission of T1D associated genes. Both classical and Bayesian statistical inference were used in the modelling and estimation. Despite the fact that this work contains five studies with different statistical models, they all concern data obtained from nationwide registries of T1D and genetics of T1D. In the analyses of T1D data, non-Mendelian transmission of T1D susceptibility alleles was not observed. In addition, non-Mendelian transmission of T1D susceptibility genes did not make a plausible explanation for the increase in T1D incidence in Finland. Instead, the Human Leucocyte Antigen associations with T1D were confirmed in the population-based analysis, which combines T1D registry information, reference sample of healthy subjects and birth cohort information of the Finnish population. Finally, a substantial familial variation in the susceptibility of T1D nephropathy was observed. The presented studies show the benefits of sophisticated statistical modelling to explore risk factors for complex diseases.

Development of regulatory T cells in the human thymus

The role of the immune system is to protect an organism against pathogens while maintaining tolerance against self. T cells are an essential component of the immune system and they develop in the thymus. The AIRE (autoimmune regulator) gene product plays an important role in T cell development, as it promotes expression of peripheral tissue antigens in the thymus. Developing T cells, thymocytes, which recognize self-antigens with high affinity are deleted. However, this deletion process is not perfect and not all autoreactive T cells are destroyed. When the distinction between self and non-self fails, tolerance breaks and the immune system attacks the host s own tissues. This results in autoimmunity. Regulatory T cells contribute to the maintenance of self-tolerance. They can actively suppress the function of autoreactive cells. Several populations of cells with regulatory properties have been described, but the best characterized population is the natural regulatory T cells (Treg cells), which develop in the thymus and express the transcription factor FOXP3. The thymic development of Treg cells in humans is the subject of this thesis. Thymocytes at different developmental stages were analyzed using flow cytometry. The CD4-CD8- double-negative (DN) thymocytes are the earliest T cell precursors in the T cell lineage. My results show that the Treg cell marker FOXP3 is up-regulated already in a subset of these DN thymocytes. FOXP3+ cells were also found among the more mature CD4+CD8+ double-positive (DP) cells and among the CD4+ and CD8+ single-positive (SP) thymocytes. The different developmental stages of the FOXP3+ thymocytes were isolated and their gene expression examined by quantitative PCR. T cell receptor (TCR) repertoire analysis was used to compare these different thymocyte populations. My data show that in humans commitment to the Treg cell lineage is an early event and suggest that the development of Treg cells follows a linear developmental pathway, FOXP3+ DN precursors evolving through the DP stage to become mature CD4+ Treg cells. Most T cells have only one kind of TCR on their cell surface, but a small fraction of cells expresses two different TCRs. My results show that the expression of two different TCRs is enriched among Treg cells. Furthermore, both receptors were capable of transmitting signals when bound by a ligand. By extrapolating flow cytometric data, it was estimated that the majority of peripheral blood Treg cells are indeed dual-specific. The high frequency of dual-specific cells among human Treg cells suggests that dual-specificity has a role in directing these cells to the Treg cell lineage. It is known that both genetic predisposition and environmental factors influence the development of autoimmunity. It is also known that the dysfunction or absence of Treg cells leads to the development of autoimmune manifestations. APECED (autoimmune polyendocrinopathy-candidiasis-ectodermal dystrophy) is a rare monogenic autoimmune disease, caused by mutations in the AIRE gene. In the absence of AIRE gene product, deletion of self-specific T cells is presumably disturbed and autoreactive T cells escape to the periphery. I examined whether Treg cells are also affected in APECED. I found that the frequency of FOXP3+ Treg cells and the level of FOXP3 expression were significantly lower in APECED patients than in controls. Additionally, when studied in cell cultures, the suppressive capacity of the patients' Treg cells was impaired. Additionally, repertoire analysis showed that the TCR repertoire of Treg cells was altered. These results suggest that AIRE contributes to the development of Treg cells in humans and the selection of Treg cells is impaired in APECED patients. In conclusion, my thesis elucidates the developmental pathway of Treg cells in humans. The differentiation of Tregs begins early during thymic development and both the cells dual-specificity and AIRE probably affect the final commitment of Treg cells.

Examination of Genetic Components Affecting Human Obesity-Related Quantitative Traits

Obesity increases the risk for several conditions, including type 2 diabetes mellitus, cardiovascular disease, hypertension, osteoarthirits and certain types of cancer. Twin- and family studies have shown that there is a major genetic component in the determination of body mass. In recent years several technological and scientific advance have been made in obesity research. For instance, novel replicated loci have been revealed by a number of genome wide association studies. This thesis aimed to investigate the association of genetic factors and obesity-related quantitative traits. The first study investigated the role of the lactase gene in anthropometric traits. We genetically defined lactose persistence by genotyping 31 720 individuals of European descent. We found that lactase persistence was significantly correlated with weight and body mass index but not with height. In the second study we performed the largest whole genome linkage scan for body mass index to date. The sample consisted of 4401 twin families and 10 535 individuals from six European countries. We found supporting evidence for two loci (3q29 and 7q36). We observed that the heritability estimate increased substantially when additional family members were removed from the analyses, which suggests reduced environmental variance in the twin sample. In the third study we assessed metabonomic, transcriptomic and genomic variation in a Finnish population cohort of 518 individuals. We formed gene expression networks to portray pathways and showed that a set of highly correlated genes of an inflammatory pathway associated with 80 serum metabolites (of 134 quantified measures). Strong association was found, for example, with several lipoprotein subclasses. We inferred causality by using genetic variation as anchors. The expression of the network genes was found to be dependent on the circulatory metabolite concentrations.

Search for genetic variants influencing human height

Human growth and attained height are determined by a combination of genetic and environmental effects and in modern Western societies > 80% of the observed variation in height is determined by genetic factors. Height is a fundamental human trait that is associated with many socioeconomic and psychosocial factors and health measures, however little is known of the identity of the specific genes that influence height variation in the general population. This thesis work aimed to identify the genetic variants that influence height in the general population by genome-wide linkage analysis utilizing large family samples. The study focused on analysis of three separate sets of families consisting of: 1) 1,417 individuals from 277 Finnish families (FinnHeight), 2) 8,450 individuals from 3,817 families from Australia and Europe (EUHeight) and 3) 9,306 individuals from 3,302 families from the United States (USHeight). The most significant finding in this study was found in the Finnish family sample where we a locus in the chromosomal region 1p21 was linked to adult height. Several regions showed evidence for linkage in the Australian, European and US families with 8q21 and 15q25 being the most significant. The region on 1p21 was followed up with further studies and we were able to show that the collagen 11-alpha-1 gene (COL11A1) residing at this location was associated with adult height. This association was also confirmed in an independent Finnish population cohort (Health 2000) consisting of 6,542 individuals. From this population sample, we estimated that homozygous males and females for this gene variant were 1.1 and 0.6 cm taller than the respective controls. In this thesis work we identified a gene variant in the COL11A1 gene that influences human height, although this variant alone explains only 0.1% of height variation in the Finnish population. We also demonstrated in this study that special stratification strategies such as performing sex-limited analyses, focusing on dizygous twin pairs, analyzing ethnic groups within a population separately and utilizing homogenous populations such as the Finns can improve the statistical power of finding QTL significantly. Also, we concluded from the results of this study that even though genetic effects explain a great proportion of height variance, it is likely that there are tens or even hundreds of genes with small individual effects underlying the genetic architecture of height.