Comparative whole genome analysis of six diagnostic brucellaphages

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.

Comprehensive deciphering prophages in genus Acetobacter on the ecology, genomic features, toxin–antitoxin system, and linkage with CRISPR-Cas system

Acetobacter is the predominant microbe in vinegar production, particularly in those natural fermentations that are achieved by complex microbial communities. Co-evolution of prophages with Acetobacter, including integration, release, and dissemination, heavily affects the genome stability and production performance of industrial strains. However, little has been discussed yet about prophages in Acetobacter. Here, prophage prediction analysis using 148 available genomes from 34 Acetobacter species was carried out. In addition, the type II toxin–antitoxin systems (TAs) and CRISPR-Cas systems encoded by prophages or the chromosome were analyzed. Totally, 12,000 prophage fragments were found, of which 350 putatively active prophages were identified in 86.5% of the selected genomes. Most of the active prophages (83.4%) belonged to the order Caudovirales dominated by the families Siphoviridae and Myroviridae prophages (71.4%). Notably, Acetobacter strains survived in complex environments that frequently carried multiple prophages compared with that in restricted habits. Acetobacter prophages showed high genome diversity and horizontal gene transfer across different bacterial species by genomic feature characterization, average nucleotide identity (ANI), and gene structure visualization analyses. About 31.14% of prophages carry type II TAS, suggesting its important role in addiction, bacterial defense, and growth-associated bioprocesses to prophages and hosts. Intriguingly, the genes coding for Cse1, Cse2, Cse3, Cse4, and Cas5e involved in type I-E and Csy4 involved in type I-F CRISPR arrays were firstly found in two prophages. Type II-C CRISPR-Cas system existed only in Acetobacter aceti, while the other Acetobacter species harbored the intact or eroded type I CRISPR-Cas systems. Totally, the results of this study provide fundamental clues for future studies on the role of prophages in the cell physiology and environmental behavior of Acetobacter.

A Proof of Principle for the Detection of Viable Brucella spp. in Raw Milk by qPCR Targeting Bacteriophages

Brucellosis is still a global health issue, and surveillance and control of this zoonotic disease in livestock remains a challenge. Human outbreaks are mainly linked to the consumption of unpasteurized dairy products. The detection of human pathogenic Brucella species in food of animal origin is time-consuming and laborious. Bacteriophages are broadly applied to the typing of Brucella isolates from pure culture. Since phages intracellularly replicate to very high numbers, they can also be used as specific indicator organisms of their host bacteria. We developed a novel real-time PCR (qPCR) assay targeting the highly conserved helicase sequence harbored in all currently known Brucella-specific lytic phages. Quality and performance tests determined a limit of detection of <1 genomic copy/µL. In raw milk artificially contaminated with Brucella microti, Iz_v phages were reliably detected after 39 h of incubation, indicating the presence of viable bacteria. The qPCR assay showed high stability in the milk matrix and significantly shortened the time to diagnosis when compared to traditional culture-based techniques. Hence, our molecular assay is a reliable and sensitive method to analyze phage titers, may help to reduce the hands-on time needed for the screening of potentially contaminated food, and reveals infection risks without bacterial isolation.

Complete Genome Sequence of Brucella abortus Phage EF4, Determined Using Long-Read Sequencing

Brucellaphage EF4 was isolated from elk feces. The 38,321-bp double-stranded DNA genome is predicted to contain 72 coding regions, 38 of which have been assigned predicted functions. This phage displays nucleotide similarity to other brucellaphages of the genus Perisivirus.

Brucellaphage EF4 was isolated from elk feces. The 38,321-bp double-stranded DNA genome is predicted to contain 72 coding regions, 38 of which have been assigned predicted functions. This phage displays nucleotide similarity to other brucellaphages of the genus Perisivirus.

Correlation of Host Range Expansion of Therapeutic Bacteriophage Sb-1 with Allele State at a Hypervariable Repeat Locus

