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.

Homogeneous model theory of metric structures

This thesis studies homogeneous classes of complete metric spaces. Over the past few decades model theory has been extended to cover a variety of nonelementary frameworks. Shelah introduced the abstact elementary classes (AEC) in the 1980s as a common framework for the study of nonelementary classes. Another direction of extension has been the development of model theory for metric structures. This thesis takes a step in the direction of combining these two by introducing an AEC-like setting for studying metric structures. To find balance between generality and the possibility to develop stability theoretic tools, we work in a homogeneous context, thus extending the usual compact approach. The homogeneous context enables the application of stability theoretic tools developed in discrete homogeneous model theory. Using these we prove categoricity transfer theorems for homogeneous metric structures with respect to isometric isomorphisms. We also show how generalized isomorphisms can be added to the class, giving a model theoretic approach to, e.g., Banach space isomorphisms or operator approximations. The novelty is the built-in treatment of these generalized isomorphisms making, e.g., stability up to perturbation the natural stability notion. With respect to these generalized isomorphisms we develop a notion of independence. It behaves well already for structures which are omega-stable up to perturbation and coincides with the one from classical homogeneous model theory over saturated enough models. We also introduce a notion of isolation and prove dominance for it.

Optimal maintenance and replacement via a semi-Markov decision model

Optimal bang-coast maintenance policies for a machine, subject to failure, are considered. The approach utilizes a semi-Markov model for the system. A simplified model for modifying the probability of machine failure with maintenance is employed. A numerical example is presented to illustrate the procedure and results.

Homogeneous hydrogenations ofα, β-unsaturated compounds using pentacyanocobaltate (II) anion as catalyst system

Hydrogenation of someα, β-unsaturated carbonyl compounds using potassium pentacyanocobaltate (II), K3Co(CN)5, as a homogeneous catalyst has been investigated. Thus, hydrogenation of 1-carvone (I), mesityl oxide (4), 2-cyclohexenone (8) and benzalacetone (6) afforded the corresponding dihydrocompounds. Hydrogenation ofβ-ionone (10) afforded a mixture of theα, β-dihydrocompounds (14) and (15). In all these cases, it was observed that the reaction proceeded to completion only in the presence of added base. Hydrogenation of 5α-androst-l-en-17β-ol-3-one acetate (19) afforded the saturated compound, 5α-androst-17β-ol-3-one (20) in 60% yield. It was found that other steroid enones and dienones were not reduced by this catalyst system.

Monte Carlo simulation of Ising chains showing several phase transitions

A Monte Carlo simulation of Ising chains with competing short-range and infiniterange interactions has been carried out. Results show that whenever the system does not enter a metastable state, variation of temperature brings about phase transitions in the Ising chain. These phase transitions, except for two sets of interaction strengths, are generally of higher order and involve changes in the long-range order while the short-range order remains unaffected.

Conformation of polypeptide-chains containing both l-residues and D-residues .2. Double-helical structures of poly-ld-peptides

Polypeptides with alternating L- and D-amino acid residues can take up stereochemically satisfactory coaxial double-helical structures, both antiparallel and parallel, which are stabilized by systematic interchain NH O hydrogen bonds. Semiempirical energy calculations over allowed regions of conformational space have yielded the characteristics of these double-helices. There are four possible types of antiparallel double-helices - A3, A4, A5 and A6, with n, the number of LD peptide units per turn, around 2.8, 3.6, 4.5 and 5.5 respectively, while for the parallel double-helices there are two types, P3 and P4, having similar helical parameters as in A3 and A4. The hydrogen-bonding scheme restricts the pitch in all the models to the narrow range of 10.0 to 11.5 Å. All these helices have large central cores whose radii increase proportionately with n. In this respect, A3 and A4 are suitable models for the structure of gramicidin A. In terms of their relative energies, antiparallel double-helices are marginally more stable than those with parallel strands. Our results indicate that the energy differences amongst the members in the antiparallel family are not significant and thus provide an explanation for the polymorphism reported for poly(γ-benzyl-LD-glutamate).

ERAP2 functional knockout in humans does not alter surface heavy chains or HLA-B27, inflammatory cytokines or endoplasmic reticulum stress markers

