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

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

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.

Cathepsin D Deficiency - Molecular and Cellular Mechanisms of Neurodegeneration

Cathepsin D (CTSD) is a lysosomal protease, the deficiency of which is fatal and associated with neurodegeneration. CTSD knock-out mice, which die at the age of four weeks, show intestinal necrosis, loss of lymphoid cells and moderate pathological changes in the brain. An active-site mutation in the CTSD gene underlies a neurodegenerative disease in newborn sheep, characterized by brain atrophy without any changes to visceral tissues. The CTSD deficiences belong to the group of neuronal ceroid-lipofuscinoses (NCLs), severe neurodegenerative lysosomal storage disorders. The aim of this thesis was to examine the molecular and cellular mechanisms behind neurodegeneration in CTSD deficiency. We found the developmental expression pattern of CTSD to resemble that of synaptophysin and the increasing expression of CTSD to coincide with the active period of myelination in the rat brain, suggesting a role for CTSD in early rat brain development. An active-site mutation underlying the congenital ovine NCL not only affected enzymatic activity, but also changed the stability, processing and transport of the mutant protein, possibly contributing to the disease pathogenesis. We also provide CTSD deficiency as a first molecular explanation for human congenital NCL, a lysosomal storage disorder, characterized by neuronal loss and demyelination in the central nervous system. Finally, we show the first evidence for synaptic abnormalities and thalamocortical changes in CTSD-deficient mice at the molecular and ultrastructural levels. Keywords: cathepsin D, congenital, cortex, lysosomal storage disorder, lysosome, mutation, neurodegeneration, neuronal ceroid-lipofuscinosis, overexpression, synapse, thalamus

Microarrays in molecular profiling of cancer - focus on head and neck squamous cell carcinoma

Microarrays have a wide range of applications in the biomedical field. From the beginning, arrays have mostly been utilized in cancer research, including classification of tumors into different subgroups and identification of clinical associations. In the microarray format, a collection of small features, such as different oligonucleotides, is attached to a solid support. The advantage of microarray technology is the ability to simultaneously measure changes in the levels of multiple biomolecules. Because many diseases, including cancer, are complex, involving an interplay between various genes and environmental factors, the detection of only a single marker molecule is usually insufficient for determining disease status. Thus, a technique that simultaneously collects information on multiple molecules allows better insights into a complex disease. Since microarrays can be custom-manufactured or obtained from a number of commercial providers, understanding data quality and comparability between different platforms is important to enable the use of the technology to areas beyond basic research. When standardized, integrated array data could ultimately help to offer a complete profile of the disease, illuminating mechanisms and genes behind disorders as well as facilitating disease diagnostics. In the first part of this work, we aimed to elucidate the comparability of gene expression measurements from different oligonucleotide and cDNA microarray platforms. We compared three different gene expression microarrays; one was a commercial oligonucleotide microarray and the others commercial and custom-made cDNA microarrays. The filtered gene expression data from the commercial platforms correlated better across experiments (r=0.78-0.86) than the expression data between the custom-made and either of the two commercial platforms (r=0.62-0.76). Although the results from different platforms correlated reasonably well, combining and comparing the measurements were not straightforward. The clone errors on the custom-made array and annotation and technical differences between the platforms introduced variability in the data. In conclusion, the different gene expression microarray platforms provided results sufficiently concordant for the research setting, but the variability represents a challenge for developing diagnostic applications for the microarrays. In the second part of the work, we performed an integrated high-resolution microarray analysis of gene copy number and expression in 38 laryngeal and oral tongue squamous cell carcinoma cell lines and primary tumors. Our aim was to pinpoint genes for which expression was impacted by changes in copy number. The data revealed that especially amplifications had a clear impact on gene expression. Across the genome, 14-32% of genes in the highly amplified regions (copy number ratio >2.5) had associated overexpression. The impact of decreased copy number on gene underexpression was less clear. Using statistical analysis across the samples, we systematically identified hundreds of genes for which an increased copy number was associated with increased expression. For example, our data implied that FADD and PPFIA1 were frequently overexpressed at the 11q13 amplicon in HNSCC. The 11q13 amplicon, including known oncogenes such as CCND1 and CTTN, is well-characterized in different type of cancers, but the roles of FADD and PPFIA1 remain obscure. Taken together, the integrated microarray analysis revealed a number of known as well as novel target genes in altered regions in HNSCC. The identified genes provide a basis for functional validation and may eventually lead to the identification of novel candidates for targeted therapy in HNSCC.

Molecular Mechanisms of Androgen Receptor Interactions

