Identification of Huge Phages from Wastewater Metagenomes

Huge phages have genomes larger than 200 kilobases, which are particularly interesting for their genetic inventory and evolution. We screened 165 wastewater metagenomes for the presence of viral sequences. After identifying over 600 potential huge phage genomes, we reduced the dataset using manual curation by excluding viral contigs that did not contain viral protein-coding genes or consisted of concatemers of several small phage genomes. This dataset showed seven fully annotated huge phage genomes. The phages grouped into distinct phylogenetic clades, likely forming new genera and families. A phylogenomic analysis between our huge phages and phages with smaller genomes, i.e., less than 200 kb, supported the hypothesis that huge phages have undergone convergent evolution. The genomes contained typical phage protein-coding genes, sequential gene cassettes for metabolic pathways, and complete inventories of tRNA genes covering all standard and rare amino acids. Our study showed a pipeline for huge phage analyses that may lead to new enzymes for therapeutic or biotechnological applications.

Ubiquitous, B12-dependent virioplankton utilizing ribonucleotide-triphosphate reductase demonstrate interseasonal dynamics and associate with a diverse range of bacterial hosts in the pelagic ocean

Through infection and lysis of their coexisting bacterial hosts, viruses impact the biogeochemical cycles sustaining globally significant pelagic oceanic ecosystems. Currently, little is known of the ecological interactions between lytic viruses and their bacterial hosts underlying these biogeochemical impacts at ecosystem scales. This study focused on populations of lytic viruses carrying the B12-dependent Class II monomeric ribonucleotide reductase (RNR) gene, ribonucleotide-triphosphate reductase (Class II RTPR), documenting seasonal changes in pelagic virioplankton and bacterioplankton using amplicon sequences of Class II RTPR and the 16S rRNA gene, respectively. Amplicon sequence libraries were analyzed using compositional data analysis tools that account for the compositional nature of these data. Both virio- and bacterioplankton communities responded to environmental changes typically seen across seasonal cycles as well as shorter term upwelling-downwelling events. Defining Class II RTPR-carrying viral populations according to major phylogenetic clades proved a more robust means of exploring virioplankton ecology than operational taxonomic units defined by percent sequence homology. Virioplankton Class II RTPR populations showed positive associations with a broad phylogenetic diversity of bacterioplankton including dominant taxa within pelagic oceanic ecosystems such as Prochlorococcus and SAR11. Temporal changes in Class II RTPR virioplankton, occurring as both free viruses and within infected cells, indicated possible viral-host pairs undergoing sustained infection and lysis cycles throughout the seasonal study. Phylogenetic relationships inferred from Class II RTPR sequences mirrored ecological patterns in virio- and bacterioplankton populations demonstrating possible genome to phenome associations for an essential viral replication gene.

Evolution of the T4 phage virion is driven by selection pressure from non-bacterial factors

Bacteriophages colonize animal and human bodies, propagating on sensitive bacteria that are symbionts, commensals, or pathogens of animals and humans. T4-like phages are dependent on abundant symbionts such as Escherichia coli, commonly present in animal and human gastrointestinal (GI) tracts. Bacteriophage T4 is one of the most complex viruses, and its intricate structure, particularly the capsid head protecting the phage genome, likely contributes substantially to the overall phage fitness in diverse environments. We investigated how individual head proteins-gp24, Hoc, and Soc-affect T4 phage survival under pressure from non-bacterial factors. We constructed a panel of T4 phage variants defective in these structural proteins: T4∆Soc, T4∆24byp24, T4∆Hoc∆Soc, T4∆Hoc∆24byp24, T4∆Soc∆24byp24, and T4∆Hoc∆Soc∆24byp24 (byp = bypass). These variants were investigated for their sensitivity to selected environmental conditions relevant to the microenvironment of the GI tract, including pH, temperature, and digestive enzymes. The simple and "primitive" structure of the phage capsid (∆24byp24) was significantly less stable at low pH and more sensitive to inactivation by digestive enzymes, and the simultaneous lack of gp24 and Soc resulted in a notable decrease in phage activity at 37°C. Gp24 was also found to be highly resistant to thermal and chemical denaturation. Thus, gp24, which was acquired relatively late in evolution, seems to play a key role in T4 withstanding environmental conditions, including those related to the animal/human GI tract, and Soc is a molecular glue that enhances this protective effect. IMPORTANCE Bacteriophages are important components of animal and human microbiota, particularly in the gastrointestinal tract, where they dominate the viral community and contribute to shaping microbial balance. However, interactions with bacterial hosts are not the only element of the equation in phage survival-phages inhabiting the GI tract are constantly exposed to increased temperature, pH fluctuations, or digestive enzymes, which raises the question of whether and how the complex structure of phage capsids contributes to their persistence in the specific microenvironment of human/animal bodies. Here we address this phage-centric perspective, identifying the role of individual head proteins in T4 phage survival in GI tract conditions. The selection pressure driving the evolution of T4-like phages could have come from the external environment that affects phage virions with increased temperature and variable pH; it is possible that in the local microenvironment along the GI tract, the phage benefits from stability-protecting proteins.

Spatiotemporal Shift of T4-Like Phage Community Structure in the Three Largest Estuaries of China

Estuaries are one of the most highly productive and economically important ecosystems at the continent-ocean interface. Estuary productivity is largely determined by the microbial community structure and activity. Viruses are major agents of microbial mortality and are key drivers of global geochemical cycles. However, the taxonomic diversity of viral communities and their spatial-temporal distribution in estuarine ecosystems have been poorly studied. In this study, we investigated the T4-like viral community composition at three major Chinese estuaries in winter and in summer. Diverse T4-like viruses, which were divided into three main clusters (Clusters I to III), were revealed. The Marine Group of Cluster III, with seven identified subgroups, was the most dominant (averaging 76.5% of the total sequences) in the Chinese estuarine ecosystems. Significant variations of T4-like viral community composition were observed among estuaries and seasons, with higher diversity occurring in winter. Among various environmental variables, temperature was a main driver of the viral communities. This study demonstrates viral assemblage diversification and seasonality in Chinese estuarine ecosystems. IMPORTANCE Viruses are ubiquitous but largely uncharacterized members of aquatic environments that cause significant mortality in microbial communities. Recent large-scale oceanic projects have greatly advanced our understanding of viral ecology in marine environments, but those studies mostly focused on oceanic regions. There have yet to be spatiotemporal studies of viral communities in estuarine ecosystems, which are unique habitats that play a significant role in global ecology and biogeochemistry. This work is the first comprehensive study that provides a detailed picture of the spatial and seasonal variation of viral communities (specifically, T4-like viral communities) in three major estuarine ecosystems in China. These findings provide much-needed knowledge regarding estuarine viral ecosystems, which currently lags in oceanic ecosystem research.

Uncovering the determinants of model Escherichia coli strain C600 susceptibility and resistance to lytic T4-like and T7-like phage

