Divergence and mosaicism among virulent soil phages of the Burkholderia cepacia complex

Description of two novel bacteriophages of the class Caudoviricetes that infect Ralstonia solanacearum and Ralstonia pseudosolanacearum

Here, we describe the isolation and characterization of two novel phages from Brazilian soil that infect Ralstonia solanacearum and Ralstonia pseudosolanacearum, which we have named "RS phage AB1 and RS phage CA1. Genome sequencing and phylogenetic analysis revealed that RS phage AB1 is a novel member of the family Peduoviridae, while RS phage CA1 could not be classified as part of any established family. Thus, we propose a new viral family, "Anamaviridae", with two subfamilies, "Kantovirinae" and "Mascarenevirinae", with the latter including RS phage CA1. We propose the species names "Cocadavirus alagoinhas" and "Acarajevirus bahia" for RS phage CA1 and RS phage AB1, respectively. These findings increase our understanding of the diversity of phages infecting plant pathogens of the genus Ralstonia.

Phage therapy to treat cystic fibrosis Burkholderia cepacia complex lung infections: perspectives and challenges

Burkholderia cepacia complex is a cause of serious lung infections in people with cystic fibrosis, exhibiting extremely high levels of antimicrobial resistance. These infections are difficult to treat and are associated with high morbidity and mortality. With a notable lack of new antibiotic classes currently in development, exploring alternative antimicrobial strategies for Burkholderia cepacia complex is crucial. One potential alternative seeing renewed interest is the use of bacteriophage (phage) therapy. This review summarises what is currently known about Burkholderia cepacia complex in cystic fibrosis, as well as challenges and insights for using phages to treat Burkholderia cepacia complex lung infections.

Comparative Metagenomic Analysis of Bacteriophages and Prophages in Gnotobiotic Mouse Models

Gnotobiotic murine models are important to understand microbiota-host interactions. Despite the role of bacteriophages as drivers for microbiome structure and function, there is no information about the structure and function of the gut virome in gnotobiotic models and the link between bacterial and bacteriophage/prophage diversity. We studied the virome of gnotobiotic murine Oligo-MM12 (12 bacterial species) and reduced Altered Schaedler Flora (ASF, three bacterial species). As reference, the virome of Specific Pathogen-Free (SPF) mice was investigated. A metagenomic approach was used to assess prophages and bacteriophages in the guts of 6-week-old female mice. We identified a positive correlation between bacteria diversity, and bacteriophages and prophages. Caudoviricetes (82.4%) were the most prominent class of phages in all samples with differing relative abundance. However, the host specificity of bacteriophages belonging to class Caudoviricetes differed depending on model bacterial diversity. We further studied the role of bacteriophages in horizontal gene transfer and microbial adaptation to the host's environment. Analysis of mobile genetic elements showed the contribution of bacteriophages to the adaptation of bacterial amino acid metabolism. Overall, our results implicate virome "dark matter" and interactions with the host system as factors for microbial community structure and function which determine host health. Taking the importance of the virome in the microbiome diversity and horizontal gene transfer, reductions in the virome might be an important factor driving losses of microbial biodiversity and the subsequent dysbiosis of the gut microbiome.

Phage Milagro: a platform for engineering a broad host range virulent phage for Burkholderia

IMPORTANCE

Burkholderia infections are a significant concern in people with CF and other immunocompromising disorders, and are difficult to treat with conventional antibiotics due to their inherent drug resistance. Bacteriophages, or bacterial viruses, are now seen as a potential alternative therapy for these infections, but most of the naturally occurring phages are temperate and have narrow host ranges, which limit their utility as therapeutics. Here we describe the temperate Burkholderia phage Milagro and our efforts to engineer this phage into a potential therapeutic by expanding the phage host range and selecting for phage mutants that are strictly virulent. This approach may be used to generate new therapeutic agents for treating intractable infections in CF patients.

Synergistic Interactions among Burkholderia cepacia Complex-Targeting Phages Reveal a Novel Therapeutic Role for Lysogenization-Capable Phages

Antimicrobial resistance is a danger to global public health and threatens many aspects of modern medicine. Bacterial species such as those of the Burkholderia cepacia complex (Bcc) cause life-threatening respiratory infections and are highly resistant to antibiotics. One promising alternative being explored to combat Bcc infections is phage therapy (PT): the use of phages to treat bacterial infections. Unfortunately, the utility of PT against many pathogenic species is limited by its prevailing paradigm: that only obligately lytic phages should be used therapeutically. It is thought that 'lysogenic' phages do not lyse all bacteria and can transfer antimicrobial resistance or virulence factors to their hosts. We argue that the tendency of a lysogenization-capable (LC) phage to form stable lysogens is not predicated exclusively on its ability to do so, and that the therapeutic suitability of a phage must be evaluated on a case-by-case basis. Concordantly, we developed several novel metrics-Efficiency of Phage Activity, Growth Reduction Coefficient, and Stable Lysogenization Frequency-and used them to evaluate eight Bcc-specific phages. Although these parameters vary considerably among Bcc phages, a strong inverse correlation (R2 = 0.67; P < 0.0001) exists between lysogen formation and antibacterial activity, indicating that certain LC phages with low frequency of stable lysogenization may be therapeutically efficacious. Moreover, we show that many LC Bcc phages interact synergistically with other phages in the first reported instance of mathematically defined polyphage synergy, and that these interactions result in the eradication of in vitro bacterial growth. Together, these findings reveal a novel therapeutic role for LC phages and challenge the current paradigm of PT. IMPORTANCE The spread of antimicrobial resistance is an imminent threat to public health around the world. Particularly concerning are species of the Burkholderia cepacia complex (Bcc), which cause life-threatening respiratory infections and are notoriously resistant to antibiotics. Phage therapy is a promising alternative being explored to combat Bcc infections and antimicrobial resistance in general, but its utility against many pathogenic species, including the Bcc, is restricted by the currently prevailing paradigm of exclusively using rare obligately lytic phages due to the perception that 'lysogenic' phages are therapeutically unsuitable. Our findings show that many lysogenization-capable phages exhibit powerful in vitro antibacterial activity both alone and through mathematically defined synergistic interactions with other phages, demonstrating a novel therapeutic role for LC phages and therefore challenging the currently prevailing paradigm of PT.

Diversity and potential function of pig gut DNA viruses

Viruses are ubiquitous in the gut of animals and play an important role in the ecology of the gut microbiome. The potential effects of these substances on the growth and development of the body are not fully known. Little is known about the effects of breeding environment on pig gut virome. Here, there are 3584 viral operational taxonomic units (vOTUs) longer than 5 kb identified by virus-enriched metagenome sequencing from 25 pig fecal samples. Only a small minority of vOTUs (11.16%) can be classified at the family level, and ∼50% of the genes could be annotated, supporting the concept of pig gut as reservoirs of substantial undescribed viral genetic diversity. The composition of pig gut virome in the six regions may be related to geography. There are only 20 viral clusters (VCs) shared among pig gut virome in six regions of Shanxi Province. These viruses rarely carry antibiotic resistance genes (ARGs). At the same time, they possess abundant auxiliary metabolic genes (AMGs) potentially involved in carbon, sulfur metabolism and cofactor biosynthesis, etc. This study has revealed the unique characteristics and potential function of pig gut DNA virome and established a foundation for the recognition of the viral roles in gut environment.

The Isolation and Characterization of a Broad Host Range Bcep22-like Podovirus JC1

