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.

An uncommon nucleotide conformation shown by molecular structure of deoxyuridine-5-phosphate and nucleic acid stereochemistry

CRYSTAL structure determinations of nucleic acid fragments have shown that several of the conformational features found in the monomeric building blocks are also manifested at the nucleic acid level. Stereochemical variations between thymine and uracil nucleotides are therefore of interest as they can provide a structural basis for some of the differences between the conformations of DNA and RNA. X-ray studies have so far not shown any major dissimilarities between these two nucleotide species although the sugar ring of deoxyribonucleotides is found to possess greater flexibility than that in ribonucleotides. We report here the molecular structure of deoxyuridine-5'-phosphate (dUMP-5') which is not a common monomer unit of DNAs as it is replaced by its 5-methyl analogue deoxythymidine-5'-phosphate (dTMP-5'). The investigation was undertaken to help determine whether or not this implied a fundamental difference between the geometries of these two molecules.

Neutron-diffraction study of structure and conformation of nucleotide uridine-5 phosphate (disodium salt)

Acta Crystallographica Section A: Foundations of Crystallography covers theoretical and fundamental aspects of the structure of matter. The journal is the prime forum for research in diffraction physics and the theory of crystallographic structure determination by diffraction methods using X-rays, neutrons and electrons. The structures include periodic and aperiodic crystals, and non-periodic disordered materials, and the corresponding Bragg, satellite and diffuse scattering, thermal motion and symmetry aspects. Spatial resolutions range from the subatomic domain in charge-density studies to nanodimensional imperfections such as dislocations and twin walls. The chemistry encompasses metals, alloys, and inorganic, organic and biological materials. Structure prediction and properties such as the theory of phase transformations are also covered.

Antiferromagnetic spin-wave states in a linear chain with nearest- and next-nearest-neighbour interactions

The spectrum of short-closed chains up to N=12 are studied by exact diagonalization to obtain the spin-wave spectrum of the Hamiltonian H=2J Sigma i=1Nsi.si+1+2J alpha Sigma i=1Nsi.si+2, -1.0chains indicate that S=0 bound magnon states lie partly below the spin-wave states with S=1 only for alpha >or=0.3 and alpha

Rate of ribonucleic acid chain growth in Mycobacterium tuberculosis H37Rv.

Two methods were employed to measure the rate of ribonucleic acid (RNA) chain growth in vivo in Mycobacterium tuberculosis H37Rv cultures growing in Sauton medium at 37 degrees C, with a generation time of 10 h. In the first, the bacteria were allowed to assimilate [3H]uracil or [3H]guanine into their RNA for short time periods. The RNA was then extracted and hydrolyzed with alkali, and the radioactivity in the resulting nucleotides and nucleosides was measured. The data obtained by this method allowed the calculation of the individual nucleotide step times during the growth of RNA chains, from which the average rate of RNA chain elongation was estimated to be about 4 nucleotides per s. The second method employed the antibiotic rifampin, which specifically inhibits the initiation of RNA synthesis without interfering with the elongation and completion of nascent RNA chains. Usint this method, the transcription time of the 16S, 23S, and 5S ribosomal RNA genes was estimated to be 7.6 min, which corresponds to a ribosomal RNA chain growth rate of 10 nucleotides per s.

Functional studies of purified transmembrane proteases, omptins, of Yersinia pestis and Salmonella enterica.

Surface proteolysis is important in migration of cells through tissue barriers. In the case of prokaryotes, surface proteolysis has been associated with invasiveness of pathogenic bacteria from the primary infection site into circulation and secondary infection sites in the host. This study addressed surface proteases of two important bacterial pathogens, Yersinia pestis which is the causative agent of the lethal systemic zoonosis, plague, and Salmonella enterica serovar Typhimurium which is an oral-faecal pathogen that annually causes millions of cases of gastoenteritis that may develop to septicaemia. Both bacterial species express an ortholog of the omptin family of transmembrane β-barrel, outer membrane proteases/adhesins. This thesis work addressed the functions of isolated plasminogen activator Pla of Y. pestis and the PgtE omptin of S. enterica. Pla and PgtE were isolated as His6-fusion proteins in denaturing conditions from recombinant Escherichia coli and activated by adding lipopolysaccharide (LPS). The structural features in LPS that enhance plasminogen activation by His6-Pla were determined, and it was found that the lack of O-specifi c chain, the presence of outer core oligosaccharide, the presence of phosphates in lipid A, as well as a low level of acylation in lipid A influence the enhancement of Pla activity by LPS. A conserved lipid A phosphate binding motif in Pla and PgtE was found important for the enhancement of enzymatic activity by LPS. The results help to explain the biological signifi cance of the genetic loss of the O-specifi c chain biosynthesis in Y. pestis as well as the variations in LPS structure upon entry of Y. pestis into the human host. Expression of Pla in Y. pestis is associated with adhesiveness to lamin of basement membranes. Here, isolated and LPS-activated His6-Pla was coated onto fluorescent microparticles. The coating conferred specifi c adhesiveness of the particles to laminin and reconstituted basement membrane, thus confi rming the intrinsic adhesive characteristics of the Pla protein. The adhesiveness is thought to direct plasmin proteolysis at tissue barriers, thus increasing tissue damage and bacterial spread. Gelatinase activity has not been previously reported in enteric bacteria. Expression of PgtE in S. enterica was associated with cleavage of porcine skin gelatin, denaturated human type I collagen, as well as DQ-gelatin. Purifi ed His6-PgtE also degraded porcine skin gelatin and human type I gelatin but did not react with DQ-gelatin, indicating that minor differences are seen in proteolysis by isolated and cell-bound PgtE. Pla was less effective in gelatin degradation. The novel gelatinase activity in S. enterica is likely to enhance bacterial dissemination during infection.

