Molecular genetics of Meckel syndrome : Ciliary genes are defective in MKS

Meckel syndrome (MKS, MIM 249000) is a severe developmental disorder that leads to death already in utero or shortly after birth. MKS diagnosis can be established by a careful ultrasound examination already at 11-14 weeks of gestation. The main features of MKS are occipital meningoencephalocele, cystic kidney dysplasia and fibrotic changes of the liver. In addition, polydactyly is frequently reported in the cases. The aim of the study was to characterize the molecular and functional defects in MKS. In this study we were able to identify two major MKS mutations in Finnish population, which cover over 90% of the cases. The first mutation is a 29 bp intronic deletion in the MKS1 gene (c.1483-7_35del) that is found in 70% of the families and the second is a C>T substitution in the coding region of CC2D2A (c.1762C>T), that is found in 20% of the MKS families. Both of these mutations result in abnormal splicing. The discovery of the disease genes has revealed that MKS is caused by primary cilia dysfunction. MKS1 gene has a conserved B9 domain, and it is found in the predicted ciliary proteome. CC2D2A protein is also found in the predicted ciliary proteome and it has a Ca2+ binding domain. The number of genes behind MKS has increased rapidly in the past years and to date, mutations have been identified in five genes (MKS1, TMEM67/MKS3, CEP290/MKS4, RPGRIP1L/MKS5 and CC2D2A/MKS6). Identification of the disease genes mutations has also revealed that MKS is an allelic disorder with other syndromes with overlapping phenotypes. Disorders that are caused by primary cilia dysfunction are collectively known as ciliopathies. Sequence analysis of all the known MKS genes in Finnish and non-Finnish families available to us, where the mutation was still unknown, revealed mutations in 14 out of the 30 families included in the study. When we collected all the reported mutations in MKS genes in different syndromes we could see that there was clearly a genotype-syndrome correlation between the mutations and the syndromes, since the same pair of mutations has never been reported in different syndromes. The basic molecular events behind MKS will not only give us information of this syndrome, but also significant novel information on early fetal development in general.

The Gerbera cDNA Microarray: A Tool for Large-Scale Identification of Genes Involved in Flower Development

During the past ten years, large-scale transcript analysis using microarrays has become a powerful tool to identify and predict functions for new genes. It allows simultaneous monitoring of the expression of thousands of genes and has become a routinely used tool in laboratories worldwide. Microarray analysis will, together with other functional genomics tools, take us closer to understanding the functions of all genes in genomes of living organisms. Flower development is a genetically regulated process which has mostly been studied in the traditional model species Arabidopsis thaliana, Antirrhinum majus and Petunia hybrida. The molecular mechanisms behind flower development in them are partly applicable in other plant systems. However, not all biological phenomena can be approached with just a few model systems. In order to understand and apply the knowledge to ecologically and economically important plants, other species also need to be studied. Sequencing of 17 000 ESTs from nine different cDNA libraries of the ornamental plant Gerbera hybrida made it possible to construct a cDNA microarray with 9000 probes. The probes of the microarray represent all different ESTs in the database. From the gerbera ESTs 20% were unique to gerbera while 373 were specific to the Asteraceae family of flowering plants. Gerbera has composite inflorescences with three different types of flowers that vary from each other morphologically. The marginal ray flowers are large, often pigmented and female, while the central disc flowers are smaller and more radially symmetrical perfect flowers. Intermediate trans flowers are similar to ray flowers but smaller in size. This feature together with the molecular tools applied to gerbera, make gerbera a unique system in comparison to the common model plants with only a single kind of flowers in their inflorescence. In the first part of this thesis, conditions for gerbera microarray analysis were optimised including experimental design, sample preparation and hybridization, as well as data analysis and verification. Moreover, in the first study, the flower and flower organ-specific genes were identified. After the reliability and reproducibility of the method were confirmed, the microarrays were utilized to investigate transcriptional differences between ray and disc flowers. This study revealed novel information about the morphological development as well as the transcriptional regulation of early stages of development in various flower types of gerbera. The most interesting finding was differential expression of MADS-box genes, suggesting the existence of flower type-specific regulatory complexes in the specification of different types of flowers. The gerbera microarray was further used to profile changes in expression during petal development. Gerbera ray flower petals are large, which makes them an ideal model to study organogenesis. Six different stages were compared and specifically analysed. Expression profiles of genes related to cell structure and growth implied that during stage two, cells divide, a process which is marked by expression of histones, cyclins and tubulins. Stage 4 was found to be a transition stage between cell division and expansion and by stage 6 cells had stopped division and instead underwent expansion. Interestingly, at the last analysed stage, stage 9, when cells did not grow any more, the highest number of upregulated genes was detected. The gerbera microarray is a fully-functioning tool for large-scale studies of flower development and correlation with real-time RT-PCR results show that it is also highly sensitive and reliable. Gene expression data presented here will be a source for gene expression mining or marker gene discovery in the future studies that will be performed in the Gerbera Laboratory. The publicly available data will also serve the plant research community world-wide.

