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.

Behavioural and physiological responses to predators of captive-bred Arctic charr: significance of genetics, learning and ontogeny

Predation forms one of the main selective forces in nature and in a vast number of prey species the behavioural responses form the main way to avoid predation. World wide numerous captive breeding programs are used to produce fish and other animal species for conservational reintroductions. However, rearing animals in the absence of predators in captivity has been shown to weaken their predator avoidance skills and lead to behavioural divergence between wild and captive-bred populations. In my thesis I studied the effects of predator odour exposures on antipredator behavioural and physiological responses of captive reared Saimaa Arctic charr. This charr population is the most endangered fish population in Finland and a sample of the remaining population has been taken to captive breeding and used for an extensive reintroduction program. Lowered responsiveness to predators is probably one of the major reasons for the poor survival probability of the charr after release into the wild. The main aims of my thesis were to explore the reasons for behavioural phenotypic variation in this charr population and whether naïve charr young could be trained to recognise their natural predators. The predator species in my thesis were burbot (Lota lota) and pikeperch (Sander lucioperca). In my thesis I showed that the captive-bred charr responded to chemical cues from burbot and pikeperch, but the magnitude of responses was linked to the predator species. The burbot odour increased the spatial odour avoidance of the charr young. On the other hand, in the pikeperch treatment charr reduced their relative swimming activity and tended to show more freezing behaviour relative to the burbot treatment. It seems evident that these different responses are related to the different hunting tactics of predator species. Furthermore, I detected wide between-family differences in antipredator responsiveness (i.e. inherited variation in antipredator behaviours) in this captive stock. Detected differences were greater in the response towards burbot than towards pikeperch. These results, in addition to predator-specific antipredator responses, suggest that there is a clear inherited component in antipredator responsiveness in Saimaa charr population and that the detected inherited differences could explain a part of the behavioural phenotypic variation in this population. In my thesis I also found out that both social learning and direct exposure to live predators enhance the antipredator responsiveness of charr young. In addition, I obtained indications that predator odour exposures (i.e. life-skills training) in alevin and fry stages can fine-tune the innate antipredator responsiveness of charr. Thus, all these methods have the potential to enhance the innate antipredator responsiveness of naïve charr young, possibly also improving the post-release survival of these trained individuals in the wild. However, the next logical phase would be to carry out large scale survival studies in the wild to test this hypothesis. Finally, the results of my thesis emphasize that possible long-term life-skills training methods should take into account not only the behavioural but also the physiological effects of training.

Molecular Tuning of Rhodopsins for High Visual Sensitivity in Different Light Environments : Variation in Absorbance Spectrum and Opsin Sequence within and between Species

