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.

Detecting recombination in 4-taxa DNA sequence alignments with Bayesian hidden Markov models and Markov chain Monte Carlo

This article presents a statistical method for detecting recombination in DNA sequence alignments, which is based on combining two probabilistic graphical models: (1) a taxon graph (phylogenetic tree) representing the relationship between the taxa, and (2) a site graph (hidden Markov model) representing interactions between different sites in the DNA sequence alignments. We adopt a Bayesian approach and sample the parameters of the model from the posterior distribution with Markov chain Monte Carlo, using a Metropolis-Hastings and Gibbs-within-Gibbs scheme. The proposed method is tested on various synthetic and real-world DNA sequence alignments, and we compare its performance with the established detection methods RECPARS, PLATO, and TOPAL, as well as with two alternative parameter estimation schemes.

Self assembly of a model amphiphilic phenylalanine peptide/polyethylene glycol block copolymer in aqueous solution

There has been great interest recently in peptide amphiphiles and block copolymers containing biomimetic peptide sequences due to applications in bionanotechnology. We investigate the self-assembly of the peptide-PEG amphiphile FFFF-PEG5000 containing the hydrophobic sequence of four phenylalanine residues conjugated to PEG of molar mass 5000. This serves as a simple model peptide amphiphile. At very low concentration, association of hydrophobic aromatic phenylalanine residues occurs, as revealed by circular dichroism and UV/vis fluorescence experiments. A critical aggregation concentration associated with the formation of hydrophobic domains is determined through pyrene fluorescence assays. At higher concentration, defined beta-sheets develop as revealed by FTIR spectroscopy and X-ray diffraction. Transmission electron microscopy reveals self-assembled straight fibril structures. These are much shorter than those observed for amyloid peptides, the finite length may be set by the end cap energy due to the hydrophobicity of phenylalanine. The combination of these techniques points to different aggregation processes depending on concentration. Hydrophobic association into irregular aggregates occurs at low concentration, well-developed beta-sheets only developing at higher concentration. Drying of FFFF-PEG5000 solutions leads to crystallization of PEG, as confirmed by polarized optical microscopy (POM), FTIR and X-ray diffraction (XRD). PEG crystallization does not disrupt local beta-sheet structure (as indicated by FTIR and XRD). However on longer lengthscales the beta-sheet fibrillar structure is perturbed because spheruilites from PEG crystallization are observed by POM. (C) 2009 Elsevier B.V. All rights reserved.

Spontaneous condensation in DNA-polystyrene-b-poly(l-lysine) polyelectrolyte block copolymer mixtures

We investigated the condensation of calf thymus DNA by amphiphilic polystyrene(m)-b-poly(l-lysine)(n) block copolymers (PSm-b- PLys(n), m, n = degree of polymerization), using small-angle X-ray scattering, polarized optical microscopy and laser scanning confocal microscopy. Microscopy studies showed that the DNA condenses in the form of fibrillar precipitates, with an irregular structure, due to electrostatic interactions between PLys and DNA. This is not modified by the presence of hydrophobic PS block. Scattering experiments show that the structure of the polyplexes corresponds to a local order of DNA rods which becomes more compact upon increasing n. It can be concluded that for DNA/ PSm-b- PLys(n) polyplexes, the balance between the PLys block length and the excess charge in the system plays an essential role in the formation of a liquid crystalline phase.

Pressure effects revealed by small angle neutron scattering on block copolymer gels

We study the effects of hydrostatic pressure (P) on aqueous solutions and gels of the block copolymer B20E610 (E, oxyethylene; B, oxybutylene; subscripts, number of repeats), by performing simultaneous small angle neutron scattering/pressure experiments. Micellar cubic gels were studied for 9.5 and 4.5 wt% B20E610 at T = 20-80 and 35-55 degrees C, respectively, while micellar isotropic solutions where Studied for 4.5 wt% B20E610 at T > 55 degrees C. We observed that the interplanar distance d(110) (cubic unit cell parameter a = root 2d(110)) decreases while the correlation length of the Cubic order (delta) increases, upon increasing P at a fixed T for 9.5 wt% B20E610. The construction of master Curves for d(110) and delta corresponding to 9.5 wt% B20E610 proved the correlation between changes in T and P. Neither d(110) and delta nor the cubic-isotropic phase transition temperature was affected by the applied pressure for 4.5 wt% B20E610. The dramatic contrast between the pressure-induced behavior observed for 9.5 and 4.5 wt% B20E610 suggests that pressure induced effects might be more effectively transmitted through samples that present wider domains of cubic structure order (9.5 wt% compared to 4.5 wt% B20E610).

Wormlike micelle formation and flow alignment of a pluronic block copolymer in aqueous solution

