Metagenomics: Read Length Matters

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.     

Prevalence of Viral Photosynthetic and Capsid Protein Genes from Cyanophages in Two Large and Deep Perialpine Lakes

Cyanophages are important components of aquatic ecosystems, but their genetic diversity has been little investigated in freshwaters. A yearlong survey was conducted in surface waters of the two largest natural perialpine lakes in France (Lake Annecy and Lake Bourget) to investigate part of this cyanophage diversity through the analysis of both structural (e.g., g20) and functional (e.g., psbA) genes. We found that these cyanophage signature genes were prevalent throughout the year but that the community compositions of g20 cyanomyoviruses were significantly different between the two lakes. In contrast, psbA-containing cyanophages seemed to be more similar between the two ecosystems. We also found that a large proportion of g20 sequences grouped with cyanomyophage isolates. psbA sequences, belonging to phages of Synechococcus spp., were characterized by distinct triplet motifs (with a novel viral triplet motif, EFE). Thus, our results show that cyanophages (i) are a diverse viral community in alpine lakes and (ii) are clearly distinct from some other freshwater and marine environments, suggesting the influence of unique biogeographic factors.     

Phylogenetic Diversity of Sequences of Cyanophage Photosynthetic Gene psbA in Marine and Freshwaters

Many cyanophage isolates which infect the marine cyanobacteria Synechococcus spp. and Prochlorococcus spp. contain a gene homologous to psbA, which codes for the D1 protein involved in photosynthesis. In the present study, cyanophage psbA gene fragments were readily amplified from freshwater and marine samples, confirming their widespread occurrence in aquatic communities. Phylogenetic analyses demonstrated that sequences from freshwaters have an evolutionary history that is distinct from that of their marine counterparts. Similarly, sequences from cyanophages infecting Prochlorococcus and Synechococcus spp. were readily discriminated, as were sequences from podoviruses and myoviruses. Viral psbA sequences from the same geographic origins clustered within different clades. For example, cyanophage psbA sequences from the Arctic Ocean fell within the Synechococcus as well as Prochlorococcus phage groups. Moreover, as psbA sequences are not confined to a single family of phages, they provide an additional genetic marker that can be used to explore the diversity and evolutionary history of cyanophages in aquatic environments.     

Diverse and dynamic populations of cyanobacterial podoviruses in the Chesapeake Bay unveiled through DNA polymerase gene sequences

Many podoviruses have been isolated which infect marine picocyanobacteria, and they may play a potentially important role in regulating the biomass and population composition of picocyanobacteria. However, little is known about the diversity and population dynamics of autochthonous cyanopodoviruses in marine environments. Using a set of newly designed PCR primers which specifically amplify the DNA pol from cyanopodoviruses, a total of 221 DNA pol sequences were retrieved from eight Chesapeake Bay virioplankton communities collected at different times and locations. All DNA pol sequences clustered with the eight known podoviruses that infect different marine picocyanobacteria, and could be divided into at least 10 different subclusters (I-X). The presence of these cyanopodovirus genotypes based on PCR-amplification of DNA pol gene sequences was supported by the existence of similar DNA pol genotypes with metagenome libraries of Chesapeake Bay virioplankton assemblages. The composition of cyanopodoviruses in the Bay also exhibited distinct winter and summer patterns which were likely related to corresponding seasonal changes in the composition of cyanobacterial populations. Our study suggests that diverse and dynamic populations of cyanopodoviruses are present in the estuarine environment. The PCR method developed in this study provides a specific and sensitive tool to explore the abundance, distribution and phylogenetic diversity of cyanopodoviruses in aquatic environments. Linking the dynamics of host and viral populations in the natural environment is critical to broader characterization of the ecological role of virioplankton within microbial communities.     

云南不同高原湖泊噬藻体psbA基因多样性

