Mining Sequential Data – in Search of Segmental Structures

Segmentation is a data mining technique yielding simplified representations of sequences of ordered points. A sequence is divided into some number of homogeneous blocks, and all points within a segment are described by a single value. The focus in this thesis is on piecewise-constant segments, where the most likely description for each segment and the most likely segmentation into some number of blocks can be computed efficiently. Representing sequences as segmentations is useful in, e.g., storage and indexing tasks in sequence databases, and segmentation can be used as a tool in learning about the structure of a given sequence. The discussion in this thesis begins with basic questions related to segmentation analysis, such as choosing the number of segments, and evaluating the obtained segmentations. Standard model selection techniques are shown to perform well for the sequence segmentation task. Segmentation evaluation is proposed with respect to a known segmentation structure. Applying segmentation on certain features of a sequence is shown to yield segmentations that are significantly close to the known underlying structure. Two extensions to the basic segmentation framework are introduced: unimodal segmentation and basis segmentation. The former is concerned with segmentations where the segment descriptions first increase and then decrease, and the latter with the interplay between different dimensions and segments in the sequence. These problems are formally defined and algorithms for solving them are provided and analyzed. Practical applications for segmentation techniques include time series and data stream analysis, text analysis, and biological sequence analysis. In this thesis segmentation applications are demonstrated in analyzing genomic sequences.

Preliminary Phonological Analysis of Denjongka of Sikkim

This thesis is a preliminary phonological description of the Tibetan-related Denjongka language of Sikkim, India. Because the language has not been much researched and the previous studies have focused on other issues than phonology, the present paper is the first of its kind. The data for this thesis was gathered in Gangtok, the capital of Sikkim, from March to May 2004. I had four language informants from four different locations in Sikkim who spoke different dialects of Denjongka. One of the informants, from whom I recorded c. 900 words and 530 sentences, was used as the main data source for the analysis. First, I will give some ethnographic background information on the people who speak Denjongka. Next, I will discuss first the segmental and then the suprasegmental phonology of the language, which were analysed much in line with American structuralism. I also used acoustic analysis enabled by the Praat-program. Eight vowel phonemes were found. The phonemic status of /E/, however, is still suspect. I present some preliminary evidence for roundedness, frontness and height assimilation among the vowels. In the interpretation adopted in this analysis, there are no diphthongs in Denjongka. Forty consonant phonemes were found: 17 plosives, 7 affricates, 5 fricatives, 5 nasals, 4 liquids and 2 approximants. Denjongka plosives and affricates have four-way aspiration/voicing distinction: voiceless aspirated, voiceless unaspirated, voiceless slightly aspirated (devoiced), and voiced unaspirated. Two voiceless nasals and two voiceless liquids were found. Two phonation types were found to be contrastive, lax/breathy and tense/creaky. Nasalisation and length in vowels are phonemic. Denjongka is an incipient tone language. Tonal phenomena, which involve mainly pitch and phonation type, are complex. Pitch is most of the time predictable from the initial consonant and the phonation type. In some cases, however, pitch is the only contrastive feature between words. The description of Denjongka in this paper differs from the traditional four-tone system, which has been used in many descriptions of Tibetan-related languages. In the four-tone system, pitch is contrastive both in the high and low register, whereas in the present analysis pitch has been established to contrast only in the high register. Lastly, the appendices include a comparative word list of the four Denjongka dialects studied in this thesis.

Turun linna kerrottuna ja kertovana tilana

The Turku castle, founded c. 1300, has changed over the centuries from a medieval defensive structure into a Renaissance palace and from a derelict jailhouse in the late 19th century into a prime example of the Medieval built heritage in Finland. Today, it is first and foremost a monument to the Medieval and Renaissance heyday of the castle. This is apparent in the architectural forms that have been carefully restored and reconstructed. It also becomes clear in all kinds of narratives, both visual (like the set of miniatures about the different stages of the construction of the castle) and textual (as during the guided tours). For the first time in the architectural history of the Turku castle, the Medieval, the Renaissance, the Modern, and the Present as architecturally constructed or reconstructed spaces can all be visited within the same hour. As a result, the monumental Turku castle may even be deemed anachronistic or inauthentic. In this study I look at the ways in which the Turku castle is, indeed, anachronistic and inauthentic. My main objective, however, is to find ways in which the anachronisms and inauthenticities are overcome in a positive way. I base my analysis of the Turku castle on three theoretical standpoints. First, I am studying the castle as space, described by Michel de Certeau as a practiced place (de Certeau 2002). Second, I am approaching the numerous narrative aspects of the castle following Paul Ricoeur s analysis of narrative as a threefold mimetic process (Ricoeur 1990). From these two theoretical settings I have summoned the concept of narrative space. The life and work at the castle are based on expectations and understandings of the historical surroundings. My third theoretical choice is to study this applied knowledge of the place as the management of blocks of knowledge in communication (Robert de Beaugrande 1980). Combining the theoretical starting points of space and narrative , I am approaching the castle as if it were an evolving set of narratives, narrated in space but also through space. Seeing e.g. the restoration teams of the mid-20th century and the present day tour guides as creative narrators, I am looking beyond the dilemma of the anachronistic spaces. What transpires is an inter-connected web of texts and spaces, tangible and intangible narratives. My analytical key to these narrative relationships is the threefold mimetic process of pre-figuration, con-figuration, and re-figuration, inspired by the writings of Paul Ricoeur (1990). This way, the past can be seen as a pool of endless possibilities to emplot place, time, and action into a narrative space. The narratives convey images of the past that may be contested by other images, and the power to narrate in the first place can be challenged and re-distributed.

