High genetic diversity in field isolates of Trypanosoma theileri assessed by analysis of cathepsin L-like sequences disclosed multiple and new genotypes infecting cattle in Thailand

In this study, we describe the first survey in Thailand of Trypanosoma theileri, a widespread and prevalent parasite of cattle that is transmitted by tabanid flies. Investigation of 210 bovine blood samples of Thai cattle from six farms by hematocrit centrifuge technique (HCT) revealed 14 samples with trypanosomes morphologically compatible to T. theileri. Additional animals were positive for T. theileri by PCR based on the Cathepsin L-like sequence (TthCATL-PCR) despite negative by HCT, indicating cryptic infections. Results revealed a prevalence of 26 +/- 15% (95% CI) of T. theileri infection. Additionally, 12 samples positive for T. theileri were detected in cattle from other 11 farms. From a total of 30 blood samples positive by HCT and/or PCR from 17 farms, seven were characterized to evaluate the genetic polymorphism of T. theileri through sequence analysis of PCR-amplified CATL DNA sequences. All CATL sequences of T. theileri from Thai cattle clustered with sequences of the previously described phylogenetic lineages TthI and TthII, supporting only two major lineages of T. theileri in cattle around the world. However, 11 of the 29 CATL sequences analyzed showed to be different, disclosing an unexpectedly large polymorphic genetic repertoire, with multiple genotypes of T. theileri not previously described in other countries circulating in Thai cattle. (C) 2011 Elsevier B.V. All rights reserved.

Sequence-specific reconstruction from fragmentary databases using seed sequences: implementation and validation on SAGE, proteome and generic sequencing data

Motivation: DNA assembly programs classically perform an all-against-all comparison of reads to identify overlaps, followed by a multiple sequence alignment and generation of a consensus sequence. If the aim is to assemble a particular segment, instead of a whole genome or transcriptome, a target-specific assembly is a more sensible approach. GenSeed is a Perl program that implements a seed-driven recursive assembly consisting of cycles comprising a similarity search, read selection and assembly. The iterative process results in a progressive extension of the original seed sequence. GenSeed was tested and validated on many applications, including the reconstruction of nuclear genes or segments, full-length transcripts, and extrachromosomal genomes. The robustness of the method was confirmed through the use of a variety of DNA and protein seeds, including short sequences derived from SAGE and proteome projects.

Trypanosoma rangeli isolates of bats from Central Brazil: Genotyping and phylogenetic analysis enable description of a new lineage using spliced-leader gene sequences

Trypanosoma rangeli infects several mammalian orders but has never confidently been described in Chiroptera, which are commonly parasitized by many trypanosome species. Here, we described trypanosomes from bats captured in Central Brazil identified as T rangeli,.T. dionisii, T cruzimarinkellei and T cruzi. Two isolates, Tra643 from Platyrrhinus lineatus and Tra1719 from Artibeus plamirostris were identified as T rangeli by morphological, biological and molecular methods, and confirmed by phylogenetic analyses. Analysis using SSU rDNA sequences clustered these bat trypanosomes together with T rangeli from other hosts, and separated them from other trypanosomes from bats. Genotyping based on length and sequence polymorphism of PCR-amplified intergenic spliced-leader gene sequences assigned Tra1719 to the lineage A whereas Tra643 was shown to be a new genotype and was assigned to the new lineage E. To our knowledge, these two isolates are the earliest T rangeli from bats and the first isolates from Central Brazil molecularly characterized. Rhodnius stali captured for this study was found infected by T rangeli and T cruzi. (c) 2008 Elsevier B.V. All rights reserved.

Phylogeny of snake trypanosomes inferred by SSU rDNA sequences, their possible transmission by phlebotomines, and taxonomic appraisal by molecular, cross-infection and morphological analysis

Blood examination by microhaematocrit and haemoculture of 459 snakes belonging to 37 species revealed 24% trypanosome prevalence in species of Viperidae (Crotalus durissus and Bothrops jararaca) and Colubridae (Pseudoboa nigra). Trypanosome cultures from C. durissus and P. nigra were behaviourally and morphologically indistinguishable. In addition, the growth and morphological features of a trypanosome from the sand fly Viannaniyia tuberculata were similar to those of snake isolates. Cross-infection experiments revealed a lack of host restriction, as snakes of 3 species were infected with the trypanosome from C. durissus. Phylogeny based on ribosomal sequences revealed that snake trypanosomes clustered together with the sand fly trypanosome, forming a new phylogenetic lineage within Trypanosoma closest to a clade of lizard trypanosomes transmitted by sand flies dagger. The clade of trypanosomes from snakes and lizards suggests an association between the evolutionary histories of these trypanosomes and their squamate hosts. Moreover, data strongly indicated that these trypanosomes are transmitted by sand flies. The flaws of the current taxonomy of snake trypanosomes are discussed, and the need for molecular parameters to be adopted is emphasized. To our knowledge, this is the first molecular phylogenetic study of snake trypanosomes.

