Mutation rates: estimating phase variation rates when fitness differences are present and their impact on population structure

Influence of phenotypic variation of Paenibacillus polymyxa E681 on growth promotion in cucumbers

The goal of the current study is to better understand how bacteria may adapt to survive under adverse environmental conditions by altering and improving their phenotypes. In this study, we report the consequences of phenotypic variation in Paenibacillus polymyxa E681 (E681), a plant growth-promoting rhizobacterium (PGPR), isolated from winter barley root that has a variety of advantageous effects on crop plants. In our previous study, two different types of bacterial cells in E681 were distinguished. We used the term F-type for the variant that doesn't produce endospores and B-type for the endospore-producing wild type. Under the circumstances of our experiment, the cucumber rhizosphere soil and the surface of the seeds produced phenotypic variance. On tryptic soy agar (TSA) plates, the B-type spontaneously converted into the F-type, but the reverse was not reversible. Intriguingly, the plant growth promotion test displayed that cucumber seedlings treated with F-type cells had characteristics resembling those of the untreated control. Whereas, growth promotion of cucumber seedlings treated with B-type depends on temperature conditions. In particular, an increased growth promotion was observed at a low temperature of 20°C. The phenotypic change from B-type to F-type did not occur at 20°C for 6 days in the growth curve analysis of E681, but it did occur on the fourth and second days at 30 and 37°C, respectively. Therefore, before using PGPR strains as a bacterial inoculant for sustainable agriculture, it is imperative to resolve phenotypic variance in these strains.

Single-cell analysis of multiple invertible promoters reveals differential inversion rates as a strong determinant of bacterial population heterogeneity

Population heterogeneity can promote bacterial fitness in response to unpredictable environmental conditions. A major mechanism of phenotypic variability in the human gut symbiont Bacteroides spp. involves the inversion of promoters that drive the expression of capsular polysaccharides, which determine the architecture of the cell surface. High-throughput single-cell sequencing reveals substantial population heterogeneity generated through combinatorial promoter inversion regulated by a broadly conserved serine recombinase. Exploiting control over population diversification, we show that populations with different initial compositions converge to a similar composition over time. Combining our data with stochastic computational modeling, we demonstrate that the differential rates of promoter inversion are a major mechanism shaping population dynamics. More broadly, our approach could be used to interrogate single-cell combinatorial phase variable states of diverse microbes including bacterial pathogens.

Adaption of Pseudomonas ogarae F113 to the Rhizosphere Environment-The AmrZ-FleQ Hub

Motility and biofilm formation are two crucial traits in the process of rhizosphere colonization by pseudomonads. The regulation of both traits requires a complex signaling network that is coordinated by the AmrZ-FleQ hub. In this review, we describe the role of this hub in the adaption to the rhizosphere. The study of the direct regulon of AmrZ and the phenotypic analyses of an amrZ mutant in Pseudomonas ogarae F113 has shown that this protein plays a crucial role in the regulation of several cellular functions, including motility, biofilm formation, iron homeostasis, and bis-(3'-5')-cyclic dimeric guanosine monophosphate (c-di-GMP) turnover, controlling the synthesis of extracellular matrix components. On the other hand, FleQ is the master regulator of flagellar synthesis in P. ogarae F113 and other pseudomonads, but its implication in the regulation of multiple traits related with environmental adaption has been shown. Genomic scale studies (ChIP-Seq and RNA-Seq) have shown that in P. ogarae F113, AmrZ and FleQ are general transcription factors that regulate multiple traits. It has also been shown that there is a common regulon shared by the two transcription factors. Moreover, these studies have shown that AmrZ and FleQ form a regulatory hub that inversely regulate traits such as motility, extracellular matrix component production, and iron homeostasis. The messenger molecule c-di-GMP plays an essential role in this hub since its production is regulated by AmrZ and it is sensed by FleQ and required for its regulatory role. This regulatory hub is functional both in culture and in the rhizosphere, indicating that the AmrZ-FleQ hub is a main player of P. ogarae F113 adaption to the rhizosphere environment.

Effects of Phenotypic Variation on Biological Properties of Endophytic Bacteria Bacillus mojavensis PS17

Simple Summary

Microorganisms play an important role in agriculture by protecting and stimulating the growth of plants. The phenotypic activities of microbial biological agents (MBA) can change under different environmental conditions. However, to adapt to these harsh conditions, genetic mutations take place in bacteria that are seen phenotypically, which might not be beneficial or less beneficial to the plants. Some adaptative mechanisms used by microorganisms, especially bacteria, to face these environmental factors lead to the appearance of subpopulations with different morphotypes that may be more adapted to survive in stressful conditions. Moreover, in favorable conditions, these subpopulations may become dominant among the overall bacterial population. In this study, Bacillus mojavensis undergoes phase variation when grown in a minimal medium, in which two colonies, opaque (morphotype I) and translucent (morphotype II), were generated. The characteristics of the generated morphotypes were determined and compared with those of their original strain. Overall, the results obtained showed that the phenotypic characteristics of morphotype I statistically differed from morphotype II. This phenomenon may be one of the factors behind the dissimilarities in the results between the laboratory and field data on the application of MBA.

Abstract

