The Vertical Distribution and Diversity of Marine Bacteriophage at a Station off Southern California

Sixty-two bacteriophages were isolated on eight indigenous bacteria from a Pacific Ocean station spanning 887-m vertical depth, on two occasions between 1999 and 2000. On the basis of 16S rRNA sequences, six hosts were tentatively identified to be in the genus Vibrio and the other two were closely related to Altermonas macleodii (W9a) and Pseudoalteromonas spp. (W13a). Restriction fragment length polymorphism (RFLP) analysis of phage genomes using AccI and HapI showed that 16 phages infecting host C4a (Vibrio) displayed 14 unique RFLP patterns. However, identical phages infecting host C4b, C6a, and C6b (all Vibrio) were obtained from both the surface layer and the hypoxic zone at 850 m. Most phage isolates from the second year had a different RFLP pattern but shared genetic similarity to the phages infecting the same host from the previous year based on a hybridization study using phage genome probes. Cluster analysis of RFLP patterns and hybridization results also indicated that phages infecting the same or genetically related hosts, in general, shared higher degrees of homology in spite of the diverse RFLP patterns. Pulsed field gel electrophoresis (PFGE) analysis of native viral genomes indicated a range in genome size from less than 40 to 200 kb, and the dominant band shifted up by about 5-10 kb in the deep samples compared to the shallow ones. Hybridization of phage genome probes with total viral community DNA from various depths suggests these isolates, or at least some of their genes, represent a detectable portion of the natural viral community and were distributed throughout the water column. Thus, the results of this study demonstrated that the genetic diversity of bacteriophage in the ocean is far greater than that of their bacterial hosts. However, host range may have contributed to the evolution of the diverse phage population in the marine environment.     

Genomic characterization of four novel bacteriophages infecting the clinical pathogen Klebsiella pneumoniae

Bacteriophages are an invaluable source of novel genetic diversity. Sequencing of phage genomes can reveal new proteins with potential uses as biotechnological and medical tools, and help unravel the diversity of biological mechanisms employed by phages to take over the host during viral infection. Aiming to expand the available collection of phage genomes, we have isolated, sequenced, and assembled the genome sequences of four phages that infect the clinical pathogen Klebsiella pneumoniae: vB_KpnP_FBKp16, vB_KpnP_FBKp27, vB_KpnM_FBKp34, and Jumbo phage vB_KpnM_FBKp24. The four phages show very low (0–13%) identity to genomic phage sequences deposited in the GenBank database. Three of the four phages encode tRNAs and have a GC content very dissimilar to that of the host. Importantly, the genome sequences of the phages reveal potentially novel DNA packaging mechanisms as well as distinct clades of tubulin spindle and nucleus shell proteins that some phages use to compartmentalize viral replication. Overall, this study contributes to uncovering previously unknown virus diversity, and provides novel candidates for phage therapy applications against antibiotic-resistant K. pneumoniae infections.     

Efficient dilution-to-extinction isolation of novel virus–host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus–host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus–host systems in culture collections and genomic databases will improve both our understanding of virus–host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.Subject terms: Bacteriophages, Microbial ecology     

Marine phage genomics: what have we learned?

Marine phages are the most abundant and diverse form of life on the planet, and their genomes have been described as the largest untapped reservoir of genomic information. To date, however, the complete genome sequences of only 17 marine phage are known. Nevertheless, these genomes have revealed some interesting features, including the presence of photosynthetic genes in cyanophage and common patterns of genomic organization. Intriguing findings are also being made from studies of the uncultivated marine viral community genome ('metavirome'). The greatest challenge in interpreting the biology of these phages, and for making comparisons with their terrestrial counterparts, is the high proportion of unidentifiable open reading frames (approximately 60%). Future studies are likely to focus on sequencing more marine phage genomes from disparate hosts and diverse environments and on further basic studies of the biology of existing marine phages.     

Analysis of High-Throughput Sequencing and Annotation Strategies for Phage Genomes

Background

Bacterial viruses (phages) play a critical role in shaping microbial populations as they influence both host mortality and horizontal gene transfer. As such, they have a significant impact on local and global ecosystem function and human health. Despite their importance, little is known about the genomic diversity harbored in phages, as methods to capture complete phage genomes have been hampered by the lack of knowledge about the target genomes, and difficulties in generating sufficient quantities of genomic DNA for sequencing. Of the approximately 550 phage genomes currently available in the public domain, fewer than 5% are marine phage.

Methodology/Principal Findings

