From qualitative data to quantitative models: analysis of the phage shock protein stress response in Escherichia coli

Background Bacteria have evolved a rich set of mechanisms for sensing and adapting to adverse conditions in their environment. These are crucial for their survival, which requires them to react to extracellular stresses such as heat shock, ethanol treatment or phage infection. Here we focus on studying the phage shock protein (Psp) stress response in Escherichia coli induced by a phage infection or other damage to the bacterial membrane. This system has not yet been theoretically modelled or analysed in silico. Results We develop a model of the Psp response system, and illustrate how such models can be constructed and analyzed in light of available sparse and qualitative information in order to generate novel biological hypotheses about their dynamical behaviour. We analyze this model using tools from Petri-net theory and study its dynamical range that is consistent with currently available knowledge by conditioning model parameters on the available data in an approximate Bayesian computation (ABC) framework. Within this ABC approach we analyze stochastic and deterministic dynamics. This analysis allows us to identify different types of behaviour and these mechanistic insights can in turn be used to design new, more detailed and time-resolved experiments. Conclusions We have developed the first mechanistic model of the Psp response in E. coli. This model allows us to predict the possible qualitative stochastic and deterministic dynamic behaviours of key molecular players in the stress response. Our inferential approach can be applied to stress response and signalling systems more generally: in the ABC framework we can condition mathematical models on qualitative data in order to delimit e.g. parameter ranges or the qualitative system dynamics in light of available end-point or qualitative information.Medical Research Council (Great Britain)Biotechnology and Biological Sciences Research Council (Great Britain)Wellcome Trust (London, England)Royal Society (Great Britain) (Wolfson Research Merit Award

Retroviruses integrate into a shared, non-palindromic DNA motif.

Many DNA-binding factors, such as transcription factors, form oligomeric complexes with structural symmetry that bind to palindromic DNA sequences1. Palindromic consensus nucleotide sequences are also found at the genomic integration sites of retroviruses2-6 and other transposable elements7-9, and it has been suggested that this palindromic consensus arises as a consequence of the structural symmetry in the integrase complex2,3. However, we show here that the palindromic consensus sequence is not present in individual integration sites of human T-cell lymphotropic virus type 1 (HTLV-1) and human immunodeficiency virus type 1 (HIV-1), but arises in the population average as a consequence of the existence of a non-palindromic nucleotide motif that occurs in approximately equal proportions on the plus strand and the minus strand of the host genome. We develop a generally applicable algorithm to sort the individual integration site sequences into plus-strand and minus-strand subpopulations, and use this to identify the integration site nucleotide motifs of five retroviruses of different genera: HTLV-1, HIV-1, murine leukaemia virus (MLV), avian sarcoma leucosis virus (ASLV) and prototype foamy virus (PFV). The results reveal a non-palindromic motif that is shared between these retroviruses

P38 and JNK have opposing effects on persistence of in vivo leukocyte migration in zebrafish

The recruitment and migration of macrophages and neutrophils is an important process during the early stages of the innate immune system in response to acute injury. Transgenic pu.1:EGFP zebrafish permit the acquisition of leukocyte migration trajectories during inflammation. Currently, these high-quality live-imaging data are mainly analysed using general statistics, for example, cell velocity. Here, we present a spatio-temporal analysis of the cell dynamics using transition matrices, which provide information of the type of cell migration. We find evidence that leukocytes exhibit types of migratory behaviour, which differ from previously described random walk processes. Dimethyl sulfoxide treatment decreased the level of persistence at early time points after wounding and ablated temporal dependencies observed in untreated embryos. We then use pharmacological inhibition of p38 and c-Jun N-terminal kinase mitogen-activated protein kinases to determine their effects on in vivo leukocyte migration patterns and discuss how they modify the characteristics of the cell migration process. In particular, we find that their respective inhibition leads to decreased and increased levels of persistent motion in leukocytes following wounding. This example shows the high level of information content, which can be gained from live-imaging data if appropriate statistical tools are used

Rôle de la réplication dans l'évolution et l'organisation du génome humain

PARIS7-Bibliothèque centrale (751132105) / SudocSudocFranceF

Human gene organization driven by the coordination of replication and transcription

In this work, we investigated a large-scale organization of the human genes with respect to putative replication origins. We developed an appropriate multiscale method to analyze the nucleotide compositional skew along the genome and found that in more than one-quarter of the genome, the skew profile presents characteristic patterns consisting of successions of N-shaped structures, designated here N-domains, bordered by putative replication origins. Our analysis of recent experimental timing data confirmed that, in a number of cases, domain borders coincide with replication initiation zones active in the early S phase, whereas the central regions replicate in the late S phase. Around the putative origins, genes are abundant and broadly expressed, and their transcription is co-oriented with replication fork progression. These features weaken progressively with the distance from putative replication origins. At the center of domains, genes are rare and expressed in few tissues. We propose that this specific organization could result from the constraints of accommodating the replication and transcription initiation processes at chromatin level, and reducing head-on collisions between the two machineries. Our findings provide a new model of gene organization in the human genome, which integrates transcription, replication, and chromatin structure as coordinated determinants of genome architecture

Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C → T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution

Model-based evolutionary analysis: the natural history of phage-shock stress response

The evolution of proteins is inseparably linked to their function. Because most biological processes involve a number of different proteins, it may become impossible to study the evolutionary properties of proteins in isolation. In the present article, we show how simple mechanistic models of biological processes can complement conventional comparative analyses of biological traits. We use the specific example of the phage-shock stress response, which has been well characterized in Escherichia coli, to elucidate patterns of gene sharing and sequence conservation across bacterial species.</jats:p