Discovering hidden structures in molecular data using a Bayesian partition model approach

Advancements in the analysis techniques have led to a rapid accumulation of biological data in databases. Such data often are in the form of sequences of observations, examples including DNA sequences and amino acid sequences of proteins. The scale and quality of the data give promises of answering various biologically relevant questions in more detail than what has been possible before. For example, one may wish to identify areas in an amino acid sequence, which are important for the function of the corresponding protein, or investigate how characteristics on the level of DNA sequence affect the adaptation of a bacterial species to its environment. Many of the interesting questions are intimately associated with the understanding of the evolutionary relationships among the items under consideration. The aim of this work is to develop novel statistical models and computational techniques to meet with the challenge of deriving meaning from the increasing amounts of data. Our main concern is on modeling the evolutionary relationships based on the observed molecular data. We operate within a Bayesian statistical framework, which allows a probabilistic quantification of the uncertainties related to a particular solution. As the basis of our modeling approach we utilize a partition model, which is used to describe the structure of data by appropriately dividing the data items into clusters of related items. Generalizations and modifications of the partition model are developed and applied to various problems. Large-scale data sets provide also a computational challenge. The models used to describe the data must be realistic enough to capture the essential features of the current modeling task but, at the same time, simple enough to make it possible to carry out the inference in practice. The partition model fulfills these two requirements. The problem-specific features can be taken into account by modifying the prior probability distributions of the model parameters. The computational efficiency stems from the ability to integrate out the parameters of the partition model analytically, which enables the use of efficient stochastic search algorithms.

Questions to Artificial Nature: A Philosophical Study of Interdisciplinary Models and their Functions in Scientific Practice

This thesis presents an interdisciplinary analysis of how models and simulations function in the production of scientific knowledge. The work is informed by three scholarly traditions: studies on models and simulations in philosophy of science, so-called micro-sociological laboratory studies within science and technology studies, and cultural-historical activity theory. Methodologically, I adopt a naturalist epistemology and combine philosophical analysis with a qualitative, empirical case study of infectious-disease modelling. This study has a dual perspective throughout the analysis: it specifies the modelling practices and examines the models as objects of research. The research questions addressed in this study are: 1) How are models constructed and what functions do they have in the production of scientific knowledge? 2) What is interdisciplinarity in model construction? 3) How do models become a general research tool and why is this process problematic? The core argument is that the mediating models as investigative instruments (cf. Morgan and Morrison 1999) take questions as a starting point, and hence their construction is intentionally guided. This argument applies the interrogative model of inquiry (e.g., Sintonen 2005; Hintikka 1981), which conceives of all knowledge acquisition as process of seeking answers to questions. The first question addresses simulation models as Artificial Nature, which is manipulated in order to answer questions that initiated the model building. This account develops further the "epistemology of simulation" (cf. Winsberg 2003) by showing the interrelatedness of researchers and their objects in the process of modelling. The second question clarifies why interdisciplinary research collaboration is demanding and difficult to maintain. The nature of the impediments to disciplinary interaction are examined by introducing the idea of object-oriented interdisciplinarity, which provides an analytical framework to study the changes in the degree of interdisciplinarity, the tools and research practices developed to support the collaboration, and the mode of collaboration in relation to the historically mutable object of research. As my interest is in the models as interdisciplinary objects, the third research problem seeks to answer my question of how we might characterise these objects, what is typical for them, and what kind of changes happen in the process of modelling. Here I examine the tension between specified, question-oriented models and more general models, and suggest that the specified models form a group of their own. I call these Tailor-made models, in opposition to the process of building a simulation platform that aims at generalisability and utility for health-policy. This tension also underlines the challenge of applying research results (or methods and tools) to discuss and solve problems in decision-making processes.

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.

Health-Related Quality of Life, Costs and Treatment of Inflammatory Rheumatic Diseases in Routine Practice : With special emphasis on biological drugs

