Gene Expression : From Microarrays to Functional Genomics

The time of the large sequencing projects has enabled unprecedented possibilities of investigating more complex aspects of living organisms. Among the high-throughput technologies based on the genomic sequences, the DNA microarrays are widely used for many purposes, including the measurement of the relative quantity of the messenger RNAs. However, the reliability of microarrays has been strongly doubted as robust analysis of the complex microarray output data has been developed only after the technology had already been spread in the community. An objective of this study consisted of increasing the performance of microarrays, and was measured by the successful validation of the results by independent techniques. To this end, emphasis has been given to the possibility of selecting candidate genes with remarkable biological significance within specific experimental design. Along with literature evidence, the re-annotation of the probes and model-based normalization algorithms were found to be beneficial when analyzing Affymetrix GeneChip data. Typically, the analysis of microarrays aims at selecting genes whose expression is significantly different in different conditions followed by grouping them in functional categories, enabling a biological interpretation of the results. Another approach investigates the global differences in the expression of functionally related groups of genes. Here, this technique has been effective in discovering patterns related to temporal changes during infection of human cells. Another aspect explored in this thesis is related to the possibility of combining independent gene expression data for creating a catalog of genes that are selectively expressed in healthy human tissues. Not all the genes present in human cells are active; some involved in basic activities (named housekeeping genes) are expressed ubiquitously. Other genes (named tissue-selective genes) provide more specific functions and they are expressed preferably in certain cell types or tissues. Defining the tissue-selective genes is also important as these genes can cause disease with phenotype in the tissues where they are expressed. The hypothesis that gene expression could be used as a measure of the relatedness of the tissues has been also proved. Microarray experiments provide long lists of candidate genes that are often difficult to interpret and prioritize. Extending the power of microarray results is possible by inferring the relationships of genes under certain conditions. Gene transcription is constantly regulated by the coordinated binding of proteins, named transcription factors, to specific portions of the its promoter sequence. In this study, the analysis of promoters from groups of candidate genes has been utilized for predicting gene networks and highlighting modules of transcription factors playing a central role in the regulation of their transcription. Specific modules have been found regulating the expression of genes selectively expressed in the hippocampus, an area of the brain having a central role in the Major Depression Disorder. Similarly, gene networks derived from microarray results have elucidated aspects of the development of the mesencephalon, another region of the brain involved in Parkinson Disease.

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 mechanisms of cancer predisposition in HNPCC/Lynch syndrome

Hereditary non-polyposis colorectal carcinoma (HNPCC; Lynch syndrome) is among the most common hereditary cancers in man and a model of cancers arising through deficient DNA mismatch repair (MMR). It is inherited in a dominant manner with predisposing germline mutations in the MMR genes, mainly MLH1, MSH2, MSH6 and PMS2. Both copies of the MMR gene need to be inactivated for cancer development. Since Lynch syndrome family members are born with one defective copy of one of the MMR genes in their germline, they only need to acquire a so called second hit to inactivate the MMR gene. Hence, they usually develop cancer at an early age. MMR gene inactivation leads to accumulation of mutations particularly in short repeat tracts, known as microsatellites, causing microsatellite instability (MSI). MSI is the hallmark of Lynch syndrome tumors, but is present in approximately 15% of sporadic tumors as well. There are several possible mechanisms of somatic inactivation (i.e. the second hit ) of MMR genes, for instance deletion of the wild-type copy, leading to loss of heterozygosity (LOH), methylation of promoter regions necessary for gene transcription, or mitotic recombination or gene conversion. In the Lynch syndrome tumors carrying germline mutations in the MMR gene, LOH was found to be the most frequent mechanism of somatic inactivation in the present study. We also studied MLH1/MSH2 deletion carriers and found that somatic mutations identical to the ones in the germline occurred frequently in colorectal cancers and were also present in extracolonic Lynch syndrome-associated tumors. Chromosome-specific marker analysis implied that gene conversion, rather than mitotic recombination or deletion of the respective gene locus accounted for wild-type inactivation. Lynch syndrome patients are predisposed to certain types of cancers, the most common ones being colorectal, endometrial and gastric cancer. Gastric cancer and uroepithelial tumors of bladder and ureter were observed to be true Lynch syndrome tumors with MMR deficiency as the driving force of tumorigenesis. Brain tumors and kidney carcinoma, on the other hand, were mostly MSS, implying the possibility of alternative routes of tumor development. These results present possible implications in clinical cancer surveillance. In about one-third of families suspected of Lynch syndrome, mutations in MMR genes are not found, and we therefore looked for alternative mechanisms of predisposition. According to our results, large genomic deletions, mainly in MSH2, and germline epimutations in MLH1, together explain a significant fraction of point mutation-negative families suspected of Lynch syndrome and are associated with characteristic clinical and family features. Our findings have important implications in the diagnosis and management of Lynch syndrome families.

