Mathematical models on the impact of noise and dyadic molecular structures on the properties of a cardiac myocyte

In cardiac myocytes (heart muscle cells), coupling of electric signal known as the action potential to contraction of the heart depends crucially on calcium-induced calcium release (CICR) in a microdomain known as the dyad. During CICR, the peak number of free calcium ions (Ca) present in the dyad is small, typically estimated to be within range 1-100. Since the free Ca ions mediate CICR, noise in Ca signaling due to the small number of free calcium ions influences Excitation-Contraction (EC) coupling gain. Noise in Ca signaling is only one noise type influencing cardiac myocytes, e.g., ion channels playing a central role in action potential propagation are stochastic machines, each of which gates more or less randomly, which produces gating noise present in membrane currents. How various noise sources influence macroscopic properties of a myocyte, how noise is attenuated and taken advantage of are largely open questions. In this thesis, the impact of noise on CICR, EC coupling and, more generally, macroscopic properties of a cardiac myocyte is investigated at multiple levels of detail using mathematical models. Complementarily to the investigation of the impact of noise on CICR, computationally-efficient yet spatially-detailed models of CICR are developed. The results of this thesis show that (1) gating noise due to the high-activity mode of L-type calcium channels playing a major role in CICR may induce early after-depolarizations associated with polymorphic tachycardia, which is a frequent precursor to sudden cardiac death in heart failure patients; (2) an increased level of voltage noise typically increases action potential duration and it skews distribution of action potential durations toward long durations in cardiac myocytes; and that (3) while a small number of Ca ions mediate CICR, Excitation-Contraction coupling is robust against this noise source, partly due to the shape of ryanodine receptor protein structures present in the cardiac dyad.

Novel conditional mouse models for targeted kidney research : Emphasis on the analysis of the biological function of nephrin

Nephrin is a transmembrane protein belonging to the immunoglobulin superfamily and is expressed primarily in the podocytes, which are highly differentiated epithelial cells needed for primary urine formation in the kidney. Mutations leading to nephrin loss abrogate podocyte morphology, and result in massive protein loss into urine and consequent early death in humans carrying specific mutations in this gene. The disease phenotype is closely replicated in respective mouse models. The purpose of this thesis was to generate novel inducible mouse-lines, which allow targeted gene deletion in a time and tissue-specific manner. A proof of principle model for succesful gene therapy for this disease was generated, which allowed podocyte specific transgene replacement to rescue gene deficient mice from perinatal lethality. Furthermore, the phenotypic consequences of nephrin restoration in the kidney and nephrin deficiency in the testis, brain and pancreas in rescued mice were investigated. A novel podocyte-specific construct was achieved by using standard cloning techniques to provide an inducible tool for in vitro and in vivo gene targeting. Using modified constructs and microinjection procedures two novel transgenic mouse-lines were generated. First, a mouse-line with doxycycline inducible expression of Cre recombinase that allows podocyte-specific gene deletion was generated. Second, a mouse-line with doxycycline inducible expression of rat nephrin, which allows podocyte-specific nephrin over-expression was made. Furthermore, it was possible to rescue nephrin deficient mice from perinatal lethality by cross-breeding them with a mouse-line with inducible rat nephrin expression that restored the missing endogenous nephrin only in the kidney after doxycycline treatment. The rescued mice were smaller, infertile, showed genital malformations and developed distinct histological abnormalities in the kidney with an altered molecular composition of the podocytes. Histological changes were also found in the testis, cerebellum and pancreas. The expression of another molecule with limited tissue expression, densin, was localized to the plasma membranes of Sertoli cells in the testis by immunofluorescence staining. Densin may be an essential adherens junction protein between Sertoli cells and developing germ cells and these junctions share similar protein assembly with kidney podocytes. This single, binary conditional construct serves as a cost- and time-efficient tool to increase the understanding of podocyte-specific key proteins in health and disease. The results verified a tightly controlled inducible podocyte-specific transgene expression in vitro and in vivo as expected. These novel mouse-lines with doxycycline inducible Cre recombinase and with rat nephrin expression will be useful for conditional gene targeting of essential podocyte proteins and to study in detail their functions in the adult mice. This is important for future diagnostic and pharmacologic development platforms.