【目的】噬藻体(cyanophage)广泛存在于自然水体生态系统中,通过侵染蓝藻进而调控蓝藻种群及群落结构,具有重要生态功能和生态地位,在控制蓝藻水华方面有巨大开发潜力。本研究旨在探究云南高原湖泊噬藻体psbA基因多样性,分析其系统进化地位,为深入了解高原湖泊生态功能、开发利用噬藻体资源奠定理论基础。【方法】以云南高原主要湖泊滇池、抚仙湖和星云湖等为研究对象,以psbA基因作为分子靶标,对湖泊水体中噬藻体遗传多样性进行研究。【结果】从不同湖泊中共获得100条环境噬藻体psbA基因序列,系统发育分析表明,湖泊的噬藻体psbA基因序列与中国东湖、中国东北稻田、日本稻田等淡水中的环境噬藻体psbA基因亲缘关系较近,与海洋环境噬藻体psbA基因亲缘关系较远;抚仙湖中的噬藻体psbA基因多样性高于滇池、星云湖和异龙湖中的噬藻体psbA基因多样性;云南高原湖泊中存在新的噬藻体类群;各湖泊秋冬季节噬藻体psbA基因遗传多样性差异不明显。【结论】云南主要高原湖泊噬藻体psbA基因遗传多样性高,与淡水环境噬藻体psbA基因亲缘关系较近,且存在独特的噬藻体类群。     

Viral metagenomics analysis of planktonic viruses in East Lake, Wuhan, China

East Lake (Lake Donghu), located in Wuhan, China, is a typical city freshwater lake that has been experiencing eutrophic conditions and algal blooming during recent years. Marine and fresh water are considered to contain a large number of viruses. However, little is known about their genetic diversity because of the limited techniques for culturing viruses. In this study, we conducted a viral metagenomic analysis using a high-throughput sequencing technique with samples collected from East Lake in Spring, Summer, Autumn, and Winter. The libraries from four samples each generated 234,669, 71,837, 12,820, and 34,236 contigs (> 90 bp each), respectively. The genetic structure of the viral community revealed a high genetic diversity covering 23 viral families, with the majority of contigs homologous to DNA viruses, including members of Myoviridae, Podoviridae, Siphoviridae, Phycodnaviridae, and Microviridae, which infect bacteria or algae, and members of Circoviridae, which infect invertebrates and vertebrates. The highest viral genetic diversity occurred in samples collected in August, then December and June, and the least diversity in March. Most contigs have low-sequence identities with known viruses. PCR detection targeting the conserved sequences of genes (g20, psbA, psbD, and DNApol) of cyanophages further confirmed that there are novel cyanophages in the East Lake. Our viral metagenomic data provide the first preliminary understanding of the virome in one freshwater lake in China and would be helpful for novel virus discovery and the control of algal blooming in the future.     

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

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

Metagenomic analysis of viral community in the Yangtze River expands known eukaryotic and prokaryotic virus diversity in freshwater

Viruses in aquatic ecosystems are characterized by extraordinary abundance and diversity. Thus far, there have been limited studies focused on viral communities in river water systems. Here, we investigated the virome of the Yangtze River Delta using viral metagenomic analysis. The compositions of viral communities from six sampling sites were analyzed and compared. By using library construction and next generation sequencing, contigs and singlet reads similar to viral sequences were classified into 17 viral families, including nine dsDNA viral families, four ssDNA viral families and four RNA viral families. Statistical analysis using Friedman test suggested that there was no significant difference among the six sampling sites (P > 0.05). The viromes in this study were all dominated by the order Caudovirales, and a group of Freshwater phage uvFW species were particularly prevalent among all the samples. The virome from Nanjing presented a unique pattern of viral community composition with a relatively high abundance of family Parvoviridae. Phylogenetic analyses based on virus hallmark genes showed that the Caudovirales order and CRESS-DNA viruses presented high genetic diversity, while viruses in the Microviridae and Parvoviridae families and the Riboviria realm were relatively conservative. Our study provides the first insight into viral community composition in large river ecosystem, revealing the diversity and stability of river water virome, contributing to the proper utilization of freshwater resource.     

Metagenomic Exploration of Viruses throughout the Indian Ocean

The characterization of global marine microbial taxonomic and functional diversity is a primary goal of the Global Ocean Sampling Expedition. As part of this study, 19 water samples were collected aboard the Sorcerer II sailing vessel from the southern Indian Ocean in an effort to more thoroughly understand the lifestyle strategies of the microbial inhabitants of this ultra-oligotrophic region. No investigations of whole virioplankton assemblages have been conducted on waters collected from the Indian Ocean or across multiple size fractions thus far. Therefore, the goals of this study were to examine the effect of size fractionation on viral consortia structure and function and understand the diversity and functional potential of the Indian Ocean virome. Five samples were selected for comprehensive metagenomic exploration; and sequencing was performed on the microbes captured on 3.0-, 0.8- and 0.1 µm membrane filters as well as the viral fraction (<0.1 µm). Phylogenetic approaches were also used to identify predicted proteins of viral origin in the larger fractions of data from all Indian Ocean samples, which were included in subsequent metagenomic analyses. Taxonomic profiling of viral sequences suggested that size fractionation of marine microbial communities enriches for specific groups of viruses within the different size classes and functional characterization further substantiated this observation. Functional analyses also revealed a relative enrichment for metabolic proteins of viral origin that potentially reflect the physiological condition of host cells in the Indian Ocean including those involved in nitrogen metabolism and oxidative phosphorylation. A novel classification method, MGTAXA, was used to assess virus-host relationships in the Indian Ocean by predicting the taxonomy of putative host genera, with Prochlorococcus, Acanthochlois and members of the SAR86 cluster comprising the most abundant predictions. This is the first study to holistically explore virioplankton dynamics across multiple size classes and provides unprecedented insight into virus diversity, metabolic potential and virus-host interactions.     