Bacteriophage JC1 is a Podoviridae phage with a C1 morphotype, isolated on host strain Burkholderia cenocepacia Van1. Phage JC1 is capable of infecting an expansive range of Burkholderia cepacia complex (Bcc) species. The JC1 genome exhibits significant similarity and synteny to Bcep22-like phages and to many Ralstonia phages. The genome of JC1 was determined to be 61,182 bp in length with a 65.4% G + C content and is predicted to encode 76 proteins and 1 tRNA gene. Unlike the other Lessieviruses, JC1 encodes a putative helicase gene in its replication module, and it is in a unique organization not found in previously analyzed phages. The JC1 genome also harbours 3 interesting moron genes, that encode a carbon storage regulator (CsrA), an N-acetyltransferase, and a phosphoadenosine phosphosulfate (PAPS) reductase. JC1 can stably lysogenize its host Van1 and integrates into the 5' end of the gene rimO. This is the first account of stable integration identified for Bcep22-like phages. JC1 has a higher global virulence index at 37 °C than at 30 °C (0.8 and 0.21, respectively); however, infection efficiency and lysogen stability are not affected by a change in temperature, and no observable temperature-sensitive switch between lytic and lysogenic lifestyle appears to exist. Although JC1 can stably lysogenize its host, it possesses some desirable characteristics for use in phage therapy. Phage JC1 has a broad host range and requires the inner core of the bacterial LPS for infection. Bacteria that mutate to evade infection by JC1 may develop a fitness disadvantage as seen in previously characterized LPS mutants lacking inner core.

Evolutionary genomics of APSE: a tailed phage that lysogenically converts the bacterium Hamiltonella defensa into a heritable protective symbiont of aphids

Background

Most phages infect free-living bacteria but a few have been identified that infect heritable symbionts of insects or other eukaryotes. Heritable symbionts are usually specialized and isolated from other bacteria with little known about the origins of associated phages. Hamiltonella defensa is a heritable bacterial symbiont of aphids that is usually infected by a tailed, double-stranded DNA phage named APSE.

Methods

We conducted comparative genomic and phylogenetic studies to determine how APSE is related to other phages and prophages.

Results

Each APSE genome was organized into four modules and two predicted functional units. Gene content and order were near-fully conserved in modules 1 and 2, which encode predicted DNA metabolism genes, and module 4, which encodes predicted virion assembly genes. Gene content of module 3, which contains predicted toxin, holin and lysozyme genes differed among haplotypes. Comparisons to other sequenced phages suggested APSE genomes are mosaics with modules 1 and 2 sharing similarities with Bordetella-Bcep-Xylostella fastidiosa-like podoviruses, module 4 sharing similarities with P22-like podoviruses, and module 3 sharing no similarities with known phages. Comparisons to other sequenced bacterial genomes identified APSE-like elements in other heritable insect symbionts (Arsenophonus spp.) and enteric bacteria in the family Morganellaceae.

Conclusions

APSEs are most closely related to phage elements in the genus Arsenophonus and other bacteria in the Morganellaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12985-021-01685-y.

Advances in Phage Therapy: Targeting the Burkholderia cepacia Complex

The increasing prevalence and worldwide distribution of multidrug-resistant bacterial pathogens is an imminent danger to public health and threatens virtually all aspects of modern medicine. Particularly concerning, yet insufficiently addressed, are the members of the Burkholderia cepacia complex (Bcc), a group of at least twenty opportunistic, hospital-transmitted, and notoriously drug-resistant species, which infect and cause morbidity in patients who are immunocompromised and those afflicted with chronic illnesses, including cystic fibrosis (CF) and chronic granulomatous disease (CGD). One potential solution to the antimicrobial resistance crisis is phage therapy-the use of phages for the treatment of bacterial infections. Although phage therapy has a long and somewhat checkered history, an impressive volume of modern research has been amassed in the past decades to show that when applied through specific, scientifically supported treatment strategies, phage therapy is highly efficacious and is a promising avenue against drug-resistant and difficult-to-treat pathogens, such as the Bcc. In this review, we discuss the clinical significance of the Bcc, the advantages of phage therapy, and the theoretical and clinical advancements made in phage therapy in general over the past decades, and apply these concepts specifically to the nascent, but growing and rapidly developing, field of Bcc phage therapy.

High genomic diversity of novel phages infecting the plant pathogen Ralstonia solanacearum, isolated in Mauritius and Reunion islands

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) is among the most important plant diseases worldwide, severely affecting a high number of crops and ornamental plants in tropical regions. Only a limited number of phages infecting R. solanacearum have been isolated over the years, despite the importance of this bacterium and the associated plant disease. The antibacterial effect or morphological traits of these R. solanacearum viruses have been well studied, but not their genomic features, which need deeper consideration. This study reports the full genome of 23 new phages infecting RSSC isolated from agricultural samples collected in Mauritius and Reunion islands, particularly affected by this plant bacterial pathogen and considered biodiversity hotspots in the Southwest Indian Ocean. The complete genomic information and phylogenetic classification is provided, revealing high genetic diversity between them and weak similarities with previous related phages. The results support our proposal of 13 new species and seven new genera of R. solanacearum phages. Our findings highlight the wide prevalence of phages of RSSC in infected agricultural settings and the underlying genetic diversity. Discoveries of this kind lead more insight into the diversity of phages in general and to optimizing their use as biocontrol agents of bacterial diseases of plants in agriculture.

Infection Kinetics and Phylogenetic Analysis of vB_EcoD_SU57, a Virulent T1-Like Drexlerviridae Coliphage

The morphology, infection kinetics, genome sequence and phylogenetic characterization of the previously isolated bacteriophage vB_EcoD_SU57 are presented. The phage vB_EcoD_SU57 was isolated on Escherichia coli strain ECOR57 from the E. coli reference collection and was shown to produce four mm clear plaques with halos. Infection kinetics, as assessed by one-step growth analyses, suggest that vB_EcoD_SU57 is a virulent phage with an adsorption rate of 8.5 × 10-10 mL × min-1, a latency period of 14 min, and a burst size of 13 PFU per bacterium. Transmission electron microscopy confirmed vB_EcoD_SU57 to be a phage that used to be classified as a Siphoviridae phage. Bioinformatics analyses showed that the genome was 46,150 base pairs long, contained 29 genes with predicted protein functions, and 51 open reading frames encoding proteins with unknown function, many of which were gathered in clusters. A putative tRNA gene was also identified. Phylogenetic analyses showed that vB_EcoD_SU57 is a Braunvirinae phage of the newly formed Drexlerviridae family and closely related to T1-like E. coli phages vB_EcoS_ACG-M12 (Guelphvirus) and Rtp (Rtpvirus) as well as the unclassified phages vB_EcoS_CEB_EC3a and ECH1.

Comparative Genome Analysis of Uropathogenic Morganella morganii Strains

Morganella morganii is an opportunistic bacterial pathogen shown to cause a wide range of clinical and community-acquired infections. This study was aimed at sequencing and comparing the genomes of three M. morganii strains isolated from the urine samples of patients with community-acquired urinary tract infections. Draft genome sequencing was conducted using the Illumina HiSeq platform. The genomes of MM 1, MM 4, and MM 190 strains have a size of 3.82–3.97 Mb and a GC content of 50.9–51%. Protein-coding sequences (CDS) represent 96.1% of the genomes, RNAs are encoded by 2.7% of genes and pseudogenes account for 1.2% of the genomes. The pan-genome containes 4,038 CDS, of which 3,279 represent core genes. Six to ten prophages and 21–33 genomic islands were identified in the genomes of MM 1, MM 4, and MM 190. More than 30 genes encode capsular biosynthesis proteins, an average of 60 genes encode motility and chemotaxis proteins, and about 70 genes are associated with fimbrial biogenesis and adhesion. We determined that all strains contained urease gene cluster ureABCEFGD and had a urease activity. Both MM 4 and MM 190 strains are capable of hemolysis and their activity correlates well with a cytotoxicity level on T-24 bladder carcinoma cells. These activities were associated with expression of RTX toxin gene hlyA, which was introduced into the genomes by a phage similar to Salmonella phage 118970_sal4.

Complete Genome of the Xanthomonas euvesicatoria Specific Bacteriophage KΦ1, Its Survival and Potential in Control of Pepper Bacterial Spot