The Myopathic Protein Myotilin in Developing Mouse and in Muscle Function

Skeletal muscle cells are highly specialised in order to accomplish their function. During development, the fusion of hundreds of immature myoblasts creates large syncytial myofibres with a highly ordered cytoplasm filled with packed myofibrils. The assembly and organisation of contractile myofibrils must be tightly controlled. Indeed, the number of proteins involved in sarcomere building is impressive, and the role of many of them has only recently begun to be elucidated. Myotilin was originally identified as a high affinity a-actinin binding protein in yeast twohybrid screen. It was then found to interact also with filamin C, actin, ZASP and FATZ-1. Human myotilin is mainly expressed in striated muscle and induces efficient actin bundling in vitro and in cells. Moreover, mutations in myotilin cause different forms of muscle disease, now collectively known as myotilinopathies. In this thesis, consisting of three publications, the work on the mouse orthologue is presented. First, the cloning and molecular characterisation of the mouse myotilin gene showed that human and mouse myotilin share high sequence homology and a similar expression pattern and gene regulation. Functional analysis of the mouse promoter revealed the myogenic factor-binding elements that are required for myotilin gene transcription. Secondly, expression of myotilin was studied during mouse embryogenesis. Surprisingly, myotilin was expressed in a wide array of tissues at some stages of development; its expression pattern became more restricted at perinatal stages and in adult life. Immunostaining of human embryos confirmed broader myotilin expression compared to the sarcomeric marker titin. Finally, in the third article, targeted deletion of myotilin gene in mice revealed that it is not essential for muscle development and function. These data altogether indicate that the mouse can be used as a model for human myotilinopathy and that loss of myotilin does not alter significantly muscle structure and function. Therefore, disease-associated mutant myotilin may act as a dominant myopathic factor.

Molecular biology and functions of epilysin (MMP-28)

Epilysin (MMP-28) is the most recently identified member of the matrix metalloproteinase (MMP) family of extracellular proteases. Together these enzymes are capable of degrading almost all components of the extracellular matrix (ECM) and are thus involved in important biological processes such as development, wound healing and immune functions, but also in pathological processes such as tumor invasion, metastasis and arthritis. MMPs do not act solely by degrading the ECM. They also regulate cell behavior by releasing growth factors and biologically active peptides from the ECM, by modulating cell surface receptors and adhesion molecules and by regulating the activity of many important mediators in inflammatory pathways. The aim of this study was to define the unique role of epilysin within the MMP-family, to elucidate how and when it is expressed and how its catalytic activity is regulated. To gain information on its essential functions and substrates, the specific aim was to characterize how epilysin affects the phenotype of epithelial cells, where it is biologically expressed. During the course of the study we found that the epilysin promoter contains a well conserved GT-box that is essential for the basic expression of this gene. Transcription factors Sp1 and Sp3 bind this sequence and could hence regulate both the basic and cell type and differentiation stage specific expression of epilysin. We cloned mouse epilysin cDNA and found that epilysin is well conserved between human and mouse genomes and that epilysin is glycosylated and activated by furin. Similarly to in human tissues, epilysin is normally expressed in a number of mouse tissues. The expression pattern differs from most other MMPs, which are expressed only in response to injury or inflammation and in pathological processes like cancer. These findings implicate that epilysin could be involved in tissue homeostasis, perhaps fine-tuning the phenotype of epithelial cells according to signals from the ECM. In view of these results, it was unexpected to find that epilysin can induce a stable epithelial to mesenchymal transition (EMT) when overexpressed in epithelial lung carcinoma cells. Transforming growth factor b (TGF-b) was recognized as a crucial mediator of this process, which was characterized by the loss of E-cadherin mediated cell-cell adhesion, elevated expression of gelatinase B and MT1-MMP and increased cell migration and invasion into collagen I gels. We also observed that epilysin is bound to the surface of epithelial cells and that this interaction is lost upon cell transformation and is susceptible to degradation by membrane type-1-MMP (MT1-MMP). The wide expression of epilysin under physiological conditions implicates that its effects on epithelial cell phenotype in vivo are not as dramatic as seen in our in vitro cell system. Nevertheless, current results indicate a possible interaction between epilysin and TGF-b also under physiological circumstances, where epilysin activity may not induce EMT but, instead, trigger less permanent changes in TGF-b signaling and cell motility. Epilysin may thus play an important role in TGF-b regulated events such as wound healing and inflammation, processes where involvement of epilysin has been indicated.

