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.

ORP1L, a new Rab7 effector and regulator of late endosome functions

Oxysterol binding protein (OSBP) homologues have been found in eukaryotic organisms ranging from yeast to humans. These evolutionary conserved proteins have in common the presence of an OSBP-related domain (ORD) which contains the fully conserved EQVSHHPP sequence motif. The ORD forms a barrel structure that binds sterols in its interior. Other domains and sequence elements found in OSBP-homologues include pleckstrin homology domains, ankyrin repeats and two phenylalanines in an acidic tract (FFAT) motifs, which target the proteins to distinct subcellular compartments. OSBP homologues have been implicated in a wide range of intracellular processes, including vesicle trafficking, lipid metabolism and cell signaling, but little is known about the functional mechanisms of these proteins. The human family of OSBP homologues consists of twelve OSBP-related proteins (ORP). This thesis work is focused on one of the family members, ORP1, of which two variants were found to be expressed tissue-specifically in humans. The shorter variant, ORP1S contains an ORD only. The N-terminally extended variant, ORP1L, comprises a pleckstrin homology domain and three ankyrin repeats in addition to the ORD. The two ORP1 variants differ in intracellular localization. ORP1S is cytosolic, while the ankyrin repeat region of ORP1L targets the protein to late endosomes/lysosomes. This part of ORP1L also has profound effects on late endosomal morphology, inducing perinuclear clustering of late endosomes. A central aim of this study was to identify molecular interactions of ORP1L on late endosomes. The morphological changes of late endosomes induced by overexpressed ORP1L implies involvement of small Rab GTPases, regulators of organelle motility, tethering, docking and/or fusion, in generation of the phenotype. A direct interaction was demonstrated between ORP1L and active Rab7. ORP1L prolongs the active state of Rab7 by stabilizing its GTP-bound form. The clustering of late endosomes/lysosomes was also shown to be linked to the minus end-directed microtubule-based dynein-dynactin motor complex through the ankyrin repeat region of ORP1L. ORP1L, Rab7 and the Rab7-interacting lysosomal protein (RILP) were found to be part of the same effector complex recruiting the dynein-dynactin complex to late endosomes, thereby promoting minus end-directed movement. The proteins were found to be physically close to each other on late endosomes and RILP was found to stabilize the ORP1L-Rab7 interaction. It is possible that ORP1L and RILP bind to each other through their C-terminal and N-terminal regions, respectively, when they are bridged by Rab7. With the results of this study we have been able to place a member of the uncharacterized OSBP-family, ORP1L, in the endocytic pathway, where it regulates motility and possibly fusion of late endosomes through interaction with the small GTPase Rab7.

Short-term effects of variable retention on epigaeic spiders and carabid beetles in Finland

Forestry has influenced forest dwelling organisms for centuries in Fennoscandia. For example, in Finland ca. 30% of the threatened species are threatened because of forestry. Nowadays forest management recommendations include practices aimed at maintaining biodiversity in harvesting, such as green-tree retention. However, the effects of these practices have been little studied. In variable retention, different numbers of trees are retained, varying from green-tree retention (at least a few live standing trees in clear-cuts) to thinning (only individual trees removed). I examined the responses of ground-dwelling spiders and carabid beetles to green-tree retention (with small and large tree groups), gap felling and thinning aimed at an uneven age structure of trees. The impacts of these harvesting methods were compared to those of clear-cutting and uncut controls. I aimed to test the hypothesis that retaining more trees positively affects populations of those species of spiders and carabids that were present before harvesting. The data come from two studies. First, spiders were collected with pitfall traps in south-central Finland in 1995 (pre-treatment) and 1998 (after-treatment) in order to examine the effects of clear-cutting, green-tree retention (with 0.01-0.02-ha sized tree groups), gap felling (with three 0.16-ha sized openings in a 1-ha stand), thinning aiming at an uneven age structure of trees and uncut control. Second, spiders and carabids were caught with pitfall traps in eastern Finland in 1998-2001 (pre-treatment and three post-treatment years) in eleven 0.09-0.55-ha sized retention-tree groups and clear-cuts adjacent to them. Original spider and carabid assemblages were better maintained after harvests that retained more trees. Thinning maintained forest spiders well. However, gap felling and large retention-tree groups maintained some forest spider and carabid species in the short-term, but negatively affected some species over time. However, use of small retention-tree groups was associated with negative effects on forest spider populations. Studies are needed on the long-term effects of variable retention on terrestrial invertebrates; especially those directed at defining appropriate retention patch size and on the importance of structural diversity provided by variable retention for invertebrate populations. However, the aims of variable retention should be specified first. For example, are retention-tree groups planned to constitute life-boats , stepping-stones or to create structural diversity? Does it suffice that some species are maintained, or do we want to preserve the most sensitive ones, and how are these best defined? Moreover, the ecological benefits and economic costs of modified logging methods should be compared to other approaches aimed at maintaining biodiversity.

