Evolutionary and convergence stability for continuous phenotypes in finite populations derived from two-allele models.

The evolution of a quantitative phenotype is often envisioned as a trait substitution sequence where mutant alleles repeatedly replace resident ones. In infinite populations, the invasion fitness of a mutant in this two-allele representation of the evolutionary process is used to characterize features about long-term phenotypic evolution, such as singular points, convergence stability (established from first-order effects of selection), branching points, and evolutionary stability (established from second-order effects of selection). Here, we try to characterize long-term phenotypic evolution in finite populations from this two-allele representation of the evolutionary process. We construct a stochastic model describing evolutionary dynamics at non-rare mutant allele frequency. We then derive stability conditions based on stationary average mutant frequencies in the presence of vanishing mutation rates. We find that the second-order stability condition obtained from second-order effects of selection is identical to convergence stability. Thus, in two-allele systems in finite populations, convergence stability is enough to characterize long-term evolution under the trait substitution sequence assumption. We perform individual-based simulations to confirm our analytic results.

Functional analysis of pyochelin-/enantiopyochelin-related genes from a pathogenicity island of Pseudomonas aeruginosa strain PA14.

Genomic islands are foreign DNA blocks inserted in so-called regions of genomic plasticity (RGP). Depending on their gene content, they are classified as pathogenicity, symbiosis, metabolic, fitness or resistance islands, although a detailed functional analysis is often lacking. Here we focused on a 34-kb pathogenicity island of Pseudomonas aeruginosa PA14 (PA14GI-6), which is inserted at RGP5 and carries genes related to those for pyochelin/enantiopyochelin biosynthesis. These enantiomeric siderophores of P. aeruginosa and certain strains of Pseudomonas protegens are assembled by a thiotemplate mechanism from salicylate and two molecules of cysteine. The biochemical function of several proteins encoded by PA14GI-6 was investigated by a series of complementation analyses using mutants affected in potential homologs. We found that PA14_54940 codes for a bifunctional salicylate synthase/salicyl-AMP ligase (for generation and activation of salicylate), that PA14_54930 specifies a dihydroaeruginoic acid (Dha) synthetase (for coupling salicylate with a cysteine-derived thiazoline ring), that PA14_54910 produces a type II thioesterase (for quality control), and that PA14_54880 encodes a serine O-acetyltransferase (for increased cysteine availability). The structure of the PA14GI-6-specified metabolite was determined by mass spectrometry, thin-layer chromatography, and HPLC as (R)-Dha, an iron chelator with antibacterial, antifungal and antitumor activity. The conservation of this genomic island in many clinical and environmental P. aeruginosa isolates of different geographical origin suggests that the ability for Dha production may confer a selective advantage to its host.

A non-circadian role for clock-genes in sleep homeostasis: a strain comparison.

BACKGROUND: We have previously reported that the expression of circadian clock-genes increases in the cerebral cortex after sleep deprivation (SD) and that the sleep rebound following SD is attenuated in mice deficient for one or more clock-genes. We hypothesized that besides generating circadian rhythms, clock-genes also play a role in the homeostatic regulation of sleep. Here we follow the time course of the forebrain changes in the expression of the clock-genes period (per)-1, per2, and of the clock-controlled gene albumin D-binding protein (dbp) during a 6 h SD and subsequent recovery sleep in three inbred strains of mice for which the homeostatic sleep rebound following SD differs. We reasoned that if clock genes are functionally implicated in sleep homeostasis then the SD-induced changes in gene expression should vary according to the genotypic differences in the sleep rebound. RESULTS: In all three strains per expression was increased when animals were kept awake but the rate of increase during the SD as well as the relative increase in per after 6 h SD were highest in the strain for which the sleep rebound was smallest; i.e., DBA/2J (D2). Moreover, whereas in the other two strains per1 and per2 reverted to control levels with recovery sleep, per2 expression specifically, remained elevated in D2 mice. dbp expression increased during the light period both during baseline and during SD although levels were reduced during the latter condition compared to baseline. In contrast to per2, dbp expression reverted to control levels with recovery sleep in D2 only, whereas in the two other strains expression remained decreased. CONCLUSION: These findings support and extend our previous findings that clock genes in the forebrain are implicated in the homeostatic regulation of sleep and suggest that sustained, high levels of per2 expression may negatively impact recovery sleep.

Biocontrol strain Pseudomonas fluorescens CHA0 and its genetically modified derivative with enhanced biocontrol capability exert comparable effects on the structure of a Sinorhizobium meliloti population in a gnotobiotic system

The impact of biocontrol strain Pseudomonas fluorescens CHA0 and of its genetically modified, antibiotic-overproducing derivative CHA0/pME3424 on a reconstructed population of the plant-beneficial Sinorhizobium meliloti bacteria was assessed in gnotobiotic systems. In sterile soil, the final density of the reconstructed S. meliloti population decreased by more than one order of magnitude in the presence of either of the Pseudomonas strains when compared to a control without addition of P. fluorescens. Moreover, there was a change in the proportion of each individual S. meliloti strain within the population. Plant tests also revealed changes in the nodulating S. meliloti population in the presence of strains CHA0 or CHA0/pME3424. In both treatments one S. meliloti strain, f43, was significantly reduced in its root nodule occupancy. Analysis of alfalfa yields showed a slight but statistically significant increase in shoot dry weight when strain CHA0 was added to the reconstructed S. meliloti population whereas no such effect was observed with CHA0/pME3424.

