Effects of antibiotic resistance alleles on bacterial evolutionary responses to viral parasites

Antibiotic resistance has wide-ranging effects on bacterial phenotypes and evolution. However, the influence of antibiotic resistance on bacterial responses to parasitic viruses remains unclear, despite the ubiquity of such viruses in nature and current interest in therapeutic applications. We experimentally investigated this by exposing various Escherichia coli genotypes, including eight antibiotic-resistant genotypes and a mutator, to different viruses (lytic bacteriophages). Across 960 populations, we measured changes in population density and sensitivity to viruses, and tested whether variation among bacterial genotypes was explained by their relative growth in the absence of parasites, or mutation rate towards phage resistance measured by fluctuation tests for each phage. We found that antibiotic resistance had relatively weak effects on adaptation to phages, although some antibiotic-resistance alleles impeded the evolution of resistance to phages via growth costs. By contrast, a mutator allele, often found in antibiotic-resistant lineages in pathogenic populations, had a relatively large positive effect on phage-resistance evolution and population density under parasitism. This suggests costs of antibiotic resistance may modify the outcome of phage therapy against pathogenic populations previously exposed to antibiotics, but the effects of any co-occurring mutator alleles are likely to be stronger.     

Carriage of λ Latent Virus Is Costly for Its Bacterial Host due to Frequent Reactivation in Monoxenic Mouse Intestine

Temperate phages, the bacterial viruses able to enter in a dormant prophage state in bacterial genomes, are present in the majority of bacterial strains for which the genome sequence is available. Although these prophages are generally considered to increase their hosts’ fitness by bringing beneficial genes, studies demonstrating such effects in ecologically relevant environments are relatively limited to few bacterial species. Here, we investigated the impact of prophage carriage in the gastrointestinal tract of monoxenic mice. Combined with mathematical modelling, these experimental results provided a quantitative estimation of key parameters governing phage-bacteria interactions within this model ecosystem. We used wild-type and mutant strains of the best known host/phage pair, Escherichia coli and phage λ. Unexpectedly, λ prophage caused a significant fitness cost for its carrier, due to an induction rate 50-fold higher than in vitro, with 1 to 2% of the prophage being induced. However, when prophage carriers were in competition with isogenic phage susceptible bacteria, the prophage indirectly benefited its carrier by killing competitors: infection of susceptible bacteria led to phage lytic development in about 80% of cases. The remaining infected bacteria were lysogenized, resulting overall in the rapid lysogenization of the susceptible lineage. Moreover, our setup enabled to demonstrate that rare events of phage gene capture by homologous recombination occurred in the intestine of monoxenic mice. To our knowledge, this study constitutes the first quantitative characterization of temperate phage-bacteria interactions in a simplified gut environment. The high prophage induction rate detected reveals DNA damage-mediated SOS response in monoxenic mouse intestine. We propose that the mammalian gut, the most densely populated bacterial ecosystem on earth, might foster bacterial evolution through high temperate phage activity.     

Malthusian Parameters as Estimators of the Fitness of Microbes: A Cautionary Tale about the Low Side of High Throughput

The maximum exponential growth rate, the Malthusian parameter (MP), is commonly used as a measure of fitness in experimental studies of adaptive evolution and of the effects of antibiotic resistance and other genes on the fitness of planktonic microbes. Thanks to automated, multi-well optical density plate readers and computers, with little hands-on effort investigators can readily obtain hundreds of estimates of MPs in less than a day. Here we compare estimates of the relative fitness of antibiotic susceptible and resistant strains of E. coli, Pseudomonas aeruginosa and Staphylococcus aureus based on MP data obtained with automated multi-well plate readers with the results from pairwise competition experiments. This leads us to question the reliability of estimates of MP obtained with these high throughput devices and the utility of these estimates of the maximum growth rates to detect fitness differences.     

Genome-driven elucidation of phage-host interplay and impact of phage resistance evolution on bacterial fitness