Staphylococci are frequent agents of health care-associated infections and include methicillin-resistant Staphylococcus aureus (MRSA), which is resistant to first-line antibiotic treatments. Bacteriophage (phage) therapy is a promising alternative antibacterial option to treat MRSA infections. S. aureus-specific phage Sb-1 has been widely used in Georgia to treat a variety of human S. aureus infections. Sb-1 has a broad host range within S. aureus, including MRSA strains, and its host range can be further expanded by adaptation to previously resistant clinical isolates. The susceptibilities of a panel of 25 genetically diverse clinical MRSA isolates to Sb-1 phage were tested, and the phage had lytic activity against 23 strains (92%). The adapted phage stock (designated Sb-1A) was tested in comparison with the parental phage (designated Sb-1P). Sb-1P had lytic activity against 78/90 strains (87%) in an expanded panel of diverse global S. aureus isolates, while eight additional strains in this panel were susceptible to Sb-1A (lytic against 86/90 strains [96%]). The Sb-1A stock was shown to be a mixed population of phage clones, including approximately 4% expanded host range mutants, designated Sb-1M. In an effort to better understand the genetic basis for this host range expansion, we sequenced the complete genomes of the parental Sb-1P and two Sb-1M mutants. Comparative genomic analysis revealed a hypervariable complex repeat structure in the Sb-1 genome that had a distinct allele that correlated with the host range expansion. This hypervariable region was previously uncharacterized in Twort-like phages and represents a novel putative host range determinant.IMPORTANCE Because of limited therapeutic options, infections caused by methicillin-resistant Staphylococcus aureus represent a serious problem in both civilian and military health care settings. Phages have potential as alternative antibacterial agents that can be used in combination with antibiotic drugs. For decades, phage Sb-1 has been used in former Soviet Union countries for antistaphylococcal treatment in humans. The therapeutic spectrum of activity of Sb-1 can be increased by selecting mutants of the phage with expanded host ranges. In this work, the host range of phage Sb-1 was expanded in the laboratory, and a hypervariable region in its genome was identified with a distinct allele state that correlated with this host range expansion. These results provide a genetic basis for better understanding the mechanisms of phage host range expansion.

Genomic characterization of three novel Basilisk-like phages infecting Bacillus anthracis

Background

In the present study, we sequenced the complete genomes of three novel bacteriophages v_B-Bak1, v_B-Bak6, v_B-Bak10 previously isolated from historical anthrax burial sites in the South Caucasus country of Georgia. We report here major trends in the molecular evolution of these phages, which we designate as “Basilisk-Like-Phages” (BLPs), and illustrate patterns in their evolution, genomic plasticity and core genome architecture.

Results

Comparative whole genome sequence analysis revealed a close evolutionary relationship between our phages and two unclassified Bacillus cereus group phages, phage Basilisk, a broad host range phage (Grose JH et al., J Vir. 2014;88(20):11846-11860) and phage PBC4, a highly host-restricted phage and close relative of Basilisk (Na H. et al. FEMS Microbiol. letters. 2016;363(12)). Genome comparisons of phages v_B-Bak1, v_B-Bak6, and v_B-Bak10 revealed significant similarity in sequence, gene content, and synteny with both Basilisk and PBC4. Transmission electron microscopy (TEM) confirmed the three phages belong to the Siphoviridae family. In contrast to the broad host range of phage Basilisk and the single-strain specificity of PBC4, our three phages displayed host specificity for Bacillus anthracis. Bacillus species including Bacillus cereus, Bacillus subtilis, Bacillus anthracoides, and Bacillus megaterium were refractory to infection.

Conclusions

Data reported here provide further insight into the shared genomic architecture, host range specificity, and molecular evolution of these rare B. cereus group phages. To date, the three phages represent the only known close relatives of the Basilisk and PBC4 phages and their shared genetic attributes and unique host specificity for B. anthracis provides additional insight into candidate host range determinants.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5056-4) contains supplementary material, which is available to authorized users.

Highly Sensitive Bacteriophage-Based Detection of Brucella abortus in Mixed Culture and Spiked Blood

