Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Computational Techniques for Haplotype Inference and for Local Alignment Significance

This thesis which consists of an introduction and four peer-reviewed original publications studies the problems of haplotype inference (haplotyping) and local alignment significance. The problems studied here belong to the broad area of bioinformatics and computational biology. The presented solutions are computationally fast and accurate, which makes them practical in high-throughput sequence data analysis. Haplotype inference is a computational problem where the goal is to estimate haplotypes from a sample of genotypes as accurately as possible. This problem is important as the direct measurement of haplotypes is difficult, whereas the genotypes are easier to quantify. Haplotypes are the key-players when studying for example the genetic causes of diseases. In this thesis, three methods are presented for the haplotype inference problem referred to as HaploParser, HIT, and BACH. HaploParser is based on a combinatorial mosaic model and hierarchical parsing that together mimic recombinations and point-mutations in a biologically plausible way. In this mosaic model, the current population is assumed to be evolved from a small founder population. Thus, the haplotypes of the current population are recombinations of the (implicit) founder haplotypes with some point--mutations. HIT (Haplotype Inference Technique) uses a hidden Markov model for haplotypes and efficient algorithms are presented to learn this model from genotype data. The model structure of HIT is analogous to the mosaic model of HaploParser with founder haplotypes. Therefore, it can be seen as a probabilistic model of recombinations and point-mutations. BACH (Bayesian Context-based Haplotyping) utilizes a context tree weighting algorithm to efficiently sum over all variable-length Markov chains to evaluate the posterior probability of a haplotype configuration. Algorithms are presented that find haplotype configurations with high posterior probability. BACH is the most accurate method presented in this thesis and has comparable performance to the best available software for haplotype inference. Local alignment significance is a computational problem where one is interested in whether the local similarities in two sequences are due to the fact that the sequences are related or just by chance. Similarity of sequences is measured by their best local alignment score and from that, a p-value is computed. This p-value is the probability of picking two sequences from the null model that have as good or better best local alignment score. Local alignment significance is used routinely for example in homology searches. In this thesis, a general framework is sketched that allows one to compute a tight upper bound for the p-value of a local pairwise alignment score. Unlike the previous methods, the presented framework is not affeced by so-called edge-effects and can handle gaps (deletions and insertions) without troublesome sampling and curve fitting.

Primordial Isocurvature Perturbation in Light of CMB and LSS Data

The increased accuracy in the cosmological observations, especially in the measurements of the comic microwave background, allow us to study the primordial perturbations in grater detail. In this thesis, we allow the possibility for a correlated isocurvature perturbations alongside the usual adiabatic perturbations. Thus far the simplest six parameter \Lambda CDM model has been able to accommodate all the observational data rather well. However, we find that the 3-year WMAP data and the 2006 Boomerang data favour a nonzero nonadiabatic contribution to the CMB angular power sprctrum. This is primordial isocurvature perturbation that is positively correlated with the primordial curvature perturbation. Compared with the adiabatic \Lambda CMD model we have four additional parameters describing the increased complexity if the primordial perturbations. Our best-fit model has a 4% nonadiabatic contribution to the CMB temperature variance and the fit is improved by \Delta\chi^2 = 9.7. We can attribute this preference for isocurvature to a feature in the peak structure of the angular power spectrum, namely, the widths of the second and third acoustic peak. Along the way, we have improved our analysis methods by identifying some issues with the parametrisation of the primordial perturbation spectra and suggesting ways to handle these. Due to the improvements, the convergence of our Markov chains is improved. The change of parametrisation has an effect on the MCMC analysis because of the change in priors. We have checked our results against this and find only marginal differences between our parametrisation.

Matrix Decomposition Methods for Data Mining : Computational Complexity and Algorithms

