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.

Genetics of T cell co-stimulatory receptors -CD28, CTLA4, ICOS and PDCD1 in immunity and transplantation

Co-stimulatory signals are essential for the activation of naïve T cells and productive immune response. Naïve T cells receive first, antigen-specific signal through T cell receptor. Co-stimulatory receptors provide the second signal which can be either activating or inhibitory. The balance between signals determines the outcome of an immune response. CD28 is crucial for T cell activation; whereas cytotoxic T lymphocyte associated antigen 4 (CTLA4) mediates critical inhibitory signal. Inducible co-stimulator (ICOS) augments cytokine expression and plays role in immunoglobulin class switching. Programmed cell death 1 (PDCD1) acts as negative regulator of T cell proliferation and cytokine responses. The co-stimulatory receptor pathways are potentially involved in self-tolerance and thus, they provide a promising therapeutic strategy for autoimmune diseases and transplantation. The genes encoding CD28, CTLA4 and ICOS are located adjacently in the chromosome region 2q33. The PDCD1 gene maps further, to the region 2q37. CTLA4 and PDCD1 are associated with the risk of a few autoimmune diseases. There is strong linkage disequilibrium (LD) on the 2q33 region; the whole gene of CD28 exists in its own LD block but CTLA4 and the 5' part of ICOS are within a same LD block. The 3' part of ICOS and PDCD1 are in their own separate LD blocks. Extended haplotypes covering the 2q33 region can be identified. This study focuses on immune related conditions like coeliac disease (CD) which is a chronic inflammatory disease with autoimmune features. Immunoglobulin A deficiency (IgAD) belongs to the group of primary antibody deficiencies characterised by reduced levels of immunoglobulins. IgAD co-occurs often with coeliac disease. Renal transplantation is needed in the end stage kidney diseases. Transplantation causes strong immune response which is tried to suppress with drugs. All these conditions are multifactorial with complex genetic background and multiple environmental factors affecting the outcome. We have screened ICOS for polymorphisms by sequencing the exon regions. We detected 11 new variants and determined their frequencies in Finnish population. We have measured linkage disequilibrium on the 2q33 region in Finnish as well as other European populations and observed conserved haplotypes. We analysed genetic association and linkage of the co-stimulatory receptor gene region aiming to study if it is a common risk locus for immune diseases. The 2q33 region was replicated to be linked to coeliac disease in Finnish population and CTLA4-ICOS haplotypes were found to be associated with CD and IgAD being the first non-HLA risk locus common for CD and immunodeficiencies. We also showed association between ICOS and the outcome of kidney transplantation. Our results suggest new evidence for CTLA4-ICOS gene region to be involved in susceptibility of coeliac disease. The earlier published contradictory association results can be explained by involvement of both CTLA4 and ICOS in disease susceptibility. The pattern of variants acting together rather than a single polymorphism may confer the disease risk. These genes may predispose also to immunodeficiencies as well as decreased graft survival and delayed graft function. Consequently, the present study indicates that like the well established HLA locus, the co-stimulatory receptor genes predispose to variety of immune disorders.

Rhizoctonia -Scots pine interactions: detection, impact on seedling performance and host defence gene response

