Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Inference on microsatellite mutation processes in the invasive mite, Varroa destructor, using reversible jump Markov chain Monte Carlo

Varroa destructor is a parasitic mite of the Eastern honeybee Apis cerana. Fifty years ago, two distinct evolutionary lineages (Korean and Japanese) invaded the Western honeybee Apis mellifera. This haplo-diploid parasite species reproduces mainly through brother sister matings, a system which largely favors the fixation of new mutations. In a worldwide sample of 225 individuals from 21 locations collected on Western honeybees and analyzed at 19 microsatellite loci, a series of de novo mutations was observed. Using historical data concerning the invasion, this original biological system has been exploited to compare three mutation models with allele size constraints for microsatellite markers: stepwise (SMM) and generalized (GSM) mutation models, and a model with mutation rate increasing exponentially with microsatellite length (ESM). Posterior probabilities of the three models have been estimated for each locus individually using reversible jump Markov Chain Monte Carlo. The relative support of each model varies widely among loci, but the GSM is the only model that always receives at least 9% support, whatever the locus. The analysis also provides robust estimates of mutation parameters for each locus and of the divergence time of the two invasive lineages (67,000 generations with a 90% credibility interval of 35,000-174,000). With an average of 10 generations per year, this divergence time fits with the last post-glacial Korea Japan land separation. (c) 2005 Elsevier Inc. All rights reserved.

Bayesian analysis of correlated evolution of discrete characters by reversible-jump Markov chain Monte Carlo

We describe a Bayesian method for investigating correlated evolution of discrete binary traits on phylogenetic trees. The method fits a continuous-time Markov model to a pair of traits, seeking the best fitting models that describe their joint evolution on a phylogeny. We employ the methodology of reversible-jump ( RJ) Markov chain Monte Carlo to search among the large number of possible models, some of which conform to independent evolution of the two traits, others to correlated evolution. The RJ Markov chain visits these models in proportion to their posterior probabilities, thereby directly estimating the support for the hypothesis of correlated evolution. In addition, the RJ Markov chain simultaneously estimates the posterior distributions of the rate parameters of the model of trait evolution. These posterior distributions can be used to test among alternative evolutionary scenarios to explain the observed data. All results are integrated over a sample of phylogenetic trees to account for phylogenetic uncertainty. We implement the method in a program called RJ Discrete and illustrate it by analyzing the question of whether mating system and advertisement of estrus by females have coevolved in the Old World monkeys and great apes.

An analysis of the impact of R&D on productivity using Bayesian model averaging with a reversible jump algorithm.

A Bayesian Model Averaging approach to the estimation of lag structures is introduced, and applied to assess the impact of R&D on agricultural productivity in the US from 1889 to 1990. Lag and structural break coefficients are estimated using a reversible jump algorithm that traverses the model space. In addition to producing estimates and standard deviations for the coe¢ cients, the probability that a given lag (or break) enters the model is estimated. The approach is extended to select models populated with Gamma distributed lags of di¤erent frequencies. Results are consistent with the hypothesis that R&D positively drives productivity. Gamma lags are found to retain their usefulness in imposing a plausible structure on lag coe¢ cients, and their role is enhanced through the use of model averaging.

Reversible phase variation in the phnE gene, which is required for phosphonate metabolism in Escherichia coli K-12

It is known that Escherichia coli K-12 is cryptic (Phn(-)) for utilization of methyl phosphonate (MePn) and that Phn(+) variants can be selected for growth on MePn as the sole P source. Variants arise from deletion via a possible slip strand mechanism of one of three direct 8-bp repeat sequences in phnE, which restores function to a component of a putative ABC type transporter. Here we show that Phn(+) variants are present at the surprisingly high frequency of >10(-2) in K-12 strains. Amplified-fragment length polymorphism analysis was used to monitor instability in phnE in various strains growing under different conditions. This revealed that, once selection for growth on MePn is removed, Phn(+) revertants reappear and accumulate at high levels through reinsertion of the 8-bp repeat element sequence. It appears that, in K-12, phnE contains a high-frequency reversible gene switch, producing phase variation which either allows ("on" form) or blocks ("off" form) MePn utilization. The switch can also block usage of other metabolizable alkyl phosphonates, including the naturally occurring 2-aminoethylphosphonate. All K-12 strains, obtained from collections, appear in the "off" form even when bearing mutations in mutS, mutD, or dnaQ which are known to enhance slip strand events between repetitive sequences. The ability to inactivate the phnE gene appears to be unique to K-12 strains since the B strain is naturally Phn(+) and lacks the inactivating 8-bp insertion in phnE, as do important pathogenic strains for which genome sequences are known and also strains isolated recently from environmental sources.

Reversible, nondegradative conversion of crystalline aromatic poly(ether ketone)s into organo-soluble poly(ether dithioketal)s