Matrix decompositions, where a given matrix is represented as a product of two other matrices, are regularly used in data mining. Most matrix decompositions have their roots in linear algebra, but the needs of data mining are not always those of linear algebra. In data mining one needs to have results that are interpretable -- and what is considered interpretable in data mining can be very different to what is considered interpretable in linear algebra. --- The purpose of this thesis is to study matrix decompositions that directly address the issue of interpretability. An example is a decomposition of binary matrices where the factor matrices are assumed to be binary and the matrix multiplication is Boolean. The restriction to binary factor matrices increases interpretability -- factor matrices are of the same type as the original matrix -- and allows the use of Boolean matrix multiplication, which is often more intuitive than normal matrix multiplication with binary matrices. Also several other decomposition methods are described, and the computational complexity of computing them is studied together with the hardness of approximating the related optimization problems. Based on these studies, algorithms for constructing the decompositions are proposed. Constructing the decompositions turns out to be computationally hard, and the proposed algorithms are mostly based on various heuristics. Nevertheless, the algorithms are shown to be capable of finding good results in empirical experiments conducted with both synthetic and real-world data.

Univariate, bivariate, and multivariate methods in corpus-based lexicography : A study of synonymy

In this dissertation, I present an overall methodological framework for studying linguistic alternations, focusing specifically on lexical variation in denoting a single meaning, that is, synonymy. As the practical example, I employ the synonymous set of the four most common Finnish verbs denoting THINK, namely ajatella, miettiä, pohtia and harkita ‘think, reflect, ponder, consider’. As a continuation to previous work, I describe in considerable detail the extension of statistical methods from dichotomous linguistic settings (e.g., Gries 2003; Bresnan et al. 2007) to polytomous ones, that is, concerning more than two possible alternative outcomes. The applied statistical methods are arranged into a succession of stages with increasing complexity, proceeding from univariate via bivariate to multivariate techniques in the end. As the central multivariate method, I argue for the use of polytomous logistic regression and demonstrate its practical implementation to the studied phenomenon, thus extending the work by Bresnan et al. (2007), who applied simple (binary) logistic regression to a dichotomous structural alternation in English. The results of the various statistical analyses confirm that a wide range of contextual features across different categories are indeed associated with the use and selection of the selected think lexemes; however, a substantial part of these features are not exemplified in current Finnish lexicographical descriptions. The multivariate analysis results indicate that the semantic classifications of syntactic argument types are on the average the most distinctive feature category, followed by overall semantic characterizations of the verb chains, and then syntactic argument types alone, with morphological features pertaining to the verb chain and extra-linguistic features relegated to the last position. In terms of overall performance of the multivariate analysis and modeling, the prediction accuracy seems to reach a ceiling at a Recall rate of roughly two-thirds of the sentences in the research corpus. The analysis of these results suggests a limit to what can be explained and determined within the immediate sentential context and applying the conventional descriptive and analytical apparatus based on currently available linguistic theories and models. The results also support Bresnan’s (2007) and others’ (e.g., Bod et al. 2003) probabilistic view of the relationship between linguistic usage and the underlying linguistic system, in which only a minority of linguistic choices are categorical, given the known context – represented as a feature cluster – that can be analytically grasped and identified. Instead, most contexts exhibit degrees of variation as to their outcomes, resulting in proportionate choices over longer stretches of usage in texts or speech.

Evaluation of methods and applications for behavioural profiling of transgenic mice

