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.

Evolutionary and conservation biology of the Finnish house sparrow

Recently it has been recognized that evolutionary aspects play a major role in conservation issues of a species. In this thesis I have combined evolutionary research with conservation studies to provide new insight into these fields. The study object of this thesis is the house sparrow, a species that has features that makes it interesting for this type of study. The house sparrow has been ubiquitous almost all over the world. Even though being still abundant, several countries have reported major declines. These declines have taken place in a relatively short time covering both urban and rural habitats. In Finland this species has declined by more than two thirds in just over two decades. In addition, as the house sparrow lives only in human inhabited areas it can also raise public awareness to conservation issues. I used both an extensive museum collection of house sparrows collected in 1980s from all over Finland as well as samples collected in 2009 from 12 of the previously collected localities. I used molecular techniques to study neutral genetic variation within and genetic differentiation between the study populations. This knowledge I then combined with data gathered on morphometric measurements. In addition I analyzed eight heavy metals from the livers of house sparrows that lived in either rural or urban areas in the 1980s and evaluated the role of heavy metal pollution as a possible cause of the declines. Even though dispersal of house sparrows is limited I found that just as the declines started in 1980s the house sparrows formed a genetically panmictic population on the scale of the whole Finland. When compared to Norway, where neutral genetic divergence has been found even with small geographic distances, I concluded that this difference would be due to contrasting landscapes. In Finland the landscape is rather homogeneous facilitating the movements of these birds and maintaining gene flow even with the low dispersal. To see whether the declines have had an effect on the neutral genetic variation of the populations I did a comparison between the historical and contemporary genetic data. I showed that even though genetic diversity has not decreased due to the drastic declines the populations have indeed become more differentiated from each other. This shows that even in a still quite abundant species the declines can have an effect on the genetic variation. It is shown that genetic diversity and differentiation may approach their new equilibriums at different rates. This emphasizes the importance of studying both of them and if the latter has increased it should be taken as a warning sign of a possible loss of genetic diversity in the future. One of the factors suggested to be responsible for the house sparrow declines is heavy metal pollution. When studying the livers of house sparrows from 1980s I discovered higher levels of heavy metal concentrations in urban than rural habitats, but the levels of the metals were comparatively low and based on that heavy metal pollution does not seem to be a direct cause for the declines in Finland. However, heavy metals are known to decrease the amount of insects in urban areas and thus in the cities heavy metals may have an indirect effect on house sparrows. Although neutral genetic variation is an important tool for conservation genetics it does not tell the whole story. Since neutral genetic variation is not affected by selection, information can be one-sided. It is possible that even neutral genetic differentiation is low, there can be substantial variation in additive genetic traits indicating local adaptation. Therefore I performed a comparison between neutral genetic differentiation and phenotypic differentiation. I discovered that two traits out of seven are likely to be under directional selection, whereas the others could be affected by random genetic drift. Bergmann s rule may be behind the observed directional selection in wing length and body mass. These results highlight the importance of estimating both neutral and adaptive genetic variation.

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

Molecular Interactions Chlorin Assemblies and NMR Separations

Photosynthesis is a chemical process in which the energy of the light quanta is transformed into chemical energy. Chlorophyll (Chl) molecules play a key role in photosynthesis; they function in the antennae systems and in the photosynthetic reaction center where the primary charge separation (CS) takes place. Bio-inspired mimicry of the CS is an essential unit in dye-sensitized solar cells. Aim of this study was to design and develop electron donor-acceptor (EDA) pairs from Chls and fullerenes (C60) or carbon nanotubes (CNT). The supramolecular approach was chosen, as long synthetic sequences required by the covalent approach lead to long reaction schemes and low yields. Here, a π-interaction between soluble CNTs and Chl was used in EDA construction. Also, a beta-face selective two-point bound Chl-C60 EDA was introduced. In addition, the photophysical properties of the supramolecular EDA dyads were analyzed. In organic chemistry, nuclear magnetic resonance (NMR) spectroscopy is the most vital analytical technique in use. Multi-dimensional NMR experiments have enabled a structural analysis of complex natural products and proteins. However, in mixture analysis NMR is still facing difficulties. In many cases overlapping signals can t be resolved even with the help of multi-dimensional experiments. In this work, an NMR tool based on simple host-guest chemistry between analytes and macromolecules was developed. Diffusion ordered NMR spectroscopy (DOSY) measures the mobilities of compounds in an NMR sample. In a liquid state NMR sample, each of the analytes has a characteristic diffusion coefficient, which is proportional to the size of the analyte. With normal DOSY experiment, provided that the diffusion coefficients of the analytes differ enough, individual spectra of analytes can be extracted. When similar sized analytes differ chemically, an additive can be introduced into the sample. Since macromolecules in a liquid state NMR sample can be considered practically stationary, even faint supramolecular interaction can change the diffusion coefficient of the analyte sufficiently for a successful resolution in DOSY. In this thesis, polyvinylpyrrolidone and polyethyleneglycol enhanced DOSY NMR techniques, which enable mixture analysis of similar in size but chemically differing natural products, are introduced.