Hypertension, kardiotoxicitet och njursvikt i samband med tyrosinkinashämmaren sunitinib

Nybildning av blodkärl från tidigare existerande kärl, angiogenes, är ett väsentligt skede vid tumörtillväxt. Denna process regleras av bland annat tillväxtfaktorer, var av den vaskulära endoteliala tillväxtfaktorn har en central roll. Hämning av angiogenes kan ske antingen extracellulärt med hjälp av humaniserade monoklonala antikroppar eller intracellulärt med hjälp av småmolekylära hämmaren. Sunitinib är en småmolekylär multikinashämmare och inhiberar flera tyrosinkinasreceptorer som påverkar tumörtillväxten och metastasutvecklingen vid cancer. Sunitinibs främsta indikationer är gastrointestinala stromacellstumörer, metastaserad njurcellscancer och neuroendokrina tumörer i bukspottskörteln. Behandling med tyrosinkinashämmare orsakar biverkningar som hypertension, kardiotoxicitet och njursvikt, vilka antas bero på de hämmande effekterna på mål som inte är väsentliga för anti-cancer-aktiviteten (”off-target” biverkningar). Bland annat AMP-aktiverat proteinkinas (AMPK), ett kinas som upprätthåller metabolisk homeostas i hjärtat, inhiberas av sunitinib och antas framkalla kardiovaskulära biverkningar. För att reducera ”off-target” biverkningar strävar man till att hitta alternativ som minskar de skadliga effekterna utan att den terapeutiska aktiviteten försvagas. Bland annat ett begränsat kaloriintag har uppvisat skyddande effekt på hjärtat via mekanismer sammankopplade till ökad resistens mot oxidativ stress, inflammation och mitokondriell dysfunktion, samt avtagande apoptos och autofagi. Detta sker delvis genom aktivering av enzymet Sirt1. Syftet med den här studien var att undersöka ifall kaloribegränsning skyddar mot kardiovaskulära och renala biverkningar inducerade av sunitinib hos råttor. Dessutom studerades vilka signalkedjor i cellen som medverkar. I studien användes 40 spontant hypertensiva råttor samt 10 normotensiva Wistar-Kyoto råttor. Försöksdjuren delades in i fem grupper beroende på behandling; I WKY kontroll, II SHR kontroll, III SHR + kaloribegränsning 70 %, IV SHR + sunitinib 3 mg/kg och V SHR + sunitinib 3 mg/kg + kaloribegränsning 70 %. Behandlingsperioden var åtta veckor. Blodtrycket mättes varje vecka med svansmanchett, urinutsöndringen undersöktes vecka 4 och vecka 8 med metabolismburar, ultraljudsundersökning av hjärtat utfördes sista veckan och blodkärlens respons till acetylkolin och natriumnitroprussid studerades i samband med avlivning. Proteinerna Sirt1 och AMPK analyserades i hjärtat med Western blotting samt förekomsten av makrofagmarkören ED1 i njurarna med immunhistokemi. Studien visade att sunitinibdosen 3 mg/kg är mycket väl tolererbar hos råttor eftersom sunitinib inte orsakade högre blodtryck, kraftigare hypertrofi eller mer omfattande njurskada jämfört med obehandlade SHR- grupper. Utgående från resultaten kan man också konstatera att kaloribegränsningen har positiva kardiovaskulära effekter.

