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.

Use of ecological information in urban planning

Urbanization leads to irreversible land-use change, which has ecological consequences such as the loss and fragmentation of green areas, and structural and functional changes in terrestrial and aquatic ecosystems. These consequences diminish ecosystem services important for human populations living in urban areas. All this results in a conflict situation: how to simultaneously meet the needs of city growth and the principles of sustainable development, and especially conserve important green areas within and around built-up areas? Urban planners and decisionmakers have an important role in this, since they must use the ecological information mainly from species and biotope inventories and biodiversity impact assessments in determining the conservation values of green areas. The main aim of this thesis was to study the use of ecological information in the urban land-use planning and decisionmaking process in the Helsinki Metropolitan Area, Finland. At first, the literature on ecological-social systems linkages related to urban planning was reviewed. Based on the review, a theoretical and conceptual framework for the research on Finnish urban setting was adapted. Secondly, factors determining the importance and effectiveness of incorporation of ecological information into the urban planning process, and the challenges related to the use of ecological information were studied. Thirdly, the importance and use of Local Ecological Knowledge in urban planning were investigated. Then, factors determining the consideration of urban green areas and related ecological information in political land-use decisionmaking were studied. Finally, in a case study illustrating the above considerations, the importance of urban stream ecosystems in the land-use planning was investigated. This thesis demonstrated that although there are several challenges in using ecological information effectively, it is considered as an increasingly important part of the basic information used in urban planning and decisionmaking process. The basic determinants for this are the recent changes in environmental legislation, but also the increasing appreciation of green areas and their conservation values by all the stakeholders. In addition, Local Ecological Knowledge in its several forms can be a source of ecological information for planners if incorporated effectively into the process. This study also showed that rare or endangered species and biotopes, and related ecological information receive priority in the urban planning process and usually pass through the decisionmaking system. Furthermore, the stream Rekolanoja case indicates that planners and residents see the value of urban stream ecosystem as increasingly important for the local health and social values, such as recreation and stress relief.

Microbial communities in Pb contaminated boreal forest soil