During the last 10-15 years interest in mouse behavioural analysis has evolved considerably. The driving force is development in molecular biological techniques that allow manipulation of the mouse genome by changing the expression of genes. Therefore, with some limitations it is possible to study how genes participate in regulation of physiological functions and to create models explaining genetic contribution to various pathological conditions. The first aim of our study was to establish a framework for behavioural phenotyping of genetically modified mice. We established comprehensive battery of tests for the initial screening of mutant mice. These included tests for exploratory and locomotor activity, emotional behaviour, sensory functions, and cognitive performance. Our interest was in the behavioural patterns of common background strains used for genetic manipulations in mice. Additionally we studied the behavioural effect of sex differences, test history, and individual housing. Our findings highlight the importance of careful consideration of genetic background for analysis of mutant mice. It was evident that some backgrounds may mask or modify the behavioural phenotype of mutants and thereby lead to false positive or negative findings. Moreover, there is no universal strain that is equally suitable for all tests, and using different backgrounds allows one to address possible phenotype modifying factors. We discovered that previous experience affected performance in several tasks. The most sensitive traits were the exploratory and emotional behaviour, as well as motor and nociceptive functions. Therefore, it may be essential to repeat some of the tests in naïve animals for assuring the phenotype. Social isolation for a long time period had strong effects on exploratory behaviour, but also on learning and memory. All experiments revealed significant interactions between strain and environmental factors (test history or housing condition) indicating genotype-dependent effects of environmental manipulations. Several mutant line analyses utilize this information. For example, we studied mice overexpressing as well as those lacking extracellular matrix protein heparin-binding growth-associated molecule (HB-GAM), and mice lacking N-syndecan (a receptor for HB-GAM). All mutant mice appeared to be fertile and healthy, without any apparent neurological or sensory defects. The lack of HB-GAM and N-syndecan, however, significantly reduced the learning capacity of the mice. On the other hand, overexpression of HB-GAM resulted in facilitated learning. Moreover, HB-GAM knockout mice displayed higher anxiety-like behaviour, whereas anxiety was reduced in HB-GAM overexpressing mice. Changes in hippocampal plasticity accompanied the behavioural phenotypes. We conclude that HB-GAM and N-syndecan are involved in the modulation of synaptic plasticity in hippocampus and play a role in regulation of anxiety- and learning-related behaviour.

Mu in vitro Transposition Technology in Functional Genetics and Genomics: Applications on Mouse and Bacteriophages

Transposons are mobile elements of genetic material that are able to move in the genomes of their host organisms using a special form of recombination called transposition. Bacteriophage Mu was the first transposon for which a cell-free in vitro transposition reaction was developed. Subsequently, the reaction has been refined and the minimal Mu in vitro reaction is useful in the generation of comprehensive libraries of mutant DNA molecules that can be used in a variety of applications. To date, the functional genetics applications of Mu in vitro technology have been subjected to either plasmids or genomic regions and entire genomes of viruses cloned on specific vectors. This study expands the use of Mu in vitro transposition in functional genetics and genomics by describing novel methods applicable to the targeted transgenesis of mouse and the whole-genome analysis of bacteriophages. The methods described here are rapid, efficient, and easily applicable to a wide variety of organisms, demonstrating the potential of the Mu transposition technology in the functional analysis of genes and genomes. First, an easy-to-use, rapid strategy to generate construct for the targeted mutagenesis of mouse genes was developed. To test the strategy, a gene encoding a neuronal K+/Cl- cotransporter was mutagenised. After a highly efficient transpositional mutagenesis, the gene fragments mutagenised were cloned into a vector backbone and transferred into bacterial cells. These constructs were screened with PCR using an effective 3D matrix system. In addition to traditional knock-out constructs, the method developed yields hypomorphic alleles that lead into reduced expression of the target gene in transgenic mice and have since been used in a follow-up study. Moreover, a scheme is devised to rapidly produce conditional alleles from the constructs produced. Next, an efficient strategy for the whole-genome analysis of bacteriophages was developed based on the transpositional mutagenesis of uncloned, infective virus genomes and their subsequent transfer into susceptible host cells. Mutant viruses able to produce viable progeny were collected and their transposon integration sites determined to map genomic regions nonessential to the viral life cycle. This method, applied here to three very different bacteriophages, PRD1, ΦYeO3 12, and PM2, does not require the target genome to be cloned and is directly applicable to all DNA and RNA viruses that have infective genomes. The method developed yielded valuable novel information on the three bacteriophages studied and whole-genome data can be complemented with concomitant studies on individual genes. Moreover, end-modified transposons constructed for this study can be used to manipulate genomes devoid of suitable restriction sites.