Phylogeography and Population Genetics of the Moor frog (Rana arvalis Nilsson) in Northern Europe

The ongoing climate change along with increasing levels of pollutants, diseases, habitat loss and fragmentation constitute global threats to the persistence of many populations, species and ecosystems. However, for the long-term persistence of local populations, one of the biggest threats is the intrinsic loss of genetic variation. In order to adapt to changes in the environment, organisms must have a sufficient supply of heritable variation in traits important for their fitness. With a loss of genetic variation, the risk of extinction will increase. For conservational practices, one should therefore understand the processes that shape the genetic population structure and also the broader (historical) phylogenetic patterning of the species in focus. In this thesis, microsatellite markers were applied to study genetic diversity and population differentiation of the protected moor frog (Rana arvalis) in Fennoscandia from both historical (evolutionary) and applied (conservation) perspectives. The results demonstrate that R. arvalis populations are highly structured over rather short geographic distances. Moreover, the results suggest that R. arvalis recolonized Fennoscandia from two directions after the last ice age. This has had implications for the genetic structuring and population differentiation, especially in the northernmost parts where the two lineages have met. Compared to more southern populations, the genetic variation decreases and the interpopulation differentiation increases dramatically towards north. This could be an outcome of serial population bottlenecking along the recolonization route. Also, current isolation and small population sizes increase the effect of drift, thus reinforcing the observed pattern. The same pattern can also be seen in island populations. However, though R. arvalis on the island of Gotland has lost most of its neutral genetic variability, our results indicate that the levels of additive genetic variation have remained high. This conforms to the conjecture that though neutral markers are widely used in conservation purposes, they may be quite uninformative about the levels of genetic variation in ecologically important traits. Finally, the evolutionary impact of the typical amphibian mating behaviour on genetic diversity was investigated. Given the short time available for larval development, it is important that mating takes place as early as possible. The genetic data and earlier capture-recapture data suggest that R. arvalis gather at mating grounds they are familiar with. However, by forming leks in random to relatedness, and having multiple paternities in single clutches, the risk of inbreeding may be minimized in this otherwise highly philopatric species.

Bioavailability and toxicity of chemicals in inland waters : the importance of prevailing water chemistry and implications for risk assessment

Physicochemical characterization of freshwater samples from Finland, Sweden, the Netherlands, and Spain revealed that water hardness and pH decreased and the quantity and quality of humic substances changed considerably in this geographical series from south to north. Since the ambient water chemistry may affect the availability of chemicals, the total aqueous concentration of a chemical may be insufficient to predict the bioconcentration, subsequent biological response, and thus risk. In addition, organisms could be affected directly by water quality characteristics. In this context the main objective of this thesis was to investigate the bioavailability of selected ecotoxicologically relevant chemicals (cadmium, benzo(a)pyrene, and pyrene) in various European surface waters and to show the importance of certain water chemistry characteristics in interpreting the bioavailability and toxicity results. The bioavailability of cadmium to Daphnia magna was examined in very soft humic lake water. Humic substances as natural ligands decreased the free and bioavailable proportion of cadmium in soft lake water. As a consequence the uptake rate and the acute toxicity decreased compared with the humic-free reference. When the hardness of humic lake water was artificially elevated, the acute toxicity of cadmium decreased, although the proportion of free cadmium increased. The decreased bioavailability of cadmium in hard water was a result of effective competition for uptake by the hardness cations, especially calcium ions. The protective role of humic substances and water hardness against cadmium toxicity was also observed in Lumbriculus variegatus, although D. magna was more sensitive to cadmium. The bioavailability of two polycyclic aromatic hydrocarbons (PAHs), pyrene and benzo(a)pyrene, was studied in European surface waters of varying water chemistry. Humic substances acted as complexing ligands with both PAHs, but the bioavailability of the more lipophilic benzo(a)pyrene to D. magna was affected more by humic substances than that of pyrene. In addition, not only the quantity of humic substances, but also their quality affected the bioavailability of benzo(a)pyrene. Nevertheless, the humic substances played a protective role in the photo-enhanced toxicity of pyrene under UV-B radiation. Water hardness had no effect on pyrene toxicity. Results indicate that the typical physicochemical characteristics of boreal freshwaters should be considered carefully in local and regional risk assessment of chemicals concerning the Fennoscandian region.