Lead contamination in the environment is of particular concern, as it is a known toxin. Until recently, however, much less attention has been given to the local contamination caused by activities at shooting ranges compared to large-scale industrial contamination. In Finland, more than 500 tons of Pb is produced each year for shotgun ammunition. The contaminant threatens various organisms, ground water and the health of human populations. However, the forest at shooting ranges usually shows no visible sign of stress compared to nearby clean environments. The aboveground biota normally reflects the belowground ecosystem. Thus, the soil microbial communities appear to bear strong resistance to contamination, despite the influence of lead. The studies forming this thesis investigated a shooting range site at Hälvälä in Southern Finland, which is heavily contaminated by lead pellets. Previously it was experimentally shown that the growth of grasses and degradation of litter are retarded. Measurements of acute toxicity of the contaminated soil or soil extracts gave conflicting results, as enchytraeid worms used as toxicity reporters were strongly affected, while reporter bacteria showed no or very minor decreases in viability. Measurements using sensitive inducible luminescent reporter bacteria suggested that the bioavailability of lead in the soil is indeed low, and this notion was supported by the very low water extractability of the lead. Nevertheless, the frequency of lead-resistant cultivable bacteria was elevated based on the isolation of cultivable strains. The bacterial and fungal diversity in heavily lead contaminated shooting sectors were compared with those of pristine sections of the shooting range area. The bacterial 16S rRNA gene and fungal ITS rRNA gene were amplified, cloned and sequenced using total DNA extracted from the soil humus layer as the template. Altogether, 917 sequenced bacterial clones and 649 sequenced fungal clones revealed a high soil microbial diversity. No effect of lead contamination was found on bacterial richness or diversity, while fungal richness and diversity significantly differed between lead contaminated and clean control areas. However, even in the case of fungi, genera that were deemed sensitive were not totally absent from the contaminated area: only their relative frequency was significantly reduced. Some operational taxonomic units (OTUs) assigned to Basidiomycota were clearly affected, and were much rarer in the lead contaminated areas. The studies of this thesis surveyed EcM sporocarps, analyzed morphotyped EcM root tips by direct sequencing, and 454-pyrosequenced fungal communities in in-growth bags. A total of 32 EcM fungi that formed conspicuous sporocarps, 27 EcM fungal OTUs from 294 root tips, and 116 EcM fungal OTUs from a total of 8 194 ITS2 454 sequences were recorded. The ordination analyses by non-parametric multidimensional scaling (NMS) indicated that Pb enrichment induced a shift in the EcM community composition. This was visible as indicative trends in the sporocarp and root tip datasets, but explicitly clear in the communities observed in the in-growth bags. The compositional shift in the EcM community was mainly attributable to an increase in the frequencies of OTUs assigned to the genus Thelephora, and to a decrease in the OTUs assigned to Pseudotomentella, Suillus and Tylospora in Pb-contaminated areas when compared to the control. The enrichment of Thelephora in contaminated areas was also observed when examining the total fungal communities in soil using DNA cloning and sequencing technology. While the compositional shifts are clear, their functional consequences for the dominant trees or soil ecosystem remain undetermined. The results indicate that at the Hälvälä shooting range, lead influences the fungal communities but not the bacterial communities. The forest ecosystem shows apparent functional redundancy, since no significant effects were seen on forest trees. Recently, by means of 454 pyrosequencing , the amount of sequences in a single analysis run can be up to one million. It has been applied in microbial ecology studies to characterize microbial communities. The handling of sequence data with traditional programs is becoming difficult and exceedingly time consuming, and novel tools are needed to handle the vast amounts of data being generated. The field of microbial ecology has recently benefited from the availability of a number of tools for describing and comparing microbial communities using robust statistical methods. However, although these programs provide methods for rapid calculation, it has become necessary to make them more amenable to larger datasets and numbers of samples from pyrosequencing. As part of this thesis, a new program was developed, MuSSA (Multi-Sample Sequence Analyser), to handle sequence data from novel high-throughput sequencing approaches in microbial community analyses. The greatest advantage of the program is that large volumes of sequence data can be manipulated, and general OTU series with a frequency value can be calculated among a large number of samples.

Carabid beetles (Coleoptera, Carabidae) as indicators of environmental change in Ranomafana National Park, Madagascar

Growing human populations and increasing exploitation of natural resources threaten nature all over the world. Tropical countries are especially vulnerable to human impact because of the high number of species, most of these endemic and still unknown. Madagascar is one of the centers of high biodiversity and renowned for its unique species. However, during the last centuries many endemic species have gone extinct and more are endangered. Because of high natural values, Madagascar is one of the global conservation priorities. The establishment of Ranomafana National Park (RNP) was intended to preserve the unique nature of Madagascar. Containing several endemic and threatened species, Ranomafana has been selected as one of UNESCO’s World Natural Heritage sites. However, due to strong human pressures the region immediately surroundings the protected area has severely degraded. Aims of this thesis were to inventory carabid fauna in RNP and evaluate their use as indicators of the environmental change. Carabid beetles were collected from protected area (secondary and primary forests) and from its degraded surrounding area. Collecting was mostly conducted by hand during years 2000-2005. Species compositions between the protected area and its surroundings were compared, and species habitat preferences and seasonal variations were studied. In total, 4498 individuals representing 127 carabid species (of which 38 are new species) were collected. Species compositions within and outside of the protected area were markedly different. Most of the species preferred forest as their primary habitat and were mainly collected from trees and bushes. Their value as indicators is based on their different habitat requirements and sensitivity to environmental variables. Some of the species were found only in the protected forest, some occupied also the degraded forests and some preferred open areas. Carabid fauna is very species rich in Ranomafana and there are still many species to be found. Most of the species are arboreal and probably cannot survive in the deforested areas outside the park. This is very likely also the case for other species. Establishment and continued protection of RNP is probably the only way to conserve this globally important area. However, new occupations and land use methods are urgently needed by the local people for improving their own lives while maintaining the forest intact.

Detection, epidemiology and host spectrum of cowpox and Borna disease virus infections