LFA-1 Integrin beta2 Chain Phosphorylation Regulates Protein Interactions and Mediates Signals in T Cells

Integrins are heterodimeric transmembrane adhesion receptors composed of alpha- and beta-subunits and they are vital for the function of multicellular organisms. Integrin-mediated adhesion is a complex process involving both affinity regulation and coupling to the actin cytoskeleton. Integrins also function as bidirectional signaling devices, regulating cell adhesion and migration after inside-out signaling, but also signal into the cell to regulate growth, differentiation and apoptosis after ligand binding. The LFA-1 integrin is exclusively expressed in leukocytes and is of fundamental importance for the function of the immune system. The LFA-1 integrins have short intracellular tails, which are devoid of catalytic activity. These cytoplasmic domains are important for integrin regulation and both the alpha and beta chain become phosphorylated. The alpha chain is constitutively phosphorylated, but the beta chain becomes phosphorylated on serine and functionally important threonine residues only after cell activation. The cytoplasmic tails of LFA-1 bind to many cytoskeletal and signaling proteins regulating numerous cell functions. However, the molecular mechanisms behind these interactions have been poorly understood. Phosphorylation of the cytoplasmic tails of the LFA-1 integrin could provide a mechanism to regulate integrin-mediated cytoskeletal interactions and take part in T cell signaling. In this study, the effects of phosphorylation of LFA-1 integrin cytoplasmic tails on different cellular functions were examined. Site-specific phosphorylation of both the alpha- and beta-chains of the LFA-1 was shown to have a role in the regulation of the LFA-1 integrin.Alpha-chain Ser1140 is needed for integrin conformational changes after chemokine- or integrin ligand-induced activation or after activation induced by active Rap1, whereas beta-chain binds to 14-3-3 proteins through the phosphorylated Thr758 and mediates cytoskeletal reorganization. Thr758 phosphorylation also acts as a molecular switch to inhibit filamin binding and allows 14-3-3 protein binding to integrin cytoplasmic domain, and it was also shown to lead to T cell adhesion, Rac-1/Cdc42 activation and expression of the T cell activation marker CD69, indicating a signaling function for Thr758 phosphorylation in T cells. Thus, phosphorylation of the cytoplasmic tails of LFA-1 plays an important role in different functions of the LFA-1 integrin in T cells. It is of vital importance to study the mechanisms and components of integrin regulation since leukocyte adhesion is involved in many functions of the immune system and defects in the regulation of LFA-1 contributes to auto-immune diseases and fundamental defects in the immune system.

Improvement in the voltage grading (axial and radial) of the generator column of a Van de Graaf generator by the use of a resistor chain

It is well known that the use of a series of resistors, connected between the equipotential rings of a Van de Graaff generator, improves the axial voltage grading of the generator. The work reported in this paper shows how the resistor chain also improves the radial voltage gradient. The electrolytic field mapping technique was adopted in the present work.

Entanglement production due to quench dynamics of an anisotropic XY chain in a transverse field

We compute concurrence and negativity as measures of two-spin entanglement generated by a power-law quench (characterized by a rate tau(-1) and an exponent alpha) which takes an anisotropic XY chain in a transverse field through a quantum critical point (QCP). We show that only spins separated by an even number of lattice spacings get entangled in such a process. Moreover, there is a critical rate of quench, tau(-1)(c), above which no two-spin entanglement is generated; the entire entanglement is multipartite. The ratio of the entanglements between consecutive even neighbors can be tuned by changing the quench rate. We also show that for large tau, the concurrence (negativity) scales as root alpha/tau(alpha/tau), and we relate this scaling behavior to defect production by the quench through a QCP.

Engineering the pentose phosphate pathway of Saccharomyces cerevisiae for production of ethanol and xylitol

