Phylogeography and Population Genetics of the Moor frog (Rana arvalis Nilsson) in Northern Europe

The ongoing climate change along with increasing levels of pollutants, diseases, habitat loss and fragmentation constitute global threats to the persistence of many populations, species and ecosystems. However, for the long-term persistence of local populations, one of the biggest threats is the intrinsic loss of genetic variation. In order to adapt to changes in the environment, organisms must have a sufficient supply of heritable variation in traits important for their fitness. With a loss of genetic variation, the risk of extinction will increase. For conservational practices, one should therefore understand the processes that shape the genetic population structure and also the broader (historical) phylogenetic patterning of the species in focus. In this thesis, microsatellite markers were applied to study genetic diversity and population differentiation of the protected moor frog (Rana arvalis) in Fennoscandia from both historical (evolutionary) and applied (conservation) perspectives. The results demonstrate that R. arvalis populations are highly structured over rather short geographic distances. Moreover, the results suggest that R. arvalis recolonized Fennoscandia from two directions after the last ice age. This has had implications for the genetic structuring and population differentiation, especially in the northernmost parts where the two lineages have met. Compared to more southern populations, the genetic variation decreases and the interpopulation differentiation increases dramatically towards north. This could be an outcome of serial population bottlenecking along the recolonization route. Also, current isolation and small population sizes increase the effect of drift, thus reinforcing the observed pattern. The same pattern can also be seen in island populations. However, though R. arvalis on the island of Gotland has lost most of its neutral genetic variability, our results indicate that the levels of additive genetic variation have remained high. This conforms to the conjecture that though neutral markers are widely used in conservation purposes, they may be quite uninformative about the levels of genetic variation in ecologically important traits. Finally, the evolutionary impact of the typical amphibian mating behaviour on genetic diversity was investigated. Given the short time available for larval development, it is important that mating takes place as early as possible. The genetic data and earlier capture-recapture data suggest that R. arvalis gather at mating grounds they are familiar with. However, by forming leks in random to relatedness, and having multiple paternities in single clutches, the risk of inbreeding may be minimized in this otherwise highly philopatric species.

Evaluation of methods and applications for behavioural profiling of transgenic mice

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

Genomic approach to organ-specific growth control

Growth is a fundamental aspect of life cycle of all organisms. Body size varies highly in most animal groups, such as mammals. Moreover, growth of a multicellular organism is not uniform enlargement of size, but different body parts and organs grow to their characteristic sizes at different times. Currently very little is known about the molecular mechanisms governing this organ-specific growth. The genome sequencing projects have provided complete genomic DNA sequences of several species over the past decade. The amount of genomic sequence information, including sequence variants within species, is constantly increasing. Based on the universal genetic code, we can make sense of this sequence information as far as it codes proteins. However, less is known about the molecular mechanisms that control expression of genes, and about the variations in gene expression that underlie many pathological states in humans. This is caused in part by lack of information about the second genetic code that consists of the binding specificities of transcription factors and the combinatorial code by which transcription factor binding sites are assembled to form tissue-specific and/or ligand-regulated enhancer elements. This thesis presents a high-throughput assay for identification of transcription factor binding specificities, which were then used to measure the DNA binding profiles of transcription factors involved in growth control. We developed ‘enhancer element locator’, a computational tool, which can be used to predict functional enhancer elements. A genome-wide prediction of human and mouse enhancer elements generated a large database of enhancer elements. This database can be used to identify target genes of signaling pathways, and to predict activated transcription factors based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c- and N-Myc genes, which has implications to organ specific growth control and tumor type specificity of oncogenes. Furthermore, we were able to locate a variation in a single nucleotide, which carries a susceptibility to colorectal cancer, to an enhancer element and propose a mechanism by which this SNP might be involved in generation of colorectal cancer.

The life cycle and genetic structure of the red alga Furcellaria lumbricalis on a salinity gradient

