Correction: genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity

Diversity and conservation of the genome architecture of phages infecting the Alphaproteobacteria

IMPORTANCE

This study reports the results of the largest analysis of genome sequences from phages that infect the Alphaproteobacteria class of bacterial hosts. We analyzed over 100 whole genome sequences of phages to construct dotplots, categorize them into genetically distinct clusters, generate a bootstrapped phylogenetic tree, compute protein orthologs, and predict packaging strategies. We determined that the phage sequences primarily cluster by the bacterial host family, phage morphotype, and genome size. We expect that the findings reported in this seminal study will facilitate future analyses that will improve our knowledge of the phages that infect these hosts.

Genomic diversity and comprehensive taxonomical classification of 61 Bacillus subtilis group member infecting bacteriophages, and the identification of ortholog taxonomic signature genes

BACKGROUND

Despite the applications of Bacillus subtilis group species in various sectors, limited information is available regarding their phages. Here, 61 B. subtilis group species-infecting phages (BSPs) were studied for their taxonomic classification considering the genome-size, genomic diversity, and the host, followed by the identification of orthologs taxonomic signature genes.

RESULTS

BSPs have widely ranging genome sizes that can be bunched into groups to demonstrate correlations to family and subfamily classifications. Comparative analysis re-confirmed the existing, BSPs-containing 14 genera and 21 species and displayed inter-genera similarities within existing subfamilies. Importantly, it also revealed the need for the creation of new taxonomic classifications, including 28 species, nine genera, and two subfamilies (New subfamily1 and New subfamily2) to accommodate inter-genera relatedness. Following pangenome analysis, no ortholog shared by all BSPs was identified, while orthologs, namely, the tail fibers/spike proteins and poly-gamma-glutamate hydrolase, that are shared by more than two-thirds of the BSPs were identified. More importantly, major capsid protein (MCP) type I, MCP type II, MCP type III and peptidoglycan binding proteins that are distinctive orthologs for Herelleviridae, Salasmaviridae, New subfamily1, and New subfamily2, respectively, were identified and analyzed which could serve as signatures to distinguish BSP members of the respective taxon.

CONCLUSIONS

In this study, we show the genomic diversity and propose a comprehensive classification of 61 BSPs, including the proposition for the creation of two new subfamilies, followed by the identification of orthologs taxonomic signature genes, potentially contributing to phage taxonomy.

Comparative Genomics of Six Lytic Bacillus subtilis Phages from the Southwest United States

Background

Despite their importance to microbial dynamics involving Bacillus subtilis, we have a limited understanding of the diversity of phages that can lyse this model organism.

Materials and Methods

Phages were isolated from soil samples collected from various sites in the southwest U.S. deserts on a wild B. subtilis strain. Their genomes were assembled, characterized, and bioinformatically compared.

Results

Six Siphoviruses with high nucleotide and amino acid similarity to each other (>80%) but very limited similarity to phages currently in GenBank were isolated. These phages have double-stranded DNA genomes (55,312 to 56,127 bp) with 86-91 putative protein coding genes, and a low GC content. Comparative genomics reveal differences in loci encoding proteins that are putatively involved in bacterial adsorption with evidence for genomic mosaicism and a possible role for small genes.

Conclusions

A comparative approach provides insights into phage evolution, including the role of indels in protein folding.

Isolation, characterization, and comparative genomic analysis of vB_PlaM_Pd22F, a new bacteriophage of the family Myoviridae

The use of phage and phage-based products for the prevention and treatment of bee disease is one of the promising natural alternatives to chemical or antibiotic treatments in beekeeping. A novel lysogenic bacteriophage, phage Pd22F (vB_PlaM_Pd22F), was isolated from Paenibacillus dendritiformis by the prophage induction method. This phage, which is capable of infecting Paenibacillus larvae and P. dendritiformis strains, was characterized by microbiological and comparative genomic analysis. Transmission electron microscopy images showed that phage Pd22F had the morphology of a myovirus. Whole-genome sequencing results showed that vB_Pla M_Pd22F has an 86,388-bp linear dsDNA genome with a GC content of 50.68%. This genome has 124 coding sequences (CDSs), 53% of which encode functionally unknown proteins and 57 of which encode proteins that show similarity to known proteins. In addition, one tRNA gene was found. The phage Pd22F genome does not contain any antimicrobial resistance genes. The similarity between the genome sequence of phage Pd22F and the whole genome sequences of other Paenibacillus phages available in the NCBI Virus Database was found to be below 50% (42%), indicating that phage Pd22F differs greatly from previously characterized phages at the DNA level. The results of comparative genomics and phylogenetic analysis revealed that Pd22F is a new phage belonging to the family Myoviridae, order Caudovirales. This is the first report of genomic and morphological characterization of a Paenibacillus dendritiformis prophage.

