Genetic background of bipolar disorder and related cognitive impairments

Bipolar disorder (BP) is a complex psychiatric disorder characterized by episodes of mania and depression. BP affects approximately 1% of the world’s population and shows no difference in lifetime prevalence between males and females. BP arises from complex interactions among genetic, developmental and environmental factors, and it is likely that several predisposing genes are involved in BP. The genetic background of BP is still poorly understood, although intensive and long-lasting research has identified several chromosomal regions and genes involved in susceptibility to BP. This thesis work aims to identify the genetic variants that influence bipolar disorder in the Finnish population by candidate gene and genome-wide linkage analyses in families with many BP cases. In addition to diagnosis-based phenotypes, neuropsychological traits that can be seen as potential endophenotypes or intermediate traits for BP were analyzed. In the first part of the thesis, we examined the role of the allelic variants of the TSNAX/DISC1 gene cluster to psychotic and bipolar spectrum disorders and found association of distinct allelic haplotypes with these two groups of disorders. The haplotype at the 5’ end of the Disrupted-in-Schizophrenia-1 gene (DISC1) was over-transmitted to males with psychotic disorder (p = 0.008; for an extended haplotype p = 0.0007 with both genders), whereas haplotypes at the 3’ end of DISC1 associated with bipolar spectrum disorder (p = 0.0002; for an extended haplotype p = 0.0001). The variants of these haplotypes also showed association with different cognitive traits. The haplotypes at the 5’ end associated with perseverations and auditory attention, while the variants at the 3’ end associated with several cognitive traits including verbal fluency and psychomotor processing speed. Second, in our complete set of BP families with 723 individuals we studied six functional candidate genes from three distinct signalling systems: serotonin-related genes (SLC6A4 and TPH2), BDNF -related genes (BDNF, CREB1 and NTRK2) and one gene related to the inflammation and cytokine system (P2RX7). We replicated association of the functional variant Val66Met of BDNF with BP and better performance in retention. The variants at the 5’ end of SLC6A4 also showed some evidence of association among males (p = 0.004), but the widely studied functional variants did not yield any significant results. A protective four-variant haplotype on P2RX7 showed evidence of association with BP and executive functions: semantic and phonemic fluency (p = 0.006 and p = 0.0003, respectively). Third, we analyzed 23 bipolar families originating from the North-Eastern region of Finland. A genome-wide scan was performed using the 6K single nucleotide polymorphism (SNP) array. We identified susceptibility loci at chromosomes 7q31 with a LOD score of 3.20 and at 9p13.1 with a LOD score of 4.02. We followed up both linkage findings in the complete set of 179 Finnish bipolar families. The finding on chromosome 9p13 was supported (maximum LOD score of 3.02), but the susceptibility gene itself remains unclarified. In the fourth part of the thesis, we wanted to test the role of the allelic variants that have associated with bipolar disorder in recent genome-wide association studies (GWAS). We could confirm findings for the DFNB31, SORCS2, SCL39A3, and DGKH genes. The best signal in this study comes from DFNB31, which remained significant after multiple testing corrections. Two variants of SORCS2 were allelic replications and presented the same signal as the haplotype analysis. However, no association was detected with the PALB2 gene, which was the most significantly associated region in the previous GWAS. Our results indicate that BP is heterogeneous and its genetic background may accordingly vary in different populations. In order to fully understand the allelic heterogeneity that underlies common diseases such as BP, complete genome sequencing for many individuals with and without the disease is required. Identification of the specific risk variants will help us better understand the pathophysiology underlying BP and will lead to the development of treatments with specific biochemical targets. In addition, it will further facilitate the identification of environmental factors that alter risk, which will potentially provide improved occupational, social and psychological advice for individuals with high risk of BP.

Identification of loci for migraine with aura