Xanthomonas euvesicatoria phage KΦ1, a member of Myoviridae family, was isolated from the rhizosphere of pepper plants showing symptoms of bacterial spot. The phage strain expressed antibacterial activity to all X. euvesicatoria strains tested and did not lyse other Xanthomonas spp., nor other less related bacterial species. The genome of KΦ1 is double-stranded DNA of 46.077 bp including 66 open reading frames and an average GC content of 62.9%, representing the first complete genome sequence published for a phage infecting xanthomonads associated with pepper or tomato. The highest genome similarity was observed between phage KΦ1 and the Xanthomonas oryzae pv. oryzae specific phage OP2. On the other hand, when compared with other members of the genus Bcep78virus, the genome similarity was lower. Forty-four (67%) predicted KΦ1 proteins shared homology with Xanthomonas phage OP2, while 20 genes (30%) were unique to KΦ1. Phage KΦ1, which is chloroform resistant and stable in different media and in the pH range 5-11, showed a high titer storage ability for at least 2 years at +4°C. Copper-hydroxide and copper-oxychloride reduced phage activity proportionally to the used concentrations and the exposure time. UV light was detrimental to the phage strain, but skim milk plus sucrose formulation extended its survival in vitro. The phages survived for at least 7 days on the surface of pepper leaves in the greenhouse, showing the ability to persist on the plant tissue without the presence of the host bacterium. Results of three repeated experiments showed that foliar applications of the unformulated KΦ1 phage suspension effectively controlled pepper bacterial spot compared to the standard treatment and the untreated control. The integration of the phage KΦ1 and copper-hydroxide treatments resulted in an increased efficacy compared to the copper-hydroxide alone.

Dietary energy drives the dynamic response of bovine rumen viral communities

BACKGROUND

Rumen microbes play a greater role in host energy acquisition than that of gut-associated microbes in monogastric animals. Although genome-enabled advancements are providing access to the vast diversity of uncultivated microbes, our understanding of variables shaping rumen microbial communities is in its infancy. Viruses have been shown to impact microbial populations through a myriad of processes, including cell lysis and reprogramming of host metabolism. However, little is known about the processes shaping the distribution of rumen viruses or how viruses may modulate microbial-driven processes in the rumen. To this end, we investigated how rumen bacterial and viral community structure and function responded in five steers fed four randomized dietary treatments in a crossover design.

RESULTS

Total digestible nutrients (TDN), a measure of dietary energy, best explained the variation in bacterial and viral communities. Additional ecological drivers of viral communities included dietary zinc content and microbial functional diversity. Using partial least squares regression, we demonstrate significant associations between the abundances of 267 viral populations and variables driving the variation in rumen viral communities. While rumen viruses were dynamic, 14 near ubiquitous viral populations were identified, suggesting the presence of a core rumen virome largely comprised of novel viruses. Moreover, analysis of virally encoded auxiliary metabolic genes (AMGs) indicates rumen viruses have glycosidic hydrolases to potentially augment the breakdown of complex carbohydrates to increase energy production. Other AMGs identified have a role in redirecting carbon to the pentose phosphate pathway and one carbon pools by folate to boost viral replication.

CONCLUSIONS

We demonstrate that rumen bacteria and viruses have differing responses and ecological drivers to dietary perturbation. Our results show that rumen viruses have implications for understanding the structuring of the previously identified core rumen microbiota and impacting microbial metabolism through a vast array of AMGs. AMGs in the rumen appear to have consequences for microbial metabolism that are largely in congruence with the current paradigm established in marine systems. This study provides a foundation for future hypotheses regarding the dynamics of viral-mediated processes in the rumen.

Impacts of Antibiotic and Bacteriophage Treatments on the Gut-Symbiont-Associated Blissus insularis (Hemiptera: Blissidae)

The Southern chinch bug, Blissus insularis, possesses specialized midgut crypts that harbor dense populations of the exocellular symbiont Burkholderia. Oral administration of antibiotics suppressed the gut symbionts in B. insularis and negatively impacted insect host fitness, as reflected by retarded development, smaller body size, and higher susceptibility to an insecticide, bifenthrin. Considering that the antibiotics probably had non-lethal but toxic effects on host fitness, attempts were conducted to reduce gut symbionts using bacteriophage treatment. Soil-lytic phages active against the cultures of specific Burkholderia ribotypes were successfully isolated using a soil enrichment protocol. Characterization of the BiBurk16MC_R phage determined its specificity to the Bi16MC_R_vitro ribotype and placed it within the family Podoviridae. Oral administration of phages to fifth-instar B. insularis, inoculated with Bi16MC_R_vitro as neonates had no deleterious effects on host fitness. However, the ingested phages failed to impact the crypt-associated Burkholderia. The observed inactivity of the phage was likely due to the blockage of the connection between the anterior and posterior midgut regions. These findings suggest that the initial colonization by Burkholderia programs the ontogeny of the midgut, providing a sheltered residence protected from microbial antagonists.

Bacteriophage application restores ethanol fermentation characteristics disrupted by Lactobacillus fermentum

BACKGROUND

Contamination of corn mash by lactic acid bacteria (LAB) reduces the efficiency of the ethanol fermentation process. The industry relies heavily on antibiotics for contamination control and there is a need to develop alternative methods. The goals of this study were to determine the diversity and abundance of bacteria contaminating commercial ethanol fermentations, and to evaluate the potential of anti-LAB bacteriophages in controlling production losses.

RESULTS

Bacterial populations in 27 corn mash samples collected from nine different commercial plants were determined by pyrosequencing of 16S rRNA amplicons. The results showed that the most abundant bacteria (>50 % of total population) in 24 of the 27 samples included LAB genera such as Lactobacillus, Streptococcus, Lactococcus, Weissella, Enterococcus, and Pediococcus. Lactobacillus was identified as the most prevalent genus at all fermentation stages in all plants, accounting for between 2.3 and 93.7 % of each population and constituting the major genus (>50 %) in nine samples from five plants and the most abundant genus in five other samples. Lactobacillus species, including L. delbrueckii, L. fermentum, L. mucosae, and L. reuteri were the most well-represented species. Two bacteriophages that target L. fermentum strains from ethanol plants, vB_LfeS_EcoSau and vB_LfeM_EcoInf (EcoSau and EcoInf), were isolated and characterized as a siphophage and a myophage, respectively. Analysis of the 31,703 bp genome of EcoSau revealed its similarity to the P335-like phage group, and the 106,701 bp genome of phage EcoInf was determined to be a novel phage type despite its distant relationship to the SPO1-like phages. Addition of phages EcoSau and EcoInf to L. fermentum-contaminated corn mash fermentation models restored the yields of ethanol and reduced levels of residual glucose, lactic acid, and acetic acid to that comparable to the infection-free control.

CONCLUSIONS

This study provides detailed insight into the microbiota contaminating commercial ethanol fermentations, and highlights the abundance of LAB, especially L. delbrueckii, L. fermentum, L. mucosae, and L. reuteri, in the process. This study suggests that phages with broad coverage of major LAB species can be applied directly to corn mash for antibiotic-free control of contamination in the ethanol fermentation industry.

Isolation and Characterization of an Aeromonas punctata Bacteriophage

An Aeromonas punctata bacteriophage, named as DH1, was isolated from East Lake, Wuhan city, China. Morphologically, phage DH1 showed a typical Myoviridae structure consisting of an isometric head (50 nm in diameter) and a visible tail. The bacteriophage had a latent period of about 90 minutes and an average burst size of about 125 PFU•Cell^-1. Restriction enzyme pattern of the bacteriophage’s genome showed that the genome is a double-stranded DNA and about 34kb in size. The sequenced genomic fragments showed highly similarities to gp04 and gp16 sequence of other Myoviridae bacteriophages at protein level.

Salmonella phages and prophages: genomics, taxonomy, and applied aspects

Since this book was originally published in 2007 there has been a significant increase in the number of Salmonella bacteriophages, particularly lytic virus, and Salmonella strains which have been fully sequenced. In addition, new insights into phage taxonomy have resulted in new phage genera, some of which have been recognized by the International Committee of Taxonomy of Viruses (ICTV). The properties of each of these genera are discussed, along with the role of phage as agents of genetic exchange, as therapeutic agents, and their involvement in phage typing.

Integration of genomic and proteomic analyses in the classification of the Siphoviridae family

Using a variety of genomic (BLASTN, ClustalW) and proteomic (Phage Proteomic Tree, CoreGenes) tools we have tackled the taxonomic status of members of the largest bacteriophage family, the Siphoviridae. In all over 400 phages were examined and we were able to propose 39 new genera, comprising 216 phage species, and add 62 species to two previously defined genera (Phic3unalikevirus; L5likevirus) grouping, in total, 390 fully sequenced phage isolates. Many of the remainders are orphans which the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) chooses not to ascribe genus status at the time being.