Visual pigments of different animal species must have evolved at some stage to match the prevailing light environments, since all visual functions depend on their ability to absorb available photons and transduce the event into a reliable neural signal. There is a large literature on correlation between the light environment and spectral sensitivity between different fish species. However, little work has been done on evolutionary adaptation between separated populations within species. More generally, little is known about the rate of evolutionary adaptation to changing spectral environments. The objective of this thesis is to illuminate the constraints under which the evolutionary tuning of visual pigments works as evident in: scope, tempo, available molecular routes, and signal/noise trade-offs. Aquatic environments offer Nature s own laboratories for research on visual pigment properties, as naturally occurring light environments offer an enormous range of variation in both spectral composition and intensity. The present thesis focuses on the visual pigments that serve dim-light vision in two groups of model species, teleost fishes and mysid crustaceans. The geographical emphasis is in the brackish Baltic Sea area with its well-known postglacial isolation history and its aquatic fauna of both marine and fresh-water origin. The absorbance spectrum of the (single) dim-light visual pigment were recorded by microspectrophotometry (MSP) in single rods of 26 fish species and single rhabdoms of 8 opossum shrimp populations of the genus Mysis inhabiting marine, brackish or freshwater environments. Additionally, spectral sensitivity was determined from six Mysis populations by electroretinogram (ERG) recording. The rod opsin gene was sequenced in individuals of four allopatric populations of the sand goby (Pomatoschistus minutus). Rod opsins of two other goby species were investigated as outgroups for comparison. Rod absorbance spectra of the Baltic subspecies or populations of the primarily marine species herring (Clupea harengus membras), sand goby (P. minutus), and flounder (Platichthys flesus) were long-wavelength-shifted compared to their marine populations. The spectral shifts are consistent with adaptation for improved quantum catch (QC) as well as improved signal-to-noise ratio (SNR) of vision in the Baltic light environment. Since the chromophore of the pigment was pure A1 in all cases, this has apparently been achieved by evolutionary tuning of the opsin visual pigment. By contrast, no opsin-based differences were evident between lake and sea populations of species of fresh-water origin, which can tune their pigment by varying chromophore ratios. A more detailed analysis of differences in absorbance spectra and opsin sequence between and within populations was conducted using the sand goby as model species. Four allopatric populations from the Baltic Sea (B), Swedish west coast (S), English Channel (E), and Adriatic Sea (A) were examined. Rod absorbance spectra, characterized by the wavelength of maximum absorbance (λmax), differed between populations and correlated with differences in the spectral light transmission of the respective water bodies. The greatest λmax shift as well as the greatest opsin sequence difference was between the Baltic and the Adriatic populations. The significant within-population variation of the Baltic λmax values (506-511 nm) was analyzed on the level of individuals and was shown to correlate well with opsin sequence substitutions. The sequences of individuals with λmax at shorter wavelengths were identical to that of the Swedish population, whereas those with λmax at longer wavelengths additionally had substitution F261F/Y in the sixth transmembrane helix of the protein. This substitution (Y261) was also present in the Baltic common gobies and is known to redshift spectra. The tuning mechanism of the long-wavelength type Baltic sand gobies is assumed to be the co-expression of F261 and Y261 in all rods to produce ≈ 5 nm redshift. The polymorphism of the Baltic sand goby population possibly indicates ambiguous selection pressures in the Baltic Sea. The visual pigments of all lake populations of the opossum shrimp (Mysis relicta) were red-shifted by 25 nm compared with all Baltic Sea populations. This is calculated to confer a significant advantage in both QC and SNR in many humus-rich lakes with reddish water. Since only A2 chromophore was present, the differences obviously reflect evolutionary tuning of the visual protein, the opsin. The changes have occurred within the ca. 9000 years that the lakes have been isolated from the Sea after the most recent glaciation. At present, it seems that the mechanism explaining the spectral differences between lake and sea populations is not an amino acid substitution at any other conventional tuning site, but the mechanism is yet to be found.

Phytoplankton ecology of subarctic lakes in Finnish Lapland

Climate is warming and it is especially seen in arctic areas, where the warming trend is expected to be greatest. Arctic freshwater ecosystems, which are a very characteristic feature of the arctic landscape, are especially sensitive to climate change. They could be used as early warning systems, but more information about the ecosystem functioning and responses are needed for proper interpretation of the observations. Phytoplankton species and assemblages could be especially suitable for climate-related studies, since they have short generation times and react rapidly to changes in the environment. In addition, phytoplankton provides a good tool for lake classifications, since different species have different requirements and tolerance ranges for various environmental factors. The use of biological indicators is especially useful in arctic areas, were many of the chemical factors commonly fall under the detection limit and therefore do not provide much information about the environment. This work brings new information about species distribution and dynamics of arctic freshwater phytoplankton in relation to environmental factors. The phytoplankton of lakes in Finnish Lapland and other European high-altitude or high-latitude areas were compared. Most lakes were oligotrophic and dominated by flagellated species belonging to chrysophytes, cryptophytes and dinoflagellates. In Finnish Lapland cryptophytes were of less importance, whereas desmids had high species richness in many of the lakes. In Pan-European scale, geographical and catchment-related factors were explaining most of the differences in species distributions between different districts, whereas lake water chemistry (especially conductivity, SiO2 and pH) was most important regionally. Seasonal and interannual variation of phytoplankton was studied in subarctic Lake Saanajärvi. Characteristic phytoplankton species in this oligotrophic, dimictic lake belonged mainly to chrysophytes and diatoms. The maximum phytoplankton biomass in Lake Saanajärvi occurs during autumn, while spring biomass is very low. During years with heavy snow cover the lake suffers from pH drop caused by melt waters, but the effects of this acid pulse are restricted to surface layers and last for a relatively short period. In addition to some chemical parameters (mainly Ca and nutrients), length of the mixing cycle and physical factors such as lake water temperature and thermal stability of water column had major impact on phytoplankton dynamics. During a year with long and strong thermal stability, the phytoplankton community developed towards an equilibrium state, with heavy dominance of only a few taxa for a longer period of time. During a year with higher windiness and less thermal stability, the species composition was more diverse and species with different functional strategies were able to occur simultaneously. The results of this work indicate that although arctic lakes in general share many common features concerning their catchment and water chemistry, large differences in biological features can be found even in a relatively small area. Most likely the lakes with very different algal flora do not respond in a similar way to differences in the environmental factors, and more information about specific arctic lake types is needed. The results also show considerable year to year differences in phytoplankton species distribution and dynamics, and these changes are most likely linked to climatic factors.

