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.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

Cellular and molecular regulation of mammalian blastocyst hatching

In mammals including humans, failure in blastocyst hatching and implantation leads to early embryonic loss and infertility. Prior to implantation, the blastocyst must hatch out of its acellular glycoprotein coat, the zona pellucida (ZP). The phenomenon of blastocyst hatching is believed to be regulated by (i) dynamic cellular components such as actin-based trophectodermal projections (TEPs), and (ii) a variety of autocrine and paracrine molecules such as growth factors, cytokines and proteases. The spatio-temporal regulation of zona lysis by blastocyst-derived cellular and molecular signaling factors is being keenly investigated. Our studies show that hamster blastocyst hatching is acelerated by growth factors such as heparin binding-epidermal growth factor and leukemia inhibitory factor and that embryo-derived, cysteine proteases including cathepsins are responsible for blastocyst hatching. Additionally, we believe that cyclooxygenase-generated prostaglandins, estradiol-17 beta mediated estrogen receptor-alpha signaling and possibly NF kappa B could be involved in peri-hatching development. Moreover, we show that TEPs are intimately involved with lysing ZP and that the TEPs potentially enrich and harbor hatching-enabling factors. These observations provide new insights into our understanding of the key cellular and molecular regulators involved in the phenomenon of mammalian blastocyst hatching, which is essential for the establishment of early pregnancy.

Molecular genetics of X-linked cone-rod dystrophy and Åland Island eye disease

Inherited retinal diseases are the most common cause of vision loss among the working population in Western countries. It is estimated that ~1 of the people worldwide suffer from vision loss due to inherited retinal diseases. The severity of these diseases varies from partial vision loss to total blindness, and at the moment no effective cure exists. To date, nearly 200 mapped loci, including 140 cloned genes for inherited retinal diseases have been identified. By a rough estimation 50% of the retinal dystrophy genes still await discovery. In this thesis we aimed to study the genetic background of two inherited retinal diseases, X-linked cone-rod dystrophy and Åland Island eye disease. X-linked cone-rod dystrophy (CORDX) is characterized by progressive loss of visual function in school age or early adulthood. Affected males show reduced visual acuity, photophobia, myopia, color vision defects, central scotomas, and variable changes in fundus. The disease is genetically heterogeneous and two disease loci, CORDX1 and CORDX2, were known prior to the present thesis work. CORDX1, located on Xp21.1-11.4, is caused by mutations in the RPGR gene. CORDX2 is located on Xq27-28 but the causative gene is still unknown. Åland Island eye disease (AIED), originally described in a family living in Åland Islands, is a congenital retinal disease characterized by decreased visual acuity, fundus hypopigmentation, nystagmus, astigmatism, red color vision defect, myopia, and defective night vision. AIED shares similarities with another retinal disease, congenital stationary night blindness (CSNB2). Mutations in the L-type calcium channel α1F-subunit gene, CACNA1F, are known to cause CSNB2, as well as AIED-like disease. The disease locus of the original AIED family maps to the same genetic interval as the CACNA1F gene, but efforts to reveal CACNA1F mutations in patients of the original AIED family have been unsuccessful. The specific aims of this study were to map the disease gene in a large Finnish family with X-linked cone-rod dystrophy and to identify the disease-causing genes in the patients of the Finnish cone-rod dystrophy family and the original AIED family. With the help of linkage and haplotype analyses, we could localize the disease gene of the Finnish cone-rod dystrophy family to the Xp11.4-Xq13.1 region, and thus establish a new genetic X-linked cone-rod dystrophy locus, CORDX3. Mutation analyses of candidate genes revealed three novel CACNA1F gene mutations: IVS28-1 GCGTC>TGG in CORDX3 patients, a 425 bp deletion, comprising exon 30 and flanking intronic regions in AIED patients, and IVS16+2T>C in an additional Finnish patient with a CSNB2-like phenotype. All three novel mutations altered splice sites of the CACNA1F gene, and resulted in defective pre-mRNA splicing suggesting altered or absent channel function as a disease mechanism. The analyses of CACNA1F mRNA also revealed novel alternative wt splice variants, which may enhance channel diversity or regulate the overall expression level of the channel. The results of our studies may be utilized in genetic counseling of the families, and they provide a basis for studies on the pathogenesis of these diseases. In the future, the knowledge of the genetic defects may be used in the identification of specific therapies for the patients.