Novel matrix metalloproteinases in intestinal inflammation and in cancer of the gastrointestinal tract

Matrix metalloproteinases (MMPs) comprise a family of 23 zinc-dependent human endopeptidases that can degrade virtually all components of the extracellular matrix (ECM). They are classified into eight subgroups according to their structure and into six subgroups based on their substrate-specificity. MMPs have been implicated in inflammation, tissue destruction, cell migration, arthritis, vascular remodeling, angiogenesis, and tumor growth and invasion. MMPs are inhibited by their natural inhibitors, tissue inhibitors of metalloproteinases (TIMPs). Different MMPs function in the same tasks depending on the tissue or cancer subtype. I investigated the role of recently discovered MMPs, especially MMPs-19 and -26, in intestinal inflammation, in intestinal and cutaneous wound healing, and in intestinal cancer. Several MMPs and TIMPs were studied to determine their exact location at tissue level and to obtain information on possible functions of MMPs in such tissues and diseases as the healthy intestine, inflammatory bowel disease (IBD), neonatal necrotizing enterocolitis (NEC), pyoderma gangrenosum (PG), and colorectal as well as pancreatic cancers. In latent celiac disease (CD), I attempted to identify markers to predict later onset of CD in children and adolescents. The main methods used were immunohistochemistry, in situ hybridization, and Taqman RT-PCR. My results show that MMP-26 is important for re-epithelialization in intestinal and cutaneous wound healing. In colon and pancreatic cancers, MMP-26 seems to be a marker of invasive potential, although it is not itself expressed at the invasive front. MMP-21 is upregulated in pancreatic cancer and may be associated with tumor differentiation. MMPs-19 and -28 are associated with normal tissue turnover in the intestine, but they disappear in tumor progression as if they were protective markers . MMP-12 is an essential protease in intestinal inflammation and tissue destruction, as seen here in NEC and in previous CD studies. In patients with type 1 diabetes (T1D), MMPs-1, -3, and -12 were upregulated in the intestinal mucosa. Furthermore, MMP-7 was strongly elevated in NEC. In a model of aberrant wound repair, PG, MMPs-8, -9, and 10 and TNFα may promote ECM destruction, while absence of MMP-1 and MMP-26 from keratinocytes retards re-epithelialization. Based on my results, I suggest MMP-26 to be considered a putative marker for poor prognosis in pancreatic and colon cancer. However, since it functions differently in various tissues and tumor subtypes, this use cannot be generalized. Furthermore, MMP-26 is a beneficial marker for wound healing if expressed by migrating epithelial cells. MMP-12 expression in latent CD patients warrants research in a larger patient population to confirm its role as a specific marker for CD in pathologically indistinct cases. MMP-7 should be considered one of the most crucial proteases in NEC-associated tissue destruction; hence, specific inhibitors of this MMP are worth investigating. In PG, TNFα inhibitors are potential therapeutic agents, as shown already in clinical trials. In conclusion, studies of several MMPs in specific diseases and in healthy tissues are needed to elucidate their roles at the tissue level. MMPs and TIMPs are not exclusively destructive or reparative in tissues. They seem to function differently in different tissues. To identify selective MMP inhibitors, we must thoroughly understand the MMP profile (degradome) and their functions in various organs not to interfere with normal reparative functions during wound repair or beneficial host-response effects during cancer initiation and growth.

Matrix metalloproteinases as biomarkers in premalignant and malignant tumors of the human skin