Redox-linked proton transfer by cytochrome c oxidase

The respiratory chain is found in the inner mitochondrial membrane of higher organisms and in the plasma membrane of many bacteria. It consists of several membrane-spanning enzymes, which conserve the energy that is liberated from the degradation of food molecules as an electrochemical proton gradient across the membrane. The proton gradient can later be utilized by the cell for different energy requiring processes, e.g. ATP production, cellular motion or active transport of ions. The difference in proton concentration between the two sides of the membrane is a result of the translocation of protons by the enzymes of the respiratory chain, from the negatively charged (N-side) to the positively charged side (P-side) of the lipid bilayer, against the proton concentration gradient. The endergonic proton transfer is driven by the flow of electrons through the enzymes of the respiratory chain, from low redox-potential electron donors to acceptors of higher potential, and ultimately to oxygen. Cytochrome c oxidase is the last enzyme in the respiratory chain and catalyzes the reduction of dioxygen to water. The redox reaction is coupled to proton transport across the membrane by a yet unresolved mechanism. Cytochrome c oxidase has two proton-conducting pathways through which protons are taken up to the interior part of the enzyme from the N-side of the membrane. The K-pathway transfers merely substrate protons, which are consumed in the process of water formation at the catalytic site. The D-pathway transfers both substrate protons and protons that are pumped to the P-side of the membrane. This thesis focuses on the role of two conserved amino acids in proton translocation by cytochrome c oxidase, glutamate 278 and tryptophan 164. Glu278 is located at the end of the D-pathway and is thought to constitute the branching point for substrate and pumped protons. In this work, it was shown that although Glu278 has an important role in the proton transfer mechanism, its presence is not an obligatory requirement. Alternative structural solutions in the area around Glu278, much like the ones present in some distantly related heme-copper oxidases, could in the absence of Glu278 support the formation of a long hydrogen-bonded water chain through which proton transfer from the D-pathway to the catalytic site is possible. The other studied amino acid, Trp164, is hydrogen bonded to the ∆-propionate of heme a3 of the catalytic site. Mutation of this amino acid showed that it may be involved in regulation of proton access to a proton acceptor, a pump site, from which the proton later is expelled to the P-side of the membrane. The ion pair that is formed by the ∆-propionate of heme a3 and arginine 473 is likely to form a gate-like structure, which regulates proton mobility to the P-side of the membrane. The same gate may also be part of an exit path through which water molecules produced at the catalytically active site are removed towards the external side of the membrane. Time-resolved optical and electrometrical experiments with the Trp164 to phenylalanine mutant revealed a so far undetected step in the proton pumping mechanism. During the A to PR transition of the catalytic cycle, a proton is transferred from Glu278 to the pump site, located somewhere in the vicinity of the ∆-propionate of heme a3. A mechanism for proton pumping by cytochrome c oxidase is proposed on the basis of the presented results and the mechanism is discussed in relation to some relevant experimental data. A common proton pumping mechanism for all members of the heme-copper oxidase family is moreover considered.

Characterization of the molecular components and function of BARE-1, Hin-Mu and Mu transposition machineries