Genomic characterization of JG068, a novel virulent podovirus active against Burkholderia cenocepacia

BACKGROUND

As is true for many other antibiotic-resistant Gram-negative pathogens, members of the Burkholderia cepacia complex (BCC) are currently being assessed for their susceptibility to phage therapy as an antimicrobial treatment. The objective of this study was to perform genomic and limited functional characterization of the novel BCC phage JG068 (vB_BceP_JG068).

RESULTS

JG068 is a podovirus that forms large, clear plaques on Burkholderia cenocepacia K56-2. Host range analysis indicates that this phage can infect environmental, clinical, and epidemic isolates of Burkholderia multivorans, B. cenocepacia, Burkholderia stabilis, and Burkholderia dolosa, likely through interaction with the host lipopolysaccharide as a receptor. The JG068 chromosome is 41,604 base pairs (bp) in length and is flanked by 216 bp short direct terminal repeats. Gene expression originates from both host and phage promoters and is in the forward direction for all 49 open reading frames. The genome sequence shows similarity to Ralstonia phage ϕRSB1, Caulobacter phage Cd1, and uncharacterized genetic loci of blood disease bacterium R229 and Burkholderia pseudomallei 1710b. CoreGenesUniqueGenes analysis indicates that JG068 belongs to the Autographivirinae subfamily and ϕKMV-like phages genus. Modules within the genome encode proteins involved in DNA-binding, morphogenesis, and lysis, but none associated with pathogenicity or lysogeny. Similar to the signal-arrest-release (SAR) endolysin of ϕKMV, inducible expression of the JG068 SAR endolysin causes lysis of Escherichia coli that is dependent on the presence of an N-terminal signal sequence. In an in vivo assay using the Galleria mellonella infection model, treatment of B. cenocepacia K56-2-infected larvae with JG068 results in a significant increase in larval survival.

CONCLUSIONS

As JG068 has a broad host range, does not encode virulence factors, is obligately lytic, and has activity against an epidemic B. cenocepacia strain in vivo, this phage is a highly promising candidate for BCC phage therapy development.

The Caulobacter crescentus phage phiCbK: genomics of a canonical phage

Background

The bacterium Caulobacter crescentus is a popular model for the study of cell cycle regulation and senescence. The large prolate siphophage phiCbK has been an important tool in C. crescentus biology, and has been studied in its own right as a model for viral morphogenesis. Although a system of some interest, to date little genomic information is available on phiCbK or its relatives.

Results

Five novel phiCbK-like C. crescentus bacteriophages, CcrMagneto, CcrSwift, CcrKarma, CcrRogue and CcrColossus, were isolated from the environment. The genomes of phage phiCbK and these five environmental phage isolates were obtained by 454 pyrosequencing. The phiCbK-like phage genomes range in size from 205 kb encoding 318 proteins (phiCbK) to 280 kb encoding 448 proteins (CcrColossus), and were found to contain nonpermuted terminal redundancies of 10 to 17 kb. A novel method of terminal ligation was developed to map genomic termini, which confirmed termini predicted by coverage analysis. This suggests that sequence coverage discontinuities may be useable as predictors of genomic termini in phage genomes. Genomic modules encoding virion morphogenesis, lysis and DNA replication proteins were identified. The phiCbK-like phages were also found to encode a number of intriguing proteins; all contain a clearly T7-like DNA polymerase, and five of the six encode a possible homolog of the C. crescentus cell cycle regulator GcrA, which may allow the phage to alter the host cell’s replicative state. The structural proteome of phage phiCbK was determined, identifying the portal, major and minor capsid proteins, the tail tape measure and possible tail fiber proteins. All six phage genomes are clearly related; phiCbK, CcrMagneto, CcrSwift, CcrKarma and CcrRogue form a group related at the DNA level, while CcrColossus is more diverged but retains significant similarity at the protein level.

Conclusions

Due to their lack of any apparent relationship to other described phages, this group is proposed as the founding cohort of a new phage type, the phiCbK-like phages. This work will serve as a foundation for future studies on morphogenesis, infection and phage-host interactions in C. crescentus.

Artificial neural networks trained to detect viral and phage structural proteins

Phages play critical roles in the survival and pathogenicity of their hosts, via lysogenic conversion factors, and in nutrient redistribution, via cell lysis. Analyses of phage- and viral-encoded genes in environmental samples provide insights into the physiological impact of viruses on microbial communities and human health. However, phage ORFs are extremely diverse of which over 70% of them are dissimilar to any genes with annotated functions in GenBank. Better identification of viruses would also aid in better detection and diagnosis of disease, in vaccine development, and generally in better understanding the physiological potential of any environment. In contrast to enzymes, viral structural protein function can be much more challenging to detect from sequence data because of low sequence conservation, few known conserved catalytic sites or sequence domains, and relatively limited experimental data. We have designed a method of predicting phage structural protein sequences that uses Artificial Neural Networks (ANNs). First, we trained ANNs to classify viral structural proteins using amino acid frequency; these correctly classify a large fraction of test cases with a high degree of specificity and sensitivity. Subsequently, we added estimates of protein isoelectric points as a feature to ANNs that classify specialized families of proteins, namely major capsid and tail proteins. As expected, these more specialized ANNs are more accurate than the structural ANNs. To experimentally validate the ANN predictions, several ORFs with no significant similarities to known sequences that are ANN-predicted structural proteins were examined by transmission electron microscopy. Some of these self-assembled into structures strongly resembling virion structures. Thus, our ANNs are new tools for identifying phage and potential prophage structural proteins that are difficult or impossible to detect by other bioinformatic analysis. The networks will be valuable when sequence is available but in vitro propagation of the phage may not be practical or possible.

Bacteriophages are extremely abundant and diverse biological entities. All phage particles are comprised of nucleic acids and structural proteins, with few other packaged proteins. Despite their simplicity and abundance, more than 70% of phage sequences in the viral Reference Sequence database encode proteins with unknown function based on FASTA annotations. As a result, the use of sequence similarity is often insufficient for detecting virus structural proteins among unknown viral sequences. Viral structural protein function is challenging to detect from sequence data because structural proteins possess few known conserved catalytic motifs and sequence domains. To address these issues we investigated the use of Artificial Neural Networks as an alternative means of predicting function. Here, we trained thousands of networks using the amino acid frequency of structural protein sequences and identified the optimal architectures with the highest accuracies. Some hypothetical protein sequences detected by our networks were expressed and visualized by TEM, and produced images that strongly resemble virion structures. Our results support the utility of our neural networks in predicting the functions of unknown viral sequences.

Comparative analysis of two phenotypically-similar but genomically-distinct Burkholderia cenocepacia-specific bacteriophages

Background

Genomic analysis of bacteriophages infecting the Burkholderia cepacia complex (BCC) is an important preliminary step in the development of a phage therapy protocol for these opportunistic pathogens. The objective of this study was to characterize KL1 (vB_BceS_KL1) and AH2 (vB_BceS_AH2), two novel Burkholderia cenocepacia-specific siphoviruses isolated from environmental samples.

Results

KL1 and AH2 exhibit several unique phenotypic similarities: they infect the same B. cenocepacia strains, they require prolonged incubation at 30°C for the formation of plaques at low titres, and they do not form plaques at similar titres following incubation at 37°C. However, despite these similarities, we have determined using whole-genome pyrosequencing that these phages show minimal relatedness to one another. The KL1 genome is 42,832 base pairs (bp) in length and is most closely related to Pseudomonas phage 73 (PA73). In contrast, the AH2 genome is 58,065 bp in length and is most closely related to Burkholderia phage BcepNazgul. Using both BLASTP and HHpred analysis, we have identified and analyzed the putative virion morphogenesis, lysis, DNA binding, and MazG proteins of these two phages. Notably, MazG homologs identified in cyanophages have been predicted to facilitate infection of stationary phase cells and may contribute to the unique plaque phenotype of KL1 and AH2.