Common migraine, i.e. migraine with (MA) or without aura (MO), is a chronic neurological disorder affecting about 10% of the Caucasian population. In MA, migraine headache is preceded by visual, sensoric and/or dysphasic reversible aura symptoms. Twin and family studies have suggested a multifactorial mode of inheritance for common migraine, and a stronger genetic component for MA than for MO. Since there is no biological or genetic marker to identify common migraine, aura symptoms provide a distinctive character to identify those suspected of suffering from migraine. The aim of this study was to identify MA susceptibility loci in well-phenotyped migraine samples with familial predisposition using different gene mapping methods. Genes coding for endothelin1 and its receptors EDNRA and ENDRB are potential candidate genes for cortical spreading depression (CSD), which is considered to be the underlying mechanism of migraine aura. The role of these genes in MA was studied in 850 Finnish migraine cases and 890 control individuals. Rare homozygous EDNRA SNPs showed nominal association with MA and with the age of onset trait (20 years). This result was also detected in the pooled analysis on 648 German MA cases and 651 control individuals when the test was adjusted for gender and sample origin. Evaluation of SNP genotyping reactions with two different DNA polymerase enzymes ensured that the genotype quality was high, and thus the discovered associations are considered reliable. The role of the 19p13 region was studied in a linkage analysis of 72 Finnish MA families. This region contains two migraine-associated genes: CACNA1A, which is associated with a predisposition to a rare Mendelian form of MA, familial hemiplegic migraine (FHM), and the insulin receptor gene (INSR) that is associated with common migraine. No evidence of linkage between the 19p13 and MA was detected. A novel visual aura locus was mapped to chromosome 9q21-q22 with significant evidence of linkage using a genome-wide linkage approach in 36 Finnish MA families. Five additional, potential loci were also detected. The 9q21-q22 region has previously been linked to occipitotemporal lobe epilepsy and MA, both of which involve prominent visual symptoms. Our result further supports a shared background for these episodic disorders.

Genetic Mechanisms Underlying Autism Spectrum Disorders

Autism is a childhood-onset developmental disorder characterized by deficits in reciprocal social interaction, verbal and non-verbal communication, and dependence on routines and rituals. It belongs to a spectrum of disorders (autism spectrum disorders, ASDs) which share core symptoms but show considerable variation in severity. The whole spectrum affects 0.6-0.7% of children worldwide, inducing a substantial public health burden and causing suffering to the affected families. Despite having a very high heritability, ASDs have shown exceptional genetic heterogeneity, which has complicated the identification of risk variants and left the etiology largely unknown. However, recent studies suggest that rare, family-specific factors contribute significantly to the genetic basis of ASDs. In this study, we investigated the role of DISC1 (Disrupted-in-schizophrenia-1) in ASDs, and identified association with markers and haplotypes previously associated with psychiatric phenotypes. We identified four polymorphic micro-RNA target sites in the 3 UTR of DISC1, and showed that hsa-miR-559 regulates DISC1 expression in vitro in an allele-specific manner. We also analyzed an extended autism pedigree with genealogical roots in Central Finland reaching back to the 17th century. To take advantage of the beneficial characteristics of population isolates to gene mapping and reduced genetic heterogeneity observed in distantly related individuals, we performed a microsatellite-based genome-wide screen for linkage and linkage disequilibrium in this pedigree. We identified a putative autism susceptibility locus on chromosome 19p13.3 and obtained further support for previously reported loci at 1q23 and 15q11-q13. To follow-up these findings, we extended our study sample from the same sub-isolate and initiated a genome-wide analysis of homozygosity and allelic sharing using high-density SNP markers. We identified a small number of haplotypes shared by different subsets of the genealogically connected cases, along with convergent biological pathways from SNP and gene expression data, which highlighted axon guidance molecules in the pathogenesis of ASDs. In conclusion, the results obtained in this thesis show that multiple distinct genetic variants are responsible for the ASD phenotype even within single pedigrees from an isolated population. We suggest that targeted resequencing of the shared haplotypes, linkage regions, and other susceptibility loci is essential to identify the causal variants. We also report a possible micro-RNA mediated regulatory mechanism, which might partially explain the wide-range neurobiological effects of the DISC1 gene.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Identification of genetic susceptibility loci for migraine