Introduction Single nucleotide polymorphisms in ERAP2 are strongly associated with ankylosing spondylitis (AS). One AS-associated single nucleotide polymorphism, rs2248374, causes a truncated ERAP2 protein that is degraded by nonsense-mediated decay. Approximately 25% of the populations of European ancestry are therefore natural ERAP2 knockouts. We investigated the effect of this associated variant on HLA class I allele presentation, surface heavy chains, endoplasmic reticulum (ER) stress markers and cytokine gene transcription in AS. Methods Patients with AS and healthy controls with either AA or GG homozygous status for rs2248374 were studied. Antibodies to CD14, CD19-ECD, HLA-A-B-C, Valpha7.2, CD161, anti-HC10 and anti-HLA-B27 were used to analyse peripheral blood mononuclear cells. Expression levels of ER stress markers (GRP78 and CHOP) and proinflammatory genes (tumour necrosis factor (TNF), IL6, IL17 and IL22) were assessed by qPCR. Results There was no significant difference in HLAclass I allele presentation or major histocompatibility class I heavy chains or ER stress markers GRP78 and CHOP or proinflammatory gene expression between genotypes for rs2248374 either between cases, between cases and controls, and between controls. Discussion Large differences were not seen in HLAB27 expression or cytokine levels between subjects with and without ERAP2 in AS cases and controls. This suggests that ERAP2 is more likely to influence AS risk through other mechanisms.

The curvature invariant for a class of homogeneous operators

For an operator T in the class B-n(), introduced by Cowen and Douglas, the simultaneous unitary equivalence class of the curvature and the covariant derivatives up to a certain order of the corresponding bundle E-T determine the unitary equivalence class of the operator T. In a subsequent paper, the authors ask if the simultaneous unitary equivalence class of the curvature and these covariant derivatives are necessary to determine the unitary equivalence class of the operator T is an element of B-n(). Here we show that some of the covariant derivatives are necessary. Our examples consist of homogeneous operators in B-n(). For homogeneous operators, the simultaneous unitary equivalence class of the curvature and all its covariant derivatives at any point w in the unit disc are determined from the simultaneous unitary equivalence class at 0. This shows that it is enough to calculate all the invariants and compare them at just one point, say 0. These calculations are then carried out in number of examples. One of our main results is that the curvature along with its covariant derivative of order (0, 1) at 0 determines the equivalence class of generic homogeneous Hermitian holomorphic vector bundles over the unit disc.

Asymptotic analysis of option pricing in a Markov modulated market

We address asymptotic analysis of option pricing in a regime switching market where the risk free interest rate, growth rate and the volatility of the stocks depend on a finite state Markov chain. We study two variations of the chain namely, when the chain is moving very fast compared to the underlying asset price and when it is moving very slow. Using quadratic hedging and asymptotic expansion, we derive corrections on the locally risk minimizing option price.

Direct probe of end-segment distribution in tethered polymer chains

End-tethered chains made of an adsorbed diblock copolymer of polystyrene (PS)-polyisoprene (PI) bearing an end-segment including a Ge atom are built by the Langmuir-Schaeffer technique. They are studied both in the dry state and in a good solvent for the PI chain using grazing incidence X-ray standing waves. The analysis of the signal provides a direct measurement of the end-segment distribution which is found to be singular and mostly localized to a plane in the dry case. In the good solvent case, end-segments are found to span the entire assembly and compare very well with results obtained by Kreer et al.

Effect of surfactants on the fluorescence spectra of water-soluble MEHPPV derivatives having grafted polyelectrolyte chains

