High-throughput identification of viral termini and packaging mechanisms in virome datasets using PhageTermVirome

Enrichment, Sequencing, and Identification of DNA Bacteriophages from Fecal Samples

Research on individual viruses and phages, as well as viral populations (viromes), is greatly expanding. Phages and viromes are increasingly suspected to have numerous impacts on the ecosystem in which they reside by interacting directly or indirectly with the other organisms present in their environment. In particular, phage communities of the gut microbiota have been associated with a wide range of diseases. However, properly investigating intestinal viromes is still very challenging, both experimentally and analytically. This chapter proposes a simple and reproducible protocol to separate and enrich DNA phage particles from fecal samples, to sequence them, and finally obtain a basic but robust bioinformatic characterization and classification of the global bacteriophage community.

Characterization of 29 newly isolated bacteriophages as a potential therapeutic agent against IMP-6-producing Klebsiella pneumoniae from clinical specimens

Carbapenemase-producing Enterobacteriaceae (CPE) are one of the most detrimental species of antibiotic-resistant bacteria globally. Phage therapy has emerged as an effective strategy for the treatment of CPE infections. In western Japan, the rise of Klebsiella pneumoniae strains harboring the pKPI-6 plasmid encoding bla IMP-6 is of increasing concern. To address this challenge, we isolated 29 phages from Japanese sewage, specifically targeting 31 K. pneumoniae strains and one Escherichia coli strain harboring the pKPI-6 plasmid. Electron microscopy analysis revealed that among the 29 isolated phages, 21 (72.4%), 5 (17.2%), and 3 (10.3%) phages belonged to myovirus, siphovirus, and podovirus morphotypes, respectively. Host range analysis showed that 18 Slopekvirus strains within the isolated phages infected 25-26 K. pneumoniae strains, indicating that most of the isolated phages have a broad host range. Notably, K. pneumoniae strain Kp21 was exclusively susceptible to phage øKp_21, whereas Kp22 exhibited susceptibility to over 20 phages. Upon administering a phage cocktail composed of 10 phages, we observed delayed emergence of phage-resistant bacteria in Kp21 but not in Kp22. Intriguingly, phage-resistant Kp21 exhibited heightened sensitivity to other bacteriophages, indicating a "trade-off" for resistance to phage øKp_21. Our proposed phage set has an adequate number of phages to combat the K. pneumoniae strain prevalent in Japan, underscoring the potential of a well-designed phage cocktail in mitigating the occurrence of phage-resistant bacteria. IMPORTANCE The emergence of Klebsiella pneumoniae harboring the bla IMP-6 plasmid poses an escalating threat in Japan. In this study, we found 29 newly isolated bacteriophages that infect K. pneumoniae strains carrying the pKPI-6 plasmid from clinical settings in western Japan. Our phages exhibited a broad host range. We applied a phage cocktail treatment composed of 10 phages against two host strains, Kp21 and Kp22, which displayed varying phage susceptibility patterns. Although the phage cocktail delayed the emergence of phage-resistant Kp21, it was unable to hinder the emergence of phage-resistant Kp22. Moreover, the phage-resistant Kp21 became sensitive to other phages that were originally non-infective to the wild-type Kp21 strains. Our study highlights the potential of a well-tailored phage cocktail in reducing the occurrence of phage-resistant bacteria.

Combined direct/indirect detection allows identification of DNA termini in diverse sequencing datasets and supports a multiple-initiation-site model for HIV plus-strand synthesis

Replication of genetic material involves the creation of characteristic termini. Determining these termini is important to refine our understanding of the mechanisms involved in maintaining the genomes of cellular organisms and viruses. Here we describe a computational approach combining direct and indirect readouts to detect termini from next-generation short-read sequencing. While a direct inference of termini can come from mapping the most prominent start positions of captured DNA fragments, this approach is insufficient in cases where the DNA termini are not captured, whether for biological or technical reasons. Thus, a complementary (indirect) approach to terminus detection can be applied, taking advantage of the imbalance in coverage between forward and reverse sequence reads near termini. A resulting metric ("strand bias") can be used to detect termini even where termini are naturally blocked from capture or ends are not captured during library preparation (e.g., in tagmentation-based protocols). Applying this analysis to datasets where known DNA termini are present, such as from linear double-stranded viral genomes, yielded distinct strand bias signals corresponding to these termini. To evaluate the potential to analyze a more complex situation, we applied the analysis to examine DNA termini present early after HIV infection in a cell culture model. We observed both the known termini expected based on standard models of HIV reverse transcription (the U5-right-end and U3-left-end termini) as well as a signal corresponding to a previously described additional initiation site for plus-strand synthesis (cPPT [central polypurine tract]). Interestingly, we also detected putative terminus signals at additional sites. The strongest of these are a set that share several characteristics with the previously characterized plus-strand initiation sites (the cPPT and 3' PPT [polypurine tract] sites): (i) an observed spike in directly captured cDNA ends, an indirect terminus signal evident in localized strand bias, (iii) a preference for location on the plus-strand, (iv) an upstream purine-rich motif, and (v) a decrease in terminus signal at late time points after infection. These characteristics are consistent in duplicate samples in two different genotypes (wild type and integrase-lacking HIV). The observation of distinct internal termini associated with multiple purine-rich regions raises a possibility that multiple internal initiations of plus-strand synthesis might contribute to HIV replication.

Chromosome folding and prophage activation reveal specific genomic architecture for intestinal bacteria

BACKGROUND

Bacteria and their viruses, bacteriophages, are the most abundant entities of the gut microbiota, a complex community of microorganisms associated with human health and disease. In this ecosystem, the interactions between these two key components are still largely unknown. In particular, the impact of the gut environment on bacteria and their associated prophages is yet to be deciphered.

RESULTS

To gain insight into the activity of lysogenic bacteriophages within the context of their host genomes, we performed proximity ligation-based sequencing (Hi-C) in both in vitro and in vivo conditions on the 12 bacterial strains of the OMM¹² synthetic bacterial community stably associated within mice gut (gnotobiotic mouse line OMM¹²). High-resolution contact maps of the chromosome 3D organization of the bacterial genomes revealed a wide diversity of architectures, differences between environments, and an overall stability over time in the gut of mice. The DNA contacts pointed at 3D signatures of prophages leading to 16 of them being predicted as functional. We also identified circularization signals and observed different 3D patterns between in vitro and in vivo conditions. Concurrent virome analysis showed that 11 of these prophages produced viral particles and that OMM¹² mice do not carry other intestinal viruses.

CONCLUSIONS

The precise identification by Hi-C of functional and active prophages within bacterial communities will unlock the study of interactions between bacteriophages and bacteria across conditions (healthy vs disease). Video Abstract.