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.

The Mechanisms of Hedgehog Signal Transduction

Studies in both vertebrates and invertebrates have identified proteins of the Hedgehog (Hh) family of secreted signaling molecules as key organizers of tissue patterning. Initially discovered in Drosophila in 1992, Hh family members have been discovered in animals with body plans as diverse as those of mammals, insects and echinoderms. In humans three related Hh genes have been identified: Sonic, Indian and Desert hedgehog (Shh, Ihh and Dhh). Transduction of the Hh signal to the cytoplasm utilizes an unusual mechanism involving consecutive repressive interactions between Hh and its receptor components, Patched (Ptc) and Smoothened (Smo). Several cytoplasmic proteins involved in Hh signal transduction are known in Drosophila, but mammalian homologs are known only for the Cubitus interruptus (Ci) transcription factor (GLI(1-3)) and for the Ci/GLI-associated protein, Suppressor of Fused (Su(fu)). In this study I analyzed the mechanisms of how the Hh receptor Ptc regulates the signal transducer Smo, and how Smo relays the Shh signal from the cell surface to the cytoplasm ultimately leading to the activation of GLI transcription factors. In Drosophila, the kinesin-like protein Costal2 (Cos2) is required for suppression of Hh target gene expression in the absence of ligand, and loss of Cos2 causes embryonic lethality. Cos2 acts by bridging Smo to the Ci. Another protein, Su(Fu) exerts a weak suppressive influence on Ci activity and loss of Su(Fu) causes subtle changes in Drosophila wing pattern. This study revealed that domains in Smo that are critical for Cos2 binding in Drosophila are dispensable for mammalian Smo function. Furthermore, by analyzing the function of Su(Fu) and the closest mouse homologs of Cos2 by protein overexpression and RNA interference I found that inhibition of the Hh response pathway in the absence of ligand does not require Cos2 activity, but instead critically depends on the activity of Su(Fu). These results indicate that a major change in the mechanism of action of a conserved signaling pathway occurred during evolution, probably through phenotypic drift made possible by the existence in some species of two parallel pathways acting between the Hh receptor and the Ci/GLI transcription factors. In a second approach to unravel Hh signaling we cloned > 90% of all human full-length protein kinase cDNAs and constructed the corresponding kinase-activity deficient mutants. Using this kinome resource as a screening tool, two kinases, MAP3K10 and DYRK2 were found to regulate Shh signaling. DYRK2 directly phosphorylated and induced the proteasome dependent degradation of the key Hh-pathway regulated transcription factor, GLI2. MAP3K10, in turn, affected GLI2 indirectly by modulating the activity of DYRK2.

Function and regulation of Cdk7

Understanding the process of cell division is crucial for modern cancer medicine due to the central role of uncontrolled cell division in this disease. Cancer involves unrestrained proliferation as a result of cells loosing normal control and being driven through the cell cycle, where they normally would be non-dividing or quiescent. Progression through the cell cycle is thought to be dependent on the sequential activation of cyclin-dependent kinases (Cdks). The full activation of Cdks requires the phosphorylation of a conserved residue (threonine-160 on human Cdk2) on the T-loop of the kinase domain. In metazoan species, a trimeric complex consisting of Cdk7, cyclin H and Mat1 has been suggested to be the T-loop kinase of several Cdks. In addition, Cdk7 have also been implicated in the regulation of transcription. Cdk7, cyclin H, and Mat1 can be found as subunits of general transcription factor TFIIH. Cdk7, in this context, phosphorylates the Carboxy-terminal domain (CTD) of the large subunit of RNA polymerase II (RNA pol II), specifically on serine-5 residues of the CTD repeat. The regulation of Cdk7 in these and other functions is not well known and the unambiguous characterization of the in vivo role of Cdk7 in both T-loop activation and CTD serine-5 phosphorylation has proved challenging. In this study, the fission yeast Cdk7-cyclin H homologous complex, Mcs6-Mcs2, is identified as the in vivo T-loop kinase of Cdk1(Cdc2). It also identifies multiple levels of regulation of Mcs6 kinase activity, i.e. association with Pmh1, a novel fission yeast protein that is the apparent homolog of metazoan Mat1, and T-loop phosphorylation of Mcs6, mediated by Csk1, a monomeric T-loop kinase with similarity to Cak1 of budding yeast. In addition, Skp1, a component of the SCF (Skp1-Cullin-F box protein) ubiquitin ligase is identified by its interactions with Mcs2 and Pmh1. The Skp1 association with Mcs2 and Pmh1 is however SCF independent and does not involve proteolytic degradation but may reflect a novel mechanism to modulate the activity or complex assembly of Mcs6. In addition to Cdk7, also Cdk8 has been shown to have CTD serine-5 kinase activity in vitro. Cdk8 is not essential in yeast but has been shown to function as a transcriptional regulator. The function of Cdk8 is unknown in flies and mammals. This prompted the investigation of murine Cdk8 and its potential role as a redundant CTD serine-5 kinase. We find that Cdk8 is required for development prior to implantation, at a time that is co-incident with a burst of Cdk8 expression during normal development. The results does not support a role of Cdk8 as a serine-5 CTD kinase in vivo but rather shows an unexpected requirement for Cdk8, early in mammalian development. The results presented in this thesis extends our current knowledge of the regulation of the cell cycle by characterizing the function of two distinct cell cycle regulating T-loop kinases, including the unambiguous identification of Mcs6, the fission yeast Cdk7 homolog, as the T-loop kinase of Cdk1. The results also indicate that the function of Mcs6 is conserved from fission yeast to human Cdk7 and suggests novel mechanisms by which the distinct functions of Cdk7 and Mcs6 could be regulated. These findings are important for our understanding of how progression of the cell cycle and proper transcription is controlled, during normal development and tissue homeostasis but also under condition where cells have escaped these control mechanisms e.g. cancer.

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.