Cereal alkylresorcinols as dietary biomarkers-absorption and occurrence in biological membranes

Several studies link the consumption of whole-grain products to a lowered risk of chronic diseases, such as certain types of cancer, type II diabetes, and cardiovascular diseases. However, the final conclusions of the exact protective mechanisms remain unclear, partly due to a lack of a suitable biomarker for the whole-grain cereals intake. Alkylresorcinols (AR) are phenolic lipids abundant in the outer parts of wheat and rye grains usually with homologues of C15:0- C25:0 alkyl chains, and are suggested to function as whole-grain biomarkers. Mammalian lignan enterolactone has also previously been studied as a potential whole-grain biomarker. In the present work a quantified gas chromatography-mass spectrometry method for the analysis of AR in plasma, erythrocytes, and lipoproteins was developed. The method was used to determine human and pig plasma AR concentrations after the intake of whole-grain wheat and rye products compared to low-fibre wheat bread diets to assess the usability of AR as biomarkers of whole-grain intake. AR plasma concentrations were compared to serum ENL concentrations. AR absorption and elimination kinetics were investigated in a pig model. AR occurrence in human erythrocyte membranes and plasma lipoproteins were determined, and the distribution of AR in blood was evaluated. Plasma AR seem to be absorbed via the lymphatic system from the small intestine, like many other lipophilic compounds. Their apparent elimination half-life is relatively short and is similar to that of tocopherols, which have a similar chemical structure. Plasma AR concentrations increased significantly after a one- to eight-week intake of whole-grain wheat and further on with whole-grain rye bread. The concentrations were also higher after habitual Finnish diet compared to diet with low-fibre bread. Inter-individual variation after a one-week intake of the same amount of bread was high, but the mean plasma AR concentrations increased with increasing AR intake. AR are incorporated into erythrocyte membranes and plasma lipoproteins, and VLDL and HDL were the main AR carriers in human plasma. Based on these studies, plasma AR could function as specific biomarkers of dietary whole-grain products. AR are exclusively found in whole-grains and are more suitable as specific biomarkers of whole-grain intake than previously investigated mammalian lignan enterolactone, that is formed from several plants other than cereals in the diet. Plasma AR C17:0/C21:0 -ratio could distinguish whether whole-grain products in the diet are mainly wheat or rye. AR could be used in epidemiological studies to determine whole-grain intake and to better assess the role of whole-grains in disease prevention.

Individual variation in in vitro efficacy of antiplatelet medication - aspirin and clopidogrel

Antiplatelet medication is known to decrease adverse effects in patients with atherothrombotic disease. However, despite ongoing antiplatelet medication considerable number of patients suffer from atherothrombotic events. The aims of the study were 1) to evaluate the individual variability in platelet functions and compare the usability of different methods in detecting it, 2) to assess variability in efficacy of antiplatelet medication with aspirin (acetylsalicylic acid) or the combination of aspirin and clopidogrel and 3) to investigate the main genetic and clinical variables as well as potential underlying mechanisms of variability in efficacy of antiplatelet medication. In comparisons of different platelet function tests in 19 healthy individuals PFA-100® correlated with traditional methods of measuring platelet function and was thus considered appropriate for testing individual variability in platelet activity. Efficacy of ongoing 100mg aspirin daily was studied in 101 patients with coronary artery disease (CAD). Aspirin response was measured with arachidonic acid (AA)-induced platelet aggregation, which reflects cyclo-oxygenase (COX)-1 dependent thromboxane (Tx) A2 formation, and PFA-100®, which evaluates platelet activation under high shear stress in the presence of collagen and epinephrine. Five percent of patients failed to show inhibition of AA-aggregation and 21% of patients had normal PFA-100® results despite aspirin and were thus considered non-responders to aspirin. Interestingly, the two methods of assessing aspirin efficacy, platelet aggregation and PFA-100®, detected different populations as being aspirin non-responders. It could be postulated that PFA-100® actually measures enhanced platelet function, which is not directly associated with TxA2 inhibition exerted by aspirin. Clopidogrel efficacy was assessed in 50 patients who received a 300mg loading dose of clopidogrel 2.5 h prior to percutaneous coronary intervention (PCI) and in 51 patients who were given a loading dose of 300mg combined with a five day treatment of 75mg clopidogrel daily mimicking ongoing treatment. Clopidogrel response was assessed with ADP-induced aggregations, due to its mechanism of action as an inhibitor of ADP-induced activation. When patients received only a loading dose of clopidogrel prior to PCI, 40% did not gain measurable inhibition of their ADP-induced platelet activity (inhibition of 10% or less). Prolongation of treatment so that all patients had reached a plateau of inhibition exerted by clopidogrel, decreased the incidence of non-responders to 20%. Polymorphisms of COX-1 and GP VI, as well as diabetes and female gender, were associated with decreased in vitro aspirin efficacy. Diabetes also impaired the in vitro efficacy of short-term clopidogrel. Decreased response to clopidogrel was associated with limited inhibition by ARMX, an antagonist of P2Y12-receptor, suggesting the reason for clopidogrel resistance to be receptor-dependent. Conclusions: Considerable numbers of CAD patients were non-responders either to aspirin, clopidogrel or both. In the future, platelet function tests may be helpful to individually select effective and safe antiplatelet medication for these patients.