Digitoitujen väri-infrailmakuvien numeerinen monikuvatulkinta metsäninventoinnissa

The use of remote sensing imagery as auxiliary data in forest inventory is based on the correlation between features extracted from the images and the ground truth. The bidirectional reflectance and radial displacement cause variation in image features located in different segments of the image but forest characteristics remaining the same. The variation has so far been diminished by different radiometric corrections. In this study the use of sun azimuth based converted image co-ordinates was examined to supplement auxiliary data extracted from digitised aerial photographs. The method was considered as an alternative for radiometric corrections. Additionally, the usefulness of multi-image interpretation of digitised aerial photographs in regression estimation of forest characteristics was studied. The state owned study area located in Leivonmäki, Central Finland and the study material consisted of five digitised and ortho-rectified colour-infrared (CIR) aerial photographs and field measurements of 388 plots, out of which 194 were relascope (Bitterlich) plots and 194 were concentric circular plots. Both the image data and the field measurements were from the year 1999. When examining the effect of the location of the image point on pixel values and texture features of Finnish forest plots in digitised CIR photographs the clearest differences were found between front-and back-lighted image halves. Inside the image half the differences between different blocks were clearly bigger on the front-lighted half than on the back-lighted half. The strength of the phenomenon varied by forest category. The differences between pixel values extracted from different image blocks were greatest in developed and mature stands and smallest in young stands. The differences between texture features were greatest in developing stands and smallest in young and mature stands. The logarithm of timber volume per hectare and the angular transformation of the proportion of broadleaved trees of the total volume were used as dependent variables in regression models. Five different converted image co-ordinates based trend surfaces were used in models in order to diminish the effect of the bidirectional reflectance. The reference model of total volume, in which the location of the image point had been ignored, resulted in RMSE of 1,268 calculated from test material. The best of the trend surfaces was the complete third order surface, which resulted in RMSE of 1,107. The reference model of the proportion of broadleaved trees resulted in RMSE of 0,4292 and the second order trend surface was the best, resulting in RMSE of 0,4270. The trend surface method is applicable, but it has to be applied by forest category and by variable. The usefulness of multi-image interpretation of digitised aerial photographs was studied by building comparable regression models using either the front-lighted image features, back-lighted image features or both. The two-image model turned out to be slightly better than the one-image models in total volume estimation. The best one-image model resulted in RMSE of 1,098 and the two-image model resulted in RMSE of 1,090. The homologous features did not improve the models of the proportion of broadleaved trees. The overall result gives motivation for further research of multi-image interpretation. The focus may be improving regression estimation and feature selection or examination of stratification used in two-phase sampling inventory techniques. Keywords: forest inventory, digitised aerial photograph, bidirectional reflectance, converted image co­ordinates, regression estimation, multi-image interpretation, pixel value, texture, trend surface

Monitoring of Fungal Growth and Degradation of Wood