The function of Bmps and Runx2 in normal tooth development and in the pathogenesis of cleidocranial dysplasia

The development of many embryonic organs is regulated by reciprocal and sequential epithelial-mesenchymal interactions. These interactions are mediated by conserved signaling pathways that are reiteratively used. Cleidocranial dysplasia (CCD) is a congenital syndrome where both bone and tooth development is affected. The syndrome is characterized by short stature, abnormal clavicles, general bone dysplasia, and supernumerary teeth. CCD is caused by mutations in RUNX2, a transcription factor that is a key regulator of osteoblast differentiation and bone formation. The first aim of this study was to analyse the expression of a family of key signal molecules, Bone morphogenetic protein (Bmp) at different stages of tooth development. Bmps have a variety of functions and they were originally discovered as signals inducing ectopic bone formation. We performed a comparative in situ hybridisation analysis of the mRNA expression of Bmp2-7 from initiation of tooth development to differentiation of dental hard tissues. The expression patterns indicated that the Bmps signal between the epithelial and mesenchymal tissues during initiation and morphogenesis of tooth development, as well as during the differentiation of odontoblasts and ameloblasts. Furthermore, they are also part of the signalling networks whereby the enamel knot regulates the patterning of tooth cusps. The second aim was to study the role of Runx2 during tooth development and thereby to gain better understanding of the pathogenesis of the tooth phenotype in CCD. We analysed the tooth phenotype of Runx2 knockout mice and examined the patterns and regulation of Runx2 gene expression.. The teeth of wild-type and Runx2 mutant mice were compared by several methods including in situ hybridisation, tissue culture, bead implantation experiments, and epithelial-mesenchymal recombination studies. Phenotypic analysis of Runx2 -/- mutant tooth development showed that teeth failed to advance beyond the bud stage. Runx2 expression was restricted to dental mesenchyme between the bud and early bell stages of tooth development and it was regulated by epithelial signals, in particular Fgfs. We searched for downstream targets of Runx2 by comparative in situ hybridisation analysis. The expression of Fgf3 was downregulated in the mesenchyme of Runx2 -/- teeth. Shh expression was absent from the enamel knot in the lower molars of Runx2 -/- and reduced in the upper molars. In conclusion, these studies showed that Runx2 regulates key epithelial-mesenchymal interactions that control advancing tooth morphogenesis and histodifferentiation of the epithelial enamel organ. In addition, in the upper molars of Runx2 mutants extra buddings occured at the palatal side of the tooth bud. We suggest that Runx2 acts as an inhibitor of successional tooth formation by preventing advancing development of the buds. Accordingly, we propose that RUNX2 haploinsuffiency in humans causes incomplete inhibition of successional tooth formation and as a result supernumerary teeth.

Identification of the Meckel syndrome gene (MKS1) exposes a novel ciliopathy