Transkingdom Analysis of the Female Reproductive Tract Reveals Bacteriophages form Communities

The female reproductive tract (FRT) microbiome plays a vital role in maintaining vaginal health. Viruses are key regulators of other microbial ecosystems, but little is known about how the FRT viruses (virome), particularly bacteriophages that comprise the phageome, impact FRT health and dysbiosis. We hypothesize that bacterial vaginosis (BV) is associated with altered FRT phageome diversity, transkingdom interplay, and bacteriophage discriminate taxa. Here, we conducted a retrospective, longitudinal analysis of vaginal swabs collected from 54 BV-positive and 46 BV-negative South African women. Bacteriome analysis revealed samples clustered into five distinct bacterial community groups (CGs), and further, bacterial alpha diversity was significantly associated with BV. Virome analysis on a subset of baseline samples showed FRT bacteriophages clustering into novel viral state types (VSTs), a viral community clustering system based on virome composition and abundance. Distinct BV bacteriophage signatures included increased alpha diversity along with discriminant Bacillus, Burkholderia, and Escherichia bacteriophages. Bacteriophage-bacteria transkingdom associations were also identified between Bacillus and Burkholderia viruses and BV-associated bacteria, providing key insights for future studies elucidating the transkingdom interactions driving BV-associated microbiome perturbations. In this cohort, bacteriophage-bacterial associations suggest complex interactions, which may play a role in the establishment and maintenance of BV.

Isolation and Characterisation of the Bundooravirus Genus and Phylogenetic Investigation of the Salasmaviridae Bacteriophages

Bacillus is a highly diverse genus containing over 200 species that can be problematic in both industrial and medical settings. This is mainly attributed to Bacillus sp. being intrinsically resistant to an array of antimicrobial compounds, hence alternative treatment options are needed. In this study, two bacteriophages, PumA1 and PumA2 were isolated and characterized. Genome nucleotide analysis identified the two phages as novel at the DNA sequence level but contained proteins similar to phi29 and other related phages. Whole genome phylogenetic investigation of 34 phi29-like phages resulted in the formation of seven clusters that aligned with recent ICTV classifications. PumA1 and PumA2 share high genetic mosaicism and form a genus with another phage named WhyPhy, more recently isolated from the United States of America. The three phages within this cluster are the only candidates to infect B. pumilus. Sequence analysis of B. pumilus phage resistant mutants revealed that PumA1 and PumA2 require polymerized and peptidoglycan bound wall teichoic acid (WTA) for their infection. Bacteriophage classification is continuously evolving with the increasing phages' sequences in public databases. Understanding phage evolution by utilizing a combination of phylogenetic approaches provides invaluable information as phages become legitimate alternatives in both human health and industrial processes.

Insights into the genome architecture and evolution of Shiga toxin encoding bacteriophages of Escherichia coli

Background

A total of 179 Shiga toxin-producing Escherichia coli (STEC) complete genomes were analyzed in terms of serotypes, prophage coding regions, and stx gene variants and their distribution. We further examined the genetic diversity of Stx-converting phage genomes (Stx phages), focusing on the lysis-lysogeny decision and lytic cassettes.

Results