Rheumatoid arthritis (RA) and other chronic inflammatory joint diseases already begin to affect patients health-related quality of life (HRQoL) in the earliest phases of these diseases. In treatment of inflammatory joint diseases, the last two decades have seen new strategies and treatment options introduced. Treatment is started at an earlier phase; combinations of disease-modifying anti-rheumatic drugs (DMARDs) and corticosteroids are used; and in refractory cases new drugs such as tumour necrosis factor (TNF) inhibitors or other biologicals can be started. In patients with new referrals to the Department of Rheumatology of the Helsinki University Central Hospital, we evaluated the 15D and the Stanford Health Assessment Questionnaire (HAQ) results at baseline and approximately 8 months after their first visit. Altogether the analysis included 295 patients with various rheumatic diseases. The mean baseline 15D score (0.822, SD 0.114) was significantly lower than for the age-matched general population (0.903, SD 0.098). Patients with osteoarthritis (OA) and spondyloarthropathies (SPA) reported the poorest HRQoL. In patients with RA and reactive arthritis (ReA) the HRQoL improved in a statistically significant manner during the 8-month follow-up. In addition, a clinically important change appeared in patients with systemic rheumatic diseases. HAQ score improved significantly in patients with RA, arthralgia and fibromyalgia, and ReA. In a study of 97 RA patients treated either with etanercept or adalimumab, we assessed their HRQoL with the RAND 36-Item Health Survey 1.0 (RAND-36) questionnaire. We also analysed changes in clinical parameters and the HAQ. With etanercept and adalimumab, the values of all domains in the RAND-36 questionnaire increased during the first 3 months. The efficacy of each in improving HRQoL was statistically significant, and the drug effects were comparable. Compared to Finnish age- and sex-matched general population values, the HRQoL of the RA patients was significantly lower at baseline and, despite the improvement, remained lower also at follow-up. Our RA patients had long-standing and severe disease that can explain the low HRQoL also at follow-up. In a pharmacoeconomic study of patients treated with infliximab we evaluated medical and work disability costs for patients with chronic inflammatory joint disease during one year before and one year after institution of infliximab treatment. Clinical and economic data for 96 patients with different arthritis diagnoses showed, in all patients, significantly improved clinical and laboratory variables. However, the medical costs increased significantly during the second period by 12 015 (95% confidence interval, 6 496 to 18 076). Only a minimal decrease in work disability costs occurred mean decrease 130 (-1 268 to 1 072). In a study involving a switch from infliximab to etanercept, we investigated the clinical outcome in 49 patients with RA. Reasons for switching were in 42% failure to respond by American College of Rheumatology (ACR) 50% criteria; in 12% adverse event; and in 46% non-medical reasons although the patients had responded to infliximab. The Disease Activity Score with 28 joints examined (DAS28) allowed us to measure patients disease activity and compare outcome between groups based on the reason for switching. In the patients in whom infliximab was switched to etanercept for nonmedical reasons, etanercept continued to suppress disease activity effectively, and 1-year drug survival for etanercept was 77% (95% CI, 62 to 97). In patients in the infliximab failure and adverse event groups, DAS28 values improved significantly during etanercept therapy. However, the 1-year drug survival of etanercept was only 43% (95% CI, 26 to 70) and 50% (95% CI, 33 to 100), respectively. Although the HRQoL of patients with inflammatory joint diseases is significantly lower than that of the general population, use of early and aggressive treatment strategies including TNF-inhibitors can improve patients HRQoL effectively. Further research is needed in finding new treatment strategies for those patients who fail to respond or lose their response to TNF-inhibitors.

X-ray and neutron scattering studies on some nanoscale structures in molecular biology

Scattering of X-rays and neutrons has been applied to the study of nanostructures with interesting biological functions. The systems studied were the protein calmodulin and its complexes, bacterial virus bacteriophage phi6, and the photosynthetic antenna complex from green sulfur bacteria, chlorosome. Information gathered using various structure determination methods has been combined to the low resolution information obtained from solution scattering. Conformational changes in calmodulin-ligand complex were studied by combining the directional information obtained from residual dipole couplings in nuclear magnetic resonance to the size information obtained from small-angle X-ray scattering from solution. The locations of non-structural protein components in a model of bacteriophage phi6, based mainly on electron microscopy, were determined by neutron scattering, deuterium labeling and contrast variation. New data are presented on the structure of the photosynthetic antenna complex of green sulfur bacteria and filamentous anoxygenic phototrophs, also known as the chlorosome. The X-ray scattering and electron cryomicroscopy results from this system are interpreted in the context of a new structural model detailed in the third paper of this dissertation. The model is found to be consistent with the results obtained from various chlorosome containing bacteria. The effect of carotenoid synthesis on the chlorosome structure and self-assembly are studied by carotenoid extraction, biosynthesis inhibition and genetic manipulation of the enzymes involved in carotenoid biosynthesis. Carotenoid composition and content are found to have a marked effect on the structural parameters and morphology of chlorosomes.

