The Mystery of Two Straight Lines in Bacterial Genome Statistics. Release 2007

In special coordinates (codon position--specific nucleotide frequencies) bacterial genomes form two straight lines in 9-dimensional space: one line for eubacterial genomes, another for archaeal genomes. All the 348 distinct bacterial genomes available in Genbank in April 2007, belong to these lines with high accuracy. The main challenge now is to explain the observed high accuracy. The new phenomenon of complementary symmetry for codon position--specific nucleotide frequencies is observed. The results of analysis of several codon usage models are presented. We demonstrate that the mean--field approximation, which is also known as context--free, or complete independence model, or Segre variety, can serve as a reasonable approximation to the real codon usage. The first two principal components of codon usage correlate strongly with genomic G+C content and the optimal growth temperature respectively. The variation of codon usage along the third component is related to the curvature of the mean-field approximation. First three eigenvalues in codon usage PCA explain 59.1%, 7.8% and 4.7% of variation. The eubacterial and archaeal genomes codon usage is clearly distributed along two third order curves with genomic G+C content as a parameter.Comment: Significantly extended version with new data for all the 348 distinct bacterial genomes available in Genbank in April 200

A Gibbs approach to Chargaff's second parity rule

Chargaff's second parity rule (CSPR) asserts that the frequencies of short polynucleotide chains are the same as those of the complementary reversed chains. Up to now, this hypothesis has only been observed empirically and there is currently no explanation for its presence in DNA strands. Here we argue that CSPR is a probabilistic consequence of the reverse complementarity between paired strands, because the Gibbs distribution associated with the chemical energy between the bonds satisfies CSPR. We develop a statistical test to study the validity of CSPR under the Gibbsian assumption and we apply it to a large set of bacterial genomes taken from the GenBank repository.Comment: 16 page

Hydrodynamic Coupling of Two Brownian Spheres to a Planar Surface

We describe direct imaging measurements of the collective and relative diffusion of two colloidal spheres near a flat plate. The bounding surface modifies the spheres' dynamics, even at separations of tens of radii. This behavior is captured by a stokeslet analysis of fluid flow driven by the spheres' and wall's no-slip boundary conditions. In particular, this analysis reveals surprising asymmetry in the normal modes for pair diffusion near a flat surface.Comment: 4 pages, 4 figure

Controllability on infinite-dimensional manifolds

Following the unified approach of A. Kriegl and P.W. Michor (1997) for a treatment of global analysis on a class of locally convex spaces known as convenient, we give a generalization of Rashevsky-Chow's theorem for control systems in regular connected manifolds modelled on convenient (infinite-dimensional) locally convex spaces which are not necessarily normable.Comment: 19 pages, 1 figur

Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonmous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa and L. lactis as their primary host. We introduce the concept of a `genome landscape,' which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such a GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.Comment: 9 Color Figures, 5 Tables, 53 Reference

Incorporating prior knowledge improves detection of differences in bacterial growth rate

BACKGROUND: Robust statistical detection of differences in the bacterial growth rate can be challenging, particularly when dealing with small differences or noisy data. The Bayesian approach provides a consistent framework for inferring model parameters and comparing hypotheses. The method captures the full uncertainty of parameter values, whilst making effective use of prior knowledge about a given system to improve estimation. RESULTS: We demonstrated the application of Bayesian analysis to bacterial growth curve comparison. Following extensive testing of the method, the analysis was applied to the large dataset of bacterial responses which are freely available at the web-resource, ComBase. Detection was found to be improved by using prior knowledge from clusters of previously analysed experimental results at similar environmental conditions. A comparison was also made to a more traditional statistical testing method, the F-test, and Bayesian analysis was found to perform more conclusively and to be capable of attributing significance to more subtle differences in growth rate. CONCLUSIONS: We have demonstrated that by making use of existing experimental knowledge, it is possible to significantly improve detection of differences in bacterial growth rate

Component-Based Real-Time Operating System for Embedded Applications

Acceptance rate: 37%, Rank (CORE): AInternational audienceAs embedded systems must constantly integrate new functionalities, their developement cycles must be based on high-level abstractions, making the software design more flexible. CBSE provides an approach to these new requirements. However, low-level services provided by operating systems are an integral part of embedded applications, furthermore deployed on resource-limited devices. Therefore, the expected benefits of CBSE must not impact on the constraints imposed by the targetted domain, such as memory footprint, energy consumption, and execution time. In this paper, we present the componentization of a legacy industry-established Real-Time Operating System, and how component-based applications are built on top of it. We use the Think framework that allows to produce flexible systems while paying for flexibility only where desired. Performed experimentions show that the induced overhead is negligeable

A method for the reconstruction of unknown non-monotonic growth functions in the chemostat

We propose an adaptive control law that allows one to identify unstable steady states of the open-loop system in the single-species chemostat model without the knowledge of the growth function. We then show how one can use this control law to trace out (reconstruct) the whole graph of the growth function. The process of tracing out the graph can be performed either continuously or step-wise. We present and compare both approaches. Even in the case of two species in competition, which is not directly accessible with our approach due to lack of controllability, feedback control improves identifiability of the non-dominant growth rate.Comment: expansion of ideas from proceedings paper (17 pages, 8 figures), proceedings paper is version v

Higher order glass-transition singularities in colloidal systems with attractive interactions

The transition from a liquid to a glass in colloidal suspensions of particles interacting through a hard core plus an attractive square-well potential is studied within the mode-coupling-theory framework. When the width of the attractive potential is much shorter than the hard-core diameter, a reentrant behavior of the liquid-glass line, and a glass-glass-transition line are found in the temperature-density plane of the model. For small well-width values, the glass-glass-transition line terminates in a third order bifurcation point, i.e. in a A_3 (cusp) singularity. On increasing the square-well width, the glass-glass line disappears, giving rise to a fourth order A_4 (swallow-tail) singularity at a critical well width. Close to the A_3 and A_4 singularities the decay of the density correlators shows stretching of huge dynamical windows, in particular logarithmic time dependence.Comment: 19 pages, 12 figures, Phys. Rev. E, in prin

Atypical AT Skew in Firmicute Genomes Results from Selection and Not from Mutation

The second parity rule states that, if there is no bias in mutation or selection, then within each strand of DNA complementary bases are present at approximately equal frequencies. In bacteria, however, there is commonly an excess of G (over C) and, to a lesser extent, T (over A) in the replicatory leading strand. The low G+C Firmicutes, such as Staphylococcus aureus, are unusual in displaying an excess of A over T on the leading strand. As mutation has been established as a major force in the generation of such skews across various bacterial taxa, this anomaly has been assumed to reflect unusual mutation biases in Firmicute genomes. Here we show that this is not the case and that mutation bias does not explain the atypical AT skew seen in S. aureus. First, recently arisen intergenic SNPs predict the classical replication-derived equilibrium enrichment of T relative to A, contrary to what is observed. Second, sites predicted to be under weak purifying selection display only weak AT skew. Third, AT skew is primarily associated with largely non-synonymous first and second codon sites and is seen with respect to their sense direction, not which replicating strand they lie on. The atypical AT skew we show to be a consequence of the strong bias for genes to be co-oriented with the replicating fork, coupled with the selective avoidance of both stop codons and costly amino acids, which tend to have T-rich codons. That intergenic sequence has more A than T, while at mutational equilibrium a preponderance of T is expected, points to a possible further unresolved selective source of skew