Uniparental DNA markers and forensic genetics in Finland

Over the years, a wide range of methods to verify identity have been developed. Molecular markers have been used for identification since the 1920s, commencing with blood types and culminating with the advent of DNA techniques in the 1980s. Identification is required by authorities in many occasions, e.g. in disputed paternity cases, identification of deceased, or crime investigation. To clarify maternal and paternal lineages, uniparental DNA markers in mtDNA and Y-chromosome can be utilized. These markers have several advantages: male specific Y-chromosome can be used to identify a male from a mixture of male and female cells, e.g. in rape cases. MtDNA is durable and has a high copy number, allowing analyses even from old or degraded samples. However, both markers are lineage-specific, not individualizing, and susceptible to genetic drift. Prior to the application of any DNA marker in forensic casework, it is of utmost importance to investigate its qualities and peculiarities in the target population. Earlier studies on the Finnish population have shown reduced variation in the Y-chromosome, but in mtDNA results have been ambiguous. The obtained results confirmed the low diversity in Y-chromosome in Finland. Detailed population analysis revealed large regional differences, and extremely reduced diversity especially in East Finland. Analysis of the qualities affecting Y-chromosomal short tandem repeat (Y-STR) variation and mutation frequencies, and search of new polymorphic markers resulted a set of Y-STRs with especially high diversity in Finland. Contrary to Y-chromosome, neither reduced diversity nor regional differences were found in mtDNA within Finland. In fact, mtDNA diversity was found similar to other European populations. The revealed peculiarities in the uniparental markers are a legacy of the Finnish population history. The obtained results challenge the traditional explanation which emphasizes relatively recent founder effects creating the observed east-west patterns. Uniparentally inherited markers, both mtDNA and Y-chromosome, are applicable for identification purposes in Finland. By adjusting the analysed Y marker set to meet the characteristics of Finnish population, Y-chromosomal diversity increases and the regional differentiation decreases, resulting increase in discrimination power and thus usefulness of Y-chromosomal analysis in forensic casework.

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 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.

Studies of the Human Transcriptome

Gene expression is one of the most critical factors influencing the phenotype of a cell. As a result of several technological advances, measuring gene expression levels has become one of the most common molecular biological measurements to study the behaviour of cells. The scientific community has produced enormous and constantly increasing collection of gene expression data from various human cells both from healthy and pathological conditions. However, while each of these studies is informative and enlighting in its own context and research setup, diverging methods and terminologies make it very challenging to integrate existing gene expression data to a more comprehensive view of human transcriptome function. On the other hand, bioinformatic science advances only through data integration and synthesis. The aim of this study was to develop biological and mathematical methods to overcome these challenges and to construct an integrated database of human transcriptome as well as to demonstrate its usage. Methods developed in this study can be divided in two distinct parts. First, the biological and medical annotation of the existing gene expression measurements needed to be encoded by systematic vocabularies. There was no single existing biomedical ontology or vocabulary suitable for this purpose. Thus, new annotation terminology was developed as a part of this work. Second part was to develop mathematical methods correcting the noise and systematic differences/errors in the data caused by various array generations. Additionally, there was a need to develop suitable computational methods for sample collection and archiving, unique sample identification, database structures, data retrieval and visualization. Bioinformatic methods were developed to analyze gene expression levels and putative functional associations of human genes by using the integrated gene expression data. Also a method to interpret individual gene expression profiles across all the healthy and pathological tissues of the reference database was developed. As a result of this work 9783 human gene expression samples measured by Affymetrix microarrays were integrated to form a unique human transcriptome resource GeneSapiens. This makes it possible to analyse expression levels of 17330 genes across 175 types of healthy and pathological human tissues. Application of this resource to interpret individual gene expression measurements allowed identification of tissue of origin with 92.0% accuracy among 44 healthy tissue types. Systematic analysis of transcriptional activity levels of 459 kinase genes was performed across 44 healthy and 55 pathological tissue types and a genome wide analysis of kinase gene co-expression networks was done. This analysis revealed biologically and medically interesting data on putative kinase gene functions in health and disease. Finally, we developed a method for alignment of gene expression profiles (AGEP) to perform analysis for individual patient samples to pinpoint gene- and pathway-specific changes in the test sample in relation to the reference transcriptome database. We also showed how large-scale gene expression data resources can be used to quantitatively characterize changes in the transcriptomic program of differentiating stem cells. Taken together, these studies indicate the power of systematic bioinformatic analyses to infer biological and medical insights from existing published datasets as well as to facilitate the interpretation of new molecular profiling data from individual patients.