As antimicrobial resistance (AMR) continues to increase, the therapeutic use of phages has re-emerged as an attractive alternative. However, knowledge of phage resistance development and bacterium-phage interaction complexity are still not fully interpreted. In this study, two lytic T4-like and T7-like phage infecting model Escherichia coli strain C600 are selected, and host genetic determinants involved in phage susceptibility and resistance are also identified using TraDIS strategy. Isolation and identification of the lytic T7-like show that though it belongs to the phage T7 family, genes encoding replication and transcription protein exhibit high differences. The TraDIS results identify a huge number of previously unidentified genes involved in phage infection, and a subset (six in susceptibility and nine in resistance) are shared under pressure of the two kinds of lytic phage. Susceptible gene wbbL has the highest value and implies the important role in phage susceptibility. Importantly, two susceptible genes QseE (QseE/QseF) and RstB (RstB/RstA), encoding the similar two-component system sensor histidine kinase (HKs), also identified. Conversely and strangely, outer membrane protein gene ompW, unlike the gene ompC encoding receptor protein of T4 phage, was shown to provide phage resistance. Overall, this study exploited a genome-wide fitness assay to uncover susceptibility and resistant genes, even the shared genes, important for the E. coli strain of both most popular high lytic T4-like and T7-like phages. This knowledge of the genetic determinants can be further used to analysis the behind function signatures to screen the potential agents to aid phage killing of MDR pathogens, which will greatly be valuable in improving the phage therapy outcome in fighting with microbial resistance.

Effects of phytoplankton, viral communities, and warming on free-living and particle-associated marine prokaryotic community structure

Free-living and particle-associated marine prokaryotes have physiological, genomic, and phylogenetic differences, yet factors influencing their temporal dynamics remain poorly constrained. In this study, we quantify the entire microbial community composition monthly over several years, including viruses, prokaryotes, phytoplankton, and total protists, from the San-Pedro Ocean Time-series using ribosomal RNA sequencing and viral metagenomics. Canonical analyses show that in addition to physicochemical factors, the double-stranded DNA viral community is the strongest factor predicting free-living prokaryotes, explaining 28% of variability, whereas the phytoplankton (via chloroplast 16S rRNA) community is strongest with particle-associated prokaryotes, explaining 31% of variability. Unexpectedly, protist community explains little variability. Our findings suggest that biotic interactions are significant determinants of the temporal dynamics of prokaryotes, and the relative importance of specific interactions varies depending on lifestyles. Also, warming influenced the prokaryotic community, which largely remained oligotrophic summer-like throughout 2014-15, with cyanobacterial populations shifting from cold-water ecotypes to warm-water ecotypes.

Novel PhoH-encoding vibriophages with lytic activity against environmental Vibrio strains

Cholera is a devastating diarrheal disease that accounts for more than 10% of children's lives worldwide, but its treatment is hampered by a rise in antibiotic resistance. One promising alternative to antibiotic therapy is the use of bacteriophages to treat antibiotic-resistant cholera infections, and control Vibrio cholera in clinical cases and in the environment, respectively. Here, we report four novel, closely related environmental myoviruses, VP4, VP6, VP18, and VP24, which we isolated from two environmental toxigenic Vibrio cholerae strains from river Kuja and Usenge beach in Kenya. High-throughput sequencing followed by bioinformatics analysis indicated that the genomes of the four bacteriophages have closely related sequences, with sizes of 148,180 bp, 148,181 bp, 148,179 bp, and 148,179 bp, and a G + C content of 36.4%. The four genomes carry the phoH gene, which is overrepresented in marine cyanophages. The isolated phages displayed a lytic activity against 15 environmental, as well as one clinical, Vibrio cholerae strains. Thus, these novel lytic vibriophages represent potential biocontrol candidates for water decontamination against pathogenic Vibrio cholerae and ought to be considered for future studies of phage therapy.

Bacterial Consumption of T4 Phages

The bacterial consumption of viruses not been reported on as of yet even though bacteria feed on almost anything. Viruses are widely distributed but have no acknowledged active biocontrol. Viral biomass undoubtedly reintegrates trophic cycles; however, the mechanisms of this phase still remain unknown. ¹³C-labelled T4 phages monitor the increase of the density of the bacterial DNA concomitant with the decrease of plaque forming units. We used ¹²C T4 phages as a control. T4 phage disappearance in wastewater sludge was found to occur mainly through predation by Aeromonadacea. Phage consumption also favours significant in situ bacterial growth. Furthermore, an isolated strain of Aeromonas was observed to grow on T4 phages as sole the source of carbon, nitrogen, and phosphorus. Bacterial species are capable of consuming bacteriophages in situ, which is likely a widespread and underestimated type of biocontrol. This assay is anticipated as a starting point for harnessing the bacterial potential in limiting the diffusion of harmful viruses within environments such as in the gut or in water.

Role of Phylogenetic Structure in the Dynamics of Coastal Viral Assemblages

Marine microbes, including viruses, are an essential part of the marine ecosystem, forming the base of the food web and driving biogeochemical cycles. Within this system, the composition of viral assemblages changes markedly with time, and some of these changes are repeatable through time; however, the extent to which these dynamics are reflected within versus among evolutionarily related groups of viruses is largely unexplored. To examine these dynamics, changes in the composition of two groups of ecologically important viruses and communities of their potential hosts were sampled every 2 weeks for 13 months at a coastal site in British Columbia, Canada. We sequenced two marker genes for viruses-the gene encoding the major capsid protein of T4-like phages and their relatives (gp23) and the RNA-dependent RNA polymerase (RdRp) gene of marnavirus-like RNA viruses-as well as marker genes for their bacterial and eukaryotic host communities, the genes encoding 16S rRNA and 18S rRNA. There were strong lagged correlations between viral diversity and community similarity of putative hosts, implying that the viruses influenced the composition of the host communities. The results showed that for both viral assemblages, the dominant clusters of phylogenetically related viruses shifted over time, and this was correlated with environmental changes. Viral clusters contained many ephemeral taxa and few persistent taxa, but within a viral assemblage, the ephemeral and persistent taxa were closely related, implying ecological dynamics within these clusters. Furthermore, these dynamics occurred in both the RNA and DNA viral assemblages surveyed, implying that this structure is common in natural viral assemblages.IMPORTANCE Viruses are major agents of microbial mortality in marine systems, yet little is known about changes in the composition of viral assemblages in relation to those of the microbial communities that they infect. Here, we sampled coastal seawater every 2 weeks for 1 year and used high-throughput sequencing of marker genes to follow changes in the composition of two groups of ecologically important viruses, as well as the communities of bacteria and protists that serve as their respective hosts. Different subsets of genetically related viruses dominated at different times. These results demonstrate that although the genetic composition of viral assemblages is highly dynamic temporally, for the most part the shuffling of genotypes occurs within a few clusters of phylogenetically related viruses. Thus, it appears that even in temperate coastal waters with large seasonal changes, the highly dynamic shuffling of viral genotypes occurs largely within a few subsets of related individuals.

Influence of Non-canonical DNA Bases on the Genomic Diversity of Tevenvirinae