Meckel syndrome (MKS, MIM 249000) is an autosomal recessive developmental disorder causing death in utero or shortly after birth. The hallmarks of the disease are cystic kidney dysplasia and fibrotic changes of the liver, occipital encephalocele with or without hydrocephalus and polydactyly. Other anomalies frequently seen in the patients are incomplete development of the male genitalia, club feet and cleft lip or palate. The clinical picture has been well characterized in the literature while the molecular pathology underlying the disease has remained unclear until now. In this study we identified the first MKS gene by utilizing the disease haplotypes in Finnish MKS families linked to the MKS1 locus on chromosome 17q23 (MKS1) locus. Subsequently, the genetic heterogeneity of MKS was established in the Finnish families. Mutations in at least four different genes can cause MKS. These genes have been mapped to the chromosomes 17q23 (MKS1), 11q13 (MKS2), 8q22 (MKS3) and 9q33 (MKS4). Two of these genes have been identified so far: The MKS1 gene (this work) and the MKS3 gene. The identified MKS1 gene was initially a novel human gene which is conserved among species. We found three different MKS mutations, one of them being the Finnish founder mutation. The information available from MKS1 orthologs in other species convinced us that the MKS1 gene is required for normal ciliogenesis. Defects of the cilial system in other human diseases and model organisms actually cause phenotypic features similar to those seen in MKS patients. The MKS3 (TMEM67) gene encodes a transmembrane protein and the gene maps to the syntenic Wpk locus in the rat, which is a model with polycystic kidney disease, agenesis of the corpus callosum and hydrocephalus. The available information from these two genes suggest that MKS1 would encode a structural component of the centriole required for normal ciliary functions, and MKS3 would be a transmembrane component most likely required for normal ciliary sensory signaling. The MKS4 locus was localized to chromosme 9q32-33 in this study by using an inbred Finnish family with two affected and two healthy children. This fourth locus contains TRIM32 gene, which is associated to another well characterized human ciliopathy, Bardet Biedl syndrome (BBS). Future studies should identify the MKS4 gene on chromosome 9q and confirm if there are more than two genes causing MKS Finnish families. The research on critical signaling pathways in organogenesis have shown that both Wnt and Hedgehog pathways are dependent on functional cilia. The MKS gene products will serve as excellent model molecules for more detailed studies of the functional role of cilia in organogenesis in more detail.

The Role Of Fumarase (FH) In Tumorigenesis

In this study, a predisposing gene for a recently characterized cancer syndrome, hereditary leiomyomatosis and renal cell cancer (HLRCC), was identified and the role of the gene was investigated in other familial cancers and in nonsyndromic tumorigenesis. HLRCC is a dominantly inherited disorder predisposing predominantly to uterine and skin leiomyomas, and also to renal cell cancer and uterine leiomyosarcoma. The disease gene was recently localized in Finnish families to 1q42-q43 by a genome-wide linkage search. Independently in the UK, a clinically similar condition, multiple cutaneous and uterine leiomyomata (MCUL), was linked to the same chromosomal region, strongly suggesting that HLRCC and MCUL are actually a single syndrome. Linkage results were confirmed by detecting loss of heterozygosity (LOH) at the disease locus in most of the patients' tumors, suggesting that this predisposing gene acts as a tumor suppressor. Through detailed investigation by genotyping of microsatellite markers and haplotype construction in Finnish and UK HLRCC/MCUL families we were able to narrow the disease locus down to 1.6 Mb. Extensive mutation screening of known and predicted transcripts in the target region resulted in identification of the HLRCC predisposing gene, fumarase (fumarate hydratase, FH). FH is a key enzyme in energy metabolism, catalyzing fumarate to malate in the tricarboxylic acid cycle (TCAC) in mitochondria. Germline alterations in FH segregating with the disease were detected in 25 of 42 HLRCC/MCUL families including whole-gene deletions, truncating small deletions/insertions and nonsense mutations, as well as substitutions or deletions of highly conserved amino acids. Biallelic inactivation was detected in almost all studied tumors of HLRCC patients. Furthermore, FH enzyme activity was reduced in the patients' normal tissues and was completely or virtually absent from tumors. Based on these findings, we extensively demonstrated that mutations in FH underlie the HLRCC/MCUL syndrome. In our studies of other familial cancers, evidence for involvement of FH defects was not found in familial prostate and breast cancers. To investigate the role of FH in sporadic tumorigenesis, we analyzed 652 lesions, including a series of 353 nonsyndromic counterparts of tumor types associated with HLRCC. Mutations in nonsyndromic tumors were rare and appeared to be limited to tumor types observed in the hereditary form of the disease. Biallelic inactivation of FH was detected in a uterine leiomyosarcoma, a cutaneous leiomyoma, a soft-tissue sarcoma, and in two uterine leiomyomas. In the uterine leiomyosarcoma and the cutaneous lesion FH mutations originated from the germline whereas the soft-tissue sarcoma harbored purely somatic changes. In uterine leiomyomas somatic mutations were detected in the two out of five tumors with LOH at the FH locus. Our findings demonstrate that FH inactivation is also involved in nonhereditary tumor development, and further support the hypothesis that FH acts as a tumor suppressor. The role of FH in predisposition to malignancies, renal cell carcinoma and leiomyosarcoma is important in the diagnosis and prevention of cancer among HLRCC patients. This study is of general clinical interest, because prior to our findings, little was known about the molecular genetics of uterine leiomyomas, the most common tumors of women.

