Analysis of high-throughput sequencing and annotation strategies for phage genomes

Droplet-based methodology for investigating bacterial population dynamics in response to phage exposure

An alarming rise in antimicrobial resistance worldwide has spurred efforts into the search for alternatives to antibiotic treatments. The use of bacteriophages, bacterial viruses harmless to humans, represents a promising approach with potential to treat bacterial infections (phage therapy). Recent advances in microscopy-based single-cell techniques have allowed researchers to develop new quantitative methodologies for assessing the interactions between bacteria and phages, especially the ability of phages to eradicate bacterial pathogen populations and to modulate growth of both commensal and pathogen populations. Here we combine droplet microfluidics with fluorescence time-lapse microscopy to characterize the growth and lysis dynamics of the bacterium Escherichia coli confined in droplets when challenged with phage. We investigated phages that promote lysis of infected E. coli cells, specifically, a phage species with DNA genome, T7 (Escherichia virus T7) and two phage species with RNA genomes, MS2 (Emesvirus zinderi) and Qβ (Qubevirus durum). Our microfluidic trapping device generated and immobilized picoliter-sized droplets, enabling stable imaging of bacterial growth and lysis in a temperature-controlled setup. Temporal information on bacterial population size was recorded for up to 25 h, allowing us to determine growth rates of bacterial populations and helping us uncover the extent and speed of phage infection. In the long-term, the development of novel microfluidic single-cell and population-level approaches will expedite research towards fundamental understanding of the genetic and molecular basis of rapid phage-induced lysis and eco-evolutionary aspects of bacteria-phage dynamics, and ultimately help identify key factors influencing the success of phage therapy.

Renaissance for Phage-Based Bacterial Control

Bacteriophages (phages) are an underutilized biological resource with vast potential for pathogen control and microbiome editing. Phage research and commercialization have increased rapidly in biomedical and agricultural industries, but adoption has been limited elsewhere. Nevertheless, converging advances in DNA sequencing, bioinformatics, microbial ecology, and synthetic biology are now poised to broaden phage applications beyond pathogen control toward the manipulation of microbial communities for defined functional improvements. Enhancements in sequencing combined with network analysis make it now feasible to identify and disrupt microbial associations to elicit desirable shifts in community structure or function, indirectly modulate species abundance, and target hub or keystone species to achieve broad functional shifts. Sequencing and bioinformatic advancements are also facilitating the use of temperate phages for safe gene delivery applications. Finally, integration of synthetic biology stands to create novel phage chassis and modular genetic components. While some fundamental, regulatory, and commercialization barriers to widespread phage use remain, many major challenges that have impeded the field now have workable solutions. Thus, a new dawn for phage-based (chemical-free) precise biocontrol and microbiome editing is on the horizon to enhance, suppress, or modulate microbial activities important for public health, food security, and more sustainable energy production and water reuse.

Genetic engineering of marine cyanophages reveals integration but not lysogeny in T7-like cyanophages

Marine cyanobacteria of the genera Synechococcus and Prochlorococcus are the most abundant photosynthetic organisms on earth, spanning vast regions of the oceans and contributing significantly to global primary production. Their viruses (cyanophages) greatly influence cyanobacterial ecology and evolution. Although many cyanophage genomes have been sequenced, insight into the functional role of cyanophage genes is limited by the lack of a cyanophage genetic engineering system. Here, we describe a simple, generalizable method for genetic engineering of cyanophages from multiple families, that we named REEP for REcombination, Enrichment and PCR screening. This method enables direct investigation of key cyanophage genes, and its simplicity makes it adaptable to other ecologically relevant host-virus systems. T7-like cyanophages often carry integrase genes and attachment sites, yet exhibit lytic infection dynamics. Here, using REEP, we investigated their ability to integrate and maintain a lysogenic life cycle. We found that these cyanophages integrate into the host genome and that the integrase and attachment site are required for integration. However, stable lysogens did not form. The frequency of integration was found to be low in both lab cultures and the oceans. These findings suggest that T7-like cyanophage integration is transient and is not part of a classical lysogenic cycle.

A rigorous assessment and comparison of enumeration methods for environmental viruses

Determining exact viral titers in a given sample is essential for many environmental and clinical applications, e.g., for studying viral ecology or application of bacteriophages for food safety. However, virus quantification is not a simple task, especially for complex environmental samples. While clonal viral isolates can be quantified with relative high accuracy using virus-specific methods, i.e., plaque assay or quantitative real-time PCR, these methods are not valid for complex and diverse environmental samples. Moreover, it is not yet known how precisely laser-based methods, i.e., epifluorescence microscopy, flow cytometry, and nanoparticle tracking analysis, quantify environmental viruses. In the present study, we compared five state-of-the-art viral quantification methods by enumerating four model viral isolates of different genome and size characteristics as well as four different environmental water samples. Although Nanoparticle tracking analysis combined with gentle staining at 30 °C could be confirmed by this study to be a reliable quantification technique for tested environmental samples, environmental samples still lack an universally applicable and accurate quantification method. Special attention has to be put on optimal sample concentrations as well as optimized sample preparations, which are specific for each method. As our results show the inefficiency when enumerating small, or single-stranded DNA or RNA viruses, the global population of viruses is presumably higher than expected.

Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome

RNA viruses in aquatic environments remain poorly studied. Here, we analysed the RNA virome from approximately 10 l water from Yangshan Deep-Water Harbour near the Yangtze River estuary in China and identified more than 4,500 distinct RNA viruses, doubling the previously known set of viruses. Phylogenomic analysis identified several major lineages, roughly, at the taxonomic ranks of class, order and family. The 719-member-strong Yangshan virus assemblage is the sister clade to the expansive class Alsuviricetes and consists of viruses with simple genomes that typically encode only RNA-dependent RNA polymerase (RdRP), capping enzyme and capsid protein. Several clades within the Yangshan assemblage independently evolved domain permutation in the RdRP. Another previously unknown clade shares ancestry with Potyviridae, the largest known plant virus family. The ‘Aquatic picorna-like viruses/Marnaviridae’ clade was greatly expanded, with more than 800 added viruses. Several RdRP-linked protein domains not previously detected in any RNA viruses were identified, such as the small ubiquitin-like modifier (SUMO) domain, phospholipase A2 and PrsW-family protease domain. Multiple viruses utilize alternative genetic codes implying protist (especially ciliate) hosts. The results reveal a vast RNA virome that includes many previously unknown groups. However, phylogenetic analysis of the RdRPs supports the previously established five-branch structure of the RNA virus evolutionary tree, with no additional phyla.