Conclusions

The nearly indistinguishable phenotypes but distinct genomes of KL1 and AH2 provide further evidence of both vast diversity and convergent evolution in the BCC-specific phage population.

Genome and proteome analysis of 7-7-1, a flagellotropic phage infecting Agrobacterium sp H13-3

Background

The flagellotropic phage 7-7-1 infects motile cells of Agrobacterium sp H13-3 by attaching to and traveling along the rotating flagellar filament to the secondary receptor at the base, where it injects its DNA into the host cell. Here we describe the complete genomic sequence of 69,391 base pairs of this unusual bacteriophage.

Methods

The sequence of the 7-7-1 genome was determined by pyro(454)sequencing to a coverage of 378-fold. It was annotated using MyRAST and a variety of internet resources. The structural proteome was analyzed by SDS-PAGE coupled electrospray ionization-tandem mass spectrometry (MS/MS).

Results

Sequence annotation and a structural proteome analysis revealed 127 open reading frames, 84 of which are unique. In six cases 7-7-1 proteins showed sequence similarity to proteins from the virulent Burkholderia myovirus BcepB1A. Unique features of the 7-7-1 genome are the physical separation of the genes encoding the small (orf100) and large (orf112) subunits of the DNA packaging complex and the apparent lack of a holin-lysin cassette. Proteomic analysis revealed the presence of 24 structural proteins, five of which were identified as baseplate (orf7), putative tail fibre (orf102), portal (orf113), major capsid (orf115) and tail sheath (orf126) proteins. In the latter case, the N-terminus was removed during capsid maturation, probably by a putative prohead protease (orf114).

Characterization of modular bacteriophage endolysins from Myoviridae phages OBP, 201φ2-1 and PVP-SE1

Peptidoglycan lytic enzymes (endolysins) induce bacterial host cell lysis in the late phase of the lytic bacteriophage replication cycle. Endolysins OBPgp279 (from Pseudomonas fluorescens phage OBP), PVP-SE1gp146 (Salmonella enterica serovar Enteritidis phage PVP-SE1) and 201φ2-1gp229 (Pseudomonas chlororaphis phage 201φ2-1) all possess a modular structure with an N-terminal cell wall binding domain and a C-terminal catalytic domain, a unique property for endolysins with a Gram-negative background. All three modular endolysins showed strong muralytic activity on the peptidoglycan of a broad range of Gram-negative bacteria, partly due to the presence of the cell wall binding domain. In the case of PVP-SE1gp146, this domain shows a binding affinity for Salmonella peptidoglycan that falls within the range of typical cell adhesion molecules (K(aff) = 1.26 × 10(6) M(-1)). Remarkably, PVP-SE1gp146 turns out to be thermoresistant up to temperatures of 90 °C, making it a potential candidate as antibacterial component in hurdle technology for food preservation. OBPgp279, on the other hand, is suggested to intrinsically destabilize the outer membrane of Pseudomonas species, thereby gaining access to their peptidoglycan and exerts an antibacterial activity of 1 logarithmic unit reduction. Addition of 0.5 mM EDTA significantly increases the antibacterial activity of the three modular endolysins up to 2-3 logarithmic units reduction. This research work offers perspectives towards elucidation of the structural differences explaining the unique biochemical and antibacterial properties of OBPgp279, PVP-SE1gp146 and 201φ2-1gp229. Furthermore, these endolysins extensively enlarge the pool of potential antibacterial compounds used against multi-drug resistant Gram-negative bacterial infections.

Cangene Gold Medal Award Lecture — Genomic analysis and modification ofBurkholderia cepaciacomplex bacteriophages1This article is based on a presentation by Dr. Karlene Lynch at the 61st Annual Meeting of the Canadian Society of Microbiologists in St. John’s, Newfoundland and Labrador, on 21 June 2011. Dr. Lynch was the recipient of the 2011 Cangene Gold Medal as the Canadian Graduate Student Microbiologist of the Year, an annual award sponsored by Cangene Corporation intended to recognize excellence in graduate research

The Burkholderia cepacia complex (Bcc) is a group of 17 Gram-negative predominantly environmental bacterial species that cause potentially fatal opportunistic infections in cystic fibrosis (CF) patients. Although its prevalence in these individuals is lower than that of Staphylococcus aureus and Pseudomonas aeruginosa , the Bcc remains a serious problem in the CF community because of the pathogenicity, transmissibility, and inherent antibiotic resistance of these organisms. An alternative treatment for Bcc infections that is currently being developed is phage therapy, the clinical use of viruses that infect bacteria. To assess the suitability of individual phage isolates for therapeutic use, the complete genome sequences of a panel of Bcc‐specific phages were determined and analyzed. These sequences encode a broad range of proteins with a gradient of relatedness to phage and bacterial gene products from Burkholderia and other genera. The majority of these phages were found not to encode virulence factors, and despite their predominantly temperate nature, a proof-of-principle experiment has shown that they may be modified to a lytic form. Both the genomic characterization and subsequent engineering of Bcc‐specific phages are fundamental to the development of an effective phage therapy strategy for these bacteria.

The promise of bacteriophage therapy for Burkholderia cepacia complex respiratory infections

In recent times, increased attention has been given to evaluating the efficacy of phage therapy, especially in scenarios where the bacterial infectious agent of interest is highly antibiotic resistant. In this regard, phage therapy is especially applicable to infections caused by the Burkholderia cepacia complex (BCC) since members of the BCC are antibiotic pan-resistant. Current studies in BCC phage therapy are unique from many other avenues of phage therapy research in that the investigation is not only comprised of phage isolation, in vitro phage characterization and assessment of in vivo infection model efficacy, but also adapting aerosol drug delivery techniques to aerosol phage formulation delivery and storage.

Genomes and characterization of phages Bcep22 and BcepIL02, founders of a novel phage type in Burkholderia cenocepacia

Within the Burkholderia cepacia complex, B. cenocepacia is the most common species associated with aggressive infections in the lungs of cystic fibrosis patients, causing disease that is often refractive to treatment by antibiotics. Phage therapy may be a potential alternative form of treatment for these infections. Here we describe the genome of the previously described therapeutic B. cenocepacia podophage BcepIL02 and its close relative, Bcep22. Phage Bcep22 was found to contain a circularly permuted genome of 63,882 bp containing 77 genes; BcepIL02 was found to be 62,714 bp and contains 76 predicted genes. Major virion-associated proteins were identified by proteomic analysis. We propose that these phages comprise the founding members of a novel podophage lineage, the Bcep22-like phages. Among the interesting features of these phages are a series of tandemly repeated putative tail fiber genes that are similar to each other and also to one or more such genes in the other phages. Both phages also contain an extremely large (ca. 4,600-amino-acid), virion-associated, multidomain protein that accounts for over 20% of the phages' coding capacity, is widely distributed among other bacterial and phage genomes, and may be involved in facilitating DNA entry in both phage and other mobile DNA elements. The phages, which were previously presumed to be virulent, show evidence of a temperate lifestyle but are apparently unable to form stable lysogens in their hosts. This ambiguity complicates determination of a phage lifestyle, a key consideration in the selection of therapeutic phages.

Genomic and functional analyses of Rhodococcus equi phages ReqiPepy6, ReqiPoco6, ReqiPine5, and ReqiDocB7

The isolation and results of genomic and functional analyses of Rhodococcus equi phages ReqiPepy6, ReqiDocB7, ReqiPine5, and ReqiPoco6 (hereafter referred to as Pepy6, DocB7, Pine5, and Poco6, respectively) are reported. Two phages, Pepy6 and Poco6, more than 75% identical, exhibited genome organization and protein sequence likeness to Lactococcus lactis phage 1706 and clostridial prophage elements. An unusually high fraction, 27%, of Pepy6 and Poco6 proteins were predicted to possess at least one transmembrane domain, a value much higher than the average of 8.5% transmembrane domain-containing proteins determined from a data set of 36,324 phage protein entries. Genome organization and protein sequence comparisons place phage Pine5 as the first nonmycobacteriophage member of the large Rosebush cluster. DocB7, which had the broadest host range among the four isolates, was not closely related to any phage or prophage in the database, and only 23 of 105 predicted encoded proteins could be assigned a functional annotation. Because of the relationship of Rhodococcus to Mycobacterium, it was anticipated that these phages should exhibit some of the features characteristic of mycobacteriophages. Traits that were identified as shared by the Rhodococcus phages and mycobacteriophages include the prevalent long-tailed morphology and the presence of genes encoding LysB-like mycolate-hydrolyzing lysis proteins. Application of DocB7 lysates to soils amended with a host strain of R. equi reduced recoverable bacterial CFU, suggesting that phage may be useful in limiting R. equi load in the environment while foals are susceptible to infection.