Poly(2-methoxy-5-[2'-ethylhexyoxy]-1,4-phenylenevinylene) (MEHPPV) derivatives with polyacrylic acid (PAA) chains grafted onto their backbone were found to be water soluble, and they exhibited a dramatic increase in their fluorescence intensity in the presence of a variety of surfactants, even at concentrations far below their critical micelle concentrations (CMC). This increase was accompanied by a blue-shift in the emission maximum. These observations are rationalized based on the postulate that the backbone conformation of the conjugated polymer is modulated upon interaction of the surfactant molecules with the polyelectrolytic tethers, which in turn results in a significant depletion of intra-chain interchromophore interactions that are known to cause red-shifted emission bands with significantly lower emission yields.

Propargyloxycarbonyl as a protecting group for the side chains of serine, threonine and tyrosine

Propargyloxycarbonyl group is used as a protecting group for the hydroxyl groups of serine, threonine and tyrosine. The propargyloxycarbonyl derivatives of these hydroxy amino acids are stable to acidic and basic reagents commonly employed in peptide synthesis. The deprotection of the O-Poc derivatives using tetrathiomolybdate does not affect commonly used protecting groups such as N-Boc, N-Cbz, N-Fmoc, methyl and benzyl esters. The di-and tripeptides synthesized using O-Poc derivatives of serine, threonine and tyrosine are stable, isolable compounds and give the hydroxy peptides in good yields when treated with tetrathiomolybdate.

Experiments on homogeneous nucleation and physicochemical properties related to atmospheric new particle formation

Aerosol particles play a role in the earth ecosystem and affect human health. A significant pathway of producing aerosol particles in the atmosphere is new particle formation, where condensable vapours nucleate and these newly formed clusters grow by condensation and coagulation. However, this phenomenon is still not fully understood. This thesis brings an insight to new particle formation from an experimental point of view. Laboratory experiments were conducted both on the nucleation process and physicochemical properties related to new particle formation. Nucleation rate measurements are used to test nucleation theories. These theories, in turn, are used to predict nucleation rates in atmospheric conditions. However, the nucleation rate measurements have proven quite difficult to conduct, as different devices can yield nucleation rates with differences of several orders of magnitude for the same substances. In this thesis, work has been done to have a greater understanding in nucleation measurements, especially those conducted in a laminar flow diffusion chamber. Systematic studies of nucleation were also made for future verification of nucleation theories. Surface tensions and densities of substances related to atmospheric new particle formation were measured. Ternary sulphuric acid + ammonia + water is a proposed candidate to participate in atmospheric nucleation. Surface tensions of an alternative candidate to nucleate in boreal forest areas, sulphuric acid + dimethylamine + water, were also measured. Binary compounds, consisting of organic acids + water are possible candidates to participate in the early growth of freshly nucleated particles. All the measured surface tensions and densities were fitted with equations, thermodynamically consistent if possible, to be easily applied to atmospheric model calculations of nucleation and subsequent evolution of particle size.

On the homogeneous and heterogeneous nucleation of some organic compounds

The conversion of a metastable phase into a thermodynamically stable phase takes place via the formation of clusters. Clusters of different sizes are formed spontaneously within the metastable mother phase, but only those larger than a certain size, called the critical size, will end up growing into a new phase. There are two types of nucleation: homogeneous, where the clusters appear in a uniform phase, and heterogeneous, when pre-existing surfaces are available and clusters form on them. The nucleation of aerosol particles from gas-phase molecules is connected not only with inorganic compounds, but also with nonvolatile organic substances found in atmosphere. The question is which ones of the myriad of organic species have the right properties and are able to participate in nucleation phenomena. This thesis discusses both homogeneous and heterogeneous nucleation, having as theoretical tool the classical nucleation theory (CNT) based on thermodynamics. Different classes of organics are investigated. The members of the first class are four dicarboxylic acids (succinic, glutaric, malonic and adipic). They can be found in both the gas and particulate phases, and represent good candidates for the aerosol formation due to their low vapor pressure and solubility. Their influence on the nucleation process has not been largely investigated in the literature and it is not fully established. The accuracy of the CNT predictions for binary water-dicarboxylic acid systems depends significantly on the good knowledge of the thermophysical properties of the organics and their aqueous solutions. A large part of the thesis is dedicated to this issue. We have shown that homogeneous and heterogeneous nucleation of succinic, glutaric and malonic acids in combination with water is unlikely to happen in atmospheric conditions. However, it seems that adipic acid could participate in the nucleation process in conditions occurring in the upper troposphere. The second class of organics is represented by n-nonane and n-propanol. Their thermophysical properties are well established, and experiments on these substances have been performed. The experimental data of binary homogeneous and heterogeneous nucleation have been compared with the theoretical predictions. Although the n-nonane - n-propanol mixture is far from being ideal, CNT seems to behave fairly well, especially when calculating the cluster composition. In the case of heterogeneous nucleation, it has been found that better characterization of the substrate - liquid interaction by means of line tension and microscopic contact angle leads to a significant improvement of the CNT prediction. Unfortunately, this can not be achieved without well defined experimental data.