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.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Molecular Dynamics Simulations of Homogeneous Nucleation

Nucleation is the first step in a phase transition where small nuclei of the new phase start appearing in the metastable old phase, such as the appearance of small liquid clusters in a supersaturated vapor. Nucleation is important in various industrial and natural processes, including atmospheric new particle formation: between 20 % to 80 % of atmospheric particle concentration is due to nucleation. These atmospheric aerosol particles have a significant effect both on climate and human health. Different simulation methods are often applied when studying things that are difficult or even impossible to measure, or when trying to distinguish between the merits of various theoretical approaches. Such simulation methods include, among others, molecular dynamics and Monte Carlo simulations. In this work molecular dynamics simulations of the homogeneous nucleation of Lennard-Jones argon have been performed. Homogeneous means that the nucleation does not occur on a pre-existing surface. The simulations include runs where the starting configuration is a supersaturated vapor and the nucleation event is observed during the simulation (direct simulations), as well as simulations of a cluster in equilibrium with a surrounding vapor (indirect simulations). The latter type are a necessity when the conditions prevent the occurrence of a nucleation event in a reasonable timeframe in the direct simulations. The effect of various temperature control schemes on the nucleation rate (the rate of appearance of clusters that are equally able to grow to macroscopic sizes and to evaporate) was studied and found to be relatively small. The method to extract the nucleation rate was also found to be of minor importance. The cluster sizes from direct and indirect simulations were used in conjunction with the nucleation theorem to calculate formation free energies for the clusters in the indirect simulations. The results agreed with density functional theory, but were higher than values from Monte Carlo simulations. The formation energies were also used to calculate surface tension for the clusters. The sizes of the clusters in the direct and indirect simulations were compared, showing that the direct simulation clusters have more atoms between the liquid-like core of the cluster and the surrounding vapor. Finally, the performance of various nucleation theories in predicting simulated nucleation rates was investigated, and the results among other things highlighted once again the inadequacy of the classical nucleation theory that is commonly employed in nucleation studies.

Molecular characterization of sediment bacterial communities affected by fish farming

Fish farming introduces nutrients, microbes and a wide variety of chemicals such as heavy metals, antifoulants and antibiotics to the surrounding environment. Introduction of antibiotics has been linked with the increased incidence of antibiotic resistant pathogenic bacteria in the farm vicinities. In this thesis molecular methods such as quantitative PCR and DNA sequencing were applied to analyze bacterial communities in sediments from fish farms and pristine locations. Altogether four farms and four pristine sites were sampled in the Baltic Sea. Two farm and two pristine locations were sampled over a surveillance period of four years. Furthermore, a new methodology was developed as a part of the study that permits amplifying single microbial genomes and capturing them according to any genetic traits, including antibiotic resistance genes. The study revealed that several resistance genes for tetracycline were found at the sediment underneath the aquaculture farms. The copy number of these genes remained elevated even at a farm that had not used any antibiotics since year 2000, six years before this study started. Similarly, an increase in the amount of mercury resistance gene merA was observed at the aquaculture sediment. The persistence of the resistance genes in absence of any selection pressure from antibiotics or heavy metals suggests that the genes may be introduced to the sediment by the farming process. This is also supported by the diversity pattern of the merA gene between farm and pristine sediments. The bacterial community-level changes in response to fish farming were very complex and no single phylogenetic groups were found that would be typical to fish farm sediments. However, the community structures had some correlation with the exposure to fish farming. Our studies suggest that the established approaches to deal with antibiotic resistance at the aquaculture, such as antibiotic cycling, are fundamentally flawed because they cannot prevent the introduction of the resistance genes and resistant bacteria to the farm area by the farming process. Further studies are required to study the entire fish farming process to identify the sources of the resistance genes and the resistant bacteria. The results also suggest that in order to prevent major microbiological changes in the surrounding aquatic environment, the farms should not be founded in shallow water where currents do not transport sedimenting matter from the farms. Finally, the technique to amplify and select microbial genomes will potentially have a considerable impact in microbial ecology and genomics.

Myotonic dystrophy type 2 (DM2) : Diagnostic methods and molecular pathology