The use of microorganism-based products in agricultural practices is gaining more interest as an alternative to chemical methods due to their non-toxic bactericidal and fungicidal properties. Various factors influence the efficacy of the microorganisms used as biological control agents in infield conditions as compared to laboratory conditions due to ecological and physiological aspects. Abiotic factors have been shown to trigger phase variations in bacterial microorganisms as a mechanism for adapting to hostile environments. In this study, we investigated the stability of the morphotype and the effects of phenotypic variation on the biological properties of Bacillus mojavensis strain PS17. B. mojavensis PS17 generated two variants (opaque and translucent) that were given the names morphotype I and II, respectively. The partial sequence of the 16S rRNA gene revealed that both morphotypes belonged to B. mojavensis. BOX and ERIC fingerprinting PCR also showed the same DNA profiles in both morphotypes. The characteristics of morphotype I did not differ from the original strain, while morphotype II showed a lower hydrolytic enzyme activity, phytohormone production, and antagonistic ability against phytopathogenic fungi. Both morphotypes demonstrated endophytic ability in tomato plants. A low growth rate of the strain PS17(II) in a minimal medium was observed in comparison to the PS17(I) strain. Furthermore, the capacity for biocontrol of B. mojavensis PS17(II) was not effective in the suppression of root rot disease in the tomato plants caused by Fusarium oxysporum f. sp. radices-lycopersici stain ZUM2407, compared to B. mojavensis PS17(I), whose inhibition was almost 47.9 ± 1.03% effective.

Mutation and Selection in Bacteria: Modelling and Calibration

Temporal evolution of a clonal bacterial population is modelled taking into account reversible mutation and selection mechanisms. For the mutation model, an efficient algorithm is proposed to verify whether experimental data can be explained by this model. The selection–mutation model has unobservable fitness parameters, and, to estimate them, we use an Approximate Bayesian Computation algorithm. The algorithms are illustrated using in vitro data for phase variable genes of Campylobacter jejuni.

Phase variation of Opa proteins of Neisseria meningitidis and the effects of bacterial transformation

Opa proteins are major proteins involved in meningococcal colonization of the nasopharynx and immune interactions. Opa proteins undergo phase variation (PV) due to the presence of the 5'-CTCTT-3' coding repeat (CR) sequence. The dynamics of PV of meningococcal Opa proteins is unknown. Opa PV, including the effect of transformation on PV, was assessed using a panel of Opa-deficient strains of Neisseria meningitidis. Analysis of Opa expression from UK disease-causing isolates was undertaken. Different opa genes demonstrated variable rates of PV, between 6.4 × 10(-4) and 6.9 × 10(-3) per cell per generation. opa genes with a longer CR tract had a higher rate of PV (r(2) = 0.77, p = 0.1212). Bacterial transformation resulted in a 180-fold increase in PV rate. The majority of opa genes in UK disease isolates (315/463, 68.0%) were in the 'on' phase, suggesting the importance of Opa proteins during invasive disease. These data provide valuable information for the first time regarding meningococcal Opa PV. The presence of Opa PV in meningococcal populations and high expression of Opa among invasive strains likely indicates the importance of this protein in bacterial colonization in the human nasopharynx. These findings have potential implications for development of vaccines derived from meningococcal outer membranes.

Mutability in Pseudomonas viridiflava as a programmed balance between antibiotic resistance and pathogenicity

Mutable bacterial cells are defective in their DNA repair system and often have a phenotype different from that of their wild-type counterparts. In human bacterial pathogens, the mutable and hypermutable phenotypes are often associated with general antibiotic resistance. Here, we quantified the occurrence of mutable cells in Pseudomonas viridiflava, a phytopathogenic bacterium in the P. syringae complex with a broad host range and capacity to live as a saprophyte. Two phenotypic variants (transparent and mucoid) were produced by this bacterium. The transparent variant had a mutator phenotype, showed general antibiotic resistance and could not induce disease on the plant species tested (bean). In contrast, the mucoid variant did not display mutability or resistance to antibiotics and was capable of inducing disease on bean. Both the transparent and mucoid variants were less fit when grown in vitro, whereas, in planta, both of the variants and wild-types attained similar population densities. Given the importance of the methyl-directed mismatch repair system (MMR) in the occurrence of mutable and hypermutable cells in human bacterial pathogens, we investigated whether mutations in mut genes were associated with mutator transparent cells in P. viridiflava. Our results showed no mutations in MMR genes in any of the P. viridiflava cells tested. Here, we report that a high mutation rate and antibiotic resistance are inversely correlated with pathogenicity in P. viridiflava, but are not associated with mutations in MMR. In addition, P. viridiflava variants differ from variants produced by other phytopathogenic bacteria in the absence of reversion to the wild-type phenotype.

Stochastic Switching of Cell Fate in Microbes

Microbes transiently differentiate into distinct, specialized cell types to generate functional diversity and cope with changing environmental conditions. Though alternate programs often entail radically different physiological and morphological states, recent single-cell studies have revealed that these crucial decisions are often left to chance. In these cases, the underlying genetic circuits leverage the intrinsic stochasticity of intracellular chemistry to drive transition between states. Understanding how these circuits transform transient gene expression fluctuations into lasting phenotypic programs will require a combination of quantitative modeling and extensive, time-resolved observation of switching events in single cells. In this article, we survey microbial cell fate decisions demonstrated to involve a random element, describe theoretical frameworks for understanding stochastic switching between states, and highlight recent advances in microfluidics that will enable characterization of key dynamic features of these circuits.

A versatile, non genetically modified organism (GMO)-based strategy for controlling low-producer mutants in Bordetella pertussis cultures using antigenic modulation

The uncontrolled presence of non-producer mutants negatively affects bioprocesses. In Bordetella pertussis cultures, avirulent mutants emerge spontaneously and accumulate. We characterized the dynamics of accumulation using high-throughput growth assays and competition experiments between virulent and avirulent (bvg(-) ) isolates. A fitness advantage of bvg(-) cells was identified as the main driver for bvg(-) accumulation under conditions of high virulence factor production. Conversely, under conditions that reduce their expression (antigenic modulation), bvg(-) takeover could be avoided. A control strategy was derived, which consists in applying modulating conditions whenever virulence factor production is not required. It has a wide range of applications, from routine laboratory operations to vaccine manufacturing, where pertussis toxin yields were increased 1.4-fold by performing early pre-culture steps in modulating conditions. Because it only requires subtle modifications of the culture medium and does not involve genetic modifications, this strategy is applicable to any B. pertussis isolate, and should facilitate regulatory acceptance of process changes for vaccine production. Strategies based on the same concept, could be derived for other industrially relevant micro-organisms. This study illustrates how a sound scientific understanding of physiological principles can be turned into a practical application for the bioprocess industry, in alignment with Quality by Design principles.