Genetic Variation and Clinical Manifestations in Inflammatory Bowel Disease

Crohn s disease (CD) and ulcerative colitis (UC), collectively known as inflammatory bowel disease (IBD), are characterised by chronic inflammation of the gastrointestinal tract. IBD prevalence in Finland is approximately 3-4 per 1000 inhabitants with a peak incidence in adolescence. The symptoms of IBD include diarrhoea, abdominal pain, fever, and weight loss. The precise aetiology of IBD is unknown but interplay of environmental risk factors and immunologic changes trigger the disease in a genetically susceptible individual. Twin and family studies have provided strong evidence for genetic factors in IBD susceptibility, and genetic factors may be more prominent in CD than UC. The first CD susceptibility gene was identified in 2001. Three common mutations R702W, G908R, and 1007fs of the CARD15/NOD2 gene are shown to associate independently with CD but the magnitude of association varies between different populations. The present study aimed at identifying mutations and genetic variations in IBD susceptibility and candidate genes. In addition, correlation to phenotype was also assessed. One of the main objectives of this study was to evaluate the role of CARD15 in a Finnish CD cohort. 271 CD patients were studied for the three common mutations and the results showed a lower mutation frequency than in other Caucasian populations. Only 16% of the patients carried one of the three mutations. Ileal location as well as stricturing and penetrating behaviour of the disease were associated with occurrence of the mutations. The whole protein coding region of CARD15 was screened for possible Finnish founder mutations. In addition to several sequence variants, five novel mutations (R38M, W355X, P727L, W907R, and R1019X) were identified in five patients. Functional consequences of these novel variants were studied in vitro, and these studies demonstrated a profound impairment of MDP response. Investigation of CARD15 mutation frequency in healthy people across three continents showed a large geographic fluctuation. No simple correlation between mutation frequency and disease incidence was seen in populations studied. The occurrence of double mutant carriers in healthy controls suggested that the penetrance of risk alleles is low. Other main objectives aimed at identifying other genetic variations that are involved in the susceptibility to IBD. We investigated the most plausible IBD candidate genes including TRAF6, SLC22A4, SLC22A5, DLG5, TLR4, TNFRSF1A, ABCB1/MDR1, IL23R, and ATG16L1. The marker for a chromosome 5 risk haplotype and the rare HLA-DRB1*0103 allele were also studied. The study cohort consisted of 699 IBD patients (240 CD and 459 UC), of which 23% had a first-degree relative with IBD. Of the several candidate genes studied, IL23R was associated with CD susceptibility, and TNFRSF1A as well as the HLA-DRB1*0103 allele with UC susceptibility. IL23R variants also showed association with the stricturing phenotype and longer disease duration in CD patients. In addition, TNFRSF1A variants were more common among familial UC and ileocolonic CD. In conclusion, the common CARD15 mutations were shown to account for 16% of CD cases in Finland. Novel CARD15 variants identified in the present study are most likely disease-causing mutations, as judged by the results of in vitro studies. The present study also confirms the IL23R association with CD susceptibility and, in addition, TNFRSF1A and HLA-DRB1*0103 allele association with UC of specific clinical phenotypes.

Ionosphere-atmosphere interaction during solar proton events