The Tevenvirinae viruses are some of the most common viruses on Earth. Representatives of this subfamily have long been used in the molecular biology studies as model organisms – since the emergence of the discipline. Tevenvirinae are promising agents for phage therapy in animals and humans, since their representatives have only lytic life cycle and many of their host bacteria are pathogens. As confirmed experimentally, some Tevenvirinae have non-canonical DNA bases. Non-canonical bases can play an essential role in the diversification of closely related viruses. The article performs a comparative and evolutionary analysis of Tevenvirinae genomes and components of Tevenvirinae genomes. A comparative analysis of these genomes and the genes associated with the synthesis of non-canonical bases allows us to conclude that non-canonical bases have a major influence on the divergence of Tevenvirinae viruses within the same habitats. Supposedly, Tevenvirinae developed a strategy for changing HGT frequency in individual populations, which was based on the accumulation of proteins for the synthesis of non-canonical bases and proteins that used those bases as substrates. Owing to this strategy, ancestors of Tevenvirinae with the highest frequency of HGT acquired genes that allowed them to exist in a certain niche, and ancestors with the lowest HGT frequency preserved the most adaptive of those genes. Given the origin and characteristics of genes associated with the synthesis of non-canonical bases in Tevenvirinae, one can assume that other phages may have similar strategies. The article demonstrates the dependence of genomic diversity of closely related Tevenvirinae on non-canonical bases.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

Isolation, characterization and analysis of bacteriophages from the haloalkaline lake Elmenteita, Kenya

As a step towards better understanding of diversity and biology of phages and their hosts in haloalkaline Lake Elmenteita, phages were isolated from sediment samples and overlying water using indigenous bacteria as hosts. 17 seemingly different phages of diverse morphotypes with different dimensions and partly exhibiting remarkably unusual ultrastructures were revealed by transmission electron microscopy. 12 clonal phage isolates were further characterized. Infection capability of the phages was optimum at 30-35°C and in alkali condition with optimum at pH 10-12. Structural protein profiles and restriction fragment length polymorphism analyses patterns were distinct for each of the phage type. Complete nucleotide sequences of phages vB-VmeM-32, vB_EauS-123 and vB_BhaS-171 genomes varied in size from 30,926-199,912 bp and G + C content of between 36.25-47.73%. A range of 56-260 potential open reading frames were identified and annotated. The results showed that the 12 phages were distinct from each other and confirmed the presence and diversity of phages in extreme environment of haloalkaline Lake Elmenteita. The phages were deposited at the German Collection of Microorganisms and Cell Cultures and three of their genomes uploaded to NCBI GenBank.

The Environment and Cyanophage Diversity: Insights From Environmental Sequencing of DNA Polymerase

Globally distributed and abundant cyanophages in the family Myoviridae have dsDNA genomes with variable gene content, including host-derived auxiliary metabolic genes (AMGs) that potentially can facilitate viral replication. However, it is not well understood how this variation in gene content interacts with environmental variables to shape cyanomyovirus communities. This project correlated the genetic repertoire of cyanomyoviruses with their phyologeny, and investigated cyanomyovirus ecotype distribution as a function of environmental conditions across locations and seasons. Reference cyanomyovirus genomes were compared for their overlap in gene content to infer phyologenetic distances, and these distances were compared to distances calculated based on DNA polymerase (gp43) gene sequences. In turn, gp43 partial gene sequences amplified from natural cyanophage communities were used to describe cyanomyovirus community composition and to assess the relationship between environmental variables. The results showed the following: (1) DNA polymerase gene phylogeny generally correlated with the similarity in gene content among reference cyanomyoviruses, and thus can be used to describe environmental cyanomyovirus communities; (2) spatial and seasonal patterns in cyanomyovirus communities were related to environmental variables; (3) salinity and temperature, combined with nutrient concentration were predictors of cyanomyovirus richness, diversity and community composition. This study shows that environmental variables shape viral communities by drawing on a diverse seed bank of viral genotypes. From these results it is evident that that viral ecotypes with their corresponding genetic repertoires underlie selection pressures. However, the mechanisms involved in selecting for specific viral genotypes remain to be fully understood.

Diversity of Active Viral Infections within the Sphagnum Microbiome

Sphagnum-dominated peatlands play an important role in global carbon storage and represent significant sources of economic and ecological value. While recent efforts to describe microbial diversity and metabolic potential of the Sphagnum microbiome have demonstrated the importance of its microbial community, little is known about the viral constituents. We used metatranscriptomics to describe the diversity and activity of viruses infecting microbes within the Sphagnum peat bog. The vegetative portions of six Sphagnum plants were obtained from a peatland in northern Minnesota, and the total RNA was extracted and sequenced. Metatranscriptomes were assembled and contigs were screened for the presence of conserved virus marker genes. Using bacteriophage capsid protein gp23 as a marker for phage diversity, we identified 33 contigs representing undocumented phages that were active in the community at the time of sampling. Similarly, RNA-dependent RNA polymerase and the nucleocytoplasmic large DNA virus (NCLDV) major capsid protein were used as markers for single-stranded RNA (ssRNA) viruses and NCLDV, respectively. In total, 114 contigs were identified as originating from undescribed ssRNA viruses, 22 of which represent nearly complete genomes. An additional 64 contigs were identified as being from NCLDVs. Finally, 7 contigs were identified as putative virophage or polinton-like viruses. We developed co-occurrence networks with these markers in relation to the expression of potential-host housekeeping gene rpb1 to predict virus-host relationships, identifying 13 groups. Together, our approach offers new tools for the identification of virus diversity and interactions in understudied clades and suggests that viruses may play a considerable role in the ecology of the Sphagnum microbiome.IMPORTANCE Sphagnum-dominated peatlands play an important role in maintaining atmospheric carbon dioxide levels by modifying conditions in the surrounding soil to favor the growth of Sphagnum over that of other plant species. This lowers the rate of decomposition and facilitates the accumulation of fixed carbon in the form of partially decomposed biomass. The unique environment produced by Sphagnum enriches for the growth of a diverse microbial consortia that benefit from and support the moss's growth, while also maintaining the hostile soil conditions. While a growing body of research has begun to characterize the microbial groups that colonize Sphagnum, little is currently known about the ecological factors that constrain community structure and define ecosystem function. Top-down population control by viruses is almost completely undescribed. This study provides insight into the significant viral influence on the Sphagnum microbiome and identifies new potential model systems to study virus-host interactions in the peatland ecosystem.

Novel T4 bacteriophages associated with black band disease in corals

Research into causative agents underlying coral disease have focused primarily on bacteria, whereas potential roles of viruses have been largely unaddressed. Bacteriophages may contribute to diseases through the lysogenic introduction of virulence genes into bacteria, or prevent diseases through lysis of bacterial pathogens. To identify candidate phages that may influence the pathogenicity of black band disease (BBD), communities of bacteria (16S rRNA) and T4-bacteriophages (gp23) were simultaneously profiled with amplicon sequencing among BBD-lesions and healthy-coral-tissue of Montipora hispida, as well as seawater (study site: the central Great Barrier Reef). Bacterial community compositions were distinct among BBD-lesions, healthy coral tissue and seawater samples, as observed in previous studies. Surprisingly, however, viral beta diversities based on both operational taxonomic unit (OTU)-compositions and overall viral community compositions of assigned taxa did not differ statistically between the BBD-lesions and healthy coral tissue. Nonetheless, relative abundances of three bacteriophage OTUs, affiliated to Cyanophage PRSM6 and Prochlorococcus phages P-SSM2, were significantly higher in BBD-lesions than in healthy tissue. These OTUs associated with BBD samples suggest the presence of bacteriophages that infect members of the cyanobacteria-dominated BBD community, and thus have potential roles in BBD pathogenicity.