Polymorphic low penetrance genes and breast cancer : The role of genes involved in metabolism of xenobiotics, estrogens and reactive oxygen species

Various endogenous and exogenous factors have been reported to increase the risk of breast cancer. Many of those are related to prolonged lifetime exposure to estrogens. Furthermore, a positive family history of breast cancer and certain benign breast diseases are known to increase the risk of breast cancer. The role of lifestyle factors, such as use of alcohol and smoking has been an area of intensive study. Alcohol has been found to increase the risk of breast cancer, whereas the role of smoking has remained obscure. A multitude of enzymes are involved in the metabolism of estrogens and xenobiotics including the carcinogens found in tobacco smoke. Many of the metabolic enzymes exhibit genetic polymorphisms that can lead to inter-individual differences in their abilities to modify hazardous substrates. Therefore, in presence of a given chemical exposure, one subgroup of women may be more susceptible to breast carcinogenesis, since they carry unfavourable forms of the polymorphic genes involved in the metabolism of the chemical. In this work, polymorphic genes encoding for cytochrome P450 (CYP) 1A1 and 1B1, N-acetyl transferase 2 (NAT2), sulfotransferase 1A1 (SULT1A1), manganese superoxide dismutase (MnSOD) and vitamin D receptor (VDR) were investigated in relation to breast cancer susceptibility in a Finnish population. CYP1A1, CYP1B1 and SULT1A1 are involved in the metabolism of both estrogens and xenobiotics, whereas NAT2 is involved only in the latter. MnSOD is an antioxidant enzyme protecting cells from oxidative damage. VDR, in turn, mediates the effects of the active form of vitamin D (1,25(OH)2D3, calcitriol) on maintenance of calcium homeostasis and it has anti-proliferative effects in many cancer cells. A 1.3-fold (95% CIs 1.01-1.73) increased risk of breast cancer was seen among women who carried the NAT2 slow acetylator genotype and a 1.5-fold (95% CI 1.1-2.0) risk was found in women with a MnSOD variant A allele containing genotypes compared to women with the NAT2 rapid acetylator genotype or to those with the MnSOD VV genotype, respectively. Instead, women with the VDR a allele containing genotypes were found to be at a decreased risk for breast cancer (OR 0.73; 95% CI 0.54-0.98) compared to women with the AA genotype. No significant overall associations were found between SULT1A1 or CYP genotypes and breast cancer risk, whereas a combination of the CYP1B1 432Val allele containing genotypes with the NAT2 slow acetylator genotypes posed a 1.5-fold (95% CI 1.03-2.24) increased risk. Moreover, NAT2 slow acetylator genotype was found to be confined to women with an advanced stage of breast cancer (stages III and IV). Further evidence for the association of xenobiotic metabolising genes with breast cancer risk was found when active smoking was taken into account. Women who smoked less than 10 cigarettes/day and carried at least one CYP1B1 432Val variant allele, were at 3.1-fold (95% CI 1.32-7.12) risk of breast cancer compared to women who smoked the same amount but did not carry the variant allele. Furthermore, the risk was significantly increased with increasing number of the CYP1B1 432Val alleles (p for trend 0.005). In addition, women who smoked less than 5 pack-years and carried the NAT2 slow acetylator genotype were at a 2.6-fold (95% CI 1.01-6.48) increased risk of breast cancer compared to women who smoked the same amount but carried the NAT2 rapid acetylator genotype. Furthermore, the combination of the CYP1B1 432Val allele and the NAT2 slow acetylator genotype increased the risk of breast cancer by 2.5-fold (95% CI 1.11-5.45) among ever smokers. Instead, the MnSOD A allele was found to be a risk factor among postmenopausal long-term smokers (>15 years of smoking) (OR 5.1; 95% CI 1.4-18.4) or among postmenopausal women who had smoked more than 10 cigarettes/day (OR 5.5; 95% CI 1.3-23.4) compared to women who had similar smoking habits but carried the MnSOD V/V genotype. Similarly, within subgroups of postmenopausal women who were using oral contraceptives, hormone replacement therapy or alcohol, women carrying the MnSOD A allele genotypes seemed to be at increased risk of breast cancer compared to women with the MnSOD V/V genotype. A positive family history of breast cancer and high parity were shown to be inversely associated with breast cancer risk among women carrying the VDR ApaI a allele or among premenopausal women carrying the SULT1A1*2 allele, respectively.