Myotonic dystrophies type 1 (DM1) and type 2 (DM2) are the most common forms of muscular dystrophy affecting adults. They are autosomal dominant diseases caused by microsatellite tri- or tetranucleotide repeat expansion mutations in transcribed but not translated gene regions. The mutant RNA accumulates in nuclei disturbing the expression of several genes. The more recently identified DM2 disease is less well known, yet more than 300 patients have been confirmed in Finland thus far, and the true number is believed to be much higher. DM1 and DM2 share some features in general clinical presentation and molecular pathology, yet they show distinctive differences, including disease severity and differential muscle and fiber type involvement. However, the molecular differences underlying DM1 and DM2 muscle pathology are not well understood. Although the primary tissue affected is muscle, both DMs show a multisystemic phenotype due to wide expression of the mutation-carrying genes. DM2 is particularly intriguing, as it shows an incredibly wide spectrum of clinical manifestations. For this reason, it constitutes a real diagnostic challenge. The core symptoms in DM2 include proximal muscle weakness, muscle pain, myotonia, cataracts, cardiac conduction defects and endocrinological disturbations; however, none of these is mandatory for the disease. Myalgic pains may be the most disabling symptom for decades, sometimes leading to incapacity for work. In addition, DM2 may cause major socio-economical consequences for the patient, if not diagnosed, due to misunderstanding and false stigmatization. In this thesis work, we have (I) improved DM2 differential diagnostics based on muscle biopsy, and (II) described abnormalities in mRNA and protein expression in DM1 and DM2 patient skeletal muscles, showing partial differences between the two diseases, which may contribute to muscle pathology in these diseases. This is the first description of histopathological differences between DM1 and DM2, which can be used in differential diagnostics. Two novel high-resolution applications of in situ -hybridization have been described, which can be used for direct visualization of the DM2 mutation in muscle biopsy sections, or mutation size determination on extended DNA-fibers. By measuring protein and mRNA expression in the samples, differential changes in expression patterns affecting contractile proteins, other structural proteins and calcium handling proteins in DM2 compared to DM1 were found. The dysregulation at mRNA level was caused by altered transciption and abnormal splicing. The findings reported here indicate that the extent of aberrant splicing is higher in DM2 compared to DM1. In addition, the described abnormalities to some extent correlate to the differences in fiber type involvement in the two disorders.

Real-Time PCR : A Molecular Approach to Investigate the Role of Intestinal Microbiota in the Pathophysiology of Irritable Bowel Syndrome

Irritable bowel syndrome (IBS) is a common multifactorial functional intestinal disorder, the pathogenesis of which is not completely understood. Increasing scientific evidence suggests that microbes are involved in the onset and maintenance of IBS symptoms. The microbiota of the human gastrointestinal (GI) tract constitutes a massive and complex ecosystem consisting mainly of obligate anaerobic microorganisms making the use of culture-based methods demanding and prone to misinterpretation. To overcome these drawbacks, an extensive panel of species- and group-specific assays for an accurate quantification of bacteria from fecal samples with real-time PCR was developed, optimized, and validated. As a result, the target bacteria were detectable at a minimum concentration range of approximately 10 000 bacterial genomes per gram of fecal sample, which corresponds to the sensitivity to detect 0.000001% subpopulations of the total fecal microbiota. The real-time PCR panel covering both commensal and pathogenic microorganisms was assessed to compare the intestinal microbiota of patients suffering from IBS with a healthy control group devoid of GI symptoms. Both the IBS and control groups showed considerable individual variation in gut microbiota composition. Sorting of the IBS patients according to the symptom subtypes (diarrhea, constipation, and alternating predominant type) revealed that lower amounts of Lactobacillus spp. were present in the samples of diarrhea predominant IBS patients, whereas constipation predominant IBS patients carried increased amounts of Veillonella spp. In the screening of intestinal pathogens, 17% of IBS samples tested positive for Staphylococcus aureus, whereas no positive cases were discovered among healthy controls. Furthermore, the methodology was applied to monitor the effects of a multispecies probiotic supplementation on GI microbiota of IBS sufferers. In the placebo-controlled double-blind probiotic intervention trial of IBS patients, each supplemented probiotic strain was detected in fecal samples. Intestinal microbiota remained stable during the trial, except for Bifidobacterium spp., which increased in the placebo group and decreased in the probiotic group. The combination of assays developed and applied in this thesis has an overall coverage of 300-400 known bacterial species, along with the number of yet unknown phylotypes. Hence, it provides good means for studying the intestinal microbiota, irrespective of the intestinal condition and health status. In particular, it allows screening and identification of microbes putatively associated with IBS. The alterations in the gut microbiota discovered here support the hypothesis that microbes are likely to contribute to the pathophysiology of IBS. The central question is whether the microbiota changes described represent the cause for, rather than the effect of, disturbed gut physiology. Therefore, more studies are needed to determine the role and importance of individual microbial species or groups in IBS. In addition, it is essential that the microbial alterations observed in this study will be confirmed using a larger set of IBS samples of different subtypes, preferably from various geographical locations.

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.

Molecular background of three lethal fetal syndromes