For decades, bacteriophages (phages) have been used for Brucella species identification in the diagnosis and epidemiology of brucellosis. Traditional Brucella phage typing is a multi-day procedure including the isolation of a pure culture, a step that can take up to three weeks. In this study, we focused on the use of brucellaphages for sensitive detection of the pathogen in clinical and other complex samples, and developed an indirect method of Brucella detection using real-time quantitative PCR monitoring of brucellaphage DNA amplification via replication on live Brucella cells. This assay allowed the detection of single bacteria (down to 1 colony-forming unit per milliliter) within 72 h without DNA extraction and purification steps. The technique was equally efficient with Brucella abortus pure culture and with mixed cultures of B. abortus and α-proteobacterial near neighbors that can be misidentified as Brucella spp., Ochrobactrum anthropi and Afipia felis. The addition of a simple short sample preparation step enabled the indirect phage-based detection of B. abortus in spiked blood, with the same high sensitivity. This indirect phage-based detection assay enables the rapid and sensitive detection of live B. abortus in mixed cultures and in blood samples, and can potentially be applied for detection in other clinical samples and other complex sample types.

Genetic Diversity of Brucella Reference and Non-reference Phages and Its Impact on Brucella-Typing

Virulent phages have been used for many years to type Brucella isolates, but until recently knowledge about the genetic makeup of these phages remains limited. In this work the host specificity and genomic sequences of the original set (deposited in 1960) of VLA Brucella reference phages Tb, Fi, Wb, Bk2, R/C, and Iz were analyzed and compared with hitherto described brucellaphages. VLA phages turned out to be different from homonymous phages in other laboratories. The host range of the phages was defined by performing plaque assays with a wide selection of Brucella strains. Propagation of the phages on different strains did not alter host specificity. Sequencing of the phages TbV, FiV, WbV, and R/CV revealed nucleotide variations when compared to same-named phages previously described by other laboratories. The phages Bk2V and IzV were sequenced for the first time. While Bk2V exhibited the same deletions as WbV, IzV possesses the largest genome of all Brucella reference phages. The duplication of a 301 bp sequence in this phage and the large deletion in Bk2V, WbV, and R/CV may be a result of recombination caused by repetitive sequences located in this DNA region. To identify new phages as potential candidates for lysotyping, the host range and Single Nucleotide Polymorphisms (SNPs) of 22 non-reference Brucella phages were determined. The phages showed lysis patterns different from those of the reference phages and thus represent novel valuable candidates in the typing set.

Analysis of the First Temperate Broad Host Range Brucellaphage (BiPBO1) Isolated from B. inopinata

Brucella species are important human and animal pathogens. Though, only little is known about mobile genetic elements of these highly pathogenic bacteria. To date, neither plasmids nor temperate phages have been described in brucellae. We analyzed genomic sequences of various reference and type strains and identified a number of putative prophages residing within the Brucella chromosomes. By induction, phage BiPBO1 was isolated from Brucella inopinata. BiPBO1 is a siphovirus that infects several Brucella species including Brucella abortus and Brucella melitensis. Integration of the phage genome occurs adjacent to a tRNA gene in chromosome 1 (chr 1). The bacterial (attB) and phage (attP) attachment sites comprise an identical sequence of 46 bp. This sequence exists in many Brucella and Ochrobactrum species. The BiPBO1 genome is composed of a 46,877 bp double-stranded DNA. Eighty-seven putative gene products were determined, of which 32 could be functionally assigned. Strongest similarities were found to a temperate phage residing in the chromosome of Ochrobactrum anthropi ATCC 49188 and to prophages identified in several families belonging to the order rhizobiales. The data suggest that horizontal gene transfer may occur between Brucella and Ochrobactrum and underpin the close relationship of these environmental and pathogenic bacteria.

On genome annotation of Brucellaphage Gadvasu (BpG): discovery of ORFans for integrated systems biology approaches

