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.

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.

Synthesis and properties of polyaromatic dendrimers possessing a repetitive amide-ester coupling sequence

A series of novel polyaromatic dendrimers that feature tris-(2-ethylamino)amine as the central core unit has been synthesized up to the third generation by employing a convergent growth strategy. The building blocks 1,3-diamino-2-hydroxypropane and 4-carboxybenzaldehyde were used for dendron construction, a process that involved the cyclic repetition of esterification, oxidation and selective amidation steps. Molecular modelling of this class of dendrimers has been used to predict potential solution state conformations employing molecular mechanics and molecular dynamic simulations. In addition, the results of preliminary metal binding studies using the first generation dendritic system are also outlined. (C) 2003 Elsevier Science Ltd. All rights reserved.

A comparison of two time-domain analysis procedures in the determination of VO2 kinetics by pseudorandom binary sequence exercise testing

The purpose of this study was to apply and compare two time-domain analysis procedures in the determination of oxygen uptake (VO2) kinetics in response to a pseudorandom binary sequence (PRBS) exercise test. PRBS exercise tests have typically been analysed in the frequency domain. However, the complex interpretation of frequency responses may have limited the application of this procedure in both sporting and clinical contexts, where a single time measurement would facilitate subject comparison. The relative potential of both a mean response time (MRT) and a peak cross-correlation time (PCCT) was investigated. This study was divided into two parts: a test-retest reliability study (part A), in which 10 healthy male subjects completed two identical PRBS exercise tests, and a comparison of the VO2 kinetics of 12 elite endurance runners (ER) and 12 elite sprinters (SR; part B). In part A, 95% limits of agreement were calculated for comparison between MRT and PCCT. The results of part A showed no significant difference between test and retest as assessed by MRT [mean (SD) 42.2 (4.2) s and 43.8 (6.9) s] or by PCCT [21.8 (3.7) s and 22.7 (4.5) s]. Measurement error (%) was lower for MRT in comparison with PCCT (16% and 25%, respectively). In part B of the study, the VO2 kinetics of ER were significantly faster than those of SR, as assessed by MRT [33.4 (3.4) s and 39.9 (7.1) s, respectively; P<0.01] and PCCT [20.9 (3.8) s and 24.8 (4.5) s; P < 0.05]. It is possible that either analysis procedure could provide a single test measurement Of VO2 kinetics; however, the greater reliability of the MRT data suggests that this method has more potential for development in the assessment Of VO2 kinetics by PRBS exercise testing.

Isolation of culturable yeasts from market wines and evaluation of the 5 center dot 8S-ITS rDNA sequence analysis for identification purposes

Aims: To test the possibility that wines available in the marketplace may contain culturable yeasts and to evaluate the 5.8S-ITS rDNA sequence analysis as adequate means for the identification of isolates. Methods and Results: As a case study, typical Greek wines were surveyed. Sequence analysis of the 5.8S-ITS rDNA was tested for its robustness in species or strain identification. Sixteen isolates could be assigned into the species Brettanomyces bruxellensis, Saccharomyces cerevisiae and Rhodotorula pinicola, whereas four isolates could not be safely identified. B. bruxellensis was the dominant species present in house wines, while non-Saccharomyces sp. were viable in aged wines of high alcohol content. Conclusions: Yeast population depends on postfermentation procedures or storage conditions. Although 5.8S-ITS rDNA sequence analysis is generally a rapid method to identify wine yeast isolates at the species level, or even below that, it may not be sufficient for some genera. Significance and Impact of the Study: This is the first report to show that commercial wines may possess diverse and potentially harmful yeast populations. The knowledge of yeasts able to reside in this niche environment is essential towards integrated quality assurance programmes. For selected species, the 5.8S-ITS rDNA sequence analysis is a rapid and accurate means.

Location, location, location: development of spatiotemporal sequence learning in infancy

We investigated infants' sensitivity to spatiotemporal structure. In Experiment 1, circles appeared in a statistically defined spatial pattern. At test 11-month-olds, but not 8-month-olds, looked longer at a novel spatial sequence. Experiment 2 presented different color/shape stimuli, but only the location sequence was violated during test; 8-month-olds preferred the novel spatial structure, but 5-month-olds did not. In Experiment 3, the locations but not color/shape pairings were constant at test; 5-month-olds showed a novelty preference. Experiment 4 examined "online learning": We recorded eye movements of 8-month-olds watching a spatiotemporal sequence. Saccade latencies to predictable locations decreased. We argue that temporal order statistics involving informative spatial relations become available to infants during the first year after birth, assisted by multiple cues.

Global attractivity in a recursive sequence

In this paper, we Study the invariant intervals, the globally attractivity of the two equilibrium points, and the oscillatory behavior of tile solutions of the difference equation x(n =) ax(n-1) - bx(n-2)/c + x(n-2), n = 1,2,......, where a, b. c > 0. (C) 2003 Elsevier Inc. All rights reserved.

Global asymptotic stability in a rational recursive sequence

In this paper, we study the global stability of the difference equation x(n) = a + bx(n-1) + cx(n-1)(2)/d - x(n-2), n = 1,2,....., where a, b greater than or equal to 0 and c, d > 0. We show that one nonnegative equilibrium point of the equation is a global attractor with a basin that is determined by the parameters, and every positive Solution of the equation in the basin exponentially converges to the attractor. (C) 2003 Elsevier Inc. All rights reserved.