Migraine is the common cause of chronic episodic headache, affecting 12%-15% of the Caucasian population (41 million Europeans and some half a million Finns), and causes considerable loss of quality of life to its sufferers, as well as being linked to increased risk for a wide range of conditions, from depression to stroke. Migraine is the 19th most severe disease in terms of disability-adjusted life years, and 9th among women. It is characterized by attacks of headache accompanied by sensitivity to external stimuli lasting 4-72 hours, and in a third of cases by neurological aura symptoms, such as loss of vision, speech or muscle function. The underlying pathophysiology, including what triggers migraine attacks and why they occur in the first place, is largely unknown. The aim of this study was to identify genetic factors associated with the hereditary susceptibility to migraine, in order to gain a better understanding of migraine mechanisms. In this thesis, we report the results of genetic linkage and association analyses on a Finnish migraine patient collection as well as migraineurs from Australia, Denmark, Germany, Iceland and the Netherlands. Altogether we studied genetic information of nearly 7,000 migraine patients and over 50,000 population-matched controls. We also developed a new migraine analysis method called the trait component analysis, which is based on individual patient responses instead of the clinical diagnosis. Using this method, we detected a number of new genetic loci for migraine, including on chromosome 17p13 (HLOD 4.65) and 10q22-q23 (female-specific HLOD 7.68) with significant evidence of linkage, along with five other loci (2p12, 8q12, 4q28-q31, 18q12-q22, and Xp22) detected with suggestive evidence of linkage. The 10q22-q23 locus was the first genetic finding in migraine to show linkage to the same locus and markers in multiple populations, with consistent detection in six different scans. Traditionally, ion channels have been thought to play a role in migraine susceptibility, but we were able to exclude any significant role for common variants in a candidate gene study of 155 ion transport genes. This was followed up by the first genome-wide association study in migraine, conducted on 2,748 migraine patients and 10,747 matched controls followed by a replication in 3,209 patients and 40,062 controls. In this study, we found interesting results with genome-wide significance, providing targets for future genetic and functional studies. Overall, we found several promising genetic loci for migraine providing a promising base for future studies in migraine.

Identifying Genetic Risk Factors in Canine Autoimmune Disorders