Mammalian Mat1 subunit of Transcription Factor IIH in gene-specific and general transcription

All protein-encoding genes in eukaryotes are transcribed into messenger RNA (mRNA) by RNA Polymerase II (RNAP II), whose activity therefore needs to be tightly controlled. An important and only partially understood level of regulation is the multiple phosphorylations of RNAP II large subunit C-terminal domain (CTD). Sequential phosphorylations regulate transcription initiation and elongation, and recruit factors involved in co-transcriptional processing of mRNA. Based largely on studies in yeast models and in vitro, the kinase activity responsible for the phosphorylation of the serine-5 (Ser5) residues of RNAP II CTD has been attributed to the Mat1/Cdk7/CycH trimer as part of Transcription Factor IIH. However, due to the lack of good mammalian genetic models, the roles of both RNAP II Ser5 phosphorylation as well as TFIIH kinase in transcription have provided ambiguous results and the in vivo kinase of Ser5 has remained elusive. The primary objective of this study was to elucidate the role of mammalian TFIIH, and specifically the Mat1 subunit in CTD phosphorylation and general RNAP II-mediated transcription. The approach utilized the Cre-LoxP system to conditionally delete murine Mat1 in cardiomyocytes and hepatocytes in vivo and and in cell culture models. The results identify the TFIIH kinase as the major mammalian Ser5 kinase and demonstrate its requirement for general transcription, noted by the use of nascent mRNA labeling. Also a role for Mat1 in regulating general mRNA turnover was identified, providing a possible rationale for earlier negative findings. A secondary objective was to identify potential gene- and tissue-specific roles of Mat1 and the TFIIH kinase through the use of tissue-specific Mat1 deletion. Mat1 was found to be required for the transcriptional function of PGC-1 in cardiomyocytes. Transriptional activation of lipogenic SREBP1 target genes following Mat1 deletion in hepatocytes revealed a repressive role for Mat1apparently mediated via co-repressor DMAP1 and the DNA methyltransferase Dnmt1. Finally, Mat1 and Cdk7 were also identified as a negative regulators of adipocyte differentiation through the inhibitory phosphorylation of Peroxisome proliferator-activated receptor (PPAR) γ. Together, these results demonstrate gene- and tissue-specific roles for the Mat1 subunit of TFIIH and open up new therapeutic possibilities in the treatment of diseases such as type II diabetes, hepatosteatosis and obesity.

Candidate gene studies on cardiovascular traits

Cardiovascular diseases (CVD) are a major cause of death and disability in Western countries and a growing health problem in the developing world. The genetic component of both coronary heart disease (CHD) and ischemic stroke events has been established in twin studies, and the traits predisposing to CVD, such as hypertension, dyslipidemias, obesity, diabetes, and smoking behavior, are all partly hereditary. Better understanding of the pathophysiology of CVD-related traits could help to target disease prevention and clinical treatment to individuals at an especially high disease risk and provide novel pharmaceutical interventions. This thesis aimed to clarify the genetic background of CVD at a population level using large Nordic population cohorts and a candidate gene approach. The first study concentrated on the allelic diversity of the thrombomodulin (THBD) gene in two Finnish cohorts, FINRISK-92 and FINRISK-97. The results from this study implied that THBD variants do not substantially contribute to CVD risk. In the second study, three other candidate genes were added to the analyses. The study investigated the epistatic effects of coagulation factor V (F5), intercellular adhesion molecule -1 (ICAM1), protein C (PROC), and THBD in the same FINRISK cohorts. The results were encouraging; we were able to identify several single SNPs and SNP combinations associating with CVD and mortality. Interestingly, THBD variants appeared in the associating SNP combinations despite the negative results from Study I, suggesting that THBD contributes to CVD through gene-gene interactions. In the third study, upstream transcription factor -1 (USF1) was analyzed in a cohort of Swedish men. USF1 was associated with metabolic syndrome, characterized by accumulation of different CVD risk factors. A putative protective and a putative risk variant were identified. A direct association with CVD was not observed. The longitudinal nature of the study also clarified the effect of USF1 variants on CVD risk factors followed in four examinations throughout adulthood. The three studies provided valuable information on the study of complex traits, highlighting the use of large study samples, the importance of replication, and the full coverage of the major allelic variants of the target genes to assure reliable findings. Although the genetic basis of coronary heart disease and ischemic stroke remains unknown, single genetic findings may facilitate the recognition of high-risk subgroups.