Genomic analysis of Acinetobacter baumannii prophages reveals remarkable diversity and suggests profound impact on bacterial virulence and fitness

The recent nomination by the World Health Organization of Acinetobacter baumannii as the number one priority pathogen for the development of new antibiotics is a direct consequence of its fast evolution of pathogenicity, and in particular of multidrug resistance. While the development of new antibiotics is critical, understanding the mechanisms behind the crescent bacterial antibiotic resistance is equally relevant. Often, resistance and other bacterial virulence elements are contained on highly mobile pieces of DNA that can easily spread to other bacteria. Prophages are one of the mediators of this form of gene transfer, and have been frequently found in bacterial genomes, often offering advantageous features to the host. Here we assess the contribution of prophages for the evolution of A. baumannii pathogenicity. We found prophages to be notably diverse and widely disseminated in A. baumannii genomes. Also remarkably, A. baumannii prophages encode for multiple putative virulence factors that may be implicated in the bacterium’s capacity to colonize host niches, evade the host immune system, subsist in unfavorable environments, and tolerate antibiotics. Overall our results point towards a significant contribution of prophages for the dissemination and evolution of pathogenicity in A. baumannii, and highlight their clinical relevance.

High Diversity of Myocyanophage in Various Aquatic Environments Revealed by High-Throughput Sequencing of Major Capsid Protein Gene With a New Set of Primers

Myocyanophages, a group of viruses infecting cyanobacteria, are abundant and play important roles in elemental cycling. Here we investigated the particle-associated viral communities retained on 0.2 μm filters and in sediment samples (representing ancient cyanophage communities) from four ocean and three lake locations, using high-throughput sequencing and a newly designed primer pair targeting a gene fragment (∼145-bp in length) encoding the cyanophage gp23 major capsid protein (MCP). Diverse viral communities were detected in all samples. The fragments of 142-, 145-, and 148-bp in length were most abundant in the amplicons, and most sequences (>92%) belonged to cyanophages. Additionally, different sequencing depths resulted in different diversity estimates of the viral community. Operational taxonomic units obtained from deep sequencing of the MCP gene covered the majority of those obtained from shallow sequencing, suggesting that deep sequencing exhibited a more complete picture of cyanophage community than shallow sequencing. Our results also revealed a wide geographic distribution of marine myocyanophages, i.e., higher dissimilarities of the myocyanophage communities corresponded with the larger distances between the sampling sites. Collectively, this study suggests that the newly designed primer pair can be effectively used to study the community and diversity of myocyanophage from different environments, and the high-throughput sequencing represents a good method to understand viral diversity.

Assessing the diversity of the g23 gene of T4-like bacteriophages from Lake Baikal with high-throughput sequencing

Based on second generation sequencing (MiSeq platform, Illumina), we determined the genetic diversity of T4-like bacteriophages of the family Myoviridae by analysing fragments of the major capsid protein gene g23 in the plankton of Lake Baikal. The sampling depth in our study was significantly higher than in those obtained by the Sanger method before. We obtained 33 701 sequences of the g23 gene fragments, 141 operational taxonomic units (OTUs) of which were identified. 86 OTUs (60.9%) had the closest relatives from lakes Bourget and Annecy, and 28 OTUs (19.8%) had the highest identity with the Baikal g23 clones, which had been previously identified in the northern and southern basins of the lake by the Sanger method. The remaining OTUs were similar to the clones from other ecosystems. We showed a high genetic diversity of T4-type bacteriophages and a genetic difference with the phage communities from other ecosystems.

Metagenomic Approaches to Assess Bacteriophages in Various Environmental Niches

Bacteriophages are ubiquitous and numerous parasites of bacteria and play a critical evolutionary role in virtually every ecosystem, yet our understanding of the extent of the diversity and role of phages remains inadequate for many ecological niches, particularly in cases in which the host is unculturable. During the past 15 years, the emergence of the field of viral metagenomics has drastically enhanced our ability to analyse the so-called viral ‘dark matter’ of the biosphere. Here, we review the evolution of viral metagenomic methodologies, as well as providing an overview of some of the most significant applications and findings in this field of research.

vConTACT: an iVirus tool to classify double-stranded DNA viruses that infect Archaea and Bacteria

Taxonomic classification of archaeal and bacterial viruses is challenging, yet also fundamental for developing a predictive understanding of microbial ecosystems. Recent identification of hundreds of thousands of new viral genomes and genome fragments, whose hosts remain unknown, requires a paradigm shift away from traditional classification approaches and towards the use of genomes for taxonomy. Here we revisited the use of genomes and their protein content as a means for developing a viral taxonomy for bacterial and archaeal viruses. A network-based analytic was evaluated and benchmarked against authority-accepted taxonomic assignments and found to be largely concordant. Exceptions were manually examined and found to represent areas of viral genome 'sequence space' that are under-sampled or prone to excessive genetic exchange. While both cases are poorly resolved by genome-based taxonomic approaches, the former will improve as viral sequence space is better sampled and the latter are uncommon. Finally, given the largely robust taxonomic capabilities of this approach, we sought to enable researchers to easily and systematically classify new viruses. Thus, we established a tool, vConTACT, as an app at iVirus, where it operates as a fast, highly scalable, user-friendly app within the free and powerful CyVerse cyberinfrastructure.

Molecular Analysis of Arthrobacter Myovirus vB_ArtM-ArV1: We Blame It on the Tail

This is the first report on a myophage that infects Arthrobacter A novel virus, vB_ArtM-ArV1 (ArV1), was isolated from soil using Arthrobacter sp. strain 68b for phage propagation. Transmission electron microscopy showed its resemblance to members of the family Myoviridae: ArV1 has an isometric head (∼74 nm in diameter) and a contractile, nonflexible tail (∼192 nm). Phylogenetic and comparative sequence analyses, however, revealed that ArV1 has more genes in common with phages from the family Siphoviridae than it does with any myovirus characterized to date. The genome of ArV1 is a linear, circularly permuted, double-stranded DNA molecule (71,200 bp) with a GC content of 61.6%. The genome includes 101 open reading frames (ORFs) yet contains no tRNA genes. More than 50% of ArV1 genes encode unique proteins that either have no reliable identity to database entries or have homologues only in Arthrobacter phages, both sipho- and myoviruses. Using bioinformatics approaches, 13 ArV1 structural genes were identified, including those coding for head, tail, tail fiber, and baseplate proteins. A further 6 ArV1 ORFs were annotated as encoding putative structural proteins based on the results of proteomic analysis. Phylogenetic analysis based on the alignment of four conserved virion proteins revealed that Arthrobacter myophages form a discrete clade that seems to occupy a position somewhat intermediate between myo- and siphoviruses. Thus, the data presented here will help to advance our understanding of genetic diversity and evolution of phages that constitute the order CaudoviralesIMPORTANCE Bacteriophages, which likely originated in the early Precambrian Era, represent the most numerous population on the planet. Approximately 95% of known phages are tailed viruses that comprise three families: Podoviridae (with short tails), Siphoviridae (with long noncontractile tails), and Myoviridae (with contractile tails). Based on the current hypothesis, myophages, which may have evolved from siphophages, are thought to have first emerged among Gram-negative bacteria, whereas they emerged only later among Gram-positive bacteria. The results of the molecular characterization of myophage vB_ArtM-ArV1 presented here conform to the aforementioned hypothesis, since, at a glance, bacteriophage vB_ArtM-ArV1 appears to be a siphovirus that possesses a seemingly functional contractile tail. Our work demonstrates that such "chimeric" myophages are of cosmopolitan nature and are likely characteristic of the ecologically important soil bacterial genus Arthrobacter.