Automatic validation of computational models using pseudo-3D spatio-temporal model checking

BACKGROUND

Computational models play an increasingly important role in systems biology for generating predictions and in synthetic biology as executable prototypes/designs. For real life (clinical) applications there is a need to scale up and build more complex spatio-temporal multiscale models; these could enable investigating how changes at small scales reflect at large scales and viceversa. Results generated by computational models can be applied to real life applications only if the models have been validated first. Traditional in silico model checking techniques only capture how non-dimensional properties (e.g. concentrations) evolve over time and are suitable for small scale systems (e.g. metabolic pathways). The validation of larger scale systems (e.g. multicellular populations) additionally requires capturing how spatial patterns and their properties change over time, which are not considered by traditional non-spatial approaches.

RESULTS

We developed and implemented a methodology for the automatic validation of computational models with respect to both their spatial and temporal properties. Stochastic biological systems are represented by abstract models which assume a linear structure of time and a pseudo-3D representation of space (2D space plus a density measure). Time series data generated by such models is provided as input to parameterised image processing modules which automatically detect and analyse spatial patterns (e.g. cell) and clusters of such patterns (e.g. cellular population). For capturing how spatial and numeric properties change over time the Probabilistic Bounded Linear Spatial Temporal Logic is introduced. Given a collection of time series data and a formal spatio-temporal specification the model checker Mudi ( http://mudi.modelchecking.org ) determines probabilistically if the formal specification holds for the computational model or not. Mudi is an approximate probabilistic model checking platform which enables users to choose between frequentist and Bayesian, estimate and statistical hypothesis testing based validation approaches. We illustrate the expressivity and efficiency of our approach based on two biological case studies namely phase variation patterning in bacterial colony growth and the chemotactic aggregation of cells.

CONCLUSIONS

The formal methodology implemented in Mudi enables the validation of computational models against spatio-temporal logic properties and is a precursor to the development and validation of more complex multidimensional and multiscale models.

Nonequilibrium population dynamics of phenotype conversion of cancer cells

Tumorigenesis is a dynamic biological process that involves distinct cancer cell subpopulations proliferating at different rates and interconverting between them. In this paper we proposed a mathematical framework of population dynamics that considers both distinctive growth rates and intercellular transitions between cancer cell populations. Our mathematical framework showed that both growth and transition influence the ratio of cancer cell subpopulations but the latter is more significant. We derived the condition that different cancer cell types can maintain distinctive subpopulations and we also explain why there always exists a stable fixed ratio after cell sorting based on putative surface markers. The cell fraction ratio can be shifted by changing either the growth rates of the subpopulations (Darwinism selection) or by environment-instructed transitions (Lamarckism induction). This insight can help us to understand the dynamics of the heterogeneity of cancer cells and lead us to new strategies to overcome cancer drug resistance.

Systematic identification and quantification of phase variation in commensal and pathogenic Escherichia coli

Bacteria have been shown to generate constant genetic variation in a process termed phase variation. We present a tool based on whole genome sequencing that allows detection and quantification of coexisting genotypes mediated by genomic inversions in bacterial cultures. We tested our method on widely used strains of Escherichia coli, and detected stable and reproducible phase variation in several invertible loci. These are shown here to be responsible for maintaining constant variation in populations grown from a single colony. Applying this tool on other bacterial strains can shed light on how pathogens adjust to hostile environments by diversifying their genomes.

Mathematical and live meningococcal models for simple sequence repeat dynamics - coherent predictions and observations

Evolvability by means of simple sequence repeat (SSR) instability is a feature under the constant influence of opposing selective pressures to expand and compress the repeat tract and is mechanistically influenced by factors that affect genetic instability. In addition to direct selection for protein expression and structural integrity, other factors that influence tract length evolution were studied. The genetic instability of SSRs that switch the expression of antibiotic resistance ON and OFF was modelled mathematically and monitored in a panel of live meningococcal strains. The mathematical model showed that the SSR length of a theoretical locus in an evolving population may be shaped by direct selection of expression status (ON or OFF), tract length dependent (α) and tract length independent factors (β). According to the model an increase in α drives the evolution towards shorter tracts. An increase in β drives the evolution towards a normal distribution of tract lengths given that an upper and a lower limit are set. Insertion and deletion biases were shown to skew allelic distributions in both directions. The meningococcal SSR model was tested in vivo by monitoring the frequency of spectinomycin resistance OFF→ON switching in a designed locus. The instability of a comprehensive panel of the homopolymeric SSRs, constituted of a range of 5-13 guanine nucleotides, was monitored in wildtype and mismatch repair deficient backgrounds. Both the repeat length itself and mismatch repair deficiency were shown to influence the genetic instability of the homopolymeric tracts. A possible insertion bias was observed in tracts ≤G10. Finally, an inverse correlation between the number of tract-encoded amino acids and growth in the presence of ON-selection illustrated a limitation to SSR expansion in an essential gene associated with the designed model locus and the protein function mediating antibiotic resistance.

Broad conditions favor the evolution of phase-variable loci

Simple sequence repeat (SSR) tracts produce stochastic on-off switching, or phase variation, in the expression of a panoply of surface molecules in many bacterial commensals and pathogens. A change to the number of repeats in a tract may alter the phase of the translational reading frame, which toggles the on-off state of the switch. Here, we construct an in silico SSR locus with mutational dynamics calibrated to those of the Haemophilus influenzae mod locus. We simulate its evolution in a regimen of two alternating environments, simultaneously varying the selection coefficient, s, and the epoch length, T. Some recent work in a simpler (two-locus) model suggested that stochastic switching in a regimen of two alternating environments may be evolutionarily favored only if the selection coefficients in the two environments are nearly equal (“symmetric”) or selection is very strong. This finding was puzzling, as it greatly restricted the conditions under which stochastic switching might evolve. Instead, we find agreement with other recent theoretical work, observing selective utility for stochastic switching if the product sT is large enough for the favored state to nearly fix in both environments. Symmetry is required neither in s nor in sT. Because we simulate finite populations and use a detailed model of the SSR locus, we are also able to examine the impact of population size and of several SSR locus parameters. Our results indicate that conditions favoring evolution and maintenance of SSR loci in bacteria are quite broad.

Bacteria experience frequent changes of environment during the infection cycle. One means to rapidly adapt is stochastic switching: a bacterial lineage will stochastically produce a variety of genotypes, so that some descendants will survive if the environment changes. Stochastic switching mediated by simple sequence repeat (SSR) loci is widespread among bacterial commensals and pathogens and influences critical interactions with host surfaces or immune effectors, thereby affecting host persistence, transmission, and virulence. Here, we use the most detailed in silico model of an SSR locus to date, with its phase variation calibrated to match the mod locus of Haemophilus influenzae. The type III restriction-modification system encoded by mod participates in the regulation of multiple other genes; thus, SSR-mediated phase variation of mod has far-reaching cis-regulatory effects. This coupling of phase-variable switching to complex phenotypic effects has been described as the “phasevarion” and is central to understanding the infection cycle of bacterial commensals and pathogens.

Phenotypic plasticity stimulated by cooperation fosters pattern diversity of bacterial colonies

Colonies of flagellated bacteria on agar plates are known to take on diverse morphologies. A diffusion-reaction model is proposed for bacterial-colony pattern formation on a surface due to time scale separation between the slow mass migration of bacteria from the point of inoculation, and the fast, but localized, dynamics of bacterial phenotypic plasticity stimulated by public-goods cooperation and phenotypic switching. By considering two switchable phenotypes in the population, the model generates pattern diversity typifying those reported by experimental studies.

Phenotypic variation in Acidovorax radicisN35 influences plant growth promotion

Acidovorax radicis N35, isolated from surface-sterilized wheat roots (Triticum aestivum), showed irreversible phenotypic variation in nutrient broth, resulting in a differing colony morphology. In addition to the wild-type form (rough colony type), a phenotypic variant form (smooth colony type) appeared at a frequency of 3.2 × 10(-3) per cell per generation on NB agar plates. In contrast to the N35 wild type, the variant N35v showed almost no cell aggregation and had lost its flagella and swarming ability. After inoculation, only the wild-type N35 significantly promoted the growth of soil-grown barley plants. After co-inoculation of axenically grown barley seedlings with differentially fluorescently labeled N35 and N35v cells, decreased competitive endophytic root colonization in the phenotypic variant N35v was observed using confocal laser scanning microscopy. In addition, 454 pyrosequencing of both phenotypes revealed almost identical genomic sequences. The only stable difference noted in the sequence of the phenotype variant N35v was a 16-nucleotide deletion identified in a gene encoding the mismatch repair protein MutL. The deletion resulted in a frameshift that revealed a new stop codon resulting in a truncated MutL protein missing a functional MutL C-terminal domain. The mutation was consistent in all investigated phenotype variant cultures and might be responsible for the observed phenotypic variation in A. radicis N35.

Plant growth promoting rhizobacteria (PGPR): the bugs to debug the root zone

Interaction of plant growth promoting rhizobacteria (PGPR) with host plants is an intricate and interdependent relationship involving not only the two partners but other biotic and abiotic factors of the rhizosphere region. Survival and establishment of PGPR in the rhizosphere is a major concern of agricultural microbiologists. Various factors that play a determining role include the composition of root exudates, properties of bacterial strain, soil status, and activities of other soil microbes. This review focuses on the different components that affect root colonization of PGPR and the underlying principles behind the success of these bugs to tide over the unfavorable conditions.

Individual-level bet hedging in the bacterium Sinorhizobium meliloti

The expression of phenotypic variability can enhance geometric mean fitness and act as a bet-hedging strategy in unpredictable environments. Metazoan bet hedging usually involves phenotypic diversification among an individual's offspring, such as differences in seed dormancy. Virtually all known microbial bet-hedging strategies, in contrast, rely on low-probability stochastic switching of a heritable phenotype by individual cells in a clonal group. This is less effective at generating within-group diversity when group size is small. Here we describe a novel microbial bet-hedging behavior that resembles individual-level metazoan bet hedging. Sinorhizobium meliloti stores carbon and energy in poly-3-hydroxybutyrate (PHB) as a contingency against carbon scarcity. We show that, when starved, dividing S. meliloti bet hedge by forming two daughter cells with different phenotypes. These have high and low PHB levels and are suited to long- and short-term starvation, respectively. The low-PHB cells have greater competitiveness for resources, whereas the high-PHB cells can survive for over a year without food, perhaps until a legume host is next available.

The spontaneous appearance rate of the yeast prion [PSI+] and its implications for the evolution of the evolvability properties of the [PSI+] system

Epigenetically inherited aggregates of the yeast prion [PSI+] cause genomewide readthrough translation that sometimes increases evolvability in certain harsh environments. The effects of natural selection on modifiers of [PSI+] appearance have been the subject of much debate. It seems likely that [PSI+] would be at least mildly deleterious in most environments, but this may be counteracted by its evolvability properties on rare occasions. Indirect selection on modifiers of [PSI+] is predicted to depend primarily on the spontaneous [PSI+] appearance rate, but this critical parameter has not previously been adequately measured. Here we measure this epimutation rate accurately and precisely as 5.8 x 10(-7) per generation, using a fluctuation test. We also determine that genetic "mimics" of [PSI+] account for up to 80% of all phenotypes involving general nonsense suppression. Using previously developed mathematical models, we can now infer that even in the absence of opportunities for adaptation, modifiers of [PSI+] are only weakly deleterious relative to genetic drift. If we assume that the spontaneous [PSI+] appearance rate is at its evolutionary optimum, then opportunities for adaptation are inferred to be rare, such that the [PSI+] system is favored only very weakly overall. But when we account for the observed increase in the [PSI+] appearance rate in response to stress, we infer much higher overall selection in favor of [PSI+] modifiers, suggesting that [PSI+]-forming ability may be a consequence of selection for evolvability.

Determinants of phase variation rate and the fitness implications of differing rates for bacterial pathogens and commensals

Phase variation (PV) of surface molecules and other phenotypes is a major adaptive strategy of pathogenic and commensal bacteria. Phase variants are produced at high frequencies and in a reversible manner by hypermutation or hypervariable methylation in specific regions of the genome. The major mechanisms of PV involve site-specific recombination, homologous recombination, simple sequence DNA repeat tracts or epigenetic modification by the dam methylase. PV rates of some of these mechanisms are subject to the influence of genome maintenance pathways such as DNA replication, recombination and repair while others are independent of these pathways. For each of these mechanisms, the rate of generation of phase variants is controlled by intrinsic and dispensable factors. These factors can impart environmental regulation on switching rates while many factors are subject to heterogeneity both within isolates of a species and between species. A major gap in our understanding is whether these environmental and epidemiological variations in PV rate have a major impact on fitness. Experimental approaches to studying the biological relevance of differing PV rates are being developed, and a recent intriguing finding is of a co-ordination of switching rates in the phase variable P-pili of uropathogenic bacteria.

Exploration of intraclonal adaptation mechanisms of Pseudomonas brassicacearum facing cadmium toxicity

Pseudomonas brassicacearum forms phenotypic variants in vitro as well as in planta during root colonization under natural conditions, leading to subpopulations (phase I and II cells) that differ in colony morphology and production of exoenzymes/secondary metabolites. The maximal concentration of cadmium allowing both variants growth was 25 μM; however, phase II cells accumulated fivefold higher Cd than phase I cells, even though both variants showed the same growth rate and kinetics, comprising a long stasis period (50 h). The whole transcriptome analysis of both variants in response to Cd was investigated using the home-made DNA microarrays. This analysis revealed completely different adaptation mechanisms developed by each variant to withstand and grow in the presence of the toxic. A re-organization of the cell wall to limit Cd entrance was noticed for phase I cells, as genes encoding levan exopolymers were downregulated at the expense of an upregulation of genes encoding alginate, and an upregulation of transporters such as cadA, and a downregulation of copper transporters. Phase II cells were unable to prevent Cd entrance and recruited genes under the control of oxyR and soxR regulation to face osmotic and oxidant stresses generated by Cd. Putrescine and spermidine metabolism appeared to play a central role in Cd tolerance. Microarray data were validated by biological analyses such as motility, oxidative stress assay, metabolite profiling with ICR-FT/MS and UPLC, capillary electrophoresis analysis of biogenic amines.

Campylobacter jejuni colonization and transmission in broiler chickens: a modelling perspective

Campylobacter jejuni is one of the most common causes of acute enteritis in the developed world. The consumption of contaminated poultry, where C. jejuni is believed to be a commensal organism, is a major risk factor. However, the dynamics of this colonization process in commercially reared chickens is still poorly understood. Quantification of these dynamics of infection at an individual level is vital to understand transmission within populations and formulate new control strategies. There are multiple potential routes of introduction of C. jejuni into a commercial flock. Introduction is followed by a rapid increase in environmental levels of C. jejuni and the level of colonization of individual broilers. Recent experimental and epidemiological evidence suggest that the celerity of this process could be masking a complex pattern of colonization and extinction of bacterial strains within individual hosts. Despite the rapidity of colonization, experimental transmission studies exhibit a highly variable and unexplained delay time in the initial stages of the process. We review past models of transmission of C. jejuni in broilers and consider simple modifications, motivated by the plausible biological mechanisms of clearance and latency, which could account for this delay. We show how simple mathematical models can be used to guide the focus of experimental studies by providing testable predictions based on our hypotheses. We conclude by suggesting that competition experiments could be used to further understand the dynamics and mechanisms underlying the colonization process. The population models for such competition processes have been extensively studied in other ecological and evolutionary contexts. However, C. jejuni can potentially adapt phenotypically through phase variation in gene expression, leading to unification of ecological and evolutionary time-scales. For a theoretician, the colonization dynamics of C. jejuni offer an experimental system to explore these 'phylodynamics', the synthesis of population dynamics and evolutionary biology.

Involvement of a quorum-sensing-regulated lipase secreted by a clinical isolate of Burkholderia glumae in severe disease symptoms in rice

Burkholderia glumae is an emerging rice pathogen in several areas around the world. Closely related Burkholderia species are important opportunistic human pathogens for specific groups of patients, such as patients with cystic fibrosis and patients with chronic granulomatous disease. Here we report that the first clinical isolate of B. glumae, strain AU6208, has retained its capability to be very pathogenic to rice. As previously reported for rice isolate B. glumae BGR1 (and also for the clinical isolate AU6208), TofI or TofR acyl homoserine lactone (AHL) quorum sensing played a pivotal role in rice virulence. We report that AHL quorum sensing in B. glumae AU6208 regulates secreted LipA lipase and toxoflavin, the phytotoxin produced by B. glumae. B. glumae AU6208 lipA mutants were no longer pathogenic to rice, indicating that the lipase is an important virulence factor. It was also established that type strain B. glumae ATCC 33617 did not produce toxoflavin and lipase and was nonpathogenic to rice. It was determined that in strain ATCC 33617 the LuxR family quorum-sensing sensor/regulator TofR was inactive. Introducing the tofR gene of B. glumae AU6208 in strain ATCC 33617 restored its ability to produce toxoflavin and the LipA lipase. This study extends the role of AHL quorum sensing in rice pathogenicity through the regulation of a lipase which was demonstrated to be a virulence factor. It is the first report of a clinical B. glumae isolate retaining strong rice pathogenicity and finally determined that B. glumae can undergo phenotypic conversion through a spontaneous mutation in the tofR regulator.

Phase variation and microevolution at homopolymeric tracts in Bordetella pertussis

BACKGROUND

Bordetella pertussis, the causative agent of whooping cough, is a highly clonal pathogen of the respiratory tract. Its lack of genetic diversity, relative to many bacterial pathogens, could limit its ability to adapt to a hostile and changing host environment. This limitation might be overcome by phase variation, as observed for other mucosal pathogens. One of the most common mechanisms of phase variation is reversible expansion or contraction of homopolymeric tracts (HPTs).

RESULTS

The genomes of B. pertussis and the two closely related species, B. bronchiseptica and B. parapertussis, were screened for homopolymeric tracts longer than expected on the basis of chance, given their nucleotide compositions. Sixty-nine such HPTs were found in total among the three genomes, 74% of which were polymorphic among the three species. Nine HPTs were genotyped in a collection of 90 geographically and temporally diverse B. pertussis strains using the polymerase chain reaction/ligase detection reaction (PCR/LDR) assay. Six HPTs were polymorphic in this collection of B. pertussis strains. Of note, one of these polymorphic HPTs was found in the fimX promoter, where a single base insertion variant was present in seven strains, all of which were isolated prior to introduction of the pertussis vaccine. Transcript abundance of fimX was found to be 3.8-fold lower in strains carrying the longer allele. HPTs in three other genes, tcfA, bapC, and BP3651, varied widely in composition across the strain collection and displayed allelic polymorphism within single cultures.

CONCLUSION

Allelic polymorphism at homopolymeric tracts is common within the B. pertussis genome. Phase variability may be an important mechanism in B. pertussis for evasion of the immune system and adaptation to different niches in the human host. High sensitivity and specificity make the PCR/LDR assay a powerful tool for investigating allelic variation at HPTs. Using this method, allelic diversity and phase variation were demonstrated at several B. pertussis loci.

Bacterial contingency loci: the role of simple sequence DNA repeats in bacterial adaptation

Bacterial pathogens face stringent challenges to their survival because of the many unpredictable, often precipitate, and dynamic changes that occur in the host environment or in the process of transmission from one host to another. Bacterial adaptation to their hosts involves either a mechanism for sensing and responding to external changes or the selection of variants that arise through mutation. Here we review how bacterial pathogens exploit localized hypermutation, through polymerase slippage of simple sequence repeats (SSRs), to generate phenotypic variation and enhanced fitness. These SSRs are located within the reading frame or in the promoter of a subset of genes, often termed contingency loci, whose functions are usually involved in direct interactions with host structures.

Intestinal versus external growth conditions change the surficial properties in a collection of environmental Escherichia coli isolates

Predicting the fate of microorganisms in the environment is increasingly warranted, especially for pathogenic strains. A major habitat of Escherichia coli, which encompasses commensal as well as pathogenic strains, is the gastrointestinal tract with conditions very different from the environment it encounters after shedding from the host or during cultivation in the laboratory. We developed two relevant growth conditions representative of intestinal (host-associated) and external (postshedding) environments to investigate the surficial properties and behaviors of a diverse subset of E. coli feedlot isolates. Surficial properties may determine an isolate's physical fate. A pronounced increase in cell hydrophobicity and concomitant biofilm mass formation was observed for isolates grown under external conditions. Isolates that exhibited the highest surface hydrophobicity also formed visible suspended aggregates under external conditions. Other than hydrophobicity, flagella-mediated motility was determinant in affecting E. coli biofilm formation under external conditions, with all four nonmotile E. coli isolates characterized as thin-biofilm formers. The majority (88%) of Ag43+ (outer membrane protein, antigen 43) isolates formed thick biofilms, whereas the majority (75%) of Ag43- isolates formed thin biofilms. The tested E. coli O157:H7 strain behaved differently from the environmental E. coli isolates: it displayed a low electrostatic charge, a small decrease in hydrophobicity upon shifts to external conditions, and very little biofilm formation. On the other hand, the commonly used laboratory strain E. coli K-12 displayed low hydrophobicity both intestinally and externally, but it formed significant biofilm mass under external conditions. Clearly, various E. coli strains manifest significant variability in surficial behavior. This variability is further modulated by growth conditions. The interacting strain-inherent and cultivation-dependent effects on surficial behavior may have broad consequences for the fate and ecology of pathogenic and commensal E. coli strains.

The role of phenotypic variation in rhizosphere Pseudomonas bacteria

Colony phase variation is a regulatory mechanism at the DNA level which usually results in high frequency, reversible switches between colonies with a different phenotype. A number of molecular mechanisms underlying phase variation are known: slipped-strand mispairing, genomic rearrangements, spontaneous mutations and epigenetic mechanisms such as differential methylation. Most examples of phenotypic variation or phase variation have been described in the context of host-pathogen interactions as mechanisms allowing pathogens to evade host immune responses. Recent reports indicate that phase variation is also relevant in competitive root colonization and biological control of phytopathogens. Many rhizospere Pseudomonas species show phenotypic variation, based on spontaneous mutation of the gacA and gacS genes. These morphological variants do not express secondary metabolites and have improved growth characteristics. The latter could contribute to efficient root colonization and success in competition, especially since (as shown for one strain) these variants were observed to revert to their wild-type form. The observation that these variants are present in rhizosphere-competent Pseudomonas bacteria suggests the existence of a conserved strategy to increase their success in the rhizosphere.

Comparative overview of the genomic and genetic differences between the pathogenic Neisseria strains and species

The availability of complete genome sequences from multiple pathogenic Neisseria strains and species has enabled a comprehensive survey of the genomic and genetic differences occurring within these species. In this review, we describe the chromosomal rearrangements that have occurred, and the genomic islands and prophages that have been identified in the various genomes. We also describe instances where specific genes are present or absent, other instances where specific genes have been inactivated, and situations where there is variation in the version of a gene that is present. We also provide an overview of mosaic genes present in these genomes, and describe the variation systems that allow the expression of particular genes to be switched ON or OFF. We have also described the presence and location of mobile non-coding elements in the various genomes. Finally, we have reviewed the incidence and properties of various extra-chromosomal elements found within these species. The overall impression is one of genomic variability and instability, resulting in increased functional flexibility within these species.

Role of RpoS and MutS in phase variation of Pseudomonas sp. PCL1171

Pseudomonas sp. strain PCL1171 undergoes reversible colony phase variation between opaque phase I and translucent phase II colonies, which is dependent on spontaneous mutations in the regulatory genes gacA and gacS. Mutation of the mutS gene and constitutive expression of rpoS increases the frequency at which gac mutants appear 1000- and 10-fold, respectively. Experiments were designed to study the relationship between gacS, rpoS and mutS. These studies showed that (i) a functional gac system is required for the expression of rpoS, (ii) RpoS suppresses the expression of mutS and therefore increases the frequency of gac mutants, and (iii) upon mutation of rpoS and gacS, the expression of mutS is increased. Mutation of gacS abolishes suppression of mutS expression in stationary growth, suggesting that additional gac-dependent factors are involved in this suppression. In conclusion, inefficient mutation repair via MutS, of which the expression is influenced by gacA/S itself and by rpoS in combination with other factors, contributes to the high frequency of mutations accumulating in gacA/S. The role of RpoS in the growth advantage of a gac mutant was analysed, and mutation of rpoS only reduced the length of the lag phase, but did not affect the growth rate, suggesting a role for both RpoS and a reduction of metabolic load in the growth advantage of a gac mutant.

Diversity in times of adversity: probabilistic strategies in microbial survival games

Population diversification strategies are ubiquitous among microbes, encompassing random phase-variation (RPV) of pathogenic bacteria, viral latency as observed in some bacteriophage and HIV, and the non-genetic diversity of bacterial stress responses. Precise conditions under which these diversification strategies confer an advantage have not been well defined. We develop a model of population growth conditioned on dynamical environmental and cellular states. Transitions among cellular states, in turn, may be biased by possibly noisy readings of the environment from cellular sensors. For various types of environmental dynamics and cellular sensor capability, we apply game-theoretic analysis to derive the evolutionarily stable strategy (ESS) for an organism and determine when that strategy is diversification. We find that: (1) RPV, effecting a sort of Parrondo paradox wherein random alternations between losing strategies produce a winning strategy, is selected when transitions between different selective environments cannot be sensed, (2) optimal RPV cell switching rates are a function of environmental lifecycle asymmetries and environmental autocorrelation, (3) probabilistic diversification upon entering a new environment is selected when sensors can detect environmental transitions but have poor precision in identifying new environments, and (4) in the presence of excess additive noise, low-pass filtering is required for evolutionary stability. We show that even when RPV is not the ESS, it may minimize growth rate variance and the risk of extinction due to 'unlucky' environmental dynamics.

Strain-specific differences in Neisseria gonorrhoeae associated with the phase variable gene repertoire

Background

There are several differences associated with the behaviour of the four main experimental Neisseria gonorrhoeae strains, FA1090, FA19, MS11, and F62. Although there is data concerning the gene complements of these strains, the reasons for the behavioural differences are currently unknown. Phase variation is a mechanism that occurs commonly within the Neisseria spp. and leads to switching of genes ON and OFF. This mechanism may provide a means for strains to express different combinations of genes, and differences in the strain-specific repertoire of phase variable genes may underlie the strain differences.

Results

By genome comparison of the four publicly available neisserial genomes a revised list of 64 genes was created that have the potential to be phase variable in N. gonorrhoeae, excluding the opa and pilC genes. Amplification and sequencing of the repeat-containing regions of these genes allowed determination of the presence of the potentially unstable repeats and the ON/OFF expression state of these genes. 35 of the 64 genes show differences in the composition or length of the repeats, of which 28 are likely to be associated with phase variation. Two genes were expressed differentially between strains causing disseminated infection and uncomplicated gonorrhoea. Further study of one of these in a range of clinical isolates showed this association to be due to sample size and is not maintained in a larger sample.

Conclusion

The results provide us with more evidence as to which genes identified through comparative genomics are indeed phase variable. The study indicates that there are large differences between these four N. gonorrhoeae strains in terms of gene expression during in vitro growth. It does not, however, identify any clear patterns by which previously reported behavioural differences can be correlated with the phase variable gene repertoire.

Involvement of genes of genome maintenance in the regulation of phase variation frequencies in Neisseria meningitidis

In Neisseria meningitidis, the reversible expression of surface antigens, i.e. phase variation, results from changes within repeated simple sequence motifs located in coding or promoter regions of the genes involved in their biosynthesis. The mutation rates of these simple sequences, which have a major influence on the generation of phenotypic diversity, can affect the fitness of the population. The aim of the present study was to investigate the involvement of genetic factors involved (mutS and dam) and not yet analysed (drg and dinB) in the regulation of phase variation frequencies of genes associated with a variety of repeat tracts. The frequency of frameshifts occurring in the polycytidine (polyC) tracts associated with siaD, spr and lgtG and in the tetranucleotide (TAAA) repeat tract associated with nadA was determined by colony immunoblotting or using the lacZ gene as a reporter. Inactivation of mutS increased the frequency of phase variation of genes presenting homopolymeric tracts of diverse length. Overexpression of dinB enhanced the instability of the homopolymeric tract associated with siaD. Investigation of the dam locus in a population of genetically distinct N. meningitidis strains revealed that 27 % of strains associated with invasive disease contained the dam gene. In all strains where a Dam function was absent, the drg gene had been inserted into the dam locus. Disruption of dam and drg in strains representative of each genotype, i.e. dam(+)/drg and dam/drg(+), did not modify phase variation frequencies. In contrast to the effects of certain genes on homopolymeric tracts, none of the genetic factors investigated affected the stability of tetranucleotide repeat tracts.

Phase variation mediated niche adaptation during prolonged experimental murine infection with Helicobacter pylori

Changes in the repeats associated with the recently redefined repertoire of 31 phase-variable genes in Helicobacter pylori were investigated following murine gastric colonization for up to one year in three unrelated H. pylori strains. Between the beginning and end of the experimental period, changes were seen in ten genes (32 %), which would alter gene expression in one or more of the three strains studied. For those genes that showed repeat length changes at the longest time points, intermediate time points showed differences between the rates of change for different functional groups of genes. Genes most likely to be associated with immediate niche fitting changed most rapidly, including phospholipase A (pldA) and LPS biosynthetic genes. Other surface proteins, which may be under adaptive immune selection, changed more slowly. Restriction-modification genes showed no particular temporal pattern. The number of genes that phase varied during adaptation to the murine gastric environment correlated inversely with their relative fitness as previously determined in this murine model of colonization. This suggests a role for these genes in determining initial fitness for colonization as well as in subsequent niche adaptation. In addition, a coding tandem repeat within a phase-variable gene which does not control actual gene expression was also investigated. This repeat was found to vary in copy number during colonization. This suggests that changes in the structures encoded by tandem repeats may also play a role in altered protein functions and/or immune evasion during H. pylori colonization.

Molecular nature of spontaneous modifications in gacS which cause colony phase variation in Pseudomonas sp. strain PCL1171

Pseudomonas sp. strain PCL1171 displays colony phase variation between opaque phase I and translucent phase II colonies, thereby regulating the production of secondary metabolites and exoenzymes. Complementation and sequence analysis of 26 phase II mutants and of 13 wild-type phase II sectors growing out of phase I colonies showed that in all these cases the phase II phenotype is caused by spontaneous mutations in gacA or/and gacS. Mutation of gac reduced both the length of the lag phase and the generation time. Isolation and sequencing of the gacS genes from the phase II bacteria revealed one insertion as well as several random point mutations, deletions, and DNA rearrangements. Most phase II colonies reverted with a high frequency, resulting in wild-type gacA and gacS genes and a phase I phenotype. Some phase II bacteria retained the phase II phenotype but changed genotypically as a result of (re)introduction of mutations in either gacA or gacS. The reversion of gacA or gacS to the wild type was not affected by mutation of recA and recB. We conclude that in Pseudomonas sp. strain PCL1171, mutations in gacA and gacS are the basis for phase variation from phase I to phase II colonies and that, since these mutations are efficiently removed, mutations in gac result in dynamic switches between the "wild-type" population and the subpopulations harboring spontaneous mutations in gacA and or gacS, thereby enabling both populations to be maintained.

Destabilization of tetranucleotide repeats in Haemophilus influenzae mutants lacking RnaseHI or the Klenow domain of PolI

A feature of Haemophilus influenzae genomes is the presence of several loci containing tracts of six or more identical tetranucleotide repeat units. These repeat tracts are unstable and mediate high frequency, reversible alterations in the expression of surface antigens. This process, termed phase variation (PV), enables H.influenzae to rapidly adapt to fluctuations in the host environment. Perturbation of lagging strand DNA synthesis is known to destabilize simple sequence repeats in yeast and Escherichia coli. By using a chromosomally located reporter construct, we demonstrated that the mutation of an H.influenzae rnhA (encoding RnaseHI) homologue increases the mutation rates of tetranucleotide repeats ∼3-fold. Additionally, deletion of the Klenow domain of DNA polymerase I (PolI) resulted in a ∼35-fold increase in tetranucleotide repeat-mediated PV rates. Deletion of the PolI 5′>3′ exonuclease domain appears to be lethal. The phenotypes of these mutants suggest that delayed or mutagenic Okazaki fragment processing destabilizes H.influenzae tetranucleotide repeat tracts.

Phase and antigenic variation in bacteria

Phase and antigenic variation result in a heterogenic phenotype of a clonal bacterial population, in which individual cells either express the phase-variable protein(s) or not, or express one of multiple antigenic forms of the protein, respectively. This form of regulation has been identified mainly, but by no means exclusively, for a wide variety of surface structures in animal pathogens and is implicated as a virulence strategy. This review provides an overview of the many bacterial proteins and structures that are under the control of phase or antigenic variation. The context is mainly within the role of the proteins and variation for pathogenesis, which reflects the main body of literature. The occurrence of phase variation in expression of genes not readily recognizable as virulence factors is highlighted as well, to illustrate that our current knowledge is incomplete. From recent genome sequence analysis, it has become clear that phase variation may be more widespread than is currently recognized, and a brief discussion is included to show how genome sequence analysis can provide novel information, as well as its limitations. The current state of knowledge of the molecular mechanisms leading to phase variation and antigenic variation are reviewed, and the way in which these mechanisms form part of the general regulatory network of the cell is addressed. Arguments both for and against a role of phase and antigenic variation in immune evasion are presented and put into new perspective by distinguishing between a role in bacterial persistence in a host and a role in facilitating evasion of cross-immunity. Finally, examples are presented to illustrate that phase-variable gene expression should be taken into account in the development of diagnostic assays and in the interpretation of experimental results and epidemiological studies.