Metagenomic analysis of a single RNA virome from the Yangshan Deep-Water Harbour in China enabled the recovery of more than 4,500 distinct RNA viruses, doubling the known set of RNA viruses to date, and provided insights into their biology.

Evaluation of Sequencing Library Preparation Protocols for Viral Metagenomic Analysis from Pristine Aquifer Groundwaters

Viral ecology of terrestrial habitats is yet-to be extensively explored, in particular the terrestrial subsurface. One problem in obtaining viral sequences from groundwater aquifer samples is the relatively low amount of virus particles. As a result, the amount of extracted DNA may not be sufficient for direct sequencing of such samples. Here we compared three DNA amplification methods to enrich viral DNA from three pristine limestone aquifer assemblages of the Hainich Critical Zone Exploratory to evaluate potential bias created by the different amplification methods as determined by viral metagenomics. Linker amplification shotgun libraries resulted in lowest redundancy among the sequencing reads and showed the highest diversity, while multiple displacement amplification produced the highest number of contigs with the longest average contig size, suggesting a combination of these two methods is suitable for the successful enrichment of viral DNA from pristine groundwater samples. In total, we identified 27,173, 5,886 and 32,613 viral contigs from the three samples from which 11.92 to 18.65% could be assigned to taxonomy using blast. Among these, members of the Caudovirales order were the most abundant group (52.20 to 69.12%) dominated by Myoviridae and Siphoviridae. Those, and the high number of unknown viral sequences, substantially expand the known virosphere.

Broad Environmental Tolerance for a Salicola Host-Phage Pair Isolated from the Cargill Solar Saltworks, Newark, CA, USA

Phages greatly influence the ecology and evolution of their bacterial hosts; however, compared to hosts, a relatively low number of phages, especially halophilic phages, have been studied. This study describes a comparative investigation of physicochemical tolerance between a strain of the halophilic bacterium, Salicola, isolated from the Cargill Saltworks (Newark, CA, USA) and its associated phage. The host grew in media between pH 6–8.5, had a salinity growth optimum of 20% total salts (ranging from 10%–30%) and an upper temperature growth limit of 48 °C. The host utilized 61 of 190 substrates tested using BIOLOG Phenotype MicroArrays. The CGφ29 phage, one of only four reported Salicola phages, is a DNA virus of the Siphoviridae family. Overall, the phage tolerated a broader range of environmental conditions than its host (salinity 0–30% total salts; pH 3–9; upper thermal limit 80 °C) and is the most thermotolerant halophilic phage ever reported. This study is the most comprehensive investigation to date of a Salicola host–phage pair and provides novel insights into extreme environmental tolerances among bacteriophages.

Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses

BACKGROUND

Viruses are key players regulating microbial ecosystems. Exploration of viral assemblages is now possible thanks to the development of metagenomics, the most powerful tool available for studying viral ecology and discovering new viruses. Unfortunately, several sources of bias lead to the misrepresentation of certain viruses within metagenomics workflows, hindering the shift from merely descriptive studies towards quantitative comparisons of communities. Therefore, benchmark studies on virus enrichment and random amplification protocols are required to better understand the sources of bias.

RESULTS

We assessed the bias introduced by viral enrichment on mock assemblages composed of seven DNA viruses, and the bias from random amplification methods on human saliva DNA viromes, using qPCR and deep sequencing, respectively. While iodixanol cushions and 0.45 μm filtration preserved the original composition of nuclease-protected viral genomes, low-force centrifugation and 0.22 μm filtration removed large viruses. Comparison of unamplified and randomly amplified saliva viromes revealed that multiple displacement amplification (MDA) induced stochastic bias from picograms of DNA template. However, the type of bias shifted to systematic using 1 ng, with only a marginal influence by amplification time. Systematic bias consisted of over-amplification of small circular genomes, and under-amplification of those with extreme GC content, a negative bias that was shared with the PCR-based sequence-independent, single-primer amplification (SISPA) method. MDA based on random priming provided by a DNA primase activity slightly outperformed those based on random hexamers and SISPA, which may reflect differences in ability to handle sequences with extreme GC content. SISPA viromes showed uneven coverage profiles, with high coverage peaks in regions with low linguistic sequence complexity. Despite misrepresentation of certain viruses after random amplification, ordination plots based on dissimilarities among contig profiles showed perfect overlapping of related amplified and unamplified saliva viromes and strong separation from unrelated saliva viromes. This result suggests that random amplification bias has a minor impact on beta diversity studies.

CONCLUSIONS

Benchmark analyses of mock and natural communities of viruses improve understanding and mitigate bias in metagenomics surveys. Bias induced by random amplification methods has only a minor impact on beta diversity studies of human saliva viromes.

Complete Genome Sequence of Ebrios, a Novel T7virus Isolated from the Ebrie Lagoon in Abidjan, Côte d'Ivoire

The lytic Escherichia coli phage Ebrios was isolated from a water sample collected in Ebrie Lagoon on the Adiopodoumé River in Abidjan (Republic of Côte d’Ivoire, West Africa). The linear genome of this Podoviridae family member contains 39,752 bp, has a G+C content of 52.9%, is composed of 53 open reading frames, and is related to the Stenotrophomonas maltophilia phage IME15.

Phage Genetic Engineering Using CRISPR⁻Cas Systems

Since their discovery over a decade ago, the class of prokaryotic immune systems known as CRISPR⁻Cas have afforded a suite of genetic tools that have revolutionized research in model organisms spanning all domains of life. CRISPR-mediated tools have also emerged for the natural targets of CRISPR⁻Cas immunity, the viruses that specifically infect bacteria, or phages. Despite their status as the most abundant biological entities on the planet, the majority of phage genes have unassigned functions. This reality underscores the need for robust genetic tools to study them. Recent reports have demonstrated that CRISPR⁻Cas systems, specifically the three major types (I, II, and III), can be harnessed to genetically engineer phages that infect diverse hosts. Here, the mechanisms of each of these systems, specific strategies used, and phage editing efficacies will be reviewed. Due to the relatively wide distribution of CRISPR⁻Cas systems across bacteria and archaea, it is anticipated that these immune systems will provide generally applicable tools that will advance the mechanistic understanding of prokaryotic viruses and accelerate the development of novel technologies based on these ubiquitous organisms.

The Human Virome Protein Cluster Database (HVPC): A Human Viral Metagenomic Database for Diversity and Function Annotation

