K. Eric Wommack

Metagenomics: read length matters.

ABSTRACT Obtaining an unbiased view of the phylogenetic composition and functional diversity within a microbial community is one central objective of metagenomic analysis. New technologies, such as 454 pyrosequencing, have dramatically reduced sequencing costs, to a level where metagenomic analysis may become a viable alternative to more-focused assessments of the phylogenetic (e.g., 16S rRNA genes) and functional diversity of microbial communities. To determine whether the short (∼100 to 200 bp) sequence reads obtained from pyrosequencing are appropriate for the phylogenetic and functional characterization of microbial communities, the results of BLAST and COG analyses were compared for long (∼750 bp) and randomly derived short reads from each of two microbial and one virioplankton metagenome libraries. Overall, BLASTX searches against the GenBank nr database found far fewer homologs within the short-sequence libraries. This was especially pronounced for a Chesapeake Bay virioplankton metagenome library. Increasing the short-read sampling depth or the length of derived short reads (up to 400 bp) did not completely resolve the discrepancy in BLASTX homolog detection. Only in cases where the long-read sequence had a close homolog (low BLAST E-score) did the derived short-read sequence also find a significant homolog. Thus, more-distant homologs of microbial and viral genes are not detected by short-read sequences. Among COG hits, derived short reads sampled at a depth of two short reads per long read missed up to 72% of the COG hits found using long reads. Noting the current limitation in computational approaches for the analysis of short sequences, the use of short-read-length libraries does not appear to be an appropriate tool for the metagenomic characterization of microbial communities.

Metagenomic Characterization of Chesapeake Bay Virioplankton

ABSTRACT Viruses are ubiquitous and abundant throughout the biosphere. In marine systems, virus-mediated processes can have significant impacts on microbial diversity and on global biogeocehmical cycling. However, viral genetic diversity remains poorly characterized. To address this shortcoming, a metagenomic library was constructed from Chesapeake Bay virioplankton. The resulting sequences constitute the largest collection of long-read double-stranded DNA (dsDNA) viral metagenome data reported to date. BLAST homology comparisons showed that Chesapeake Bay virioplankton contained a high proportion of unknown (homologous only to environmental sequences) and novel (no significant homolog) sequences. This analysis suggests that dsDNA viruses are likely one of the largest reservoirs of unknown genetic diversity in the biosphere. The taxonomic origin of BLAST homologs to viral library sequences agreed well with reported abundances of cooccurring bacterial subphyla within the estuary and indicated that cyanophages were abundant. However, the low proportion of Siphophage homologs contradicts a previous assertion that this family comprises most bacteriophage diversity. Identification and analyses of cyanobacterial homologs of the psbA gene illustrated the value of metagenomic studies of virioplankton. The phylogeny of inferred PsbA protein sequences suggested that Chesapeake Bay cyanophage strains are endemic in that environment. The ratio of psbA homologous sequences to total cyanophage sequences in the metagenome indicated that the psbA gene may be nearly universal in Chesapeake Bay cyanophage genomes. Furthermore, the low frequency of psbD homologs in the library supports the prediction that Chesapeake Bay cyanophage populations are dominated by Podoviridae.

Assembly of viral metagenomes from yellowstone hot springs.

ABSTRACT Thermophilic viruses were reported decades ago; however, knowledge of their diversity, biology, and ecological impact is limited. Previous research on thermophilic viruses focused on cultivated strains. This study examined metagenomic profiles of viruses directly isolated from two mildly alkaline hot springs, Bear Paw (74°C) and Octopus (93°C). Using a new method for constructing libraries from picograms of DNA, nearly 30 Mb of viral DNA sequence was determined. In contrast to previous studies, sequences were assembled at 50% and 95% identity, creating composite contigs up to 35 kb and facilitating analysis of the inherent heterogeneity in the populations. Lowering the assembly identity reduced the estimated number of viral types from 1,440 and 1,310 to 548 and 283, respectively. Surprisingly, the diversity of viral species in these springs approaches that in moderate-temperature environments. While most known thermophilic viruses have a chronic, nonlytic infection lifestyle, analysis of coding sequences suggests lytic viruses are more common in geothermal environments than previously thought. The 50% assembly included one contig with high similarity and perfect synteny to nine genes from Pyrobaculum spherical virus (PSV). In fact, nearly all the genes of the 28-kb genome of PSV have apparent homologs in the metagenomes. Similarities to thermoacidophilic viruses isolated on other continents were limited to specific open reading frames but were equally strong. Nearly 25% of the reads showed significant similarity between the hot springs, suggesting a common subterranean source. To our knowledge, this is the first application of metagenomics to viruses of geothermal origin.

Abundance and Diversity of Viruses in Six Delaware Soils

ABSTRACT The importance of viruses in marine microbial ecology has been established over the past decade. Specifically, viruses influence bacterial abundance and community composition through lysis and alter bacterial genetic diversity through transduction and lysogenic conversion. By contrast, the abundance and distribution of viruses in soils are almost completely unknown. This study describes the abundance and diversity of autochthonous viruses in six Delaware soils: two agricultural soils, two coastal plain forest soils, and two piedmont forest soils. Viral abundance was measured using epifluorescence microscopy, while viral diversity was assessed from morphological data obtained through transmission electron microscopy. Extracted soil virus communities were dominated by bacteriophages that demonstrated a wide range of capsid diameters (20 nm to 160 nm) and morphologies, including filamentous forms and phages with elongated capsids. The reciprocal Simpsons index suggests that forest soils harbor more diverse assemblages of viruses, particularly in terms of morphological distribution. Repeated extractions of virus-like particles (VLPs) from soils indicated that the initial round of extraction removes approximately 70% of extractable viruses. Higher VLP abundances were observed in forest soils (1.31 × 109 to 4.17 × 109 g−1 dry weight) than in agricultural soils (8.7 × 108 to 1.1 × 109 g−1 dry weight). Soil VLP abundance was significantly correlated to moisture content (r = 0.988) but not to soil texture. Land use (agricultural or forested) was significantly correlated to both bacterial (r = 0.885) and viral (r = 0.812) abundances, as were soil organic matter and water content. Thus, land use is a significant factor influencing viral abundance and diversity in soils.