To advance the study of phage biology through comparative genomic approaches we used marine cyanophage as a model system. We compared DNA preparation methodologies (DNA extraction directly from either phage lysates or CsCl purified phage particles), and sequencing strategies that utilize either Sanger sequencing of a linker amplification shotgun library (LASL) or of a whole genome shotgun library (WGSL), or 454 pyrosequencing methods. We demonstrate that genomic DNA sample preparation directly from a phage lysate, combined with 454 pyrosequencing, is best suited for phage genome sequencing at scale, as this method is capable of capturing complete continuous genomes with high accuracy. In addition, we describe an automated annotation informatics pipeline that delivers high-quality annotation and yields few false positives and negatives in ORF calling.

Conclusions/Significance

These DNA preparation, sequencing and annotation strategies enable a high-throughput approach to the burgeoning field of phage genomics.     

Diversity and distribution of single-stranded DNA phages in the North Atlantic Ocean

Knowledge of marine phages is highly biased toward double-stranded DNA (dsDNA) phages; however, recent metagenomic surveys have also identified single-stranded DNA (ssDNA) phages in the oceans. Here, we describe two complete ssDNA phage genomes that were reconstructed from a viral metagenome from 80 m depth at the Bermuda Atlantic Time-series Study (BATS) site in the northwestern Sargasso Sea and examine their spatial and temporal distributions. Both genomes (SARssφ1 and SARssφ2) exhibited similarity to known phages of the Microviridae family in terms of size, GC content, genome organization and protein sequence. PCR amplification of the replication initiation protein (Rep) gene revealed narrow and distinct depth distributions for the newly described ssDNA phages within the upper 200 m of the water column at the BATS site. Comparison of Rep gene sequences obtained from the BATS site over time revealed changes in the diversity of ssDNA phages over monthly time scales, although some nearly identical sequences were recovered from samples collected 4 years apart. Examination of ssDNA phage diversity along transects through the North Atlantic Ocean revealed a positive correlation between genetic distance and geographic distance between sampling sites. Together, the data suggest fundamental differences between the distribution of these ssDNA phages and the distribution of known marine dsDNA phages, possibly because of differences in host range, host distribution, virion stability, or viral evolution mechanisms and rates. Future work needs to elucidate the host ranges for oceanic ssDNA phages and determine their ecological roles in the marine ecosystem.     

Marine phage genomics

Marine phages are the most abundant biological entities in the oceans. They play important roles in carbon cycling through marine food webs, gene transfer by transduction and conversion of hosts by lysogeny. The handful of marine phage genomes that have been sequenced to date, along with prophages in marine bacterial genomes, and partial sequencing of uncultivated phages are yielding glimpses of the tremendous diversity and physiological potential of the marine phage community. Common gene modules in diverse phages are providing the information necessary to make evolutionary comparisons. Finally, deciphering phage genomes is providing clues about the adaptive response of phages and their hosts to environmental cues.     

Quantitative molecular assessment of chimerism across tissues in marmosets and tamarins

Background

Periodic spacing of A-tracts (short runs of A or T) with the DNA helical period of ~10?C11?bp is characteristic of intrinsically bent DNA. In eukaryotes, the DNA bending is related to chromatin structure and nucleosome positioning. However, the physiological role of strong sequence periodicity detected in many prokaryotic genomes is not clear.

Results

We developed measures of intensity and persistency of DNA curvature-related sequence periodicity and applied them to prokaryotic chromosomes and phages. The results indicate that strong periodic signals present in chromosomes are generally absent in phage genomes. Moreover, chromosomes containing prophages are less likely to possess a persistent periodic signal than chromosomes with no prophages.

Conclusions

Absence of DNA curvature-related sequence periodicity in phages could arise from constraints associated with DNA packaging in the viral capsid. Lack of prophages in chromosomes with persistent periodic signal suggests that the sequence periodicity and concomitant DNA curvature could play a role in protecting the chromosomes from integration of phage DNA.     

Conserved genomes of ΦKMV-like bacteriophages (T7 supergroup) active on Pseudomonas aeruginosa

The T7-like ΦKMV bacteriophage active on Pseudomonas aeruginosa was previously isolated by us and shown to have DNA resistant to many endonucleases. A loss of sensitive sites might be a consequence of a long ΦKMV evolution on different hosts. To elucidate, whether this trait is shared by other similar phages, several new ΦKMV-like phages were isolated from different sources and compared. All studied ΦKMV-like phages formed three groups, insignificantly differing in the number and localization of endonuclease-sensitive DNA sites. This confirms that the present-day phages of this species have highly conserved genomes. Mutational “restoration” of the lost sites may be restricted by a lethal effect. The ΦKMV-like phages were shown for the first time to increase the rate of in vitro accumulation of giant KZ-like phages of P. aeruginosa. This effect is characteristic only of ΦKMV-like phages.     

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.     