This study identified the molecular defects underlying three lethal fetal syndromes. Lethal Congenital Contracture Syndrome 1 (LCCS1, MIM 253310) and Lethal Arthrogryposis with Anterior Horn Cell Disease (LAAHD, MIM 611890) are fetal motor neuron diseases. They affect the nerve cells that control voluntary muscle movement, and eventually result in severe atrophy of spinal cord motor neurons and fetal immobility. Both LCCS1 and LAAHD are caused by mutations in the GLE1 gene, which encodes for a multifunctional protein involved in posttranscriptional mRNA processing. LCCS2 and LCCS3, two syndromes that are clinically similar to LCCS1, are caused by defective proteins involved in the synthesis of inositol hexakisphosphate (IP6), an essential cofactor of GLE1. This suggests a common mechanism behind these fetal motor neuron diseases, and along with accumulating evidence from genetic studies of more late-onset motor neuron diseases such as Spinal muscular atrophy (SMA) and Amyotrophic lateral sclerosis (ALS), implicates mRNA processing as a common mechanism in motor neuron disease pathogenesis. We also studied gle1-/- zebrafish in order to investigate whether they would be a good model for studying the pathogenesis of LCCS1 and LAAHD. Mutant zebrafish exhibit cell death in their central nervous system at two days post fertilization, and the distribution of mRNA within the cells of mutant zebrafish differs from controls, encouraging further studies. The third lethal fetal syndrome is described in this study for the first time. Cocoon syndrome (MIM 613630) was discovered in a Finnish family with two affected individuals. Its hallmarks are the encasement of the limbs under the skin, and severe craniofacial abnormalities, including the lack of skull bones. We showed that Cocoon syndrome is caused by a mutation in the gene encoding the conserved helix-loop-helix ubiquitous kinase CHUK, also known as IκB kinase α (IKKα). The mutation results in the complete lack of CHUK protein expression. CHUK is a subunit of the IκB kinase enzyme that inhibits NF-κB transcription factors, but in addition, it has an essential, independent role in controlling keratinocyte differentiation, as well as informing morphogenetic events such as limb and skeletal patterning. CHUK also acts as a tumor suppressor, and is frequently inactivated in cancer. This study has brought significant new information about the molecular background of these three lethal fetal syndromes, as well as provided knowledge about the prerequisites of normal human development.

Molecular and Clinical Characteristics of Pituitary Adenoma Predisposition (PAP)

Pituitary adenomas are common benign neoplasms. Although most of them are sporadic, a minority occurs in familial settings. Heterozygous germline mutations in the aryl hydrocarbon receptor interacting protein (AIP) gene were found to underlie familial pituitary adenomas, a condition designated as pituitary adenoma predisposition (PAP). PAP confers incomplete penetrance of mostly growth hormone (GH) secreting adenomas in young patients, who often lack a family history of pituitary adenomas. This thesis work aimed to clarify the molecular and clinical characteristics of PAP. Applying the multiplex ligation-dependent probe amplification assay (MLPA), we found large genomic AIP deletions to account for a subset of PAP. Therefore, MLPA could be considered in PAP suspected patients with no AIP mutations found with conventional sequencing. We generated an Aip mouse model to examine pituitary tumorigenesis in vivo. The heterozygous Aip mutation conferred complete penetrance of pituitary adenomas that were mostly GH-secreting, rendering the phenotype of the Aip mouse similar to that of PAP patients. We suggest that AIP may function as a candidate gatekeeper gene in somatotrophs. To clarify molecular mechanisms of tumorigenesis, we elucidated the expression of AIP-related molecules in human and mouse pituitary tumors. The expression of aryl hydrocarbon receptor nuclear translocator (ARNT) was reduced in mouse Aip-deficient adenomas, and similar ARNT reduction was also evident in human AIP mutation positive adenomas. This suggests that in addition to participating in the hypoxia pathway, estrogen receptor signaling and xenobiotic response pathways, ARNT may play a role in AIP-related tumorigenesis. We also studied the characteristics and the response to therapy of PAP patients and found them to have an aggressive disease phenotype with young age at onset. Therefore, improvement in treatment outcomes of PAP patients would require their efficient identification and earlier diagnosis of the pituitary adenomas. The possible role of the RET proto-oncogene in tumorigenesis of familial AIP mutation negative pituitary adenomas was evaluated, but none of the found RET germline variants were considered pathogenic. Surprisingly, RET immunohistochemistry suggested possible underexpression of RET in AIP mutation positive pituitary adenomas an observation that merits further investigation.

Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.

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.

Detection and molecular epidemiology of tick-borne encephalitis virus infections

Pdf-file, link above