Risk factors for open angle glaucoma; a clinical and molecular genetic study

Glaucoma, optic neuropathy with excavation in the optic nerve head and corresponding visual field defect, is one of the leading causes for blindness worldwide. However, visual disability can often be avoided or delayed if the disease is diagnosed at an early stage. Therefore, recognising the risk factors for development and progression of glaucoma may prevent further damage. The purpose of the present study was to evaluate factors associated with visual disability caused by glaucoma and the genetic features of two risk factors, exfoliation syndrome (ES) and a positive family history of glaucoma. The present study material consisted of three study groups 1) deceased glaucoma patients from the Ekenäs practice 2) glaucoma families from the Ekenäs region and 3) population based families with and without exfoliation syndrome from Kökar Island. For the retrospective study, 106 patients with open angle glaucoma (OAG) were identified. At the last visit, 17 patients were visually impaired. Blindness induced by glaucoma was found in one or both eyes in 16 patients and in both eyes in six patients. The cumulative incidence of glaucoma caused blindness for one eye was 6% at 5 years, 9% at 10 years, and 15% at 15 years from initialising the treatment. The factors associated with blindness caused by glaucoma were an advanced stage of glaucoma at diagnosis, fluctuation in intraocular pressure during treatment, the presence of exfoliation syndrome, and poor patient compliance. A cross-sectional population based study performed in 1960-1962 on Kökar Island and the same population was followed until 2002. In total 965 subjects (530 over 50 years) have been examined at least once. The prevalence of exfoliation syndrome (ES) was 18% among subjects older than 50 years. Seventy-five of all 78 ES-positives belonged to the same extended pedigree. According to the segregation and family analysis, exfoliation syndrome seemed to be inherited as an autosomal dominant trait with reduced penetrance. The penetrance was more reduced for males, but the risk for glaucoma was higher in males than in females. To find the gene or genes associated with exfoliation syndrome, a genome wide scan was performed for 64 members (28 ES affected and 36 controls) of the Kökar pedigree. A promising result was found: the highest two-point LOD score of 3.45 (θ=0.04) in chromosome18q12.1-21.33. The presence of mutations in glaucoma genes TIGR/MYOC (myocilin) and OPTN (optineurin) was analysed in eight glaucoma families from the Ekenäs region. An inheritance pattern resembling autosomal dominant mode was detected in all these families. Primary open angle glaucoma or exfoliation glaucoma was found in 35% of 136 family members and 28% were suspected to have glaucoma. No mutations were detected in these families.

Molecular mediators linking stroke and carotid artery disease