Population genetics of the invasive plant species Impatiens glandulifera in Southern Finland

Biological invasions affect biodiversity worldwide, and, consequently, the invaded ecosystems may suffer from significant losses in economic and cultural values. Impatiens glandulifera Royle (Balsaminaceae) is an invasive annual herb, native to the western Himalayas and introduced into Europe in the 19th century as a garden ornamental plant. The massive invasion of I. glandulifera is due to its high reproductive output, rapid growth and its ability to outcompete native species. In Finland, the first observations regarding the presence of I. glandulifera date from the year 1947, and today it is considered a serious problem in riparian habitats. The aim of this master’s thesis research is to reveal the population genetic structure of I. glandulifera in Finland and to find out whether there have been one or multiple invasions in Finland. The study focuses on investigating the origin of I. glandulifera in Southern Finland, by comparing plant samples from the Helsinki region with those from its native region and other regions of invasion. Samples from four populations in Helsinki and from the United Kingdom, Canada, India and Pakistan were collected and genotyped using 11 microsatellite markers. The genetic analyses were evaluated using the programs Arlequin and Structure. The results of the genetic analyses suggested that I. glandulifera has been introduced to Finland more than once. Multiple introductions are supported by the higher level of genetic diversity detected within and among Finnish populations than would be expected for a single introduction. Results of the Bayesian Structure analysis divided the four Finnish populations into four clusters. This geographical structure was further supported by pairwise Fst values among populations. The causes and potential consequences of such multiple introductions of I. glandulifera in Finland and further perspectives are discussed.

Protein Kinase C and Anaplastic Lymphoma Kinase Targeted Compounds

The protein kinases (PKs) belong to the largest single family of enzymes, phosphotransferases, which catalyze the phosphorylation of other enzymes and proteins and function primarily in signal transduction. Consequently, PKs regulate cell mechanisms such as growth, differentiation, and proliferation. Dysfunction of these cellular mechanisms may lead to cancer, a major predicament in health care. Even though there is a range of clinically available cancer-fighting drugs, increasing number of cancer cases and setbacks such as drug resistance, constantly keep cancer research active. At the commencement of this study an isophthalic acid derivative had been suggested to bind to the regulatory domain of protein kinase C (PKC). In order to investigate the biological effects and structure-activity relationships (SARs) of this new chemical entity, a library of compounds was synthesized. The best compounds induced apoptosis in human leukemia HL-60 cells and were not cytotoxic in Swiss 3T3 fibroblasts. In addition, the best apoptosis inducers were neither cytotoxic nor mutagenic. Furthermore, results from binding affinity assays of PKC isoforms revealed the pharmacophores of these isophthalic acid derivatives. The best inhibition constants of the tested compounds were measured to 210 nM for PKCα and to 530 nM for PKCδ. Among natural compounds targeting the regulatory domain of PKC, the target of bistramide A has been a matter of debate. It was initially found to activate PKCδ; however, actin was recently reported as the main target. In order to clarify and to further study the biological effects of bistramide A, the total syntheses of the natural compound and two isomers were performed. Biological assays of the compounds revealed accumulation of 4n polyploid cells as the primary mode of action and the compounds showed similar overall antiproliferative activities. However, each compound showed a distinct distribution of antimitotic effect presumably via actin binding, proapoptotic effect presumably via PKCδ, and pro-differentiation effect as evidenced by CD11b expression. Furthermore, it was shown that the antimitotic and proapoptotic effects of bistramide A were not secondary effects of actin binding but independent effects. The third aim in this study was to synthesize a library of a new class of urea-based type II inhibitors targeted at the kinase domain of anaplastic lymphoma kinase (ALK). The best compounds in this library showed IC50 values as low as 390 nM for ALK while the initial low cellular activities were successfully increased even by more than 70 times for NPM-ALK- positive BaF3 cells. More importantly, selective antiproliferative activity on ALK-positive cell lines was achieved; while the best compound affected the BaF3 and SU-DHL-1 cells with IC50 values of 0.5 and 0.8 μM, respectively, they were less toxic to the NPM-ALK-negative human leukemic cells U937 (IC50 = 3.2 μM) and BaF3 parental cells (IC50 = 5.4 μM). Furthermore, SAR studies of the synthesized compounds revealed functional groups and positions of the scaffold, which enhanced the enzymatic and cellular activities.

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.