Quantitative analysis of the bacteriophage Qβ infection cycle

Droplet-based methodology for investigating bacterial population dynamics in response to phage exposure

An alarming rise in antimicrobial resistance worldwide has spurred efforts into the search for alternatives to antibiotic treatments. The use of bacteriophages, bacterial viruses harmless to humans, represents a promising approach with potential to treat bacterial infections (phage therapy). Recent advances in microscopy-based single-cell techniques have allowed researchers to develop new quantitative methodologies for assessing the interactions between bacteria and phages, especially the ability of phages to eradicate bacterial pathogen populations and to modulate growth of both commensal and pathogen populations. Here we combine droplet microfluidics with fluorescence time-lapse microscopy to characterize the growth and lysis dynamics of the bacterium Escherichia coli confined in droplets when challenged with phage. We investigated phages that promote lysis of infected E. coli cells, specifically, a phage species with DNA genome, T7 (Escherichia virus T7) and two phage species with RNA genomes, MS2 (Emesvirus zinderi) and Qβ (Qubevirus durum). Our microfluidic trapping device generated and immobilized picoliter-sized droplets, enabling stable imaging of bacterial growth and lysis in a temperature-controlled setup. Temporal information on bacterial population size was recorded for up to 25 h, allowing us to determine growth rates of bacterial populations and helping us uncover the extent and speed of phage infection. In the long-term, the development of novel microfluidic single-cell and population-level approaches will expedite research towards fundamental understanding of the genetic and molecular basis of rapid phage-induced lysis and eco-evolutionary aspects of bacteria-phage dynamics, and ultimately help identify key factors influencing the success of phage therapy.

Isolation and Molecular Characterization of Two Novel Lytic Bacteriophages for the Biocontrol of Escherichia coli in Uterine Infections: In Vitro and Ex Vivo Preliminary Studies in Veterinary Medicine

E. coli is one of the etiological agents responsible for pyometra in female dogs, with conventional treatment involving ovariohysterectomy. Here, we report the isolation and full characterization of two novel lytic phages, viz. vB_EcoM_Uniso11 (ph0011) and vB_EcoM_Uniso21 (ph0021). Both phages belong to the order Caudovirales and present myovirus-like morphotypes, with phage ph0011 being classified as Myoviridae genus Asteriusvirus and phage ph0021 being classified as Myoviridae genus Tequatrovirus, based on their complete genome sequences. The 348,288 bp phage ph0011 and 165,222 bp phage ph0021 genomes do not encode toxins, integrases or antimicrobial resistance genes neither depolymerases related sequences. Both phages were shown to be effective against at least twelve E. coli clinical isolates in in vitro antibacterial activity assays. Based on their features, both phages have potential for controlling pyometra infections caused by E. coli. Phage ph0011 (reduction of 4.24 log CFU/mL) was more effective than phage ph0021 (reduction of 1.90 log CFU/mL) after 12 h of incubation at MOI 1000. As a cocktail, the two phages were highly effective in reducing the bacterial load (reduction of 5.57 log CFU/mL) at MOI 100, after 12 h of treatment. Both phages were structurally and functionally stabilized in vaginal egg formulations.

Patterns of virus growth across the diversity of life

Although viruses in their natural habitats add up to less than 10% of the biomass, they contribute more than 90% of the genome sequences [1]. These viral sequences or 'viromes' encode viruses that populate the Earth's oceans [2, 3] and terrestrial environments [4, 5], where their infections impact life across diverse ecological niches and scales [6, 7], including humans [8-10]. Most viruses have yet to be isolated and cultured [11-13], and surprisingly few efforts have explored what analysis of available data might reveal about their nature. Here, we compiled and analyzed seven decades of one-step growth and other data for viruses from six major families, including their infections of archaeal, bacterial and eukaryotic hosts [14-191]. We found that the use of host cell biomass for virus production was highest for archaea at 10%, followed by bacteria at 1% and eukarya at 0.01%, highlighting the degree to which viruses of archaea and bacteria exploit their host cells. For individual host cells, the yield of virus progeny spanned a relatively narrow range (10-1000 infectious particles per cell) compared with the million-fold difference in size between the smallest and largest cells. Furthermore, healthy and infected host cells were remarkably similar in the time they needed to multiply themselves or their virus progeny. Specifically, the doubling time of healthy cells and the delay time for virus release from infected cells were not only correlated (r = 0.71, p < 10-10, n = 101); they also spanned the same range from tens of minutes to about a week. These results have implications for better understanding the growth, spread and persistence of viruses in complex natural habitats that abound with diverse hosts, including humans and their associated microbes.