Valko- ja ruskolahosienet tunnetaan luonnossa tehokkaimpina puun ja karikkeen lignoselluloosan lahottajina. Valkolahosienet pystyvät hajottamaan kaikkia puun osia: ligniiniä, selluloosaa ja hemiselluloosaa. Selektiivisesti ligniiniä hajottavat sienet lahottavat puusta suhteessa enemmän vaikeasti hajoavaa ligniiniä kuin selluloosaa tai hemiselluloosaa, jolloin jäljelle jää valkoista ja miltei puhdasta selluloosaa. Bioteknisissä sovelluksissa juuri selektiviiviset valkolahottajat ovat kiinnostavia. Niiden avulla voidaan puuhaketta esikäsitellä esimerkiksi paperinvalmistuksessa haitallisen ligniinin poistamiseksi. Ruskolahosienet ovat huomattavia puun, puutavaran ja puisten rakenteiden lahottajia, kuten tässä työssä käytetty Gloeophyllum trabeum (saunasieni ) ja Poria (Postia) placenta (istukkakääpä). Ruskolahosienet hajottavat puusta hemiselluloosan lisäksi selluloosaa, jolloin jää jäljelle ruskea ja jauhomaiseksi mureneva ligniini. Ruskolahosienet muovaavat ligniiniä jonkin verran. Kahden ruskolahosienen G. trabeumin ja P. placentan lisäksi tutkittiin valkolahosieniä, joista Ceriporiopsis subvermispora (karstakääpä) ja harvinainen Physisporinus rivulosus -sieni (talikääpä) hajottavat ligniiniä erittäin selektiivisesti. Phanerochaete chrysosporium on kaikkialla paljon tutkittu sieni, ja Phlebia radiata valkolahosientä (rusorypykkä) on tutkittu paljon mikrobiologian osastolla. Lisäksi tutkittiin Phlebia tremellosa -sienten (hytyrypykkä) ligninolyyttisten entsyymien tuottoa ja 14C-leimatun synteettisen ligniinin (DHP) hajotusta. P. radiata ja P. tremellosa -sienten on todettu aiemmin hajottavan ligniiniä selektiivisesti. Työssä selvitettiin miten sienten kasvua voi mitata, miten vertailukelpoisia eri mittaamismenetelmillä saadut tulokset ovat ja ilmenevätkö sienten aktiivisimmat kasvuvaiheet samaan aikaan eri menetelmillä mitattuna. Tärkeimmät tulokset olivat seuraavat havainnot: (i) P. radiata ja P. tremellosa -sienikannat tuottivat ligniini- ja mangaaniperoksidaasientsyymejä (LiP ja MnP) sekä lakkaasia, ja sienistä puhdistettiin 2-3 LiP- ja P. radiatasta yksi MnP-entsyymi; (ii) P. tremellosa -sienet hajottivat leimattua synteettistä ligniiniä (DHP) yhtä hyvin kuin paljon tutkitut P. chrysosporium ja P. radiata -sienet; (iii) puu, sienen luonnollinen kasvualusta, lisäsi valkolaho- ja ruskolahosienten demetoksylaatiota [O14CH3]-leimatusta ligniinin malliyhdisteestä 14CO2:ksi ilman puuta olleeseen alustaan verrattuna; (iv) demetoksylaatio (14CO2:n tuotto) oli normaalissa ilma-atmosfäärissä useimmiten parempi happeen verrattuna; (v) hapessa paras 14CO2:n tuotto saatiin puupalakasvatuksissa, joihin oli lisätty ravinnetyppeä tai typen lisäksi glukoosia sekä valkolaho- että ruskolahosienillä; (vi) ilmassa 14CO2:n tuotto oli puulla voimakkainta valkolahosienillä ilman lisäravinteita, kun taas G. trabeum -sienellä se oli yhtä hyvä eri alustoissa; (vii) biomassan muodostuminen rihmastojen ergosterolipitoisuuksista mitattuna oli ruskolahosienillä parempi kuin valkolahosienillä; (viii) ja biomassojen huippupitoisuudet olivat 6:lla sienellä eri suuruisia ja niiden maksimimäärien ajankohdat vaihtelivat viiden viikon kasvatusten kuluessa. Mikrobiologian osastolla Viikissä eristetty ja paljon tutkittu P. radiata -valkolahosieni oli mukana kaikissa tehdyissä kokeissa. Sienen LiP-aktiivisuus ja 14CO2:n tuotto 14C-rengas-leimatusta synteettisestä ligniinistä (DHP) korreloivat erittäin hyvin. Biomassan muodostuminen ergosterolilla määritettynä tuki hyvin entsyymiaktiivisuusmittauksilla ja isotooppikasvatuksilla saatuja tuloksia.

Nutrition and bone metabolism : In vivo effects of inorganic phosphate on rats and in vitro effects of bioactive tripeptides on human osteoblasts