Animal model and molecular interactions of Cln5

Neuronal ceroid lipofuscinoses (NCLs) are a family of inherited pediatric neurodegenerative disorders, leading to retinal degeneration, death of selective neuronal populations and accumulation of autofluorscent ceroid-lipopigments. The clinical manifestations are generally similar in all forms. The Finnish variant late infantile neuronal ceroid lipofuscinosis (vLINCLFin) is a form of NCL, especially enriched in the Finnish population. The aim of this thesis was to analyse the brain pathology of vLINCLFin utilising the novel Cln5-/- mouse model. Gene expression profiling of the brains of already symptomatic Cln5-/- mice revealed that inflammation, neurodegeneration and defects in myelinization are the major characteristics of the later stages of the disease. Histological characterization of the brain pathology confirmed that the thalamocortical system is affected in Cln5-/- mice, similarly to the other NCL mouse models. However, whereas the brain pathology in all other analyzed NCL mice initiate in the thalamus and spread only months later to the cortex, we observed that the sequence of events is uniquely reversed in Cln5-/- mice; beginning in the cortex and spreading to the thalamus only months later. We could also show that even though neurodegeneration is inititated in the cortex, reactive gliosis and loss of myelin are evident in specific nuclei of the thalamus already in the 1 month old brain. To obtain a deeper insight into the disturbed metabolic pathways, we performed gene expression profiling of presymptomatic mouse brains. We validated these findings with immunohistological analyses, and could show that cytoskeleton and myelin were affected in Cln5-/- mice. Comparison of gene expression profiling results of Cln5-/- and Cln1-/- mice, further highlighted that these two NCL models share a common defective pathway, leading to disturbances in the neuronal growth cone and cytoskeleton. Encouraged by the evidence of this defected pathway, we analyzed the molecular interactions of NCL-proteins and observed that Cln5 and Cln1/Ppt1 proteins interact with each other. Furthermore, we demonstrated that Cln5 and Cln1/Ppt1 share an interaction partner, the F1-ATP synthase, potentially linking both vLINCLFIN and INCL diseases to disturbed lipid metabolism. In addition, Cln5 was shown to interact with other NCL proteins; Cln2, Cln3, Cln6 and Cln8, implicating a central role for Cln5 in the NCL pathophysiology. This study is the first to describe the brain pathology and gene expression changes in the Cln5-/- mouse. Together the findings presented in this thesis represent novel information of the disease processes and the molecular mechanisms behind vLINCLFin and have highlighted that vLINCLFin forms a very important model to analyze the pathophysiology of NCL diseases.

Molecular and Cellular Mechanisms Behind Juvenile Neuronal Ceroid Lipofuscinosis (JNCL, Batten disease)