Carotid artery disease is the most prevalent etiologic precursor of ischemic stroke, which is a major health hazard and the second most common cause of death in the world. If a patient presents with a symptomatic high-grade (>70%) stenosis in the internal carotid artery, the treatment of choice is carotid endarterectomy. However, the natural course of radiologically equivalent carotid lesions may be clinically quite diverse, and the reason for that is unknown. It would be of utmost importance to develop molecular markers that predict the symptomatic phenotype of an atherosclerotic carotid plaque (CP) and help to differentiate vulnerable lesions from stable ones. The aim of this study was to investigate the morphologic and molecular factors that associate with stroke-prone CPs. In addition to immunohistochemistry, DNA microarrays were utilized to identify molecular markers that would differentiate between symptomatic and asymptomatic CPs. Endothelial adhesion molecule expression (ICAM-1, VCAM-1, P-selectin, and E-selectin) did not differ between symptomatic and asymptomatic patients. Denudation of endothelial cells was associated with symptom-generating carotid lesions, but in studies on the mechanism of decay of endothelial cells, markers of apoptosis (TUNEL, activated caspase 3) were found to be decreased in the endothelium of symptomatic lesions. Furthermore, markers of endothelial apoptosis were directly associated with those of cell proliferation (Ki-67) in all plaques. FasL expression was significantly increased on the endothelium of symptomatic CPs. DNA microarray analysis revealed prominent induction of specific genes in symptomatic CPs, including those subserving iron and heme metabolism, namely HO-1, and hemoglobin scavenger receptor CD163. HO-1 and CD163 proteins were also increased in symptomatic CPs and associated with intraplaque iron deposits, which, however, did not correlate with symptom status itself. ADRP, the gene for adipophilin, was also overexpressed in symptomatic CPs. Adipophilin expression was markedly increased in ulcerated CPs and colocalized with extravasated red blood cells and cholesterol crystals. Taken together, the phenotypic characteristics and the numerous possible molecular mediators of the destabilization of carotid plaques provide potential platforms for future research. The denudation of the endothelial lining observed in symptomatic CPs may lead to direct thromboembolism and maintain harmful oxidative and inflammatory processes, predispose to plaque microhemorrhages, and contribute to lipid accumulation into the plaque, thereby making it vulnerable to rupture.

Migraine and glutamate: Modulators of glutamatergic signalling as potential treatments of neuropsychiatric disorders

Migraine is a complex neurological disorder with a well-documented genetic basis. Migraine is a product of allelic variation in genes of neurological, vascular and hormonal origin interacting with environmental triggers. Presentation can include attacks of head pain with symptoms of nausea, emesis, photophobia, phonophobia, and occasionally, visual sensory disturbances, known as aura. Migraine pain is difficult to ignore, associated with a deep sense of malaise and manifests as a throbbing, pulsatile headache, localized to one side of the head that intensifies with physical activity and that can last from 4-72 hours. Migraine is diagnosed according to criteria developed by the International Headache Society (IHS) and is subdivided into two main types based on the occurrence of aura symptoms that may be present in the early stages of the headache: migraine with aura (MA) and migraine without aura (MO). The majority (about 70%) of migraineurs are diagnosed with the MO subtype whilst the remaining 30% experience MA accompanied by neurological symptoms that manifest as fully reversible, visual, sensory and/or dysphasic speech disturbances in conjunction with their headache. Glutamate is the primary excitatory neurotransmitter in the central nervous system (CNS) and over-excitation of glutamate receptors is regarded as a contributing factor, through various mechanisms, to the pathology of migraine. In this chapter we present an overview of the pathophysiology and co-morbidity of migraine with other psychiatric disorders and discuss the role of the glutamatergic system in migraine, its molecular components as potential drug targets, in addition to the current treatments and progress of modulators of glutamatergic signaling.

Quantum Interference in Coherent Molecular Conductance

Coherent electronic transport through individual molecules is crucially sensitive to quantum interference. We investigate the zero-bias and zero-temperature conductance through pi-conjugated annulene molecules weakly coupled to two leads for different source-drain configurations, finding an important reduction for certain transmission channels and for particular geometries as a consequence of destructive quantum interference between states with definite momenta. When translational symmetry is broken by an external perturbation we find an abrupt increase of the conductance through those channels. Previous studies concentrated on the effect at the Fermi energy, where this effect is very small. By analyzing the effect of symmetry breaking on the main transmission channels we find a much larger response thus leading to the possibility of a larger switching of the conductance through single molecules.

Interactions and turnover of the muscular dystrophy protein myotilin