Improving the genetic representation of rare taxa within complex microbial communities using DNA normalization methods

Complex microbial communities typically contain a large number of low abundance species, which collectively, comprise a considerable proportion of the community. This ‘rare biosphere’ has been speculated to contain keystone species and act as a repository of genomic diversity to facilitate community adaptation. Many environmental microbes are currently resistant to cultivation, and can only be accessed via culture‐independent approaches. To enhance our understanding of the role of the rare biosphere, we aimed to improve their metagenomic representation using DNA normalization methods, and assess normalization success via shotgun DNA sequencing. A synthetic metagenome was constructed from the genomic DNA of five bacterial species, pooled in a defined ratio spanning three orders of magnitude. The synthetic metagenome was fractionated and thermally renatured, allowing the most abundant sequences to hybridize. Double‐stranded DNA was removed either by hydroxyapatite chromatography, or by a duplex‐specific nuclease (DSN). The chromatographic method failed to enrich for the genomes present in low starting abundance, whereas the DSN method resulted in all genomes reaching near equimolar abundance. The representation of the rarest member was increased by approximately 450‐fold. De novo assembly of the normalized metagenome enabled up to 18.0% of genes from the rarest organism to be assembled, in contrast to the un‐normalized sample, where genes were not able to be assembled at the same sequencing depth. This study has demonstrated that the application of normalization methods to metagenomic samples is a powerful tool to enrich for sequences from rare taxa, which will shed further light on their ecological niches.     

Fingerprinting Diazotroph Communities in the Chesapeake Bay by Using a DNA Macroarray

Investigations of the distribution and diversity of nitrogen-fixing microorganisms in natural environments have often relied on PCR amplification and sequence analysis of a portion of one of the key enzymes in nitrogen fixation, dinitrogenase reductase, encoded by nifH. Recent work has suggested that DNA macroarrays provide semiquantitative fingerprints of diversity within mixtures of nifH amplicons (G. F. Steward, B. D. Jenkins, B. B. Ward, and J. P. Zehr, Appl. Environ. Microbiol. 70:1455-1465, 2004). Here we report the application of macroarrays for a study in the Chesapeake Bay. Samples from different locations in the bay yielded distinct fingerprints. Analysis of replicates and samples from different locations by cluster analysis showed that replicates clustered together, whereas different samples formed distinct clusters. There was a correspondence between the hybridization pattern observed and that predicted from the distribution of sequence types in a corresponding clone library. Some discrepancies between the methods were observed which are likely a result of the high nifH sequence diversity in the Chesapeake Bay and the limited number of sequences represented on this version of the array. Analyses of sequences in the clone library indicate that the Chesapeake Bay harbors unique, phylogenetically diverse diazotrophs. The macroarray hybridization patterns suggest that there are spatially variable communities of diazotrophs, which have been confirmed by quantitative PCR methods (S. M. Short, B. D. Jenkins, and J. P. Zehr, Appl. Environ. Microbiol., in press). The results show that DNA macroarrays have great potential for mapping the spatial and temporal variability of functional gene diversity in the environment.     

Seasonal Dynamics and Metagenomic Characterization of Estuarine Viriobenthos Assemblages by Randomly Amplified Polymorphic DNA PCR