Things Are Getting Hairy: Enterobacteria Bacteriophage vB_PcaM_CBB

Enterobacteria phage vB_PcaM_CBB is a "jumbo" phage belonging to the family Myoviridae. It possesses highly atypical whisker-like structures along the length of its contractile tail. It has a broad host range with the capability of infecting species of the genera Erwinia, Pectobacterium, and Cronobacter. With a genome of 355,922 bp, excluding a predicted terminal repeat of 22,456 bp, phage CBB is the third largest phage sequenced to date. Its genome was predicted to encode 554 ORFs with 33 tRNAs. Based on prediction and proteome analysis of the virions, 29% of its predicted ORFs could be functionally assigned. Protein comparison shows that CBB shares between 33-38% of its proteins with Cronobacter phage GAP32, coliphages PBECO4 and 121Q as well as Klebsiella phage vB_KleM_Rak2. This work presents a detailed and comparative analysis of vB_PcaM_CBB of a highly atypical jumbo myoviridae phage, contributing to a better understanding of phage diversity and biology.

Characterization of LysPBC4, a novel Bacillus cereus-specific endolysin of bacteriophage PBC4

Bacillus cereus is a spore-forming, Gram-positive bacterium and is a major food-borne pathogen. A B. cereus-specific bacteriophage PBC4 was isolated from the soil of a stock farm, and its genome was analyzed. PBC4 belongs to the Siphoviridae family and has a genome consisting of 80 647-bp-long double-stranded DNA, including 123 genes and two tRNAs. LysPBC4, the endolysin of PBC4, has an enzymatically active domain (EAD) on its N-terminal region and a putative cell wall-binding domain (CBD) on its C-terminal region, respectively. Although the phage PBC4 showed a very limited host range, LysPBC4 could lyse all of the B. cereus strains tested. However, LysPBC4 did not kill other bacteria such as B. subtilis or Listeria, indicating that the endolysin has specific lytic activity against the B. cereus group species. Furthermore, LysPBC4_CBD fused with enhanced green fluorescent protein (EGFP) could decorate limited strains of B. cereus group, suggesting that the LysPBC4_CBD may be a promising material for specific detection of B. cereus.

Composition and Interactions among Bacterial, Microeukaryotic, and T4-like Viral Assemblages in Lakes from Both Polar Zones

In this study we assess global biogeography and correlation patterns among three components of microbial life: bacteria, microeukaryotes, and T4-like myoviruses. In addition to environmental and biogeographical considerations, we have focused our study on samples from high-latitude pristine lakes from both poles, since these simple island-like ecosystems represent ideal ecological models to probe the relationships among microbial components and with the environment. Bacterial assemblages were dominated by members of the same groups found to dominate freshwater ecosystems elsewhere, and microeukaryotic assemblages were dominated by photosynthetic microalgae. Despite inter-lake variations in community composition, the overall percentages of OTUs shared among sites was remarkable, indicating that many microeukaryotic, bacterial, and viral OTUs are globally-distributed. We observed an intriguing negative correlation between bacterial and microeukaryotic diversity values. Notably, our analyses show significant global correlations between bacterial and microeukaryotic community structures, and between the phylogenetic compositions of bacterial and T4-like virus assemblages. Overall, environmental filtering emerged as the main factor driving community structures.

Biogeography of bacteriophages at four hydrothermal vent sites in the Antarctic based on g23 sequence diversity

In this study, which was carried out within the ChEsSO consortium project (Chemosynthetically driven ecosystems south of the Polar Front), we sampled two hydrothermal vent sites on the East Scotia Ridge, Scotia Sea, one in the Kemp Caldera, South Sandwich Arc and one in the Bransfield Strait, north-west of the Antarctic Peninsula, which exhibit strong differences in their chemical characteristics. We compared a subset of their bacteriophage population by Sanger- and 454-sequencing of g23, which codes for the major capsid protein of T4likeviruses. We found that the sites differ vastly in their bacteriophage diversity, which reflects the differences in the chemical conditions and therefore putatively the differences in microbial hosts living at these sites. Comparing phage diversity in the vent samples to other aquatic samples, the vent samples formed a distinct separate cluster, which also included the non-vent control samples that were taken several hundred meters above the vent chimneys. This indicates that the influence of the vents on the microbial population and therefore also the bacteriophage population extends much further than anticipated.

The complete genome, structural proteome, comparative genomics and phylogenetic analysis of a broad host lytic bacteriophage ϕD3 infecting pectinolytic Dickeya spp

Plant necrotrophic Dickeya spp. are among the top ten most devastating bacterial plant pathogens able to infect a number of different plant species worldwide including economically important crops. Little is known of the lytic bacteriophages infecting Dickeya spp. A broad host lytic bacteriophage ϕD3 belonging to the family Myoviridae and order Caudovirales has been isolated in our previous study. This report provides detailed information of its annotated genome, structural proteome and phylogenetic relationships with known lytic bacteriophages infecting species of the Enterobacteriaceae family.

Using signature genes as tools to assess environmental viral ecology and diversity

Viruses (including bacteriophages) are the most abundant biological entities on the planet. As such, they are thought to have a major impact on all aspects of microbial community structure and function. Despite this critical role in ecosystem processes, the study of virus/phage diversity has lagged far behind parallel studies of the bacterial and eukaryotic kingdoms, largely due to the absence of any universal phylogenetic marker. Here we review the development and use of signature genes to investigate viral diversity, as a viable strategy for data sets of specific virus groups. Genes that have been used include those encoding structural proteins, such as portal protein, major capsid protein, and tail sheath protein, auxiliary metabolism genes, such as psbA, psbB,and phoH, and several polymerase genes. These marker genes have been used in combination with PCR-based fingerprinting and/or sequencing strategies to investigate spatial, temporal, and seasonal variations and diversity in a wide range of habitats.

Genomic characteristics and environmental distributions of the uncultivated Far-T4 phages