Reconstructing viral genomes from the environment using fosmid clones: the case of haloviruses

Background

Metaviriomes, the viral genomes present in an environment, have been studied by direct sequencing of the viral DNA or by cloning in small insert libraries. The short reads generated by both approaches make it very difficult to assemble and annotate such flexible genomic entities. Many environmental viruses belong to unknown groups or prey on uncultured and little known cellular lineages, and hence might not be present in databases.

Methodology and Principal Findings

Here we have used a different approach, the cloning of viral DNA into fosmids before sequencing, to obtain natural contigs that are close to the size of a viral genome. We have studied a relatively low diversity extreme environment: saturated NaCl brines, which simplifies the analysis and interpretation of the data. Forty-two different viral genomes were retrieved, and some of these were almost complete, and could be tentatively identified as head-tail phages (Caudovirales).

Conclusions and Significance

We found a cluster of phage genomes that most likely infect Haloquadratum walsbyi, the square archaeon and major component of the community in these hypersaline habitats. The identity of the prey could be confirmed by the presence of CRISPR spacer sequences shared by the virus and one of the available strain genomes. Other viral clusters detected appeared to prey on the Nanohaloarchaea and on the bacterium Salinibacter ruber, covering most of the diversity of microbes found in this type of environment. This approach appears then as a viable alternative to describe metaviriomes in a much more detailed and reliable way than by the more common approaches based on direct sequencing. An example of transfer of a CRISPR cluster including repeats and spacers was accidentally found supporting the dynamic nature and frequent transfer of this peculiar prokaryotic mechanism of cell protection.     

Preliminary Characterization of Bacteriophages Infecting the Thermophilic Actinomycete Thermomonospora

Bacteriophages lysing strains of Thermomonospora alba and T. fusca were isolated, following specific enrichment, from vegetable composts. Four Thermomonospora phages were distinguished by plaque morphology and host range. Electron microscopy of phage particles, termperature inactivation profiles, and electrophoretic analyses of major virion proteins and genomic DNA were used in the comparison and initial characterization of these phages. The four phages studied possessed polyhedral heads and long tails; genomes were linear double-stranded DNA molecules, 35 to 45 kilobases in length, which probably contain cohesive ends. Transfection of Thermomonospora protoplasts with purified genomic DNA from one of the phages was demonstrated.     

Codon Bias is a Major Factor Explaining Phage Evolution in Translationally Biased Hosts

The size and diversity of bacteriophage populations require methodologies to quantitatively study the landscape of phage differences. Statistical approaches are confronted with small genome sizes forbidding significant single-phage analysis, and comparative methods analyzing full phage genomes represent an alternative but they are of difficult interpretation due to lateral gene transfer, which creates a mosaic spectrum of related phage species. Based on a large-scale codon bias analysis of 116 DNA phages hosted by 11 translationally biased bacteria belonging to different phylogenetic families, we observe that phage genomes are almost always under codon selective pressure imposed by translationally biased hosts, and we propose a classification of phages with translationally biased hosts which is based on adaptation patterns. We introduce a computational method for comparing phages sharing homologous proteins, possibly accepted by different hosts. We observe that throughout phages, independently from the host, capsid genes appear to be the most affected by host translational bias. For coliphages, genes involved in virion morphogenesis, host interaction and ssDNA binding are also affected by adaptive pressure. Adaptation affects long and small phages in a significant way. We analyze in more detail the Microviridae phage space to illustrate the potentiality of the approach. The small number of directions in adaptation observed in phages grouped around ϕX174 is discussed in the light of functional bias. The adaptation analysis of the set of Microviridae phages defined around ϕMH2K shows that phage classification based on adaptation does not reflect bacterial phylogeny. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Single-cell enabled comparative genomics of a deep ocean SAR11 bathytype