We show that most STEC isolates belong to non-O157 serotypes (73 %), regardless the sources and geographical regions. While the majority of STEC genomes contain a single stx gene (61 %), strains containing two (35 %), three (3 %) and four (1 %) stx genes were also found, being stx2 the most prevalent gene variant. Their location is exclusively found in intact prophage regions, indicating that they are phage-borne. We further demonstrate that Stx phages can be grouped into four clusters (A, B, C and D), three subclusters (A1, A2 and A3) and one singleton, based on their shared gene content. This cluster distribution is in good agreement with their predicted virion morphologies. Stx phage genomes are highly diverse with a vast number of 1,838 gene phamilies (phams) of related sequences (of which 677 are orphams i.e. unique genes) and, although having high mosaicism, they are generally organized into three major transcripts. While the mechanisms that guide lysis–lysogeny decision are complex, there is a strong selective pressure to maintain the stx genes location close to the lytic cassette composed of predicted SAR-endolysin and pin-holin lytic proteins. The evolution of STEC Stx phages seems to be strongly related to acquiring genetic material, probably from horizontal gene transfer events.

Conclusions

This work provides novel insights on the genetic structure of Stx phages, showing a high genetic diversity throughout the genomes, where the various lysis-lysogeny regulatory systems are in contrast with an uncommon, but conserved, lytic system always adjacent to stx genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07685-0.

Two novel bacteriophage genera from a groundwater reservoir highlight subsurface environments as underexplored biotopes in bacteriophage ecology

Although bacteriophages are central entities in bacterial ecology and population dynamics, there is currently no literature on the genomes of bacteriophages isolated from groundwater. Using a collection of bacterial isolates from an aquifer as hosts, this study isolated, sequenced and characterised two bacteriophages native to the groundwater reservoir. Host phylogenetic analyses revealed that the phages targeted B. mycoides and a novel Pseudomonas species. These results suggest that both bacteriophages represent new genera, highlighting that groundwater reservoirs, and probably other subsurface environments as well, are underexplored biotopes in terms of the presence and ecology of bacteriophages.

Genome Sequences of Two Bacillus Phages Isolated from Indonesia

Here, the genome sequences of two soil bacteriophages isolated from a red chili plantation in Indonesia are presented. The genome of vB_BspS_SplendidRed (42,859 bp) is highly similar to Bacillus phage Ray17 from the United States, while vB_BspM_MarvelLand (156,945 bp) is highly similar to Bacillus phage BC01 from South Korea.

Here, the genome sequences of two soil bacteriophages isolated from a red chili plantation in Indonesia are presented. The genome of vB_BspS_SplendidRed (42,859 bp) is highly similar to Bacillus phage Ray17 from the United States, while vB_BspM_MarvelLand (156,945 bp) is highly similar to Bacillus phage BC01 from South Korea.

Characterization and comparative genomic analysis of virulent and temperate Bacillus megaterium bacteriophages

Next-Generation Sequencing (NGS) technologies provide unique possibilities for the comprehensive assessment of the environmental diversity of bacteriophages. Several Bacillus bacteriophages have been isolated, but very few Bacillus megaterium bacteriophages have been characterized. In this study, we describe the biological characteristics, whole genome sequences, and annotations for two new isolates of the B. megaterium bacteriophages (BM5 and BM10), which were isolated from Egyptian soil samples. Growth analyses indicated that the phages BM5 and BM10 have a shorter latent period (25 and 30 min, respectively) and a smaller burst size (103 and 117 PFU, respectively), in comparison to what is typical for Bacillus phages. The genome sizes of the phages BM5 and BM10 were 165,031 bp and 165,213 bp, respectively, with modular organization. Bioinformatic analyses of these genomes enabled the assignment of putative functions to 97 and 65 putative ORFs, respectively. Comparative analysis of the BM5 and BM10 genome structures, in conjunction with other B. megaterium bacteriophages, revealed relatively high levels of sequence and organizational identity. Both genomic comparisons and phylogenetic analyses support the conclusion that the sequenced phages (BM5 and BM10) belong to different sub-clusters (L5 and L7, respectively), within the L-cluster, and display different lifestyles (lysogenic and lytic, respectively). Moreover, sequenced phages encode proteins associated with Bacillus pathogenesis. In addition, BM5 does not contain any tRNA sequences, whereas BM10 genome codes for 17 tRNAs.

A CRISPR-Cas9-Based Toolkit for Fast and Precise In Vivo Genetic Engineering of Bacillus subtilis Phages