Direct enumeration and genetic analyses indicate that aquatic sediments harbor abundant and diverse viral communities. Thus far, synecological analysis of estuarine sediment viral diversity over an annual cycle has not been reported. This oversight is due in large part to a lack of molecular genetic approaches for assessing viral diversity within a large collection of environmental samples. Here, randomly amplified polymorphic DNA PCR (RAPD-PCR) was used to examine viral genotypic diversity within Chesapeake Bay sediments. Using a single 10-mer oligonucleotide primer for all samples, RAPD-PCR analysis of sediment viral assemblages yielded unique banding patterns across spatial and temporal scales, with the occurrence of specific bands varying among the sample set. Cluster analysis of RAPD-PCR amplicon banding patterns indicated that sediment viral assemblages changed with season and to a lesser extent with geographic location. Sequence analysis of RAPD-PCR amplicons revealed that 76% of sediment viral sequences were not homologous to any sequence in the GenBank nonredundant protein database. Of the GenBank sequence homologs, the majority belonged to viruses within the Podoviridae (24%) and Myoviridae (22%) viral families, which agrees with the previously observed frequencies of these morphological families in Chesapeake Bay sediments. Furthermore, the majority of the sediment viral sequences homologous to GenBank nonredundant protein sequences were phages or prophages (57%). Hence, RAPD-PCR proved to be a reliable and useful approach for characterization of viral assemblages and the genetic diversity of viruses within aquatic sediments.Large numbers of viruses, an estimated abundance greater than 10³¹ viruses worldwide (11, 26), have been found in a variety of environments, including seawater (38), freshwater (19), sediments (25, 28), and soils (34). Viruses are not only abundant but also likely to significantly influence the population dynamics and genotypic composition of their bacterial host populations (29, 33). Process-level investigations of viral activity in sediments have shown that viruses are an active component of sediment microbial communities (23). Glud and Middelboe (23) found that bacterial growth rates and viral production increased in parallel with respiration, suggesting that viruses are active members of benthic microbial communities. Previous studies have shown that sediment viral abundance exceeds coexisting bacterial abundance by 10- to 1,000-fold (15, 17, 25), creating the potential for viral processes to influence the microbial ecology of aquatic sediments. However, with the exception of small-scale metagenomic investigations (4, 8), there exists little information on the genetic content of viriobenthos assemblages or how the composition of these assemblages changes over ecological gradients.Despite the high abundances of viruses in nature, the lack of a shared genetic marker creates a difficult problem when attempting to examine viral genetic diversity in environmental samples (31). Gene g20 encodes a multifunctional protein within the collar between the capsid and tail in T4-like bacteriophages and has been of significant importance in examining the genetic diversity of cyanomyoviruses (22, 24, 32). As well, others have been able to evaluate the diversity of unidentified aquatic picornavirus-like viruses using the RNA-dependent RNA polymerase gene (13). Other studies have attempted to examine phage genetic diversity based on the DNA polymerase gene (6, 21). Unfortunately, not all known phages contain these specific genes; hence, their use as universal markers is markedly inadequate. Thus, molecular methods that do not rely on polymorphism analysis of a single gene product must be used to circumvent these limitations.Recently, metagenomic approaches (i.e., sequencing of random genomic DNA fragments from whole microbial assemblages) have been used to examine genetic diversity within viral (18) and prokaryotic (10) assemblages. For sediment environments, metagenomic analysis has revealed that the viriobenthos is perhaps the most diverse of all viral assemblages, having been estimated to contain more than 10,000 genotypes per kg of sediment (4). Viral assemblages within a wide range of environments including marine (2, 8) and estuarine (3) waters, soils (20), stromatolites (16), and equine (9) and human feces (5, 40) have been examined. Overall, these studies have shown that a relatively low proportion (∼30%) of viral metagenome sequences are similar to sequences found in the nonredundant GenBank database (nr database), but the probability of detecting significant BLAST homologs increases twofold when queries against other viral metagenome sequence libraries are included (3). Thus, the function of most viral genes is currently unknown; however, these genes are broadly distributed among viruses.While large-scale metagenomics offers unprecedented resolution of the diversity and composition of a viral assemblage, the significant costs and computational requirements preclude routine application in a large collection of environmental samples. Recently, Winget and Wommack (36) introduced a new, low-cost, high-throughput means for genetic analysis of viral diversity utilizing random amplified polymorphic DNA PCR (RAPD-PCR). In this approach, a single 10-bp oligonucleotide serves as both the forward and reverse primers in a single thermocycler reaction. Target sequences in the template DNA are randomly selected; thus, development of a RAPD-PCR assay requires no prior information on the DNA coding content within the sample or organism—a significant advantage considering the largely unknown nature of most viral genes.In this study, we assess the potential of RAPD-PCR as a tool to examine genotype-scale compositional changes in the Chesapeake Bay viriobenthos and to explore the genetic diversity of viruses within Chesapeake Bay sediments. To our knowledge, this is the first study to use RAPD-PCR for evaluating sediment viral diversity and documenting compositional changes in viriobenthos assemblages over time and geographic location.     