The incidence of non-melanoma skin cancer is increasing worldwide. Basal cell carcinoma followed by squamous cell carcinoma and malignant melanoma are the most frequent skin tumors. Immunosuppressed patients have an increased risk of neoplasia, of which non-melanoma skin cancer is the most common. Matrix metalloproteinases (MMPs) are proteolytic enzymes that collectively are capable of degrading virtually all components of the extracellular matrix. MMPs can also process substrates distinct from extracellular matrix proteins and influence cell proliferation, differentiation, angiogenesis, and apoptosis. MMP activity is regulated by their natural inhibitors, tissue inhibitors of metallopro-teinases (TIMPs). In this study, the expression patterns of MMPs, TIMPs, and certain cancer-related molecules were investigated in premalignant and malignant lesions of the human skin. As methods were used immunohistochemisty, in situ hybridization, and reverse transcriptase polymerase chain reaction (RT-PCR) from the cell cultures. Our aim was to evaluate the expression pattern of MMPs in extramammary Paget's disease in order to find markers for more advanced tumors, as well as to shed light on the origin of this rare neoplasm. Novel MMPs -21, -26, and -28 were studied in melanoma cell culture, in primary cutaneous melanomas, and their sentinel nodes. The MMP expression profile in keratoacanthomas and well-differentiated squamous cell carcinomas was analyzed to find markers to differentiate benign keratinocyte hyperproliferation from malignantly transformed cells. Squamous cell carcinomas of immunosuppressed organ transplant recipients were compared to squamous cell carcinomas of matched immunocompetent controls to investigate the factors explaining their more aggressive nature. We found that MMP-7 and -19 proteins are abundant in extramammary Paget's disease and that their presence may predict an underlying adenocarcinoma in these patients. In melanomas, MMP-21 was upregulated in early phases of melanoma progression, but disappeared from the more aggressive tumors with lymph node metastases. The presence of MMP-13 in primary melanomas and lymph node metastases may relate to more aggressive disease. In keratoacanthomas, the expression of MMP-7 and -9 is rare and therefore should raise a suspicion of well-differentiated squamous cell carcinomas. Furthermore, MMP-19 and p16 were observed in benign keratinocyte hyperproliferation of keratoacanthomas, whereas they were generally lost from malignant keratinocytes of SCCs. MMP-26 staining was significantly stronger in squamous cell carcinomas and Bowen s disease samples of organ transplant recipients and it may contribute to the more aggressive nature of squamous cell carcinomas in immunosuppressed patients. In addition, the staining for MMP-9 was significantly stronger in macrophages surrounding the tumors of the immunocompetent group and in neutrophils of those patients on cyclosporin medication. In conclusion, based on our studies, MMP-7 and -19 might serve as biomarkers for more aggressive extramammary Paget's disease and MMP-21 for malignant transformation of melanocytes. MMP -7, -9, and -26, however, could play an important role in the pathobiology of keratinocyte derived malignancies.

Basement membrane and provisional matrix components (collagen type IV, laminins and von Willebrand factor) in rheumatoid arthritis and Sjögren s syndrome

The aim of this study was twofold- Firstly, to determine the composition of the type IV collagen which are the major components of the basement membrane (BM), in the synovial lining of the rheumatoid arthritis (RA) patient and in the BM in the labial salivary gland of the Sjögrens syndrome (SS) patient. Secondly, this thesis aimed to investigate the role of the BM component laminin α4 and laminin α5 in the migration of neutrophils from the blood vessels thorough the synovial lining layer into synovial fluid and the presence of vWF in the microvasculature of labial salivary gland in SS. Our studies showed that certain α chains type IV collagen are low in RA compared to control synovial linings, while laminin α5 exhibited a pattern of low expression regions at the synovial lining interface towards the joint cavity and fluid. Also, high numbers of macrophage-like lining cells containing MMP-9 were found in the lining. MMP-9 was also found in the synovial fluid. Collagen α1/2 (IV) mRNA was found to be present in high amount compared to the other α(IV) chains and also showed intense labelling in immunohistochemical staining in normal and SS patients. In healthy glands α5(IV) and α6(IV) chains were found to be continuous around ducts but discontinuous around acini. The α5(IV) and α6(IV) mRNAs were present in LSG explants and HSG cell line, while in SS these chains seemed to be absent or appear only in patches around the ductal BM and tended to be absent around acini in immunohistochemical staining, indicating that their synthesis and/or degradation seemed to be locally regulated around acinar cells. The provisional matrix component vWF serves as a marker of vascular damage. Microvasculature in SS showed signs of focal damage which in turn might impair arteriolar feeding, capillary transudation and venular drainage of blood. However, capillary density was not decreased but rather increased, perhaps as a result of angiogenesis compensatory to microvascular damage. Microvascular involvement of LSG may contribute to the pathogenesis of this syndrome. This twofold approach allows us to understand the intricate relation between the ECM components and the immunopathological changes that occur during the pathogenesis of these inflammatory rheumatic disease processes. Also notably this study highlights the importance of maintaining a healthy ECM to prevent the progression or possibly allow reversal of the disease to a considerable level. Furthermore, it can be speculated that a healthy BM could quarantine the inflamed region or in case of cancer cells barricade the movement of malignant cells thereby preventing further spread to the surrounding areas. This understanding can be further applied to design appropriate drugs which act specifically to maintain a proper BM/BM like intercellular matrix composition.