Human virome, including those of bacteria (bacteriophages) have received an increasing attention recently, owing to the rapid developments in human microbiome research and the awareness of the far-reaching influence of microbiomes on health and disease. Nevertheless, human viromes are still underrepresented in literature making viruses a virtually untapped resource of diversity, functional and physiological information. Here we present the human virome protein cluster database as an effort to improve functional annotation and characterization of human viromes. The database was built out of hundreds of virome datasets from six different body sites. We also show the utility of this database through its use for the characterization of three bronchoalveolar lavage (BAL) viromes from one healthy control in addition to one moderate and one severe chronic obstructive pulmonary disease (COPD) patients. The use of the database allowed for a better functional annotation, which were otherwise poorly characterized when limited to annotation using sequences from full-length viral genomes. In addition, our BAL samples gave a first insight into viral communities of COPD patients and confirm a state of dysbiosis for viruses that increases with disease progression. Moreover, they shed light on the potential role of phages in the horizontal gene transfer of bacterial virulence factors, a phenomenon that highlights a possible contribution of phages to etiopathology.

The Number of Target Molecules of the Amplification Step Limits Accuracy and Sensitivity in Ultradeep-Sequencing Viral Population Studies

The invention of next-generation sequencing (NGS) techniques marked the coming of a new era in the detection of the genetic diversity of intrahost viral populations. A good understanding of the genetic structure of these populations requires, first, the ability to identify the different isolates or variants and, second, the ability to accurately quantify them. However, the initial amplification step of NGS studies can impose potential quantitative biases, modifying the variant relative frequencies. In particular, the number of target molecules (NTM) used during the amplification step is vastly overlooked although of primary importance, as it sets the limit of the accuracy and sensitivity of the sequencing procedure. In the present article, we investigated quantitative biases in an NGS study of populations of a multipartite single-stranded DNA (ssDNA) virus at different steps of the procedure. We studied 20 independent populations of the ssDNA virus faba bean necrotic stunt virus (FBNSV) in two host plants, Vicia faba and Medicago truncatula FBNSV is a multipartite virus composed of eight genomic segments, whose specific and host-dependent relative frequencies are defined as the "genome formula." Our results show a significant distortion of the FBNSV genome formula after the amplification and sequencing steps. We also quantified the genetic bottleneck occurring at the amplification step by documenting the NTM of two genomic segments of FBNSV. We argue that the NTM must be documented and carefully considered when determining the sensitivity and accuracy of data from NGS studies.IMPORTANCE The advent of next-generation sequencing (NGS) techniques now enables study of the genetic diversity of viral populations. A good understanding of the genetic structure of these populations first requires the ability to identify the different isolates or variants and second requires the ability to accurately quantify them. Prior to sequencing, viral genomes need to be amplified, a step that potentially imposes quantitative biases and modifies the viral population structure. In particular, the number of target molecules (NTM) used during the amplification step is of primary importance, as it sets the limit of the accuracy and sensitivity of the sequencing procedure. In this work, we used 20 replicated populations of the multipartite faba bean necrotic stunt virus (FBNSV) to estimate the various limitations of ultradeep-sequencing studies performed on intrahost viral populations. We report quantitative biases during rolling-circle amplification and the NTM of two genomic segments of FBNSV.

Genome analysis of three novel lytic Vibrio coralliilyticus phages isolated from seawater, Okinawa, Japan

Three novel Vibrio phages were isolated from seawater in Okinawa. The Vibrio phage RYC infected Vibrio coralliilyticus SWA 07, while Vibrio phages CKB-S1 and CKB-S2 infected the coral pathogen V. coralliilyticus P1 (LMG 23696). The Vibrio phages CKB-S1 and CKB-S2 displayed head-tail structures whereas the Vibrio phage RYC showed a tailless non-enveloped capsid. All these Vibrio phages contained linear and double-stranded DNA. The whole genome sequencing revealed that Vibrio phage RYC has a larger genome size compared to Vibrio phages CKB-S1 and CKB-S2, and six tRNAs genes were found only in Vibrio phage RYC. Genome-wide comparison showed that Vibrio phage CKB-S1 was closely related, but was not identical, to Vibrio parahaemolyticus phages VP16T and VP16C. Meanwhile, the Vibrio phages RYC and CKB-S2 did not show high genome-wide similarity to any phages. These results suggest that the Vibrio phages CKB-S1, CKB-S2 and RYC are novel phages, which need further exploration, especially for their potential applications in phage therapy.

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.

The genomic content and context of auxiliary metabolic genes in marine cyanomyoviruses

Viruses of marine cyanobacteria frequently contain auxiliary metabolic genes (AMGs) that augment host metabolism during infection, but little is known about their adaptive significance. We analyzed the distribution and genomic context of 33 AMGs across 60 cyanomyovirus genomes. Similarity in AMG content among cyanomyoviruses was only weakly correlated with phylogenetic relatedness; however, AMG content was generally conserved within the same operational taxonomic unit (OTU). A virus' AMG repertoire was also correlated with its isolation host and environment (coastal versus open ocean). A new analytical method based on shared co-linear blocks revealed that variation in the genomic location of an AMG was negatively correlated with its frequency across the genomes. We propose that rare AMGs are more frequently gained or lost as a result of fluctuating selection pressures, whereas common AMGs are associated with stable selection pressures. Finally, we describe a unique cyanomyovirus (S-CAM7) that lacks many AMGs including the photosynthesis gene psbA.

Genome analysis and gene nblA identification of Microcystis aeruginosa myovirus (MaMV-DC) reveal the evidence for horizontal gene transfer events between cyanomyovirus and host

The genome sequence, genetic characterization and nblA gene function of Microcystis aeruginosa myovirus isolated from Lake Dianchi in China (MaMV-DC) have been analysed. The genome DNA is 169 223 bp long, with 170 predicted protein-coding genes (001L–170L) and a tRNA gene. About one-sixth of these genes have homologues in the host cyanobacteria M. aeruginosa. The genome carries a gene homologous to host nblA, which encodes a protein involved in the degradation of cyanobacterial phycobilisome. Its expression during MaMV-DC infection was confirmed by reverse transcriptase PCR and Western blot detection and abundant expression was companied by the significant decline of phycocyanin content and massive release of progeny MaMV-DC. In addition, expressing MaMV-DC nblA reduced the phycocyanin peak and the phycocyanin to chlorophyll ratio in model cyanobacteria. These results confirm that horizontal gene transfer events have occurred between cyanobacterial host and cyanomyovirus and suggest that MaMV-DC carrying host-derived genes (such as 005L, that codes for NblA) is responsible for more efficient expression of cyanophage genes and release of progeny cyanophage. This study provides novel insight into the horizontal gene transfer in cyanophage and the interactions between cyanophage and their host.