Rhizoctonia spp. are ubiquitous soil inhabiting fungi that enter into pathogenic or symbiotic associations with plants. In general Rhizoctonia spp. are regarded as plant pathogenic fungi and many cause root rot and other plant diseases which results in considerable economic losses both in agriculture and forestry. Many Rhizoctonia strains enter into symbiotic mycorrhizal associations with orchids and some hypovirulent strains are promising biocontrol candidates in preventing host plant infection by pathogenic Rhizoctonia strains. This work focuses on uni- and binucleate Rhizoctonia (respectively UNR and BNR) strains belonging to the teleomorphic genus Ceratobasidium, but multinucleate Rhizoctonia (MNR) belonging to teleomorphic genus Thanatephorus and ectomycorrhizal fungal species, such as Suillus bovinus, were also included in DNA probe development work. Strain specific probes were developed to target rDNA ITS (internal transcribed spacer) sequences (ITS1, 5.8S and ITS2) and applied in Southern dot blot and liquid hybridization assays. Liquid hybridization was more sensitive and the size of the hybridized PCR products could be detected simultaneously, but the advantage in Southern hybridization was that sample DNA could be used without additional PCR amplification. The impacts of four Finnish BNR Ceratorhiza sp. strains 251, 266, 268 and 269 were investigated on Scot pine (Pinus sylvestris) seedling growth, and the infection biology and infection levels were microscopically examined following tryphan blue staining of infected roots. All BNR strains enhanced early seedling growth and affected the root architecture, while the infection levels remained low. The fungal infection was restricted to the outer cortical regions of long roots and typical monilioid cells detected with strain 268. The interactions of pathogenic UNR Ceratobasidium bicorne strain 1983-111/1N, and endophytic BNR Ceratorhiza sp. strain 268 were studied in single or dual inoculated Scots pine roots. The fungal infection levels and host defence-gene activity of nine transcripts [phenylalanine ammonia lyase (pal1), silbene synthase (STS), chalcone synthase (CHS), short-root specific peroxidase (Psyp1), antimicrobial peptide gene (Sp-AMP), rapidly elicited defence-related gene (PsACRE), germin-like protein (PsGER1), CuZn- superoxide dismutase (SOD), and dehydrin-like protein (dhy-like)] were measured from differentially treated and un-treated control roots by quantitative real time PCR (qRT-PCR). The infection level of pathogenic UNR was restricted in BNR- pre-inoculated Scots pine roots, while UNR was more competitive in simultaneous dual infection. The STS transcript was highly up-regulated in all treated roots, while CHS, pal1, and Psyp1 transcripts were more moderately activated. No significant activity of Sp-AMP, PsACRE, PsGER1, SOD, or dhy-like transcripts were detected compared to control roots. The integrated experiments presented, provide tools to assist in the future detection of these fungi in the environment and to understand the host infection biology and defence, and relationships between these interacting fungi in roots and soils. This study further confirms the complexity of the Rhizoctonia group both phylogenetically and in their infection biology and plant host specificity. The knowledge obtained could be applied in integrated forestry nursery management programmes.

Pathogenicity, functional significance and clinical phenotype of mismatch repair gene MSH2 variants found in cancer patients