Neurodegenerative disorders are chronic, progressive, and often fatal disorders of the nervous system caused by dysfunction, and ultimately, death of neuronal cells. The underlying mechanisms of neurodegeneration are poorly understood, and monogenic disorders can be utilised as disease models to elucidate the pathogenesis. Juvenile neuronal ceroid-lipofuscinosis (JNCL, Batten disease) is a recessively inherited lysosomal storage disorder with progressive neurodegeneration and accumulation of autofluorescent storage material in most tissues. It is caused by mutations in the CLN3 gene; however, the exact function of the corresponding CLN3 protein, as well as the molecular mechanisms of JNCL pathogenesis have remained elusive. JNCL disease exclusively affects the central nervous system leaving other organs unaffected, and therefore it is of a particular importance to conduct studies in brain tissue and neuronal cells. The aim of this thesis project was to elucidate the molecular and cell biological mechanisms underlying JNCL. This was the first study to describe the endogenous Cln3 protein, and it was shown that Cln3 localised to neuronal cells in the mouse brain. At a subcellular level, endogenous Cln3 was localised to the presynaptic terminals and to the synaptosome compartment, but not to the synaptic vesicles. Studies with the CLN3-deficient cells demonstrated an impaired endocytic membrane trafficking, and established an interconnection between CLN3, microtubulus-binding Hook1 and Rab proteins. This novel data was not only important in characterising the roles of CLN3 in cells, but also provided significant information delineating the versatile role of the Rab proteins. To identify affected cellular pathways in JNCL, global gene expression profiling of the knock-out mouse Cln3-/- neurons was performed and systematically analysed; this revealed a slight dysfunction of the mitochondria, cytoskeletal abnormality in the microtubule plus-end, and an impaired recovery from depolarizing stimulus when specific N-type Ca2+ channels were inhibited, thus leading to a prolonged time of higher intracellular calcium. All these defective pathways are interrelated, and may together be sufficient to initiate the neurodegenerative process. Results of this thesis also suggest that in neuronal cells, CLN3 most likely functions at endocytic vesicles at the presynaptic terminal, potentially involved in the regulation of the calcium-mediated synaptic transmission.

Molecular profiling of malignant pleural mesothelioma and pulmonary fibrosis

Malignant mesothelioma (MM) is a rare, usually incurable, disease mainly caused by former exposure to asbestos. Even though MM has a strong etiological link, genetic factors may play a role, since not all cases can be linked to former asbestos exposure. This thesis focuses on lung diseases, mainly malignant mesothelioma (MM), and idiopathic pulmonary fibrosis (IPF), which resembles asbestosis. The specific asbestos-related pathways associated with malignant as well as non-malignant lung diseases, still need to be clarified. Since most patients diagnosed with MM or asbestosis/fibrosis have a dismal prognosis and few therapeutic options are available, early diagnosis and better understanding of the disease pathogenesis are of the utmost importance. The first objective of this thesis was to identify asbestos specific differentially expressed genes. This was approached by using high-resolution gene expression arrays, and three different human lung cell lines, as well as with three different bioinformatics approaches. Since the first study aimed to elucidate potential early changes, the second study was used to screen DNA copy number changes in MM tumour samples. This was performed using genome wide microarrays for identification of DNA copy number changes characterstic for MM. Study III focused on the role of gremlin in the regulation of bone morphogenetic protein (BMPs) in IPF. Further studies were conducted in asbestos-exposed cell cultures as well as in an asbestos-induced mouse model. Furthermore, GATA-6 was studied in MM and metastatic pleural adenocarcinoma. The GATA transcription factors are important during embryonic development, but their role in cancer is still unclear. GATA-6 is a co-factor/target of thyroid transcription factor 1 (TTF-1), which is used in differential diagnostics of pleural MM and adenocarcinoma. Bioinformatics probed the genes and biological processes ordered in terms of significance, clusters, and highly enriched chromosomal regions. The study revealed several already identified targets, produced new ideas about genes which are central for asbestos exposure, as well as provided supplementary data for researchers to check their own novel findings or ideas. The analysis revealed DNA copy number changes characteristic for MM tumors. The most common regions of loss were detected in 1p, 3p, 6q, 9p, 13, 14, and 22, and gains at 17q. The histological features in asbestosis and IPF are very similar, wherefore IPF can be studied in asbestos models. The BMP antagonist gremlin was up-regulated by asbestos exposure in human epithelial cell lines, which was also observed in Study I. The transforming growth factor (TGF) -β and BMP expression and signaling activities were measured from murine and human fibrotic lungs. BMP-7 signaling was down-regulated in response to up-regulation of gremlin, and restoration of BMP-7 signaling prevented progression of fibrosis in mice. Therefore, the study suggests that the restoration of BMP-7 signaling in fibrotic lung could potentially aid in the treatment of IPF patients. Study IV revealed that GATA-6 was strongly expressed in the majority of the MM cases, and correlated statistically significant with longer survival in subgroups of MM.