Molecular basis of colorectal cancer predisposition

Colorectal cancer (CRC) is one of the most frequent malignancies in Western countries. Inherited factors have been suggested to be involved in 35% of CRCs. The hereditary CRC syndromes explain only ~6% of all CRCs, indicating that a large proportion of the inherited susceptibility is still unexplained. Much of the remaining genetic predisposition for CRC is probably due to undiscovered low-penetrance variations. This study was conducted to identify germline and somatic changes that contribute to CRC predisposition and tumorigenesis. MLH1 and MSH2, that underlie Hereditary non-polyposis colorectal cancer (HNPCC) are considered to be tumor suppressor genes; the first hit is inherited in the germline and somatic inactivation of the wild type allele is required for tumor initiation. In a recent study, frequent loss of the mutant allele in HNPCC tumors was detected and a new model, arguing against the two-hit hypothesis, was proposed for somatic HNPCC tumorigenesis. We tested this hypothesis by conducting LOH analysis on 25 colorectal HNPCC tumors with a known germline mutation in the MLH1 or MSH2 genes. LOH was detected in 56% of the tumors. All the losses targeted the wild type allele supporting the classical two-hit model for HNPCC tumorigenesis. The variants 3020insC, R702W and G908R in NOD2 predispose to Crohn s disease. Contribution of NOD2 to CRC predisposition has been examined in several case-control series, with conflicting results. We have previously shown that 3020insC does not predispose to CRC in Finnish CRC patients. To expand our previous study the variants R702W and G908R were genotyped in a population-based series of 1042 Finnish CRC patients and 508 healthy controls. Association analyses did not show significant evidence for association of the variants with CRC. Single nucleotide polymorphism (SNP) rs6983267 at chromosome 8q24 was the first CRC susceptibility variant identified through genome-wide association studies. To characterize the role of rs6983267 in CRC predisposition in the Finnish population, we genotyped the SNP in the case-control material of 1042 cases and 1012 controls and showed that G allele of rs6983267 is associated with the increased risk of CRC (OR 1.22; P=0.0018). Examination of allelic imbalance in the tumors heterozygous for rs6983267 revealed that copy number increase affected 22% of the tumors and interestingly, it favored the G allele. By utilizing a computer algorithm, Enhancer Element Locator (EEL), an evolutionary conserved regulatory motif containing rs6983267 was identified. The SNP affected the binding site of TCF4, a transcription factor that mediates Wnt signaling in cells, and has proven to be crucial in colorectal neoplasia. The preferential binding of TCF4 to the risk allele G was showed in vitro and in vivo. The element drove lacZ marker gene expression in mouse embryos in a pattern that is consistent with genes regulated by the Wnt signaling pathway. These results suggest that rs6983267 at 8q24 exerts its effect in CRC predisposition by regulating gene expression. The most obvious target gene for the enhancer element is MYC, residing ~335 kb downstream, however further studies are required to establish the transcriptional target(s) of the predicted enhancer element.

Latent TGF-β binding proteins -3 and -4 : transcriptional control and extracellular matrix targeting