Assessing Illumina technology for the high-throughput sequencing of bacteriophage genomes

Bacteriophages are the most abundant biological entities on the planet, playing crucial roles in the shaping of bacterial populations. Phages have smaller genomes than their bacterial hosts, yet there are currently fewer fully sequenced phage than bacterial genomes. We assessed the suitability of Illumina technology for high-throughput sequencing and subsequent assembly of phage genomes. In silico datasets reveal that 30× coverage is sufficient to correctly assemble the complete genome of ~98.5% of known phages, with experimental data confirming that the majority of phage genomes can be assembled at 30× coverage. Furthermore, in silico data demonstrate it is possible to co-sequence multiple phages from different hosts, without introducing assembly errors.

Genetic characterization of ØVC8 lytic phage for Vibrio cholerae O1

Background

Epidemics and pandemics of cholera, a diarrheal disease, are attributed to Vibrio cholerae serogroups O1 and O139. In recent years, specific lytic phages of V. cholerae have been proposed to be important factors in the cyclic occurrence of cholera in endemic areas. However, the role and potential participation of lytic phages during long interepidemic periods of cholera in non-endemic regions have not yet been described. The purpose of this study was to isolate and characterize specific lytic phages of V. cholerae O1 strains.

Methods

Sixteen phages were isolated from wastewater samples collected at the Endhó Dam in Hidalgo State, Mexico, concentrated with PEG/NaCl, and purified by density gradient. The lytic activity of the purified phages was tested using different V. cholerae O1 and O139 strains. Phage morphology was visualized by transmission electron microscopy (TEM), and phage genome sequencing was performed using the Genome Analyzer IIx System. Genome assembly and bioinformatics analysis were performed using a set of high-throughput programs. Phage structural proteins were analyzed by mass spectrometry.

Results

Sixteen phages with lytic and lysogenic activity were isolated; only phage ØVC8 showed specific lytic activity against V. cholerae O1 strains. TEM images of ØVC8 revealed a phage with a short tail and an isometric head. The ØVC8 genome comprises linear double-stranded DNA of 39,422 bp with 50.8 % G + C. Of the 48 annotated ORFs, 16 exhibit homology with sequences of known function and several conserved domains. Bioinformatics analysis showed multiple conserved domains, including an Ig domain, suggesting that ØVC8 might adhere to different mucus substrates such as the human intestinal epithelium. The results suggest that ØVC8 genome utilize the “single-stranded cohesive ends” packaging strategy of the lambda-like group. The two structural proteins sequenced and analyzed are proteins of known function.

Conclusions

ØVC8 is a lytic phage with specific activity against V. cholerae O1 strains and is grouped as a member of the VP2-like phage subfamily. The encoding of an Ig domain by ØVC8 makes this phage a good candidate for use in phage therapy and an alternative tool for monitoring V. cholerae populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0490-x) contains supplementary material, which is available to authorized users.

Comparative Genomic and Phylogenomic Analyses Reveal a Conserved Core Genome Shared by Estuarine and Oceanic Cyanopodoviruses

Podoviruses are among the major viral groups that infect marine picocyanobacteria Prochlorococcus and Synechococcus. Here, we reported the genome sequences of five Synechococcus podoviruses isolated from the estuarine environment, and performed comparative genomic and phylogenomic analyses based on a total of 20 cyanopodovirus genomes. The genomes of all the known marine cyanopodoviruses are highly syntenic. A pan-genome of 349 clustered orthologous groups was determined, among which 15 were core genes. These core genes make up nearly half of each genome in length, reflecting the high level of genome conservation among this cyanophage type. The whole genome phylogenies based on concatenated core genes and gene content were highly consistent and confirmed the separation of two discrete marine cyanopodovirus clusters MPP-A and MPP-B. The genomes within cluster MPP-B grouped into subclusters mainly corresponding to Prochlorococcus or Synechococcus host types. Auxiliary metabolic genes tend to occur in a specific phylogenetic group of these cyanopodoviruses. All the MPP-B phages analyzed here encode the photosynthesis gene psbA, which are absent in all the MPP-A genomes thus far. Interestingly, all the MPP-B and two MPP-A Synechococcus podoviruses encode the thymidylate synthase gene thyX, while at the same genome locus all the MPP-B Prochlorococcus podoviruses encode the transaldolase gene talC. Both genes are hypothesized to have the potential to facilitate the biosynthesis of deoxynucleotide for phage replication. Inheritance of specific functional genes could be important to the evolution and ecological fitness of certain cyanophage genotypes. Our analyses demonstrate that cyanopodoviruses of estuarine and oceanic origins share a conserved core genome and suggest that accessory genes may be related to environmental adaptation.

Unraveling the genome structure of cyanobacterial podovirus A-4L with long direct terminal repeats

The freshwater cyanobacterial virus (cyanophage) A-4L, a podovirus, can infect the model cyanobacterium Anabaena sp. strain PCC 7120 resulting in a high burst size and forming concentric plaques on its lawns. The complete genome sequence of A-4L was determined by the combination of high-throughput sequencing, terminal transferase-mediated polymerase chain reaction and restriction mapping. It contains 41,750 bp with 810 bp direct terminal repeats and 38 potential open reading frames. As compared with other cyanobacterial podoviruses in diverse ecosystems, the A-4L has the longest terminal repeat and shares similar genome organizations with freshwater members. Furthermore, phylogenetic analysis based on concatenated sequences of eight core proteins indicated that freshwater cyanobacterial podoviruses were clustered together and distinct from marine counterparts, suggesting a clear divergence in the cyanobacterial podovirus lineage between freshwater and marine ecosystems. Our findings uncover the unique genome structure of A-4L which contains long direct terminal repeats, and create the first model system to address knowledge gaps in understanding cyanobacterial virus-host interactions at the molecular level.

Targeted diversity generation by intraterrestrial archaea and archaeal viruses

In the evolutionary arms race between microbes, their parasites, and their neighbours, the capacity for rapid protein diversification is a potent weapon. Diversity-generating retroelements (DGRs) use mutagenic reverse transcription and retrohoming to generate myriad variants of a target gene. Originally discovered in pathogens, these retroelements have been identified in bacteria and their viruses, but never in archaea. Here we report the discovery of intact DGRs in two distinct intraterrestrial archaeal systems: a novel virus that appears to infect archaea in the marine subsurface, and, separately, two uncultivated nanoarchaea from the terrestrial subsurface. The viral DGR system targets putative tail fibre ligand-binding domains, potentially generating >10¹⁸ protein variants. The two single-cell nanoarchaeal genomes each possess ≥4 distinct DGRs. Against an expected background of low genome-wide mutation rates, these results demonstrate a previously unsuspected potential for rapid, targeted sequence diversification in intraterrestrial archaea and their viruses.