Molecular Alterations in Asbestos-Related Lung Cancer

Background: Asbestos is a well known cancer-causing mineral fibre, which has a synergistic effect on lung cancer risk in combination with tobacco smoking. Several in vitro and in vivo experiments have demonstrated that asbestos can evoke chromosomal damage and cause alterations as well as gene expression changes. Lung tumours, in general, have very complex karyotypes with several recurrently gained and lost chromosomal regions and this has made it difficult to identify specific molecular changes related primarily to asbestos exposure. The main aim of these studies has been to characterize asbestos-related lung cancer at a molecular level. Methods: Samples from asbestos-exposed and non-exposed lung cancer patients were studied using array comparative genomic hybridization (aCGH) and fluorescent in situ hybridization (FISH) to detect copy number alterations (CNA) as well as microsatellite analysis to detect allelic imbalance (AI). In addition, asbestos-exposed cell lines were studied using gene expression microarrays. Results: Eighteen chromosomal regions showing differential copy number in the lung tumours of asbestos-exposed patients compared to those of non-exposed patients were identified. The most significant differences were detected at 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2 and 19p13.1-p13.3 (p<0.005). The alterations at 2p and 9q were validated and characterized in detail using AI and FISH analysis in a larger study population. Furthermore, in vitro studies were performed to examine the early gene expression changes induced by asbestos in three different lung cell lines. The results revealed specific asbestos-associated gene expression profiles and biological processes as well as chromosomal regions enriched with genes believed to contribute to the common asbestos-related responses in the cell lines. Interestingly, the most significant region enriched with asbestos-response genes was identified at 2p22, close to the previously identified region showing asbestos-related CNA in lung tumours. Additionally, in this thesis, the dysregulated biological processes (Gene Ontology terms) detected in the cell line experiment were compared to dysregulated processes identified in patient samples in a later study (Ruosaari et al., 2008a). Commonly affected processes such as those related to protein ubiquitination, ion transport and surprisingly sensory perception of smell were identified. Conclusions: The identification of specific CNA and dysregulated biological processes shed some light on the underlying genes acting as mediators in asbestos-related lung carcinogenesis. It is postulated that the combination of several asbestos-specific molecular alterations could be used to develop a diagnostic method for the identification of asbestos-related lung cancer.

Genetic aspects of outcome in kidney transplantation: cytokine and thrombosis associated candidate genes and gene expression biomarkers

Kidney transplantation (Tx) is the treatment of choice for end stage renal disease. Immunosuppressive medications are given to prevent an immunological rejection of the transplant. However, immunosuppressive drugs increase e.g. the risk of infection, cancer or nephrotoxicity. A major genetic contributors to immunological acceptance of the graft are human leukocyte antigen (HLA) genes. Also other non-HLA gene polymorphisms may predict the future risk of complications before Tx, possibly enabling individualised immunotherapy. Graft function after Tx is monitored using non-specific clinical symptoms and laboratory markers. The definitive diagnosis of graft rejection however relies on a biopsy of the graft. In the acute rejection (AR) diagnostics there is a need for an alternative to biopsy that would be an easily repeatable and simple method for regular use. Frequent surveillance of acute or subclinical rejection (SCR) may improve long-term function. In this thesis, associations between cytokine and thrombosis associated candidate genes and the outcome of kidney Tx were studied. Cytotoxic and co-stimulatory T lymphocyte molecule gene expression biomarkers for the diagnosis of the AR and the SCR were also investigated. We found that polymorphisms in the cytokine genes tumor necrosis factor and interleukin 10 (IL10) of the recipients were associated with AR. In addition, certain IL10 gene polymorphisms of the donors were associated with the incidence of cytomegalovirus infection and occurrence of later infection in a subpopulation of recipients. Further, polymorphisms in genes related to the risk of thrombosis and those of certain cytokines were not associated with the occurrence of thrombosis, infarction, AR or graft survival. In the study of biomarkers for AR, whole blood samples were prospectively collected from adult kidney Tx patients. With real-time quantitative PCR (RT-QPCR) gene expression quantities of CD154 and ICOS differentiated the patients with AR from those without, but not from the patients with other causes of graft dysfunction. Biomarkers for SCR were studied in paediatric kidney Tx patients. We used RT-QPCR to quantify the gene expression of immunological candidate genes in a low-density array format. In addition, we used RT-QPCR to validate the results of the microarray analysis. No gene marker differentiated patients with SCR from those without SCR. This research demonstrates the lack of robust markers among polymorphisms or biomarkers in investigated genes that could be included in routine analysis in a clinical laboratory. In genetic studies, kidney Tx can be regarded as a complex trait, i.e. several environmental and genetic factors may determine its outcome. A number of currently unknown genetic factors probably influence the results of Tx.