Nutrition affects bone health throughout life. To optimize peak bone mass development and maintenance, it is important to pay attention to the dietary factors that enhance and impair bone metabolism. In this study, the in vivo effects of inorganic dietary phosphate and the in vitro effects of bioactive tripeptides, IPP, VPP and LKP were investigated. Dietary phosphate intake is increased through the use of convenience foods and soft drinks rich in phosphate-containing food additives. Our results show that increased dietary phosphate intake hinders mineral deposition in cortical bone and diminishes bone mineral density (BMD) in the aged skeleton in a rodent model (Study I). In the growing skeleton (Study II), increased phosphate intake was observed to reduce bone material and structural properties, leading to diminished bone strength. Studies I and II revealed that a low Ca:P ratio has negative effects on the mature and growing rat skeleton even when calcium intake is sufficient. High dietary protein intake is beneficial for bone health. Protein is essential for bone turnover and matrix formation. In addition, hydrolysis of proteins in the gastrointestinal tract produces short peptides that possess a biological function beyond that of being tissue building blocks. The effects of three bioactive tripeptides, IPP, VPP and LKP, were assessed in short- and long-term in vitro experiments. Short-term treatment (24 h) with tripeptide IPP, VPP or LKP influenced osteoblast gene expression (Study III). IPP in particular, regulates genes associated with cell differentiation, cell growth and cell signal transduction. The upregulation of these genes indicates that IPP enhances osteoblast proliferation and differentiation. Long-term treatment with IPP enhanced osteoblast gene expression in favour of bone formation and increased mineralization (Study IV). The in vivo effects of IPP on osteoblast differentiation might differ since eating frequency drives food consumption, and protein degradation products, such as bioactive peptides, are available periodically, not continuously as in this study. To sum up, Studies I and II raise concern about the appropriate amount of dietary phosphate to support bone health as excess is harmful. Studies III and IV in turn, support findings of the beneficial effects of dietary protein on bone and provide a mechanistic explanation since cell proliferation and osteoblast function were improved by treatment with bioactive tripeptide IPP.

Organized Nanostructures of Thermoresponsive Poly(N-isopropylacrylamide) Block Copolymers Obtained Through Controlled RAFT Polymerization

Kontrolloidut radikaalipolymerointimenetelmät, kuten RAFT-polymerointi, ovat moderni tapa valmistaa polymeerejä säädellysti. RAFT-polymeroinnilla polymeerien ketjunpituutta, moolimassajakaumaa, mikrorakennetta (taktisuus, järjestys), koostumusta ja funktionaalisuutta kyetään hallitsemaan. Siten menetelmällä voidaan valmistaa uudenlaisia polymeeriarkkitektuureja, kuten blokki- ja tähtipolymeerejä, sekä hybridimateriaaleja ja biokonjugaatteja. Polymeeristen rakennuspalikoiden itsejärjestyminen, missä huolellisesti syntetisoidut polymeerit järjestyvät halutulla tavalla nanoskaalassa, on suosittu tutkimuskohde materiaalitieteessä. On huomattava, että blokkipolymeerien itsejärjestyminen on vielä suhteellisen nuori tutkimusaihe. Tämän hetkiset polymeeriset nanomateriaalit ovat suhteellisen yksinkertaisia luonnon luomuksiin verrattuina, tarjoten jatkuvasti uusia mahdollisuuksia seuraavan sukupolven polymeereille. Tässä työssä RAFT-polymeroinnilla syntetisoitiin amfifiilisiä di- ja triblokkikopolymeerejä sekä tutkittiin niiden järjestymistä nanorakenteiksi. Kaikissa blokkikopolymeereissä käytettiin lämpöherkkää poly(N-isopropyyliakryyliamidia). Siten polymeerit ja tutkitut materiaalit reagoivat lämpötilanmuutokseen ympäristössä eli ovat ns. ympäristöherkkiä. Työssä tutkittiin taktisuuden kontrollointia N-isopropyyliakryyliamidin RAFT-polymeroinnissa. Polymeerin taktisuutta sekä ketjunpituutta ja blokkijärjestystä säätämällä voitiin hallita polymeerin itsejärjestymistä vesiliuoksessa. Amfifiiliset polymeerit järjestyivät laimeissa vesiliuoksissa erilaisiksi misellirakenteiksi, muodostaen ns. mikrosäiliöitä. Tällaisilla polymeereillä odotetaan olevan sovelluksia esim. lääkeainevapautuksessa. Amfifiilejä käytetään myös esimerkiksi apuaineina pinnoitteissa ja kosmetiikassa. Kiinteässä tilassa tutkitut triblokkikopolymeerit muodostivat teoreettisesti ennustettuja morfologioita. Lämpöherkän materiaalin hydrogeelit toimivat suodatinmembraanina nanokokoluokassa. RAFT-polymeroinnilla syntetisoituja polymeereja voidaan sellaisenaan käyttää kultananopartikkeleiden päällystämiseen. Kultananopartikkelit ovat erittäin kiinostavia mm. niiden stabiilisuuden ja ainutlaatuisten pintaominaisuuksien vuoksi. Kun amfifiilisiä polymeerejä kiinnitettiin kultapartikkelin pinnalle, sen liuos- ja optisia ominaisuuksia voitiin säädellä pH:n ja lämpötilan avulla. Tällaisilla kultananopartikkeleilla on sovelluksia mm. diagnostiikassa, sensoreina ja solukuvauksessa.