Autoimmune diseases are more common in dogs than in humans and are already threatening the future of some highly predisposed dog breeds. Susceptibility to autoimmune diseases is controlled by environmental and genetic factors, especially the major histocompatibility complex (MHC) gene region. Dogs show a similar physiology, disease presentation and clinical response as humans, making them an excellent disease model for autoimmune diseases common to both species. The genetic background of canine autoimmune disorders is largely unknown, but recent annotation of the dog genome and subsequent development of new genomic tools offer a unique opportunity to map novel autoimmune genes in various breeds. Many autoimmune disorders show breed-specific enrichment, supporting a strong genetic background. Furthermore, the presence of hundreds of breeds as genetic isolates facilitates gene mapping in complex autoimmune disorders. Identification of novel predisposing genes establishes breeds as models and may reveal novel candidate genes for the corresponding human disorders. Genetic studies will eventually shed light on common biological functions and interactions between genes and the environment. This study aimed to identify genetic risk factors in various autoimmune disorders, including systemic lupus erythematosus (SLE)-related diseases, comprising immune-mediated rheumatic disease (IMRD) and steroid-responsive meningitis arteritis (SMRA) as well as Addison s disease (AD) in Nova Scotia Duck Tolling Retrievers (NSDTRs) and chronic superficial keratitis (CSK) in German Shepherd dogs (GSDs). We used two different approaches to identify genetic risk factors. Firstly, a candidate gene approach was applied to test the potential association of MHC class II, also known as a dog leukocyte antigen (DLA) in canine species. Secondly, a genome-wide association study (GWAS) was performed to identify novel risk loci for SLE-related disease and AD in NSDTRs. We identified DLA risk haplotypes for an IMRD subphenotype of SLE-related disease, AD and CSK, but not in SMRA, and show that the MHC class II gene region is a major genetic risk factor in canine autoimmune diseases. An elevated risk was found for IMRD in dogs that carried the DLA-DRB1*00601/DQA1*005011/DQB1*02001 haplotype (OR = 2.0, 99% CI = 1.03-3.95, p = 0.01) and for ANA-positive IMRD dogs (OR = 2.3, 99% CI = 1.07-5.04, p-value 0.007). We also found that DLA-DRB1*01502/DQA*00601/DQB1*02301 haplotype was significantly associated with AD in NSDTRs (OR = 2.1, CI = 1.0-4.4, P = 0.044) and the DLA-DRB1*01501/DQA1*00601/DQB1*00301 haplotype with the CSK in GSDs (OR=2.67, CI=1.17-6.44, p= 0.02). In addition, we found that homozygosity for the risk haplotype increases the risk for each disease phenotype and that an overall homozygosity for the DLA region predisposes to CSK and AD. Our results have enabled the development of genetic tests to improve breeding practices by avoiding the production of puppies homozygous for risk haplotypes. We also performed the first successful GWAS for a complex disease in dogs. With less than 100 cases and 100 controls, we identified five risk loci for SLE-related disease and AD and found strong candidate genes involved in a novel T-cell activation pathway. We show that an inbred dog population has fewer risk factors, but each of them has a stronger genetic risk. Ongoing studies aim to identify the causative mutations and bring new knowledge to help diagnostics, treatment and understanding of the aetiology of SLE-related diseases.

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.

Genetic and epigenetic mechanisms of tumor predisposition in hereditary non-polyposis colorectal carcinoma and sporadic cancers

Individuals with inherited deficiency in DNA mismatch repair(MMR) (Lynch syndrome) LS are predisposed to different cancers in a non-random fashion. Endometrial cancer (EC) is the most common extracolonic malignancy in LS. LS represents the best characterized form of hereditary nonpolyposis colorectal carcinoma (HNPCC). Other forms of familial non-polyposis colon cancer exist, including familial colorectal cancer type X (FCCX). This syndrome resembles LS, but MMR gene defects are excluded and the predisposition genes are unknown so far. To address why different organs are differently susceptible to cancer development, we examined molecular similarities and differences in selected cancers whose frequency varies in LS individuals. Tumors that are common (colorectal, endometrial, gastric) and less common (brain, urological) in LS were characterized for MMR protein expression, microsatellite instability (MSI), and by altered DNA methylation. We also studied samples of histologically normal endometrium, endometrial hyperplasia,and cancer for molecular alterations to identify potential markers that could predict malignant transformation in LS and sporadic cases. Our results suggest that brain and kidney tumors follow a different pathway for cancer development than the most common LS related cancers.Our results suggest also that MMR defects are detectable in endometrial tissues from a proportion of LS mutation carriers prior to endometrial cancer development. Traditionally (complex) atypical hyperplasia has been considered critical for progression to malignancy. Our results suggest that complex hyperplasia without atypia is equally important as a precursor lesion of malignancy. Tumor profiles from Egypt were compared with colorectal tumors from Finland to evaluate if there are differences specific to the ethnic origin (East vs.West). Results showed for the first time a distinct genetic and epigenetic signature in the Egyptian CRC marked by high methylation of microsatellite stable tumors associated with advanced stage, and low frequency of Wnt signaling activation, suggesting a novel pathway. DNA samples from FCCX families were studied with genome wide linkage analysis using microsatellite markers. Selected genes from the linked areas were tested for possible mutations that could explain predisposition to a large number of colon adenomas and carcinomas seen in these families. Based on the results from the linkage analysis, a number of areas with tentative linkage were identified in family 20. We narrowed down these areas by additional microsatellite markers to found a mutation in the BMPR1A gene. Sequencing of an additional 17 FCCX families resulted in a BMPR1A mutation frequency of 2/18 families (11%). Clarification of the mechanisms of the differential tumor susceptibility in LS increases the understanding of gene and organ specific targets of MMR deficiency. While it is generally accepted that widespread MMR deficiency and consequent microsatellite instability (MSI) drives tumorigenesis in LS, the timing of molecular alterations is controversial. In particular, it is important to know that alterations may occur several years before cancer formation, at stages that are still histologically regarded as normal. Identification of molecular markers that could predict the risk of malignant transformation may be used to improve surveillance and cancer prevention in genetically predisposed individuals. Significant fractions of families with colorectal and/or endometrial cancer presently lack molecular definition altogether. Our findings expand the phenotypic spectrum of BMPR1A mutations and, for the first time, link FCCX families to the germline mutation of a specific gene. In particular, our observations encourage screening of additional families with FCCX for BMPR1A mutation, which is necessary in obtaining a reliable estimate of the share of BMPR1A-associated cases among all FCCX families worldwide. Clinically, the identification of predisposing mutations enables targeted cancer prevention in proven mutation carriers and thereby reduces cancer morbidity and mortality in the respective families.