Molecular Modelling of Estrogen-Producing Enzymes CYP450 Aromatase and 17beta-Hydroxysteroid Dehydrogenase Type 1

Breast cancer is the most common cancer in women in the western countries. Approximately two-thirds of breast cancer tumours are hormone dependent, requiring estrogens to grow. Estrogens are formed in the human body via a multistep route starting from cholesterol. The final steps in the biosynthesis include the CYP450 aromatase enzyme, converting the male hormones androgens (preferred substrate androstenedione ASD) into estrogens(estrone E1), and the 17beta-HSD1 enzyme, converting the biologically less active E1 into the active hormone 17beta-hydroxyestradiol E2. E2 is bound to the nuclear estrogen receptors causing a cascade of biochemical reactions leading to cell proliferation in normal tissue, and to tumour growth in cancer tissue. Aromatase and 17beta-HSD1 are expressed in or near the breast tumour, locally providing the tissue with estrogens. One approach in treating hormone dependent breast tumours is to block the local estrogen production by inhibiting these two enzymes. Aromatase inhibitors are already on the market in treating breast cancer, despite the lack of an experimentally solved structure. The structure of 17beta-HSD1, on the other hand, has been solved, but no commercial drugs have emerged from the drug discovery projects reported in the literature. Computer-assisted molecular modelling is an invaluable tool in modern drug design projects. Modelling techniques can be used to generate a model of the target protein and to design novel inhibitors for them even if the target protein structure is unknown. Molecular modelling has applications in predicting the activities of theoretical inhibitors and in finding possible active inhibitors from a compound database. Inhibitor binding at atomic level can also be studied with molecular modelling. To clarify the interactions between the aromatase enzyme and its substrate and inhibitors, we generated a homology model based on a mammalian CYP450 enzyme, rabbit progesterone 21-hydroxylase CYP2C5. The model was carefully validated using molecular dynamics simulations (MDS) with and without the natural substrate ASD. Binding orientation of the inhibitors was based on the hypothesis that the inhibitors coordinate to the heme iron, and were studied using MDS. The inhibitors were dietary phytoestrogens, which have been shown to reduce the risk for breast cancer. To further validate the model, the interactions of a commercial breast cancer drug were studied with MDS and ligand–protein docking. In the case of 17beta-HSD1, a 3D QSAR model was generated on the basis of MDS of an enzyme complex with active inhibitor and ligand–protein docking, employing a compound library synthesised in our laboratory. Furthermore, four pharmacophore hypotheses with and without a bound substrate or an inhibitor were developed and used in screening a commercial database of drug-like compounds. The homology model of aromatase showed stable behaviour in MDS and was capable of explaining most of the results from mutagenesis studies. We were able to identify the active site residues contributing to the inhibitor binding, and explain differences in coordination geometry corresponding to the inhibitory activity. Interactions between the inhibitors and aromatase were in agreement with the mutagenesis studies reported for aromatase. Simulations of 17beta-HSD1 with inhibitors revealed an inhibitor binding mode with hydrogen bond interactions previously not reported, and a hydrophobic pocket capable of accommodating a bulky side chain. Pharmacophore hypothesis generation, followed by virtual screening, was able to identify several compounds that can be used in lead compound generation. The visualisation of the interaction fields from the QSAR model and the pharmacophores provided us with novel ideas for inhibitor development in our drug discovery project.

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.

Methicillin-resistant Staphylococcus aureus in Finland: recent changes in the epidemiology, long-term facility aspects, and phenotypic and molecular detection of isolates