An approach to palaeoseismicity in the Olkiluoto (sea) area during the early Holocene

Olkiluoto Island is situated in the northern Baltic Sea, near the southwestern coast of Finland, and is the proposed location of a spent nuclear fuel repository. This study examined Holocene palaeoseismicity in the Olkiluoto area and in the surrounding sea areas by computer simulations together with acoustic-seismic, sedimentological and dating methods. The most abundant rock type on the island is migmatic mica gneiss, intruded by tonalites, granodiorites and granites. The surrounding Baltic Sea seabed consists of Palaeoproterozoic crystalline bedrock, which is to a great extent covered by younger Mesoproterozoic sedimentary rocks. The area contains several ancient deep-seated fracture zones that divide it into bedrock blocks. The response of bedrock at the Olkiluoto site was modelled considering four future ice-age scenarios. Each scenario produced shear displacements of fractures with different times of occurrence and varying recovery rates. Generally, the larger the maximum ice load, the larger were the permanent shear displacements. For a basic case, the maximum shear displacements were a few centimetres at the proposed nuclear waste repository level, at proximately 500 m b.s.l. High-resolution, low-frequency echo-sounding was used to examine the Holocene submarine sedimentary structures and possible direct and indirect indicators of palaeoseismic activity in the northern Baltic Sea. Echo-sounding profiles of Holocene submarine sediments revealed slides and slumps, normal faults, debris flows and turbidite-type structures. The profiles also showed pockmarks and other structures related to gas or groundwater seepages, which might be related to fracture zone activation. Evidence of postglacial reactivation in the study area was derived from the spatial occurrence of some of the structures, especial the faults and the seepages, in the vicinity of some old bedrock fracture zones. Palaeoseismic event(s) (a single or several events) in the Olkiluoto area were dated and the palaeoenvironment was characterized using palaeomagnetic, biostratigraphical and lithostratigraphical methods, enhancing the reliability of the chronology. Combined lithostratigraphy, biostratigraphy and palaeomagnetic stratigraphy revealed an age estimation of 10 650 to 10 200 cal. years BP for the palaeoseismic event(s). All Holocene sediment faults in the northern Baltic Sea occur at the same stratigraphical level, the age of which is estimated at 10 700 cal. years BP (9500 radiocarbon years BP). Their movement is suggested to have been triggered by palaeoseismic event(s) when the Late Weichselian ice sheet was retreating from the site and bedrock stresses were released along the bedrock fracture zones. Since no younger or repeated traces of seismic events were found, it corroborates the suggestion that the major seismic activity occurred within a short time during and after the last deglaciation. The origin of the gas/groundwater seepages remains unclear. Their reflections in the echo-sounding profiles imply that part of the gas is derived from the organic-bearing Litorina and modern gyttja clays. However, at least some of the gas is derived from the bedrock. Additional information could be gained by pore water analysis from the pockmarks. Information on postglacial fault activation and possible gas and/or fluid discharges under high hydraulic heads has relevance in evaluating the safety assessment of a planned spent nuclear fuel repository in the region.

Conformal Field Theory Methods for Variants of Schramm-Loewner Evolutions

This thesis consists of an introduction, four research articles and an appendix. The thesis studies relations between two different approaches to continuum limit of models of two dimensional statistical mechanics at criticality. The approach of conformal field theory (CFT) could be thought of as the algebraic classification of some basic objects in these models. It has been succesfully used by physicists since 1980's. The other approach, Schramm-Loewner evolutions (SLEs), is a recently introduced set of mathematical methods to study random curves or interfaces occurring in the continuum limit of the models. The first and second included articles argue on basis of statistical mechanics what would be a plausible relation between SLEs and conformal field theory. The first article studies multiple SLEs, several random curves simultaneously in a domain. The proposed definition is compatible with a natural commutation requirement suggested by Dubédat. The curves of multiple SLE may form different topological configurations, ``pure geometries''. We conjecture a relation between the topological configurations and CFT concepts of conformal blocks and operator product expansions. Example applications of multiple SLEs include crossing probabilities for percolation and Ising model. The second article studies SLE variants that represent models with boundary conditions implemented by primary fields. The most well known of these, SLE(kappa, rho), is shown to be simple in terms of the Coulomb gas formalism of CFT. In the third article the space of local martingales for variants of SLE is shown to carry a representation of Virasoro algebra. Finding this structure is guided by the relation of SLEs and CFTs in general, but the result is established in a straightforward fashion. This article, too, emphasizes multiple SLEs and proposes a possible way of treating pure geometries in terms of Coulomb gas. The fourth article states results of applications of the Virasoro structure to the open questions of SLE reversibility and duality. Proofs of the stated results are provided in the appendix. The objective is an indirect computation of certain polynomial expected values. Provided that these expected values exist, in generic cases they are shown to possess the desired properties, thus giving support for both reversibility and duality.