Molecular Pathways of Disturbed Sleep and Depression: Studies on Adenosine and Gene Expression Patterns

Background: Adenosine is a potent sleep-promoting substance, and one of its targets is the basal forebrain. Fairly little is known about its mechanism of action in the basal forebrain and about the receptor subtype mediating its regulating effects on sleep homeostasis. Homeostatic deficiency might be one of the causes of the profoundly disturbed sleep pattern in major depressive disorder, which could explain the reduced amounts of delta-activity-rich stages 3 and 4. Since major depression has a relatively high heritability, and on the other hand adenosine regulates sleep homeostasis and might also be involved in mood modulation, adenosine-related genes should be considered for their possible contribution to a predisposition for depression and disturbed sleep in humans. Depression is a complex disorder likely involving the abnormal functioning of several genes. Novel target genes which could serve as the possible common substrates for depression and comorbid disturbed sleep should be identified. In this way specific brain areas related to sleep regulation should be studied by using animal model of depression which represents more homogenous phenotype as compared to humans. It is also important to study these brain areas during the development of depressive-like features to understand how early changes could facilitate pathophysiological changes in depression. Aims and methods: We aimed to find out whether, in the basal forebrain, adenosine induces recovery non-rapid eye movement (NREM) sleep after prolonged waking through the A1 or/and A2A receptor subtype. A1 and A2A receptor antagonists were perfused into the rat basal forebrain during 3 h of sleep deprivation, and the amount of NREM sleep and delta power during recovery NREM sleep were analyzed. We then explored whether polymorphisms in genes related to the metabolism, transport and signaling of adenosine could predispose to depression accompanied by signs of disturbed sleep. DNA from 1423 individuals representative of the Finnish population and including controls and cases with depression, depression accompanied by early morning awakenings and depression accompanied by fatigue, was used in the study to investigate the possible association between polymorphisms from adenosine-related genes and cases. Finally to find common molecular substrates of depression and disturbed sleep, gene expression changes were investigated in specific brain areas in the rat clomipramine model of depression. We focused on the basal forebrain of 3-week old clomipramine-treated rats which develop depressive-like symptoms later in adulthood and on the hypothalamus of adult female clomipramine-treated rats. Results: Blocking of the A1 receptor during sleep deprivation resulted in a reduction of the recovery NREM sleep amount and delta power, whereas A2A receptor antagonism had no effect. Polymorphisms in adenosine-related genes SLC29A3 (equilibrative nucleoside transporter type 3) in women and SLC28A1 (concentrative nucleoside transporter type 1) in men associated with depression alone as well as when accompanied by early morning awakenings and fatigue. In Study III the basal forebrain of postnatal rats treated with clomipramine displayed disturbances in gamma-aminobutyric acid (GABA) receptor type A signaling, in synaptic transmission and possible epigenetic changes. CREB1 was identified as a common transcription denominator which also mediates epigenetic regulation. In the hypothalamus the major changes included the expression of genes in GABA-A receptor pathway, K+ channel-related, glutamatergic and mitochondrial genes, as well as an overexpression of genes related to RNA and mRNA processing. Conclusions: Adenosine plays an important role in sleep homeostasis by promoting recovery NREM sleep via the A1 receptor subtype in the basal forebrain. Also adenosine levels might contribute to the risk of depression with disturbed sleep, since the genes encoding nucleoside transporters showed the strongest associations with depression alone and when accompanied by signs of disturbed sleep in both women and men. Sleep and mood abnormalities in major depressive disorder could be a consequence of multiple changes at the transcriptional level, GABA-A receptor signaling and synaptic transmission in sleep-related basal forebrain and the hypothalamus.

Genetics of Cardiovascular Disease: A Candidate Gene Study of USF1

