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.

Lipid-Containing Icosahedral dsDNA Bacteriophages: Entry, Exit and Structure

In this thesis three icosahedral lipid-containing double-stranded (ds) deoxyribonucleic acid (DNA) bacteriophages have been studied: PRD1, Bam35 and P23-77. The work focuses on the entry, exit and structure of the viruses. PRD1 is the type member of the Tectiviridae family, infecting a variety of Gram-negative bacteria. The PRD1 receptor binding complex, consisting of the penton protein P31, the spike protein P5 and the receptor binding protein P2 recognizes a specific receptor on the host surface. In this study we found that the transmembrane protein P16 has an important stabilization function as the fourth member of the receptor binding complex and protein P16 may have a role in the formation of a tubular membrane structure, which is needed in the ejection of the genome into the cell. Phage Bam35 (Tectiviridae), which infects Gram-positive hosts, has been earlier found to resemble PRD1 in morphology and genome organization The uncharacterized early and late events in the Bam35 life cycle were studied by electrochemical methods. Physiological changes in the beginning of the infection were found to be similar in both lysogenic and nonlysogenic cell lines, Bam35 inducing a temporal decrease of membrane voltage and K+ efflux. At the end of the infection cycle physiological changes were observed only in the nonlysogenic cell line. The strong K+ efflux 40 min after infection and the induced premature cell lysis propose that Bam35 has a similar holin-endolysin lysis system to that of PRD1. Thermophilic icosahedral dsDNA Thermus phages P23-65H, P23-72 and P23-77 have been proposed to belong to the Tectiviridae family. In this study these phages were compared to each other. Analysis of structural protein patterns and stability revealed these phages to be very similar but not identical. The most stable of the studied viruses, P23-77, was further analyzed in more detail. Cryo-electron microscopy and three-dimensional image reconstruction was used to determine the structure of virus to 14 Å resolution. Results of thin layer chromatography for neutral lipids together with analysis of the three dimensional reconstruction of P23-77 virus particle revealed the presence of an internal lipid membrane. The overall capsid architecture of P23-77 is similar to PRD1 and Bam35, but most closely it resembles the structure of the capsid of archaeal virus SH1. This complicates the classification of dsDNA, internal lipid-containing icosahedral viruses.

Phonological Mean Length of Utterance as an indicator of typical and impaired phonological development in children acquiring Finnish

The aim of this study was to examine the applicability of the Phonological Mean Length of Utterance (pMLU) method to the data of children acquiring Finnish, for both typically developing children and children with a Specific Language Impairment (SLI). Study I examined typically developing children at the end of the one-word stage (N=17, mean age 1;8), and Study II analysed children s (N=5) productions in a follow-up study with four assessment points (ages 2;0, 2;6, 3;0, 3;6). Study III was carried out in the form of a review article that examined recent research on the phonological development of children acquiring Finnish and compared the results with general trends and cross-linguistic findings in phonological development. Study IV included children with SLI (N=4, mean age 4;10) and age-matched peers. The analyses in Studies I, II and IV were made using the quantitative pMLU method. In the pMLU method, pMLU values are counted for both the words that the children targeted (so-called target words) and the words produced by the children. When the child s average pMLU value was divided with the average target word pMLU value, it is possible to examine that child s accuracy in producing the words with the Whole-Word Proximity (PWP) value. In addition, the number of entirely correctly produced words is counted to obtain the Whole-Word Correctness (PWC) value. Qualitative analyses were carried out in order to examine how the children s phoneme inventories and deficiencies in phonotactics would explain the observed pMLU, PWP and PWC values. The results showed that the pMLU values for children acquiring Finnish were relatively high already at the end of the one-word stage (Study I). The values were found to reflect the characteristics of the ambient language. Typological features that lead to cross-linguistic differences in pMLU values were also observed in the review article (Study III), which noted that in the course of phonological acquisition there are a large number of language-specific phenomena and processes. Study II indicated that overall the children s phonological development during the follow-up period was reflected in the pMLU, PWP and PWC values, although the method showed limitations in detecting qualitative differences between the children. Correct vowels were not scored in the pMLU counts, which led to some misleadingly high pMLU and PWP results: vowel errors were only reflected in the PWC values. Typically developing children in Study II reached the highest possible pMLU results already around age 3;6. At the same time, the differences between the children with SLI and age-matched peers in the pMLU values were very prominent (Study IV). The values for the children with SLI were similar to the ones reported for two-year-old children. Qualitative analyses revealed that the phonologies of the children with SLI largely resembled the ones of younger, typically developing children. However, unusual errors were also witnessed (e.g., vowel errors, omissions of word-initial stops, consonants added to the initial position in words beginning with a vowel). This dissertation provides an application of a new tool for quantitative phonological assessment and analysis in children acquiring Finnish. The preliminary results suggest that, with some modifications, the pMLU method can be used to assess children s phonological development and that it has some advantages compared to the earlier, segment-oriented approaches. Qualitative analyses complemented the pMLU s observations on the children s phonologies. More research is needed in order to verify the levels of the pMLU, PWP and PWC values in children acquiring Finnish.