Paleomagnetic and rock magnetic study with emphasis on the Precambrian intrusions and impact structures in Fennoscandia and South Africa

The importance of supercontinents in our understanding of the geological evolution of the planet Earth has been recently emphasized. The role of paleomagnetism in reconstructing lithospheric blocks in their ancient paleopositions is vital. Paleomagnetism is the only quantitative tool for providing ancient latitudes and azimuthal orientations of continents. It also yields information of content of the geomagnetic field in the past. In order to obtain a continuous record on the positions of continents, dated intrusive rocks are required in temporal progression. This is not always possible due to pulse-like occurrences of dykes. In this work we demonstrate that studies of meteorite impact-related rocks may fill some gaps in the paleomagnetic record. This dissertation is based on paleomagnetic and rock magnetic data obtained from samples of the Jänisjärvi impact structure (Russian Karelia, most recent 40Ar-39Ar age of 682 Ma), the Salla diabase dyke (North Finland, U-Pb 1122 Ma), the Valaam monzodioritic sill (Russian Karelia, U-Pb 1458 Ma), and the Vredefort impact structure (South Africa, 2023 Ma). The paleomagnetic study of Jänisjärvi samples was made in order to obtain a pole for Baltica, which lacks paleomagnetic data from 750 to ca. 600 Ma. The position of Baltica at ca. 700 Ma is relevant in order to verify whether the supercontinent Rodinia was already fragmented. The paleomagnetic study of the Salla dyke was conducted to examine the position of Baltica at the onset of supercontinent Rodinia's formation. The virtual geomagnetic pole (VGP) from Salla dyke provides hints that the Mesoproterozoic Baltica - Laurentia unity in the Hudsonland (Columbia, Nuna) supercontinent assembly may have lasted until 1.12 Ga. Moreover, the new VGP of Salla dyke provides new constraint on the timing of the rotation of Baltica relative to Laurentia (e.g. Gower et al., 1990). A paleomagnetic study of the Valaam sill was carried out in order to shed light into the question of existence of Baltica-Laurentia unity in the supercontinent Hudsonland. Combined with results from dyke complex of the Lake Ladoga region (Schehrbakova et al., 2008) a new robust paleomagnetic pole for Baltica is obtained. This pole places Baltica on a latitude of 10°. This low latitude location is supported also by Mesoproterozoic 1.5 1.3 Ga red-bed sedimentation (for example the Satakunta sandstone). The Vredefort impactite samples provide a well dated (2.02 Ga) pole for the Kaapvaal Craton. Rock magnetic data reveal unusually high Koenigsberger ratios (Q values) in all studied lithologies of the Vredefort dome. The high Q values are now first time also seen in samples from the Johannesburg Dome (ca. 120 km away) where there is no impact evidence. Thus, a direct causative link of high Q values to the Vredefort impact event can be ruled out.

Computational Methods for Reconstruction and Analysis of Genome-Scale Metabolic Networks

Metabolism is the cellular subsystem responsible for generation of energy from nutrients and production of building blocks for larger macromolecules. Computational and statistical modeling of metabolism is vital to many disciplines including bioengineering, the study of diseases, drug target identification, and understanding the evolution of metabolism. In this thesis, we propose efficient computational methods for metabolic modeling. The techniques presented are targeted particularly at the analysis of large metabolic models encompassing the whole metabolism of one or several organisms. We concentrate on three major themes of metabolic modeling: metabolic pathway analysis, metabolic reconstruction and the study of evolution of metabolism. In the first part of this thesis, we study metabolic pathway analysis. We propose a novel modeling framework called gapless modeling to study biochemically viable metabolic networks and pathways. In addition, we investigate the utilization of atom-level information on metabolism to improve the quality of pathway analyses. We describe efficient algorithms for discovering both gapless and atom-level metabolic pathways, and conduct experiments with large-scale metabolic networks. The presented gapless approach offers a compromise in terms of complexity and feasibility between the previous graph-theoretic and stoichiometric approaches to metabolic modeling. Gapless pathway analysis shows that microbial metabolic networks are not as robust to random damage as suggested by previous studies. Furthermore the amino acid biosynthesis pathways of the fungal species Trichoderma reesei discovered from atom-level data are shown to closely correspond to those of Saccharomyces cerevisiae. In the second part, we propose computational methods for metabolic reconstruction in the gapless modeling framework. We study the task of reconstructing a metabolic network that does not suffer from connectivity problems. Such problems often limit the usability of reconstructed models, and typically require a significant amount of manual postprocessing. We formulate gapless metabolic reconstruction as an optimization problem and propose an efficient divide-and-conquer strategy to solve it with real-world instances. We also describe computational techniques for solving problems stemming from ambiguities in metabolite naming. These techniques have been implemented in a web-based sofware ReMatch intended for reconstruction of models for 13C metabolic flux analysis. In the third part, we extend our scope from single to multiple metabolic networks and propose an algorithm for inferring gapless metabolic networks of ancestral species from phylogenetic data. Experimenting with 16 fungal species, we show that the method is able to generate results that are easily interpretable and that provide hypotheses about the evolution of metabolism.