Identifying active phage lysins through functional viral metagenomics

Recent metagenomic sequencing studies of uncultured viral populations have provided novel insights into the ecology of environmental bacteriophage. At the same time, viral metagenomes could also represent a potential source of recombinant proteins with biotechnological value. In order to identify such proteins, a novel two-step screening technique was devised for cloning phage lytic enzymes from uncultured viral DNA. This plasmid-based approach first involves a primary screen in which transformed Escherichia coli clones that demonstrate colony lysis following exposure to inducing agent are identified. This effect, which can be due to the expression of membrane-permeabilizing phage holins, is discerned by the development a hemolytic effect in surrounding blood agar. In a secondary step, the clones identified in the primary screen are overlaid with autoclaved Gram-negative bacteria (specifically Pseudomonas aeruginosa) to assay directly for recombinant expression of lytic enzymes, which are often encoded proximally to holins in phage genomes. As proof-of-principle, the method was applied to a viral metagenomic library constructed from mixed animal feces, and 26 actively expressed lytic enzymes were cloned. These proteins include both Gram-positive-like and Gram-negative-like enzymes, as well as several atypical lysins whose predicted structures are less common among known phage. Overall, this study represents one of the first functional screens of a viral metagenomic population, and it provides a general approach for characterizing lysins from uncultured phage.

Genetic and phenotypic diversity in Burkholderia: contributions by prophage and phage-like elements

BACKGROUND

Burkholderia species exhibit enormous phenotypic diversity, ranging from the nonpathogenic, soil- and water-inhabiting Burkholderia thailandensis to the virulent, host-adapted mammalian pathogen B. mallei. Genomic diversity is evident within Burkholderia species as well. Individual isolates of Burkholderia pseudomallei and B. thailandensis, for example, carry a variety of strain-specific genomic islands (GIs), including putative pathogenicity and metabolic islands, prophage-like islands, and prophages. These GIs may provide some strains with a competitive advantage in the environment and/or in the host relative to other strains.

RESULTS

Here we present the results of analysis of 37 prophages, putative prophages, and prophage-like elements from six different Burkholderia species. Five of these were spontaneously induced to form bacteriophage particles from B. pseudomallei and B. thailandensis strains and were isolated and fully sequenced; 24 were computationally predicted in sequenced Burkholderia genomes; and eight are previously characterized prophages or prophage-like elements. The results reveal numerous differences in both genome structure and gene content among elements derived from different species as well as from strains within species, due in part to the incorporation of additional DNA, or 'morons' into the prophage genomes. Implications for pathogenicity are also discussed. Lastly, RNAseq analysis of gene expression showed that many of the genes in varphi1026b that appear to contribute to phage and lysogen fitness were expressed independently of the phage structural and replication genes.

CONCLUSIONS

This study provides the first estimate of the relative contribution of prophages to the vast phenotypic diversity found among the Burkholderiae.

Efficacy of bacteriophage therapy in a model of Burkholderia cenocepacia pulmonary infection

The therapeutic potential of bacteriophages (phages) in a mouse model of acute Burkholderia cenocepacia pulmonary infection was assessed. Phage treatment was administered by either intranasal inhalation or intraperitoneal injection. Bacterial density, macrophage inflammatory protein 2 (MIP-2), and tumor necrosis factor alpha (TNF-alpha) levels were significantly reduced in lungs of mice treated with intraperitoneal phages (P < .05). No significant differences in lung bacterial density or MIP-2 levels were found between untreated mice and mice treated with intranasal phages, intraperitoneal ultraviolet-inactivated phages, or intraperitoneal lambda phage control mice. Mock-infected mice treated with phage showed no significant increase in lung MIP-2 or TNF-alpha levels compared with mock-infected/mock-treated mice. We have demonstrated the efficacy of phage therapy in an acute B. cenocepacia lung infection model. Systemic phage administration was more effective than inhalational administration, suggesting that circulating phages have better access to bacteria in lungs than do topical phages.

Genomic and biological analysis of phage Xfas53 and related prophages of Xylella fastidiosa

We report the plaque propagation and genomic analysis of Xfas53, a temperate phage of Xylella fastidiosa. Xfas53 was isolated from supernatants of X. fastidiosa strain 53 and forms plaques on the sequenced strain Temecula. Xfas53 forms short-tailed virions, morphologically similar to podophage P22. The 36.7-kb genome is predicted to encode 45 proteins. The Xfas53 terminase and structural genes are related at a protein and gene order level to P22. The left arm of the Xfas53 genome has over 90% nucleotide identity to multiple prophage elements of the sequenced X. fastidiosa strains. This arm encodes proteins involved in DNA metabolism, integration, and lysogenic control. In contrast to Xfas53, each of these prophages encodes head and DNA packaging proteins related to the siphophage lambda and tail morphogenesis proteins related to those of myophage P2. Therefore, it appears that Xfas53 was formed by recombination between a widespread family of X. fastidiosa P2-related prophage elements and a podophage distantly related to phage P22. The lysis cassette of Xfas53 is predicted to encode a pinholin, a signal anchor and release (SAR) endolysin, and Rz and Rz1 equivalents. The holin gene encodes a pinholin and appears to be subject to an unprecedented degree of negative regulation at both the level of expression, with rho-independent transcriptional termination and RNA structure-dependent translational repression, and the level of holin function, with two upstream translational starts predicted to encode antiholin products. A notable feature of Xfas53 and related prophages is the presence of 220- to 390-nucleotide degenerate tandem direct repeats encoding putative DNA binding proteins. Additionally, each phage encodes at least two BroN domain-containing proteins possibly involved in lysogenic control. Xfas53 exhibits unusually slow adsorption kinetics, possibly an adaptation to the confined niche of its slow-growing host.

Isolation and genomic characterization of the first phage infecting Iodobacteria: ϕPLPE, a myovirus having a novel set of features

The aquatic phage ϕPLPE infects a bacterium of the genus Iodobacter that are common inhabitants of rivers, streams and canals that produce violacein-like pigments. Our characterization of ϕPLPE reveals it to be a small, contractile-tailed phage whose 47.5 kb genome sequence is phylogenetically distant from all previously characterized phages. The genome has a generally modular organization (e.g. replication/recombination, structure/morphogenesis, lysis/lysogeny) and approximately half of its 84 open reading frames have no known homologues. It behaves as a virulent phage under the host growth conditions we have employed and, with the exception of an anti-repressor (ant) homologue, the genome lacks all the genes associated with a lysogenic lifestyle. Thus, either ϕPLPE was once a temperate phage that has lost most of its lysogeny cassette or it is a virulent phage that acquired an ant-like gene presumably for some function other than the control of lysogeny. The ϕPLPE genome has few bacterial gene homologues with the interesting exception of a putative acylhydrolase (acylase). This function has been implicated in bacterial quorum sensing since it degrades homoserine-lactone signalling molecules and can disrupt or modulate quorum signalling from either the emitter or its competitors. ϕPLPE may be an example of a phage co-opting components of the bacterial quorum-sensing apparatus to its own advantage.

Classification of Myoviridae bacteriophages using protein sequence similarity

Background

We advocate unifying classical and genomic classification of bacteriophages by integration of proteomic data and physicochemical parameters. Our previous application of this approach to the entirely sequenced members of the Podoviridae fully supported the current phage classification of the International Committee on Taxonomy of Viruses (ICTV). It appears that horizontal gene transfer generally does not totally obliterate evolutionary relationships between phages.