Viral metagenomics (viromics) is a tremendous tool to reveal viral taxonomic and functional diversity across ecosystems ranging from the human gut to the world's oceans. As with microbes however, there appear vast swaths of “dark matter” yet to be documented for viruses, even among relatively well-studied viral types. Here, we use viromics to explore the “Far-T4 phages” sequence space, a neighbor clade from the well-studied T4-like phages that was first detected through PCR study in seawater and subsequently identified in freshwater lakes through 454-sequenced viromes. To advance the description of these viruses beyond this single marker gene, we explore Far-T4 genome fragments assembled from two deeply-sequenced freshwater viromes. Single gene phylogenetic trees confirm that the Far-T4 phages are divergent from the T4-like phages, genome fragments reveal largely collinear genome organizations, and both data led to the delineation of five Far-T4 clades. Three-dimensional models of major capsid proteins are consistent with a T4-like structure, and highlight a highly conserved core flanked by variable insertions. Finally, we contextualize these now better characterized Far-T4 phages by re-analyzing 196 previously published viromes. These suggest that Far-T4 are common in freshwater and seawater as only four of 82 aquatic viromes lacked Far-T4-like sequences. Variability in representation across the five newly identified clades suggests clade-specific niche differentiation may be occurring across the different biomes, though the underlying mechanism remains unidentified. While complete genome assembly from complex communities and the lack of host linkage information still bottleneck virus discovery through viromes, these findings exemplify the power of metagenomics approaches to assess the diversity, evolutionary history, and genomic characteristics of novel uncultivated phages.

Composite conserved promoter-terminator motifs (PeSLs) that mediate modular shuffling in the diverse T4-like myoviruses

The diverse T4-like phages (Tquatrovirinae) infect a wide array of gram-negative bacterial hosts. The genome architecture of these phages is generally well conserved, most of the phylogenetically variable genes being grouped together in a series hyperplastic regions (HPRs) that are interspersed among large blocks of conserved core genes. Recent evidence from a pair of closely related T4-like phages has suggested that small, composite terminator/promoter sequences (promoterearly stem loop [PeSLs]) were implicated in mediating the high levels of genetic plasticity by indels occurring within the HPRs. Here, we present the genome sequence analysis of two T4-like phages, PST (168 kb, 272 open reading frames [ORFs]) and nt-1 (248 kb, 405 ORFs). These two phages were chosen for comparative sequence analysis because, although they are closely related to phages that have been previously sequenced (T4 and KVP40, respectively), they have different host ranges. In each case, one member of the pair infects a bacterial strain that is a human pathogen, whereas the other phage's host is a nonpathogen. Despite belonging to phylogenetically distant branches of the T4-likes, these pairs of phage have diverged from each other in part by a mechanism apparently involving PeSL-mediated recombination. This analysis confirms a role of PeSL sequences in the generation of genomic diversity by serving as a point of genetic exchange between otherwise unrelated sequences within the HPRs. Finally, the palette of divergent genes swapped by PeSL-mediated homologous recombination is discussed in the context of the PeSLs' potentially important role in facilitating phage adaption to new hosts and environments.

Characterization of a novel Pseudomonas aeruginosa bacteriophage, KPP25, of the family Podoviridae

Pseudomonas aeruginosa phages belonging to the family Podoviridae are one of the well-characterized phage groups. In this study, a novel P. aeruginosa phage, KPP25, was isolated and characterized. Phage KPP25's morphology was indicative of the family Podoviridae; however, analyses of the whole genome and the virion proteins suggested that it did not belong to any of the known podophage genera. Based on these analyses, phage KPP25 appears to be a novel podophage infecting P. aeruginosa.

The genome, proteome and phylogenetic analysis of Sinorhizobium meliloti phage ΦM12, the founder of a new group of T4-superfamily phages

Phage ΦM12 is an important transducing phage of the nitrogen-fixing rhizobial bacterium Sinorhizobium meliloti. Here we report the genome, phylogenetic analysis, and proteome of ΦM12, the first report of the genome and proteome of a rhizobium-infecting T4-superfamily phage. The structural genes of ΦM12 are most similar to T4-superfamily phages of cyanobacteria. ΦM12 is the first reported T4-superfamily phage to lack genes encoding class I ribonucleotide reductase (RNR) and exonuclease dexA, and to possess a class II coenzyme B12-dependent RNR. ΦM12's novel collection of genes establishes it as the founder of a new group of T4-superfamily phages, fusing features of cyanophages and phages of enteric bacteria.

Isolation, growth and genome of the Rhodothermus RM378 thermophilic bacteriophage

Several bacteriophages that infect different strains of the thermophilic bacterium Rhodothermus marinus were isolated and their infection pattern was studied. One phage, named RM378 was cultivated and characterized. The RM378 genome was also sequenced and analyzed. The phage was grouped as a member of the Myoviridae family with A2 morphology. It had a moderately elongated head, with dimensions of 85 and 95 nm between opposite apices and a 150 nm long tail, attached with a connector to the head. RM378 showed a virulent behavior that followed a lytic cycle of infection. It routinely gave lysates with 10(11) pfu/ml, and sometimes reached titers as high as 10(13) pfu/ml. The titer remained stable up to 65 °C but the phage lost viability when incubated at higher temperatures. Heating for 30 min at 96 °C lowered the titer by 10(4). The RM378 genome consisted of ds DNA of 129.908 bp with a GC ratio of 42.0% and contained about 120 ORFs. A few structural proteins, such as the major head protein corresponding to the gp23 in T4, could be identified. Only 29 gene products as probable homologs to other proteins of known function could be predicted, with most showing only low similarity to known proteins in other bacteriophages. These and other studies based on sequence analysis of a large number of phage genomes showed RM378 to be distantly related to all other known T4-like phages.

Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists

Characterizing ecological relationships between viruses, bacteria and protists in the ocean are critical to understanding ecosystem function, yet these relationships are infrequently investigated together. We evaluated these relationships through microbial association network analysis of samples collected approximately monthly from March 2008 to January 2011 in the surface ocean (0-5 m) at the San Pedro Ocean Time series station. Bacterial, T4-like myoviral and protistan communities were described by Automated Ribosomal Intergenic Spacer Analysis and terminal restriction fragment length polymorphism of the gene encoding the major capsid protein (g23) and 18S ribosomal DNA, respectively. Concurrent shifts in community structure suggested similar timing of responses to environmental and biological parameters. We linked T4-like myoviral, bacterial and protistan operational taxonomic units by local similarity correlations, which were then visualized as association networks. Network links (correlations) potentially represent synergistic and antagonistic relationships such as viral lysis, grazing, competition or other interactions. We found that virus-bacteria relationships were more cross-linked than protist-bacteria relationships, suggestive of increased taxonomic specificity in virus-bacteria relationships. We also found that 80% of bacterial-protist and 74% of bacterial-viral correlations were positive, with the latter suggesting that at monthly and seasonal timescales, viruses may be following their hosts more often than controlling host abundance.

Strong seasonality and interannual recurrence in marine myovirus communities

The temporal community dynamics and persistence of different viral types in the marine environment are still mostly obscure. Polymorphism of the major capsid protein gene, g23, was used to investigate the community composition dynamics of T4-like myoviruses in a North Atlantic fjord for a period of 2 years. A total of 160 unique operational taxonomic units (OTUs) were identified by terminal restriction fragment length polymorphism (TRFLP) of the gene g23. Three major community profiles were identified (winter-spring, summer, and autumn), which resulted in a clear seasonal succession pattern. These seasonal transitions were recurrent over the 2 years and significantly correlated with progression of seawater temperature, Synechococcus abundance, and turbidity. The appearance of the autumn viral communities was concomitant with the occurrence of prominent Synechococcus blooms. As a whole, we found a highly dynamic T4-like viral community with strong seasonality and recurrence patterns. These communities were unexpectedly dominated by a group of persistently abundant viruses.

Phage cluster relationships identified through single gene analysis