Diversity-generating retroelements (DGRs) are genetic elements that introduce sequence variation within target genes in bacteria and their viruses. Here, Paul et al. report the discovery of DGRs in an archaeal virus and in two archaea from marine and terrestrial subsurface environments, respectively.

Life-style and genome structure of marine Pseudoalteromonas siphovirus B8b isolated from the northwestern Mediterranean Sea

Marine viruses (phages) alter bacterial diversity and evolution with impacts on marine biogeochemical cycles, and yet few well-developed model systems limit opportunities for hypothesis testing. Here we isolate phage B8b from the Mediterranean Sea using Pseudoalteromonas sp. QC-44 as a host and characterize it using myriad techniques. Morphologically, phage B8b was classified as a member of the Siphoviridae family. One-step growth analyses showed that this siphovirus had a latent period of 70 min and released 172 new viral particles per cell. Host range analysis against 89 bacterial host strains revealed that phage B8b infected 3 Pseudoalteromonas strains (52 tested, >99.9% 16S rRNA gene nucleotide identity) and 1 non-Pseudoaltermonas strain belonging to Alteromonas sp. (37 strains from 6 genera tested), which helps bound the phylogenetic distance possible in a phage-mediated horizontal gene transfer event. The Pseudoalteromonas phage B8b genome size was 42.7 kb, with clear structural and replication modules where the former were delineated leveraging identification of 16 structural genes by virion structural proteomics, only 4 of which had any similarity to known structural proteins. In nature, this phage was common in coastal marine environments in both photic and aphotic layers (found in 26.5% of available viral metagenomes), but not abundant in any sample (average per sample abundance was 0.65% of the reads). Together these data improve our understanding of siphoviruses in nature, and provide foundational information for a new ‘rare virosphere’ phage–host model system.

Purifying the impure: sequencing metagenomes and metatranscriptomes from complex animal-associated samples

The accessibility of high-throughput sequencing has revolutionized many fields of biology. In order to better understand host-associated viral and microbial communities, a comprehensive workflow for DNA and RNA extraction was developed. The workflow concurrently generates viral and microbial metagenomes, as well as metatranscriptomes, from a single sample for next-generation sequencing. The coupling of these approaches provides an overview of both the taxonomical characteristics and the community encoded functions. The presented methods use Cystic Fibrosis (CF) sputum, a problematic sample type, because it is exceptionally viscous and contains high amount of mucins, free neutrophil DNA, and other unknown contaminants. The protocols described here target these problems and successfully recover viral and microbial DNA with minimal human DNA contamination. To complement the metagenomics studies, a metatranscriptomics protocol was optimized to recover both microbial and host mRNA that contains relatively few ribosomal RNA (rRNA) sequences. An overview of the data characteristics is presented to serve as a reference for assessing the success of the methods. Additional CF sputum samples were also collected to (i) evaluate the consistency of the microbiome profiles across seven consecutive days within a single patient, and (ii) compare the consistency of metagenomic approach to a 16S ribosomal RNA gene-based sequencing. The results showed that daily fluctuation of microbial profiles without antibiotic perturbation was minimal and the taxonomy profiles of the common CF-associated bacteria were highly similar between the 16S rDNA libraries and metagenomes generated from the hypotonic lysis (HL)-derived DNA. However, the differences between 16S rDNA taxonomical profiles generated from total DNA and HL-derived DNA suggest that hypotonic lysis and the washing steps benefit in not only removing the human-derived DNA, but also microbial-derived extracellular DNA that may misrepresent the actual microbial profiles.

Unraveling the unseen players in the ocean - a field guide to water chemistry and marine microbiology

Here we introduce a series of thoroughly tested and well standardized research protocols adapted for use in remote marine environments. The sampling protocols include the assessment of resources available to the microbial community (dissolved organic carbon, particulate organic matter, inorganic nutrients), and a comprehensive description of the viral and bacterial communities (via direct viral and microbial counts, enumeration of autofluorescent microbes, and construction of viral and microbial metagenomes). We use a combination of methods, which represent a dispersed field of scientific disciplines comprising already established protocols and some of the most recent techniques developed. Especially metagenomic sequencing techniques used for viral and bacterial community characterization, have been established only in recent years, and are thus still subjected to constant improvement. This has led to a variety of sampling and sample processing procedures currently in use. The set of methods presented here provides an up to date approach to collect and process environmental samples. Parameters addressed with these protocols yield the minimum on information essential to characterize and understand the underlying mechanisms of viral and microbial community dynamics. It gives easy to follow guidelines to conduct comprehensive surveys and discusses critical steps and potential caveats pertinent to each technique.

Ribonucleotide reductases reveal novel viral diversity and predict biological and ecological features of unknown marine viruses

Virioplankton play a crucial role in aquatic ecosystems as top-down regulators of bacterial populations and agents of horizontal gene transfer and nutrient cycling. However, the biology and ecology of virioplankton populations in the environment remain poorly understood. Ribonucleotide reductases (RNRs) are ancient enzymes that reduce ribonucleotides to deoxyribonucleotides and thus prime DNA synthesis. Composed of three classes according to O2 reactivity, RNRs can be predictive of the physiological conditions surrounding DNA synthesis. RNRs are universal among cellular life, common within viral genomes and virioplankton shotgun metagenomes (viromes), and estimated to occur within >90% of the dsDNA virioplankton sampled in this study. RNRs occur across diverse viral groups, including all three morphological families of tailed phages, making these genes attractive for studies of viral diversity. Differing patterns in virioplankton diversity were clear from RNRs sampled across a broad oceanic transect. The most abundant RNRs belonged to novel lineages of podoviruses infecting α-proteobacteria, a bacterial class critical to oceanic carbon cycling. RNR class was predictive of phage morphology among cyanophages and RNR distribution frequencies among cyanophages were largely consistent with the predictions of the "kill the winner-cost of resistance" model. RNRs were also identified for the first time to our knowledge within ssDNA viromes. These data indicate that RNR polymorphism provides a means of connecting the biological and ecological features of virioplankton populations.

Complete genome sequence of invertebrate iridovirus IIV-25 isolated from a blackfly larva