Molecular basis for the development of diagnostic methods and specific treatment modalities for idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is an interstitial lung disease with unknown aetiology and poor prognosis. IPF is characterized by alveolar epithelial damage that leads tissue remodelling and ultimately to the loss of normal lung architecture and function. Treatment has been focused on anti-inflammatory therapies, but due to their poor efficacy new therapeutic modalities are being sought. There is a need for early diagnosis and also for differential diagnostic markers for IPF and other interstitial lung diseases. The study utilized patient material obtained from bronchoalveolar lavage (BAL), diagnostic biopsies or lung transplantation. Human pulmonary fibroblast cell cultures were propagated and asbestos-induced pulmonary fibrosis in mice was used as an experimental animal model of IPF. The possible markers for IPF were scanned by immunohistochemistry, RT-PCR, ELISA and western blot. Matrix metalloproteinases (MMPs) are proteolytic enzymes that participate in tissue remodelling. Microarray studies have introduced potential markers that could serve as additional tools for the assessment of IPF and one of the most promising was MMP 7. MMP-7 protein levels were measured in the BAL fluid of patients with idiopathic interstitial lung diseases or idiopathic cough. MMP-7 was however similarly elevated in the BAL fluid of all these disorders and thus cannot be used as a differential diagnostic marker for IPF. Activation of transforming growth factor (TGF)-ß is considered to be a key element in the progression of IPF. Bone morphogenetic proteins (BMP) are negative regulators of intracellular TGF-ß signalling and BMP-4 signalling is in turn negatively regulated by gremlin. Gremlin was found to be highly upregulated in the IPF lungs and IPF fibroblasts. Gremlin was detected in the thickened IPF parenchyma and endothelium of small capillaries, whereas in non-specific interstitial pneumonia it localized predominantly in the alveolar epithelium. Parenchymal gremlin immunoreactivity might indicate IPF-type interstitial pneumonia. Gremlin mRNA levels were higher in patients with end-stage fibrosis suggesting that gremlin might be a marker for more advanced disease. Characterization of the fibroblastic foci in the IPF lungs showed that immunoreactivity to platelet-derived growth factor (PDGF) receptor-α and PDGF receptor-β was elevated in IPF parenchyma, but the fibroblastic foci showed only minor immunoreactivity to the PDGF receptors or the antioxidant peroxiredoxin II. Ki67 positive cells were also observed predominantly outside the fibroblastic foci, suggesting that the fibroblastic foci may not be composed of actively proliferating cells. When inhibition of profibrotic PDGF-signalling by imatinib mesylate was assessed, imatinib mesylate reduced asbestos-induced pulmonary fibrosis in mice as well as human pulmonary fibroblast migration in vitro but it had no effect on the lung inflammation.

Enzymatic cross-linking of β-casein and its impact on digestibility and allergenicity