On comparing two sequences of numbers and its applications to clustering analysis

A conceptual problem that appears in different contexts of clustering analysis is that of measuring the degree of compatibility between two sequences of numbers. This problem is usually addressed by means of numerical indexes referred to as sequence correlation indexes. This paper elaborates on why some specific sequence correlation indexes may not be good choices depending on the application scenario in hand. A variant of the Product-Moment correlation coefficient and a weighted formulation for the Goodman-Kruskal and Kendall`s indexes are derived that may be more appropriate for some particular application scenarios. The proposed and existing indexes are analyzed from different perspectives, such as their sensitivity to the ranks and magnitudes of the sequences under evaluation, among other relevant aspects of the problem. The results help suggesting scenarios within the context of clustering analysis that are possibly more appropriate for the application of each index. (C) 2008 Elsevier Inc. All rights reserved.

Using stochastic volatility models to analyse weekly ozone averages in Mexico City

In this paper we make use of some stochastic volatility models to analyse the behaviour of a weekly ozone average measurements series. The models considered here have been used previously in problems related to financial time series. Two models are considered and their parameters are estimated using a Bayesian approach based on Markov chain Monte Carlo (MCMC) methods. Both models are applied to the data provided by the monitoring network of the Metropolitan Area of Mexico City. The selection of the best model for that specific data set is performed using the Deviance Information Criterion and the Conditional Predictive Ordinate method.

On the cutting stock problem under stochastic demand

This paper addresses the one-dimensional cutting stock problem when demand is a random variable. The problem is formulated as a two-stage stochastic nonlinear program with recourse. The first stage decision variables are the number of objects to be cut according to a cutting pattern. The second stage decision variables are the number of holding or backordering items due to the decisions made in the first stage. The problem`s objective is to minimize the total expected cost incurred in both stages, due to waste and holding or backordering penalties. A Simplex-based method with column generation is proposed for solving a linear relaxation of the resulting optimization problem. The proposed method is evaluated by using two well-known measures of uncertainty effects in stochastic programming: the value of stochastic solution-VSS-and the expected value of perfect information-EVPI. The optimal two-stage solution is shown to be more effective than the alternative wait-and-see and expected value approaches, even under small variations in the parameters of the problem.

A new scale-invariant ratio and finite-size scaling for the stochastic susceptible-infected-recovered model

The critical behavior of the stochastic susceptible-infected-recovered model on a square lattice is obtained by numerical simulations and finite-size scaling. The order parameter as well as the distribution in the number of recovered individuals is determined as a function of the infection rate for several values of the system size. The analysis around criticality is obtained by exploring the close relationship between the present model and standard percolation theory. The quantity UP, equal to the ratio U between the second moment and the squared first moment of the size distribution multiplied by the order parameter P, is shown to have, for a square system, a universal value 1.0167(1) that is the same for site and bond percolation, confirming further that the SIR model is also in the percolation class.

Stochastic lattice gas model describing the dynamics of the SIRS epidemic process

We study a stochastic process describing the onset of spreading dynamics of an epidemic in a population composed of individuals of three classes: susceptible (S), infected (I), and recovered (R). The stochastic process is defined by local rules and involves the following cyclic process: S -> I -> R -> S (SIRS). The open process S -> I -> R (SIR) is studied as a particular case of the SIRS process. The epidemic process is analyzed at different levels of description: by a stochastic lattice gas model and by a birth and death process. By means of Monte Carlo simulations and dynamical mean-field approximations we show that the SIRS stochastic lattice gas model exhibit a line of critical points separating the two phases: an absorbing phase where the lattice is completely full of S individuals and an active phase where S, I and R individuals coexist, which may or may not present population cycles. The critical line, that corresponds to the onset of epidemic spreading, is shown to belong in the directed percolation universality class. By considering the birth and death process we analyze the role of noise in stabilizing the oscillations. (C) 2009 Elsevier B.V. All rights reserved.

CRITICAL BEHAVIOR AND THRESHOLD OF COEXISTENCE OF A PREDATOR-PREY STOCHASTIC MODEL IN A 2D LATTICE

We investigate the critical behavior of a stochastic lattice model describing a predator-prey system. By means of Monte Carlo procedure we simulate the model defined on a regular square lattice and determine the threshold of species coexistence, that is, the critical phase boundaries related to the transition between an active state, where both species coexist and an absorbing state where one of the species is extinct. A finite size scaling analysis is employed to determine the order parameter, order parameter fluctuations, correlation length and the critical exponents. Our numerical results for the critical exponents agree with those of the directed percolation universality class. We also check the validity of the hyperscaling relation and present the data collapse curves.

Glassy states in the stochastic Potts model