The life cycle and genetic diversity of the red alga Furcellaria lumbricalis (Hudson) Lamouroux were investigated in 15 populations in northern Europe. The occurrence of different life cycle phases and seasonality of reproduction were studied in four brackish populations in the northern Baltic Sea. Furthermore, a new method, based on genome screening with ISSR markers combined with a restriction-ligation method, was developed to discover microsatellite markers for population genetic analyses. The mitochondrial DNA cox2-3 spacer sequence and four microsatellite markers were used to examine the genetic diversity and differentiation of red algal populations in northern Europe. In addition, clonality and small-scale genetic structure of one Irish and four Baltic Sea populations were studied with microsatellite markers. It was discovered that at the low salinities of the northern Baltic Sea, only tetrasporophytes and males were present in the populations of F. lumbricalis and that winter was the main season for tetrasporangial production. Furthermore, the population occurring at the lowest salinity (3.6 practical salinity units, psu) did not produce spores. The size of the tetraspores was smaller in the Baltic Sea populations than that in the Irish population, and there were more deformed spores in the Baltic Sea populations than in the Irish populations. Studies with microsatellite markers indicated that clonality is a common phenomenon in the Baltic Sea populations of F. lumbricalis, although the proportion of clonal individuals varied among populations. Some genetic divergence occurred within locations both in Ireland and in the northern Baltic Sea. Even though no carpogonia were detected in the field samples during the study, the microsatellite data indicated that sexual reproduction occurs at least occasionally in the northern Baltic Sea. The genetic diversity of F. lumbricalis was highest in Brittany, France. Since no variation was discovered in the mtDNA cox2-3 spacer sequence, which is generally regarded as an informative phylogeographic marker in red algae, it can be assumed that the studied populations probably share the same origin.

Immune response to insulin and changes in the gut immune system in children with or at risk for type 1 diabetes

Type 1 diabetes (T1D) is considered to be an autoimmune disease. The cause of T1D is the destruction of insulin-producing β-cells in the pancreatic islets. The autoimmune nature of T1D is characterized by the presence of autoreactive T-cells and autoantibodies against β-cell molecules. Insulin is the only β-cell-specific autoantigen associated with T1D but the insulin autoantibodies (IAAs) are difficult to measure with proper sensitivity. T-cell assays for detection of autoreactive T-cells, such as insulin-specific T-cells, have also proven to be difficult to perform. The genetic risk of T1D is associated with the HLA gene region but the environmental factors also play an important role. The most studied environmental risk factors of T1D are enteroviruses and cow's milk which both affect the immune system through the gut. One hypothesis is that the insulin-specific immune response develops against bovine insulin in cow's milk during early infancy and later spreads to include human insulin. The aims of this study were to determine whether the separation of immunoglobulin (Ig)G from plasma would improve the sensitivity of the IAA assay and how insulin treatment affects the cellular immune response to insulin in newly diagnosed patients. Furthermore, the effect of insulin concentration in mother's breast milk on the development of antibodies to dietary insulin in the child was examined. Small intestinal biopsies were also obtained from children with T1D to characterize any immunological changes associated with T1D in the gut. The isolation of the IgG fraction from the plasma of T1D patients negative for plasma IAA led to detectable IAA levels that exceeded those in the control children. Thus the isolation of IgG may improve the sensitivity of the IAA assay. The effect of insulin treatment on insulin-specific T-cells was studied by culturing peripheral blood mononuclear cells with insulin. The insulin stimulation induced increased expression of regulatory T-cell markers, such as Foxp3, in those patients treated with insulin than in patients examined before initiating insulin treatment. This finding suggests that insulin treatment in patients with T1D stimulates regulatory T-cells in vivo and this may partly explain the difficulties in measuring autoantigen-specific T-cell responses in recently diagnosed patients. The stimulation of regulatory T-cells by insulin treatment may also explain the remission period often seen after initiating insulin treatment. In the third study we showed that insulin concentration in mother's breast milk correlates inversely with the levels of bovine insulin-specific antibodies in those infants who were exposed to cow's milk proteins in their diet, suggesting that human insulin in breast milk induces tolerance to dietary bovine insulin. However, in infants who later developed T1D-associated autoantibodies, the insulin concentration in their mother's breast milk was increased. This finding may indicate that in those children prone to β-cell autoimmunity, breast milk insulin does not promote tolerance to insulin. In the small intestinal biopsies the presence of several immunological markers were quantified with the RT-PCR. From these markers the expression of the interleukin (IL)-18 cytokine was significantly increased in the gut in patients with T1D compared with children with celiac disease or control children. The increased IL-18 expression lends further support for the hypothesis that the gut immune system is involved in the pathogenesis of T1D.

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.