Cardiovascular diseases (CVD) are major contributors to morbidity and mortality worldwide. Several interacting environmental, biochemical, and genetic risk factors can increase disease susceptibility. While some of the genes involved in the etiology of CVD are known, many are yet to be discovered. During the last few decades, scientists have searched for these genes with genome-wide linkage and association methods, and with more targeted candidate gene studies. This thesis investigates variation within the upstream transcription factor 1 (USF1) gene locus in relation to CVD risk factors, atherosclerosis, and incidence and prevalence of CVD. This candidate gene was first identified in Finnish families ascertained for familial combined hyperlipidemia, a common dyslipidemia predisposing to coronary heart disease. The gene is a ubiquitously expressed transcription factor regulating expression of several genes from lipid and glucose metabolism, inflammation, and endothelial function. First, we examined association between USF1 variants and several CVD risk factors, such as lipid phenotypes, body composition measures, and metabolic syndrome, in two prospective population cohorts. Our data suggested that USF1 contributes to these CVD risk factors at the population level. Notably, the associations with quantitative measurements were mostly detected among study subjects with CVD or metabolic syndrome, suggesting complex interactions between USF1 effects and the pathophysiological state of an individual. Second, we investigated how variation at the USF1 locus contributes to atherosclerotic lesions of the coronary arteries and abdominal aorta. For this, we used two study samples of middle-aged men with detailed measurements of atherosclerosis obtained in autopsy. USF1 variation significantly associated with areas of several types of lesions, especially with calcification of the arteries. Next, we tested what effect the USF1 risk variants have on sudden cardiac death and incidence of CVD. The atherosclerosis-associated risk variant increased the risk of sudden cardiac death of the same study subjects. Furthermore, USF1 alleles associated with incidence of CVD in the Finnish population follow-up cohorts. These associations were especially prominent among women, suggesting a sex specific effect, which has also been detected in subsequent studies. Finally, as some of the low-yield DNA samples of the Finnish follow-up study cohort needed to be whole-genome amplified (WGA) prior to genotyping, we evaluated whether the produced WGA genotypes were of good quality. Although the samples giving genotype discrepancies could not be detected before genotyping with standard laboratory quality control methods, our results suggested that enhanced quality control at the time of the genotyping could identify such samples. In addition, combining two WGA reactions into one pooled DNA sample for genotyping markedly reduced the number of discrepancies and samples showing them. In conclusion, USF1 seems to have a role in the etiology of CVD. Additional studies are warranted to identify functional variants and to study interactions between USF1 and other genetic or environmental factors. This USF1 study, and other studies with low DNA yield of some samples, can benefit from whole genome amplification of the low-yield samples prior to genotyping. Careful quality control procedures are, however, needed in WGA genotyping.

Characterization of the TRIM37 gene and mutations underlying mulibrey nanism

Mulibrey nanism is a hereditary developmental disorder, characterized by prenatal onset growth failure without postnatal catch-up growth, distinctive craniofacial features, progressive cardiopathy and failure of sexual maturation. In addition, the patients develop insulin resistance syndrome and type 2 diabetes and they have an increased risk of developing tumors. The TRIM37 gene that underlies mulibrey nanism encodes for a member of the tripartite motif (TRIM) protein family. The physiological function of TRIM37 and the pathogenetic mechanisms leading from TRIM37 dysfunction to the mulibrey nanism phenotype are unknown. However, TRIM37 localizes at least partially to peroxisomes, and possesses ubiquitin E3-ligase activity. Thus, it may mediate ubiquitin dependent protein degradation, suggesting that accumulation of yet unknown substrate proteins may underlie the disease pathogenesis. In this study, the TRIM37 gene was characterized in detail. A transcription initiation window, with several separate transcription start sites, was identified and the putative promoter region immediately upstream from the transcription initiation window was shown to possess basal promoter activity. Further, several alternative splice variants of the gene were identified, including a highly expressed testis specific variant, encoding for an identical protein product with the main transcript. Expression of TRIM37 mRNA was detected in several different tissues, with highest expression seen in testis and in brain, when the expression patterns of the two major transcripts in different human tissues were studied by quantitative real-time PCR. Several mulibrey nanism patients were studied and thirteen novel mutations in TRIM37 were found, including three mutations (p.Gly322Val, p.Cys109Ser, p.Glu271_Ser287), that are likely to express mutant TRIM37 proteins. These mutations were further shown to alter the subcellular localization of the mutant proteins. Most of the mulibrey nanism associated mutations however, lead to premature termination codons and degradation of mRNA. All the TRIM37 mutations identified to date predict loss-of-function alleles, and thus no phenotype-genotype correlation is seen among the patients. In order to understand the pathogenetic mechanisms underlying mulibrey nanism, an animal model for the disorder is needed. For the development of a Trim37 knock-out mouse, the mouse Trim37 gene was characterized. Alternative splice variants, were identified, including a testis specific variant predicting a longer protein product. Further, a strictly tissue and cell-specific pattern of Trim37 expression was observed in developing and adult mouse tissues, when studied by immunohistochemical methods. This distribution of Trim37 expression in mouse tissues is in agreement with the clinical findings in human mulibrey nanism patients. This thesis work gives new tools for the diagnostics of mulibrey nanism as well as for studying the molecular pathogenesis behind this interesting disorder.