Extracellular matrix (ECM) is a complex network of various proteins and proteoglycans which provides tissues with structural strength and resilience. By harvesting signaling molecules like growth factors ECM has the capacity to control cellular functions including proliferation, differentiation and cell survival. Latent transforming growth factor β (TGF-β) binding proteins (LTBPs) associate fibrillar structures of the ECM and mediate the efficient secretion and ECM deposition of latent TGF-β. The current work was conducted to determine the regulatory regions of LTBP-3 and -4 genes to gain insight into their tissue-specific expression which also has impact on TGF-β biology. Furthermore, the current research aimed at defining the ECM targeting of the N-terminal variants of LTBP-4 (LTBP-4S and -4L), which is required to understand their functions in tissues and to gain insight into conditions in which TGF-β is activated. To characterize the regulatory regions of LTBP-3 and -4 genes in silico and functional promoter analysis techniques were employed. It was found that the expression of LTBP-4S and -4L are under control of two independent promoters. This finding was in accordance with the observed expression patterns of LTBP-4S and -4L in human tissues. All promoter regions characterized in this study were TATAless, GC-rich and highly conserved between human and mouse species. Putative binding sites for Sp1 and GATA family of transcription factors were recognized in all of these regulatory regions. It is possible that these transcription factors control the basal expression of LTBP-3 and -4 genes. Smad binding element was found within the LTBP-3 and -4S promoter regions, but it was not present in LTBP-4L promoter. Although this element important for TGF-β signaling was present in LTBP-4S promoter, TGF-β did not induce its transcriptional activity. LTBP-3 promoter activity and mRNA expression instead were stimulated by TGF-β1 in osteosarcoma cells. It was found that the stimulatory effect of TGF-β was mediated by Smad and Erk MAPK signaling pathways. The current work explored the ECM targeting of LTBP-4S and identified binding partners of this protein. It was found that the N-terminal end of LTBP-4S possesses fibronectin (FN) binding sites which are critical for its ECM targeting. FN deficient fibroblasts incorporated LTBP-4S into their ECM only after addition of exogenous FN. Furthermore, LTBP-4S was found to have heparin binding regions, of which the C-terminal binding site mediated fibroblast adhesion. Soluble heparin prevented the ECM association of LTBP-4S in fibroblast cultures. In the current work it was observed that there are significant differences in the secretion, processing and ECM targeting of LTBP-4S and -4L. Interestingly, it was observed that most of the secreted LTBP-4L was associated with latent TGF-β1, whereas LTBP-4S was mainly secreted as a free form from CHO cells. This thesis provides information on transcriptional regulation of LTBP-3 and -4 genes, which is required for the deeper understanding of their tissue-specific functions. Further, the current work elucidates the structural variability of LTBPs, which appears to have impact on secretion and ECM targeting of TGF-β. These findings may advance understanding the abnormal activation of TGF-β which is associated with connective tissue disorders and cancer.

Genomic microarrays in chromosomal analysis of leukemia

Chromosomal alterations in leukemia have been shown to have prognostic and predictive significance and are also important minimal residual disease (MRD) markers in the follow-up of leukemia patients. Although specific oncogenes and tumor suppressors have been discovered in some of the chromosomal alterations, the role and target genes of many alterations in leukemia remain unknown. In addition, a number of leukemia patients have a normal karyotype by standard cytogenetics, but have variability in clinical course and are often molecularly heterogeneous. Cytogenetic methods traditionally used in leukemia analysis and diagnostics; G-banding, various fluorescence in situ hybridization (FISH) techniques, and chromosomal comparative genomic hybridization (cCGH), have enormously increased knowledge about the leukemia genome, but have limitations in resolution or in genomic coverage. In the last decade, the development of microarray comparative genomic hybridization (array-CGH, aCGH) for DNA copy number analysis and the SNP microarray (SNP-array) method for simultaneous copy number and loss of heterozygosity (LOH) analysis has enabled investigation of chromosomal and gene alterations genome-wide with high resolution and high throughput. In these studies, genetic alterations were analyzed in acute myeloid leukemia (AML) and chronic lymphocytic leukemia (CLL). The aim was to screen and characterize genomic alterations that could play role in leukemia pathogenesis by using aCGH and SNP-arrays. One of the most important goals was to screen cryptic alterations in karyotypically normal leukemia patients. In addition, chromosomal changes were evaluated to narrow the target regions, to find new markers, and to obtain tumor suppressor and oncogene candidates. The work presented here shows the capability of aCGH to detect submicroscopic copy number alterations in leukemia, with information about breakpoints and genes involved in the alterations, and that genome-wide microarray analyses with aCGH and SNP-array are advantageous methods in the research and diagnosis of leukemia. The most important findings were the cryptic changes detected with aCGH in karyotypically normal AML and CLL, characterization of amplified genes in 11q marker chromosomes, detection of deletion-based mechanisms of MLL-ARHGEF12 fusion gene formation, and detection of LOH without copy number alteration in karyotypically normal AML. These alterations harbor candidate oncogenes and tumor suppressors for further studies.

Stanniocalcin-1 in cell stress and differentiation

