A phylogenetic mixture model for detecting pattern-heterogeneity in gene sequence or character-state data

We describe a general likelihood-based 'mixture model' for inferring phylogenetic trees from gene-sequence or other character-state data. The model accommodates cases in which different sites in the alignment evolve in qualitatively distinct ways, but does not require prior knowledge of these patterns or partitioning of the data. We call this qualitative variability in the pattern of evolution across sites "pattern-heterogeneity" to distinguish it from both a homogenous process of evolution and from one characterized principally by differences in rates of evolution. We present studies to show that the model correctly retrieves the signals of pattern-heterogeneity from simulated gene-sequence data, and we apply the method to protein-coding genes and to a ribosomal 12S data set. The mixture model outperforms conventional partitioning in both these data sets. We implement the mixture model such that it can simultaneously detect rate- and pattern-heterogeneity. The model simplifies to a homogeneous model or a rate- variability model as special cases, and therefore always performs at least as well as these two approaches, and often considerably improves upon them. We make the model available within a Bayesian Markov-chain Monte Carlo framework for phylogenetic inference, as an easy-to-use computer program.

Mapping quantitative trait loci for flag leave senescence as a yield determinant in winter wheat under optimal and drought stressed environments

The timing of flag leaf senescence (FLS) is an important determinant of yield under stress and optimal environments. A doubled haploid population derived from crossing the photo period-sensitive variety Beaver,with the photo period-insensitive variety Soissons, varied significantly for this trait, measured as the percent green flag leaf area remaining at 14 days and 35 days after anthesis. This trait also showed a significantly positive correlation with yield under variable environmental regimes. QTL analysis based on a genetic map derived from 48 doubled haploid lines using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers, revealed the genetic control of this trait. The coincidence of QTL for senescence on chromosomes 2B and 2D under drought-stressed and optimal environments, respectively, indicate a complex genetic mechanism of this trait involving the re-mobilisation of resources from the source to the sink during senescence.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Analysis of the antigen combining site: Correlations between length and sequence composition of the hypervariable loops and the nature of the antigen

It has long been suggested that the overall shape of the antigen combining site (ACS) of antibodies is correlated with the nature of the antigen. For example, deep pockets are characteristic of antibodies that bind haptens, grooves indicate peptide binders, while antibodies that bind to proteins have relatively flat combining sites. In. 1996, MacCallum, Martin and Thornton used a fractal shape descriptor and showed a strong correlation of the shape of the binding region with the general nature of the antigen. However, the shape of the ACS is determined primarily by the lengths of the six complementarity-determining regions (CDRs). Here, we make a direct correlation between the lengths of the CDRs and the nature of the antigen. In addition, we show significant differences in the residue composition of the CDRs of antibodies that bind to different antigen classes. As well as helping us to understand the process of antigen recognition, autoimmune disease and cross-reactivity these results are of direct application in the design of antibody phage libraries and modification of affinity. (C) 2003 Elsevier Science Ltd. All rights reserved.

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.

Inter-simple sequence repeat and aggressiveness analyses revealed high genetic diversity, recombination and long-range dispersal in Fusarium culmorum

Inter-simple sequence repeat (ISSR) analysis and aggressiveness assays were used to investigate genetic variability within a global collection of Fusarium culmorum isolates. A set of four ISSR primers were tested, of which three primers amplified a total of 37 bands out of which 30 (81%) were polymorphic. The intraspecific diversity was high, ranging from four to 28 different ISSR genotypes for F. culmorum depending on the primer. The combined analysis of ISSR data revealed 59 different genotypes clustered into seven distinct clades amongst 75 isolates of F. culmorum examined. All the isolates were assayed to test their aggressiveness on a winter wheat cv. 'Armada'. A significant quantitative variation for aggressiveness was found among the isolates. The ISSR and aggressiveness variation existed on a macro- as well as micro-geographical scale. The data suggested a long-range dispersal of F. culmorum and indicated that this fungus may have been introduced into Canada from Europe. In addition to the high level of intraspecific diversity observed in F. culmorum, the index of multilocus association calculated using ISSR data indicated that reproduction in F. culmorum cannot be exclusively clonal and recombination is likely to occur.

Population subdivision and molecular sequence variation: Theory and analysis of Drosophila ananassae data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to Mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter Theta to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of Theta, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.

Recognition of polyimide sequence information by a molecular tweezer

Specific monomer sequences in aromatic copolyimides are recognized through their -stacking and hydrogen-bonding interactions with a sterically and electronically complementary molecular tweezer. These interactions enable the tweezer molecule to read monomer sequences comprising up to 27 aromatic rings by multiple adjacent binding to neighboring sites on the polymer chain.

Recognition of sequence-information in synthetic copolymer chains by a conformationally-constrained tweezer molecule

A novel type of tweezer molecule containing electron-rich 2-pyrenyloxy arms has been designed to exploit intramolecular hydrogen bonding in stabilising a preferred conformation for supramolecular complexation to complementary sequences in aromatic copolyimides. This tweezer-conformation is demonstrated by single-crystal X-ray analyses of the tweezer molecule itself and of its complex with an aromatic diimide model-compound. In terms of its ability to bind selectively to polyimide chains, the new tweezer molecule shows very high sensitivity to sequence effects. Thus, even low concentrations of tweezer relative to diimide units (<2.5 mol%) are sufficient to produce dramatic, sequence-related splittings of the pyromellitimide proton NMR resonances. These induced resonance-shifts arise from ring-current shielding of pyromellitimide protons by the pyrenyloxy arms of the tweezer-molecule, and the magnitude of such shielding is a function of the tweezer-binding constant for any particular monomer sequence. Recognition of both short-range and long-range sequences is observed, the latter arising from cumulative ring-current shielding of diimide protons by tweezer molecules binding at multiple adjacent sites on the copolymer chain.