L-rhamnose-1-dehydrogenase gene and L-rhamnose catabolism in the yeast Pichia stipitis

The purpose of this work was to identify some of the genes of the catabolic route of L-rhamnose in the yeast Pichia stipitis. There are at least two distinctly different pathways for L-rhamnose catabolism. The one described in bacteria has phosphorylated intermediates and the enzymes and the genes of this route have been described. The pathway described in yeast does not have phosphorylated intermediates. The intermediates and the enzymes of this pathway are known but none of the genes have been identified. The work was started by purifying the L-rhamnose dehydrogenase, which oxidates L-rhamnose to rhamnonic acid-gamma-lactone. NAD is used as a cofactor in this reaction. A DEAE ion exchange column was used for purification. The active fraction was further purified using a non-denaturing PAGE and the active protein identified by zymogram staining. In the last step the protein was separated in a SDS-PAGE, the protein band trypsinated and analysed by MALDI-TOF MS. This resulted in the identification of the corresponding gene, RHA1, which was then, after a codon change, expressed in Saccharomyces cerevisiae. Also C- or N-terminal histidine tags were added but as the activity of the enzyme was lost or strongly reduced these were not used. The kinetic properties of the protein were analysed in the cell extract. Substrate specifity was tested with different sugars; L-rhamnose, L-lyxose and L-mannose were oxidated by the enzyme. Vmax values were 180 nkat/mg, 160 nkat/mg and 72 nkat/mg, respectively. The highest affinity was towards L-rhamnose, the Km value being 0.9 mM. Lower affinities were obtained with L-lyxose, Km 4.3 mM, and L-mannose Km 25 mM. Northern analysis was done to study the transcription of RHA1 with different carbon sources. Transcription was observed only on L-rhamnose suggesting that RHA1 expression is L-rhamnose induced. A RHA1 deletion cassette for P. stipitis was constructed but the cassette had integrated randomly and not targeted to delete the RHA1 gene. Enzyme assays for L-lactaldehyde dehydrogenase were done similarly to L-rhamnose dehydrogenase assays. NAD is used as a cofactor also in this reaction where L-lactaldehyde is oxidised to L-lactate. The observed enzyme activities were very low and the activity was lost during the purification procedures.

Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

Regulation of Growth by Drosophila FoxO Transcription Factor

Forkhead box class O (FoxO) transcription factors are members of the forkhead box transcription factor superfamily, with orthologues in various species such as human, worm and fly. FoxO proteins are key regulators of growth, metabolism, stress resistance and, consequently, life span. FoxOs integrate signals from different pathways, e.g. the growth controlling Insulin-TOR signaling pathway and the stress induced JNK and Hippo signaling pathways. FoxO proteins have evolved to guide the cellular response to varying energy and stress conditions by inducing the expression of genes involved in the regulation of growth and metabolism. This work has aimed to deepen the understanding of how FoxO executes its biological functions. A particular emphasis has been laid to its role in growth control. Specifically, evidence is presented indicating that FoxO restricts tissue growth in a situation when TOR signaling is high. This finding can have implications in a human condition called Tuberous sclerosis, manifested by multiple benign tumors. Further, it is shown that FoxO directly binds to the promoter and regulates the expression of a Drosophila Adenylate cyclase gene, ac76e, which in turn modulates the fly s development and growth systemically. These results strengthen FoxOs position among central size regulators as it is able to operate at the level of individual cells as well as in the whole organism. Finally, an attempt to reveal the regulatory network upstream of FoxO has been carried out. Several putative FoxO activity regulators were identified in an RNAi screen of Drosophila kinases and phosphatases. The results underscore that FoxO is regulated through an elaborate network, ensuring the correct execution of key cellular processes in metabolism and response to stress. Overall, the evidence provided in this study strengthens our view of FoxO as a key integrator of growth and stress signals.