Transposable elements, transposons, are discrete DNA segments that are able to move or copy themselves from one locus to another within or between their host genome(s) without a requirement for DNA homology. They are abundant residents in virtually all the genomes studied, for instance, the genomic portion of TEs is approximately 3% in Saccharomyces cerevisiae, 45% in humans, and apparently more than 70% in some plant genomes such as maize and barley. Transposons plays essential role in genome evolution, in lateral transfer of antibiotic resistance genes among bacteria and in life cycle of certain viruses such as HIV-1 and bacteriophage Mu. Despite the diversity of transposable elements they all use a fundamentally similar mechanism called transpositional DNA recombination (transposition) for the movement within and between the genomes of their host organisms. The DNA breakage and joining reactions that underlie their transposition are chemically similar in virtually all known transposition systems. The similarity of the reactions is also reflected in the structure and function of the catalyzing enzymes, transposases and integrases. The transposition reactions take place within the context of a transposition machinery, which can be particularly complex, as in the case of the VLP (virus like particle) machinery of retroelements, which in vivo contains RNA or cDNA and a number of element encoded structural and catalytic proteins. Yet, the minimal core machinery required for transposition comprises a multimer of transposase or integrase proteins and their binding sites at the element DNA ends only. Although the chemistry of DNA transposition is fairly well characterized, the components and function of the transposition machinery have been investigated in detail for only a small group of elements. This work focuses on the identification, characterization, and functional studies of the molecular components of the transposition machineries of BARE-1, Hin-Mu and Mu. For BARE-1 and Hin-Mu transpositional activity has not been shown previously, whereas bacteriophage Mu is a general model of transposition. For BARE-1, which is a retroelement of barley (Hordeum vulgare), the protein and DNA components of the functional VLP machinery were identified from cell extracts. In the case of Hin-Mu, which is a Mu-like prophage in Haemophilus influenzae Rd genome, the components of the core machinery (transposase and its binding sites) were characterized and their functionality was studied by using an in vitro methodology developed for Mu. The function of Mu core machinery was studied for its ability to use various DNA substrates: Hin-Mu end specific DNA substrates and Mu end specific hairpin substrates. The hairpin processing reaction by MuA was characterized in detail. New information was gained of all three machineries. The components or their activity required for functional BARE-1 VLP machinery and retrotransposon life cycle were present in vivo and VLP-like structures could be detected. The Hin-Mu core machinery components were identified and shown to be functional. The components of the Mu and Hin-Mu core machineries were partially interchangeable, reflecting both evolutionary conservation and flexibility within the core machineries. The Mu core machinery displayed surprising flexibility in substrate usage, as it was able to utilize Hin-Mu end specific DNA substrates and to process Mu end DNA hairpin substrates. This flexibility may be evolutionarily and mechanistically important.

Missing-In-Metastasis (MIM) regulates cell morphology by promoting plasma membrane and actin cytoskeleton dynamics

The cells of multicellular organisms have differentiated to carry out specific functions that are often accompanied by distinct cell morphology. The actin cytoskeleton is one of the key regulators of cell shape subsequently controlling multiple cellular events including cell migration, cell division, endo- and exocytosis. A large set of actin regulating proteins has evolved to achieve and tightly coordinate this wide range of functions. Some actin regulator proteins have so-called house keeping roles and are essential for all eukaryotic cells, but some have evolved to meet the requirements of more specialized cell-types found in higher organisms enabling complex functions of differentiated organs, such as liver, kidney and brain. Often processes mediated by the actin cytoskeleton, like formation of cellular protrusions during cell migration, are intimately linked to plasma membrane remodeling. Thus, a close cooperation between these two cellular compartments is necessary, yet not much is known about the underlying molecular mechanisms. This study focused on a vertebrate-specific protein called missing-in-metastasis (MIM), which was originally characterized as a metastasis suppressor of bladder cancer. We demonstrated that MIM regulates the dynamics of actin cytoskeleton via its WH2 domain, and is expressed in a cell-type specific manner. Interestingly, further examination showed that the IM-domain of MIM displays a novel membrane tubulation activity, which induces formation of filopodia in cells. Following studies demonstrated that this membrane deformation activity is crucial for cell protrusions driven by MIM. In mammals, there are five members of IM-domain protein family. Functions and expression patterns of these family members have remained poorly characterized. To understand the physiological functions of MIM, we generated MIM knockout mice. MIM-deficient mice display no apparent developmental defects, but instead suffer from progressive renal disease and increased susceptibility to tumors. This indicates that MIM plays a role in the maintenance of specific physiological functions associated with distinct cell morphologies. Taken together, these studies implicate MIM both in the regulation of the actin cytoskeleton and the plasma membrane. Our results thus suggest that members of MIM/IRSp53 protein family coordinate the actin cytoskeleton:plasma membrane interface to control cell and tissue morphogenesis in multicellular organisms.