Stanniocalcin-1 (STC-1) is a 56 kD homodimeric protein which was originally identified in bony fish, where it regulates calcium/phosphate homeostasis and protects against toxic hypercalcemia. STC-1 was considered unique to fish until the cloning of cDNA for human STC-1 in 1995 and mouse Stc-1 in 1996. STC-1 is conserved through evolution with human and salmon STC-1 sharing 60% identity and 80% similarity. The surprisingly high homology between mammalian and fish STC-1 and the protective actions of STC-1 in terminally differentiated neurons, originally reported by my colleagues, prompted me to further study the role of STC-1 in cell stress and differentiation. One purpose was to determine whether there is an inter-relationship between terminally differentiated cells and STC-1 expression. The study revealed an accumulation of STC-1 in mature megakaryocytes and adipocytes, i.e. postmitotic cells with limited or lost proliferative capacity. Still proliferating uninduced cells were negative for STC-1 mRNA and protein, whereas differentiating cells accumulated STC-1 in their cytoplasm. Interestingly, in liposarcomas the grade inversely correlated with STC-1 expression. Another aim was to study how STC-1 gene expression is regulated. Given that IL-6 is a cytokine with neuroprotective actions, by unknown mechanisms, we examined whether IL-6 regulates STC-1 gene expression. Treatment of human neural Paju cells with IL-6 induced a dose-dependent upregulation of STC-1 mRNA levels. This induction of STC-1 expression by IL-6 occurred mainly through the MAPK signaling pathway. Furthermore, I studied the role of IL-6-mediated STC-1 expression as a mechanism of cytoprotection conferred by hypoxic preconditioning (HOPC) in brain and heart. My findings show that Stc-1 was upregulated in brain after hypoxia treatment. In the brain of IL-6 deficient mice, however, no upregulation of Stc-1 expression was evident. After induced brain injury the STC-1 response in brains of IL-6 transgenic mice, with IL-6 overexpression in astroglial cells, was stronger than in brains of WT mice. These results indicate that IL-6-mediated expression of STC-1 is one molecular mechanism of HOPC-induced tolerance to brain ischemia. The protection conferred by HOPC in heart occurs during a bimodal time course comprising early and delayed preconditioning. Interestingly, my results showed that the expression of Stc-1 in heart was upregulated in a biphasic manner during HOPC. IL-6 deficient mice did not, however, show a similar biphasic manner of Stc-1 upregulation as did WT mice. Instead, only an early upregulation of Stc-1 expression was evident. The results suggest that the upregulation of Stc-1 during the delayed preconditioning is IL-6-dependent. The upregulated expression of Stc-1 during the early preconditioning, however, is only partly IL-6-dependent and possibly also directly mediated by HIF-1. These findings suggest that STC-1 is a pro-survival protein for terminally differentiated cells and that STC-1 expression may in fact be regulated by stress. In addition, I show that STC-1 gene upregulation, mediated in part by IL-6, is a new mechanism of protection conferred by HOPC in brain and heart. Because of its importance for fundamental biological processes, such as differentiation and cytoprotection, STC-1 may have therapeutic implications for management of stroke, neurodegenerative diseases, cancer, and obesity.

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.

Novel insights on functions of myotilin/palladin family members