Functional Analysis of MrpL55, Vig and Mat1 Acting at the Crossroad of Cell Cycle, Transcription and Metabolism Regulation

Cell proliferation, transcription and metabolism are regulated by complex partly overlapping signaling networks involving proteins in various subcellular compartments. The objective of this study was to increase our knowledge on such regulatory networks and their interrelationships through analysis of MrpL55, Vig, and Mat1 representing three gene products implicated in regulation of cell cycle, transcription, and metabolism. Genome-wide and biochemical in vitro studies have previously revealed MrpL55 as a component of the large subunit of the mitochondrial ribosome and demonstrated a possible role for the protein in cell cycle regulation. Vig has been implicated in heterochromatin formation and identified as a constituent of the RNAi-induced silencing complex (RISC) involved in cell cycle regulation and RNAi-directed transcriptional gene silencing (TGS) coupled to RNA polymerase II (RNAPII) transcription. Mat1 has been characterized as a regulatory subunit of cyclin-dependent kinase 7 (Cdk7) complex phosphorylating and regulating critical targets involved in cell cycle progression, energy metabolism and transcription by RNAPII. The first part of the study explored whether mRpL55 is required for cell viability or involved in a regulation of energy metabolism and cell proliferation. The results revealed a dynamic requirement of the essential Drosophila mRpL55 gene during development and suggested a function of MrpL55 in cell cycle control either at the G1/S or G2/M transition prior to cell differentiation. This first in vivo characterization of a metazoan-specific constituent of the large subunit of mitochondrial ribosome also demonstrated forth compelling evidence of the interconnection of nuclear and mitochondrial genomes as well as complex functions of the evolutionarily young metazoan-specific mitochondrial ribosomal proteins. In studies on the Drosophila RISC complex regulation, it was noted that Vig, a protein involved in heterochromatin formation, unlike other analyzed RISC associated proteins Argonaute2 and R2D2, is dynamically phosphorylated in a dsRNA-independent manner. Vig displays similarity with a known in vivo substrate for protein kinase C (PKC), human chromatin remodeling factor Ki-1/57, and is efficiently phosphorylated by PKC on multiple sites in vitro. These results suggest that function of the RISC complex protein Vig in RNAi-directed TGS and chromatin modification may be regulated through dsRNA-independent phosphorylation by PKC. In the third part of this study the role of Mat1 in regulating RNAPII transcription was investigated using cultured murine immortal fibroblasts with a conditional allele of Mat1. The results demonstrated that phosphorylation of the carboxy-terminal domain (CTD) of the large subunit of RNAPII in the heptapeptide YSPTSPS repeat in Mat-/- cells was over 10-fold reduced on Serine-5 and subsequently on Serine-2. Occupancy of the hypophosphorylated RNAPII in gene bodies was detectably decreased, whereas capping, splicing, histone methylation and mRNA levels were generally not affected. However, a subset of transcripts in absence of Mat1 was repressed and associated with decreased occupancy of RNAPII at promoters as well as defective capping. The results identify the Cdk7-CycH-Mat1 kinase submodule of TFIIH as a stimulatory non-essential regulator of transcriptional elongation and a genespecific essential factor for stable binding of RNAPII at the promoter region and capping. The results of these studies suggest important roles for both MrpL55 and Mat1 in cell cycle progression and their possible interplay at the G2/M stage in undifferentiated cells. The identified function of Mat1 and of TFIIH kinase complex in gene-specific transcriptional repression is challenging for further studies in regard to a possible link to Vig and RISC-mediated transcriptional gene silencing.