Protein modification via enzymatic cross-linking is an attractive way for altering food structure so as to create products with increased quality and nutritional value. These modifications are expected to affect not only the structure and physico-chemical properties of proteins but also their physiological characteristics, such as digestibility in the GI-tract and allergenicity. Protein cross-linking enzymes such as transglutaminases are currently commercially available, but also other types of cross-linking enzymes are being explored intensively. In this study, enzymatic cross-linking of β-casein, the most abundant bovine milk protein, was studied. Enzymatic cross-linking reactions were performed by fungal Trichoderma reesei tyrosinase (TrTyr) and the performance of the enzyme was compared to that of transglutaminase from Streptoverticillium mobaraense (Tgase). Enzymatic cross-linking reactions were followed by different analytical techniques, such as size exclusion chromatography -Ultra violet/Visible multi angle light scattering (SEC-UV/Vis-MALLS), phosphorus nuclear magnetic resonance spectroscopy (31P-NMR), atomic force (AFM) and matrix-assisted laser desorption/ionisation-time of flight mass spectrometry (MALDI-TOF MS). The research results showed that in both cases cross-linking of β-casein resulted in the formation of high molecular mass (MM ca. 1 350 kg mol-1), disk-shaped nanoparticles when the highest enzyme dosage and longest incubation times were used. According to SEC-UV/Vis-MALLS data, commercial β-casein was cross-linked almost completely when TrTyr and Tgase were used as cross-linking enzymes. In the case of TrTyr, high degree of cross-linking was confirmed by 31P-NMR where it was shown that 91 % of the tyrosine side-chains were involved in the cross-linking. The impact of enzymatic cross-linking of β-casein on in vitro digestibility by pepsin was followed by various analytical techniques. The research results demonstrated that enzymatically cross-linked β-casein was stable under the acidic conditions present in the stomach. Furthermore, it was found that cross-linked β-casein was more resistant to pepsin digestion when compared to that of non modified β-casein. The effects of enzymatic cross-linking of β-casein on allergenicity were also studied by different biochemical test methods. On the basis of the research results, enzymatic cross-linking decreased allergenicity of native β-casein by 14 % when cross-linked by TrTyr and by 6 % after treatment by Tgase. It can be concluded that in addition to the basic understanding of the reaction mechanism of TrTyr on protein matrix, the research results obtained in this study can have high impact on various applications like food, cosmetic, medical, textile and packing sectors.

Inference on Cointegration in Vector Autoregressive Models (summary section only)

In the thesis we consider inference for cointegration in vector autoregressive (VAR) models. The thesis consists of an introduction and four papers. The first paper proposes a new test for cointegration in VAR models that is directly based on the eigenvalues of the least squares (LS) estimate of the autoregressive matrix. In the second paper we compare a small sample correction for the likelihood ratio (LR) test of cointegrating rank and the bootstrap. The simulation experiments show that the bootstrap works very well in practice and dominates the correction factor. The tests are applied to international stock prices data, and the .nite sample performance of the tests are investigated by simulating the data. The third paper studies the demand for money in Sweden 1970—2000 using the I(2) model. In the fourth paper we re-examine the evidence of cointegration between international stock prices. The paper shows that some of the previous empirical results can be explained by the small-sample bias and size distortion of Johansen’s LR tests for cointegration. In all papers we work with two data sets. The first data set is a Swedish money demand data set with observations on the money stock, the consumer price index, gross domestic product (GDP), the short-term interest rate and the long-term interest rate. The data are quarterly and the sample period is 1970(1)—2000(1). The second data set consists of month-end stock market index observations for Finland, France, Germany, Sweden, the United Kingdom and the United States from 1980(1) to 1997(2). Both data sets are typical of the sample sizes encountered in economic data, and the applications illustrate the usefulness of the models and tests discussed in the thesis.

Occupational chronic solvent encephalopathy in Finland 1995 - 2007 : incidence and diagnostic methods