Molecular Genetics of Age-related Macular Degeneration

Age-related macular degeneration (AMD; OMIM # 603075) is an eye disease of the elderly, signs of which appear after the age of 50. In the Western world it is a leading cause of permanent visual loss with a prevalence of 8.5% in persons under 54 years of age and of 37% in persons over 75 years of age. Early forms of AMD may be asymptomatic, but in the late forms usually a central scotoma in the visual field follows severely complicating daily tasks. Smoking, age, and genetic predisposition are known risk factors for AMD. Until recently no true susceptibility genes had been identified though the composition of drusen deposits, the hallmarks of AMD, has suggested that the complement system might play a role in the pathogenesis of AMD. When four groups reported in March 2005, that, on chromosome 1q32, a Y402H variant in the complement factor H (CFH) gene confers risk for AMD in independent Caucasian samples, a new period in the field of genetic research of AMD started. CFH is a key regulator of the complement system. Thus, it is logical to speculate, that it plays a role in the pathogenesis of AMD. We performed a case-control association study to analyse whether the CFH Y402H variant contain a risk for AMD in the Finnish population. Although the population of Finland represents a genetic isolate, the CFH Y402H polymorphism was associated with AMD also in our patient sample with similar risk allele frequencies as in the other Caucasian populations. We further evaluated the effects of this variant, but no association between lesion subtype (predominantly classic, minimally classic or occult lesion) or lesion size of neovascular AMD and the CFH Y402H variant was detected. Neither did the variant have an effect on the photodynamic therapy (PDT) outcome. The patients that respond to PDT carried the risk genotype as frequently as those who did not respond, and no difference was found in the number of PDT sessions needed in patients with or without the risk genotypes of CFH Y402H. Functional analyses, however, showed that the binding of C-reactive protein (CRP) to CFH was significantly reduced in patients with the risk genotype of Y402H. In the past two years, the LOC387715/ high-temperature requirement factor A1 (HTRA1) locus on 10q26 has also been repeatedly associated with AMD in several populations. The recent discovery of the LOC387715 protein on the mitochondrial outer membrane suggests that the LOC387715 gene, not HTRA1, is the true predisposing gene in this region, although its biological function is still unknown. In our Finnish patient material, patients with AMD carried the A69S risk genotype of LOC387715 more frequently than the controls. Also, for the first time, an interaction between the CFH Y402H and the LOC387715 A69S variants was found. The most recently detected susceptibilty gene of AMD, the complement component 3 (C3) gene, encodes the central component of the complement system, C3. In our Finnish sample, an additive gene effect for the C3 locus was detected, though weaker than the effects for the two main loci, CFH and LOC387715. Instead, the hemicentin-1 or the elongation of very long chain fatty acids-like 4 genes that have also been suggested as candidate genes for AMD did not carry a risk for AMD in the Finnish population. This was the first series of molecular genetic study of AMD in Finland. We showed that two common risk variants, CFH Y402H and LOC387715 A69S, represent a high risk of AMD also in the isolated Finnish population, and furthermore, that they had a statistical interaction. It was demonstrated that the CFH Y402H risk genotype affects the binding of CFH to CRP thus suggesting that complement indeed plays an important role in the pathogenesis of AMD.

Ion channel gene variants predisposing to severe human ventricular arrhythmias

Congenital long QT syndrome (LQTS) with an estimated prevalence of 1:2000-1:10 000 manifests with prolonged QT interval on electrocardiogram and risk for ventricular arrhythmias and sudden death. Several ion channel genes and hundreds of mutations in these genes have been identified to underlie the disorder. In Finland, four LQTS founder mutations of potassium channel genes account for up to 40-70% of genetic spectrum of LQTS. Acquired LQTS has similar clinical manifestations, but often arises from usage of QT-prolonging medication or electrolyte disturbances. A prolonged QT interval is associated with increased morbidity and mortality not only in clinical LQTS but also in patients with ischemic heart disease and in the general population. The principal aim of this study was to estimate the actual prevalence of LQTS founder mutations in Finland and to calculate their effect on QT interval in the Finnish background population. Using a large population-based sample of over 6000 Finnish individuals from the Health 2000 Survey, we identified LQTS founder mutations KCNQ1 G589D (n=8), KCNQ1 IVS7-2A>G (n=1), KCNH2 L552S (n=2), and KCNH2 R176W (n=16) in 27 study participants. This resulted in a weighted prevalence estimate of 0.4% for LQTS in Finland. Using a linear regression model, the founder mutations resulted in a 22- to 50-ms prolongation of the age-, sex-, and heart rate-adjusted QT interval. Collectively, these data suggest that one of 250 individuals in Finland may be genetically predisposed to ventricular arrhythmias arising from the four LQTS founder mutations. A KCNE1 D85N minor allele with a frequency of 1.4% was associated with a 10-ms prolongation in adjusted QT interval and could thus identify individuals at increased risk of ventricular arrhythmias at the population level. In addition, the previously reported associations of KCNH2 K897T, KCNH2 rs3807375, and NOS1AP rs2880058 with QT interval duration were confirmed in the present study. In a separate study, LQTS founder mutations were identified in a subgroup of acquired LQTS, providing further evidence that congenital LQTS gene mutations may underlie acquired LQTS. Catecholaminergic polymorphic ventricular tachycardia (CPVT) is characterized by exercise-induced ventricular arrhythmias in a structurally normal heart and results from defects in the cardiac Ca2+ signaling proteins, mainly ryanodine receptor type 2 (RyR2). In a patient population of typical CPVT, RyR2 mutations were identifiable in 25% (4/16) of patients, implying that noncoding variants or other genes are involved in CPVT pathogenesis. A 1.1 kb RyR2 exon 3 deletion was identified in two patients independently, suggesting that this region may provide a new target for RyR2-related molecular genetic studies. Two novel RyR2 mutations showing a gain-of-function defect in vitro were identified in three victims of sudden cardiac death. Extended pedigree analyses revealed some surviving mutation carriers with mild structural abnormalities of the heart and resting ventricular arrhythmias suggesting that not all RyR2 mutations lead to a typical CPVT phenotype, underscoring the relevance of tailored risk stratification of a RyR2 mutation carrier.

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

The many ways to delimit species: hairs, genes and surface chemistry

Species identification forms the basis for understanding the diversity of the living world, but it is also a prerequisite for understanding many evolutionary patterns and processes. The most promising approach for correctly delimiting and identifying species is to integrate many types of information in the same study. Our aim was to test how cuticular hydro- carbons, traditional morphometrics, genetic polymorphisms in nuclear markers (allozymes and DNA microsatellites) and DNA barcoding (partial mitochondrial COI gene) perform in delimiting species. As an example, we used two closely related Formica ants, F. fusca and F. lemani, sampled from a sympatric population in the northern part of their distribu- tion. Morphological characters vary and overlap in different parts of their distribution areas, but cuticular hydrocarbons include a strong taxonomic signal and our aim is to test the degree to which morphological and genetic data correspond to the chemical data. In the morphological analysis, species were best separated by the combined number of hairs on pro- notum and mesonotum, but individual workers overlapped in hair numbers, as previously noted by several authors. Nests of the two species were separated but not clustered according to species in a Principal Component Analysis made on nuclear genetic data. However, model-based Bayesian clustering resulted in perfect separation of the species and gave no indication of hybridization. Furthermore, F. lemani and F. fusca did not share any mitochondrial haplotypes, and the species were perfectly separated in a phylogenetic tree. We conclude that F. fusca and F. lemani are valid species that can be separated in our study area relatively well with all methods employed. However, the unusually small genetic differen- tiation in nuclear markers (FST = 0.12) shows that they are closely related, and occasional hybridization between F. fusca and F. lemani cannot be ruled out.