Molecular Characterization of Viruses Causing the Cassava Brown Streak Disease Epidemic in Eastern Africa

Cassava brown streak disease (CBSD) was described for the first time in Tanganyika (now Tanzania) about seven decades ago. Tanganyika (now Tanzania) about seven decades ago. It was endemic in the lowland areas of East Africa and inland parts of Malawi and caused by Cassava brown streak virus (CBSV; genus Ipomovirus; Potyviridae). However, in 1990s CBSD was observed at high altitude areas in Uganda. The causes for spread to new locations were not known.The present work was thus initiated to generate information on genetic variability, clarify the taxonomy of the virus or viruses associated with CBSD in Eastern Africa as well as to understand the evolutionary forces acting on their genes. It also sought to develop a molecular based diagnostic tool for detection of CBSD-associated virus isolates. Comparison of the CP-encoding sequences of CBSD-associated virus isolates collected from Uganda and north-western Tanzania in 2007 and the partial sequences available in Genbank revealed occurrence of two genetically distinct groups of isolates. Two isolates were selected to represent the two groups. The complete genomes of isolates MLB3 (TZ:Mlb3:07) and Kor6 (TZ:Kor6:08) obtained from North-Western (Kagera) and North-Eastern (Tanga) Tanzania, respectively, were sequenced. The genomes were 9069 and 8995 nucleotides (nt), respectively. They translated into polyproteins that were predicted to yield ten mature proteins after cleavage. Nine proteins were typical in the family Potyviridae, namely P1, P3, 6K1, CI, 6K2, VPg, NIa-Pro, NIb and CP, but the viruses did not contain HC-Pro. Interestingly, genomes of both isolates contained a Maf/HAM1-like sequence (HAM1h; 678 nucleotides, 25 kDa) recombined between the NIb and CP domains in the 3’-proximal part of the genomes. HAM1h was also identified in Euphorbia ringspot virus (EuRSV) whose sequence was in GenBank. The HAM1 gene is widely spread in both prokaryotes and eukaryotes. In yeast (Saccharomyces cerevisiae) it is known to be a nucleoside triphosphate (NTP) pyrophosphatase. Novel information was obtained on the structural variation at the N-termini of polyproteins of viruses in the genus Ipomovirus. Cucumber vein yellowing virus (CVYV) and Squash vein yellowing virus (SqVYV) contain a duplicated P1 (P1a and P1b) but lack the HC-Pro. On the other hand, Sweet potato mild mottle virus (SPMMV), has a single but large P1 and has HC-Pro. Both virus isolates (TZ:Mlb3:07 & TZ:Kor6:08) characterized in this study contained a single P1 and lacked the HC-Pro which indicates unique evolution in the family Potyviridae. Comparison of 12 complete genomes of CBSD-associated viruses which included two genomes characterized in this study, revealed genetic identity of 69.0–70.3% (nt) and amino acid (aa) identities of 73.6–74.4% at polyprotein level. Comparison was also made among 68 complete CP sequences, which indicated 69.0-70.3 and 73.6-74.4 % identity at nt and aa levels, respectively. The genetic variation was large enough for dermacation of CBSD-associated virus isolates into two distinct species. The name CBSV was retained for isolates that were related to CBSV isolates available in database whereas the new virus described for the first time in this study was named Ugandan cassava brown streak virus (UCBSV) by the International Committee on Virus Taxonomy (ICTV). The isolates TZ:Mlb3:07 and TZ:Kor6:08 belong to UCBSV and CBSV, respectively. The isolates of CBSV and UCBSV were 79.3-95.5% and 86.3-99.3 % identitical at nt level, respectively, suggesting more variation amongst CBSV isolates. The main sources of variation in plant viruses are mutations and recombination. Signals for recombination events were detected in 50% of isolates of each virus. Recombination events were detected in coding and non-coding (3’-UTR) sequences except in the 5’UTR and P3. There was no evidence for recombination between isolates of CBSV and UCBSV. The non-synonomous (dN) to synonomous (dS) nucleotide substitution ratio (ω) for the HAM1h and CP domains of both viruses were ≤ 0.184 suggesting that most sites of these proteins were evolving under strong purifying selection. However, there were individual amino acid sites that were submitted to adaptive evolution. For instance, adaptive evolution was detected in the HAM1h of UCBSV (n=15) where 12 aa sites were under positive selection (P< 0.05) but not in CBSV (n=12). The CP of CBSV (n=23) contained 12 aa sites (p<0.01) while only 5 aa sites in the CP gene of UCBSV were predicted to be submitted to positive selection pressure (p<0.01). The advantages offered by the aa sites under positive selection could not be established but occurrence of such sites in the terminal ends of UCBSV-HAMIh, for example, was interpreted as a requirement for proteolysis during polyprotein processing. Two different primer pairs that simultaneously detect UCBSV and CBSV isolates were developed in this study. They were used successfully to study distribution of CBSV, UCBSV and their mixed infections in Tanzania and Uganda. It was established that the two viruses co-infect cassava and that incidences of co-infection could be as high as 50% around Lake Victoria on the Tanzanian side. Furthermore, it was revealed for the first time that both UCBSV and CBSV were widely distributed in Eastern Africa. The primer pair was also used to confirm infection in a close relative of cassava, Manihot glaziovii (Müller Arg.) with CBSV. DNA barcoding of M. glaziovii was done by sequencing the matK gene. Two out of seven M. glaziovii from the coastal areas of Korogwe and Kibaha in north eastern Tanzania were shown to be infected by CBSV but not UCBSV isolates. Detection in M. glaziovii has an implication in control and management of CBSD as it is likely to serve as virus reservoir. This study has contributed to the understanding of evolution of CBSV and UCBSV, which cause CBSD epidemic in Eastern Africa. The detection tools developed in this work will be useful in plant breeding, verification of the phytosanitary status of materials in regional and international movement of germplasm, and in all diagnostic activities related to management of CBSD. Whereas there are still many issues to be resolved such as the function and biological significance of HAM1h and its origin, this work has laid a foundation upon which the studies on these aspects can be based.