The occurrence of occupational chronic solvent encephalopathy (CSE) seems to decrease, but still every year reveals new cases. To prevent CSE and early retirement of solvent-exposed workers, actions should focus on early CSE detection and diagnosis. Identifying the work tasks and solvent exposure associated with high risk for CSE is crucial. Clinical and exposure data of all the 128 cases diagnosed with CSE as an occupational disease in Finland during 1995-2007 was collected from the patient records at the Finnish Institute of Occupational Health (FIOH) in Helsinki. The data on the number of exposed workers in Finland were gathered from the Finnish Job-exposure Matrix (FINJEM) and the number of employed from the national workforce survey. We analyzed the work tasks and solvent exposure of CSE patients and the findings in brain magnetic resonance imaging (MRI), quantitative electroencephalography (QEEG), and event-related potentials (ERP). The annual number of new cases diminished from 18 to 3, and the incidence of CSE decreased from 8.6 to 1.2 / million employed per year. The highest incidence of CSE was in workers with their main exposure to aromatic hydrocarbons; during 1995-2006 the incidence decreased from 1.2 to 0.3 / 1 000 exposed workers per year. The work tasks with the highest incidence of CSE were floor layers and lacquerers, wooden surface finishers, and industrial, metal, or car painters. Among 71 CSE patients, brain MRI revealed atrophy or white matter hyperintensities or both in 38% of the cases. Atrophy which was associated with duration of exposure was most frequently located in the cerebellum and in the frontal or parietal brain areas. QEEG in a group of 47 patients revealed increased power of the theta band in the frontal brain area. In a group of 86 patients, the P300 amplitude of auditory ERP was decreased, but at individual level, all the amplitude values were classified as normal. In 11 CSE patients and 13 age-matched controls, ERP elicited by a multimodal paradigm including an auditory, a visual detection, and a recognition memory task under single and dual-task conditions corroborated the decrease of auditory P300 amplitude in CSE patients in single-task condition. In dual-task conditions, the auditory P300 component was, more often in patients than in controls, unrecognizable. Due to the paucity and non-specificity of the findings, brain MRI serves mainly for differential diagnostics in CSE. QEEG and auditory P300 are insensitive at individual level and not useful in the clinical diagnostics of CSE. A multimodal ERP paradigm may, however, provide a more sensitive method to diagnose slight cognitive disturbances such as CSE.

The decline of malaria in Finland – the impact of the vector and social variables

Background: Malaria was prevalent in Finland in the 18th century. It declined slowly without deliberate counter-measures and the last indigenous case was reported in 1954. In the present analysis of indigenous malaria in Finland, an effort was made to construct a data set on annual malaria cases of maximum temporal length to be able to evaluate the significance of different factors assumed to affect malaria trends. Methods: To analyse the long-term trend malaria statistics were collected from 1750–2008. During that time, malaria frequency decreased from about 20,000 – 50,000 per 1,000,000 people to less than 1 per 1,000,000 people. To assess the cause of the decline, a correlation analysis was performed between malaria frequency per million people and temperature data, animal husbandry, consolidation of land by redistribution and household size. Results: Anopheles messeae and Anopheles beklemishevi exist only as larvae in June and most of July. The females seek an overwintering place in August. Those that overwinter together with humans may act as vectors. They have to stay in their overwintering place from September to May because of the cold climate. The temperatures between June and July determine the number of malaria cases during the following transmission season. This did not, however, have an impact on the longterm trend of malaria. The change in animal husbandry and reclamation of wetlands may also be excluded as a possible cause for the decline of malaria. The long-term social changes, such as land consolidation and decreasing household size, showed a strong correlation with the decline of Plasmodium. Conclusion: The indigenous malaria in Finland faded out evenly in the whole country during 200 years with limited or no counter-measures or medication. It appears that malaria in Finland was basically a social disease and that malaria trends were strongly linked to changes in human behaviour. Decreasing household size caused fewer interactions between families and accordingly decreasing recolonization possibilities for Plasmodium. The permanent drop of the household size was the precondition for a permanent eradication of malaria.