When considering the interactions between bacteriophages and their host, the issue of phage-resistance emergence is a key element in understanding the ecological impact of phages on the bacterial population. It is also an essential parameter for the implementation of phage therapy to combat antibiotic-resistant pathogens. This study investigates the phenotypic and genetic responses of five Pseudomonas aeruginosa strains (PAO1, A5803, AA43, CHA, and PAK) to the infection by seven phages with distinct evolutionary backgrounds and recognised receptors (LPS/T4P). Emerging phage-insensitivity was generally accompanied by self and cross-resistance mechanisms. Significant differences were observed between the reference PAO1 responses compared to other clinical representatives. LPS-dependent phage infections in clinical strains selected for mutations in the “global regulatory” and “other” genes, rather than in the LPS-synthesis clusters detected in PAO1 clones. Reduced fitness, as proxied by the growth rate, was correlated with large deletion (20–500 kbp) and phage carrier state. Multi-phage resistance was significantly correlated with a reduced growth rate but only in the PAO1 population. In addition, we observed that the presence of prophages decreased the lytic phage maintenance seemingly protecting the host against carrier state and occasional lytic phage propagation, thus preventing a significant reduction in bacterial growth rate.Subject terms: Bacteriophages, Biodiversity     

Effects of sequential and simultaneous applications of bacteriophages on populations of Pseudomonas aeruginosa in vitro and in wax moth larvae

Interest in using bacteriophages to treat bacterial infections (phage therapy) is growing, but there have been few experiments comparing the effects of different treatment strategies on both bacterial densities and resistance evolution. While it is established that multiphage therapy is typically more effective than the application of a single phage type, it is not clear if it is best to apply phages simultaneously or sequentially. We tried single- and multiphage therapy against Pseudomonas aeruginosa PAO1 in vitro, using different combinations of phages either simultaneously or sequentially. Across different phage combinations, simultaneous application was consistently equal or superior to sequential application in terms of reducing bacterial population density, and there was no difference (on average) in terms of minimizing resistance. Phage-resistant bacteria emerged in all experimental treatments and incurred significant fitness costs, expressed as reduced growth rate in the absence of phages. Finally, phage therapy increased the life span of wax moth larvae infected with P. aeruginosa, and a phage cocktail was the most effective short-term treatment. When the ratio of phages to bacteria was very high, phage cocktails cured otherwise lethal infections. These results suggest that while adding all available phages simultaneously tends to be the most successful short-term strategy, there are sequential strategies that are equally effective and potentially better over longer time scales.     

A Window of Opportunity to Control the Bacterial Pathogen Pseudomonas aeruginosa Combining Antibiotics and Phages

The evolution of antibiotic resistance in bacteria is a global concern and the use of bacteriophages alone or in combined therapies is attracting increasing attention as an alternative. Evolutionary theory predicts that the probability of bacterial resistance to both phages and antibiotics will be lower than to either separately, due for example to fitness costs or to trade-offs between phage resistance mechanisms and bacterial growth. In this study, we assess the population impacts of either individual or combined treatments of a bacteriophage and streptomycin on the nosocomial pathogen Pseudomonas aeruginosa. We show that combining phage and antibiotics substantially increases bacterial control compared to either separately, and that there is a specific time delay in antibiotic introduction independent of antibiotic dose, that minimizes both bacterial density and resistance to either antibiotics or phage. These results have implications for optimal combined therapeutic approaches.     

Construction of two ureolytic model organisms for the study of microbially induced calcium carbonate precipitation

Two bacterial strains, Pseudomonas aeruginosa MJK1 and Escherichia coli MJK2, were constructed that both express green fluorescent protein (GFP) and carry out ureolysis. These two novel model organisms are useful for studying bacterial carbonate mineral precipitation processes and specifically ureolysis-driven microbially induced calcium carbonate precipitation (MICP). The strains were constructed by adding plasmid-borne urease genes (ureABC, ureD and ureFG) to the strains P. aeruginosa AH298 and E. coli AF504gfp, both of which already carried unstable GFP derivatives. The ureolytic activities of the two new strains were compared to the common, non-GFP expressing, model organism Sporosarcina pasteurii in planktonic culture under standard laboratory growth conditions. It was found that the engineered strains exhibited a lower ureolysis rate per cell but were able to grow faster and to a higher population density under the conditions of this study. Both engineered strains were successfully grown as biofilms in capillary flow cell reactors and ureolysis-induced calcium carbonate mineral precipitation was observed microscopically. The undisturbed spatiotemporal distribution of biomass and calcium carbonate minerals were successfully resolved in 3D using confocal laser scanning microscopy. Observations of this nature were not possible previously because no obligate urease producer that expresses GFP had been available. Future observations using these organisms will allow researchers to further improve engineered application of MICP as well as study natural mineralization processes in model systems.     