Evolutionary Genetics in the WIld : from Populations to Individuals

Predicting evolutionary outcomes and reconstructing past evolutionary transitions are among the main goals of evolutionary biology. Ultimately, understanding the mechanisms of evolutionary change will also provide answers to the timely question of whether and how organisms will adapt to changing environmental conditions. In this thesis, I have investigated the relative roles of natural selection, random genetic drift and genetic correlations in the evolution of complex traits at different levels of organisation from populations to individuals. I have shown that natural selection has been the driving force behind body shape divergence of marine and freshwater threespine stickleback (Gasterosteus aculeatus) populations, while genetic drift may have played a significant role in the more fine scale divergence among isolated freshwater populations. These results are concurrent with the patterns that have emerged in the published studies comparing the relative importance of natural selection and genetic drift as explanations for population divergence in different traits and taxa. I have also shown that body shape and armour divergence among threespine stickleback populations is likely to be biased by the patterns of genetic variation and covariation. Body shape and armour variation along the most likely direction of evolution the direction of maximum genetic variance reflects the general patterns of variation observed wild populations across the distribution range of the threespine stickleback. Conversely, it appears that genetic correlations between the sexes have not imposed significant constraints on the evolution of sexual dimorphism in threespine stickleback body shape and armour. I have demonstrated that the patterns of evolution seen in the wild can be experimentally recreated to tease out the effects of different selection agents in detail. In addition, I have shown how important it is to take into account the correlative nature of traits, when making interpretations about the effects of natural selection on individual traits. Overall, this thesis provides a demonstration of how considering the relative roles of different mechanism of evolutionary change at different levels of organisation can aid in an emergence of a comprehensive picture of how adaptive divergence in wild populations occurs.

The effects of eutrophication on alternative reproductive tactics in threespine sticklebacks

When organisms compete for mates and fertilisations, the process of sexual selection drives the evolution of traits that increase reproductive success. The traits targeted by selection, and the extent to which they change, are constrained by the local environment. Sexual selection due to female mate choice can be undermined by alternative reproductive tactics (ARTs), which refers to discontinuous variation in traits or behaviours used in reproduction. As human activities are rapidly changing our planet, this raises the question how ARTs will be affected. Fish show a bewildering diversity of ARTs, which make them good model organisms to answer these questions. One example of human-induced environmental change, which is affecting aquatic ecosystems around the world, is eutrophication, the over-enrichment of water bodies with nutrients. One of its effects is decreased underwater visibility due to increases in both turbidity and vegetation density. The aims of this thesis were to investigate the effects increased turbidity and vegetation density have on an ART in sticklebacks, a fish common to marine and fresh water bodies of the Northern hemisphere. I furthermore investigated how this affected sexual selection for male size, a trait commonly under selection. I used a combination of behavioural observations in microcosms, where I manipulated underwater visibility, with collection of genetic material to reconstruct parentage of broods, and thus identify sneak fertilisations. The results show that turbidity might have weak negative effects on the frequency of sneaking behaviour, although this behaviour was rather infrequent in these experiments, which complicates firm conclusions. In dense vegetation the number of sneak fertilisations decreased slightly, as fewer nesting males sneaked, while the number of non-nesting males sneaking remained constant. The paternity analyses revealed that a significantly smaller fraction of eggs was sneak fertilised under dense vegetation. Furthermore, amongst the nesting males that sneaked, the amount of eggs sneak fertilised correlated positively with courtship success. A reduction in sneaking by these males under dense vegetation equalised the distribution of fertilisation success, in turn contributing to a decrease in the opportunity for selection. Under dense vegetation significantly more males built nests, which has also been observed in previous field studies. In a separate experiment we addressed if such changes in the proportion of nesters and non-nesters, without changes in visibility, affected the incidence of sneak fertilisation. My results show this was not the case, likely because sneaking is an opportunistic tactic shown by both nesters and non-nesters. Non-nesters did sneak proportionately more when there were many of them, which could be due to changes in the cost-benefit ratio of sneaking. As nesters can only attack one intruder at a time, the costs and risks per sneaker will decrease as the number of sneakers increases. The defensive behaviours shown by the nesters before spawning shifted to a more aggressive form of nest defence. This could be because less aggressive behaviours lose their effectiveness when the number of intruders increases. It could also indicate that the risks associated with aggressive behaviours decrease when there are fewer fellow nesters, as other studies indicate nesters are competitive and aggressive individuals. Under turbid conditions I did not detect changes in the opportunity for selection, based on fertilisation success, nor was male size under significant selection under clear or turbid conditions. More thorough analyses under densely vegetated conditions across the nesting, courtship and fertilisation stages revealed a decrease in the opportunity for selection across all stages. A reduction in sneaking by nesters contributed to this. During the nesting stage, but not during later stages, body size was under significant directional selection under sparse, but not dense vegetation. This illustrates the importance of considering all selection stages to get a complete picture of how environmental changes affect sexual selection. Leaving out certain stages or subgroups can result in incomplete or misleading results.