New molecular strategies for prevention of fibroproliferative vascular disease : An experimental approach to restenosis

Vascular intimal hyperplasia is a major complication following angioplasty. The hallmark feature of this disorder is accumulation of dedifferentiated smooth muscle cells (SMCs) to the luminal side of the injured artery, cellular proliferation, migration, and synthesis of extracellular matrix. This finally results in intimal hyperplasia, which is currently considered an untreatable condition. According to current knowledge, a major part of neointimal cells derive from circulating precursor cells. This has outdated the traditional in vitro cell culture methods of studying neointimal cell migration and proliferation using cultured medial SMCs. Somatostatin and some of its analogs with different selectivity for the five somatostatin receptors (sst1 through sst5) have been shown to have vasculoprotective properties in animal studies. However, clinical trials using analogs selective for sst2/sst3/sst5 to prevent restenosis after percutaneous transluminal coronary angioplasty (PTCA) have failed to show any major benefits. Sirolimus is a cell cycle inhibitor that has been suggested to act synergistically with the protein-tyrosine kinase inhibitor imatinib to inhibit intimal hyperplasia in rat already at well-tolerated submaximal oral doses. The mechanisms behind this synergy and its long-term efficacy are not known. The aim of this study was to set up an ex vivo vascular explant culture model to measure neointimal cell activity without excluding the participation of circulating progenitor cells. Furthermore, two novel potential vasculoprotective treatment strategies were evaluated in detail in rat models of intimal hyperplasia and in the ex vivo explant model: sst1/sst4-selective somatostatin receptor analogs and combination treatment with sirolimus and imatinib. This study shows how whole vessel explants can be used to study the kinetics of neointimal cells and their progenitors, and to evaluate the anti-migratory and anti-proliferative properties of potential vasculoprotective compounds. It also shows how the influx of neointimal progenitor cells occurs already during the first days after vascular injury, how the contribution of cell migration is more important in the injury response than cell proliferation, and how the adventitia actively contribute in vascular repair. The vasculoprotective effect of somatostatin is mediated preferentially through sst4, and through inhibition of cell migration rather than of proliferation, which may explain why sst2/sst3/sst5-selective analogs have failed in clinical trials. Furthermore, a brief early oral treatment with the combination of sirolimus and imatinib at submaximal doses results in long-term synergistic suppression of intimal hyperplasia. The synergy is a result of inhibition of post-operative thrombocytosis and leukocytosis, inhibition of neointimal cell migration to the injury-site, and maintenance of cell integrity by inhibition of apoptosis and SMC dedifferentiation. In conclusion, the influx of progenitor cells already during the first days after injury and the high neointimal cell migratory activity underlines the importance of early therapeutic intervention with anti-migratory compounds to prevent neointimal hyperplasia. Sst4-selective analogs and the combination therapy with sirolimus and imatinib represent potential targets for the development of such vasculoprotective therapies.

Water-ethanol mixtures by molecular dynamics and x-ray Compton scattering