Functions of human papillomavirus E5 oncogene in epithelial cells and in the onset of cervical cancer

Human papillomaviruses (HPVs) are the causal agents of cervical cancer, which is the second most common cancer among women worldwide. Cellular transformation and carcinogenesis depend on the activities of viral E5, E6 and E7 proteins. Alterations in cell-cell contacts and in communication between epithelial cells take place during cervical carcinogenesis, leading to changes in cell morphology, increased cell motility and finally invasion. The aim of this thesis was to study genome-wide effects of the HPV type 16 (HPV-16) E5 protein on the expression of host cell messenger RNAs (mRNAs) and microRNAs by applying microarray technology. The results showed that the HPV-16 E5 protein alters several cellular pathways involved in cellular adhesion, motility and proliferation as well as in the extracellular matrix. The E5 protein was observed to enhance wound healing of epithelial cell monolayers by increasing cell motility in vivo. HPV-16 E5-induced alterations in the expression of cellular microRNAs and their target genes seem to favour increased proliferation and tumorigenesis. E5 was also shown to affect the expression of adherens junction proteins in HaCaT epithelial keratinocytes. In addition, a study of a membrane cytoskeletal cross-linker protein, ezrin, revealed that when activated, it localizes to adherens junctions. The results suggest that ezrin distribution to forming adherens junctions is due to Rac1 activity in epithelial cells. These studies reveal for the first time the holistic effects of HPV-16 E5 protein in promoting precancerous events in epithelial cells. The results contribute to identifyinging novel markers for cervical precancerous stages and to predicting disease behaviour.

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.

Studies on Causes of Multiple Sclerosis : From Genes to Transcriptome