Studies of the star-forming structures in the dense interstellar medium: a view by dust extinction

New stars in galaxies form in dense, molecular clouds of the interstellar medium. Measuring how the mass is distributed in these clouds is of crucial importance for the current theories of star formation. This is because several open issues in them, such as the strength of different mechanism regulating star formation and the origin of stellar masses, can be addressed using detailed information on the cloud structure. Unfortunately, quantifying the mass distribution in molecular clouds accurately over a wide spatial and dynamical range is a fundamental problem in the modern astrophysics. This thesis presents studies examining the structure of dense molecular clouds and the distribution of mass in them, with the emphasis on nearby clouds that are sites of low-mass star formation. In particular, this thesis concentrates on investigating the mass distributions using the near infrared dust extinction mapping technique. In this technique, the gas column densities towards molecular clouds are determined by examining radiation from the stars that shine through the clouds. In addition, the thesis examines the feasibility of using a similar technique to derive the masses of molecular clouds in nearby external galaxies. The papers presented in this thesis demonstrate how the near infrared dust extinction mapping technique can be used to extract detailed information on the mass distribution in nearby molecular clouds. Furthermore, such information is used to examine characteristics crucial for the star formation in the clouds. Regarding the use of extinction mapping technique in nearby galaxies, the papers of this thesis show that deriving the masses of molecular clouds using the technique suffers from strong biases. However, it is shown that some structural properties can still be examined with the technique.

Inhomogeneous Structures in Holographic Superfluids I

We begin an investigation of inhomogeneous structures in holographic superfluids. As a first example, we study domain wall like defects in the 3+1 dimensional Einstein-Maxwell-Higgs theory, which was developed as a dual model for a holographic superconductor. In [1], we reported on such "dark solitons" in holographic superfluids. In this work, we present an extensive numerical study of their properties, working in the probe limit. We construct dark solitons for two possible condensing operators, and find that both of them share common features with their standard superfluid counterparts. However, both are characterized by two distinct coherence length scales (one for order parameter, one for charge condensate). We study the relative charge depletion factor and find that solitons in the two different condensates have very distinct depletion characteristics. We also study quasiparticle excitations above the holographic superfluid, and find that the scale of the excitations is comparable to the soliton coherence length scales.

Methods and applications for improving parameter prediction models for stand structures in Finland

This thesis report attempts to improve the models for predicting forest stand structure for practical use, e.g. forest management planning (FMP) purposes in Finland. Comparisons were made between Weibull and Johnson s SB distribution and alternative regression estimation methods. Data used for preliminary studies was local but the final models were based on representative data. Models were validated mainly in terms of bias and RMSE in the main stand characteristics (e.g. volume) using independent data. The bivariate SBB distribution model was used to mimic realistic variations in tree dimensions by including within-diameter-class height variation. Using the traditional method, diameter distribution with the expected height resulted in reduced height variation, whereas the alternative bivariate method utilized the error-term of the height model. The lack of models for FMP was covered to some extent by the models for peatland and juvenile stands. The validation of these models showed that the more sophisticated regression estimation methods provided slightly improved accuracy. A flexible prediction and application for stand structure consisted of seemingly unrelated regression models for eight stand characteristics, the parameters of three optional distributions and Näslund s height curve. The cross-model covariance structure was used for linear prediction application, in which the expected values of the models were calibrated with the known stand characteristics. This provided a framework to validate the optional distributions and the optional set of stand characteristics. Height distribution is recommended for the earliest state of stands because of its continuous feature. From the mean height of about 4 m, Weibull dbh-frequency distribution is recommended in young stands if the input variables consist of arithmetic stand characteristics. In advanced stands, basal area-dbh distribution models are recommended. Näslund s height curve proved useful. Some efficient transformations of stand characteristics are introduced, e.g. the shape index, which combined the basal area, the stem number and the median diameter. Shape index enabled SB model for peatland stands to detect large variation in stand densities. This model also demonstrated reasonable behaviour for stands in mineral soils.