Population Dynamics of Chesapeake Bay Virioplankton: Total-Community Analysis by Pulsed-Field Gel Electrophoresis

Recognition of viruses as the most abundant component of aquatic microbial communities has stimulated investigations of the impact of viruses on bacterio- and phytoplankton host communities. From results of field studies to date, it is concluded that in most aquatic environments, a reduction in the number of bacteria on a daily basis is caused by viral infection. However, the modest amount of in situ virus-mediated mortality may be less significant than viral infection serving to maintain clonal diversity in the host communities directly, through gene transmission (i.e., transduction), and indirectly, by elimination of numerically dominant host species. If the latter mechanism for controlling community diversity prevails, then the overall structure of aquatic viral communities would be expected to change as well over short seasonal and spatial scales. To determine whether this occurs, pulsed-field gel electrophoresis (PFGE) was used to monitor the population dynamics of Chesapeake Bay virioplankton for an annual cycle (1 year). Virioplankton in water samples collected at six stations along a transect running the length of the bay were concentrated 100-fold by ultrafiltration. Viruses were further concentrated by ultracentrifugation, and the concentrated samples were embedded in agarose. PFGE analysis of virus DNA in the agarose plugs yielded several distinct bands, ranging from 50 to 300 kb. Principal-component and cluster analyses of the virus PFGE fingerprints indicated that changes in virioplankton community structure were correlated with time, geographical location, and extent of water column stratification. From the results of this study, it is concluded that, based on the dynamic nature of the Chesapeake Bay virioplankton community structure, the clonal diversity of bacterio- and phytoplankton host communities is an important component of the virus community.     

Endospores of halophilic bacteria of the family Bacillaceae isolated from non-saline Japanese soil may be transported by Kosa event (Asian dust storm)

Background

Natural microbial communities are extremely complex and dynamic systems in terms of their population structure and functions. However, little is known about the in situ functions of the microbial communities.

Results

This study describes the application of proteomic approaches (metaproteomics) to observe expressed protein profiles of natural microbial communities (metaproteomes). The technique was validated using a constructed community and subsequently used to analyze Chesapeake Bay microbial community (0.2 to 3.0 μm) metaproteomes. Chesapeake Bay metaproteomes contained proteins from pI 4–8 with apparent molecular masses between 10–80 kDa. Replicated middle Bay metaproteomes shared ~92% of all detected spots, but only shared 30% and 70% of common protein spots with upper and lower Bay metaproteomes. MALDI-TOF analysis of highly expressed proteins produced no significant matches to known proteins. Three Chesapeake Bay proteins were tentatively identified by LC-MS/MS sequencing coupled with MS-BLAST searching. The proteins identified were of marine microbial origin and correlated with abundant Chesapeake Bay microbial lineages, Bacteroides and α-proteobacteria.

Conclusion

Our results represent the first metaproteomic study of aquatic microbial assemblages and demonstrate the potential of metaproteomic approaches to link metagenomic data, taxonomic diversity, functional diversity and biological processes in natural environments.     

Phylogenetic Diversity of Marine Cyanophage Isolates and Natural Virus Communities as Revealed by Sequences of Viral Capsid Assembly Protein Gene g20

In order to characterize the genetic diversity and phylogenetic affiliations of marine cyanophage isolates and natural cyanophage assemblages, oligonucleotide primers CPS1 and CPS8 were designed to specifically amplify ca. 592-bp fragments of the gene for viral capsid assembly protein g20. Phylogenetic analysis of isolated cyanophages revealed that the marine cyanophages were highly diverse yet more closely related to each other than to enteric coliphage T4. Genetically related marine cyanophage isolates were widely distributed without significant geographic segregation (i.e., no correlation between genetic variation and geographic distance). Cloning and sequencing analysis of six natural virus concentrates from estuarine and oligotrophic offshore environments revealed nine phylogenetic groups in a total of 114 different g20 homologs, with up to six clusters and 29 genotypes encountered in a single sample. The composition and structure of natural cyanophage communities in the estuary and open-ocean samples were different from each other, with unique phylogenetic clusters found for each environment. Changes in clonal diversity were also observed from the surface waters to the deep chlorophyll maximum layer in the open ocean. Only three clusters contained known cyanophage isolates, while the identities of the other six clusters remain unknown. Whether or not these unidentified groups are composed of bacteriophages that infect different Synechococcus groups or other closely related cyanobacteria remains to be determined. The high genetic diversity of marine cyanophage assemblages revealed by the g20 sequences suggests that marine viruses can potentially play important roles in regulating microbial genetic diversity.     