Principles of sequence-recognition in aromatic polyimides

Pyrene-based molecular tweezers show sequence-specific binding to aromatic polyimides through sterically-controlled donor-acceptor pi-stacking and hydrogen bonding; H-1 NMR spectra of tweezer-complexes with polyimides having different sequence-restrictions show conclusively that the detection of long range sequence-information results from multiple tweezer-binding at adjacent imide residues.

Design of a turn-linker-turn foldamer by incorporating meta-amino benzoic acid in the middle of a helix forming hexapeptide sequence: A helix breaking approach

Single crystal X-ray diffraction studies reveal that the incorporation of meta-amino benzoic acid in the middle of a helix forming hexapeptide sequence such as in peptide I Boc-Ile(1)-Aib(2)-Val(3)-m-ABA(4)-Ile(5)-Aib(6)-Leu(7)-OMe (Aib: alpha-amino isobutyric acid: m-ABA: meta-amino benzoic acid) breaks the helix propagation to produce a turn-linker-turn (T-L-T) foldamer in the solid state. In the crystalline state two conformational isomers of peptide I self-assemble in antiparallel fashion through intermolecular hydrogen bonds and aromatic pi-pi interactions to form a molecular duplex. The duplexes are further interconnected through intermolecular hydrogen bonds to form a layer of peptides. The layers are stacked one on top of the other through van der Waals interactions to form hydrophilic channels filled with solvent methanol. (C) 2009 Elsevier B.V. All rights reserved.

Multiple lyotropic polymorphism of a PEG-peptide diblock copolymer in aqueous solution

Nematic and hexagonal columnar liquid crystal phase formation by a PEG-peptide conjugate is reported. The results are relevant to peptide-polymer Conjugates and bionanomaterial self-assembly (with relevance to PEGylated peptides used in therapeutic applications). The use of modified fragments of the amyloid beta peptide is especially interesting with respect to amyloid fibrillization and its control.

Polymorphism in [Co(SCN)(4)(ppz-H)(2)] (ppz, piperazine)

A one-pot reaction of [Co(NO3)(2)center dot 6H(2)O and piperazine] with NH4SCN/NaSCN in water-methanol (1:1) solvent leads to two polymorphs of [Co(SCN)(4)(ppz-H)(2)] (ppz, piperazine) (I and II). X-ray crystal structure reveals both have same space group but the differences in the alignment of pendant SCN- leads to two polymorphs. In I, trifurcated N-H...S hydrogen bonding plays a prominent role in crystal packing leading to S...S interactions between SCN fragments but in II, no such trifurcation arises and thereby the crystal packing occurs through hydrogen bonding interactions only leading to a distinctly different network topology. TG/DSC and FT-IR study reveal they are enantiotropically related. (c) 2007 Elsevier Ltd. All rights reserved.

Effect of sequence distribution on the morphology, crystallization, melting, and biodegradation of poly(epsilon-caprolactone-co-epsilon-caprolactam) copolymers

Two types of poly(epsilon-caprolactone (CLo)-co-poly(epsilon-caprolactam (CLa)) copolymers were prepared by catalyzed hydrolytic ring-opening polymerization. Both cyclic comonomers were added simultaneously in the reaction medium for the First type or materials where copolymers have a random distribution of counits, as evidenced by H-1 and C-13 NMR. For the second type of copolymers, the cyclic comonomers were added sequentially, yielding diblock poly(ester-amides). The materials were characterized by differential scanning calorimetry (DSC), wide- and small-angle X-ray scattering (WAXS and SAXS), and transmission and scanning electron microscopies (TEM and SEM). Their biodegradation in compost was also studied. All copolymers were found to be miscible by the absence of structure in the melt. TEM revealed that all samples exhibited a crystalline lamellar morphology. DSC and WAXS showed that in a wide composition range (CLo contents from 6 to 55%) only the CLa units were capable of crystallization in the random copolymers. The block copolymer samples only experience a small reduction of crystallization and melting temperature with composition, and this was attributed to a dilution effect caused by the miscible noncrystalline CLo units. The comparison between block and random copolymers provided a unique opportunity to distinguish the dilution effect of the CLo units on the crystallization and melting of the polyamide phase from the chemical composition effect in the random copolymers case, where the CLa sequences are interrupted statistically by the CLo units, making the crystallization of the polyamide strongly composition dependent. Finally, the enzymatic degradation of the copolymers in composted soil indicate a synergistic behavior where much faster degradation was obtained for random copolymers witha CLo content larger than 30% than for neat PCL.

Concomitant polymorphism of an antiferromagnetically coupled dicopper(II,II) complex with single strand helical assembly: synthesis, structure, DSC, magnetic and heterogeneous catalytic studies

Two concomitant polymorphic coordination complexes (dark blue - I and black - II) with the formula (Cu2C44H60N4O4) have been synthesized and characterized crystallographically. Magnetic measurements show the presence of a strong antiferromagnetic interaction and the 2J value corresponds extremely well to the theoretically calculated one, indicating the fact that it follows nicely the magneto-structural relationship. Immobilization of the copper(II) complex I on a 2D-hexagonal mesoporous silica showed good catalytic efficiency in the liquid phase partial oxidation of olefins in the presence of TBHP as an oxidant. (C) 2009 Elsevier Ltd. All rights reserved.