Simulation of the M13 life cycle I: Assembly of a genetically-structured deterministic chemical kinetic simulation

To expand the quantitative, systems level understanding and foster the expansion of the biotechnological applications of the filamentous bacteriophage M13, we have unified the accumulated quantitative information on M13 biology into a genetically-structured, experimentally-based computational simulation of the entire phage life cycle. The deterministic chemical kinetic simulation explicitly includes the molecular details of DNA replication, mRNA transcription, protein translation and particle assembly, as well as the competing protein-protein and protein-nucleic acid interactions that control the timing and extent of phage production. The simulation reproduces the holistic behavior of M13, closely matching experimentally reported values of the intracellular levels of phage species and the timing of events in the M13 life cycle. The computational model provides a quantitative description of phage biology, highlights gaps in the present understanding of M13, and offers a framework for exploring alternative mechanisms of regulation in the context of the complete M13 life cycle.

Virus wars: using one virus to block the spread of another

The failure of traditional interventions to block and cure HIV infections has led to novel proposals that involve treating infections with therapeutic viruses-infectious viruses that specifically inhibit HIV propagation in the host. Early efforts in evaluating these proposals have been limited chiefly to mathematical models of dynamics, for lack of suitable empirical systems. Here we propose, develop and analyze an empirical system of a therapeutic virus that protects a host cell population against a lethal virus. The empirical system uses E. coli bacteria as the host cell population, an RNA phage as the lethal virus and a filamentous phage as the therapeutic virus. Basic dynamic properties are established for each virus alone and then together. Observed dynamics broadly agree with those predicted by a computer simulation model, although some differences are noted. Two cases of dynamics are contrasted, differing in whether the therapeutic virus is introduced before the lethal virus or after the lethal virus. The therapeutic virus increases in both cases but by different mechanisms. With the therapeutic virus introduced first, it spreads infectiously without any appreciable change in host dynamics. With the therapeutic virus introduced second, host abundance is depressed at the time therapy is applied; following an initial period of therapeutic virus spread by infection, the subsequent rise of protection is through reproduction by hosts already protected. This latter outcome is due to inheritance of the therapeutic virus state when the protected cell divides. Overall, the work establishes the feasibility and robustness to details of a viral interference using a therapeutic virus.

Contribution of silent mutations to thermal adaptation of RNA bacteriophage Qβ

Changes in protein function and other biological properties, such as RNA structure, are crucial for adaptation of organisms to novel or inhibitory environments. To investigate how mutations that do not alter amino acid sequence may be positively selected, we performed a thermal adaptation experiment using the single-stranded RNA bacteriophage Qβ in which the culture temperature was increased from 37.2°C to 41.2°C and finally to an inhibitory temperature of 43.6°C in a stepwise manner in three independent lines. Whole-genome analysis revealed 31 mutations, including 14 mutations that did not result in amino acid sequence alterations, in this thermal adaptation. Eight of the 31 mutations were observed in all three lines. Reconstruction and fitness analyses of Qβ strains containing only mutations observed in all three lines indicated that five mutations that did not result in amino acid sequence changes but increased the amplification ratio appeared in the course of adaptation to growth at 41.2°C. Moreover, these mutations provided a suitable genetic background for subsequent mutations, altering the fitness contribution from deleterious to beneficial. These results clearly showed that mutations that do not alter the amino acid sequence play important roles in adaptation of this single-stranded RNA virus to elevated temperature.

IMPORTANCE

Recent studies using whole-genome analysis technology suggested the importance of mutations that do not alter the amino acid sequence for adaptation of organisms to novel environmental conditions. It is necessary to investigate how these mutations may be positively selected and to determine to what degree such mutations that do not alter amino acid sequences contribute to adaptive evolution. Here, we report the roles of these silent mutations in thermal adaptation of RNA bacteriophage Qβ based on experimental evolution during which Qβ showed adaptation to growth at an inhibitory temperature. Intriguingly, four synonymous mutations and one mutation in the untranslated region that spread widely in the Qβ population during the adaptation process at moderately high temperature provided a suitable genetic background to alter the fitness contribution of subsequent mutations from deleterious to beneficial at a higher temperature.