The striated muscle sarcomere is a force generating and transducing unit as well as an important sensor of extracellular cues and a coordinator of cellular signals. The borders of individual sarcomeres are formed by the Z-disks. The Z-disk component myotilin interacts with Z-disk core structural proteins and with regulators of signaling cascades. Missense mutations in the gene encoding myotilin cause dominantly inherited muscle disorders, myotilinopathies, by an unknown mechanism. In this thesis the functions of myotilin were further characterized to clarify the molecular biological basis and the pathogenetic mechanisms of inherited muscle disorders, mainly caused by mutated myotilin. Myotilin has an important function in the assembly and maintenance of the Z-disks probably through its actin-organizing properties. Our results show that the Ig-domains of myotilin are needed for both binding and bundling actin and define the Ig domains as actin-binding modules. The disease-causing mutations appear not to change the interplay between actin and myotilin. Interactions between Z-disk proteins regulate muscle functions and disruption of these interactions results in muscle disorders. Mutations in Z-disk components myotilin, ZASP/Cypher and FATZ-2 (calsarcin-1/myozenin-2) are associated with myopathies. We showed that proteins from the myotilin and FATZ families interact via a novel and unique type of class III PDZ binding motif with the PDZ domains of ZASP and other Enigma family members and that the interactions can be modulated by phosphorylation. The morphological findings typical of myotilinopathies include Z-disk alterations and aggregation of dense filamentous material. The causes and mechanisms of protein aggregation in myotilinopathy patients are unknown, but impaired degradation might explain in part the abnormal protein accumulation. We showed that myotilin is degraded by the calcium-dependent, non-lysosomal cysteine protease calpain and by the proteasome pathway, and that wild type and mutant myotilin differ in their sensitivity to degradation. These studies identify the first functional difference between mutated and wild type myotilin. Furthermore, if degradation of myotilin is disturbed, it accumulates in cells in a manner resembling that seen in myotilinopathy patients. Based on the results, we propose a model where mutant myotilin escapes proteolytic breakdown and forms protein aggregates, leading to disruption of myofibrils and muscular dystrophy. In conclusion, the main results of this study demonstrate that myotilin is a Z-disk structural protein interacting with several Z-disk components. The turnover of myotilin is regulated by calpain and the ubiquitin proteasome system and mutations in myotilin seem to affect the degradation of myotilin, leading to protein accumulations in cells. These findings are important for understanding myotilin-linked muscle diseases and designing treatments for these disorders.

Prognostic molecular factors and algorithms in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common of the non-Hodgkin lymphomas. As DLBCL is characterized by heterogeneous clinical and biological features, its prognosis varies. To date, the International Prognostic Index has been the strongest predictor of outcome for DLBCL patients. However, no biological characters of the disease are taken into account. Gene expression profiling studies have identified two major cell-of-origin phenotypes in DLBCL with different prognoses, the favourable germinal centre B-cell-like (GCB) and the unfavourable activated B-cell-like (ABC) phenotypes. However, results of the prognostic impact of the immunohistochemically defined GCB and non-GCB distinction are controversial. Furthermore, since the addition of the CD20 antibody rituximab to chemotherapy has been established as the standard treatment of DLBCL, all molecular markers need to be evaluated in the post-rituximab era. In this study, we aimed to evaluate the predictive value of immunohistochemically defined cell-of-origin classification in DLBCL patients. The GCB and non-GCB phenotypes were defined according to the Hans algorithm (CD10, BCL6 and MUM1/IRF4) among 90 immunochemotherapy- and 104 chemotherapy-treated DLBCL patients. In the chemotherapy group, we observed a significant difference in survival between GCB and non-GCB patients, with a good and a poor prognosis, respectively. However, in the rituximab group, no prognostic value of the GCB phenotype was observed. Likewise, among 29 high-risk de novo DLBCL patients receiving high-dose chemotherapy and autologous stem cell transplantation, the survival of non-GCB patients was improved, but no difference in outcome was seen between GCB and non-GCB subgroups. Since the results suggested that the Hans algorithm was not applicable in immunochemotherapy-treated DLBCL patients, we aimed to further focus on algorithms based on ABC markers. We examined the modified activated B-cell-like algorithm based (MUM1/IRF4 and FOXP1), as well as a previously reported Muris algorithm (BCL2, CD10 and MUM1/IRF4) among 88 DLBCL patients uniformly treated with immunochemotherapy. Both algorithms distinguished the unfavourable ABC-like subgroup with a significantly inferior failure-free survival relative to the GCB-like DLBCL patients. Similarly, the results of the individual predictive molecular markers transcription factor FOXP1 and anti-apoptotic protein BCL2 have been inconsistent and should be assessed in immunochemotherapy-treated DLBCL patients. The markers were evaluated in a cohort of 117 patients treated with rituximab and chemotherapy. FOXP1 expression could not distinguish between patients, with favourable and those with poor outcomes. In contrast, BCL2-negative DLBCL patients had significantly superior survival relative to BCL2-positive patients. Our results indicate that the immunohistochemically defined cell-of-origin classification in DLBCL has a prognostic impact in the immunochemotherapy era, when the identifying algorithms are based on ABC-associated markers. We also propose that BCL2 negativity is predictive of a favourable outcome. Further investigational efforts are, however, warranted to identify the molecular features of DLBCL that could enable individualized cancer therapy in routine patient care.