Results

CoreGenes/CoreExtractor proteome comparison techniques applied to 102 Myoviridae suggest the establishment of three subfamilies (Peduovirinae, Teequatrovirinae, the Spounavirinae) and eight new independent genera (Bcep781, BcepMu, FelixO1, HAP1, Bzx1, PB1, phiCD119, and phiKZ-like viruses). The Peduovirinae subfamily, derived from the P2-related phages, is composed of two distinct genera: the "P2-like viruses", and the "HP1-like viruses". At present, the more complex Teequatrovirinae subfamily has two genera, the "T4-like" and "KVP40-like viruses". In the genus "T4-like viruses" proper, four groups sharing >70% proteins are distinguished: T4-type, 44RR-type, RB43-type, and RB49-type viruses. The Spounavirinae contain the "SPO1-"and "Twort-like viruses."

Conclusion

The hierarchical clustering of these groupings provide biologically significant subdivisions, which are consistent with our previous analysis of the Podoviridae.

Experimental bacteriophage therapy increases survival of Galleria mellonella larvae infected with clinically relevant strains of the Burkholderia cepacia complex

The Burkholderia cepacia complex (BCC) is a group of bacterial pathogens that are highly antibiotic resistant and associated with debilitating respiratory infections. Although bacteriophages of the BCC have been isolated and characterized, no studies have yet examined phage therapy against the BCC in vivo. In a caterpillar infection model, we show that BCC phage therapy is an alternative treatment possibility and is highly effective under specific conditions.

Preparation of a phage DNA fragment library for whole genome shotgun sequencing

The most efficient method to determine the genomic sequence of a dsDNA phage is to use a whole genome shotgun approach (WGSA). Preparation of a library where each genomic fragment has an equal chance of being represented is critical to the success of the WGSA. For many phages, there are regions of the genome likely to be under-represented in the shotgun library, which results in more gaps in the shotgun assembly than predicted by the Poisson distribution. However, as phage genomes are relatively small, this increased number of gaps does not present an insurmountable impediment to using the WGSA. This chapter will focus on construction of a high-quality random library and sequence analysis of this library in a 96-well format. Techniques are described for the mechanical fragmentation of genomic DNA into 2 kb average size fragments, preparation of the fragmented DNA for shotgun cloning, and advice on the choice of cloning vector for library preparation. Protocols for deepwell block culture, plasmid isolation, and sequencing in 96-well format are given. The rationale for determining the total number of random clones from a library to sequence for a 50 and 150 kb genome is explained. The steps involved in going from hundreds of shotgun sequencing traces to generating contigs will be outlined as well as how to close gaps in the sequence by primer walking on phage DNA and PCR-generated templates. Finally, examples will be given of how biological information about the phage genomic termini can be derived by analysis of the organization of individual clones in the shotgun sequence assembly. Specific examples are given for the circularly permuted termini of pac type phages, the direct terminal repeats found in most T7-like phages, variable host DNA at either end as in the Mu-like phages, and the 5' and 3' overhanging ends of cos type phages. The end result of these steps is the entire DNA sequence of a novel phage, ready for gene prediction.

Isolation of new Stenotrophomonas bacteriophages and genomic characterization of temperate phage S1

Twenty-two phages that infect Stenotrophomonas species were isolated through sewage enrichment and prophage induction. Of them, S1, S3, and S4 were selected due to their wide host ranges compared to those of the other phages. S1 and S4 are temperate siphoviruses, while S3 is a virulent myovirus. The genomes of S3 and S4, about 33 and 200 kb, were resistant to restriction digestion. The lytic cycles lasted 30 min for S3 and about 75 min for S1 and S4. The burst size for S3 was 100 virions/cell, while S1 and S4 produced about 75 virus particles/cell. The frequency of bacteriophage-insensitive host mutants, calculated by dividing the number of surviving colonies by the bacterial titer of a parallel, uninfected culture, ranged between 10(-5) and 10(-6) for S3 and 10(-3) and 10(-4) for S1 and S4. The 40,287-bp genome of S1 contains 48 open reading frames (ORFs) and 12-bp 5' protruding cohesive ends. By using a combination of bioinformatics and experimental evidence, functions were ascribed to 21 ORFs. The morphogenetic and lysis modules are well-conserved, but no lysis-lysogeny switch or DNA replication gene clusters were recognized. Two major clusters of genes with respect to transcriptional orientation were observed. Interspersed among them were lysogenic conversion genes encoding phosphoadenosine phosphosulfate reductase and GspM, a protein involved in the general secretion system II. The attP site of S1 may be located within a gene that presents over 75% homology to a Stenotrophomonas chromosomal determinant.

Unifying classical and molecular taxonomic classification: analysis of the Podoviridae using BLASTP-based tools

We defined phage genera by measuring genome relationships by the numbers of shared homologous/orthologous proteins. We used BLAST-based tools (CoreExtractor.vbs and CoreGenes) to analyze 55 fully sequenced bacteriophage genomes from the NCBI and EBI databases. This approach was first applied to the T7-related phages. Using a cut-off score of 40% homologous proteins, we identified three genera within the T7-related phages, redefined the phi29-related phages, and introduced five novel genera. The T7- and phi29-related phages were given subfamily status and named "Autographivirinae" and "Picovirinae", respectively. Our results confirm and refine the ICTV phage classification, enable elimination of errors in public databases, and provide a straightforward tool for the molecular classification of new phage genomes.

The Widespread Evolutionary Significance of Viruses

In the last 30 years, the study of virus evolution has undergone a transformation. Originally concerned with disease and its emergence, virus evolution had not been well integrated into the general study of evolution. This chapter reviews the developments that have brought us to this new appreciation for the general significance of virus evolution to all life. We now know that viruses numerically dominate all habitats of life, especially the oceans. Theoretical developments in the 1970s regarding quasispecies, error rates, and error thresholds have yielded many practical insights into virus–host dynamics. The human diseases of HIV-1 and hepatitis C virus cannot be understood without this evolutionary framework. Yet recent developments with poliovirus demonstrate that viral fitness can be the result of a consortia, not one fittest type, a basic Darwinian concept in evolutionary biology. Darwinian principles do apply to viruses, such as with Fisher population genetics, but other features, such as reticulated and quasispecies-based evolution distinguish virus evolution from classical studies. The available phylogenetic tools have greatly aided our analysis of virus evolution, but these methods struggle to characterize the role of virus populations. Missing from many of these considerations has been the major role played by persisting viruses in stable virus evolution and disease emergence. In many cases, extreme stability is seen with persisting RNA viruses. Indeed, examples are known in which it is the persistently infected host that has better survival. We have also recently come to appreciate the vast diversity of phage (DNA viruses) of prokaryotes as a system that evolves by genetic exchanges across vast populations (Chapter 10). This has been proposed to be the “big bang” of biological evolution. In the large DNA viruses of aquatic microbes we see surprisingly large, complex and diverse viruses. With both prokaryotic and eukaryotic DNA viruses, recombination is the main engine of virus evolution, and virus host co-evolution is common, although not uniform. Viral emergence appears to be an unending phenomenon and we can currently witness a selective sweep by retroviruses that infect and become endogenized in koala bears.

The genome of epsilon15, a serotype-converting, Group E1 Salmonella enterica-specific bacteriophage

The genome sequence of the Salmonella enterica serovar Anatum-specific, serotype-converting bacteriophage epsilon15 has been completed. The nonredundant genome contains 39,671 bp and 51 putative genes. It most closely resembles the genome of phiV10, an Escherichia coli O157:H7-specific temperate phage, with which it shares 36 related genes. More distant relatives include the Burkholderia cepacia-specific phage, BcepC6B (8 similar genes), the Bordetella bronchiseptica-specific phage, BPP-1 (8 similar genes) and the Photobacterium profundum prophage, P Pphipr1 (6 similar genes). epsilon15 gene identifications based on homologies with known gene families include the terminase small and large subunits, integrase, endolysin, two holins, two DNA methylase enzymes (one adenine-specific and one cytosine-specific) and a RecT-like enzyme. Genes identified experimentally include those coding for the serotype conversion proteins, the tail fiber, the major capsid protein and the major repressor. epsilon15's attP site and the Salmonella attB site with which it interacts during lysogenization have also been determined.