Members of the family Iridoviridae are animal viruses that infect only invertebrates and poikilothermic vertebrates. Invertebrate iridovirus 25 (IIV-25) was originally isolated from the larva of a blackfly (Simulium spp., order Diptera) found in the Ystwyth river near Aberystwyth, Wales. IIV-25 virions are icosahedral, have a diameter of ~130 nm, and contain a dsDNA genome of 204.8 kbp, with a G+C content of 30.32 %, that codes for 177 proteins. Here, we describe the complete genome sequence of this virus and its annotation. This is the fifth genome sequence of an invertebrate iridovirus reported.

Bacteriophage-based synthetic biology for the study of infectious diseases

Since their discovery, bacteriophages have contributed enormously to our understanding of molecular biology as model systems. Furthermore, bacteriophages have provided many tools that have advanced the fields of genetic engineering and synthetic biology. Here, we discuss bacteriophage-based technologies and their application to the study of infectious diseases. New strategies for engineering genomes have the potential to accelerate the design of novel phages as therapies, diagnostics, and tools. Though almost a century has elapsed since their discovery, bacteriophages continue to have a major impact on modern biological sciences, especially with the growth of multidrug-resistant bacteria and interest in the microbiome.

Complete genome sequence of invertebrate iridovirus IIV22A, a variant of IIV22, isolated originally from a blackfly larva

Members of the family Iridoviridae are animal viruses that infect only invertebrates and poikilothermic vertebrates. The invertebrate iridoviruses 22 (IIV22) and 25 (IIV25) were originally isolated from a single sample of blackfly larva (Simulium spp., order Diptera) collected from the Ystwyth river near Aberystwyth, Wales. Recently, the genomes of IIV22 (197.7 kbp) and IIV25 (204.8 kbp) were sequenced and reported. Here, we describe the complete genome sequence of IIV22A, a variant that was isolated from the same pool of virions collected from the blackfly larva from which the IIV22 virion genome originated. The IIV22A genome, 196.5 kbp, is smaller than IIV22. Nevertheless, it contains 7 supplementary putative ORFs. Its analysis enables evaluation of the degree of genomic polymorphisms within an IIV isolate. Despite the occurrence of this IIV variant with IIV22 and IIV25 in a single blackfly larva and the features of their DNA polymerase, we found no evidence of lateral genetic transfers between the genomes of these two IIV species.

Genome sequence of a crustacean iridovirus, IIV31, isolated from the pill bug, Armadillidium vulgare

Members of the family Iridoviridae are animal viruses that infect only invertebrates and poikilothermic vertebrates. The invertebrate iridovirus 31 (IIV31) was originally isolated from adult pill bugs, Armadillidium vulgare (class Crustacea, order Isopoda, suborder Oniscidea), found in southern California on the campus of the University of California, Riverside, USA. IIV31 virions are icosahedral, have a diameter of about 135 nm, and contain a dsDNA genome 220.222 kbp in length, with 35.09 mol % G+C content and 203 ORFs. Here, we describe the complete genome sequence of this virus and its annotation. This is the eighth genome sequence of an IIV reported.

Caught in the middle with multiple displacement amplification: the myth of pooling for avoiding multiple displacement amplification bias in a metagenome

BACKGROUND

Shotgun metagenomics has become an important tool for investigating the ecology of microorganisms. Underlying these investigations is the assumption that metagenome sequence data accurately estimates the census of microbial populations. Multiple displacement amplification (MDA) of microbial community DNA is often used in cases where it is difficult to obtain enough DNA for sequencing; however, MDA can result in amplification biases that may impact subsequent estimates of population census from metagenome data. Some have posited that pooling replicate MDA reactions negates these biases and restores the accuracy of population analyses. This assumption has not been empirically tested.

RESULTS

Using mock viral communities, we examined the influence of pooling on population-scale analyses. In pooled and single reaction MDA treatments, sequence coverage of viral populations was highly variable and coverage patterns across viral genomes were nearly identical, indicating that initial priming biases were reproducible and that pooling did not alleviate biases. In contrast, control unamplified sequence libraries showed relatively even coverage across phage genomes.

CONCLUSIONS

MDA should be avoided for metagenomic investigations that require quantitative estimates of microbial taxa and gene functional groups. While MDA is an indispensable technique in applications such as single-cell genomics, amplification biases cannot be overcome by combining replicate MDA reactions. Alternative library preparation techniques should be utilized for quantitative microbial ecology studies utilizing metagenomic sequencing approaches.

Complete genome sequence of invertebrate iridovirus IIV30 isolated from the corn earworm, Helicoverpa zea

Members of the family Iridoviridae are animal viruses that infect only invertebrates and poikilothermic vertebrates. The invertebrate iridovirus 30 (IIV30) was originally isolated from a larva of the corn earworm, Helicoverpa zea (order lepidoptera, Family Noctuidae) in western Australia. The IIV30 virions are icosahedral, have a diameter of about 130nm, and contain a dsDNA genome of 198.5kbp with 28.11% in GC content and 177 coding sequences. Here we describe its complete genome sequence and annotate the genes for which we could assign a putative function. This is the sixth genome sequence of an invertebrate iridovirus reported.

Metabolic reprogramming by viruses in the sunlit and dark ocean

Background

Marine ecosystem function is largely determined by matter and energy transformations mediated by microbial community interaction networks. Viral infection modulates network properties through mortality, gene transfer and metabolic reprogramming.

Results

Here we explore the nature and extent of viral metabolic reprogramming throughout the Pacific Ocean depth continuum. We describe 35 marine viral gene families with potential to reprogram metabolic flux through central metabolic pathways recovered from Pacific Ocean waters. Four of these families have been previously reported but 31 are novel. These known and new carbon pathway auxiliary metabolic genes were recovered from a total of 22 viral metagenomes in which viral auxiliary metabolic genes were differentiated from low-level cellular DNA inputs based on small subunit ribosomal RNA gene content, taxonomy, fragment recruitment and genomic context information. Auxiliary metabolic gene distribution patterns reveal that marine viruses target overlapping, but relatively distinct pathways in sunlit and dark ocean waters to redirect host carbon flux towards energy production and viral genome replication under low nutrient, niche-differentiated conditions throughout the depth continuum.

Conclusions

Given half of ocean microbes are infected by viruses at any given time, these findings of broad viral metabolic reprogramming suggest the need for renewed consideration of viruses in global ocean carbon models.

Complete genome sequence of invertebrate iridescent virus 22 isolated from a blackfly larva

Members of the family Iridoviridae are animal viruses that infect only invertebrates and poikilothermic vertebrates. Invertebrate iridescent virus 22 (IIV-22) was originally isolated from the larva of a blackfly (Simulium sp., order Diptera) found in the Ystwyth river, near Aberystwyth, Wales, UK. IIV-22 virions are icosahedral, with a diameter of about 130 nm and contain a dsDNA genome that is 197.7 kb in length, has a G+C content of 28.05 mol% and contains 167 coding sequences. Here, we describe the complete genome sequence of this virus and its annotation. This is the fourth genome sequence of an invertebrate iridovirus to be reported.

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.