Hybrid Multiscale Finite Volume method for two-phase flow in porous media

We present a novel hybrid (or multiphysics) algorithm, which couples pore-scale and Darcy descriptions of two-phase flow in porous media. The flow at the pore-scale is described by the Navier?Stokes equations, and the Volume of Fluid (VOF) method is used to model the evolution of the fluid?fluid interface. An extension of the Multiscale Finite Volume (MsFV) method is employed to construct the Darcy-scale problem. First, a set of local interpolators for pressure and velocity is constructed by solving the Navier?Stokes equations; then, a coarse mass-conservation problem is constructed by averaging the pore-scale velocity over the cells of a coarse grid, which act as control volumes; finally, a conservative pore-scale velocity field is reconstructed and used to advect the fluid?fluid interface. The method relies on the localization assumptions used to compute the interpolators (which are quite straightforward extensions of the standard MsFV) and on the postulate that the coarse-scale fluxes are proportional to the coarse-pressure differences. By numerical simulations of two-phase problems, we demonstrate that these assumptions provide hybrid solutions that are in good agreement with reference pore-scale solutions and are able to model the transition from stable to unstable flow regimes. Our hybrid method can naturally take advantage of several adaptive strategies and allows considering pore-scale fluxes only in some regions, while Darcy fluxes are used in the rest of the domain. Moreover, since the method relies on the assumption that the relationship between coarse-scale fluxes and pressure differences is local, it can be used as a numerical tool to investigate the limits of validity of Darcy's law and to understand the link between pore-scale quantities and their corresponding Darcy-scale variables.

Random mating with a finite number of matings.

Random mating is the null model central to population genetics. One assumption behind random mating is that individuals mate an infinite number of times. This is obviously unrealistic. Here we show that when each female mates a finite number of times, the effective size of the population is substantially decreased.

Oxygen-sensing reporter strain of Pseudomonas fluorescens for monitoring the distribution of low-oxygen habitats in soil.

The root-colonizing bacterium Pseudomonas fluorescens CHA0 was used to construct an oxygen-responsive biosensor. An anaerobically inducible promoter of Pseudomonas aeruginosa, which depends on the FNR (fumarate and nitrate reductase regulation)-like transcriptional regulator ANR (anaerobic regulation of arginine deiminase and nitrate reductase pathways), was fused to the structural lacZ gene of Escherichia coli. By inserting the reporter fusion into the chromosomal attTn7 site of P. fluorescens CHA0 by using a mini-Tn7 transposon, the reporter strain, CHA900, was obtained. Grown in glutamate-yeast extract medium in an oxystat at defined oxygen levels, the biosensor CHA900 responded to a decrease in oxygen concentration from 210 x 10(2) Pa to 2 x 10(2) Pa of O(2) by a nearly 100-fold increase in beta-galactosidase activity. Half-maximal induction of the reporter occurred at about 5 x 10(2) Pa. This dose response closely resembles that found for E. coli promoters which are activated by the FNR protein. In a carbon-free buffer or in bulk soil, the biosensor CHA900 still responded to a decrease in oxygen concentration, although here induction was about 10 times lower and the low oxygen response was gradually lost within 3 days. Introduced into a barley-soil microcosm, the biosensor could report decreasing oxygen concentrations in the rhizosphere for a 6-day period. When the water content in the microcosm was raised from 60% to 85% of field capacity, expression of the reporter gene was elevated about twofold above a basal level after 2 days of incubation, suggesting that a water content of 85% caused mild anoxia. Increased compaction of the soil was shown to have a faster and more dramatic effect on the expression of the oxygen reporter than soil water content alone, indicating that factors other than the water-filled pore space influenced the oxygen status of the soil. These experiments illustrate the utility of the biosensor for detecting low oxygen concentrations in the rhizosphere and other soil habitats.

Predominance of Nonculturable Cells of the Biocontrol Strain Pseudomonas fluorescens CHA0 in the Surface Horizon of Large Outdoor Lysimeters.