Short-term plant-decomposer feedbacks in grassland plants

Plant species differ in their effects on ecosystem productivity and it is recognised that these effects are partly due to plant species-specific influences on soil processes. Until recently, however, not much attention was given to the potential role played by soil biota in these species-specific effects. While soil decomposers are responsible for governing the availability of nutrients for plant production, they simultaneously depend on the amount of carbon provided by plants. Litter and rhizodeposition constitute the two basal resources that plants provide to soil decomposer food webs. While it has been shown that both of these can have effects on soil decomposer communities that differ among plant species, the putative significance of these effects for plant nitrogen (N) acquisition is currently understudied. My PhD work aimed at clarifying whether the species-specific influences of three temperate grassland plants on the soil microfood-web, through rhizodeposition and litter, can feed back to plant N uptake. The methods and approach used (15N labelling of plant litter in microcosm experiments) revealed to be an effective combination of tools in studying these feedbacks. Plant effects on soil organisms were shown to differ significantly between plant species and the effects could be followed across several trophic levels. The labelling of litter further permitted the evaluation of plant acquisition of N derived from soil organic matter. The results show that the structure of the soil microfood-web can have a significant role in plant N acquisition when the structure is experimentally manipulated, such as when comparing systems consisting of microbes to those consisting of microbes and their grazers. However, despite this, the results indicate that differences in N uptake from soil organic matter between different plant species are not related to the effects these species exert on the structure of the soil microfood-web. Rather, these differences in N uptake seem to be determined by other species-specific traits of live plants and their litter. My results thus indicate that different resources provided by different plant species may not induce species-specific decomposer feedbacks on plant N uptake from soil organic matter. This further suggests that the species-specific plant effects on soil decomposer communities may not, at least in the short term, have significant consequences on plant production.

Genomic approach to organ-specific growth control

Growth is a fundamental aspect of life cycle of all organisms. Body size varies highly in most animal groups, such as mammals. Moreover, growth of a multicellular organism is not uniform enlargement of size, but different body parts and organs grow to their characteristic sizes at different times. Currently very little is known about the molecular mechanisms governing this organ-specific growth. The genome sequencing projects have provided complete genomic DNA sequences of several species over the past decade. The amount of genomic sequence information, including sequence variants within species, is constantly increasing. Based on the universal genetic code, we can make sense of this sequence information as far as it codes proteins. However, less is known about the molecular mechanisms that control expression of genes, and about the variations in gene expression that underlie many pathological states in humans. This is caused in part by lack of information about the second genetic code that consists of the binding specificities of transcription factors and the combinatorial code by which transcription factor binding sites are assembled to form tissue-specific and/or ligand-regulated enhancer elements. This thesis presents a high-throughput assay for identification of transcription factor binding specificities, which were then used to measure the DNA binding profiles of transcription factors involved in growth control. We developed ‘enhancer element locator’, a computational tool, which can be used to predict functional enhancer elements. A genome-wide prediction of human and mouse enhancer elements generated a large database of enhancer elements. This database can be used to identify target genes of signaling pathways, and to predict activated transcription factors based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c- and N-Myc genes, which has implications to organ specific growth control and tumor type specificity of oncogenes. Furthermore, we were able to locate a variation in a single nucleotide, which carries a susceptibility to colorectal cancer, to an enhancer element and propose a mechanism by which this SNP might be involved in generation of colorectal cancer.

LFA-1 Integrin beta2 Chain Phosphorylation Regulates Protein Interactions and Mediates Signals in T Cells