Background

Phylogenetic comparison of bacteriophages requires whole genome approaches such as dotplot analysis, genome pairwise maps, and gene content analysis. Currently mycobacteriophages, a highly studied phage group, are categorized into related clusters based on the comparative analysis of whole genome sequences. With the recent explosion of phage isolation, a simple method for phage cluster prediction would facilitate analysis of crude or complex samples without whole genome isolation and sequencing. The hypothesis of this study was that mycobacteriophage-cluster prediction is possible using comparison of a single, ubiquitous, semi-conserved gene. Tape Measure Protein (TMP) was selected to test the hypothesis because it is typically the longest gene in mycobacteriophage genomes and because regions within the TMP gene are conserved.

Results

A single gene, TMP, identified the known Mycobacteriophage clusters and subclusters using a Gepard dotplot comparison or a phylogenetic tree constructed from global alignment and maximum likelihood comparisons. Gepard analysis of 247 mycobacteriophage TMP sequences appropriately recovered 98.8% of the subcluster assignments that were made by whole-genome comparison. Subcluster-specific primers within TMP allow for PCR determination of the mycobacteriophage subcluster from DNA samples. Using the single-gene comparison approach for siphovirus coliphages, phage groupings by TMP comparison reflected relationships observed in a whole genome dotplot comparison and confirm the potential utility of this approach to another widely studied group of phages.

Conclusions

TMP sequence comparison and PCR results support the hypothesis that a single gene can be used for distinguishing phage cluster and subcluster assignments. TMP single-gene analysis can quickly and accurately aid in mycobacteriophage classification.

Diversity of the major capsid genes (g23) of T4-like bacteriophages in the eutrophic Lake Kotokel in East Siberia, Russia

Numerous studies revealed high diversity of T4-like bacteriophages in various environments, but so far, little is known about T4-like virus diversity in freshwater bodies, particularly in eutrophic lakes. The present study was aimed at elucidating molecular diversity of T4-like bacteriophages in eutrophic Lake Kotokel located near Lake Baikal by partial sequencing of the major capsid genes (g23) of T4-like bacteriophages. The majority of g23 fragments from Lake Kotokel were most similar to those from freshwater lakes and paddy fields. Despite the proximity and direct water connection between Lake Kotokel and Lake Baikal, g23 sequence assemblages from two lakes were different. UniFrac analysis showed that uncultured T4-like viruses from Lake Kotokel tended to cluster with those from the distant lake of the same trophic status. This fact suggested that the trophic conditions affected the formation of viral populations, particularly of T4-like viruses, in freshwater environments.

Short-term observations of marine bacterial and viral communities: patterns, connections and resilience

Observation of short-term temporal variation in bacterial and viral communities is important for understanding patterns of aquatic microbial diversity. We collected surface seawater once daily for 38 consecutive days with seven more samples interspersed over 40 more days at one location ∼2 km from Santa Catalina Island, California. Bacterial communities were analyzed by automated ribosomal intergenic spacer analysis (ARISA) and viral communities were analyzed by terminal restriction fragment length polymorphism (TRFLP) of the conserved T4-like myoviral gene encoding the major capsid protein (g23). Common bacterial and viral taxa were consistently dominant, and relatively few displayed dramatic increases/decreases or 'boom/bust' patterns that might be expected from dynamic predator-prey interactions. Association network analysis showed most significant covariations (associations) occurred among bacterial taxa or among viral taxa and there were several modular (highly-interconnected) associations (P≤0.005). Associations observed between bacteria and viruses (P≤0.005) occurred with a median time lag of 2 days. Regression of all pairwise Bray-Curtis similarities between samples indicated a rate of bacterial community change that slows from 2.1%-0.18% per day over a week to 2 months; the rate stays around 0.4% per day for viruses. Our interpretation is that, over the scale of days, individual bacterial and viral OTUs can be dynamic and patterned; resulting in statistical associations regarded as potential ecological interactions. However, over the scale of weeks, average bacterial community variation is slower, suggesting that there is strong community-level ecological resilience, that is, a tendency to converge towards a 'mean' microbial community set by longer-term controlling factors.

Pseudomonas aeruginosa keratitis in mice: effects of topical bacteriophage KPP12 administration

The therapeutic effects of bacteriophage (phage) KPP12 in Pseudomonas aeruginosa keratitis were investigated in mice. Morphological analysis showed that phage KPP12 is a member of the family Myoviridae, morphotype A1, and DNA sequence analysis revealed that phage KPP12 is similar to PB1-like viruses. Analysis of the phage KPP12 genome did not identify any genes related to drug resistance, pathogenicity or lysogenicity, and so phage KPP12 may be a good candidate for therapeutic. KPP12 showed a broad host range for P. aeruginosa strains isolated from clinical ophthalmic infections. Inoculation of the scarified cornea with P. aeruginosa caused severe keratitis and eventual corneal perforation. Subsequent single-dose administration of KPP12 eye-drops significantly improved disease outcome, and preserved the structural integrity and transparency of the infected cornea. KPP12 treatment resulted in the suppression of neutrophil infiltration and greatly enhanced bacterial clearance in the infected cornea. These results indicate that bacteriophage eye-drops may be a novel adjunctive or alternative therapeutic agent for the treatment of infectious keratitis secondary to antibiotic-resistant bacteria.

Ocean viruses: rigorously evaluating the metagenomic sample-to-sequence pipeline

As new environments are studied, viruses consistently emerge as important and prominent players in natural and man-made ecosystems. However, much of what we know is built both upon the foundation of the culturable minority and using methods that are often insufficiently ground-truthed. Here, we review the modern culture-independent viral metagenomic sample-to-sequence pipeline and how next-generation sequencing techniques are drastically altering our ability to systematically and rigorously evaluate them. Together, a series of studies quantitatively evaluate existing and new methods that allow-even for ultra-low DNA samples-the generation of replicable, near-quantitative datasets that maximize inter-comparability and biological inference.

Evaluation of two approaches for assessing the genetic similarity of virioplankton populations as defined by genome size

Viral production estimates show that virioplankton communities turn over rapidly in aquatic ecosystems. Thus, it is likely that the genetic identity of viral populations comprising the virioplankton also change over temporal and spatial scales, reflecting shifts in viral-host interactions. However, there are few approaches that can provide data on the genotypic identity of viral populations at low cost and with the sample throughput necessary to assess dynamic changes in the virioplankton. This study examined two of these approaches-T4-like major capsid protein (g23) gene polymorphism and randomly amplified polymorphic DNA-PCR (RAPD-PCR) fingerprinting-to ask how well each technique could track differences in virioplankton populations over time and geographic location. Seasonal changes in overall virioplankton composition were apparent from pulsed-field gel electrophoresis (PFGE) analysis. T4-like phages containing similar g23 proteins were found within both small- and large-genome populations, including populations from different geographic locations and times. The surprising occurrence of T4-like g23 within small genomic groups (23 to 64 kb) indicated that the genome size range of T4-like phages may be broader than previously believed. In contrast, RAPD-PCR fingerprinting detected high genotypic similarity within PFGE bands from the same location, time, and genome size class without the requirement for DNA sequencing. Unlike g23 polymorphism, RAPD-PCR fingerprints showed a greater temporal than geographic variation. Thus, while polymorphism in a viral signature gene, such as g23, can be a powerful tool for inferring evolutionary relationships, the degree to which this approach can capture fine-scale variability within virioplankton populations is less clear.