Water-ethanol mixtures are commonly used in industry and house holds. However, quite surprisingly their molecular-level structure is still not completely understood. In particular, there is evidence that the local intermolecular geometries depend significantly on the concentration. The aim of this study was to gain information on the molecular-level structures of water-ethanol mixtures by two computational methods. The methods are classical molecular dynamics (MD), where the movement of molecules can be studied, and x-ray Compton scattering, in which the scattering cross section is sensitive to the electron momentum density. Firstly, the water-ethanol mixtures were studied with MD simulations, with the mixture concentration ranging from 0 to 100%. For the simulations well-established force fields were used for the water and ethanol molecules (TIP4P and OPLS-AA, respectively). Moreover, two models were used for ethanol, rigid and non-rigid. In the rigid model the intramolecular bond lengths are fixed, whereas in the non-rigid model the lengths are determined by harmonic potentials. Secondly, mixtures with three different concentrations employing both ethanol models were studied by calculating the experimentally observable x-ray quantity, the Compton profile. In the MD simulations a slight underestimation in the density was observed as compared to experiment. Furthermore, a positive excess of hydrogen bonding with water molecules and a negative one with ethanol was quantified. Also, the mixture was found more structured when the ethanol concentration was higher. Negligible differences in the results were found between the two ethanol models. In contrast, in the Compton scattering results a notable difference between the ethanol models was observed. For the rigid model the Compton profiles were similar for all the concentrations, but for the non-rigid model they were distinct. This leads to two possibilities of how the mixing occurs. Either the mixing is similar in all concentrations (as suggested by the rigid model) or the mixing changes for different concentrations (as suggested by the non-rigid model). Either way, this study shows that the choice of the force field is essential in the microscopic structure formation in the MD simulations. When the sources of uncertainty in the calculated Compton profiles were analyzed, it was found that more statistics needs to be collected to reduce the statistical uncertainty in the final results. The obtained Compton scattering results can be considered somewhat preliminary, but clearly indicative of the behaviour of the water-ethanol mixtures when the force field is modified. The next step is to collect more statistics and compare the results with experimental data to decide which ethanol model describes the mixture better. This way, valuable information on the microscopic structure of water-ethanol mixtures can be found. In addition, information on the force fields in the MD simulations and on the ability of the MD simulations to reproduce the microscopic structure of binary liquids is obtained.

Discovering hidden structures in molecular data using a Bayesian partition model approach

Advancements in the analysis techniques have led to a rapid accumulation of biological data in databases. Such data often are in the form of sequences of observations, examples including DNA sequences and amino acid sequences of proteins. The scale and quality of the data give promises of answering various biologically relevant questions in more detail than what has been possible before. For example, one may wish to identify areas in an amino acid sequence, which are important for the function of the corresponding protein, or investigate how characteristics on the level of DNA sequence affect the adaptation of a bacterial species to its environment. Many of the interesting questions are intimately associated with the understanding of the evolutionary relationships among the items under consideration. The aim of this work is to develop novel statistical models and computational techniques to meet with the challenge of deriving meaning from the increasing amounts of data. Our main concern is on modeling the evolutionary relationships based on the observed molecular data. We operate within a Bayesian statistical framework, which allows a probabilistic quantification of the uncertainties related to a particular solution. As the basis of our modeling approach we utilize a partition model, which is used to describe the structure of data by appropriately dividing the data items into clusters of related items. Generalizations and modifications of the partition model are developed and applied to various problems. Large-scale data sets provide also a computational challenge. The models used to describe the data must be realistic enough to capture the essential features of the current modeling task but, at the same time, simple enough to make it possible to carry out the inference in practice. The partition model fulfills these two requirements. The problem-specific features can be taken into account by modifying the prior probability distributions of the model parameters. The computational efficiency stems from the ability to integrate out the parameters of the partition model analytically, which enables the use of efficient stochastic search algorithms.

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

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

Mechanisms and molecular regulation of mammalian tooth replacement