The self-assembly into wormlike micelles of a poly(ethylene oxide)-b-poly(propylene oxide)-b-poly(ethylene oxide) triblock copolymer Pluronic P84 in aqueous salt solution (2 M NaCl) has been studied by rheology, small-angle X-ray and neutron scattering (SAXS/SANS), and light scattering. Measurements of the flow curves by controlled stress rheometry indicated phase separation under flow. SAXS on solutions subjected to capillary flow showed alignment of micelles at intermediate shear rates, although loss of alignment was observed for high shear rates. For dilute solutions, SAXS and static light scattering data on unaligned samples could be superposed over three decades in scattering vector, providing unique information on the wormlike micelle structure over several length scales. SANS data provided information on even shorter length scales, in particular, concerning "blob" scattering from the micelle corona. The data could be modeled based on a system of semiflexible self-avoiding cylinders with a circular cross-section, as described by the wormlike chain model with excluded volume interactions. The micelle structure was compared at two temperatures close to the cloud point (47 degrees C). The micellar radius was found not to vary with temperature in this region, although the contour length increased with increasing temperature, whereas the Kuhn length decreased. These variations result in an increase of the low-concentration radius of gyration with increasing temperature. This was consistent with dynamic light scattering results, and, applying theoretical results from the literature, this is in agreement with an increase in endcap energy due to changes in hydration of the poly(ethylene oxide) blocks as the temperature is increased.

Modeling the structure of amorphous MoS3: a neutron diffraction and reverse Monte Carlo study

A model for the structure of amorphous molybdenum trisulfide, a-MoS3, has been created using reverse Monte Carlo methods. This model, which consists of chains Of MoS6 units sharing three sulfurs with each of its two neighbors and forming alternate long, nonbonded, and short, bonded, Mo-Mo separations, is a good fit to the neutron diffraction data and is chemically and physically realistic. The paper identifies the limitations of previous models based on Mo-3 triangular clusters in accounting for the available experimental data.

Ellipsometric study of adsorption on nanopatterned block copolymer substrates

We report ellipsometrically obtained adsorption isotherms for a carefully chosen test liquid on block copolymer films of Kraton G1650, compared with adsorption isotherms on homogeneous films of the constituent polymers. Standard atomic force microscopy images imply the outer surface of Kraton G1650 is chemically patterned on the nanoscale, but this could instead be a reflection of structure buried beneath a 10 nm layer of the lower energy component. Our test liquid was chosen on the basis that it did not dissolve in either component and in addition that it was nonwetting on the lower energy polymer while forming thick adsorbed films on pure substrates of the higher energy component. Our ellipsometry data for Kraton G1650 rule out the presence of segregation by the lower energy constituent to the outer surface, implying a mixed surface consistent with Cassie's law. We discuss implications of our findings and related work for the outer surface structures of block copolymer films.

Structure, rheology and shear alignment of Pluronic block copolymer mixtures

The structure and flow behaviour of binary mixtures of Pluronic block copolymers P85 and P123 is investigated by small-angle scattering, rheometry and mobility tests. Micelle dimensions are probed by dynamic light scattering. The micelle hydrodynamic radius for the 50/50 mixture is larger than that for either P85 or P123 alone, Clue to the formation of mixed micelles with a higher association number. The phase diagram for 50/50 mixtures contains regions Of Cubic and hexagonal phases similar to those for the parent homopolymers, however the region of stability of the cubic phase is enhanced at low temperature and concentrations above 40 wt%. This is ascribed to favourable packing of the mixed micelles containing core blocks with two different chain lengths, but similar corona chain lengths. The shear flow alignment of face-centred cubic and hexagonal phases is probed by in situ small-angle X-ray or neutron scattering with simultaneous rheology. The hexagonal phase can be aligned using steady shear in a Couette geometry, however the high modulus Cubic phase cannot be aligned well in this way. This requires the application of oscillatory shear or compression. (C) 2008 Elsevier Inc. All rights reserved.

In situ formation of gold nanoparticles with a thermoresponsive block copolymer corona

Gold nanoparticles with a polymer coating exhibiting large and reversible thermoresponsiveness are prepared via a one-pot synthesis method using narrow polydispersity thermoresponsive block copolymers. (C) 2007 Elsevier B.V. All rights reserved.

Finding the smallest eigenvalue by the inverse Monte Carlo method with refinement

Finding the smallest eigenvalue of a given square matrix A of order n is computationally very intensive problem. The most popular method for this problem is the Inverse Power Method which uses LU-decomposition and forward and backward solving of the factored system at every iteration step. An alternative to this method is the Resolvent Monte Carlo method which uses representation of the resolvent matrix [I -qA](-m) as a series and then performs Monte Carlo iterations (random walks) on the elements of the matrix. This leads to great savings in computations, but the method has many restrictions and a very slow convergence. In this paper we propose a method that includes fast Monte Carlo procedure for finding the inverse matrix, refinement procedure to improve approximation of the inverse if necessary, and Monte Carlo power iterations to compute the smallest eigenvalue. We provide not only theoretical estimations about accuracy and convergence but also results from numerical tests performed on a number of test matrices.