Staphylococcus aureus is one of the most important bacteria that cause disease in humans, and methicillin-resistant S. aureus (MRSA) has become the most commonly identified antibiotic-resistant pathogen in many parts of the world. MRSA rates have been stable for many years in the Nordic countries and the Netherlands with a low MRSA prevalence in Europe, but in the recent decades, MRSA rates have increased in those low-prevalence countries as well. MRSA has been established as a major hospital pathogen, but has also been found increasingly in long-term facilities (LTF) and in communities of persons with no connections to the health-care setting. In Finland, the annual number of MRSA isolates reported to the National Infectious Disease Register (NIDR) has constantly increased, especially outside the Helsinki metropolitan area. Molecular typing has revealed numerous outbreak strains of MRSA, some of which have previously been associated with community acquisition. In this work, data on MRSA cases notified to the NIDR and on MRSA strain types identified with pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), and staphylococcal cassette chromosome mec (SCCmec) typing at the National Reference Laboratory (NRL) in Finland from 1997 to 2004 were analyzed. An increasing trend in MRSA incidence in Finland from 1997 to 2004 was shown. In addition, non-multi-drug resistant (NMDR) MRSA isolates, especially those resistant only to methicillin/oxacillin, showed an emerging trend. The predominant MRSA strains changed over time and place, but two internationally spread epidemic strains of MRSA, FIN-16 and FIN-21, were related to the increase detected most recently. Those strains were also one cause of the strikingly increasing invasive MRSA findings. The rise of MRSA strains with SCCmec types IV or V, possible community-acquired MRSA was also detected. With questionnaires, the diagnostic methods used for MRSA identification in Finnish microbiology laboratories and the number of MRSA screening specimens studied were reviewed. Surveys, which focused on the MRSA situation in long-term facilities in 2001 and on the background information of MRSA-positive persons in 2001-2003, were also carried out. The rates of MRSA and screening practices varied widely across geographic regions. Part of the NMDR MRSA strains could remain undetected in some laboratories because of insufficient diagnostic techniques used. The increasing proportion of elderly population carrying MRSA suggests that MRSA is an emerging problem in Finnish long-term facilities. Among the patients, 50% of the specimens were taken on a clinical basis, 43% on a screening basis after exposure to MRSA, 3% on a screening basis because of hospital contact abroad, and 4% for other reasons. In response to an outbreak of MRSA possessing a new genotype that occurred in a health care ward and in an associated nursing home of a small municipality in Northern Finland in autumn 2003, a point-prevalence survey was performed six months later. In the same study, the molecular epidemiology of MRSA and methicillin-sensitive S. aureus (MSSA) strains were also assessed, the results to the national strain collection compared, and the difficulties of MRSA screening with low-level oxacillin-resistant isolates encountered. The original MRSA outbreak in LTF, which consisted of isolates possessing a nationally new PFGE profile (FIN-22) and internationally rare MLST type (ST-27), was confined. Another previously unrecognized MRSA strain was found with additional screening, possibly indicating that current routine MRSA screening methods may be insufficiently sensitive for strains possessing low-level oxacillin resistance. Most of the MSSA strains found were genotypically related to the epidemic MRSA strains, but only a few of them had received the SCCmec element, and all those strains possessed the new SCCmec type V. In the second largest nursing home in Finland, the colonization of S. aureus and MRSA, and the role of screening sites along with broth enrichment culture on the sensitivity to detect S. aureus were studied. Combining the use of enrichment broth and perineal swabbing, in addition to nostrils and skin lesions swabbing, may be an alternative for throat swabs in the nursing home setting, especially when residents are uncooperative. Finally, in order to evaluate adequate phenotypic and genotypic methods needed for reliable laboratory diagnostics of MRSA, oxacillin disk diffusion and MIC tests to the cefoxitin disk diffusion method at both +35°C and +30°C, both with or without an addition of sodium chloride (NaCl) to the Müller Hinton test medium, and in-house PCR to two commercial molecular methods (the GenoType® MRSA test and the EVIGENETM MRSA Detection test) with different bacterial species in addition to S. aureus were compared. The cefoxitin disk diffusion method was superior to that of oxacillin disk diffusion and to the MIC tests in predicting mecA-mediated resistance in S. aureus when incubating at +35°C with or without the addition of NaCl to the test medium. Both the Geno Type® MRSA and EVIGENETM MRSA Detection tests are usable, accurate, cost-effective, and sufficiently fast methods for rapid MRSA confirmation from a pure culture.

Molecular methods for the epidemiological analysis of methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae

Methicillin-resistant Staphylococcus aureus (MRSA) and Streptococcus pneumoniae are major health problems worldwide, both found in symptomless carriage but also causing even life-threatening infections. The aim of this thesis was to characterise MRSA and S. pneumoniae in detail by using several molecular typing methods for various epidemiological purposes: clonality analysis, epidemiological surveillance, outbreak investigation, and virulence factor analysis. The characteristics of MRSA isolates from the strain collection of the Finnish National Infectious Disease Register (NIDR) and pneumococcal isolates collected from military recruits and children with acute otitis media (AOM) were analysed using various typing techniques. Antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), spa typing, staphylococcal cassette chromosome mec (SCCmec) typing, and the detection of Panton-Valentine leukocidin (PVL) genes were performed for MRSA isolates. Pneumococcal isolates were analysed using antimicrobial susceptibility testing, serotyping, MLST, and by detecting pilus islet 1 (PI-1) and 2 (PI-2) genes. Several international community- and hospital-associated MRSA clones were recognised in Finland. The genetic diversity among MRSA FIN-4 isolates and among FIN-16 isolates was low. Overall, MRSA blood isolates from 1997 to 2006 were genetically diverse. spa typing was found to be a highly discriminatory, rapid and accurate typing method and it also qualifies as the primary typing method in countries with a long history of PFGE-based MRSA strain nomenclature. However, additional typing by another method, e.g. PFGE, is needed in certain situations to be able to provide adequate discrimination for epidemiological surveillance and outbreak investigation. An outbreak of pneumonia was associated with one pneumococcal strain among military recruits, previously healthy young men living in a crowded setting. The pneumococcal carriage rate after the outbreak was found to be exceptionally high. PI-1 genes were detected at a rather low prevalence among pneumococcal isolates from children with AOM. However, the study demonstrated that PI-1 has existed among pneumococcal isolates prior to pneumococcal conjugate vaccine and the increased antimicrobial resistance era. Moreover, PI-1 was found to associate with the serotype rather than the genotype. This study adds to our understanding of the molecular epidemiology of MRSA strains in Finland and the importance of an appropriate genotyping method to be able to perform high-level laboratory-based surveillance of MRSA. Epidemiological and molecular analyses of S. pneumoniae add to our knowledge of the characteristics of pneumococcal strains in Finland.

Molecular diagnosis of developmental disorders

ABSTRACT Idiopathic developmental disorders (DDs) affect ~1% of the population worldwide. This being a considerable amount, efforts are being made to elucidate the disease mechanisms. One or several genetic factors cause 30-40% of DDs, and only 10% are caused by environmental factors. The remaining 50% of DD patients go undiagnosed, mostly due to a lack of diagnostic techniques. The cause in most undiagnosed cases is though to be a genetic factor or a combination of genetic and environmental factors. Despite the surge of new technologies entering the market, their implementation into diagnostic laboratories is hampered by costs, lack of information about the expected diagnostic yield, and the wide range of selection. This study evaluates new microarray methods in diagnosing idiopathic DDs, providing information about their added diagnostic value. Study I analysed 150 patients by array comparative genomic hybridization (array CGH, 44K and 244K), with a subsequent 18% diagnostic yield. These results are supported by other studies, indicating an enourmous added diagnostic value of array CGH, compared with conventional cytogenetic analysis. Nevertheless, 80% of the patients remained undiagnosed in Study I. In an effort to diagnose more patients, in Study IV the resolution was increased from 8.9 Kb of the 244K CGH array to 0.7 Kb, by using a single-nucleotide polymorphism (SNP) array. However, no additional pathogenic changes were detected in the 35 patients assessed, and thus, for diagnostic purposes, an array platform with ca 9 Kb resolution appears adequate. The recent vast increase in reports of detected aberrations and associated phenotypes has enabled characterization of several new syndromes first based on a common aberration and thereafter by delineation of common clinical characteristics. In Study II, a familial deletion at 9q22.2q22.32 with variable penetrance was described. Despite several reports of aberrations in the adjacent area at 9q associated with Gorlin syndrome, the patients in this family had a unique phenotype and did not present with the syndrome. In Study III, a familial duplication of chromosome 6p22.2 was described. The duplication caused increased expression of an important enzyme of the γ-aminobutyric acid (GABA) degradation pathway, causing oxidative stress of the brain, and thus, very likely, the mild mental retardation of these patients. These two case studies attempted to pinpoint candidate genes and to resolve the pathogenic mechanism causing the clinical characteristics of the patients. Presenting rare genetic and clinical findings to the international science and medical community enables interpretation of similar findings in other patients. The added value of molecular karyotyping in patients with idiopathic DD is evident. As a first line of testing, arrays with a median resolution of at least 9 Kb should be considered and further characterization of detected aberrations undertaken when possible. Diagnostic whole-exome sequencing may be the best option for patients who remain undiagnosed after high-resolution array analysis.