The Emerging Phenotype: Ecological Genetics of Color, Form and Sex in the Common Frog

One of the main aims of evolutionary biology is to explain why organisms vary phenotypically as they do. Proximately, this variation arises from genetic differences and from environmental influences, the latter of which is referred to as phenotypic plasticity. Phenotypic plasticity is thus a central concept in evolutionary biology, and understanding its relative importance in causing the phenotypic variation and differentiation is important, for instance in anticipating the consequences of human induced environmental changes. The aim of this thesis was to study geographic variation and local adaptation, as well as sex ratios and environmental sex reversal, in the common frog (Rana temporaria). These themes cover three different aspects of phenotypic plasticity, which emerges as the central concept for the thesis. The first two chapters address geographic variation and local adaptation in two potentially thermally adaptive traits, namely the degree of melanism and the relative leg length. The results show that although there is an increasing latitudinal trend in the degree of melanism in wild populations across Scandinavian Peninsula, this cline has no direct genetic basis and is thus environmentally induced. The second chapter demonstrates that although there is no linear, latitudinally ordered phenotypic trend in relative leg length that would be expected under Allen s rule an ecogeographical rule linking extremity length to climatic conditions there seems to be such a trend at the genetic level, hidden under environmental effects. The first two chapters thus view phenotypic plasticity through its ecological role and evolution, and demonstrate that it can both give rise to phenotypic variation and hide evolutionary patterns in studies that focus solely on phenotypes. The last three chapters relate to phenotypic plasticity through its ecological and evolutionary role in sex determination, and consequent effects on population sex ratio, genetic recombination and the evolution of sex chromosomes. The results show that while sex ratios are strongly female biased and there is evidence of environmental sex reversals, these reversals are unlikely to have caused the sex ratio skew, at least directly. The results demonstrate that environmental sex reversal can have an effect on the evolution of sex chromosomes, as the recombination patterns between them seem to be controlled by phenotypic, rather than genetic, sex. This potentially allows Y chromosomes to recombine, lending support for the recent hypothesis suggesting that sex-reversal may play an important role on the rejuvenation of Y chromosomes.

Ecology, population genetics and conservation of the African violet

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

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.

Phylogeography and hybrid swarms: history of brackish water bivalve diversity in North European marginal seas