We have studied by numerical simulations the relaxation of the stochastic seven-state Potts model after a quench from a high temperature down to a temperature below the first-order transition. For quench temperatures just below the transition temperature the phase ordering occurs by simple coarsening under the action of surface tension. For sufficient low temperatures however the straightening of the interface between domains drives the system toward a metastable disordered state, identified as a glassy state. Escaping from this state occurs, if the quench temperature is nonzero, by a thermal activated dynamics that eventually drives the system toward the equilibrium state. (C) 2009 Elsevier B.V. All rights reserved.

The stochastic nature of predator-prey cycles

We study by numerical simulations the time correlation function of a stochastic lattice model describing the dynamics of coexistence of two interacting biological species that present time cycles in the number of species individuals. Its asymptotic behavior is shown to decrease in time as a sinusoidal exponential function from which we extract the dominant eigenvalue of the evolution operator related to the stochastic dynamics showing that it is complex with the imaginary part being the frequency of the population cycles. The transition from the oscillatory to the nonoscillatory behavior occurs when the asymptotic behavior of the time correlation function becomes a pure exponential, that is, when the real part of the complex eigenvalue equals a real eigenvalue. We also show that the amplitude of the undamped oscillations increases with the square root of the area of the habitat as ordinary random fluctuations. (C) 2009 Elsevier B.V. All rights reserved.

Mean Motion in Stochastic Plasmas With a Space-Dependent Diffusion Coefficient

Radial transport in the tokamap, which has been proposed as a simple model for the motion in a stochastic plasma, is investigated. A theory for previous numerical findings is presented. The new results are stimulated by the fact that the radial diffusion coefficients is space-dependent. The space-dependence of the transport coefficient has several interesting effects which have not been elucidated so far. Among the new findings are the analytical predictions for the scaling of the mean radial displacement with time and the relation between the Fokker-Planck diffusion coefficient and the diffusion coefficient from the mean square displacement. The applicability to other systems is also discussed. (c) 2009 WILEY-VCH GmbH & Co. KGaA, Weinheim

GEOCHEMISTRY AND U-PB ZIRCON AGE OF THE INTERNAL OSSA-MORENA ZONE OPHIOLITE SEQUENCES: A REMNANT OF RHEIC OCEAN IN SW IBERIA

The Early Paleozoic geodynamic evolution in SW Iberia is believed to have been dominated by the opening of the Rheic Ocean. The Rheic Ocean is generally accepted to have resulted from the drift of peri-Gondwanan terranes such as Avalonia from the northern margin of Gondwana during Late Cambrian-Early Ordovician times. The closure of the Rheic Ocean was the final result of a continent-continent collision between Gondwana and Laurussia that produced the Variscan orogen. The Ossa-Morena Zone is a peri-Gondwana terrane, which preserves spread fragments of ophiolites - the Internal Ossa-Morena Zones Ophiolite Sequences (IOMZOS). The final patchwork of the IOMZOS shows a complete oceanic lithospheric sequence with geochemical characteristics similar to the ocean-floor basalts, without any orogenic fingerprint and/or crustal contamination. The IOMZOS were obducted and imbricated with high pressure lithologies. Based on structural, petrological and whole-rock geochemical data, the authors argue that the IOMZOS represent fragments of the oceanic lithosphere from the Rheic Ocean. Zircon SHRIMP U-Pb geochronological data on metagabbros point to an age of ca. 480 Ma for IOMZOS, providing evidence of a well-developed ocean in SW Iberia during this period, reinforcing the interpretation of the Rheic Ocean as a wide ocean among the peri-Gondwanan terranes during Early Ordovician times.

STABILITY ANALYSIS WITH APPLICATIONS OF A TWO-DIMENSIONAL DYNAMICAL SYSTEM ARISING FROM A STOCHASTIC MODEL FOR AN ASSET MARKET

We analyze the stability properties of equilibrium solutions and periodicity of orbits in a two-dimensional dynamical system whose orbits mimic the evolution of the price of an asset and the excess demand for that asset. The construction of the system is grounded upon a heterogeneous interacting agent model for a single risky asset market. An advantage of this construction procedure is that the resulting dynamical system becomes a macroscopic market model which mirrors the market quantities and qualities that would typically be taken into account solely at the microscopic level of modeling. The system`s parameters correspond to: (a) the proportion of speculators in a market; (b) the traders` speculative trend; (c) the degree of heterogeneity of idiosyncratic evaluations of the market agents with respect to the asset`s fundamental value; and (d) the strength of the feedback of the population excess demand on the asset price update increment. This correspondence allows us to employ our results in order to infer plausible causes for the emergence of price and demand fluctuations in a real asset market. The employment of dynamical systems for studying evolution of stochastic models of socio-economic phenomena is quite usual in the area of heterogeneous interacting agent models. However, in the vast majority of the cases present in the literature, these dynamical systems are one-dimensional. Our work is among the few in the area that construct and study analytically a two-dimensional dynamical system and apply it for explanation of socio-economic phenomena.