Crystalline aromatic poly(ether ketone)s Such as PEEK and PEK may be cleanly and reversibly derivatized by dithioketalization of the carbonyl groups With 1,2-ethanedithiol or 1,3-propanedithiol under strong acid conditions. The resulting 1,3-dithiolane and 1,3-dithiane polymers are hydrolytically stable, amorphous, and readily soluble in organic solvents such as chloroform and THF and are thus (unlike their parent polymers) easily characterized by gel permeation chromatography (GPC). GPC analysis of a range of derivatized PEEK samples using light-scattering detection revealed, in some instances, a bimodal molecular weight distribution with a small but potentially significant (and previously undetected) very high-molecular-weight fraction.

A pressure-jump time-resolved x-ray diffraction study of cubic-cubic transition kinetics in monoolein

In the past two decades, the geometric pathways involved in the transformations between inverse bicontinuous cubic phases in amphiphilic systems have been extensively theoretically modeled. However, little experimental data exists on the cubic-cubic transformation in pure lipid systems. We have used pressure-jump time-resolved X-ray diffraction to investigate the transition between the gyroid Q(II)(G) and double-diamond Q(II)(D) phases in mixtures of 1-monoolein in 30 wt% water. We find for this system that the cubic-cubic transition occurs without any detectable intermediate structures. In addition, we have determined the kinetics of the transition, in both the forward and reverse directions, as a function of pressure-jump amplitude, temperature, and water content. A recently developed model allows (at least in principle) the calculation of the activation energy for lipid phase transitions from such data. The analysis is applicable only if kinetic reproducibility is achieved, at least within one sample, and achievement of such kinetic reproducibility is shown here, by carrying out prolonged pressure-cycling. The rate of transformation shows clear and consistent trends with pressure-jump amplitude, temperature, and water content, all of which are shown to be in agreement with the effect of the shift in the position of the cubic-cubic phase boundary following a change in the thermodynamic parameters.

Tin and tin-titanium as catalyst components for reversible hydrogen storage of sodium aluminium hydride

This paper is concerned with the effects of adding tin and/or titanium dopant to sodium aluminium hydride for both dehydrogenation and re-hydrogenation reactions during their reversible storage of molecular hydrogen. Temperature programmed decomposition (TPD) measurements show that the dehydrogenation kinetics of NaAlH4 are significantly enhanced upon doping the material with 2 mol% of tributyltin hydride, Sn(Bu)(3)H but the tin catalyst dopant is shown to be inferior than titanium. On the other hand, in this preliminary work, a significant synergetic catalytic effect is clearly revealed in material co-doped with both titanium and tin catalysts which shows the highest reversible rates of dehydrogenation and re-hydrogenation (after their hydrogen depletion). The re-hydrogenation rates of depleted Sn/Ti/NaAlH4 evaluated at both 9.5 and 140 bars hydrogen are also found to be favourable compared to the Ti/NaAlH4, which clearly suggest the importance of the catalyst choice. Basing on these results some mechanistic insights for the catalytic reversible dehydrogenation and re-hydrogenation processes of Sn/Ti/NaAlH4 are therefore made. (C) 2006 Elsevier Ltd. All rights reserved.

Pressure-jump X-ray studies of liquid crystal transitions in lipids

In this paper, we give an overview of our studies by static and time-resolved X-ray diffraction of inverse cubic phases and phase transitions in lipids. In 1, we briefly discuss the lyotropic phase behaviour of lipids, focusing attention on non-lamellar structures, and their geometric/topological relationship to fusion processes in lipid membranes. Possible pathways for transitions between different cubic phases are also outlined. In 2, we discuss the effects of hydrostatic pressure on lipid membranes and lipid phase transitions, and describe how the parameters required to predict the pressure dependence of lipid phase transition temperatures can be conveniently measured. We review some earlier results of inverse bicontinuous cubic phases from our laboratory, showing effects such as pressure-induced formation and swelling. In 3, we describe the technique of pressure-jump synchrotron X-ray diffraction. We present results that have been obtained from the lipid system 1:2 dilauroylphosphatidylcholine/lauric acid for cubic-inverse hexagonal, cubic-cubic and lamellar-cubic transitions. The rate of transition was found to increase with the amplitude of the pressure-jump and with increasing temperature. Evidence for intermediate structures occurring transiently during the transitions was also obtained. In 4, we describe an IDL-based 'AXCESS' software package being developed in our laboratory to permit batch processing and analysis of the large X-ray datasets produced by pressure-jump synchrotron experiments. In 5, we present some recent results on the fluid lamellar-Pn3m cubic phase transition of the single-chain lipid 1-monoelaidin, which we have studied both by pressure-jump and temperature-jump X-ray diffraction. Finally, in 6, we give a few indicators of future directions of this research. We anticipate that the most useful technical advance will be the development of pressure-jump apparatus on the microsecond time-scale, which will involve the use of a stack of piezoelectric pressure actuators. The pressure-jump technique is not restricted to lipid phase transitions, but can be used to study a wide range of soft matter transitions, ranging from protein unfolding and DNA unwinding and transitions, to phase transitions in thermotropic liquid crystals, surfactants and block copolymers.

Reversible ageing effects in cryogenically cooled infrared filter radiometers