The persistence of the biocontrol agent Pseudomonas fluorescens CHA0 in the surface horizon of 12 large outdoor lysimeters planted with winter wheat, Phacelia tanacetifolia followed by spring wheat, or maize was monitored for 1 year. Soil was inoculated with a spontaneous rifampin-resistant mutant (CHA0-Rif) of CHA0, and the strain was studied by using colony counts, Kogure's direct viable counts, and total counts (immunofluorescence). The number of culturable cells of the inoculant decreased progressively from 8 to 2 log CFU/g of soil or lower. However, culturable cells of CHA0-Rif accounted for less than 1% of the total cells of the inoculant 8 months after release in autumn. Since viable but nonculturable cells represented less than a quarter of the latter, most cells of CHA0-Rif in soil were thus inactive-dormant or dead at that time. Nonculturable cells of the inoculant were predominant also in the surface horizon of the lysimeters inoculated in the spring, and a significant fraction of them were viable. Results suggest that the occurrence of nonculturable cells of CHA0-Rif was influenced by climatic factors (water availability and soil temperature) and the abundance of roots in soil. The fact that the inoculant persisted as mixed populations of cells of different physiological states, in which nonculturable cells were predominant, needs to be taken into account when assessing the autecology of wild-type or genetically modified pseudomonads released into the soil ecosystem.

The Complete Chromosomal Organization of the Reference Strain of the Leishmania Genome Project, L. major ;Friedlin'.

In recent years, analysis of the genomes of many organisms has received increasing international attention. The bulk of the effort to date has centred on the Human Genome Project and analysis of model organisms such as yeast, Drosophila and Caenorhabditis elegans. More recently, the revolution in genome sequencing and gene identification has begun to impact on infectious disease organisms. Initially, much of the effort was concentrated on prokaryotes, but small eukaryotic genomes, including the protozoan parasites Plasmodium, Toxoplasma and trypanosomatids (Leishmania, Trypanosoma brucei and T. cruzi), as well as some multicellular organisms, such as Brugia and Schistosoma, are benefiting from the technological advances of the genome era. These advances promise a radical new approach to the development of novel diagnostic tools, chemotherapeutic targets and vaccines for infectious disease organisms, as well as to the more detailed analysis of cell biology and function.Several networks or consortia linking laboratories around the world have been established to support these parasite genome projects[1] (for more information, see http://www.ebi.ac.uk/ parasites/paratable.html). Five of these networks were supported by an initiative launched in 1994 by the Specific Programme for Research and Tropical Diseases (TDR) of the WHO[2, 3, 4, 5, 6]. The Leishmania Genome Network (LGN) is one of these[3]. Its activities are reported at http://www.ebi.ac.uk/parasites/leish.html, and its current aim is to map and sequence the genome of Leishmania by the year 2002. All the mapping, hybridization and sequence data are also publicly available from LeishDB, an AceDB-based genome database (http://www.ebi.ac.uk/parasites/LGN/leissssoft.html).

Development and experimental validation of a finite element model of total ankle replacement.

Total ankle replacement remains a less satisfactory solution compared to other joint replacements. The goal of this study was to develop and validate a finite element model of total ankle replacement, for future testing of hypotheses related to clinical issues. To validate the finite element model, an experimental setup was specifically developed and applied on 8 cadaveric tibias. A non-cemented press fit tibial component of a mobile bearing prosthesis was inserted into the tibias. Two extreme anterior and posterior positions of the mobile bearing insert were considered, as well as a centered one. An axial force of 2kN was applied for each insert position. Strains were measured on the bone surface using digital image correlation. Tibias were CT scanned before implantation, after implantation, and after mechanical tests and removal of the prosthesis. The finite element model replicated the experimental setup. The first CT was used to build the geometry and evaluate the mechanical properties of the tibias. The second CT was used to set the implant position. The third CT was used to assess the bone-implant interface conditions. The coefficient of determination (R-squared) between the measured and predicted strains was 0.91. Predicted bone strains were maximal around the implant keel, especially at the anterior and posterior ends. The finite element model presented here is validated for future tests using more physiological loading conditions.

Thermomechanical analysis of strain localization in a ductile detachment zone

At the latitude of the Thor-Odin dome (British Columbia) the Columbia River Detachment defines the eastern margin of the Shuswap metamorphic core complex and localizes in a 1 km thick muscovite-bearing quartzite mylonite. We present a combined Ar-40/Ar-39, (micro) structural, and oxygen isotope study of the deformation history in the detachment and evaluate the spatial and temporal relationships between microstructure formation and localization of strain. High-precision Ar-40/Ar-39 geochronology from different levels in the mylonite delineates a pattern of increasingly younger (49.0 to 47.9 Ma) deformation ages in deeper levels of the mylonitic footwall. The correlation of Ar-40/Ar-39 ages with decreasing deformation temperatures (similar to 550 degrees-400 degrees C) in the top 200 m of the mylonite indicates that deformation migrated downward from the contact with the hanging wall. Strain localization was diachronous in progressively deeper levels of the footwall and was likely controlled by fluid-assisted strain hardening due to advective heat removal and contemporaneous reaction weakening due to dissolution-reprecipitation of white mica. The observed constant high-stress microstructures across the entire detachment indicate that flow stress was buffered by the interplay of strain rate and temperature, where high strain rates at elevated temperature produced the same microstructure as lower strain rates under decreasing temperature conditions. The combined data suggest that the complex interplay among temporally nonuniform rates of footwall exhumation, heat advection, and embrittlement by meteoric fluids strongly determines the thermomechanical behavior of extensional detachments.