A(2) expression and assembly regulates lysis in Qβ infections

The capsids of ssRNA phages comprise a single copy of an ~45 kDa maturation protein that serves to recognize the conjugative pilus as receptor, to protect the ends of the viral RNA and also to escort the genomic RNA into the host cytoplasm. In the Alloleviviridae, represented by the canonical phage Qβ, the maturation protein A(2) also causes lysis. This is achieved by inhibiting the activity of MurA, which catalyses the first committed step of murein biosynthesis. Previously, it was shown that Qβ virions, with a single copy of A(2), inhibit MurA activity. This led to a model for lysis timing in which, during phage infection, A(2) is not active as a MurA inhibitor until assembled into virion particles, thus preventing premature lysis before a sufficient yield of viable progeny has accumulated. Here we report that MurA inactivates purified Qβ particles, casting doubt on the notion that A(2) must assemble into particles prior to MurA inhibition. Furthermore, quantification of A(2) protein induced from a plasmid indicated that lysis is entrained when the amount of the lysis protein is approximately equimolar to that of cellular MurA. Qβ por mutants, isolated as suppressors that overcome a murA(rat) mutation that reduces the affinity of MurA for A(2), were shown to be missense mutations in A(2) that increase the translation of the maturation protein. Because of the increased production of A(2), the por mutants have an attenuated infection cycle and reduced burst size, indicating that a delicate balance between assembled and unassembled A(2) levels regulates lysis timing.

The three faces of riboviral spontaneous mutation: spectrum, mode of genome replication, and mutation rate

Riboviruses (RNA viruses without DNA replication intermediates) are the most abundant pathogens infecting animals and plants. Only a few riboviral infections can be controlled with antiviral drugs, mainly because of the rapid appearance of resistance mutations. Little reliable information is available concerning i) kinds and relative frequencies of mutations (the mutational spectrum), ii) mode of genome replication and mutation accumulation, and iii) rates of spontaneous mutation. To illuminate these issues, we developed a model in vivo system based on phage Qß infecting its natural host, Escherichia coli. The Qß RT gene encoding the Read-Through protein was used as a mutation reporter. To reduce uncertainties in mutation frequencies due to selection, the experimental Qß populations were established after a single cycle of infection and selection against RT⁻ mutants during phage growth was ameliorated by plasmid-based RT complementation in trans. The dynamics of Qß genome replication were confirmed to reflect the linear process of iterative copying (the stamping-machine mode). A total of 32 RT mutants were detected among 7,517 Qß isolates. Sequencing analysis of 45 RT mutations revealed a spectrum dominated by 39 transitions, plus 4 transversions and 2 indels. A clear template•primer mismatch bias was observed: A•C>C•A>U•G>G•U> transversion mismatches. The average mutation rate per base replication was ≈9.1×10⁻⁶ for base substitutions and ≈2.3×10⁻⁷ for indels. The estimated mutation rate per genome replication, μ_g, was ≈0.04 (or, per phage generation, ≈0.08), although secondary RT mutations arose during the growth of some RT mutants at a rate about 7-fold higher, signaling the possible impact of transitory bouts of hypermutation. These results are contrasted with those previously reported for other riboviruses to depict the current state of the art in riboviral mutagenesis.

Viral disease is a subject of major concern in public health. Diseases produced by riboviruses (RNA viruses sensu stricto) represent a special urgency, because these viruses display an exceptional capability to generate resistance mutations against antiviral drugs. Unfortunately, little is known about the rate and nature of spontaneous mutation in riboviruses. Thus, characterization of their mutation process may be helpful in the development of improved ways to counteract riboviral diseases. In this study, we investigated the mutation process in vivo of a model ribovirus, the bacteriophage Qß, focusing on three key aspects: i) the kinds and relative frequencies of mutations, ii) the mode of genome replication, and iii) the rate of spontaneous mutation. Our results, combined with other information about riboviral mutagenesis, depict a ribovirus mutation spectrum largely dominated by transitions, a predominantly linear mode of genome replication, and a mutation rate per genome replication on the order of 0.04 for bacteriophages and plant viruses but perhaps an order of magnitude higher for mammalian riboviruses.

The incorporation of the A2 protein to produce novel Qβ virus-like particles using cell-free protein synthesis