Algorithms for 13C metabolic flux analysis

The metabolism of an organism consists of a network of biochemical reactions that transform small molecules, or metabolites, into others in order to produce energy and building blocks for essential macromolecules. The goal of metabolic flux analysis is to uncover the rates, or the fluxes, of those biochemical reactions. In a steady state, the sum of the fluxes that produce an internal metabolite is equal to the sum of the fluxes that consume the same molecule. Thus the steady state imposes linear balance constraints to the fluxes. In general, the balance constraints imposed by the steady state are not sufficient to uncover all the fluxes of a metabolic network. The fluxes through cycles and alternative pathways between the same source and target metabolites remain unknown. More information about the fluxes can be obtained from isotopic labelling experiments, where a cell population is fed with labelled nutrients, such as glucose that contains 13C atoms. Labels are then transferred by biochemical reactions to other metabolites. The relative abundances of different labelling patterns in internal metabolites depend on the fluxes of pathways producing them. Thus, the relative abundances of different labelling patterns contain information about the fluxes that cannot be uncovered from the balance constraints derived from the steady state. The field of research that estimates the fluxes utilizing the measured constraints to the relative abundances of different labelling patterns induced by 13C labelled nutrients is called 13C metabolic flux analysis. There exist two approaches of 13C metabolic flux analysis. In the optimization approach, a non-linear optimization task, where candidate fluxes are iteratively generated until they fit to the measured abundances of different labelling patterns, is constructed. In the direct approach, linear balance constraints given by the steady state are augmented with linear constraints derived from the abundances of different labelling patterns of metabolites. Thus, mathematically involved non-linear optimization methods that can get stuck to the local optima can be avoided. On the other hand, the direct approach may require more measurement data than the optimization approach to obtain the same flux information. Furthermore, the optimization framework can easily be applied regardless of the labelling measurement technology and with all network topologies. In this thesis we present a formal computational framework for direct 13C metabolic flux analysis. The aim of our study is to construct as many linear constraints to the fluxes from the 13C labelling measurements using only computational methods that avoid non-linear techniques and are independent from the type of measurement data, the labelling of external nutrients and the topology of the metabolic network. The presented framework is the first representative of the direct approach for 13C metabolic flux analysis that is free from restricting assumptions made about these parameters.In our framework, measurement data is first propagated from the measured metabolites to other metabolites. The propagation is facilitated by the flow analysis of metabolite fragments in the network. Then new linear constraints to the fluxes are derived from the propagated data by applying the techniques of linear algebra.Based on the results of the fragment flow analysis, we also present an experiment planning method that selects sets of metabolites whose relative abundances of different labelling patterns are most useful for 13C metabolic flux analysis. Furthermore, we give computational tools to process raw 13C labelling data produced by tandem mass spectrometry to a form suitable for 13C metabolic flux analysis.

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.

Structure and function of GDNF receptor alpha splice variants