In most non-mammalian vertebrates, such as fish and reptiles, teeth are replaced continuously. However, tooth replacement in most mammals, including human, takes place only once and further renewal is apparently inhibited. It is not known how tooth replacement is genetically regulated, and little is known on the physiological mechanism and evolutionary reduction of tooth replacement in mammals. In this study I have attempted to address these questions. In a rare human condition cleidocranial dysplasia, caused by a mutation in a Runt domain transcription factor Runx2, tooth replacement is continued. Runx2 mutant mice were used to investigate the molecular mechanisms of Runx2 function. Microarray analysis from dissected embryonic day 14 Runx2 mutant and wild type dental mesenchymes revealed many downstream targets of Runx2, which were validated using in situ hybridization and tissue culture methods. Wnt signaling inhibitor Dkk1 was identified as a candidate target, and in tissue culture conditions it was shown that Dkk1 is induced by FGF4 and this induction is Runx2 dependent. These experiments demonstrated a connection between Runx2, FGF and Wnt signaling in tooth development and possibly also in tooth replacement. The role of Wnt signaling in tooth replacement was further investigated by using a transgenic mouse model where Wnt signaling mediator β-catenin is continuously stabilized in dental epithelium. This stabilization led to activated Wnt signaling and to the formation of multiple enamel knots. In vitro and transplantation experiments were performed to examine the process of extra tooth formation. We showed that new teeth were continuously generated and that new teeth form from pre-existing teeth. A morphodynamic activator-inhibitor model was used to simulate enamel knot formation. By increasing the intrinsic production rate of the activator (β-catenin), the multiple enamel knot phenotype was reproduced by computer simulations. It was thus concluded that β-catenin acts as an upstream activator of enamel knots, closely linking Wnt signaling to the regulation of tooth renewal. As mice do not normally replace teeth, we used other model animals to investigate the physiological and genetic mechanisms of tooth replacement. Sorex araneus, the common shrew was earlier reported to have non-functional tooth replacement in all antemolar tooth positions. We showed by histological and gene expression studies that there is tooth replacement only in one position, the premolar 4 and that the deciduous tooth is diminished in size and disappears during embryogenesis without becoming functional. The growth rates of deciduous and permanent premolar 4 were measured and it was shown by competence inference that the early initiation of the replacement tooth in relation to the developmental stage of the deciduous tooth led to the inhibition of deciduous tooth morphogenesis. It was concluded that the evolutionary loss of deciduous teeth may involve the early activation of replacement teeth, which in turn suppress their predecessors. Mustela putorius furo, the ferret, has a dentition that resembles that of the human as ferrets have teeth that belong to all four tooth families, and all the antemolar teeth are replaced once. To investigate the replacement mechanism, histological serial sections from different embryonic stages were analyzed. It was noticed that tooth replacement is a process which involves the growth and detachment of the dental lamina from the lingual cervical loop of the deciduous tooth. Detachment of the deciduous tooth leads to a free successional dental lamina, which grows deeper into the mesenchyme, and later buds the replacement tooth. A careful 3D analysis of serial histological sections was performed and it was shown that replacement teeth are initiated from the successional dental lamina and not from the epithelium of the deciduous tooth. The molecular regulation of tooth replacement was studied and it was shown by examination of expression patterns of candidate regulatory genes that BMP/Wnt inhibitor Sostdc1 was strongly expressed in the buccal aspect of the dental lamina, and in the intersection between the detaching deciduous tooth and the successional dental lamina, suggesting a role for Sostdc1 in the process of detachment. Shh was expressed in the enamel knot and in the inner enamel epithelium in both generations of teeth supporting the view that the morphogenesis of both generations of teeth is regulated by similar mechanisms. In summary, histological and molecular studies on different model animals and transgenic mouse models were used to investigate tooth replacement. This thesis work has significantly contributed to the knowledge on the physiological mechanisms and molecular regulation of tooth replacement and its evolutionary suppression in mammals.

Mutations as molecular tools : The metabolic-rate dependent molecular clock and DNA barcoding of allied species