Relationship between fitness and heterogeneity in exponentially growing microbial populations

Despite major environmental and genetic differences, microbial metabolic networks are known to generate consistent physiological outcomes across vastly different organisms. This remarkable robustness suggests that, at least in bacteria, metabolic activity may be guided by universal principles. The constrained optimization of evolutionarily motivated objective functions, such as the growth rate, has emerged as the key theoretical assumption for the study of bacterial metabolism. While conceptually and practically useful in many situations, the idea that certain functions are optimized is hard to validate in data. Moreover, it is not always clear how optimality can be reconciled with the high degree of single-cell variability observed in experiments within microbial populations. To shed light on these issues, we develop an inverse modeling framework that connects the fitness of a population of cells (represented by the mean single-cell growth rate) to the underlying metabolic variability through the maximum entropy inference of the distribution of metabolic phenotypes from data. While no clear objective function emerges, we find that, as the medium gets richer, the fitness and inferred variability for Escherichia coli populations follow and slowly approach the theoretically optimal bound defined by minimal reduction of variability at given fitness. These results suggest that bacterial metabolism may be crucially shaped by a population-level trade-off between growth and heterogeneity.     

Differential freshwater flagellate community response to bacterial food quality with a focus on Limnohabitans bacteria

Different bacterial strains can have different value as food for heterotrophic nanoflagellates (HNF), thus modulating HNF growth and community composition. We examined the influence of prey food quality using four Limnohabitans strains, one Polynucleobacter strain and one freshwater actinobacterial strain on growth (growth rate, length of lag phase and growth efficiency) and community composition of a natural HNF community from a freshwater reservoir. Pyrosequencing of eukaryotic small subunit rRNA amplicons was used to assess time-course changes in HNF community composition. All four Limnohabitans strains and the Polynucleobacter strain yielded significant HNF community growth while the actinobacterial strain did not although it was detected in HNF food vacuoles. Notably, even within the Limnohabitans strains we found significant prey-related differences in HNF growth parameters, which could not be related only to size of the bacterial prey. Sequence data characterizing the HNF communities showed also that different bacterial prey items induced highly significant differences in community composition of flagellates. Generally, Stramenopiles dominated the communities and phylotypes closely related to Pedospumella (Chrysophyceae) were most abundant bacterivorous flagellates rapidly reacting to addition of the bacterial prey of high food quality.     

Phage formation in Staphylococcus muscae cultures. VIII. Effect of the protein factor and aspartic acid on virus synthesis with various bacterial strains

1. Four strains of Staphylococcus muscae have been isolated which differ in their growth rates and phage syntheses in Fildes' synthetic medium. 2. Two of the strains when singly infected cannot release phage in Fildes' synthetic medium unless a substance present in certain acid-hydrolyzed proteins is added to the medium. One of these strains also requires other substance(s) present in acid-hydrolyzed proteins in order to grow in Fildes' medium. 3. The two strains which do not require the addition of the phage-stimulating factor have been found either to synthesize this substance, or one similar to it. One of these strains will not grow in Fildes' medium unless substance(s) present in acid-hydrolyzed proteins is added to the medium. 4. The purified acid-hydrolyzed protein factor necessary for virus liberation does not affect the multiplication rate of uninfected S. muscae cells in Fildes' synthetic medium. 5. The substance is not needed for the adsorption or the invasion of the host cell by the virus. In the absence of the factor, the virus is adsorbed to the cell and "kills" it. 6. An analysis carried out by means of the one-step growth curve technique has indicated that the substance is not concerned simply with the mechanism of virus release, but is necessary for some initial stage in virus synthesis. 7. With one bacterial strain not requiring the AHPF, aspartic acid had to be present at least during the minimum latent period for the cell to form virus. 8. In the absence of aspartic acid, the virus was adsorbed to the cell and killed it, but no virus was released from singly infected bacteria. 9. If the cells were grown in a medium containing aspartic acid and then resuspended in the medium minus aspartic acid, no virus was released, although such cells contained at least two times the amount of aspartic acid necessary for the burst size in the complete medium. 10. Aspartic acid, a constituent of the virus particle, appears from an analysis of one-step growth curves to take part in the initial phase of phage synthesis. 11. The effect of amino acids on virus formation is discussed in relation to the time sequence of virus protein and desoxyribonucleic acid synthesis.     