Several orthopoxviruses (OPV) and Borna disease virus (BDV) are enveloped, zoonotic viruses with a wide geographical distribution. OPV antibodies cross-react, and former smallpox vaccination has therefore protected human populations from another OPV infection, rodent-borne cowpox virus (CPXV). Cowpox in humans and cats usually manifests as a mild, self-limiting dermatitis and constitutional symptoms, but it can be severe and even life-threatening in the immunocompromised. Classical Borna disease is a progressive meningoencephalomyelitis in horses and sheep known in central Europe for centuries. Nowadays the virus or its close relative infects humans and also several other species in central Europe and elsewhere, but the existence of human Borna disease with its suspected neuropsychiatric symptoms is controversial. The epidemiology of BDV is largely unknown, and the present situation is even more intriguing following the recent detection of several-million-year-old, endogenized BDV genes in primate and various other vertebrate genomes. The aims of this study were to elucidate the importance of CPXV and BDV in Finland and in possible host species, and particularly to 1) establish relevant methods for the detection of CPXV and other OPVs as well as BDV in Finland, 2) determine whether CPXV and BDV exist in Finland, 3) discover how common OPV immunity is in different age groups in Finland, 4) characterize possible disease cases and clarify their epidemiological context, 5) establish the hosts and possible reservoir species of these viruses and their geographical distribution in wild rodents, and 6) elucidate the infection kinetics of BDV in the bank vole. An indirect immunofluorescence assay and avidity measurement were established for the detection, timing and verification of OPV or BDV antibodies in thousands of blood samples from humans, horses, ruminants, lynxes, gallinaceous birds, dogs, cats and rodents. The mostly vaccine-derived OPV seroprevalence was found to decrease gradually according to the year of birth of the sampled human subjects from 100% to 10% in those born after 1977. On the other hand, OPV antibodies indicating natural contact with CPXV or other OPVs were commonly found in domestic and wild animals: the horse, cow, lynx, dog, cat and, with a prevalence occasionally even as high as 92%, in wild rodents, including some previously undetected species and new regions. Antibodies to BDV were detected in humans, horses, a dog, cats, and for the first time in wild rodents, such as bank voles (Myodes glareolus). Because of the controversy within the human Borna disease field, extra verification methods were established for BDV antibody findings: recombinant nucleocapsid and phosphoproteins were produced in Escherichia coli and in a baculovirus system, and peptide arrays were additionally applied. With these verification assays, Finnish human, equine, feline and rodent BDV infections were confirmed. Taken together, wide host spectra were evident for both OPV and BDV infections based on the antibody findings, and OPV infections were found to be geographically broadly distributed. PCR amplification methods were utilised for hundreds of blood and tissue samples. The methods included conventional, nested and real-time PCRs with or without the reverse transcription step and detecting four or two genes of OPVs and BDV, respectively. OPV DNA could be amplified from two human patients and three bank voles, whereas no BDV RNA was detected in naturally infected individuals. Based on the phylogenetic analyses, the Finnish OPV sequences were closely related although not identical to a Russian CPXV isolate, and clearly different from other CPXV strains. Moreover, the Finnish sequences only equalled each other, but the short amplicons obtained from German rodents were identical to monkeypox virus, in addition to German CPXV variants. This reflects the close relationship of all OPVs. In summary, RNA of the Finnish BDV variant could not be detected with the available PCR methods, but OPV DNA infrequently could. The OPV species infecting the patients of this study was proven to be CPXV, which is most probably also responsible for the rodent infections. Multiple cell lines and some newborn rodents were utilised in the isolation of CPXV and BDV from patient and wildlife samples. CPXV could be isolated from a child with severe, generalised cowpox. BDV isolation attempts from rodents were unsuccessful in this study. However, in parallel studies, a transient BDV infection of cells inoculated with equine brain material was detected, and BDV antigens discovered in archival animal brains using established immunohistology. Thus, based on several independent methods, both CPXV and BDV (or a closely related agent) were shown to be present in Finland. Bank voles could be productively infected with BDV. This experimental infection did not result in notable pathological findings or symptoms, despite the intense spread of the virus in the central and peripheral nervous system. Infected voles commonly excreted the virus in urine and faeces, which emphasises their possible role as a BDV reservoir. Moreover, BDV RNA was regularly reverse transcribed into DNA in bank voles, which was detected by amplifying DNA by PCR without reverse transcription, and verified with nuclease treatments. This finding indicates that BDV genes could be endogenized during an acute infection. Although further transmission studies are needed, this experimental infection demonstrated that the bank vole can function as a potential BDV reservoir. In summary, multiple methods were established and applied in large panels to detect two zoonoses novel to Finland: cowpox virus and Borna disease virus. Moreover, new information was obtained on their geographical distribution, host spectrum, epidemiology and infection kinetics.