Brucellaphage Gadvasu (BpG) is a lytic phage infecting Brucella spp. Brucellaphages contain dsDNA as genetic material and are short-tailed particles with host-specificity. Here, we report the challenges on annotation in the complete genome sequence of BpG when compared with that of a recent broad host-range brucellaphage Pr, an original reference genome. The extracted DNA was subjected to genome sequencing with Illumina technology and assembled using SSAKE/Velvet. A significant number of genes were found to be similar between the phages with sequence analysis revealing conserved open reading frames that correspond to 33 gene ontology classifiers, transcriptional terminators and a few putative transcriptional promoters. The analyses revealed that the genome constitutes 1269 contigs and 275 genes encoding 260 proteins. The sequence comparison from the reference data indicated that the genome shares an approximately 70 % nucleotide similarity and differs mainly in the region encoding proteins. We bring this commentary providing an overview of how this exemplar genome can allow us to understand these known unknown regions in brucellaphages.

Whole genome sequence comparison of ten diagnostic brucellaphages propagated on two Brucella abortus hosts

Background

Recently the genome sequences of two brucellaphages, isolated in Georgia (Tb) and Mexico (Pr) were reported revealing pronounced sequence homogeneity and the presence of two major indels discriminating the two phages. Subsequent genome sequencing of six diagnostic brucellaphages: Tbilisi (Tb), Firenze (Fz), Weybridge (Wb), S708, Berkeley (Bk) and R/C phages identified three major genetic groups. However, the propensity for fine-scale genetic variability of diverse brucellaphages grown on multiple hosts within a single Brucella species remains unknown.

Methods

We sequenced the complete genomes of ten brucellaphages following initial propagation on B. abortus strain 141 and after subsequent propagation on B. abortus strain S19. All brucellaphages were isolated and propagated at the Eliava Institute including the original Tb phage. Genomic libraries were quantified using the Qbit and sheared on the Covaris M220. QC for fragmentation was performed on the BioAnalyzer 2100. DNA libraries were prepared using an Illumina Paired-End protocol and sequenced on the Illumina MiSeq. Sequence analysis was performed using Geneious and MEGA software.

Results

Comparative whole genome sequence analysis revealed genetic homogeneity consistent with previously published data as well as multiple nucleotide variations. Genomic changes as a result of passages were observed in similar genes and predominantly occurred at identical sites in separate phages. Multiple instances of within-sample genetic heterogeneity were observed often at shared genomics positions across phages. Positive selection was detected in the tail collar protein gene. We also identified a Staphylothermus marinus F1-like CRISPR spacer and sequences orthologous to both prophage antirepressors of Brucella spp. and intergenic sequences encoded by Ochrobactrum anthropi.

Conclusion

We surveyed whole genome level diversity in phage lytic for B. abortus as they are propagated on alternate vaccine strains within the species. Our data extend previous results indicating select variable hotspots and broad genomic homogeneity as well as multiple common polymorphisms and within-sample variation. These data also provide additional genomes for future reference in comparative studies involving the molecular evolution and host specificity of brucellaphages.

Isolation and Genomic Characterization of the T4-Like Bacteriophage PM2 Infecting Pectobacterium carotovorum subsp. carotovorum

In order to control Pectobacterium carotovorum subsp. carotovorum, a novel virulent bacteriophage PM2 was isolated. Bacteriophage PM2 can infect 48% of P. carotovorum subsp. carotovorum and 78% of P. carotovorum subsp. brasilliensis but none of atrosepticum, betavasculorum, odoriferum and wasabiae isolates had been infected with PM2. PM2 phage belongs to the family Myoviridae, and contains a large head and contractile tail. It has a 170,286 base pair genome that encodes 291 open reading frames (ORFs) and 12 tRNAs. Most ORFs in bacteriophage PM2 share a high level of homology with T4-like phages including IME08, RB69, and JS98. Phylogenetic analysis based on the amino acid sequence of terminase large subunits confirmed that PM2 is classified as a T4-like phage. It contains no integrase- or no repressor-coding genes related to the lysogenic cycle, and lifestyle prediction using PHACT software suggested that PM2 is a virulent bacteriophage.