Randomly amplified polymorphic DNA PCR as a tool for assessment of marine viral richness

Recent discoveries have uncovered considerable genetic diversity among aquatic viruses and raised questions about the variability of this diversity within and between environments. Studies of the temporal and spatial dynamics of aquatic viral assemblages have been hindered by the lack of a common genetic marker among viruses for rapid diversity assessments. Randomly amplified polymorphic DNA (RAPD) PCR bypasses this obstacle by sampling at the genetic level without requiring viral isolation or previous sequence knowledge. In this study, the utility of RAPD-PCR for assessing DNA viral richness within Chesapeake Bay water samples was evaluated. RAPD-PCR using single 10-mer oligonucleotide primers successfully produced amplicons from a variety of viral samples, and banding patterns were highly reproducible, indicating that each band likely represents a single amplicon originating from viral template DNA. In agreement with observations from other community profiling techniques, resulting RAPD-PCR banding patterns revealed more temporal than spatial variability in Chesapeake Bay virioplankton assemblages. High-quality hybridization probes and sequence information were also easily generated from single RAPD-PCR products or whole reactions. Thus, the RAPD-PCR technique appears to be practical and efficient for routine use in high-resolution viral diversity studies by providing assemblage comparisons through fingerprinting, probing, or sequence information.     

Shotgun metagenomics indicates novel family A DNA polymerases predominate within marine virioplankton

Virioplankton have a significant role in marine ecosystems, yet we know little of the predominant biological characteristics of aquatic viruses that influence the flow of nutrients and energy through microbial communities. Family A DNA polymerases, critical to DNA replication and repair in prokaryotes, are found in many tailed bacteriophages. The essential role of DNA polymerase in viral replication makes it a useful target for connecting viral diversity with an important biological feature of viruses. Capturing the full diversity of this polymorphic gene by targeted approaches has been difficult; thus, full-length DNA polymerase genes were assembled out of virioplankton shotgun metagenomic sequence libraries (viromes). Within the viromes novel DNA polymerases were common and found in both double-stranded (ds) DNA and single-stranded (ss) DNA libraries. Finding DNA polymerase genes in ssDNA viral libraries was unexpected, as no such genes have been previously reported from ssDNA phage. Surprisingly, the most common virioplankton DNA polymerases were related to a siphovirus infecting an α-proteobacterial symbiont of a marine sponge and not the podoviral T7-like polymerases seen in many other studies. Amino acids predictive of catalytic efficiency and fidelity linked perfectly to the environmental clades, indicating that most DNA polymerase-carrying virioplankton utilize a lower efficiency, higher fidelity enzyme. Comparisons with previously reported, PCR-amplified DNA polymerase sequences indicated that the most common virioplankton metagenomic DNA polymerases formed a new group that included siphoviruses. These data indicate that slower-replicating, lytic or lysogenic phage populations rather than fast-replicating, highly lytic phages may predominate within the virioplankton.     

Amplification of uncultured single-stranded DNA viruses from rice paddy soil

Viruses are known to be the most numerous biological entities in soil; however, little is known about their diversity in this environment. In order to explore the genetic diversity of soil viruses, we isolated viruses by centrifugation and sequential filtration before performing a metagenomic investigation. We adopted multiple-displacement amplification (MDA), an isothermal whole-genome amplification method with phi29 polymerase and random hexamers, to amplify viral DNA and construct clone libraries for metagenome sequencing. By the MDA method, the diversity of both single-stranded DNA (ssDNA) viruses and double-stranded DNA viruses could be investigated at the same time. On the contrary, by eliminating the denaturing step in the MDA reaction, only ssDNA viral diversity could be explored selectively. Irrespective of the denaturing step, more than 60% of the soil metagenome sequences did not show significant hits (E-value criterion, 0.001) with previously reported viral sequences. Those hits that were considered to be significant were also distantly related to known ssDNA viruses (average amino acid similarity, approximately 34%). Phylogenetic analysis showed that replication-related proteins (which were the most frequently detected proteins) related to those of ssDNA viruses obtained from the metagenomic sequences were diverse and novel. Putative circular genome components of ssDNA viruses that are unrelated to known viruses were assembled from the metagenomic sequences. In conclusion, ssDNA viral diversity in soil is more complex than previously thought. Soil is therefore a rich pool of previously unknown ssDNA viruses.