Effects of large-scale human land use on Capercaillie (Tetrao urogallus L.) populations in Finland

The Capercaillie (Tetrao urogallus L.) is often used as a focal species for landscape ecological studies: the minimum size for its lekking area is 300 ha, and the annual home range for an individual may cover 30 80 km2. In Finland, Capercaillie populations have decreased by approximately 40 85%, with the declines likely to have started in the 1940s. Although the declines have partly stabilized from the 1990s onwards, it is obvious that the negative population trend was at least partly caused by changes in human land use. The aim of this thesis was to study the connections between human land use and Capercaillie populations in Finland, using several spatial and temporal scales. First, the effect of forest age structure on Capercaillie population trends was studied in 18 forestry board districts in Finland, during 1965 1988. Second, the abundances of Capercaillie and Moose (Alces alces L.) were compared in terms of several land-use variables on a scale of 50 × 50 km grids and in five regions in Finland. Third, the effects of forest cover and fine-grain forest fragmentation on Capercaillie lekking area persistence were studied in three study locations in Finland, on 1000 and 3000 m spatial scales surrounding the leks. The analyses considering lekking areas were performed with two definitions for forest: > 60 and > 152 m3ha 1 of timber volume. The results show that patterns and processes at large spatial scales strongly influence Capercaillie in Finland. In particular, in southwestern and eastern Finland, high forest cover and low human impact were found to be beneficial for this species. Forest cover (> 60 m3ha 1 of timber) surrounding the lekking sites positively affected lekking area persistence only at the larger landscape scale (3000 m radius). The effects of older forest classes were hard to assess due to scarcity of older forests in several study areas. Young and middle-aged forest classes were common in the vicinity of areas with high Capercaillie abundances especially in northern Finland. The increase in the amount of younger forest classes did not provide a good explanation for Capercaillie population decline in 1965 1988. In addition, there was no significant connection between mature forests (> 152 m3ha 1 of timber) and lekking area persistence in Finland. It seems that in present-day Finnish landscapes, area covered with old forest is either too scarce to efficiently explain the abundance of Capercaillie and the persistence of the lekking areas, or the effect of forest age is only important when considering smaller spatial scales than the ones studied in this thesis. In conclusion, larger spatial scales should be considered for assessing the future Capercaillie management. According to the proposed multi-level planning, the first priority should be to secure the large, regional-scale forest cover, and the second priority should be to maintain fine-grained, heterogeneous structure within the separate forest patches. A management unit covering hundreds of hectares, or even tens or hundreds of square kilometers, should be covered, which requires regional-level land-use planning and co-operation between forest owners.

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.

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.

Impact of the Human Bacterial Environment on Mycobacteriosis and Allergy