DNA ja siinä sijaitsevat geenit ohjaavat kaikkea solujen toimintaa. DNA-molekyyleihin kuitenkin kertyy mutaatioita sekä ympäristön vaikutuksen, että solujen oman toiminnan tuloksena. Mikäli virheitä ei korjata, saattaa tuloksena olla solun muuttuminen syöpäsoluksi. Soluilla onkin käytössä useita DNA-virheiden korjausmekanismeja, joista yksi on ns. mismatch repair (MMR). MMR vastaa DNA:n kahdentumisessa syntyvien virheiden korjauksesta. Periytyvät mutaatiot geeneissä, jotka vastaavat MMR-proteiinien rakentamisesta, aiheuttavat ongelmia DNA:n korjauksessa ja altistavat kantajansa periytyvälle ei-polypoottiselle paksusuolisyöpäoireyhtymälle (hereditary nonpolyposis colorectal cancer, HNPCC). Yleisimmin mutatoituneet MMR-geenit ovat MLH1 ja MSH2. HNPCC periytyy vallitsevasti, eli jo toiselta vanhemmalta peritty geenivirhe altistaa syövälle. MMR-geenivirheen kantaja sairastuu syöpään elämänsä aikana suurella todennäköisyydellä, ja sairastumisikä on vain noin 40 vuotta. Syövälle altistavan geenivirheen löytäminen mutaation kantajilta on hyvin tärkeää, sillä säännöllinen seuranta mahdollistaa kehittymässä olevan kasvaimen havaitsemisen ja poistamisen jo aikaisessa vaiheessa. Tämän on osoitettu alentavan syöpäkuolleisuutta merkittävästi. Varma tieto altistuksen alkuperästä on tärkeä myös niille syöpäsuvun jäsenille, jotka eivät kanna kyseistä mutaatiota. Syövälle altistavien mutaatioiden ohella MMR-geeneistä löydetään säännöllisesti muutoksia, jotka ovat normaalia henkilöiden välistä geneettistä vaihtelua, eikä niiden oleteta lisäävän syöpäaltistusta. Altistavien mutaatioiden erottaminen näistä neutraaleista variaatioista on vaikeaa, mutta välttämätöntä altistuneiden tehokkaan seurannan varmistamiseksi. Tässä väitöskirjassa tutkittiin 18:a MSH2 -geenin mutaatiota. Mutaatiot oli löydetty perheistä, joissa esiintyi paljon syöpiä, mutta niiden vaikutus DNA:n korjaustehoon ja syöpäaltistukseen oli epäselvä. Työssä tutkittiin kunkin mutaation vaikutusta MSH2-proteiinin normaaliin toimintaan, ja tuloksia verrattiin potilaiden ja sukujen kliinisiin tietoihin. Tutkituista mutaatiosta 12 aiheutti puutteita MMR-korjauksessa. Nämä mutaatiot tulkittiin syövälle altistaviksi. Analyyseissä normaalisti toimineet 4 mutaatiota eivät todennäköisesti ole syynä syövän syntyyn kyseisillä perheillä. Tulkinta jätettiin avoimeksi 2 mutaation kohdalla. Tutkimuksesta hyötyivät suoraan kuvattujen mutaatioiden kantajaperheet, joiden geenivirheen syöpäaltistuksesta saatiin tietoa, mahdollistaen perinnöllisyysneuvonnan ja seurannan kohdentamisen sitä tarvitseville. Työ selvensi myös mekanismeja, joilla mutatoitunut MSH2-proteiini voi menettää toimintakykynsä.

Immune response to insulin and changes in the gut immune system in children with or at risk for type 1 diabetes

Type 1 diabetes (T1D) is considered to be an autoimmune disease. The cause of T1D is the destruction of insulin-producing β-cells in the pancreatic islets. The autoimmune nature of T1D is characterized by the presence of autoreactive T-cells and autoantibodies against β-cell molecules. Insulin is the only β-cell-specific autoantigen associated with T1D but the insulin autoantibodies (IAAs) are difficult to measure with proper sensitivity. T-cell assays for detection of autoreactive T-cells, such as insulin-specific T-cells, have also proven to be difficult to perform. The genetic risk of T1D is associated with the HLA gene region but the environmental factors also play an important role. The most studied environmental risk factors of T1D are enteroviruses and cow's milk which both affect the immune system through the gut. One hypothesis is that the insulin-specific immune response develops against bovine insulin in cow's milk during early infancy and later spreads to include human insulin. The aims of this study were to determine whether the separation of immunoglobulin (Ig)G from plasma would improve the sensitivity of the IAA assay and how insulin treatment affects the cellular immune response to insulin in newly diagnosed patients. Furthermore, the effect of insulin concentration in mother's breast milk on the development of antibodies to dietary insulin in the child was examined. Small intestinal biopsies were also obtained from children with T1D to characterize any immunological changes associated with T1D in the gut. The isolation of the IgG fraction from the plasma of T1D patients negative for plasma IAA led to detectable IAA levels that exceeded those in the control children. Thus the isolation of IgG may improve the sensitivity of the IAA assay. The effect of insulin treatment on insulin-specific T-cells was studied by culturing peripheral blood mononuclear cells with insulin. The insulin stimulation induced increased expression of regulatory T-cell markers, such as Foxp3, in those patients treated with insulin than in patients examined before initiating insulin treatment. This finding suggests that insulin treatment in patients with T1D stimulates regulatory T-cells in vivo and this may partly explain the difficulties in measuring autoantigen-specific T-cell responses in recently diagnosed patients. The stimulation of regulatory T-cells by insulin treatment may also explain the remission period often seen after initiating insulin treatment. In the third study we showed that insulin concentration in mother's breast milk correlates inversely with the levels of bovine insulin-specific antibodies in those infants who were exposed to cow's milk proteins in their diet, suggesting that human insulin in breast milk induces tolerance to dietary bovine insulin. However, in infants who later developed T1D-associated autoantibodies, the insulin concentration in their mother's breast milk was increased. This finding may indicate that in those children prone to β-cell autoimmunity, breast milk insulin does not promote tolerance to insulin. In the small intestinal biopsies the presence of several immunological markers were quantified with the RT-PCR. From these markers the expression of the interleukin (IL)-18 cytokine was significantly increased in the gut in patients with T1D compared with children with celiac disease or control children. The increased IL-18 expression lends further support for the hypothesis that the gut immune system is involved in the pathogenesis of T1D.