Bacterioplankton of the SAR11 clade are the most abundant microorganisms in marine systems, usually representing 25% or more of the total bacterial cells in seawater worldwide. SAR11 is divided into subclades with distinct spatiotemporal distributions (ecotypes), some of which appear to be specific to deep water. Here we examine the genomic basis for deep ocean distribution of one SAR11 bathytype (depth-specific ecotype), subclade Ic. Four single-cell Ic genomes, with estimated completeness of 55%–86%, were isolated from 770 m at station ALOHA and compared with eight SAR11 surface genomes and metagenomic datasets. Subclade Ic genomes dominated metagenomic fragment recruitment below the euphotic zone. They had similar COG distributions, high local synteny and shared a large number (69%) of orthologous clusters with SAR11 surface genomes, yet were distinct at the 16S rRNA gene and amino-acid level, and formed a separate, monophyletic group in phylogenetic trees. Subclade Ic genomes were enriched in genes associated with membrane/cell wall/envelope biosynthesis and showed evidence of unique phage defenses. The majority of subclade Ic-specfic genes were hypothetical, and some were highly abundant in deep ocean metagenomic data, potentially masking mechanisms for niche differentiation. However, the evidence suggests these organisms have a similar metabolism to their surface counterparts, and that subclade Ic adaptations to the deep ocean do not involve large variations in gene content, but rather more subtle differences previously observed deep ocean genomic data, like preferential amino-acid substitutions, larger coding regions among SAR11 clade orthologs, larger intergenic regions and larger estimated average genome size.     

Genomic Diversity of Phages Infecting Probiotic Strains of Lactobacillus paracasei

Strains of the Lactobacillus casei group have been extensively studied because some are used as probiotics in foods. Conversely, their phages have received much less attention. We analyzed the complete genome sequences of five L. paracasei temperate phages: C_L1, C_L2, iLp84, iLp1308, and iA2. Only phage iA2 could not replicate in an indicator strain. The genome lengths ranged from 34,155 bp (iA2) to 39,474 bp (C_L1). Phages iA2 and iLp1308 (34,176 bp) possess the smallest genomes reported, thus far, for phages of the L. casei group. The GC contents of the five phage genomes ranged from 44.8 to 45.6%. As observed with many other phages, their genomes were organized as follows: genes coding for DNA packaging, morphogenesis, lysis, lysogeny, and replication. Phages C_L1, C_L2, and iLp1308 are highly related to each other. Phage iLp84 was also related to these three phages, but the similarities were limited to gene products involved in DNA packaging and structural proteins. Genomic fragments of phages C_L1, C_L2, iLp1308, and iLp84 were found in several genomes of L. casei strains. Prophage iA2 is unrelated to these four phages, but almost all of its genome was found in at least four L. casei strains. Overall, these phages are distinct from previously characterized Lactobacillus phages. Our results highlight the diversity of L. casei phages and indicate frequent DNA exchanges between phages and their hosts.     

Analysis of six prophages in Lactococcus lactis IL1403: different genetic structure of temperate and virulent phage populations

We report the genetic organisation of six prophages present in the genome of Lactococcus lactis IL1403. The three larger prophages (36–42 kb), belong to the already described P335 group of temperate phages, whereas the three smaller ones (13–15 kb) are most probably satellites relying on helper phage(s) for multiplication. These data give a new insight into the genetic structure of lactococcal phage populations. P335 temperate phages have variable genomes, sharing homology over only 10–33% of their length. In contrast, virulent phages have highly similar genomes sharing homology over >90% of their length. Further analysis of genetic structure in all known groups of phages active on other bacterial hosts such as Escherichia coli, Bacillus subtilis, Mycobacterium and Streptococcus thermophilus confirmed the existence of two types of genetic structure related to the phage way of life. This might reflect different intensities of horizontal DNA exchange: low among purely virulent phages and high among temperate phages and their lytic homologues. We suggest that the constraints on genetic exchange among purely virulent phages reflect their optimal genetic organisation, adapted to a more specialised and extreme form of parasitism than temperate/lytic phages.     

Evaluation of a transposase protocol for rapid generation of shotgun high-throughput sequencing libraries from nanogram quantities of DNA

Construction of DNA fragment libraries for next-generation sequencing can prove challenging, especially for samples with low DNA yield. Protocols devised to circumvent the problems associated with low starting quantities of DNA can result in amplification biases that skew the distribution of genomes in metagenomic data. Moreover, sample throughput can be slow, as current library construction techniques are time-consuming. This study evaluated Nextera, a new transposon-based method that is designed for quick production of DNA fragment libraries from a small quantity of DNA. The sequence read distribution across nine phage genomes in a mock viral assemblage met predictions for six of the least-abundant phages; however, the rank order of the most abundant phages differed slightly from predictions. De novo genome assemblies from Nextera libraries provided long contigs spanning over half of the phage genome; in four cases where full-length genome sequences were available for comparison, consensus sequences were found to match over 99% of the genome with near-perfect identity. Analysis of areas of low and high sequence coverage within phage genomes indicated that GC content may influence coverage of sequences from Nextera libraries. Comparisons of phage genomes prepared using both Nextera and a standard 454 FLX Titanium library preparation protocol suggested that the coverage biases according to GC content observed within the Nextera libraries were largely attributable to bias in the Nextera protocol rather than to the 454 sequencing technology. Nevertheless, given suitable sequence coverage, the Nextera protocol produced high-quality data for genomic studies. For metagenomics analyses, effects of GC amplification bias would need to be considered; however, the library preparation standardization that Nextera provides should benefit comparative metagenomic analyses.     

