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.

Sequence diversity and differential expression of major phenylpropanoid-flavonoid biosynthetic genes among three mango varieties

Background: Mango fruits contain a broad spectrum of phenolic compounds which impart potential health benefits; their biosynthesis is catalysed by enzymes in the phenylpropanoid-flavonoid (PF) pathway. The aim of this study was to reveal the variability in genes involved in the PF pathway in three different mango varieties Mangifera indica L., a member of the family Anacardiaceae: Kensington Pride (KP), Irwin (IW) and Nam Doc Mai (NDM) and to determine associations with gene expression and mango flavonoid profiles. Results: A close evolutionary relationship between mango genes and those from the woody species poplar of the Salicaceae family (Populus trichocarpa) and grape of the Vitaceae family (Vitis vinifera), was revealed through phylogenetic analysis of PF pathway genes. We discovered 145 SNPs in total within coding sequences with an average frequency of one SNP every 316bp. Variety IW had the highest SNP frequency (one SNP every 258bp) while KP and NDM had similar frequencies (one SNP every 369bp and 360bp, respectively). The position in the PF pathway appeared to influence the extent of genetic diversity of the encoded enzymes. The entry point enzymes phenylalanine lyase (PAL), cinnamate 4-mono-oxygenase (C4H) and chalcone synthase (CHS) had low levels of SNP diversity in their coding sequences, whereas anthocyanidin reductase (ANR) showed the highest SNP frequency followed by flavonoid 3'-hydroxylase (F3'H). Quantitative PCR revealed characteristic patterns of gene expression that differed between mango peel and flesh, and between varieties. Conclusions: The combination of mango expressed sequence tags and availability of well-established reference PF biosynthetic genes from other plant species allowed the identification of coding sequences of genes that may lead to the formation of important flavonoid compounds in mango fruits and facilitated characterisation of single nucleotide polymorphisms between varieties. We discovered an association between the extent of sequence variation and position in the pathway for up-stream genes. The high expression of PAL, C4H and CHS genes in mango peel compared to flesh is associated with high amounts of total phenolic contents in peels, which suggest that these genes have an influence on total flavonoid levels in mango fruit peel and flesh. In addition, the particularly high expression levels of ANR in KP and NDM peels compared to IW peel and the significant accumulation of its product epicatechin gallate (ECG) in those extracts reflects the rate-limiting role of ANR on ECG biosynthesis in mango. © 2015 Hoang et al.

The shifting influence of drought and heat stress for crops in Northeast Australia

Characterization of drought environment types (ETs) has proven useful for breeding crops for drought-prone regions. Here we consider how changes in climate and atmospheric carbon dioxide (CO2) concentrations will affect drought ET frequencies in sorghum and wheat systems of Northeast Australia. We also modify APSIM (the Agricultural Production Systems Simulator) to incorporate extreme heat effects on grain number and weight, and then evaluate changes in the occurrence of heat-induced yield losses of more than 10, as well as the co-occurrence of drought and heat. More than six million simulations spanning representative locations, soil types, management systems, and 33 climate projections led to three key findings. First, the projected frequency of drought decreased slightly for most climate projections for both sorghum and wheat, but for different reasons. In sorghum, warming exacerbated drought stresses by raising the atmospheric vapor pressure deficit and reducing transpiration efficiency (TE), but an increase in TE due to elevated CO2 more than offset these effects. In wheat, warming reduced drought stress during spring by hastening development through winter and reducing exposure to terminal drought. Elevated CO2 increased TE but also raised radiation use efficiency and overall growth rates and water use, thereby offsetting much of the drought reduction from warming. Second, adding explicit effects of heat on grain number and grain size often switched projected yield impacts from positive to negative. Finally, although average yield losses associated with drought will remain generally higher than for heat stress for the next half century, the relative importance of heat is steadily growing. This trend, as well as the likely high degree of genetic variability in heat tolerance, suggests that more emphasis on heat tolerance is warranted in breeding programs. At the same time, work on drought tolerance should continue with an emphasis on drought that co-occurs with extreme heat. This article is protected by copyright. All rights reserved.

Isotope shifts and hyperfine structure in the 555.8-nm S-1(0) -> P-3(1) line of Yb