Phage and bacteria support mutual diversity in a narrowing staircase of coexistence

The competitive exclusion principle states that phage diversity M should not exceed bacterial diversity N. By analyzing the steady-state solutions of multistrain equations, we find a new constraint: the diversity N of bacteria living on the same resources is constrained to be M or M+1 in terms of the diversity of their phage predators. We quantify how the parameter space of coexistence exponentially decreases with diversity. For diversity to grow, an open or evolving ecosystem needs to climb a narrowing ‘diversity staircase'' by alternatingly adding new bacteria and phages. The unfolding coevolutionary arms race will typically favor high growth rate, but a phage that infects two bacterial strains differently can occasionally eliminate the fastest growing bacteria. This context-dependent fitness allows abrupt resetting of the ‘Red-Queen''s race'' and constrains the local diversity.     

<Emphasis Type="Italic">In silico</Emphasis> comparison of bacterial strains using mutual information

Fast-sequencing throughput methods have increased the number of completely sequenced bacterial genomes to about 400 by December 2006, with the number increasing rapidly. These include several strains. In silico methods of comparative genomics are of use in categorizing and phylogenetically sorting these bacteria. Various word-based tools have been used for quantifying the similarities and differences between entire genomes. The simple di-nucleotide frequency comparison, codon specificity and k-mer repeat detection are among some of the well-known methods.In this paper, we show that the Mutual Information function, which is a measure of correlations and a concept from Information Theory, is very effective in determining the similarities and differences among genome sequences of various strains of bacteria such as the plant pathogen Xylella fastidiosa, marine Cyanobacteria Prochlorococcus marinus or animal and human pathogens such as species of Ehrlichia and Legionella. The short-range three-base periodicity, small sequence repeats and long-range correlations taken together constitute a genome signature that can be used as a technique for identifying new bacterial strains with the help of strains already catalogued in the database.There have been several applications of using the Mutual Information function as a measure of correlations in genomics but this is the first whole genome analysis done to detect strain similarities and differences.     

A Pipeline for Screening Small Molecules with Growth Inhibitory Activity against Burkholderia cenocepacia

Infections with the bacteria Burkholderia cepacia complex (Bcc) are very difficult to eradicate in cystic fibrosis patients due the intrinsic resistance of Bcc to most available antibiotics and the emergence of multiple antibiotic resistant strains during antibiotic treatment. In this work, we used a whole-cell based assay to screen a diverse collection of small molecules for growth inhibitors of a relevant strain of Bcc, B. cenocepacia K56-2. The primary screen used bacterial growth in 96-well plate format and identified 206 primary actives among 30,259 compounds. From 100 compounds with no previous record of antibacterial activity secondary screening and data mining selected a total of Bce bioactives that were further analyzed. An experimental pipeline, evaluating in vitro antibacterial and antibiofilm activity, toxicity and in vivo antibacterial activity using C. elegans was used for prioritizing compounds with better chances to be further investigated as potential Bcc antibacterial drugs. This high throughput screen, along with the in vitro and in vivo analysis highlights the utility of this experimental method to quickly identify bioactives as a starting point of antibacterial drug discovery.     

Isolation and Partial Characterization of a Virulent Bacteriophage IHQ1 Specific for Aeromonas punctata from Stream Water