Seasonality and monthly dynamics of marine myovirus communities

Marine myoviruses (i.e. bacteriophages with a contractile tail sheath) are numerically abundant and genetically diverse. We developed a terminal restriction fragment length polymorphism assay (TRFLP) for g23, the conserved gene encoding the major capsid protein, to investigate T4-like myovirus communities at USC's Microbial Observatory at the San Pedro Ocean Time-series (SPOT), where we previously reported bacterial seasonality. Between 71 and 154 operational taxonomic units (OTUs) were observed monthly over 3 years. Roughly 25% of OTUs were detected in 31 or more months. T4-like myoviral community structure varied seasonally with some OTUs peaking repeatedly in spring-summer and others in fall-winter, while moderately abundant OTUs persisted year-round. Recurring community structure was demonstrated using discriminant function analysis (DFA, selecting taxa that best predict months) and average Bray-Curtis similarity. DFA showed communities from adjacent months or 12 months apart were positively auto-correlated, while communities 3-7 months apart were negatively auto-correlated. Bray-Curtis similarity was highest between adjacent months - with a local maximum at 12-month and local minima at 6- and 18- to 20-month lags. The T4-like virus community at SPOT exhibited seasonality, yet the somewhat unexpected persistence of moderately abundant OTUs and predictability of the community add new twists to existing conceptual models of marine viruses.

Assessing the diversity and specificity of two freshwater viral communities through metagenomics

Transitions between saline and fresh waters have been shown to be infrequent for microorganisms. Based on host-specific interactions, the presence of specific clades among hosts suggests the existence of freshwater-specific viral clades. Yet, little is known about the composition and diversity of the temperate freshwater viral communities, and even if freshwater lakes and marine waters harbor distinct clades for particular viral sub-families, this distinction remains to be demonstrated on a community scale.

To help identify the characteristics and potential specificities of freshwater viral communities, such communities from two lakes differing by their ecological parameters were studied through metagenomics. Both the cluster richness and the species richness of the Lake Bourget virome were significantly higher that those of the Lake Pavin, highlighting a trend similar to the one observed for microorganisms (i.e. the specie richness observed in mesotrophic lakes is greater than the one observed in oligotrophic lakes). Using 29 previously published viromes, the cluster richness was shown to vary between different environment types and appeared significantly higher in marine ecosystems than in other biomes. Furthermore, significant genetic similarity between viral communities of related environments was highlighted as freshwater, marine and hypersaline environments were separated from each other despite the vast geographical distances between sample locations within each of these biomes. An automated phylogeny procedure was then applied to marker genes of the major families of single-stranded (Microviridae, Circoviridae, Nanoviridae) and double-stranded (Caudovirales) DNA viruses. These phylogenetic analyses all spotlighted a very broad diversity and previously unknown clades undetectable by PCR analysis, clades that gathered sequences from the two lakes. Thus, the two freshwater viromes appear closely related, despite the significant ecological differences between the two lakes. Furthermore, freshwater viral communities appear genetically distinct from other aquatic ecosystems, demonstrating the specificity of freshwater viruses at a community scale for the first time.

A novel cyanophage with a cyanobacterial nonbleaching protein A gene in the genome

A cyanophage, PaV-LD, has been isolated from harmful filamentous cyanobacterium Planktothrix agardhii in Lake Donghu, a shallow freshwater lake in China. Here, we present the cyanophage's genomic organization and major structural proteins. The genome is a 95,299-bp-long, linear double-stranded DNA and contains 142 potential genes. BLAST searches revealed 29 proteins of known function in cyanophages, cyanobacteria, or bacteria. Thirteen major structural proteins ranging in size from 27 kDa to 172 kDa were identified by SDS-PAGE and mass-spectrometric analysis. The genome lacks major genes that are necessary to the tail structure, and the tailless PaV-LD has been confirmed by an electron microscopy comparison with other tail cyanophages and phages. Phylogenetic analysis of the major capsid proteins also reveals an independent branch of PaV-LD that is quite different from other known tail cyanophages and phages. Moreover, the unique genome carries a nonbleaching protein A (NblA) gene (open reading frame [ORF] 022L), which is present in all phycobilisome-containing organisms and mediates phycobilisome degradation. Western blot detection confirmed that 022L was expressed after PaV-LD infection in the host filamentous cyanobacterium. In addition, its appearance was companied by a significant decline of phycocyanobilin content and a color change of the cyanobacterial cells from blue-green to yellow-green. The biological function of PaV-LD nblA was further confirmed by expression in a model cyanobacterium via an integration platform, by spectroscopic analysis and electron microscopy observation. The data indicate that PaV-LD is an exceptional cyanophage of filamentous cyanobacteria, and this novel cyanophage will also provide us with a new vision of the cyanophage-host interactions.

Genome of low-temperature T4-related bacteriophage vB_EcoM-VR7

The complete genome sequence of the T4-related low-temperature Escherichia coli bacteriophage vB_EcoM-VR7 was determined. The genome sequence is 169,285 bp long, with a G+C content of 40.3%. Overall, 95% of the phage genome is coding. It encodes 293 putative protein-encoding open reading frames (ORFs) and tRNA(Met). More than half (59%) of the genomic DNA lacks significant homology with the DNA of T4, but once translated, 72% of the vB_EcoM-VR7 genome (211 ORFs) encodes protein homologues of T4 genes. Overall, 46 vB_EcoM-VR7 ORFs have no homologues in T4 but are derived from other T4-related phages, nine ORFs show similarities to bacterial or non-T4-related phage genes, and 27 ORFs are unique to vB_EcoM-VR7. This phage lacks several T4 enzymes involved in host DNA degradation; however, there is extensive representation of the DNA replication, recombination and repair enzymes as well as the viral capsid and tail structural genes.

The gp38 adhesins of the T4 superfamily: a complex modular determinant of the phage's host specificity

The tail fiber adhesins are the primary determinants of host range in the T4-type bacteriophages. Among the indispensable virion components, the sequences of the long tail fiber genes and their associated adhesins are among the most variable. The predominant form of the adhesin in the T4-type phages is not even the version of the gene encoded by T4, the archetype of the superfamily, but rather a small unrelated protein (gp38) encoded by closely related phages such as T2 and T6. This gp38 adhesin has a modular design: its N-terminal attachment domain binds at the tip of the tail fiber, whereas the C-terminal specificity domain determines its host receptor affinity. This specificity domain has a series of four hypervariable segments (HVSs) that are separated by a set of highly conserved glycine-rich motifs (GRMs) that apparently form the domain’s conserved structural core. The role of gp38’s various components was examined by a comparative analysis of a large series of gp38 adhesins from T-even superfamily phages with differing host specificities. A deletion analysis revealed that the individual HVSs and GRMs are essential to the T6 adhesin’s function and suggests that these different components all act in synergy to mediate adsorption. The evolutionary advantages of the modular design of the adhesin involving both conserved structural elements and multiple independent and easily interchanged specificity determinants are discussed.