Childhood-onset mitochondrial diseases : with emphasis on molecular diagnostics of neurological disorders

Childhood-onset mitochondrial diseases comprise a heterogeneous group of disorders, which may manifest with almost any symptom and affect any tissue or organ. Due to challenging diagnostics, most children still lack a specific aetiological diagnosis. The aim of this thesis was to find molecular causes for childhood-onset mitochondrial disorders in Finland. We identified the underlying cause for 25 children, and found three new diseases, which had not been diagnosed in Finland before. These diseases caused severe progressive infantile-onset encephalomyopathies, and were due to defects in mitochondrial DNA (mtDNA) maintenance. Furthermore, the thesis provides the molecular background of Finnish patients with ‘leukoencephalopathy with brain stem and spinal cord involvement and elevated brain lactate’ (LBSL). A new phenotype was identified to be due to mutations in Twinkle, resembling ‘infantile onset spinocerebellar ataxia’ (IOSCA). These mutations caused mtDNA depletion in the liver, thus confirming the essential role of Twinkle in mtDNA maintenance, and expanding the molecular background of mtDNA depletion syndromes. The major aetiology for infantile mitochondrial myopathy in Finland was discovered to be due to mutations in thymidine kinase 2 (TK2). A novel mutation with Finnish ancestry was identified, and a genotype-phenotype correlation with mutation-specific distribution of tissue involvement was found, thus proving that deficient TK2 may cause multi-tissue depletion and impair neuronal function. This work established the molecular diagnosis and advanced the knowledge of phenotypes among paediatric patients with polymerase gamma (POLG) mutations. The patients showed severe early-onset encephalopathy with intractable epilepsy. POLG mutations are not a prevalent cause of children’s ataxias, although ataxia is a major presenting symptom among adults. Our findings indicate that POLG mutations should be investigated even if typical MRI, histochemical or biochemical abnormalities are lacking. LBSL patients showed considerable variation in phenotype despite identical mutations. A common, most likely European, ancestry, and a relative high carrier frequency of these mutations in Finland were discovered; suggesting that LBSL may be a quite common leukoencephalopathy in other populations as well. The results suggest that MRI findings are so unique that the diagnosis of LBSL is possible to make without genetic studies. This thesis work has resulted in identification of new mitochondrial disorders in Finland, enhancing the understanding of the clinical variability and the importance of tissue-specificity of these disorders. In addition to providing specific diagnosis to the patients, these findings give light to the underlying pathogenetic mechanisms of childhood-onset mitochondrial disorders.