Among the most striking natural phenomena affecting ozone are solar proton events (SPE), during which high-energy protons precipitate into the middle atmosphere in the polar regions. Ionisation caused by the protons results in changes in the lower ionosphere, and in production of neutral odd nitrogen and odd hydrogen species which then destroy ozone in well-known catalytic chemical reaction chains. Large SPEs are able to decrease the ozone concentration of upper stratosphere and mesosphere, but are not expected to significantly affect the ozone layer at 15--30~km altitude. In this work we have used the Sodankylä Ion and Neutral Chemistry Model (SIC) in studies of the short-term effects caused by SPEs. The model results were found to be in a good agreement with ionospheric observations from incoherent scatter radars, riometers, and VLF radio receivers as well as with measurements from the GOMOS/Envisat satellite instrument. For the first time, GOMOS was able to observe the SPE effects on odd nitrogen and ozone in the winter polar region. Ozone observations from GOMOS were validated against those from MIPAS/Envisat instrument, and a good agreement was found throughout the middle atmosphere. For the case of the SPE of October/November 2003, long-term ozone depletion was observed in the upper stratosphere. The depletion was further enhanced by the descent of odd nitrogen from the mesosphere inside the polar vortex, until the recovery occurred in late December. During the event, substantial diurnal variation of ozone depletion was seen in the mesosphere, caused mainly by the the strong diurnal cycle of the odd hydrogen species. In the lower ionosphere, SPEs increase the electron density which is very low in normal conditions. Therefore, SPEs make radar observations easier. In the case of the SPE of October, 1989, we studied the sunset transition of negative charge from electrons to ions, a long-standing problem. The observed phenomenon, which is controlled by the amount of solar radiation, was successfully explained by considering twilight changes in both the rate of photodetachment of negative ions and concentrations of minor neutral species. Changes in the magnetic field of the Earth control the extent of SPE-affected area. For the SPE of November 2001, the results indicated that for low and middle levels of geomagnetic disturbance the estimated cosmic radio noise absorption levels based on a magnetic field model are in a good agreement with ionospheric observations. For high levels of disturbance, the model overestimates the stretching of the geomagnetic field and the geographical extent of SPE-affected area. This work shows the importance of ionosphere-atmosphere interaction for SPE studies. By using both ionospheric and atmospheric observations, we have been able to cover for the most part the whole chain of SPE-triggered processes, from proton-induced ionisation to depletion of ozone.

Quantum Chemical Studies on Tropospheric Nucleation Mechanisms Involving Sulfuric Acid

Nucleation is the first step of the process by which gas molecules in the atmosphere condense to form liquid or solid particles. Despite the importance of atmospheric new-particle formation for both climate and health-related issues, little information exists on its precise molecular-level mechanisms. In this thesis, potential nucleation mechanisms involving sulfuric acid together with either water and ammonia or reactive biogenic molecules are studied using quantum chemical methods. Quantum chemistry calculations are based on the numerical solution of Schrödinger's equation for a system of atoms and electrons subject to various sets of approximations, the precise details of which give rise to a large number of model chemistries. A comparison of several different model chemistries indicates that the computational method must be chosen with care if accurate results for sulfuric acid - water - ammonia clusters are desired. Specifically, binding energies are incorrectly predicted by some popular density functionals, and vibrational anharmonicity must be accounted for if quantitatively reliable formation free energies are desired. The calculations reported in this thesis show that a combination of different high-level energy corrections and advanced thermochemical analysis can quantitatively replicate experimental results concerning the hydration of sulfuric acid. The role of ammonia in sulfuric acid - water nucleation was revealed by a series of calculations on molecular clusters of increasing size with respect to all three co-ordinates; sulfuric acid, water and ammonia. As indicated by experimental measurements, ammonia significantly assists the growth of clusters in the sulfuric acid - co-ordinate. The calculations presented in this thesis predict that in atmospheric conditions, this effect becomes important as the number of acid molecules increases from two to three. On the other hand, small molecular clusters are unlikely to contain more than one ammonia molecule per sulfuric acid. This implies that the average NH3:H2SO4 mole ratio of small molecular clusters in atmospheric conditions is likely to be between 1:3 and 1:1. Calculations on charged clusters confirm the experimental result that the HSO4- ion is much more strongly hydrated than neutral sulfuric acid. Preliminary calculations on HSO4- NH3 clusters indicate that ammonia is likely to play at most a minor role in ion-induced nucleation in the sulfuric acid - water system. Calculations of thermodynamic and kinetic parameters for the reaction of stabilized Criegee Intermediates with sulfuric acid demonstrate that quantum chemistry is a powerful tool for investigating chemically complicated nucleation mechanisms. The calculations indicate that if the biogenic Criegee Intermediates have sufficiently long lifetimes in atmospheric conditions, the studied reaction may be an important source of nucleation precursors.