Mutation and recombination are the fundamental processes leading to genetic variation in natural populations. This variation forms the raw material for evolution through natural selection and drift. Therefore, studying mutation rates may reveal information about evolutionary histories as well as phylogenetic interrelationships of organisms. In this thesis two molecular tools, DNA barcoding and the molecular clock were examined. In the first part, the efficiency of mutations to delineate closely related species was tested and the implications for conservation practices were assessed. The second part investigated the proposition that a constant mutation rate exists within invertebrates, in form of a metabolic-rate dependent molecular clock, which can be applied to accurately date speciation events. DNA barcoding aspires to be an efficient technique to not only distinguish between species but also reveal population-level variation solely relying on mutations found on a short stretch of a single gene. In this thesis barcoding was applied to discriminate between Hylochares populations from Russian Karelia and new Hylochares findings from the greater Helsinki region in Finland. Although barcoding failed to delineate the two reproductively isolated groups, their distinct morphological features and differing life-history traits led to their classification as two closely related, although separate species. The lack of genetic differentiation appears to be due to a recent divergence event not yet reflected in the beetles molecular make-up. Thus, the Russian Hylochares was described as a new species. The Finnish species, previously considered as locally extinct, was recognized as endangered. Even if, due to their identical genetic make-up, the populations had been regarded as conspecific, conservation strategies based on prior knowledge from Russia would not have guaranteed the survival of the Finnish beetle. Therefore, new conservation actions based on detailed studies of the biology and life-history of the Finnish Hylochares were conducted to protect this endemic rarity in Finland. The idea behind the strict molecular clock is that mutation rates are constant over evolutionary time and may thus be used to infer species divergence dates. However, one of the most recent theories argues that a strict clock does not tick per unit of time but that it has a constant substitution rate per unit of mass-specific metabolic energy. Therefore, according to this hypothesis, molecular clocks have to be recalibrated taking body size and temperature into account. This thesis tested the temperature effect on mutation rates in equally sized invertebrates. For the first dataset (family Eucnemidae, Coleoptera) the phylogenetic interrelationships and evolutionary history of the genus Arrhipis had to be inferred before the influence of temperature on substitution rates could be studied. Further, a second, larger invertebrate dataset (family Syrphidae, Diptera) was employed. Several methodological approaches, a number of genes and multiple molecular clock models revealed that there was no consistent relationship between temperature and mutation rate for the taxa under study. Thus, the body size effect, observed in vertebrates but controversial for invertebrates, rather than temperature may be the underlying driving force behind the metabolic-rate dependent molecular clock. Therefore, the metabolic-rate dependent molecular clock does not hold for the here studied invertebrate groups. This thesis emphasizes that molecular techniques relying on mutation rates have to be applied with caution. Whereas they may work satisfactorily under certain conditions for specific taxa, they may fail for others. The molecular clock as well as DNA barcoding should incorporate all the information and data available to obtain comprehensive estimations of the existing biodiversity and its evolutionary history.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Molecular mechanisms of hypertension-induced heart failure : Experimental studies with special emphasis on local renin-angiotensin system, cardiac metabolism and levosimendan

Hypertension is a major risk factor for stroke, ischaemic heart disease, and the development of heart failure. Hypertension-induced heart failure is usually preceded by the development of left ventricular hypertrophy (LVH), which represents an adaptive and compensatory response to the increased cardiac workload. Biomechanical stress and neurohumoral activation are the most important triggers of pathologic hypertrophy and the transition of cardiac hypertrophy to heart failure. Non-clinical and clinical studies have also revealed derangements of energy metabolism in hypertensive heart failure. The goal of this study was to investigate in experimental models the molecular mechanisms and signalling pathways involved in hypertension-induced heart failure with special emphasis on local renin-angiotensin-aldosterone system (RAAS), cardiac metabolism, and calcium sensitizers, a novel class of inotropic agents used currently in the treatment of acute decompensated heart failure. Two different animal models of hypertensive heart failure were used in the present study, i.e. hypertensive and salt-sensitive Dahl/Rapp rats on a high salt diet (a salt-sensitive model of hypertensive heart failure) and double transgenic rats (dTGR) harboring human renin and human angiotensinogen genes (a transgenic model of hypertensive heart failure with increased local RAAS activity). The influence of angiotensin II (Ang II) on cardiac substrate utilization and cardiac metabolomic profile was investigated by using gas chromatography coupled to time-of-flight mass spectrometry to detect 247 intermediary metabolites. It was found that Ang II could alter cardiac metabolomics both in normotensive and hypertensive rats in an Ang II receptor type 1 (AT1)-dependent manner. A distinct substrate use from fatty acid oxidation towards glycolysis was found in dTGR. Altered cardiac substrate utilization in dTGR was associated with mitochondrial dysfunction. Cardiac expression of the redox-sensitive metabolic sensor sirtuin1 (SIRT1) was increased in dTGR. Resveratrol supplementation prevented cardiovascular mortality and ameliorated Ang II-induced cardiac remodeling in dTGR via blood pressure-dependent pathways and mechanisms linked to increased mitochondrial biogenesis. Resveratrol dose-dependently increased SIRT1 activity in vitro. Oral levosimendan treatment was also found to improve survival and systolic function in dTGR via blood pressure-independent mechanisms, and ameliorate Ang II-induced coronary and cardiomyocyte damage. Finally, using Dahl/Rapp rats it was demonstrated that oral levosimendan as well as the AT1 receptor antagonist valsartan improved survival and prevented cardiac remodeling. The beneficial effects of levosimendan were associated with improved diastolic function without significantly improved systolic changes. These positive effects were potentiated when the drug combination was administered. In conclusion, the present study points to an important role for local RAAS in the pathophysiology of hypertension-induced heart failure as well as its involvement as a regulator of cardiac substrate utilization and mitochondrial function. Our findings suggest a therapeutic role for natural polyphenol resveratrol and calcium sensitizer, levosimendan, and the novel drug combination of valsartan and levosimendan, in prevention of hypertension-induced heart failure. The present study also provides a better understanding of the pathophysiology of hypertension-induced heart failure, and may help identify potential targets for novel therapeutic interventions.