U12-type Spliceosome: Localization and Effects of Splicing Efficiency on Gene Expression

The removal of non-coding sequences, introns, is an essential part of messenger RNA processing. In most metazoan organisms, the U12-type spliceosome processes a subset of introns containing highly conserved recognition sequences. U12-type introns constitute less than 0,5% of all introns and reside preferentially in genes related to information processing functions, as opposed to genes encoding for metabolic enzymes. It has previously been shown that the excision of U12-type introns is inefficient compared to that of U2-type introns, supporting the model that these introns could provide a rate-limiting control for gene expression. The low efficiency of U12-type splicing is believed to have important consequences to gene expression by limiting the production of mature mRNAs from genes containing U12-type introns. The inefficiency of U12-type splicing has been attributed to the low abundance of the components of the U12-type spliceosome in cells, but this hypothesis has not been proven. The aim of the first part of this work was to study the effect of the abundance of the spliceosomal snRNA components on splicing. Cells with a low abundance of the U12-type spliceosome were found to inefficiently process U12-type introns encoded by a transfected construct, but the expression levels of endogenous genes were not found to be affected by the abundance of the U12-type spliceosome. However, significant levels of endogenous unspliced U12-type intron-containing pre-mRNAs were detected in cells. Together these results support the idea that U12-type splicing may limit gene expression in some situations. The inefficiency of U12-type splicing has also promoted the idea that the U12-type spliceosome may control gene expression, limiting the mRNA levels of some U12-type intron-containing genes. While the identities of the primary target genes that contain U12-type introns are relatively well known, little has previously been known about the downstream genes and pathways potentially affected by the efficiency of U12-type intron processing. Here, the effects of U12-type splicing efficiency on a whole organism were studied in a Drosophila line with a mutation in an essential U12-type spliceosome component. Genes containing U12-type introns showed variable gene-specific responses to the splicing defect, which points to variation in the susceptibility of different genes to changes in splicing efficiency. Surprisingly, microarray screening revealed that metabolic genes were enriched among downstream effects, and that the phenotype could largely be attributed to one U12-type intron-containing mitochondrial gene. Gene expression control by the U12-type spliceosome could thus have widespread effects on metabolic functions in the organism. The subcellular localization of the U12-type spliceosome components was studied as a response to a recent dispute on the localization of the U12-type spliceosome. All components studied were found to be nuclear indicating that the processing of U12-type introns occurs within the nucleus, thus clarifying a question central to the field. The results suggest that the U12-type spliceosome can limit the expression of genes that contain U12-type introns in a gene-specific manner. Through its limiting role in pre-mRNA processing, the U12-type splicing activity can affect specific genetic pathways, which in the case of Drosophila are involved in metabolic functions.

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.

Microbe-induced Innate Immune Responses in Human Dendritic Cells