Poikkijuovaisen luuranko- ja sydänlihaksen supistumisyksikkö, sarkomeeri, koostuu tarkoin järjestyneistä aktiini- ja myosiinisäikeistä. Rakenne eroaa muista solutyypeistä, joissa aktiinisäikeistö muovautuu jatkuvasti ja sen järjestyminen säätelee solun muotoa, solujakautumista, soluliikettä ja solunsisäisten organellien kuljetusta. Myotilin, palladin ja myopalladin kuuluvat proteiiniperheeseen, jonka yhteispiirteenä ovat immunoglobuliinin kaltaiset (Igl) domeenit. Proteiinit liittyvät aktiinitukirankaan ja niiden arvellaan toimivan solutukirangan rakenne-elementteinä ja säätelijöinä. Myotilinia ja myopalladinia ilmennetään poikkijuovaisessa lihaksessa. Sen sijaan palladinin eri silmukointimuotoja tavataan monissa kudostyypeissä kuten hermostossa, ja eri muodoilla saattaa olla solutyypistä riippuvia tehtäviä. Poikkijuovaisessa lihaksessa kaikki perheen jäsenet sijaitsevat aktiinisäikeitä yhdistävässä Z-levyssä ja ne sitovat Z-levyn rakenneproteiinia, -aktiniinia. Myotilingeenin pistemutaatiot johtavat periytyviin lihastauteihin, kun taas palladinin mutaatioiden on kuvattu liittyvän periytyvään haimasyöpään ja lisääntyneeseen sydäninfarktin riskiin. Tässä tutkimuksessa selvitettin myotilinin ja pallainin toimintaa. Kokeissa löydettiin uusia palladinin 90-92kDa alatyyppiin sitoutuvia proteiineja. Yksi niistä on aktiinidynamiikkaa säätelevä profilin. Profilinilla on kahdenlaisia tehtäviä; se edesauttaa aktiinisäikeiden muodostumista, mutta se voi myös eristää yksittäisiä aktiinimolekyylejä ja edistää säikeiden hajoamista. Solutasolla palladinin ja profilinin sijainti on yhtenevä runsaasti aktiinia sisältävillä solujen reuna-alueilla. Palladinin ja profilinin sidos on heikko ja hyvin dynaaminen, joka sopii palladinin tehtävään aktiinisäideiden muodostumisen koordinoijana. Toinen palladinin sitoutumiskumppani on aktiinisäikeitä yhteensitova -aktiniini. -Aktiniini liittää solutukirangan solukalvon proteiineihin ja ankkuroi solunsisäisiä viestintämolekyylejä. Sitoutumista välittävä alue on hyvin samankaltainen palladinissa ja myotilinissa. Luurankolihaksen liiallinen toistuva venytys muuttaa Z-levyjen rakennetta ja muotoa. Prosessin aikana syntyy uusia aktiinifilamenttejä sisältäviä tiivistymiä ja lopulta uusia sarkomeereja. Löydöstemme perusteella myotilinin uudelleenjärjestyminen noudattaa aktiinin muutoksia. Tämä viittaa siihen, että myotilin liittää yhteen uudismuodostuvia aktiinisäikeitä ja vakauttaa niitä. Myotilin saattaa myös ankkuroida viesti- tai rakennemolekyylejä, joiden tehtävänä on edesauttaa Z-levyjen uudismuodostusta. Tulostemme perusteella arvelemme, että myotilin toimii Z-levyjen rakenteen vakaajana ja aktiinisäikeiden säätelijänä. Palladinin puute johtaa sikiöaikaiseen kuolemaan hiirillä, mutta myotilinin puutoksella ei ole samanlaisia vaikutuksia. Tuotettujen myotilin poistogeenisten hiirten todetiin syntyvän ja kehittyvän normaalisti eikä niillä esiintynyt rakenteellisia tai toiminnallisia häiriöitä. Toisaalta aiemmissa kokeissa, joissa hiirille on siirretty ihmisen lihastautia aikaansaava myotilingeeni, nähdään samankaltaisia kuin sairailla ihmisillä. Näin ollen muuntunut myotilin näyttä olevan lihaksen toiminnalle haitallisempi kuin myotilinin puute. Myotilinin ja palladinin yhteisvaikutusta selvittääksemme risteytimme myotilin poistegeenisen hiiren ja hiirilinjan, joka ilmentää puutteellisesti palladinin 200 kDa muotoa. Puutteellisesti 200 kDa palladinia ilmentävien hiirten sydänlihaksessa todettiin vähäisiä hienorakenteen muutoksia, mutta risteytetyillä hiirillä tavattiin rakenteellisia ja toiminnallisia muutoksia myös luurankolihaksessa. Tulosten perusteella voidaan todeta, että palladinin 200 kDa muoto säätelee sydänlihassolujen rakennetta. Luurankolihaksessa sen sijaan myotilinilla ja palladinilla näyttäisi olevan päällekkäisiä tehtäviä.

Expression Analysis of Host-Induced Genes and Proteins of Potato Pathogen Pectobacterium atrosepticum