Aeromonas punctata is the causative agent of septicemia, diarrhea, wound infections, meningitis, peritonitis, and infections of the joints, bones and eyes. Bacteriophages are often considered alternative agents for controlling bacterial infection and contamination. In this study, we described the isolation and preliminary characterization of bacteriophage IHQ1 (family Myoviridae) active against the Gram-negative bacterial strain A. punctata. This virulent bacteriophage was isolated from stream water sample. Genome analysis indicated that phage IHQ1 was a double-stranded DNA virus with an approximate genome size of 25–28 kb. The initial characterization of this newly isolated phage showed that it has a narrow host range and infects only A. punctata as it failed to infect seven other clinically isolated pathogenic strains, i.e., methicillin-resistant Staphylococcus aureus 6403, MRSA 17644, Acinetobacter 33408, Acinetobacter 1172, Pseudomonas aeruginosa 22250, P. aeruginosa 11219, and Escherichia coli. Proteomic pattern of phage IHQ1, generated by SDS-PAGE using purified phage particles, showed three major and three minor protein bands with molecular weights ranging from 25 to 70 kDa. The adsorption rate of phage IHQ1 to the host bacterium was also determined, which was significantly enhanced by the addition of 10 mM CaCl₂. From the single-step growth experiment, it was inferred that the latent time period of phage IHQ1 was 24 min and a burst size of 626 phages per cell. Moreover, the pH and thermal stability of phage IHQ1 were also investigated. The maximum stability of the phage was observed at optimal pH 7.0, and it was totally unstable at extreme acidic pH 3; however, it was comparatively stable at alkaline pH 11.0. At 37°C the phage showed maximum number of plaques, and the viability was almost 100%. The existence of Aeromonas bacteriophage is very promising for the eradication of this opportunistic pathogen and also for future applications such as the design of new detection and phage typing (diagnosis) methods. The specificity of the bacteriophage for A. punctata makes it an attractive candidate for phage therapy of A. punctata infections.     

A Bacteriophage-lysing Strain of Staphylococcus Employed in the Manufacture of Dry Sausage

A bacteriophage of a certain Staphylococcus (a strain of Staphylococcus lactis) employed in the manufacture of dry sausage has been characterized. The host range of this bacteriophage is wide. In addition to the original host, 15 other strains (out of 40 strains tested) were found to support reproduction of the phage. The sensitive strains represented Staphylococcus saprophyticus and different types of S. lactis.

The growth rate of the bacterial host did not influence the rates of phage adsorption, nor the maximal reproduction rate of new particles. With increasing bacterial growth rate, the “lag” observed before phage reproduction started was distinctly decreased. This phase was shorter with the original host strain than with other sensitive strains.

Resistant cultures of the original host strain were easily obtained. These cultures grew as rapidly and gave as good yields of cell mass as the original phage-sensitive host. However, phage resistance was frequently lost.

     

Plausibility of inferred ancestral phenotypes and the evaluation of alternative models of limb evolution in scincid lizards

Phylogenetic approaches to inferring ancestral character states are becoming increasingly sophisticated; however, the potential remains for available methods to yield strongly supported but inaccurate ancestral state estimates. The consistency of ancestral states inferred for two or more characters affords a useful criterion for evaluating ancestral trait reconstructions. Ancestral state estimates for multiple characters that entail plausible phenotypes when considered together may reasonably be assumed to be reliable. However, the accuracy of inferred ancestral states for one or more characters may be questionable where combined reconstructions imply implausible phenotypes for a proportion of internal nodes. This criterion for assessing reconstructed ancestral states is applied here in evaluating inferences of ancestral limb morphology in the scincid lizard clade Lerista. Ancestral numbers of digits for the manus and pes inferred assuming the models that best fit the data entail ancestral digit configurations for many nodes that differ fundamentally from configurations observed among known species. However, when an alternative model is assumed for the pes, inferred ancestral digit configurations are invariably represented among observed phenotypes. This indicates that a suboptimal model for the pes (and not the model providing the best fit to the data) yields accurate ancestral state estimates.     

Phage selection for bacterial cheats leads to population decline

While predators and parasites are known for their effects on bacterial population biology, their impact on the dynamics of bacterial social evolution remains largely unclear. Siderophores are iron-chelating molecules that are key to the survival of certain bacterial species in iron-limited environments, but their production can be subject to cheating by non-producing genotypes. In a selection experiment conducted over approximately 20 bacterial generations and involving 140 populations of the pathogenic bacterium Pseudomonas aeruginosa PAO1, we assessed the impact of a lytic phage on competition between siderophore producers and non-producers. We show that the presence of lytic phages favours the non-producing genotype in competition, regardless of whether iron use relies on siderophores. Interestingly, phage pressure resulted in higher siderophore production, which constitutes a cost to the producers and may explain why they were outcompeted by non-producers. By the end of the experiment, however, cheating load reduced the fitness of mixed populations relative to producer monocultures, and only monocultures of producers managed to grow in the presence of phage in situations where siderophores were necessary to access iron. These results suggest that public goods production may be modulated in the presence of natural enemies with consequences for the evolution of social strategies.     