Two types of antigen-presenting cells (APCs), macrophages and dendritic cells (DCs), function at the interface of innate and adaptive immunity. Through recognition of conserved microbial patterns, they are able to detect the invading pathogens. This leads to activation of signal transduction pathways that in turn induce gene expression of various molecules required for immune responses and eventually pathogen clearance. Cytokines are among the genes induced upon detection of microbes. They play an important role in regulating host immune responses during microbial infection. Chemotactic cytokines, chemokines, are involved in migratory events of immune cells. Cytokines also promote the differentiation of distinct T cell responses. Because of the multiple roles of cytokines in the immune system, the cytokine network needs to be tightly regulated. In this work, the induction of innate immune responses was studied using human primary macrophages or DCs as cell models. Salmonella enterica serovar Typhimurium served as a model for an intracellular bacterium, whereas Sendai virus was used in virus experiments. The starting point of this study was that DCs of mouse origin had recently been characterized as host cells for Salmonella. However, only little was known about the immune responses initiated in Salmonella-infected human DCs. Thus, cellular responses of macrophages and DCs, in particular the pattern of cytokine production, to Salmonella infection were compared. Salmonella-induced macrophages and DCs were found to produce multiple cytokines including interferon (IFN) -gamma, which is conventionally produced by T and natural killer (NK) cells. Both macrophages and DCs also promoted the intracellular survival of the bacterium. Phenotypic maturation of DCs as characterized by upregulation of costimulatory and human leukocyte antigen (HLA) molecules, and production of CCL19 chemokine, were also detected upon infection with Salmonella. Another focus of this PhD work was to unravel the regulatory events controlling the expression of cytokine genes encoding for CCL19 and type III IFNs, which are central to DC biology. We found that the promoters of CCL19 and type III IFNs contain similar regulatory elements that bind nuclear factor kappaB (NF-kappaB) and interferon regulatory factors (IRFs), which could mediate transcriptional activation of the genes. The regulation of type III IFNs in virus infection resembled that of type I IFNs a cytokine class traditionally regarded as antiviral. The induction of type I and type III IFNs was also observed in response to bacterial infection. Taken together, this work identifies new details about the interaction of Salmonella with its phagocytic host cells of human origin. In addition, studies provide information on the regulatory events controlling the expression of CCL19 and the most recently identified IFN family genes, type III IFN genes.

Type III Secretion System of Phytopathogenic Bacterium Pseudomonas syringae:From Gene to Function

The type III secretion system (T3SS) is an essential requirement for the virulence of many Gram-negative bacteria which infect plants, animals and men. Pathogens use the T3SS to deliver effector proteins from the bacterial cytoplasm to the eukaryotic host cells, where the effectors subvert host defenses. The best candidates for directing effector protein traffic are the bacterial type III-associated appendages, called needles or pili. In plant pathogenic bacteria, the best characterized example of a T3SS-associated appendage is the HrpA pilus of the plant pathogen Pseudomonas syringae pv. tomato DC3000. The components of the T3SS in plant pathogens are encoded by a cluster of hrp (hypersensitive reaction and pathogenicity) genes. Two major classes of T3SS-secreted proteins are: harpin proteins such as HrpZ which are exported into extracellular space, and avirulence (Avr) proteins such as AvrPto which are translocated directly to the plant cytoplasm. This study deals with the structural and functional characterization of the T3SS-associated HrpA pilus and the T3SS-secreted harpins. By insertional mutagenesis analysis of HrpA, we located the optimal epitope insertion site in the amino-terminus of HrpA, and revealed the potential application of the HrpA pilus as a carrier of antigenic determinants for vaccination. By pulse-expression of proteins combined with immuno-electron microscopy, we discovered the Hrp pilus assembly strategy as addition of HrpA subunits to the distal end of the growing pilus, and we showed for the first time that secretion of HrpZ occurs at the tip of the pilus. The pilus thus functions as a conduit delivering proteins to the extracellular milieu. By using phage-display and scanning-insertion mutagenesis methods we identified a conserved HrpZ-binding peptide and localized the peptide-binding site to the central domain of HrpZ. We also found that the HrpZ specifically interacts with a host bean protein. Taken together, the current results provide deeper insight into the molecular mechanism of T3SS-associated pilus assembly and effector protein translocation, which will be helpful for further studies on the pathogenic mechanisms of Gram-negative bacteria and for developing new strategies to prevent bacterial infection.

