A quasispecies approach to the evolution of sexual replication in unicellular organisms

This study develops a simplified model describing the evolutionary dynamics of a population composed of obligate sexually and asexually reproducing, unicellular organisms. The model assumes that the organisms have diploid genomes consisting of two chromosomes, and that the sexual organisms replicate by first dividing into haploid intermediates, which then combine with other haploids, followed by the normal mitotic division of the resulting diploid into two new daughter cells. We assume that the fitness landscape of the diploids is analogous to the single-fitness-peak approach often used in single-chromosome studies. That is, we assume a master chromosome that becomes defective with just one point mutation. The diploid fitness then depends on whether the genome has zero, one, or two copies of the master chromosome. We also assume that only pairs of haploids with a master chromosome are capable of combining so as to produce sexual diploid cells, and that this process is described by second-order kinetics. We find that, in a range of intermediate values of the replication fidelity, sexually reproducing cells can outcompete asexual ones, provided the initial abundance of sexual cells is above some threshold value. The range of values where sexual reproduction outcompetes asexual reproduction increases with decreasing replication rate and increasing population density. We critically evaluate a common approach, based on a group selection perspective, used to study the competition between populations and show its flaws in addressing the evolution of sex problem.

ENERGY OF GLOBAL FRAMES

The energy of a unit vector field X on a closed Riemannian manifold M is defined as the energy of the section into T(1) M determined by X. For odd-dimensional spheres, the energy functional has an infimum for each dimension 2k + 1 which is not attained by any non-singular vector field for k > 1. For k = 1, Hopf vector fields are the unique minima. In this paper we show that for any closed Riemannian manifold, the energy of a frame defined on the manifold, possibly except on a finite subset, admits a lower bound in terms of the total scalar curvature of the manifold. In particular, for odd-dimensional spheres this lower bound is attained by a family of frames defined on the sphere minus one point and consisting of vector fields parallel along geodesics.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

A Floating-point Extended Kalman Filter Implementation for Autonomous Mobile Robots

Localization and Mapping are two of the most important capabilities for autonomous mobile robots and have been receiving considerable attention from the scientific computing community over the last 10 years. One of the most efficient methods to address these problems is based on the use of the Extended Kalman Filter (EKF). The EKF simultaneously estimates a model of the environment (map) and the position of the robot based on odometric and exteroceptive sensor information. As this algorithm demands a considerable amount of computation, it is usually executed on high end PCs coupled to the robot. In this work we present an FPGA-based architecture for the EKF algorithm that is capable of processing two-dimensional maps containing up to 1.8 k features at real time (14 Hz), a three-fold improvement over a Pentium M 1.6 GHz, and a 13-fold improvement over an ARM920T 200 MHz. The proposed architecture also consumes only 1.3% of the Pentium and 12.3% of the ARM energy per feature.

Influence of super-ohmic dissipation on a disordered quantum critical point

We investigate the combined influence of quenched randomness and dissipation on a quantum critical point with O(N) order-parameter symmetry. Utilizing a strong-disorder renormalization group, we determine the critical behavior in one space dimension exactly. For super-ohmic dissipation, we find a Kosterlitz-Thouless type transition with conventional (power-law) dynamical scaling. The dynamical critical exponent depends on the spectral density of the dissipative baths. We also discuss the Griffiths singularities, and we determine observables.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

NIELSEN COINCIDENCE THEORY OF FIBRE-PRESERVING MAPS AND DOLD`S FIXED POINT INDEX

Let M -> B, N -> B be fibrations and f(1), f(2): M -> N be a pair of fibre-preserving maps. Using normal bordism techniques we define an invariant which is an obstruction to deforming the pair f(1), f(2) over B to a coincidence free pair of maps. In the special case where the two fibrations axe the same and one of the maps is the identity, a weak version of our omega-invariant turns out to equal Dold`s fixed point index of fibre-preserving maps. The concepts of Reidemeister classes and Nielsen coincidence classes over B are developed. As an illustration we compute e.g. the minimal number of coincidence components for all homotopy classes of maps between S(1)-bundles over S(1) as well as their Nielsen and Reidemeister numbers.

Isoelectric Point Determination for Glossoscolex paulistus Extracellular Hemoglobin: Oligomeric Stability in Acidic pH and Relevance to Protein-Surfactant Interactions

The extracellular hemoglobin from Glossoscolex paulistus (HbGp) has a molecular mass of 3.6 M Da, It has a high oligomeric stability at pH 7.0 and low autoxidation rates, as compared to vertebrate hemoglobins. In this work, fluorescence and light scattering experiments were performed with the three oxidation forms of HbGp exposed to acidic pH. Our focus is on the HbGp stability at acidic pH and also on the determination of the isoelectric point (pI) of the protein. Our results show that the protein in the cyanomet form is more stable than in the other two forms, in the whole range. Our zeta-potential data are consistent with light scattering results. Average values apt obtained by different techniques were 5.6 +/- 0.5, 5.4 +/- 0.2 and 5.2 +/- 0.5 for the oxy, met, and cyanomet forms. Dynamic light scattering (DLS) experiments have shown that, at pH 6.0, the aggregation (oligomeric) state of oxy-, met- and cyanomet-HbGp remains the same as that at 7.0. The interaction between the oxy-HbGp and ionic surfactants at pH 5.0 and 6.0 was also monitored in the present study. At pH 5,0, below the protein pI, the effects of sodium dodecyl sulfate (SDS) and cetyltrimethylammonium chloride (CTAC) are inverted when compared to pH 7.0. For CTAC, in acid pH 5.0, no precipitation is observed, while for SDS an intense light scattering appears due to a precipitation process. HbGp interacts strongly with the cationic surfactant at pH 7.0 and with the anionic one at pH 5.0. This effect is due to the predominance, in the protein surface, of residues presenting opposite charges to the surfactant headgroups. This information can be relevant for the development of extracellular hemoglobin-based artificial blood substitutes.