Virus-like particles (VLPs) have been employed for a number of nanometric applications because they self-assemble, exhibit a high degree of symmetry, and can be genetically and chemically modified. However, high symmetry does not allow for a single unique modification site on the VLP. Here, we demonstrate the co-expression of the cytotoxic A2 protein and the coat protein of the bacteriophage Qβ to form a nearly monodispersed population of novel VLPs. Cell-free protein synthesis allows for direct access and optimization of protein-synthesis and VLP-assembly. The A2 is shown to be incorporated at high efficiency, approaching a theoretical maximum of one A2 per VLP. This work demonstrates de novo production of a novel VLP, which contains a unique site that has the potential for future nanometric engineering applications.

Ongoing phenotypic and genomic changes in experimental coevolution of RNA bacteriophage Qβ and Escherichia coli

According to the Red Queen hypothesis or arms race dynamics, coevolution drives continuous adaptation and counter-adaptation. Experimental models under simplified environments consisting of bacteria and bacteriophages have been used to analyze the ongoing process of coevolution, but the analysis of both parasites and their hosts in ongoing adaptation and counter-adaptation remained to be performed at the levels of population dynamics and molecular evolution to understand how the phenotypes and genotypes of coevolving parasite–host pairs change through the arms race. Copropagation experiments with Escherichia coli and the lytic RNA bacteriophage Qβ in a spatially unstructured environment revealed coexistence for 54 days (equivalent to 163–165 replication generations of Qβ) and fitness analysis indicated that they were in an arms race. E. coli first adapted by developing partial resistance to infection and later increasing specific growth rate. The phage counter-adapted by improving release efficiency with a change in host specificity and decrease in virulence. Whole-genome analysis indicated that the phage accumulated 7.5 mutations, mainly in the A2 gene, 3.4-fold faster than in Qβ propagated alone. E. coli showed fixation of two mutations (in traQ and csdA) faster than in sole E. coli experimental evolution. These observations suggest that the virus and its host can coexist in an evolutionary arms race, despite a difference in genome mutability (i.e., mutations per genome per replication) of approximately one to three orders of magnitude.

To examine the ongoing changes driven by host–parasite interactions, we have constructed a coevolution model consisting of Escherichia coli and the lytic RNA bacteriophage Qβ (Qβ) in a spatially unstructured environment. In coevolution through 54 daily copropagations of the parasite and its host, E. coli first evolved partial resistance to infection and later accelerated its specific growth rate, while the phage counter-adapted by improving release efficiency with a change in host specificity and a decrease in virulence. Whole-genome analysis of E. coli and Qβ revealed accelerated molecular evolution in comparison with Qβ propagation in this study and E. coli sole passage reported previously. The results of the present study indicated that, despite the large difference in mutability of their genomes (approximately one to three orders of magnitude difference), a host with larger genome size (4.6 Mbp) and a lower spontaneous mutation rate (5.4×10⁻¹⁰ per bp per replication) and a parasite with a smaller genome size (4,217 bases) and a higher mutation rate (1.5×10⁻³ to 1.5×10⁻⁵ per base per replication) were capable of changing their phenotypes to coexist in an arms race.

Identification of two forms of Q{beta} replicase with different thermal stabilities but identical RNA replication activity

The enzyme Qβ replicase is an RNA-dependent RNA polymerase, which plays a central role in infection by the simple single-stranded RNA virus bacteriophage Qβ. This enzyme has been used in a number of applications because of its unique activity in amplifying RNA from an RNA template. Determination of the thermal stability of Qβ replicase is important to gain an understanding of its function and potential applications, but data reported to date have been contradictory. Here, we provide evidence that these previous inconsistencies were due to the heterogeneous forms of the replicase with different stabilities. We purified two forms of replicase expressed in Escherichia coli, which differed in their thermal stability but showed identical RNA replication activity. Furthermore, we found that the replicase undergoes conversion between these forms due to oxidation, and the Cys-533 residue in the catalytic β subunit and Cys-82 residue in the EF-Tu subunit of the replicase are essential prerequisites for this conversion to occur. These results strongly suggest that the thermal stable replicase contains the intersubunit disulfide bond between these cysteines. The established strategies for isolating and purifying a thermally stable replicase should increase the usefulness of Qβ replicase in various applications, and the data regarding thermal stability obtained in this study may yield insight into the precise mechanism of infection by bacteriophage Qβ.