Rz/Rz1 lysis gene equivalents in phages of Gram-negative hosts

Under usual laboratory conditions, lysis by bacteriophage lambda requires only the holin and endolysin genes, but not the Rz and Rz1 genes, of the lysis cassette. Defects in Rz or Rz1 block lysis only in the presence of high concentrations of divalent cations. The lambda Rz and Rz1 lysis genes are remarkable in that Rz1, encoding an outer membrane lipoprotein, is completely embedded in the +1 register within Rz, which itself encodes an integral inner membrane protein. While Rz and Rz1 equivalents have been identified in T7 and P2, most phages, including such well-studied classic phages as T4, P1, T1, Mu and SP6, lack annotated Rz/Rz1 equivalents. Here we report that a search strategy based primarily on gene arrangement and membrane localization signals rather than sequence similarity has revealed that Rz/Rz1 equivalents are nearly ubiquitous among phages of Gram-negative hosts, with 120 of 137 phages possessing genes that fit the search criteria. In the case of T4, a deletion of a non-overlapping gene pair pseT.2 and pseT.3 identified as Rz/Rz1 equivalents resulted in the same divalent cation-dependent lysis phenotype. Remarkably, in T1 and six other phages, Rz/Rz1 pairs were not found but a single gene encoding an outer membrane lipoprotein with a C-terminal transmembrane domain capable of integration into the inner membrane was identified. These proteins were named "spanins," since their protein products are predicted to span the periplasm providing a physical connection between the inner and outer membranes. The T1 spanin gene was shown to complement the lambda Rz-Rz1- lysis defect, indicating that spanins function as Rz/Rz1 equivalents. The widespread presence of Rz/Rz1 or their spanin equivalents in phages of Gram-negative hosts suggests a strong selective advantage and that their role in the ecology of these phages is greater than that inferred from the mild laboratory phenotype.

Role of phages in the pathogenesis of Burkholderia, or 'Where are the toxin genes in Burkholderia phages?'

Most bacteria of the genus Burkholderia are soil- and rhizosphere-associated, and rhizosphere associated, noted for their metabolic plasticity in the utilization of a wide range of organic compounds as carbon sources. Many Burkholderia species are also opportunistic human and plant pathogens, and the distinction between environmental, plant, and human pathogens is not always clear. Burkholderia phages are not uncommon and multiple cryptic prophages are identifiable in the sequenced Burkholderia genomes. Phages have played a crucial role in the transmission of virulence factors among many important pathogens; however, the data do not yet support a significant correlation between phages and pathogenicity in the Burkholderia. This may be due to the role of Burkholderia as a 'versaphile' such that selection is occurring in several niches, including as a pathogen and in the context of environmental survival.

The complete genome sequence of Yersinia pseudotuberculosis IP31758, the causative agent of Far East scarlet-like fever

The first reported Far East scarlet-like fever (FESLF) epidemic swept the Pacific coastal region of Russia in the late 1950s. Symptoms of the severe infection included erythematous skin rash and desquamation, exanthema, hyperhemic tongue, and a toxic shock syndrome. The term FESLF was coined for the infection because it shares clinical presentations with scarlet fever caused by group A streptococci. The causative agent was later identified as Yersinia pseudotuberculosis, although the range of morbidities was vastly different from classical pseudotuberculosis symptoms. To understand the origin and emergence of the peculiar clinical features of FESLF, we have sequenced the genome of the FESLF-causing strain Y. pseudotuberculosis IP31758 and compared it with that of another Y. pseudotuberculosis strain, IP32953, which causes classical gastrointestinal symptoms. The unique gene pool of Y pseudotuberculosis IP31758 accounts for more than 260 strain-specific genes and introduces individual physiological capabilities and virulence determinants, with a significant proportion horizontally acquired that likely originated from Enterobacteriaceae and other soil-dwelling bacteria that persist in the same ecological niche. The mobile genome pool includes two novel plasmids phylogenetically unrelated to all currently reported Yersinia plasmids. An icm/dot type IVB secretion system, shared only with the intracellular persisting pathogens of the order Legionellales, was found on the larger plasmid and could contribute to scarlatinoid fever symptoms in patients due to the introduction of immunomodulatory and immunosuppressive capabilities. We determined the common and unique traits resulting from genome evolution and speciation within the genus Yersinia and drew a more accurate species border between Y. pseudotuberculosis and Y. pestis. In contrast to the lack of genetic diversity observed in the evolutionary young descending Y. pestis lineage, the population genetics of Y. pseudotuberculosis is more heterogenous. Both Y. pseudotuberculosis strains IP31758 and the previously sequenced Y. pseudotuberculosis strain IP32953 have evolved by the acquisition of specific plasmids and by the horizontal acquisition and incorporation of different genetic information into the chromosome, which all together or independently seems to potentially impact the phenotypic adaptation of these two strains.

Curated list of prokaryote viruses with fully sequenced genomes

Genome sequencing is of enormous importance for classification of prokaryote viruses and for understanding the evolution of these viruses. This survey covers 284 sequenced viruses for which a full description has been published and for which the morphology is known. This corresponds to 219 (4%) of tailed and 75 (36%) of tailless viruses of prokaryotes. The number of sequenced tailless viruses almost doubles if viruses of unknown morphology are counted. The sequences are from representatives of 15 virus families and three groups without family status, including eight taxa of archaeal viruses. Tailed phages, especially those with large genomes and hosts other than enterobacteria or lactococci, mycobacteria and pseudomonads, are vastly under investigated.

Salmonella phages and prophages--genomics and practical aspects

Numerous bacteriophages specific to Salmonella have been isolated or identified as part of host genome sequencing projects. Phylogenetic analysis of the sequenced phages, based on related protein content using CoreGenes, reveals that these viruses fall into five groupings (P27-like, P2-like, lambdoid, P22-like, and T7-like) and three outliers (epsilon15, KS7, and Felix O1). The P27 group is only represented by ST64B; the P2 group contains Fels-2, SopEphi, and PSP3; the lambdoid Salmonella phages include Gifsy-1, Gifsy-2, and Fels-1. The P22-like viruses include epsilon34, ES18, P22, ST104, and ST64T. The only member of the T7-like group is SP6. The properties of each of these phages are discussed, along with their role as agents of genetic exchange and as therapeutic agents and their involvement in phage typing.

Plasticity of the gene functions for DNA replication in the T4-like phages

We have completely sequenced and annotated the genomes of several relatives of the bacteriophage T4, including three coliphages (RB43, RB49 and RB69), three Aeromonas salmonicida phages (44RR2.8t, 25 and 31) and one Aeromonas hydrophila phage (Aeh1). In addition, we have partially sequenced and annotated the T4-like genomes of coliphage RB16 (a close relative of RB43), A. salmonicida phage 65, Acinetobacter johnsonii phage 133 and Vibrio natriegens phage nt-1. Each of these phage genomes exhibited a unique sequence that distinguished it from its relatives, although there were examples of genomes that are very similar to each other. As a group the phages compared here diverge from one another by several criteria, including (a) host range, (b) genome size in the range between approximately 160 kb and approximately 250 kb, (c) content and genetic organization of their T4-like genes for DNA metabolism, (d) mutational drift of the predicted T4-like gene products and their regulatory sites and (e) content of open-reading frames that have no counterparts in T4 or other known organisms (novel ORFs). We have observed a number of DNA rearrangements of the T4 genome type, some exhibiting proximity to putative homing endonuclease genes. Also, we cite and discuss examples of sequence divergence in the predicted sites for protein-protein and protein-nucleic acid interactions of homologues of the T4 DNA replication proteins, with emphasis on the diversity in sequence, molecular form and regulation of the phage-encoded DNA polymerase, gp43. Five of the sequenced phage genomes are predicted to encode split forms of this polymerase. Our studies suggest that the modular construction and plasticity of the T4 genome type and several of its replication proteins may offer resilience to mutation, including DNA rearrangements, and facilitate the adaptation of T4-like phages to different bacterial hosts in nature.