This study addressed the large-scale molecular zoogeography in two brackish water bivalve molluscs, Macoma balthica and Cerastoderma glaucum, and genetic signatures of the postglacial colonization of Northern Europe by them. The traditional view poses that M. balthica in the Baltic, White and Barents seas (i.e. marginal seas) represent direct postglacial descendants of the adjacent Northeast Atlantic populations, but this has recently been challenged by observations of close genetic affinities between these marginal populations and those of the Northeast Pacific. The primary aim of the thesis was to verify, quantify and characterize the Pacific genetic contribution across North European populations of M. balthica and to resolve the phylogeographic histories of the two bivalve taxa in range-wide studies using information from mitochondrial DNA (mtDNA) and nuclear allozyme polymorphisms. The presence of recent Pacific genetic influence in M. balthica of the Baltic, White and Barents seas, along with an Atlantic element, was confirmed by mtDNA sequence data. On a broader temporal and geographical scale, altogether four independent trans-Arctic invasions of Macoma from the Pacific since the Miocene seem to have been involved in generating the current North Atlantic lineage diversity. The latest trans-Arctic invasion that affected the current Baltic, White and Barents Sea populations probably took place in the early post-glacial. The nuclear genetic compositions of these marginal sea populations are intermediate between those of pure Pacific and Atlantic subspecies. In the marginal sea populations of mixed ancestry (Barents, White and Northern Baltic seas), the Pacific and Atlantic components are now randomly associated in the genomes of individual clams, which indicates both pervasive historical interbreeding between the previously long-isolated lineages (subspecies), and current isolation of these populations from the adjacent pure Atlantic populations. These mixed populations can be characterized as self-supporting hybrid swarms, and they arguably represent the most extensive marine animal hybrid swarms so far documented. Each of the three swarms still has a distinct genetic composition, and the relative Pacific contributions vary from 30 to 90 % in local populations. This diversity highlights the potential of introgressive hybridization to rapidly give rise to new evolutionarily and ecologically significant units in the marine realm. In the south of the Danish straits and in the Southern Baltic Sea, a broad genetic transition zone links the pure North Sea subspecies M. balthica rubra to the inner Baltic hybrid swarm, which has about 60 % of Pacific contribution in its genome. This transition zone has no regular smooth clinal structure, but its populations show strong genotypic disequilibria typical of a hybrid zone maintained by the interplay of selection and gene flow by dispersing pelagic larvae. The structure of the genetic transition is partly in line with features of Baltic water circulation and salinity stratification, with greater penetration of Atlantic genes on the Baltic south coast and in deeper water populations. In all, the scenarios of historical isolation and secondary contact that arise from the phylogeographic studies of both Macoma and Cerastoderma shed light to the more general but enigmatic patterns seen in marine phylogeography, where deep genetic breaks are often seen in species with high dispersal potential.

Molecular details of phage phi6 RNA-dependent RNA synthesis

For most RNA viruses RNA-dependent RNA polymerases (RdRPs) encoded by the virus are responsible for the entire RNA metabolism. Thus, RdRPs are critical components in the viral life cycle. However, it is not fully understood how these important enzymes function during viral replication. Double-stranded RNA (dsRNA) viruses perform the synthesis of their RNA genome within a proteinacous viral particle containing an RdRP as a minor constituent. The phi6 bacteriophage is the best-studied dsRNA virus, providing an excellent background for studies of its RNA synthesis. The purified recombinant phi6 RdRP is highly active in vitro and it possesses both RNA replication and transcription activities. The crystal structure of the phi6 polymerase, solved in complex with a number of ligands, provides a working model for detailed in vitro studies of RNA-dependent RNA polymerization. In this thesis, the primer-independent initiation of the phi6 RdRP was studied in vitro using biochemical and structural methods. A C-terminal, four-amino-acid-long loop protruding into the central cavity of the phi6 RdRP has been suggested to stabilize the incoming nucleotides of the initiation complex formation through stacking interactions. A similar structural element has been found from several other viral RdRPs. In this thesis, this so-called initiation platform loop was subjected to site-directed mutagenesis to address its role in the initiation. It was found that the initiation mode of the mutants is primer-dependent, requiring either an oligonucleotide primer or a back-priming initiation mechanism for the RNA synthesis. The crystal structure of a mutant RdRP with altered initiation platform revealed a set of contacts important for primer-independent initiation. Since phi6 RdRP is structurally and functionally homologous to several viral RdRPs, among them the hepatitis C virus RdRP, these results provide further general insight to understand primer-independent initiation. In this study it is demonstrated that manganese phasing could be used as a practical tool for solving structures of large proteins with a bound manganese ion. The phi6 RdRP was used as a case study to obtain phases for crystallographic analysis. Manganese ions are naturally bound to the phi6 RdRP at the palm domain of the enzyme. In a crystallographic experiment, X-ray diffraction data from a phi6 RdRP crystal were collected at a wavelength of 1.89 Å, which is the K edge of manganese. With this data an automatically built model of the core region of the protein could be obtained. Finally, in this work terminal nucleotidyl transferase (TNTase) activity of the phi6 RdRP was documented in the isolated polymerase as well as in the viral particle. This is the first time that such an activity has been reported in a polymerase of a dsRNA virus. The phi6 RdRP used uridine triphosphates as the sole substrate in a TNTase reaction but could accept several heterologous templates. The RdRP was able to add one or a few non-templated nucleotides to the 3' end of the single- or double-stranded RNA substrate. Based on the results on particle-mediated TNTase activity and previous structural information of the polymerase, a model for termination of the RNA-dependent RNA synthesis is suggested in this thesis.

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

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

Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.