My work describes two sectors of the human bacterial environment: 1. The sources of exposure to infectious non-tuberculous mycobacteria. 2. Bacteria in dust, reflecting the airborne bacterial exposure in environments protecting from or predisposing to allergic disorders. Non-tuberculous mycobacteria (NTM) transmit to humans and animals from the environment. Infection by NTM in Finland has increased during the past decade beyond that by Mycobacterium tuberculosis. Among the farm animals, porcine mycobacteriosis is the predominant NTM disease in Finland. Symptoms of mycobacteriosis are found in 0.34 % of slaughtered pigs. Soil and drinking water are suspected as sources for humans and bedding materials for pigs. To achieve quantitative data on the sources of human and porcine NTM exposure, methods for quantitation of environmental NTM are needed. We developed a quantitative real-time PCR method, utilizing primers targeted at the 16S rRNA gene of the genus of Mycobacterium. With this method, I found in Finnish sphagnum peat, sandy soils and mud high contents of mycobacterial DNA, 106 to 107 genome equivalents per gram. A similar result was obtained by a method based on the Mycobacterium-specific hybridization of 16S rRNA. Since rRNA is found mainly in live cells, this result shows that the DNA detected by qPCR mainly represented live mycobacteria. Next, I investigated the occurrence of environmental mycobacteria in the bedding materials obtained from 5 pig farms with high prevalence (>4 %) of mycobacteriosis. When I used for quantification the same qPCR methods as for the soils, I found that piggery samples contained non-mycobacterial DNA that was amplified in spite of several mismatches with the primers. I therefore improved the qPCR assay by designing Mycobacterium-specific detection probes. Using the probe qPCR assay, I found 105 to 107 genome equivalents of mycobacterial DNA in unused bedding materials and up to 1000 fold more in the bedding collected after use in the piggery. This result shows that there was a source of mycobacteria in the bedding materials purchased by the piggery and that mycobacteria increased in the bedding materials during use in the piggery. Allergic diseases have reached epidemic proportions in urbanized countries. At the same time, childhood in rural environment or simple living conditions appears to protect against allergic disorders. Exposure to immunoreactive microbial components in rural environments seems to prevent allergies. I searched for differences in the bacterial communities of two indoor dusts, an urban house dust shown to possess immunoreactivity of the TH2-type and a farm barn dust with TH1-activity. The immunoreactivities of the dusts were revealed by my collaborators, in vitro in human dendritic cells and in vivo in mouse. The dusts accumulated >10 years in the respiratory zone (>1.5 m above floor), thus reflecting the long-term content of airborne bacteria at the two sites. I investigated these dusts by cloning and sequencing of bacterial 16S rRNA genes from dust contained DNA. From the TH2-active urban house dust, I isolated 139 16S rRNA gene clones. The most prevalent genera among the clones were Corynebacterium (5 species, 34 clones), Streptococcus (8 species, 33 clones), Staphylococcus (5 species, 9 clones) and Finegoldia (1 species, 9 clones). Almost all of these species are known as colonizers of the human skin and oral cavity. Species of Corynebacterium and Streptococcus have been reported to contain anti-inflammatory lipoarabinomannans and immunmoreactive beta-glucans respectively. Streptococcus mitis, found in the urban house dust is known as an inducer of TH2 polarized immunity, characteristic of allergic disorders. I isolated 152 DNA clones from the TH1-active farm barn dust and found species quite different from those found from the urban house dust. Among others, I found DNA clones representing Bacillus licheniformis, Acinetobacter lwoffii and Lactobacillus each of which was recently reported to possess anti-allergy immunoreactivity. Moreover, the farm barn dust contained dramatically higher bacterial diversity than the urban house dust. Exposure to this dust thus stimulated the human dendritic cells by multiple microbial components. Such stimulation was reported to promote TH1 immunity. The biodiversity in dust may thus be connected to its immunoreactivity. Furthermore, the bacterial biomass in the farm barn dust consisted of live intact bacteria mainly. In the urban house dust only ~1 % of the biomass appeared as intact bacteria, as judged by microscoping. Fragmented microbes may possess bioactivity different from that of intact cells. This was recently shown for moulds. If this is also valid for bacteria, the different immunoreactivities of the two dusts may be explained by the intactness of dustborne bacteria. Based on these results, we offer three factors potentially contributing to the polarized immunoreactivities of the two dusts: (i) the species-composition, (ii) the biodiversity and (iii) the intactness of the dustborne bacterial biomass. The risk of childhood atopic diseases is 4-fold lower in the Russian compared with the Finnish Karelia. This difference across the country border is not explainable by different geo-climatic factors or genetic susceptibilities of the two populations. Instead, the explanation must be lifestyle-related. It has already been reported that the microbiological quality of drinking water differs on the two sides of the borders. In collaboration with allergists, I investigated dusts collected from homes in the Russian Karelia and in the Finnish Karelia. I found that bacterial 16S rRNA genes cloned from the Russian Karelian dusts (10 homes, 234 clones) predominantly represented Gram-positive taxa (the phyla Actinobacteria and Firmicutes, 67%). The Russian Karelian dusts contained nine-fold more of muramic acid (60 to 70 ng mg-1) than the Finnish Karelian dusts (3 to 11 ng mg-1). Among the DNA clones isolated from the Finnish side (n=231), Gram-negative taxa (40%) outnumbered the Gram-positives (34%). Out of the 465 DNA clones isolated from the Karelian dusts, 242 were assigned to cultured validly described bacterial species. In Russian Karelia, animal-associated species e.g. Staphylococcus and Macrococcus were numerous (27 clones, 14 unique species). This finding may connect to the difference in the prevalence of allergy, as childhood contacts with pets and farm animals have been connected with low allergy risk. Plant-associated bacteria and plant-borne 16S rRNA genes (chloroplast) were frequent among the DNA clones isolated from the Finnish Karelia, indicating components originating from plants. In conclusion, my work revealed three major differences between the bacterial communtites in the Russian and in the Finnish Karelian homes: (i) the high prevalence of Gram-positive bacteria on the Russian side and of Gram-negative bacteria on the Finnish side and (ii) the rich presence of animal-associated bacteria on the Russian side whereas (iii) plant-associated bacteria prevailed on the Finnish side. One or several of these factors may connect to the differences in the prevalence of allergy.