Phages across the biosphere: contrasts of viruses in soil and aquatic environments

Despite the predominance of aquatic environments on the planet Earth, microbial abundance and diversity within soil environments exceed that of the aquatic realm. Most of what we know of viral ecology within natural systems has come through investigations of aquatic environments. However, the aquatic-bias in viral ecology is beginning to change as the cultivation-independent approaches, which revealed the extraordinary abundance and diversity of viruses within aquatic systems, are now being applied to soils. This review briefly summarizes recent investigations of viral abundance and diversity in soil environments.

Sampling Natural Viral Communities from Soil for Culture-Independent Analyses

ABSTRACT An essential first step in investigations of viruses in soil is the evaluation of viral recovery methods suitable for subsequent culture-independent analyses. In this study, four elution buffers (10% beef extract, 250 mM glycine buffer, 10 mM sodium pyrophosphate, and 1% potassium citrate) and three enumeration techniques (plaque assay, epifluorescence microscopy [EFM], and transmission electron microscopy [TEM]) were compared to determine the best method of extracting autochthonous bacteriophages from two Delaware agricultural soils. Beef extract and glycine buffer were the most effective in eluting viable phages inoculated into soils (up to 29% recovery); however, extraction efficiency varied significantly with phage strain. Potassium citrate eluted the highest numbers of virus-like particles from both soils based on enumerations by EFM (mean, 5.3 × 108 g of dry soil−1), but specific soil-eluant combinations posed significant problems to enumeration by EFM. Observations of virus-like particles under TEM gave confidence that the particles were, in fact, phages, but TEM enumerations yielded measurements of phage abundance (mean, 1.5×108 g of dry soil−1) that were about five times lower. Clearly, the measurement of phage abundance in soils varies with both the extraction and enumeration methodology; thus, it is important to assess multiple extraction and enumeration approaches prior to undertaking ecological studies of phages in a particular soil.

VIROME: a standard operating procedure for analysis of viral metagenome sequences

One consistent finding among studies using shotgun metagenomics to analyze whole viral communities is that most viral sequences show no significant homology to known sequences. Thus, bioinformatic analyses based on sequence collections such as GenBank nr, which are largely comprised of sequences from known organisms, tend to ignore a majority of sequences within most shotgun viral metagenome libraries. Here we describe a bioinformatic pipeline, the Viral Informatics Resource for Metagenome Exploration (VIROME), that emphasizes the classification of viral metagenome sequences (predicted open-reading frames) based on homology search results against both known and environmental sequences. Functional and taxonomic information is derived from five annotated sequence databases which are linked to the UniRef 100 database. Environmental classifications are obtained from hits against a custom database, MetaGenomes On-Line, which contains 49 million predicted environmental peptides. Each predicted viral metagenomic ORF run through the VIROME pipeline is placed into one of seven ORF classes, thus, every sequence receives a meaningful annotation. Additionally, the pipeline includes quality control measures to remove contaminating and poor quality sequence and assesses the potential amount of cellular DNA contamination in a viral metagenome library by screening for rRNA genes. Access to the VIROME pipeline and analysis results are provided through a web-application interface that is dynamically linked to a relational back-end database. The VIROME web-application interface is designed to allow users flexibility in retrieving sequences (reads, ORFs, predicted peptides) and search results for focused secondary analyses.

Diverse and Unique Picocyanobacteria in Chesapeake Bay, Revealed by 16S-23S rRNA Internal Transcribed Spacer Sequences

ABSTRACT rRNA internal transcribed spacer phylogeny showed that Chesapeake Bay is populated with diverse Synechococcus strains, including members of the poorly studied marine cluster B. Marine cluster B prevailed in the upper bay, while marine cluster A was common in the lower bay. Interestingly, marine cluster B Synechococcus included phycocyanin- and phycoerythrin-rich strains.

Lysogenic virus-host interactions predominate at deep-sea diffuse-flow hydrothermal vents.

The consequences of viral infection within microbial communities are dependent on the nature of the viral life cycle. Among the possible outcomes is the substantial influence of temperate viruses on the phenotypes of lysogenic prokaryotes through various forms of genetic exchange. To date, no marine microbial ecosystem has consistently shown a predisposition for containing significant numbers of inducible temperate viruses. Here, we show that deep-sea diffuse-flow hydrothermal vent waters display a consistently high incidence of lysogenic hosts and harbor substantial populations of temperate viruses. Genetic fingerprinting and initial metagenomic analyses indicate that temperate viruses in vent waters appear to be a less diverse subset of the larger virioplankton community and that these viral populations contain an extraordinarily high frequency of novel genes. Thus, it appears likely that temperate viruses are key players in the ecology of prokaryotes within the extreme geothermal ecosystems of the deep sea.

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

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