From deep sequencing to viral tagging: recent advances in viral metagenomics

Culture-independent high-throughput sequencing has provided unprecedented insights into microbial ecology, particularly for Earth's most ubiquitous and diverse inhabitants - the viruses. A plethora of methods now exist for amplifying the vanishingly small amounts of nucleic acids in natural viral communities in order to sequence them, and sequencing depth is now so great that viral genomes can be detected and assembled even amid large concentrations of non-viral DNA. Complementing these advances in amplification and sequencing is the ability to physically link fluorescently labeled viruses to their host cells via high-throughput flow sorting. Sequencing of such isolated virus-host pairs facilitates cultivation-independent exploration of the natural host range of viruses. Within the next decade, as these technologies become widespread, we can expect to see a systematic expansion of our knowledge of viruses and their hosts.

The Pacific Ocean virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology

Bacteria and their viruses (phage) are fundamental drivers of many ecosystem processes including global biogeochemistry and horizontal gene transfer. While databases and resources for studying function in uncultured bacterial communities are relatively advanced, many fewer exist for their viral counterparts. The issue is largely technical in that the majority (often 90%) of viral sequences are functionally 'unknown' making viruses a virtually untapped resource of functional and physiological information. Here, we provide a community resource that organizes this unknown sequence space into 27 K high confidence protein clusters using 32 viral metagenomes from four biogeographic regions in the Pacific Ocean that vary by season, depth, and proximity to land, and include some of the first deep pelagic ocean viral metagenomes. These protein clusters more than double currently available viral protein clusters, including those from environmental datasets. Further, a protein cluster guided analysis of functional diversity revealed that richness decreased (i) from deep to surface waters, (ii) from winter to summer, (iii) and with distance from shore in surface waters only. These data provide a framework from which to draw on for future metadata-enabled functional inquiries of the vast viral unknown.

Preparation of metagenomic libraries from naturally occurring marine viruses

Microbes are now well recognized as major drivers of the biogeochemical cycling that fuels the Earth, and their viruses (phages) are known to be abundant and important in microbial mortality, horizontal gene transfer, and modulating microbial metabolic output. Investigation of environmental phages has been frustrated by an inability to culture the vast majority of naturally occurring diversity coupled with the lack of robust, quantitative, culture-independent methods for studying this uncultured majority. However, for double-stranded DNA phages, a quantitative viral metagenomic sample-to-sequence workflow now exists. Here, we review these advances with special emphasis on the technical details of preparing DNA sequencing libraries for metagenomic sequencing from environmentally relevant low-input DNA samples. Library preparation steps broadly involve manipulating the sample DNA by fragmentation, end repair and adaptor ligation, size fractionation, and amplification. One critical area of future research and development is parallel advances for alternate nucleic acid types such as single-stranded DNA and RNA viruses that are also abundant in nature. Combinations of recent advances in fragmentation (e.g., acoustic shearing and tagmentation), ligation reactions (adaptor-to-template ratio reference table availability), size fractionation (non-gel-sizing), and amplification (linear amplification for deep sequencing and linker amplification protocols) enhance our ability to generate quantitatively representative metagenomic datasets from low-input DNA samples. Such datasets are already providing new insights into the role of viruses in marine systems and will continue to do so as new environments are explored and synergies and paradigms emerge from large-scale comparative analyses.

Genomic analysis of cold-active Colwelliaphage 9A and psychrophilic phage-host interactions

The 104 kb genome of cold-active bacteriophage 9A, which replicates in the marine psychrophilic gamma-proteobacterium Colwellia psychrerythraea strain 34H (between -12 and 8 °C), was sequenced and analyzed to investigate elements of molecular adaptation to low temperature and phage-host interactions in the cold. Most characterized ORFs indicated closest similarity to gamma-proteobacteria and their phages, though no single module provided definitive phylogenetic grouping. A subset of primary structural features linked to psychrophily suggested that the majority of annotated phage proteins were not psychrophilic; those that were, primarily serve phage-specific functions and may also contribute to 9A's restricted temperature range for replication as compared to host. Comparative analyses suggest ribonucleotide reductase genes were acquired laterally from host. Neither restriction modification nor the CRISPR-Cas system appeared to be the predominant phage defense mechanism of Cp34H or other cold-adapted bacteria; we hypothesize that psychrophilic hosts rely more on the use of extracellular polymeric material to block cell surface receptors recognized by phages. The relative dearth of evidence for genome-specific defenses, genetic transfer events or auxiliary metabolic genes suggest that the 9A-Cp34H system may be less tightly coupled than are other genomically characterized marine phage-host systems, with possible implications for phage specificity under different environmental conditions.

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.

Genomic and functional analysis of Vibrio phage SIO-2 reveals novel insights into ecology and evolution of marine siphoviruses

We report on a genomic and functional analysis of a novel marine siphovirus, the Vibrio phage SIO-2. This phage is lytic for related Vibrio species of great ecological interest including the broadly antagonistic bacterium Vibrio sp. SWAT3 as well as notable members of the Harveyi clade (V.harveyi ATTC BAA-1116 and V.campbellii ATCC 25920). Vibrio phage SIO-2 has a circularly permuted genome of 80598 bp, which displays unusual features. This genome is larger than that of most known siphoviruses and only 38 of the 116 predicted proteins had homologues in databases. Another divergence is manifest by the origin of core genes, most of which share robust similarities with unrelated viruses and bacteria spanning a wide range of phyla. These core genes are arranged in the same order as in most bacteriophages but they are unusually interspaced at two places with insertions of DNA comprising a high density of uncharacterized genes. The acquisition of these DNA inserts is associated with morphological variation of SIO-2 capsid, which assembles as a large (80 nm) shell with a novel T=12 symmetry. These atypical structural features confer on SIO-2 a remarkable stability to a variety of physical, chemical and environmental factors. Given this high level of functional and genomic novelty, SIO-2 emerges as a model of considerable interest in ecological and evolutionary studies.

Evaluation of methods to concentrate and purify ocean virus communities through comparative, replicated metagenomics