Phages are currently under discussion as a solution for the antibiotic crisis, as they may cure diseases caused by multi-drug-resistant pathogens. However, knowledge of phage biology and genetics is limited, which impedes risk assessment of therapeutic applications. In order to enable advances in phage genetic research, the aim of this work was to create a toolkit for simple and fast genetic engineering of phages recruiting Bacillus subtilis as host system. The model organism B. subtilis represents a non-pathogenic surrogate of its harmful relatives, such as Bacillus anthracis or Bacillus cereus. This toolkit comprises the application CutSPR, a bioinformatic tool for rapid primer design, and facilitates the cloning of specific CRISPR-Cas9-based mutagenesis plasmids. The employment of the prophage-free and super-competent B. subtilis TS01 strain enables an easy and fast introduction of specific constructs for in vivo phage mutagenesis. Clean gene deletions and a functional clean gene insertion into the genome of the model phage vB_BsuP-Goe1 served as proof of concept and demonstrate reliability and high efficiency. The here presented toolkit allows comprehensive investigation of the diverse phage genetic pool, a better understanding of phage biology, and safe phage applications.

Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity

The bacteriophage population is large, dynamic, ancient, and genetically diverse. Limited genomic information shows that phage genomes are mosaic, and the genetic architecture of phage populations remains ill-defined. To understand the population structure of phages infecting a single host strain, we isolated, sequenced, and compared 627 phages of Mycobacterium smegmatis. Their genetic diversity is considerable, and there are 28 distinct genomic types (clusters) with related nucleotide sequences. However, amino acid sequence comparisons show pervasive genomic mosaicism, and quantification of inter-cluster and intra-cluster relatedness reveals a continuum of genetic diversity, albeit with uneven representation of different phages. Furthermore, rarefaction analysis shows that the mycobacteriophage population is not closed, and there is a constant influx of genes from other sources. Phage isolation and analysis was performed by a large consortium of academic institutions, illustrating the substantial benefits of a disseminated, structured program involving large numbers of freshman undergraduates in scientific discovery.

DOI:http://dx.doi.org/10.7554/eLife.06416.001

Viruses are unable to replicate independently. To generate copies of itself, a virus must instead invade a target cell and commandeer that cell's replication machinery. Different viruses are able to invade different types of cell, and a group of viruses known as bacteriophages (or phages for short) replicate within bacteria. The enormous number and diversity of phages in the world means that they play an important role in virtually every ecosystem.

Despite their importance, relatively little is known about how different phage populations are related to each other and how they evolved. Many phages contain their genetic information in the form of strands of DNA. Using genetic sequencing to find out where and how different genes are encoded in the DNA can reveal information about how different viruses are related to each other. These relationships are particularly complicated in phages, as they can exchange genes with other viruses and microbes.

Previous studies comparing the genomes—the complete DNA sequence—of reasonably small numbers of phages that infect the Mycobacterium group of bacteria have found that the phages can be sorted into ‘clusters’ based on similarities in their genes and where these are encoded in their DNA. However, the number of phages investigated so far has been too small to conclude how different clusters are related. Are the clusters separate, or do they form a ‘continuum’ with different genes and DNA sequences shared between different clusters?

Here, Pope, Bowman, Russell et al. compare the individual genomes of 627 bacteriophages that infect the bacterial species Mycobacterium smegmatis. This is by far the largest number of phage genomes analyzed from a single host species. The large number of genomes analyzed allowed a much clearer understanding of the complexity and diversity of these phages to be obtained. The isolation, sequencing and analysis of the hundreds of M. smegmatis bacteriophage genomes was performed by an integrated research and education program, called the Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) program. This enabled thousands of undergraduate students from different institutions to contribute to the phage discovery and sequencing project, and co-author the report. SEA-PHAGES therefore shows that it is possible to successfully incorporate genuine scientific research into an undergraduate course, and that doing so can benefit both the students and researchers involved.

The results show that while the genomes could be categorized into 28 clusters, the genomes are not completely unrelated. Instead, a spread of diversity is seen, as genes and groups of genes are shared between different clusters. Pope, Bowman, Russell et al. further reveal that the phage population is in a constant state of change, and continuously acquires genes from other microorganisms and viruses.

DOI:http://dx.doi.org/10.7554/eLife.06416.002