Genetic diversity in cultured and wild marine cyanomyoviruses reveals phosphorus stress as a strong selective agent

Viruses that infect marine cyanobacteria–cyanophages–often carry genes with orthologs in their cyanobacterial hosts, and the frequency of these genes can vary with habitat. To explore habitat-influenced genomic diversity more deeply, we used the genomes of 28 cultured cyanomyoviruses as references to identify phage genes in three ocean habitats. Only about 6–11% of genes were consistently observed in the wild, revealing high gene-content variability in these populations. Numerous shared phage/host genes differed in relative frequency between environments, including genes related to phosphorous acquisition, photorespiration, photosynthesis and the pentose phosphate pathway, possibly reflecting environmental selection for these genes in cyanomyovirus genomes. The strongest emergent signal was related to phosphorous availability; a higher fraction of genomes from relatively low-phosphorus environments–the Sargasso and Mediterranean Sea–contained host-like phosphorus assimilation genes compared with those from the N. Pacific Gyre. These genes are known to be upregulated when the host is phosphorous starved, a response mediated by pho box motifs in phage genomes that bind a host regulatory protein. Eleven cyanomyoviruses have predicted pho boxes upstream of the phosphate-acquisition genes pstS and phoA; eight of these have a conserved cyanophage-specific gene (PhCOG173) between the pho box and pstS. PhCOG173 is also found upstream of other shared phage/host genes, suggesting a unique regulatory role. Pho boxes are found upstream of high light-inducible (hli) genes in cyanomyoviruses, suggesting that this motif may have a broader role than regulating phosphorous-stress responses in infected hosts or that these hlis are involved in the phosphorous-stress response.     

A New Microviridae Phage Isolated from a Failed Biotechnological Process Driven by Escherichia coli

Bacteriophages are present in every environment that supports bacterial growth, including manmade ecological niches. Virulent phages may even slow or, in more severe cases, interrupt bioprocesses driven by bacteria. Escherichia coli is one of the most widely used bacteria for large-scale bioprocesses; however, literature describing phage-host interactions in this industrial context is sparse. Here, we describe phage MED1 isolated from a failed industrial process. Phage MED1 (Microviridae family, with a single-stranded DNA [ssDNA] genome) is highly similar to the archetypal phage phiX174, sharing >95% identity between their genomic sequences. Whole-genome phylogenetic analysis of 52 microvirus genomes from public databases revealed three genotypes (alpha3, G4, and phiX174). Phage MED1 belongs to the phiX174 group. We analyzed the distribution of single nucleotide variants in MED1 and 18 other phiX174-like genomes and found that there are more missense mutations in genes G, B, and E than in the other genes of these genomes. Gene G encodes the spike protein, involved in host attachment. The evolution of this protein likely results from the selective pressure on phages to rapidly adapt to the molecular diversity found at the surface of their hosts.     

Genetic Variation of Lactobacillus delbrueckii subsp. lactis Bacteriophages Isolated from Cheese Processing Plants in Finland

The genomes of four Lactobacillus delbrueckii subsp. lactis bacteriophages were characterized by restriction endonuclease mapping, Southern hybridization, and heteroduplex analysis. The phages were isolated from different cheese processing plants in Finland between 1950 and 1972. All four phages had a small isometric head and a long noncontractile tail. Two different types of genome (double-stranded DNA) organization existed among the different phages, the pac type and the cos type, corresponding to alternative types of phage DNA packaging. Three phages belonged to the pac type, and a fourth was a cos-type phage. The pac-type phages were genetically closely related. In the genomes of the pac-type phages, three putative insertion/deletions (0.7 to 0.8 kb, 1.0 kb, and 1.5 kb) and one other region (0.9 kb) containing clustered base substitutions were discovered and localized. At the phenotype level, three main differences were observed among the pac-type phages. These concerned two minor structural proteins and the efficiency of phage DNA packaging. The genomes of the pac-type phages showed only weak homology with that of the cos-type phage. Phage-related DNA, probably a defective prophage, was located in the chromosome of the host strain sensitive to the cos-type phage. This DNA exhibited homology under stringent conditions to the pac-type phages.