Slow protein fluctuations explain the emergence of growth phenotypes and persistence in clonal bacterial populations

One of the most challenging problems in microbiology is to understand how a small fraction of microbes that resists killing by antibiotics can emerge in a population of genetically identical cells, the phenomenon known as persistence or drug tolerance. Its characteristic signature is the biphasic kill curve, whereby microbes exposed to a bactericidal agent are initially killed very rapidly but then much more slowly. Here we relate this problem to the more general problem of understanding the emergence of distinct growth phenotypes in clonal populations. We address the problem mathematically by adopting the framework of the phenomenon of so-called weak ergodicity breaking, well known in dynamical physical systems, which we extend to the biological context. We show analytically and by direct stochastic simulations that distinct growth phenotypes can emerge as a consequence of slow-down of stochastic fluctuations in the expression of a gene controlling growth rate. In the regime of fast gene transcription, the system is ergodic, the growth rate distribution is unimodal, and accounts for one phenotype only. In contrast, at slow transcription and fast translation, weakly non-ergodic components emerge, the population distribution of growth rates becomes bimodal, and two distinct growth phenotypes are identified. When coupled to the well-established growth rate dependence of antibiotic killing, this model describes the observed fast and slow killing phases, and reproduces much of the phenomenology of bacterial persistence. The model has major implications for efforts to develop control strategies for persistent infections.Comment: 26 pages, 7 figure

Interrogation of global mutagenesis data with a genome scale model of Neisseria meningitidis to assess gene fitness in vitro and in sera

BACKGROUND: Neisseria meningitidis is an important human commensal and pathogen that causes several thousand deaths each year, mostly in young children. How the pathogen replicates and causes disease in the host is largely unknown, particularly the role of metabolism in colonization and disease. Completed genome sequences are available for several strains but our understanding of how these data relate to phenotype remains limited. RESULTS: To investigate the metabolism of N. meningitidis we generated and selected a representative Tn5 library on rich medium, a minimal defined medium and in human serum to identify genes essential for growth under these conditions. To relate these data to a systems-wide understanding of the pathogen's biology we constructed a genome-scale metabolic network: Nmb_iTM560. This model was able to distinguish essential and non-essential genes as predicted by the global mutagenesis. These essentiality data, the library and the Nmb_iTM560 model are powerful and widely applicable resources for the study of meningococcal metabolism and physiology. We demonstrate the utility of these resources by predicting and demonstrating metabolic requirements on minimal medium such as a requirement for PEP carboxylase, and by describing the nutritional and biochemical status of N. meningitidis when grown in serum, including a requirement for both the synthesis and transport of amino acids. CONCLUSIONS: This study describes the application of a genome scale transposon library combined with an experimentally validated genome-scale metabolic network of N. meningitidis to identify essential genes and provide novel insight to the pathogen's metabolism both in vitro and during infection

Modeling and simulation of the main metabolism in Escherichia coli and its several single-gene knockout mutants with experimental verification

Peer reviewedPublisher PD

Effect of the spray volume adjustment model on the efficiency of fungicides and residues in processing tomato

This study compared the effects of a proportionate spray volume (PSV) adjustment model and a fixed model (300 L/ha) on the infestation of processing tomato with potato late blight (Phytophthora infestans (Mont.) de Bary) (PLB) and azoxystrobin and chlorothalonil residues in fruits in three consecutive seasons. The fungicides were applied in alternating system with or without two spreader adjuvants. The proportionate spray volume adjustment model was based on the number of leaves on plants and spray volume index. The modified Quick, Easy, Cheap, Effective, Rugged, and Safe (QuEChERS) method was optimized and validated for extraction of azoxystrobin and chlorothalonil residue. Gas chromatography with a nitrogen and phosphorus detector and an electron capture detector were used for the analysis of fungicides. The results showed that higher fungicidal residues were connected with lower infestation of tomato with PLB. PSV adjustment model resulted in lower infestation of tomato than the fixed model (300 L/ha) when fungicides were applied at half the dose without adjuvants. Higher expected spray interception into the tomato canopy with the PSV system was recognized as the reasons of better control of PLB. The spreader adjuvants did not have positive effect on the biological efficacy of spray volume application systems. The results suggest that PSV adjustment model can be used to determine the spray volume for fungicide application for processing tomato crop

The Genetic Requirements for Fast and Slow Growth in Mycobacteria

Mycobacterium tuberculosis infects a third of the world's population. Primary tuberculosis involving active fast bacterial replication is often followed by asymptomatic latent tuberculosis, which is characterised by slow or non-replicating bacteria. Reactivation of the latent infection involving a switch back to active bacterial replication can lead to post-primary transmissible tuberculosis. Mycobacterial mechanisms involved in slow growth or switching growth rate provide rational targets for the development of new drugs against persistent mycobacterial infection. Using chemostat culture to control growth rate, we screened a transposon mutant library by Transposon site hybridization (TraSH) selection to define the genetic requirements for slow and fast growth of Mycobacterium bovis (BCG) and for the requirements of switching growth rate. We identified 84 genes that are exclusively required for slow growth (69 hours doubling time) and 256 genes required for switching from slow to fast growth. To validate these findings we performed experiments using individual M. tuberculosis and M. bovis BCG knock out mutants. We have demonstrated that growth rate control is a carefully orchestrated process which requires a distinct set of genes encoding several virulence determinants, gene regulators, and metabolic enzymes. The mce1 locus appears to be a component of the switch to slow growth rate, which is consistent with the proposed role in virulence of M. tuberculosis. These results suggest novel perspectives for unravelling the mechanisms involved in the switch between acute and persistent TB infections and provide a means to study aspects of this important phenomenon in vitro

Acorn: A grid computing system for constraint based modeling and visualization of the genome scale metabolic reaction networks via a web interface

Constraint-based approaches facilitate the prediction of cellular metabolic capabilities, based, in turn on predictions of the repertoire of enzymes encoded in the genome. Recently, genome annotations have been used to reconstruct genome scale metabolic reaction networks for numerous species, including Homo sapiens, which allow simulations that provide valuable insights into topics, including predictions of gene essentiality of pathogens, interpretation of genetic polymorphism in metabolic disease syndromes and suggestions for novel approaches to microbial metabolic engineering. These constraint-based simulations are being integrated with the functional genomics portals, an activity that requires efficient implementation of the constraint-based simulations in the web-based environment