Association between molecular variation in the phosphoglucose isomerase (Pgi) locus and migration propensity in the Glanville fritillary butterfly

Habitat fragmentation produces patches of suitable habitat surrounded by unfavourable matrix habitat. A species may persist in such a fragmented landscape in an equilibrium between the extinctions and recolonizations of local populations, thus forming a metapopulation. Migration between local populations is necessary for the long-term persistence of a metapopulation. The Glanville fritillary butterfly (Melitaea cinxia) forms a metapopulation in the Åland islands in Finland. There is migration between the populations, the extent of which is affected by several environmental factors and variation in the phenotype of individual butterflies. Different allelic forms of the glycolytic enzyme phosphoglucose isomerase (Pgi) has been identified as a possible genetic factor influencing flight performance and migration rate in this species. The frequency of a certain Pgi allele, Pgi-f, follows the same pattern in relation to population age and connectivity as migration propensity. Furthermore, variation in flight metabolic performance, which is likely to affect migration propensity, has been linked to genetic variation in Pgi or a closely linked locus. The aim of this study was to investigate the association between Pgi genotype and the migration propensity in the Glanville fritillary both at the individual and population levels using a statistical modelling approach. A mark-release-recapture (MRR) study was conducted in a habitat patch network of M. cinxia in Åland to collect data on the movements of individual butterflies. Larval samples from the study area were also collected for population level examinations. Each butterfly and larva was genotyped at the Pgi locus. The MRR data was parameterised with two mathematical models of migration: the Virtual Migration Model (VM) and the spatially explicit diffusion model. VM model predicted and observed numbers of emigrants from populations with high and low frequencies of Pgi-f were compared. Posterior predictive data sets were simulated based on the parameters of the diffusion model. Lack-of-fit of observed values to the model predicted values of several descriptors of movements were detected, and the effect of Pgi genotype on the deviations was assessed by randomizations including the genotype information. This study revealed a possible difference in the effect of Pgi genotype on migration propensity between the two sexes in the Glanville fritillary. The females with and males without the Pgi-f allele moved more between habitat patches, which is probably related to differences in the function of flight in the two sexes. Females may use their high flight capacity to migrate between habitat patches to find suitable oviposition sites, whereas males may use it to acquire mates by keeping a territory and fighting off other intruding males, possibly causing them to emigrate. The results were consistent across different movement descriptors and at the individual and population levels. The effect of Pgi is likely to be dependent on the structure of the landscape and the prevailing environmental conditions.