We apply our technique of using a Rb-stabilized ring-cavity resonator to measure the frequencies of various spectral components in the 555.8-nm 1S0-->3P1 line of Yb. We determine the isotope shifts with 60 kHz precision, which is an order-of-magnitude improvement over the best previous measurement on this line. There are two overlapping transitions, 171Yb(1/2-->3/2) and 173Yb(5/2-->3/2), which we resolve by applying a magnetic field. We thus obtain the hyperfine constants in the 3P1 state of the odd isotopes with a significantly improved precision. Knowledge of isotope shifts and hyperfine structure should prove useful for high-precision calculations in Yb necessary to interpret ongoing experiments testing parity and time-reversal symmetry violation in the laws of physics.

Viscoelastic Properties of Nanocrystalline Films of Semiconducting Chalcogenides at Liquid/Liquid Interface

The interfacial shear rheological properties of a continuous single-crystalline film of CuS and a 3D particulate gel of CdS nanoparticles (3−5 nm in diameter) formed at toluene−water interfaces have been studied. The ultrathin films (50 nm in thickness) are formed in situ in the shear cell through a reaction at the toluene−water interface between a metal−organic compound in the organic layer and an appropriate reagent for sulfidation in the aqueous layer. Linear viscoelastic spectra of the nanofilms reveal solid-like rheological behavior with the storage modulus higher than the loss modulus over the range of angular frequencies probed. Large strain amplitude sweep measurements on the CdS nanofilms formed at different reactant concentrations suggest that they form a weakly flocculated gel. Under steady shear, the films exhibit a yield stress, followed by a steady shear thinning at high shear rates. The viscoelastic and flow behavior of these films that are in common with those of many 3D “soft” materials like gels, foams, and concentrated colloidal suspensions can be described by the “soft” glassy rheology model.

The shifting influence of drought and heat stress for crops in Northeast Australia

Characterization of drought environment types (ETs) has proven useful for breeding crops for drought-prone regions. Here we consider how changes in climate and atmospheric carbon dioxide (CO2) concentrations will affect drought ET frequencies in sorghum and wheat systems of Northeast Australia. We also modify APSIM (the Agricultural Production Systems Simulator) to incorporate extreme heat effects on grain number and weight, and then evaluate changes in the occurrence of heat-induced yield losses of more than 10%, as well as the co-occurrence of drought and heat. More than six million simulations spanning representative locations, soil types, management systems, and 33 climate projections led to three key findings. First, the projected frequency of drought decreased slightly for most climate projections for both sorghum and wheat, but for different reasons. In sorghum, warming exacerbated drought stresses by raising the atmospheric vapor pressure deficit and reducing transpiration efficiency (TE), but an increase in TE due to elevated CO2 more than offset these effects. In wheat, warming reduced drought stress during spring by hastening development through winter and reducing exposure to terminal drought. Elevated CO2 increased TE but also raised radiation use efficiency and overall growth rates and water use, thereby offsetting much of the drought reduction from warming. Second, adding explicit effects of heat on grain number and grain size often switched projected yield impacts from positive to negative. Finally, although average yield losses associated with drought will remain generally higher than for heat stress for the next half century, the relative importance of heat is steadily growing. This trend, as well as the likely high degree of genetic variability in heat tolerance, suggests that more emphasis on heat tolerance is warranted in breeding programs. At the same time, work on drought tolerance should continue with an emphasis on drought that co-occurs with extreme heat. This article is protected by copyright. All rights reserved.

Characterization of maize growing environments in eastern and southern Africa using the APSIM model

Maize grown in eastern and southern Africa experiences random occurrences of drought. This uncertainty creates difficulty in developing superior varieties and their agronomy. Characterisation of drought types and their frequencies could help in better defining selection environments for improving resistance to drought. We used the well tested APSIM maize model to characterise major drought stress patterns and their frequencies across six countries of the region including Ethiopia, Kenya, Tanzania, Malawi, Mozambique and Zimbabwe. The database thus generated covered 35 sites, 17 to 86 years of daily climate records, 3 varieties and 3 planting densities from a total of 11,174 simulations. The analysis identified four major drought environment types including those characterised by low-stress which occurred in 42% of the years, mid-season drought occurring in 15% of the years, late-terminal stress which occurred in 22% of the years and early-terminal drought occurring in 21% of the years. These frequencies varied in relation to sites, genotypes and management. The simulations showed that early terminal stress could result in a yield reduction of 70% compared with low-stress environmental types. The study presents the importance of environmental characterization in contributing to maize improvement in eastern and southern Africa.

Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

Molecular vibrations and conformation of primary dithiocarbamate ester. Infrared and NMR studies on S-methyl dithiocarbamate