Integrins are heterodimeric transmembrane adhesion receptors composed of alpha- and beta-subunits and they are vital for the function of multicellular organisms. Integrin-mediated adhesion is a complex process involving both affinity regulation and coupling to the actin cytoskeleton. Integrins also function as bidirectional signaling devices, regulating cell adhesion and migration after inside-out signaling, but also signal into the cell to regulate growth, differentiation and apoptosis after ligand binding. The LFA-1 integrin is exclusively expressed in leukocytes and is of fundamental importance for the function of the immune system. The LFA-1 integrins have short intracellular tails, which are devoid of catalytic activity. These cytoplasmic domains are important for integrin regulation and both the alpha and beta chain become phosphorylated. The alpha chain is constitutively phosphorylated, but the beta chain becomes phosphorylated on serine and functionally important threonine residues only after cell activation. The cytoplasmic tails of LFA-1 bind to many cytoskeletal and signaling proteins regulating numerous cell functions. However, the molecular mechanisms behind these interactions have been poorly understood. Phosphorylation of the cytoplasmic tails of the LFA-1 integrin could provide a mechanism to regulate integrin-mediated cytoskeletal interactions and take part in T cell signaling. In this study, the effects of phosphorylation of LFA-1 integrin cytoplasmic tails on different cellular functions were examined. Site-specific phosphorylation of both the alpha- and beta-chains of the LFA-1 was shown to have a role in the regulation of the LFA-1 integrin.Alpha-chain Ser1140 is needed for integrin conformational changes after chemokine- or integrin ligand-induced activation or after activation induced by active Rap1, whereas beta-chain binds to 14-3-3 proteins through the phosphorylated Thr758 and mediates cytoskeletal reorganization. Thr758 phosphorylation also acts as a molecular switch to inhibit filamin binding and allows 14-3-3 protein binding to integrin cytoplasmic domain, and it was also shown to lead to T cell adhesion, Rac-1/Cdc42 activation and expression of the T cell activation marker CD69, indicating a signaling function for Thr758 phosphorylation in T cells. Thus, phosphorylation of the cytoplasmic tails of LFA-1 plays an important role in different functions of the LFA-1 integrin in T cells. It is of vital importance to study the mechanisms and components of integrin regulation since leukocyte adhesion is involved in many functions of the immune system and defects in the regulation of LFA-1 contributes to auto-immune diseases and fundamental defects in the immune system.

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

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

Environmental effects of thermal and radioactive discharges from nuclear power plants in the boreal brackish-water conditions of the northern Baltic Sea