Multiple sclerosis (MS) is the most common cause of neurological disability in young adults, affecting more than two million people worldwide. It manifests as a chronic inflammation in the central nervous system (CNS) and causes demyelination and neurodegeneration. Depending on the location of the demyelinated plaques and axonal loss, a variety of symptoms can be observed including deficits in vision, coordination, balance and movement. With a typical age of onset at 20-40 years, the social and economic impacts of MS on lives of the patients and their families are considerable. Unfortunately the current treatments are relatively inefficient and the development of more effective treatments has been impeded by our limited understanding of the causes and pathogenesis of MS. Risk of MS is higher in biological relatives of MS patients than in the general population. Twin and adoption studies have shown that familial clustering of MS is explained by shared genetic factors rather than by shared familial environment. While the involvement of the human leukocyte antigen (HLA) genes was first discovered four decades ago, additional genetic risk factors have only recently been identified through genome-wide association studies (GWAS). Current evidence suggests that MS is a highly polygenic disease with perhaps hundreds of common variants with relatively modest effects contributing to susceptibility. Despite extensive research, the majority of these risk factors still remain to be identified. In this thesis the aim was to identify novel genes and pathways involved in MS. Using genome-wide microarray technology, gene expression levels in peripheral blood mononuclear cells (PBMC) from 12 MS patients and 15 controls were profiled and more than 600 genes with altered expression in MS were identified. Three of five selected findings, DEFA1A3, LILRA4 and TNFRSF25, were successfully replicated in an independent sample. Increased expression of DEFA1A3 in MS is a particularly interesting observation, because its elevated levels have previously been reported also in several other autoimmune diseases. A systematic review of seven microarray studies was then performed leading to identification of 229 genes, in which either decreased or increased expression in MS had been reported in at least two studies. In general there was relatively little overlap across the experiments: 11 of the 229 genes had been reported in three studies and only HSPA1A in four studies. Nevertheless, these 229 genes were associated with several immunological pathways including interleukin pathways related to type 2 and type 17 helper T cells and regulatory T cells. However, whether these pathways are involved in causing MS or related to secondary processes activated after disease onset remains to be investigated. The 229 genes were also compared with loci identified in published MS GWASs. Single nucleotide polymorphisms (SNP) in 17 of the 229 loci had been reported to be associated with MS with P-value less than 0.0001 including variants in CXCR4 and SAPS2, which were the only loci where evidence for correlation between the associated variant and gene expression was found. The CXCR4 variant was further tested for association with MS in a large case-control sample and the previously reported suggestive association was replicated (P-value is 0.0004). Finally, common genetic variants in candidate genes, which had been selected on the basis of showing association with other autoimmune diseases (MYO9B) or showing differential expression in MS in our study (DEFA1A3, LILRA4 and TNFRSF25), were tested for association with MS, but no evidence of association was found. In conclusion, through a systematic review of genome-wide expression studies in MS we have identified several promising candidate genes and pathways for future studies. In addition, we have replicated a previously suggested association of a SNP variant upstream of CXCR4 with MS. Keywords: autoimmune disease, common variant, CXCR4, DEFA1A3, HSPA1A,gene expression, genetic association, GWAS, MS, multiple sclerosis, systematic review

LIN28B, LIN28A, KISS1, and KISS1R in idiopathic central precocious puberty

Abstract Background Pubertal timing is a strongly heritable trait, but no single puberty gene has been identified. Thus, the genetic background of idiopathic central precocious puberty (ICPP) is poorly understood. Overall, the genetic modulation of pubertal onset most likely arises from the additive effect of multiple genes, but also monogenic causes of ICPP probably exist, as cases of familial ICPP have been reported. Mutations in KISS1 and KISSR, coding for kisspeptin and its receptor, involved in GnRH secretion and puberty onset, have been suggested causative for monogenic ICPP. Variation in LIN28B was associated with timing of puberty in genome-wide association (GWA) studies. LIN28B is a human ortholog of the gene that controls, through microRNAs, developmental timing in C. elegans. In addition, Lin28a transgenic mice manifest the puberty phenotypes identified in the human GWAS. Thus, both LIN28B and LIN28A may have a role in pubertal development and are good candidate genes for monogenic ICPP. Methods Thirty girls with ICPP were included in the study. ICPP was defined by pubertal onset before 8 yrs of age, and a pubertal LH response to GnRH testing. The coding regions of LIN28B, LIN28A, KISS1, and KISS1R were sequenced. The missense change in LIN28B was also screened in 132 control subjects. Results No rare variants were detected in KISS1 or KISS1R in the 30 subjects with ICPP. In LIN28B, one missense change, His199Arg, was found in one subject with ICPP. However, this variant was also detected in one of the 132 controls. No variation in LIN28A was found. Conclusions We did not find any evidence that mutations in LIN28B or LIN28A would underlie ICPP. In addition, we confirmed that mutations in KISS1 and KISS1R are not a common cause for ICPP.