Pectobacterium atrosepticum on Gram-negatiivinen bakteeri, joka aiheuttaa perunan tyvi- ja märkämätää. P. atrosepticum bakteerin optimilämpötila on melko alhainen ja se on yleinen lauhkeilla alueilla. Tyvimätä leviää pääasiassa siemenperunan välityksellä ja siksi se on ongelma erityisesti siemenperunan tuotannossa. P. atrosepticum kannan SCRI1043 genomi on julkaistu ja sitä tutkitaan malliorganismina märkä- ja tyvimädän taudinaiheuttamisen ymmärtämiseksi. Tämä opportunistinen taudinaiheuttaja voi elää isäntäkasvissa kuukausia piilevänä, aiheuttamatta näkyviä oireita. Suotuisissa olosuhteissa bakteerit alkavat jakautua ja tuottaa kasvin kudoksia hajottavia entsyymejä. Mädäntyvä kasvimassa tarjoaa ravinteita bakteerien kasvuun ja mahdollistaa isäntäkasvin asuttamisen. Soluseiniä hajottavien entsyymien merkitys taudinaiheuttamisessa on hyvin tunnettu, mutta oireettomasta jaksosta ja taudin alkuvaiheista tiedätään vain vähän. Bakteerin genomi sisältää monia toksiineja, adhesiineja, hemolysiineja ja muita proteiineja, joilla saattaa olla merkitys taudinaiheuttamisessa. Tässä työssä käytettiin proteomiikkaa ja mikrosiruanalysiä P. atrosepticum bakteerin erittyvien proteiinien ja geeniekspression tutkimiseen. Proteiinit, jotka eritetään ulos bakteerista, toimivat todennäköisesti taudinaiheuttamisessa, koska ne ovat suorassa kontaktissa isäntäkasvin kanssa. Analyysit suoritettiin olosuhteissa, jotka muistuttavat kasvin soluvälitilaa: matala pH, vähän ravinteita ja matala lämpötila. Isäntäkasvin läsnäolon vaikutusta proteiinien tuottoon ja geeniekspressioon tutkittiin lisäämällä perunauutetta kasvatusalustaan. Tutkimuksessa tunnistettiin P. atrosepticum bakteerin monia jo tunnettuja ja mahdollisesti taudinaiheuttamiseen liittyviä proteiineja. Perunauute lisäsi hiljattain tunnistetun, proteiinien eritysreittiä (tyyppi VI sekreetio, T6SS) koodaavien geenien ilmentymistä. Lisäksi bakteerin havaittiin erittävän useita T6SS:n liittyviä proteiineja kasvualustaan, johon oli lisätty perunauutetta. T6SS:n merkitys bakteereille on vielä epäselvä ja sen vaikutuksesta taudinaiheuttamiseen on julkaistu ristiriitaisia tuloksia. Märkä- ja tyvimädän ymmärtäminen molekulaarisella tasolla luo pohjan tautien kontrollointiin tähtäävään soveltavaan tutkimukseen. Tämä tutkimus lisää tietoa kasvi-patogeeni- interaktiosta ja sitä voidaan tulevaisuudessa käyttää hyväksi esimerkiksi diagnostiikassa, resistenttien perunalajikkeiden jalostuksessa tai viljely- ja varastointiolosuhteiden parantamisessa.

Genetic diversity of Atlantic salmon (Salmo salar) and European whitefish (Coregonus lavaretus) in the Baltic Sea

We described the patterns and extent of microsatellite DNA variation in historical and present-day Atlantic salmon (Salmo salar L.) stocks in the Baltic Sea and neighbouring areas, and in European whitefish (Coregonus lavaretus) ecotypes, populations and run-timing types in Finland. Moreover, the amount and pattern of genetic diversity in historical salmon populations before human impact were described, and the proportion of diversity maintained in the present hatchery stocks evaluated. Salmon populations in the Baltic Sea were, on average, significantly less variable than eastern Atlantic populations, and the diversity of landlocked populations (Lakes Vänern, Saimaa, Onega and Ladoga) was in turn significantly lower than that of anadromous salmon populations in the Baltic Sea populations. Within the Baltic Sea, the anadromous populations of Atlantic salmon formed three clear groups, corresponding to the northern (Gulf of Bothnia), eastern (Gulf of Finland and eastern Baltic Main Basin) and southern (western Baltic Main Basin) regions. Based on microsatellite data, three salmon population groups in the Baltic Sea were considered potentially different colonization lineages. In short- and long-term breeding programmes of Atlantic salmon, the average observed rate of loss of alleles was 4.9% and 2.0% per generation and the average rate of loss of heterozygosity was 1.4% and 1% per generation, respectively. When comparing the genetic parameters of stocks before and after hatchery breeding of several successive generations (Rivers Iijoki and Oulujoki), statistically significant changes in allele frequencies were common, while large wild stock in the Teno River has remained temporally very stable over 56 years. Despite the observed losses of genetic diversity in broodstock breeding, a large proportion of the genetic resources of the extirpated stocks are still conserved in the broodstocks. Genetic differentiation among European whitefish ecotypes was generally low, thus giving support to the hypothesis of one native European whitefish species in Fennoscandia. Among the ecotypes, the northern, large sparsely rakered, bottom-dwelling whitefish was the most unique. The known genetic differences in quantitative traits have thus either developed independently of potential phylogenetic lineages, or the lineages have mixed and the quantitative traits of the ecotypes, like gill-raker number, have later changed according to environment and selection pressures. Overall, genetic distances between the anadromous whitefish populations along the Finnish coast, especially in the Bothnian Bay area, were small. Wild whitefish populations studied had slightly higher allelic diversity than hatchery-reared populations in corresponding rivers.