In this paper we report the observation of drifts in the responsivity of cryogenically cooled InSb detector-based infrared filter radiometers which have very strong wavelength dependence. These drifts can result in the increase or decrease of the response of the filter radiometers by over 5%. The origin of these variations was investigated and was shown to arise due to a thin film of ice formed on the multi-layer bandpass filter used to define the spectral response of the filter radiometer. The thin layer of ice interacts with the characteristics of the filter (which itself consists of a number of thin layers) and modifies the filter spectral transmission thus modifying the response of the filter radiometer of which the filter is part of. These observations are particularly relevant to space instruments which use infrared filter radiometers for earth observation. Debris from the spacecraft engines is known to accumulate on cold surfaces of instruments carried on board. The deposition of this debris on cold filters can modify the spectral response of the instruments, which use these filters to define a spectral response. Crown Copyright (c) 2004 Published by Elsevier B.V. All rights reserved.

Reversible sphere-to-lamellar wetting transition at the interface of a diblock copolymer system

We use ellipsometry to investigate a transition in the morphology of a sphere-forming diblock copolymer thin-film system. At an interface the diblock morphology may differ from the bulk when the interfacial tension favours wetting of the minority domain, thereby inducing a sphere-to-lamella transition. In a small, favourable window in energetics, one may observe this transition simply by adjusting the temperature. Ellipsometry is ideally suited to the study of the transition because the additional interface created by the wetting layer affects the polarisation of light reflected from the sample. Here we study thin films of poly(butadiene-ethylene oxide) (PB-PEO), which order to form PEO minority spheres in a PB matrix. As temperature is varied, the reversible transition from a partially wetting layer of PEO spheres to a full wetting layer at the substrate is investigated.

Stochastic volatility jump-diffusions for European equity index dynamics

Major research on equity index dynamics has investigated only US indices (usually the S&P 500) and has provided contradictory results. In this paper a clarification and extension of that previous research is given. We find that European equity indices have quite different dynamics from the S&P 500. Each of the European indices considered may be satisfactorily modelled using either an affine model with price and volatility jumps or a GARCH volatility process without jumps. The S&P 500 dynamics are much more difficult to capture in a jump-diffusion framework.

Reversible addition of CO to coordinatively unsaturated high-spin iron(II) complexes

Several new coordinatively unsaturated iron(II) complexes of the types [Fe(EN-iPr)X2] (E = P, S, Se; X = Cl, Br) and [Fe(ON-iPr)2X]X containing bidentate EN ligands based on N-(2-pyridinyl)aminophosphines as well as oxo, thio, and seleno derivatives thereof were prepared and characterized by NMR spectroscopy and X-ray crystallography. Mössbauer spectroscopy and magnetization studies confirmed their high-spin nature with magnetic moments very close to 4.9 μB, reflecting the expected four unpaired d-electrons in all these compounds. Stable low-spin carbonyl complexes of the types [Fe(PN-iPr)2(CO)X]X (X = Cl, Br) and cis-CO,cis-Br-[Fe(PN-iPr)(CO)2X2] (X = Br) were obtained by reacting cis-Fe(CO)4X2 with the stronger PN donor ligands, but not with the weaker EN donor ligands (E = O, S, Se). Furthermore, the reactivity of [Fe(PN-iPr)X2] toward CO was investigated by IR spectroscopy. Whereas at room temperature no reaction took place, at −50 °C [Fe(PN-iPr)X2] added readily CO to form, depending on the nature of X, the mono- and dicarbonyl complexes [Fe(PN-iPr)(X)2(CO)] (X = Cl) and [Fe(PN-iPr)(CO)2X2] (X = Cl, Br), respectively. In the case of X = Br, two isomeric dicarbonyl complexes, namely, cis-CO,trans-Br-[Fe(PN-iPr)(CO)2Br2] (major species) and cis-CO,cis-Br-[Fe(PN-iPr)(CO)2Br2] (minor species), are formed. The addition of CO to [Fe(PN-iPr)X2] was investigated in detail by means of DFT/B3LYP calculations. This study strongly supports the experimental findings that at low temperature two isomeric low-spin dicarbonyl complexes are formed. For kinetic reasons cis,trans-[Fe(PN-iPr)(CO)2Br2] releases CO at elevated temperature, re-forming [Fe(PN-iPr)Br2], while the corresponding cis,cis isomer is stable under these conditions.

Spectral analysis of diffusions with jump boundary

In this paper we consider one-dimensional diffusions with constant coefficients in a finite interval with jump boundary and a certain deterministic jump distribution. We use coupling methods in order to identify the spectral gap in the case of a large drift and prove that there is a threshold drift above which the bottom of the spectrum no longer depends on the drift. As a corollary to our result we are able to answer two questions concerning elliptic eigenvalue problems with non-local boundary conditions formulated previously by Iddo Ben-Ari and Ross Pinsky.

On the spectral gap of Brownian motion with jump boundary

In this paper we consider the Brownian motion with jump boundary and present a new proof of a recent result of Li, Leung and Rakesh concerning the exact convergence rate in the one-dimensional case. Our methods are dierent and mainly probabilistic relying on coupling methods adapted to the special situation under investigation. Moreover we answer a question raised by Ben-Ari and Pinsky concerning the dependence of the spectral gap from the jump distribution in a multi-dimensional setting.