The i.r. spectra of a primary dithiocarbamate ester namely, S-methyl dithiocarbamate (SMDTC) and its N-dideuterated compound have been measured between 4000 and 30 cm−1. Spectra in solution and at liquid nitrogen temperature have also been obtained. Assignment of all the fundamentals has been proposed and supported from a full normal coordinate analysis. The band assignments for SMDTC have been compared with those of related molecules and the characteristic bands of primary thioamides are derived. Conformational flexibility of SMDTC has been examined by i.r. and proton NMR spectroscopy. The hindered rotation around the C---N bond has been studied by a complete line shape analysis. The magnitude of ---NH2 and ---CH3 torsional barriers is also estimated from vibrational frequencies and force constants.

Phosphine resistance in India is characterised by a dihydrolipoamide dehydrogenase variant that is otherwise unobserved in eukaryotes

Phosphine (PH3) fumigation is the primary method worldwide for controlling insect pests of stored commodities. Over-reliance on phosphine, however, has led to the emergence of strong resistance. Detailed genetic studies previously identified two loci, rph1 and rph2, that interact synergistically to create a strong resistance phenotype. We compared the genetics of phosphine resistance in strains of Rhyzopertha dominica and Tribolium castaneum from India and Australia, countries having similar pest species but widely differing in pest management practices. Sequencing analysis of the rph2 locus, dihydrolipoamide dehydrogenase (dld), identified two structurally equivalent variants, Proline49>Serine (P49S) in one R. dominica strain and P45S in three strains of T. castaneum from India. These variants of the DLD protein likely affect FAD cofactor interaction with the enzyme. A survey of insects from storage facilities across southern India revealed that the P45/49S variant is distributed throughout the region at very high frequencies, in up to 94% of R. dominica and 97% of T. castaneum in the state of Tamil Nadu. The abundance of the P45/49S variant in insect populations contrasted sharply with the evolutionary record in which the variant was absent from eukaryotic DLD sequences. This suggests that the variant is unlikely to provide a strong selective advantage in the absence of phosphine fumigation.

Temperature dependence of 35Cl NQR in 3,4-Dichlorophenol

NQR frequencies in 3,4-dichlorophenol are investigated in the temperature range 77 K to room temperature. Two resonances have been observed throughout the temperature range, corresponding to the two chemically inequivalent chlorine sites. Using Bayer's theory and Brown's method torsional frequencies and their temperature dependence in this range are estimated.

Nonlinear mode coupling of surface acoustic waves on an isotropic solid

The nonlinear mode coupling between two co-directional quasi-harmonic Rayleigh surface waves on an isotropic solid is analysed using the method of multiple scales. This procedure yields a system of six semi-linear hyperbolic partial differential equations with the same principal part governing the slow variations in the (complex) amplitudes of the two fundamental, the two second harmonic and the two combination frequency waves at the second stage of the perturbation expansion. A numerical solution of these equations for excitation by monochromatic signals at two arbitrary frequencies, indicates that there is a continuous transfer of energy back and forth among the fundamental, second harmonic and combination frequency waves due to mode coupling. The mode coupling tends to be more pronounced as the frequencies of the interacting waves approach each other.

A problem in micropolar elasticity

In this paper a three-dimensional analysis for statics and dynamics of a class of simply supported rectangular plates made up of micropolar elastic material is presented. The solution is in the form of series, in which each term is explicitly determined. For free vibrations, the frequencies are obtained by the solution of a closed form characteristic equation.

Temperature-dependence of 35Cl N.Q.R. in sodium chloroacetate

N.q.r. in sodium chloroacetate has been investigated at temperatures from 77 K to room temperature (c.300 K). A single line has been observed throughout this temperature range. Torsional frequencies of the molecule have been calculated in the above temperature range from Bayer's theory. Also the temperature coefficient of the torsional frequencies has been calculated by Brown's method.

Vibration of generally orthotropic skew plates

Vibration problem of generally orthotropic plates with particular attention to plates of skew geometry is studied. The formulation is based on orthotropic plate theory with arbitrary orientation of the principal axes of orthotropy. The boundary conditions considered are combinations of simply supported, clamped, and free-edge conditions. Approximate solution for frequencies and modes is obtained by the Ritz method using products of appropriate beam characteristic functions as admissible functions. The variation of frequencies and modes with orientation of the axes of orthotropy is examined for different skew angles and boundary conditions. Features such as "crossings" and "quasi-degeneracies" of the frequency curves are found to occur with variation of the orientation of the axes of orthotropy for a given geometry of the skew plate. It is also found that for each combination of skew angle and side ratio, a particular orientation of the axes gives the highest value for the fundamental frequency of the plate.