Molecular genetics of Usher syndrome -inherited deafness and blindness

Usher syndrome (USH) is an inherited blindness and deafness disorder with variable vestibular dysfunction. The syndrome is divided into three subtypes according to the progression and severity of clinical symptoms. The gene mutated in Usher syndrome type 3 (USH3), clarin 1 (CLRN1), was identified in Finland in 2001 and two mutations were identified in Finnish patients at that time. Prior to this thesis study, the two CLRN1 gene mutations were the only USH mutations identified in Finnish USH patients. To further clarify the Finnish USH mutation spectrum, all nine USH genes were studied. Seven mutations were identified: one was a previously known mutation in CLRN1, four were novel mutations in myosin VIIa (MYO7A) and two were a novel and a previously known mutation in usherin (USH2A). Another aim of this thesis research was to further study the structure and function of the CLRN1 gene, and to clarify the effects of mutations on protein function. The search for new splice variants resulted in the identification of eight novel splice variants in addition to the three splice variants that were already known prior to this study. Studies of the possible promoter regions for these splice variants showed the most active region included the 1000 bases upstream of the translation start site in the first exon of the main three exon splice variant. The 232 aa CLRN1 protein encoded by the main (three-exon) splice variant was transported to the plasma membrane when expressed in cultured cells. Western blot studies suggested that CLRN1 forms dimers and multimers. The CLRN1 mutant proteins studied were retained in the endoplasmic reticulum (ER) and some of the USH3 mutations caused CLRN1 to be unstable. During this study, two novel CLRN1 sequence alterations were identified and their pathogenicity was studied with cell culture protein expression. Previous studies with mice had shown that Clrn1 is expressed in mouse cochlear hair cells and spiral ganglion cells, but the expression profile in mouse retina remained unknown. The Clrn1 knockout mice display cochlear cell disruption/death, but do not have a retinal phenotype. The zebrafish, Danio rerio, clrn1 was found to be expressed in hair cells associated with hearing and balance. Clrn1 expression was also found in the inner nuclear layer (INL), photoreceptor layer and retinal pigment epithelium layer (RPE) of the zebrafish retina. When Clrn1 production was knocked down with injected morpholino oligonucleotides (MO) targeting Clrn1 translation or correct splicing, the zebrafish larvae showed symptoms similar to USH3 patients. These larvae had balance/hearing problems and reduced response to visual stimuli. The knowledge this thesis research has provided about the mutations in USH genes and the Finnish USH mutation spectrum are important in USH patient diagnostics. The extended information about the structure and function of CLRN1 is a step further in exploring USH3 pathogenesis caused by mutated CLRN1 as well as a step in finding a cure for the disease.

Frequency-comb-referenced molecular spectroscopy in the mid-infrared region

A simple method for absolute frequency measurements of molecular transitions in the mid-infrared region is reported. The method is based on a cw singly-resonant optical parametric oscillator (SRO), which is tunable from 3.2 to 3.45 µm. The mid- infrared frequency of the SRO is referenced to an optical frequency comb through its pump and signal beams. Sub-Doppler spectroscopy and absolute frequency measurement of the P(7) transition of the ν3 band of CH4 are demonstrated.