During recent decades, thermal and radioactive discharges from nuclear power plants into the aquatic environment have become the subject of lively debate as an ecological concern. The target of this thesis was to summarize the large quantity of results obtained in extensive monitoring programmes and studies carried out in recipient sea areas off the Finnish nuclear power plants at Loviisa and Olkiluoto during more than four decades. The Loviisa NPP is located on the coast of the Gulf of Finland and Olkiluoto NPP on that of the Bothnian Sea. The state of the Gulf of Finland is clearly more eutrophic; the nutrient concentrations in the surface water are about 1½ 2 times higher at Loviisa than at Olkiluoto, and the total phosphorus concentrations still increased in both areas (even doubled at Loviisa) between the early 1970s and 2000. Thus, it is a challenge to distinguish the local effects of thermal discharges from the general eutrophication process of the Gulf of Finland. The salinity is generally low in the brackish-water conditions of the northern Baltic Sea, being however about 1 higher at Olkiluoto than at Loviisa (the salinity of surface water varying at the latter from near to 0 in early spring to 4 6 in late autumn). Thus, many marine and fresh-water organisms live in the Loviisa area close to their limit of existence, which makes the biota sensitive to any additional stress. The characteristics of the discharge areas of the two sites differ from each other in many respects: the discharge area at Loviisa is a semi-enclosed bay in the inner archipelago, where the exchange of water is limited, while the discharge area at Olkiluoto is more open, and the exchange of water with the open Bothnian Sea is more effective. The effects of the cooling water discharged from the power plants on the temperatures in the sea were most obvious in winter. The formation of a permanent ice cover in the discharge areas has been delayed in early winter, and the break-up of the ice occurs earlier in spring. The prolonging of the growing season and the disturbance of the overwintering time, in conditions where the biota has adjusted to a distinct rest period in winter, have been the most significant biological effects of the thermal pollution. The soft-bottom macrofauna at Loviisa has deteriorated to the point of almost total extinction at many sampling stations during the past 40 years. A similar decline has been reported for the whole eastern Gulf of Finland. However, the local eutrophication process seems to have contributed into the decline of the zoobenthos in the discharge area at Loviisa. Thermal discharges have increased the production of organic matter, which again has led to more organic bottom deposits. These have in turn increased the tendency of the isolated deeps to a depletion of oxygen, and this has further caused strong remobilization of phosphorus from the bottom sediments. Phytoplankton primary production and primary production capacity doubled in the whole area between the late 1960s and the late 1990s, but started to decrease a little at the beginning of this century. The focus of the production shifted from spring to mid- and late summer. The general rise in the level of primary production was mainly due to the increase in nutrient concentrations over the whole Gulf of Finland, but the thermal discharge contributed to a stronger increase of production in the discharge area compared to that in the intake area. The eutrophication of littoral vegetation in the discharge area has been the most obvious, unambiguous and significant biological effect of the heated water. Myriophyllum spicatum, Potamogeton perfoliatus and Potamogeton pectinatus, and vigorous growths of numerous filamentous algae as their epiphytes have strongly increased in the vicinity of the cooling water outlet, where they have formed dense populations in the littoral zone in late summer. However, the strongest increase of phytobenthos has extended only to a distance of about 1 km from the outlet, i.e., the changes in vegetation have been largest in those areas that remain ice-free in winter. Similar trends were also discernible at Olkiluoto, but to a clearly smaller extent, which was due to the definitely weaker level of background eutrophy and nutrient concentrations in the Bothnian Sea, and the differing local hydrographical and biological factors prevailing in the Olkiluoto area. The level of primary production has also increased at Olkiluoto, but has remained at a clearly lower level than at Loviisa. In spite of the analogous changes observed in the macrozoobenthos, the benthic fauna has remained strong and diversified in the Olkiluoto area. Small amounts of local discharge nuclides were regularly detected in environmental samples taken from the discharge areas: tritium in seawater samples, and activation products, such as 60Co, 58Co, 54Mn, 110mAg, 51Cr, in suspended particulate matter, bottom sediments and in several indicator organisms (e.g., periphyton and Fucus vesiculosus) that effectively accumulate radioactive substances from the medium. The tritium discharges and the consequent detection frequency and concentrations of tritium in seawater were higher at Loviisa, but the concentrations of the activation products were higher at Olkiluoto, where traces of local discharge nuclides were also observed over a clearly wider area, due to the better exchange of water than at Loviisa, where local discharge nuclides were only detected outside Hästholmsfjärden Bay quite rarely and in smaller amounts. At the farthest, an insignificant trace amount (0.2 Bq kg-1 d.w.) of 60Co originating from Olkiluoto was detected in Fucus at a distance of 137 km from the power plant. Discharge nuclides from the local nuclear power plants were almost exclusively detected at the lower trophic levels of the ecosystems. Traces of local discharge nuclides were very seldom detected in fish, and even then only in very low quantities. As a consequence of the reduced discharges, the concentrations of local discharge nuclides in the environment have decreased noticeably in recent years at both Loviisa and Olkiluoto. Although the concentrations in environmental samples, and above all, the discharge data, are presented as seemingly large numbers, the radiation doses caused by them to the population and to the biota are very low, practically insignificant. The effects of the thermal discharges have been more significant, at least to the wildlife in the discharge areas of the cooling water, although the area of impact has been relatively small. The results show that the nutrient level and the exchange of water in the discharge area of a nuclear power plant are of crucial importance.

A systematic-ecological approach to Baltic Sea ice studies of algae and protists