Transcriptome and glycome profiling of human hematopoietic CD133+ and CD34+ cells

New blood cells are continuously provided by self-renewing multipotent hematopoietic stem cells (HSC). The capacity of HSCs to regenerate the hematopoietic system is utilized in the treatment of patients with hematological malignancies. HSCs can be enriched using an antibody-based recognition of CD34 or CD133 glycoproteins on the cell surface. The CD133+ and CD34+ cells may have partly different roles in hematopoiesis. Furthermore, each cell has a glycome typical for that cell type. Knowledge of HSC glycobiology can be used to design therapeutic cells with improved cell proliferation or homing properties. The present studies characterize the global gene expression profile of human cord blood-derived CD133+ and CD34+ cells, and demonstrate the differences between CD133+ and CD34+ cell populations that may have an impact in transplantation when CD133+ and CD34+ selected cells are used. In addition, these studies unravel the glycome profile of primitive hematopoietic cells and reveal the transcriptional regulation of N-glycan biosynthesis in CD133+ and CD34+ cells. The gene expression profile of CD133+ cells represents 690 differentially expressed transcripts between CD133+ cells and CD133- cells. CD34+ cells have 620 transcripts differentially expressed when compared to CD34- cells. The integrated CD133+/CD34+ cell gene expression profiles proffer novel transcripts to specify HSCs. Furthermore, the differences between the gene expression profiles of CD133+ and CD34+ cells indicate differences in the transcriptional regulation of CD133+ and CD34+ cells. CD133+ cells express a lower number of hematopoietic lineage differentiation marker genes than CD34+ cells. The expression profiles suggest a more primitive nature of CD133+ cells. Moreover, CD133+ cells have characteristic glycome that differ from the glycome of CD133- cells. High mannose-type and biantennary complex-type N-glycans are enriched in CD133+ cells. N-glycosylation-related gene expression pattern of CD133+ cells identify the key genes regulating the CD133+ cell-specific glycosylation including the overexpression of MGAT2 and underexpression of MGAT4. The putative role of MAN1C1 in the increase of unprocessed high mannose-type N-glycans in CD133+ cells is also discussed. These studies provide new information on the characteristics of HSCs. Improved understanding of HSC biology can be used to design therapeutic cells with improved cell proliferation and homing properties. As a result, HSC engineering could further their clinical use.