Viruses have global impact through mortality, nutrient cycling and horizontal gene transfer, yet their study is limited by complex methodologies with little validation. Here, we use triplicate metagenomes to compare common aquatic viral concentration and purification methods across four combinations as follows: (i) tangential flow filtration (TFF) and DNase + CsCl, (ii) FeCl₃ precipitation and DNase, (iii) FeCl₃ precipitation and DNase + CsCl and (iv) FeCl₃ precipitation and DNase + sucrose. Taxonomic data (30% of reads) suggested that purification methods were statistically indistinguishable at any taxonomic level while concentration methods were significantly different at family and genus levels. Specifically, TFF-concentrated viral metagenomes had significantly fewer abundant viral types (Podoviridae and Phycodnaviridae) and more variability among Myoviridae than FeCl₃-precipitated viral metagenomes. More comprehensive analyses using protein clusters (66% of reads) and k-mers (100% of reads) showed 50–53% of these data were common to all four methods, and revealed trace bacterial DNA contamination in TFF-concentrated metagenomes and one of three replicates concentrated using FeCl₃ and purified by DNase alone. Shared k-mer analyses also revealed that polymerases used in amplification impact the resulting metagenomes, with TaKaRa enriching for ‘rare’ reads relative to PfuTurbo. Together these results provide empirical data for making experimental design decisions in culture-independent viral ecology studies.

Towards quantitative metagenomics of wild viruses and other ultra-low concentration DNA samples: a rigorous assessment and optimization of the linker amplification method

Metagenomics generates and tests hypotheses about dynamics and mechanistic drivers in wild populations, yet commonly suffers from insufficient (< 1 ng) starting genomic material for sequencing. Current solutions for amplifying sufficient DNA for metagenomics analyses include linear amplification for deep sequencing (LADS), which requires more DNA than is normally available, linker-amplified shotgun libraries (LASLs), which is prohibitively low throughput, and whole-genome amplification, which is significantly biased and thus non-quantitative. Here, we adapt the LASL approach to next generation sequencing by offering an alternate polymerase for challenging samples, developing a more efficient sizing step, integrating a 'reconditioning PCR' step to increase yield and minimize late-cycle PCR artefacts, and empirically documenting the quantitative capability of the optimized method with both laboratory isolate and wild community viral DNA. Our optimized linker amplification method requires as little as 1 pg of DNA and is the most precise and accurate available, with G + C content amplification biases less than 1.5-fold, even for complex samples as diverse as a wild virus community. While optimized here for 454 sequencing, this linker amplification method can be used to prepare metagenomics libraries for sequencing with next-generation platforms, including Illumina and Ion Torrent, the first of which we tested and present data for here.

A novel virus genome discovered in an extreme environment suggests recombination between unrelated groups of RNA and DNA viruses

BACKGROUND

Viruses are known to be the most abundant organisms on earth, yet little is known about their collective origin and evolutionary history. With exceptionally high rates of genetic mutation and mosaicism, it is not currently possible to resolve deep evolutionary histories of the known major virus groups. Metagenomics offers a potential means of establishing a more comprehensive view of viral evolution as vast amounts of new sequence data becomes available for comparative analysis.

RESULTS

Bioinformatic analysis of viral metagenomic sequences derived from a hot, acidic lake revealed a circular, putatively single-stranded DNA virus encoding a major capsid protein similar to those found only in single-stranded RNA viruses. The presence and circular configuration of the complete virus genome was confirmed by inverse PCR amplification from native DNA extracted from lake sediment. The virus genome appears to be the result of a RNA-DNA recombination event between two ostensibly unrelated virus groups. Environmental sequence databases were examined for homologous genes arranged in similar configurations and three similar putative virus genomes from marine environments were identified. This result indicates the existence of a widespread but previously undetected group of viruses.

CONCLUSIONS

This unique viral genome carries implications for theories of virus emergence and evolution, as no mechanism for interviral RNA-DNA recombination has yet been identified, and only scant evidence exists that genetic exchange occurs between such distinct virus lineages.

REVIEWERS

This article was reviewed by EK, MK (nominated by PF) and AM. For the full reviews, please go to the Reviewers' comments section.

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

Background

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

Results

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

Conclusions

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

The P-SSP7 cyanophage has a linear genome with direct terminal repeats

P-SSP7 is a T7-like phage that infects the cyanobacterium Prochlorococcus MED4. MED4 is a member of the high-light-adapted Prochlorococcus ecotypes that are abundant in the surface oceans and contribute significantly to primary production. P-SSP7 has become a model system for the investigation of T7-like phages that infect Prochlorococcus. It was classified as T7-like based on genome content and organization. However, because its genome assembled as a circular molecule, it was thought to be circularly permuted and to lack the direct terminal repeats found in other T7-like phages. Here we sequenced the ends of the P-SSP7 genome and found that the genome map is linear and contains a 206 bp repeat at both genome ends. Furthermore, we found that a 728 bp region of the genome originally placed downstream of the last ORF is actually located upstream of the first ORF on the genome map. These findings suggest that P-SSP7 is likely to use the direct terminal repeats for genome replication and packaging in a similar manner to other T7-like phages. Moreover, these results highlight the importance of experimentally verifying the ends of phage genomes, and will facilitate the use of P-SSP7 as a model for the correct assembly and end determination of the many T7-like phages isolated from the marine environment that are currently being sequenced.

Draft genome sequence of four coccolithoviruses: Emiliania huxleyi virus EhV-88, EhV-201, EhV-207, and EhV-208

The Coccolithoviridae are a group of viruses which infect the marine coccolithophorid microalga Emiliania huxleyi. The Emiliania huxleyi viruses (known as EhVs) described herein have 160- to 180-nm diameter icosahedral structures, have genomes of approximately 400 kbp, and consist of more than 450 predicted coding sequences (CDSs). Here, we describe the genomic features of four newly sequenced coccolithoviruses (EhV-88, EhV-201, EhV-207, and EhV-208) together with their draft genome sequences and their annotations, highlighting the homology and heterogeneity of these genomes to the EhV-86 model reference genome.

Grinder: a versatile amplicon and shotgun sequence simulator

We introduce Grinder (http://sourceforge.net/projects/biogrinder/), an open-source bioinformatic tool to simulate amplicon and shotgun (genomic, metagenomic, transcriptomic and metatranscriptomic) datasets from reference sequences. This is the first tool to simulate amplicon datasets (e.g. 16S rRNA) widely used by microbial ecologists. Grinder can create sequence libraries with a specific community structure, α and β diversities and experimental biases (e.g. chimeras, gene copy number variation) for commonly used sequencing platforms. This versatility allows the creation of simple to complex read datasets necessary for hypothesis testing when developing bioinformatic software, benchmarking existing tools or designing sequence-based experiments. Grinder is particularly useful for simulating clinical or environmental microbial communities and complements the use of in vitro mock communities.