The seasonal occurrence of sea ice that annually covers almost half the Baltic Sea area provides a unique habitat for halo- and cold temperature-tolerant extremophiles. Baltic Sea ice biology has more than 100 years of tradition that began with the floristic observation of species by the early pioneers using light microscopic techniques that were the only thing available at the time. Since the discovery of life within sea ice, more technologies have become available for taxonomy. Electron microscopy and genetic evidence have been used to identify sea ice biota revealing increased numbers of taxa. Meanwhile ecologists have used light microscopic cell enumeration in addition to the chemical and physical properties of sea ice in attempts to explain the food web structure of sea ice and its functions. Thus, during the Baltic winter, the sea ice hosts more abundant and diverse microbial communities than the water column beneath it. These communities are typically dominated by autotrophic diatoms together with a diverse assortment of dinoflagellates, auto- and heterotrophic flagellates, ciliates, metazoan rotifers and bacteria, which are mostly responsible for the recycling of nutrients. This thesis comprises ecological and systematic studies. In addition to the results of the previous studies carried out on landfast ice, the data presented here provide new insight into the spatial distribution of pelagial sea ice, which has remained largely unexplored. The studies reveal spatial heterogeneity in the pelagial sea ice of the Gulf of Bothnia. There were mismatches in chlorophyll-a concentrations and in photosynthetic efficiencies of the communities studied. The temporal succession was followed and experimental studies performed investigating the community responses towards increased or decreased light in landfast ice in the Gulf of Finland. The systematic studies carried out with established dinoflagellate cultures revealed a new resting cyst belonging to common sea ice dinoflagellate, Scrippsiella hangoei (Schiller) Larsen 1995. The cyst can be used to explain the overwintering of this species during prolonged periods of darkness. The dissimilarities and similarities in the material isolated from the sea ice called for description of a new subspecies Heterocapsa arctica ssp. frigida. The cells obtained in the cultured material were unlike those of the previously described species, necessitating description of ssp. frigida. As a result of its own unique habitus, the subspecies had been noted by Finnish taxonomists during the past three decades and thus its annual occurrence and geographical distribution in the Baltic Sea. This illustrates how combining ecology and systematics increases our understanding of organisms.

Causes and consequences of variation in nestling immune function

Defence against pathogens is a vital need of all living organisms that has led to the evolution of complex immune mechanisms. However, although immunocompetence the ability to resist pathogens and control infection has in recent decades become a focus for research in evolutionary ecology, the variation in immune function observed in natural populations is relatively little understood. This thesis examines sources of this variation (environmental, genetic and maternal effects) during the nestling stage and its fitness consequences in wild populations of passerines: the blue tit (Cyanistes caeruleus) and the collared flycatcher (Ficedula albicollis). A developing organism may face a dilemma as to whether to allocate limited resources to growth or to immune defences. The optimal level of investment in immunity is shaped inherently by specific requirements of the environment. If the probability of contracting infection is low, maintaining high growth rates even at the expense of immune function may be advantageous for nestlings, as body mass is usually a good predictor of post-fledging survival. In experiments with blue tits and haematophagous hen fleas (Ceratophyllus gallinae) using two methods, methionine supplementation (to manipulate nestlings resource allocation to cellular immune function) and food supplementation (to increase resource availability), I confirmed that there is a trade-off between growth and immunity and that the abundance of ectoparasites is an environmental factor affecting allocation of resources to immune function. A cross-fostering experiment also revealed that environmental heterogeneity in terms of abundance of ectoparasites may contribute to maintaining additive genetic variation in immunity and other traits. Animal model analysis of extensive data collected from the population of collared flycatchers on Gotland (Sweden) allowed examination of the narrow-sense heritability of PHA-response the most commonly used index of cellular immunocompetence in avian studies. PHA-response is not heritable in this population, but is subject to a non-heritable origin (presumably maternal) effect. However, experimental manipulation of yolk androgen levels indicates that the mechanism of the maternal effect in PHA-response is not in ovo deposition of androgens. The relationship between PHA-response and recruitment was studied for over 1300 collared flycatcher nestlings. Multivariate selection analysis shows that it is body mass, not PHA-response, that is under direct selection. PHA-response appears to be related to recruitment because of its positive relationship with body mass. These results imply that either PHA-response fails to capture the immune mechanisms that are relevant for defence against pathogens encountered by fledglings or that the selection pressure from parasites is not as strong as commonly assumed.