FGFR1 regulated gene-expression, cell proliferation and differentiation in the developing midbrain and hindbrain

The neuroectodermal tissue close to the midbrain hindbrain boundary (MHB) is an important secondary organizer in the developing neural tube. This so-called isthmic organizer (IsO) regulates cellular survival, patterning and proliferation in the midbrain (Mb) and rhombomere 1 (R1) of the hindbrain. Signaling molecules of the IsO, such as fibroblast growth factor 8 (FGF8) and WNT1 are expressed in distinct bands of cells around the MHB. It has been previously shown that FGF-receptor 1 (FGFR1) is required for the normal development of this brain region in the mouse embryo. In the present study, we have compared the gene expression profiles of wild-type and Fgfr1 mutant embryos. We show that the loss of Fgfr1 results in the downregulation of several genes expressed close to the MHB and in the disappearance of gene expression gradients in the midbrain and R1. Our microarray screen identified several previously uncharacterized genes which may participate in the development of midbrain R1 region. Our results also show altered neurogenesis in the midbrain and R1 of the Fgfr1 mutants. Interestingly, the neuronal progenitors in midbrain and R1 show different responses to the loss of signaling through FGFR1. As Wnt1 expression at the MHB region requires the FGF signaling pathway, WNT target genes, including Drapc1, were also identified in our screen. The microarray data analysis also suggested that the cells next to the midbrain hindbrain boundary express distinct cell cycle regulators. We showed that the cells close to the border appeared to have unique features. These cells proliferate less rapidly than the surrounding cells. Unlike the cells further away from the boundary, these cells express Fgfr1 but not the other FGF receptors. The slowly proliferating boundary cells are necessary for development of the characteristic isthmic constriction. They may also contribute to compartmentalization of this brain region.

Applications of gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other teleost fish species

Evolutionary genetics incorporates traditional population genetics and studies of the origins of genetic variation by mutation and recombination, and the molecular evolution of genomes. Among the primary forces that have potential to affect the genetic variation within and among populations, including those that may lead to adaptation and speciation, are genetic drift, gene flow, mutations and natural selection. The main challenges in knowing the genetic basis of evolutionary changes is to distinguish the adaptive selection forces that cause existent DNA sequence variants and also to identify the nucleotide differences responsible for the observed phenotypic variation. To understand the effects of various forces, interpretation of gene sequence variation has been the principal basis of many evolutionary genetic studies. The main aim of this thesis was to assess different forms of teleost gene sequence polymorphisms in evolutionary genetic studies of Atlantic salmon (Salmo salar) and other species. Firstly, the level of Darwinian adaptive evolution affected coding regions of the growth hormone (GH) gene during the teleost evolution was investigated based on the sequence data existing in public databases. Secondly, a target gene approach was used to identify within population variation in the growth hormone 1 (GH1) gene in salmon. Then, a new strategy for single nucleotide polymorphisms (SNPs) discovery in salmonid fishes was introduced, and, finally, the usefulness of a limited number of SNP markers as molecular tools in several applications of population genetics in Atlantic salmon was assessed. This thesis showed that the gene sequences in databases can be utilized to perform comparative studies of molecular evolution, and some putative evidence of the existence of Darwinian selection during the teleost GH evolution was presented. In addition, existent sequence data was exploited to investigate GH1 gene variation within Atlantic salmon populations throughout its range. Purifying selection is suggested to be the predominant evolutionary force controlling the genetic variation of this gene in salmon, and some support for gene flow between continents was also observed. The novel approach to SNP discovery in species with duplicated genome fragments introduced here proved to be an effective method, and this may have several applications in evolutionary genetics with different species - e.g. when developing gene-targeted markers to investigate quantitative genetic variation. The thesis also demonstrated that only a few SNPs performed highly similar signals in some of the population genetic analyses when compared with the microsatellite markers. This may have useful applications when estimating genetic diversity in genes having a potential role in ecological and conservation issues, or when using hard biological samples in genetic studies as SNPs can be applied with relatively highly degraded DNA.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.