Genomic and Gene-Expression Comparisons among Phage-Resistant Type-IV Pilus Mutants of Pseudomonas syringae pathovar phaseolicola

Pseudomonas syringae pv. phaseolicola (Pph) is a significant bacterial pathogen of agricultural crops, and phage Φ6 and other members of the dsRNA virus family Cystoviridae undergo lytic (virulent) infection of Pph, using the type IV pilus as the initial site of cellular attachment. Despite the popularity of Pph/phage Φ6 as a model system in evolutionary biology, Pph resistance to phage Φ6 remains poorly characterized. To investigate differences between phage Φ6 resistant Pph strains, we examined genomic and gene expression variation among three bacterial genotypes that differ in the number of type IV pili expressed per cell: ordinary (wild-type), non-piliated, and super-piliated. Genome sequencing of non-piliated and super-piliated Pph identified few mutations that separate these genotypes from wild type Pph–and none present in genes known to be directly involved in type IV pilus expression. Expression analysis revealed that 81.1% of gene ontology (GO) terms up-regulated in the non-piliated strain were down-regulated in the super-piliated strain. This differential expression is particularly prevalent in genes associated with respiration—specifically genes in the tricarboxylic acid cycle (TCA) cycle, aerobic respiration, and acetyl-CoA metabolism. The expression patterns of the TCA pathway appear to be generally up and down-regulated, in non-piliated and super-piliated Pph respectively. As pilus retraction is mediated by an ATP motor, loss of retraction ability might lead to a lower energy draw on the bacterial cell, leading to a different energy balance than wild type. The lower metabolic rate of the super-piliated strain is potentially a result of its loss of ability to retract.     

Phage loss and the breakdown of a defensive symbiosis in aphids

Terrestrial arthropods are often infected with heritable bacterial symbionts, which may themselves be infected by bacteriophages. However, what role, if any, bacteriophages play in the regulation and maintenance of insect–bacteria symbioses is largely unknown. Infection of the aphid Acyrthosiphon pisum by the bacterial symbiont Hamiltonella defensa confers protection against parasitoid wasps, but only when H. defensa is itself infected by the phage A. pisum secondary endosymbiont (APSE). Here, we use a controlled genetic background and correlation-based assays to show that loss of APSE is associated with up to sevenfold increases in the intra-aphid abundance of H. defensa. APSE loss is also associated with severe deleterious effects on aphid fitness: aphids infected with H. defensa lacking APSE have a significantly delayed onset of reproduction, lower weight at adulthood and half as many total offspring as aphids infected with phage-harbouring H. defensa, indicating that phage loss can rapidly lead to the breakdown of the defensive symbiosis. Our results overall indicate that bacteriophages play critical roles in both aphid defence and the maintenance of heritable symbiosis.     

Effects of bacterial secondary symbionts on host plant use in pea aphids

Aphids possess several facultative bacterial symbionts that have important effects on their hosts'' biology. These have been most closely studied in the pea aphid (Acyrthosiphon pisum), a species that feeds on multiple host plants. Whether secondary symbionts influence host plant utilization is unclear. We report the fitness consequences of introducing different strains of the symbiont Hamiltonella defensa into three aphid clones collected on Lathyrus pratensis that naturally lack symbionts, and of removing symbionts from 20 natural aphid–bacterial associations. Infection decreased fitness on Lathyrus but not on Vicia faba, a plant on which most pea aphids readily feed. This may explain the unusually low prevalence of symbionts in aphids collected on Lathyrus. There was no effect of presence of symbiont on performance of the aphids on the host plants of the clones from which the H. defensa strains were isolated. Removing the symbiont from natural aphid–bacterial associations led to an average approximate 20 per cent reduction in fecundity, both on the natural host plant and on V. faba, suggesting general rather than plant-species-specific effects of the symbiont. Throughout, we find significant genetic variation among aphid clones. The results provide no evidence that secondary symbionts have a major direct role in facilitating aphid utilization of particular host plant species.