The androgen receptor (AR) mediates the effects of the male sex-steroid hormones (androgens), testosterone and 5?-dihydrotestosterone. Androgens are critical in the development and maintenance of male sexual characteristics. AR is a member of the steroid receptor ligand-inducible transcription factor family. The steroid receptor family is a subgroup of the nuclear receptor superfamily that also includes receptors for the active forms of vitamin A, vitamin D3, and thyroid hormones. Like all nuclear receptors, AR has a conserved modular structure consisting of a non-conserved amino-terminal domain (NTD), containing the intrinsic activation function 1, a highly conserved DNA-binding domain, and a conserved ligand-binding domain (LBD) that harbors the activation function 2. Each of these domains plays an important role in receptor function and signaling, either via intra- and inter-receptor interactions, interactions with specific DNA sequences, termed hormone response elements, or via functional interactions with domain-specific proteins, termed coregulators (coactivators and corepressors). Upon binding androgens, AR acquires a new conformational state, translocates to the nucleus, binds to androgen response elements, homodimerizes and recruits sequence-specific coregulatory factors and the basal transcription machinery. This set of events is required to activate gene transcription (expression). Gene transcription is a strictly modulated process that governs cell growth, cell homeostasis, cell function and cell death. Disruptions of AR transcriptional activity caused by receptor mutations and/or altered coregulator interactions are linked to a wide spectrum of androgen insensitivity syndromes, and to the pathogenesis of prostate cancer (CaP). The treatment of CaP usually involves androgen depletion therapy (ADT). ADT achieves significant clinical responses during the early stages of the disease. However, under the selective pressure of androgen withdrawal, androgen-dependent CaP can progress to an androgen-independent CaP. Androgen-independent CaP is invariably a more aggressive and untreatable form of the disease. Advancing our understanding of the molecular mechanisms behind the switch in androgen-dependency would improve our success of treating CaP and other AR related illnesses. This study evaluates how clinically identified AR mutations affect the receptor s transcriptional activity. We reveal that a potential molecular abnormality in androgen insensitivity syndrome and CaP patients is caused by disruptions of the important intra-receptor NTD/LBD interaction. We demonstrate that the same AR LBD mutations can also disrupt the recruitment of the p160 coactivator protein GRIP1. Our investigations reveal that 30% of patients with advanced, untreated local CaP have somatic mutations that may lead to increases in AR activity. We report that somatic mutations that activate AR may lead to early relapse in ADT. Our results demonstrate that the types of ADT a CaP patient receives may cause a clustering of mutations to a particular region of the receptor. Furthermore, the mutations that arise before and during ADT do not always result in a receptor that is more active, indicating that coregulator interactions play a pivotal role in the progression of androgen-independent CaP. To improve CaP therapy, it is necessary to identify critical coregulators of AR. We screened a HeLa cell cDNA library and identified small carboxyl-terminal domain phosphatase 2 (SCP2). SCP2 is a protein phosphatase that directly interacts with the AR NTD and represses AR activity. We demonstrated that reducing the endogenous cellular levels of SCP2 causes more AR to load on to the prostate specific antigen (PSA) gene promoter and enhancer regions. Additionally, under the same conditions, more RNA polymerase II was recruited to the PSA promoter region and overall there was an increase in androgen-dependent transcription of the PSA gene, revealing that SCP2 could play a role in the pathogenesis of CaP.

Molecular Background of Common Dyslipidemias

The leading cause of death in the Western world continues to be coronary heart disease (CHD). At the root of the disease process is dyslipidemia an aberration in the relevant amounts of circulating blood lipids. Cholesterol builds up in the arterial wall and following rupture of these plaques, myocardial infarction or stroke can occur. Heart disease runs in families and a number of hereditary forms are known. The leading cause of adult dyslipidemia presently however is overweight and obesity. This thesis work presents an investigation of the molecular genetics of common, hereditary dyslipidemia and the tightly related condition of obesity. Familial combined hyperlipidemia (FCHL) is the most common hereditary dyslipidemia in man with an estimated population prevalence of 1-6%. This complex disease is characterized by elevated levels of serum total cholesterol, triglycerides or both and is observed in about 20% of individuals with premature CHD. Our group identified the disease to be associated with genetic variation in the USF1 transcription factor gene. USF1 has a key role in regulating other genes that control lipid and glucose metabolism as well as the inflammatory response all central processes in the progression of atherosclerosis and CHD. The first two works of this thesis aimed at understanding how these USF1 variants result in increased disease risk. Among the many, non-coding single-nucleotide polymorphisms (SNPs) that associated with the disease, one was found to have a functional effect. The risk-enhancing allele of this SNP seems to eradicate the ability of the important hormone insulin to induce the expression of USF1 in peripheral tissues. The resultant changes in the expression of numerous USF1 target genes over time probably enhance and accelerate the atherogenic processes. Dyslipidemias often represent an outcome of obesity and in the final work of this thesis we wanted to address the metabolic pathways related to acquired obesity. It is recognized that active processes in adipose tissue play an important role in the development of dyslipidemia, insulin resistance and other pathological conditions associated with obesity. To minimize the confounding effects of genetic differences present in most human studies, we investigated a rare collection of identical twins that differed significantly in the amount of body fat. In the obese, but otherwise healthy young adults, several notable changes were observed. In addition to chronic inflammation, the adipose tissue of the obese co-twins was characterized by a marked (47%) decrease in amount of mitochondrial DNA (mtDNA) a change associated with mitochondrial dysfunction. The catabolism of branched chain amino acids (BCAAs) was identified as the most down-regulated process in the obese co-twins. A concordant increase in the serum level of these insulin secretagogues was identified. This hyperaminoacidemia may provide the feed-back signal from insulin resistant adipose tissue to the pancreas to ensure an appropriately augmented secretory response. The down regulation of BCAA catabolism correlated closely with liver fat accumulation and insulin. The single most up-regulated gene (5.9 fold) in the obese co-twins was osteopontin (SPP1) a cytokine involved in macrophage recruitment to adipose tissue. SPP1 is here implicated as an important player in the development of insulin resistance. These studies of exceptional study samples provide better understanding of the underlying pathology in common dyslipidemias and other obesity associated diseases important for future improvement of intervention strategies and treatments to combat atherosclerosis and coronary heart disease.