The growth factors of the glial cell line-derived neurotrophic factor (GDNF) family consisting of GDNF, neurturin (NRTN), artemin (ARTN) and persephin (PSPN), are involved in the development, differentiation and maintenance of many types of neurons. They also have important functions outside the nervous system in the development of kidney, testis and thyroid gland. Each of these GFLs preferentially binds to one of the glycosylphosphatidylinositol (GPI)-anchored GDNF family receptors α (GFRα). GDNF binds to GFRα1, NRTN to GFRα2, ARTN to GFRα3 and PSPN to GFRα4. The GFLs in the complex with their cognate GFRα receptors all bind to and signal through the receptor tyrosine kinase RET. Alternative splicing of the mouse GFRα4 gene yields three splice isoforms. These had been described as putative GPI-anchored, transmembrane and soluble forms. My goal was to characterise the function of the different forms of mouse GFRα4. I firstly found that the putative GPI-anchored GFRα4 (GFRα4-GPI) is glycosylated, membrane-bound, GPI-anchored and interacts with PSPN and RET. We also showed that mouse GFRα4-GPI mediates PSPN-induced phosphorylation of RET, promotes PSPN-dependent neuronal differentiation of the rat pheochromocytoma cell line PC6-3 and PSPN-dependent survival of cerebellar granule neurons (CGN). However, although this receptor can mediate PSPN-signalling and activate RET, GFRα4-GPI does not recruit RET into lipid rafts. The recruitment of RET into lipid rafts has previously been thought to be a crucial event for GDNF- and GFL-mediated signalling via RET. I secondly demonstrated that the putative transmembrane GFRα4 (GFRα4-TM) is indeed a real transmembrane GFRα4 protein. Although it has a weak binding capacity for PSPN, it can not mediate PSPN-dependent phosphorylation of RET, neuronal differentiation or survival. These data show that GFRα4-TM is inactive as a receptor for PSPN. Surprisingly, GFRα4-TM can negatively regulate PSPN-mediated signalling via GFRα4-GPI. GFRα4-TM interacts with GFRα4-GPI and blocks PSPN-induced phosphorylation of RET, neuronal differentiation as well as survival. Taken together, our data show that GFRα4-TM may act as a dominant negative inhibitor of PSPN-mediated signaling. The most exciting part of my work was the finding that the putative soluble GFRα4 (GFRα4-sol) can form homodimers and function as an agonist of the RET receptor. In the absence of PSPN, GFRα4-sol can promote the phosphorylation of RET, trigger the activation of the PI-3K/AKT pathway, induce neuronal differentiation and support the survival of CGN. Our findings are in line with a recent publication showing the GFRα4-sol might contribute to the inherited cancer syndrome multiple endocrine neoplasia type 2. Our data provide an explanation to how GFRα4-sol may cause or modify the disease. Mammalian GFRα4 receptors all lack the first Cys-rich domain which is present in other GFRα receptors. In the final part of my work I have studied the function of this particular domain. I created a truncated GFRα1 construct lacking the first Cys-rich domain. Using binding assays in both cellular and cell-free systems, phosphorylation assays with RET, as well as neurite outgrowth assays, we found that the first Cys-rich domain contributes to an optimal function of GFRα1, by stabilizing the interaction between GDNF and GFRα1.

The viral coat protein is regulated by HSP70 and HSP40 in Potato virus A infection

Plus-stranded (plus) RNA viruses multiply within a cellular environment as tightly integrated units and rely on the genetic information carried within their genomes for multiplication and, hence, persistence. The minimal genomes of plus RNA viruses are unable to encode the molecular machineries that are required for virus multiplication. This sets requisites for the virus, which must form compatible interactions with host components during multiplication to successfully utilize primary metabolites as building blocks or metabolic energy, and to divert the protein synthesis machinery for production of viral proteins. In fact, the emerging picture of a virus-infected cell displays tight integration with the virus, from simple host and virus protein interactions through to major changes in the physiological state of the host cell. This study set out to develop a method for the identification of host components, mainly host proteins, that interact with proteins of Potato virus A (PVA; Potyvirus) during infection. This goal was approached by developing affinity-tag based methods for the purification of viral proteins complexed with associated host proteins from infected plants. Using this method, host membrane-associated viral ribonucleoprotein (RNP) complexes were obtained, and several host and viral proteins could be identified as components of these complexes. One of the host proteins identified using this strategy was a member of the heat shock protein 70 (HSP70) family, and this protein was chosen for further analysis. To enable the analysis of viral gene expression, a second method was developed based on Agrobacterium-mediated virus genome delivery into plant cells, and detection of virally expressed Renilla luciferase (RLUC) as a quantitative measure of viral gene expression. Using this method, it was observed that down-regulation of HSP70 caused a PVA coat protein (CP)-mediated defect associated with replication. Further experimentation suggested that CP can inhibit viral gene expression and that a distinct translational activity coupled to replication, referred to as replication-associated translation (RAT), exists. Unlike translation of replication-deficient viral RNA, RAT was dependent on HSP70 and its co-chaperone CPIP. HSP70 and CPIP together regulated CP turnover by promoting its modification by ubiquitin. Based on these results, an HSP70 and CPIP-driven mechanism that functions to regulate CP during viral RNA replication and/or translation is proposed, possibly to prevent premature particle assembly caused by CP association with viral RNA.