The baker s yeast Saccharomyces cerevisiae has a long tradition in alcohol production from D-glucose of e.g. starch. However, without genetic modifications it is unable to utilise the 5-carbon sugars D-xylose and L arabinose present in plant biomass. In this study, one key metabolic step of the catabolic D-xylose pathway in recombinant D-xylose-utilising S. cerevisiae strains was studied. This step, carried out by xylulokinase (XK), was shown to be rate-limiting, because overexpression of the xylulokinase-encoding gene XKS1 increased both the specific ethanol production rate and the yield from D xylose. In addition, less of the unwanted side product xylitol was produced. Recombinant D-xylose-utilizing S. cerevisiae strains have been constructed by expressing the genes coding for the first two enzymes of the pathway, D-xylose reductase (XR) and xylitol dehydrogenase (XDH) from the D-xylose-utilising yeast Pichia stipitis. In this study, the ability of endogenous genes of S. cerevisiae to enable D-xylose utilisation was evaluated. Overexpression of the GRE3 gene coding for an unspecific aldose reductase and the ScXYL2 gene coding for a xylitol dehydrogenase homologue enabled growth on D-xylose in aerobic conditions. However, the strain with GRE3 and ScXYL2 had a lower growth rate and accumulated more xylitol compared to the strain with the corresponding enzymes from P. stipitis. Use of the strictly NADPH-dependent Gre3p instead of the P. stipitis XR able to utilise both NADH and NADPH leads to a more severe redox imbalance. In a S. cerevisiae strain not engineered for D-xylose utilisation the presence of D-xylose increased xylitol dehydrogenase activity and the expression of the genes SOR1 or SOR2 coding for sorbitol dehydrogenase. Thus, D-xylose utilisation by S. cerevisiae with activities encoded by ScXYL2 or possibly SOR1 or SOR2, and GRE3 is feasible, but requires efficient redox balance engineering. Compared to D-xylose, D-glucose is a cheap and readily available substrate and thus an attractive alternative for xylitol manufacture. In this study, the pentose phosphate pathway (PPP) of S. cerevisiae was engineered for production of xylitol from D-glucose. Xylitol was formed from D-xylulose 5-phosphate in strains lacking transketolase activity and expressing the gene coding for XDH from P. stipitis. In addition to xylitol, ribitol, D-ribose and D-ribulose were also formed. Deletion of the xylulokinase-encoding gene increased xylitol production, whereas the expression of DOG1 coding for sugar phosphate phosphatase increased ribitol, D-ribose and D-ribulose production. Strains lacking phosphoglucose isomerase (Pgi1p) activity were shown to produce 5 carbon compounds through PPP when DOG1 was overexpressed. Expression of genes encoding glyceraldehyde 3-phosphate dehydrogenase of Bacillus subtilis, GapB, or NAD-dependent glutamate dehydrogenase Gdh2p of S. cerevisiae, altered the cellular redox balance and enhanced growth of pgi1 strains on D glucose, but co-expression with DOG1 reduced growth on higher D-glucose concentrations. Strains lacking both transketolase and phosphoglucose isomerase activities tolerated only low D-glucose concentrations, but the yield of 5-carbon sugars and sugar alcohols on D-glucose was about 50% (w/w).

Synthesis, characterization, thermal degradation, and comparative chain dynamics studies of weak-link polysulfide polymers

Synthesis, spectroscopic and thermal characterization of two new classes of polysulfide polymers: poly[1(phenoxymethyl) ethylene polysulfide] (PPMEP), and poly [1-(phenoxy) ethylene polysulfide] (PPEP) is presented. The direct pyrolysis mass spectrometry (DP-MS) technique, used to study the thermal degradation behavior of these polysulfide polymers, indicated that the polymers underwent degradation through the weak-links scission. The thermal stability of the polysulfide polymers decreased as the ``rank'' (number of sulfur atoms in the polysulfide linkage; n=1, 2, 4) increased. The main-chain flexibility of these polysulfide polymers in terms of their C-13 NMR spinlattice relaxation time (T-1) measurements on the backbone methine (-CH-) and methylene (-CH2-) carbons are reported here for the first time. A comparative study of the solution chain dynamics indicated that it increased as ``rank'' of the polysulfide linkages decreased as well as by introducing side chain spacers such as, ether (-O-) or methyleneoxy (-CH2O-) groups.

Structure of disodium deoxyuridine 5'-phosphate pentahydrate

Disodium deoxyuridine 5'-nhosDhate pentahvdrate, Na2(C9H l INEOsP). 5 H20, Call 11N208 P2-. 2Na +. 5 H20, crystallizes in the monoclinic space group P2: with a = 7.250 (4), b = 35.45 (2), c = 7.132 (4)/~, fl = 102.2 (4) °, Z = 4. The Cu Ka intensity data were collected photographically and estimated visually. The structure was obtained by the minimum-function method and difference syntheses and refined to an R of 0.089. In both molecules the uracil base has an anti conformation (2cN = 57.1 and 59.9 °) with respect to the sugar. The deoxyribose moiety of molecule B shows a typical C(l')-exo puckering, with C(I') displaced by 0.52 /k from the best plane. The furanose ring conformation of molecule A can be described as C(2')-endo,C(l')-exo. Both the molecules have an unusual trans-gauche conformation about the exocyclic C(4')-C(5') bond with (~0oo = 171.1, 172.2°; ~0oc = -64.7, -65.9°).

Conformational studies of the triphenylphosphazenyl side chain in cyclophosphazenes. II. Crystal and molecular structure of 2,4,4,6,6,8,8-heptachloro-2-triphenylphosphazenylcyclotetra(phosphazene), N4P4Cl7(NPPh3)

Abstract is not available.