Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2002, p. 1576-1584, Vol. 68, No. 4
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.4.1576-1584.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Department of Marine Sciences, University of Georgia, Athens, Georgia 30602,1 Center of Marine Biotechnology, University of Maryland Biotechnology Institute, Baltimore, Maryland 21202,2 Department of Microbiology, The University of Tennessee, Knoxville, Tennessee 379963
Received 21 November 2001/ Accepted 25 January 2002
|
|
|---|
|
|
|---|
Viruses also potentially influence the genetic diversity of microbial communities, mainly via virus selection pressure and virus-mediated genetic exchanges such as transduction, transformation, and lysogenic conversion (6, 9, 27). Cyanophage isolates are diverse in terms both of their morphology (20, 27) and their genetic fingerprints (11, 32). However, little is known about the genetic diversity and phyletic linkages of cyanophages in natural marine environments. Recently, the genes for viral capsid assembly g20 proteins from three marine cyanophages that infect marine Synechococcus spp. were sequenced and their sequences were similar to that of the g20-encoding gene of bacteriophage T4, a gene responsible for the initiation of head assembly, DNA packing, and head-tail junction (7). Several conserved regions were found to be unique among the g20 sequences of three marine cyanomyoviruses (7). A set of PCR primers (CPS1 and CPS2) was designed based on the conserved regions of these cyanophages and used to specifically amplify a 165-bp g20 fragment from cyanomyoviruses in natural seawater (7). This primer set was subsequently applied in PCR-DGGE (denaturing gradient gel electrophoresis) to study cyanophage diversity along a transect in the Atlantic Ocean (30, 31). High genetic diversity of cyanophage g20 sequences was found in both surface waters and the deep layer of the euphotic zone. Natural cyanophage community structure changed significantly from surface to depth, and maximum diversity was always coincident with maximum Synechococcus sp. abundance (31).
It is difficult to resolve the complex interactions between cyanophage and Synechococcus populations without knowing their phylogenetic distribution and phyletic variation in different environments. While natural marine Synechococcus populations have been well characterized phylogenetically (13, 14, 25), little is known about the phyletic relatedness and diversity within natural cyanophage communities in different marine environments. Our earlier study has shown that the short g20 gene sequence flanked by primers CPS1 and CPS2 (7) yielded less robust phylogenetic inference compared to that flanked by primers CPS1 and CPS8 (Y. Zhong, F. Chen, and R. E. Hodson, Abstr. 100th Gen. Meet. Am. Soc. Microbiol. abstr. 100, p. 465, 2000).
In this study, primers CPS1 and CPS8 were used to amplify the g20 gene fragment (ca. 592 bp) from cyanophage isolates and natural viral assemblages. The phylogenetic diversity of cyanophage communities in various marine environments was investigated based on the partial g20 gene sequences recovered from natural viral communities.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Specificity test of PCR amplification with primers CPS1 and CPS2
|
|
View this table: [in a new window] |
TABLE 2. Stations where VCs were collected
|
Clone library construction.
The PCR amplicon from each VC was purified with the Wizard PCR prep DNA purification system (Promega). The purified products were ligated into the pGEM-T Easy cloning vector (Promega) and then transformed into JM109 competent cells (Promega) in accordance with the manufacturer's instructions. Positive clones (white colonies) were picked randomly and transferred onto a new agar plate. Clones were numbered 1 to 80 and assigned the prefix SE (for Skidaway Estuarine), GS (for Gulf Stream), or SS (for Sargasso Sea).
Sequencing.
About 35 clones from each clone library were randomly picked, and the plasmid inserts were PCR amplified with vector-specific primers T7 and SP6 along with Expand High Fidelity DNA polymerase (Roche). PCR amplification involved a 3-min initial denaturation, followed by 35 cycles of 94°C for 15s, 45°C for 15s, 73°C for 1 min, and a final 4-min extension at 73°C. PCR products of the expected sizes were purified as described above and sequenced bidirectionally with primers T7 and SP6 using the BigDye Terminator Cycle Sequencing Kit (Applied Biosystems Inc.) on an ABI310 or ABI3700 automated DNA sequencer (Applied Biosystems Inc.). For cyanophage isolates, the PCR products were purified and sequenced directly on both strands with CPS1 and CPS8, respectively.
Phylogenetic analysis.
Pairs of sequences were aligned and corrected manually with Sequencher 3.0 software (Gene Codes Corporation, Ann Arbor, Mich.), and the consensus sequences were further analyzed. Sequence alignment and similarity analysis were performed with the Genetics Computer Group (Madison, Wis.) package and checked visually for chimeric artifacts. Phylogenetic reconstruction was accomplished with the phylogeny inference package (PHYLIP version 3.5 [5]). Evolutionary distances were calculated by the Jukes-Cantor method with the DNADIST program. Phylogenetic trees were constructed from distance matrices by the neighbor-joining method, which was implemented with the NEIGHBOR program. The node reproducibility for tree topology was estimated by bootstrap analysis, which included 100 replicate data sets.
Nucleotide sequence accession numbers.
The nucleotide sequences determined in this study have been submitted to GenBank and assigned accession numbers AY027938 to AY028078.
|
|
|---|
Phylogenetic diversity of cyanophage isolates.
The sequence alignment of the g20 gene fragments from 11 cyanophage isolates demonstrated that the region amplified by primers CPS1 and CPS8 was highly conserved and therefore suitable for inference of genetic relatedness among cyanophages (data not shown). Nucleotide similarity analysis revealed that cyanophage isolates shared 64.7 to 96.1% sequence similarity with each other, while enteric coliphage T4 was only 50.4 to 55.2% similar to any cyanophage isolates (Table 3). Based on the phylogenetic affiliation among cyanophage isolates and clones from six natural virus communities (shown later; refer to Fig. 3), these isolates were assigned to three clusters, with most of the isolates falling into clusters I and III and only one isolate (31B) falling into cluster II. S-BnM1 and P17 were distant from any isolates or clones; therefore, they were not grouped into any clusters (Fig. 1).
|
View this table: [in a new window] |
TABLE 3. Percent nucleotide sequence similarity between cyanophage isolates and phage T4
|
![]() View larger version (33K): [in a new window] |
FIG. 3. Neighbor-joining tree showing the phylogenetic affiliation of cyanophage isolates (boldface type) and representative clones from all six of the natural virus communities studied. The tree was constructed on the basis of a 176-amino-acid sequence alignment with T4 as the outgroup. Each value in parentheses is the number of different nucleotide sequences in the same cluster and same community as the representative clone. Clusters A through F and I through III were assigned on the basis of phylogenetic relatedness. Bootstrap values of less than 50 were not shown. The scale bar indicates 0.1 substitution per site.
|
![]() View larger version (13K): [in a new window] |
FIG. 1. Phylogenetic affiliation of cyanophage isolates based on 552-bp nucleotide sequences between primers CPS1 and CPS8. Coliphage T4 was used as an outgroup. Values at tree branches are bootstrap values with 100 replicates. The scale bar is equivalent to 0.1 replacement per site. Only bootstrap values larger than 50 are shown.
|
![]() ![]() View larger version (64K): [in a new window] |
FIG. 2. Phylogenetic relationships among the clones amplified with primers CPS1 and CPS8 from concentrated natural virus communities at Skidaway Estuarine surface SE1 (a), Gulf Stream surface GS26 (b), Gulf Stream DCM GS27 (c), Sargasso Sea surface SS48 (d), Sargasso Sea DCM SS47 (e), and Sargasso Sea DCM SS40 (f). The neighbor-joining tree was constructed based on 552-bp nucleotide sequence alignment with T4 as the outgroup. Only bootstrap values of greater than 50 are shown. Each value in parentheses is the number of clones bearing identical nucleotide sequences in a community. The clusters were assigned on the basis of Fig. 3.
|
(ii) Gulf Stream surface sample (GS26).
The Gulf Stream surface sample (GS26) community exhibited much less diversity in the cyanophage g20 gene pool than the Skidaway Estuarine surface community did. Among the 37 clones sequenced, 13 different g20 sequences were found and assigned to three clusters. GS2608 and GS2624 were in clusters III and B, respectively, while all of the other sequences fell into cluster E with 95.8 to 99.8% pairwise nucleotide similarity (Fig. 2b).
(iii) Gulf Stream DCM sample (GS27).
The cyanophage community at the DCM of the Gulf Stream appeared to be more diverse than that at the surface (GS26). Among the 33 clones analyzed, 23 different g20 gene sequences were encountered. The DCM layer contained the same three clusters as did the surface sample (clusters III, B, and E) plus three other clusters (C, D, and II) that were not present on the surface (Fig. 2c). However, unlike the surface sample (GS26), where cluster E was dominant, clusters B and E were the two major clusters at the DCM of the Gulf Stream.
(iv) Sargasso Sea surface sample (SS48).
The g20 gene diversity and phylogenetic structure of the Sargasso Sea surface sample (SS48) were similar to those of the surface sample in the Gulf Stream (GS26). Fifteen different g20 sequences were found in the 35 clones analyzed. Of the three phyletic clusters in this community, clusters B and E were both present and E was the dominant cluster (Fig. 2d). The higher number of repeat sequences and fewer phyletic groups in GS26 and SS48 indicated a less diversified g20 gene pool in these oceanic surface waters.
(v) Sargasso Sea DCM sample (SS47 and SS40).
Two samples in the DCM layer of the Sargasso Sea were analyzed and exhibited similar phylogenetic patterns. Among the 36 and 31 clones studied in samples SS47 and SS40, 24 and 26 different g20 sequences were identified, respectively. Four common phyletic clusters (II, C, E, and B) were shared and the tree topologies within clusters II, E, and B were very similar between the two communities (Fig. 2e and f). Compared to Sargasso Sea surface sample SS48, the two DCM samples (SS47 and SS40) both contained more clusters. The dominant cluster shifted from cluster E in the surface to cluster B in the DCM layer, as was observed in the Gulf Stream. However, the two DCM samples in the Sargasso Sea were also differed from each other. For example, clone SS4716, found in the SS47 sample, did not cluster with any other clones or cyanophage isolates whereas SS40 had two more clusters than SS47.
Comprehensive results.
Among 207 g20 clones randomly picked from six natural viral communities, 114 different g20 homologs were encountered. No clone had a nucleotide sequence identical to that of any of the cyanophage isolates in this study, indicating that cyanophages are very diverse in natural environments. The distribution of clones and g20 sequences within each community is summarized in Table 4. Nine distinct clusters are revealed in the comprehensive phylogenetic tree in Fig. 3 based on the amino acid sequence alignment of g20. Clusters I, II, and III comprised all of the 11 known isolates and 32 different natural clones (nucleotide sequences). Cluster I included five isolates and 15 different sequences that were exclusively from the Skidaway Estuarine sample. The pairwise sequence similarity levels in cluster I ranged from 75.2 to 98.0%. Cluster II included an isolate from the Sargasso Sea and eight different sequences from three oceanic DCM communities in the Gulf Stream and the Sargasso Sea. The sequence identity levels in this cluster ranged from 73.2 to 99.5%. Cluster III consisted of six isolates and nine different clones from the Skidaway Estuary, the Gulf Stream, and the Sargasso Sea, and there was 67.0 to 98.8% sequence similarity in this cluster.
|
View this table: [in a new window] |
TABLE 4. Distribution of g20 gene clones and nucleotide sequences in the six virus communities studied
|
|
|
|---|
Phylogenetic analysis with isolated cyanophages indicated that marine cyanophages were genetically divergent yet more closely related to each other than to bacteriophage T4 (Fig. 1). Similar results were reported for marine microalgal viruses, which formed a unique group distant from other double-stranded DNA viruses infecting eukaryotes (2). Similarity analysis of g20 sequences further suggested that no correlation was present between genetic variation and geographic distance (Tables 1 and 3). Isolates from the same or adjacent water bodies could share a high level of g20 identity, e.g., 90.4% between 27A and 32A and 96.1% between P79 and P81. However, isolates from different water bodies could also share a high level of similarity, e.g., 96.1% between S-WHM1 and P77. On the other hand, cyanophage isolates from the same region can be as diverse as those from different oceans. For example, there is 66.2% similarity between 31B and 44B and 68.3% similarity between 31B and P17. Our results obtained with cyanophage isolates support the suggestion that genetically related marine viruses are widely distributed in the ocean without significant geographic segregation (10).
Natural cyanophage communities from both estuarine and oceanic environments exhibited even greater genetic diversity than did the cyanophage isolates we studied. Among the nine phylogenetic clusters identified from the g20 sequences, three clusters contained known cyanophage isolates while the other six clusters contained only natural clones. Even within a single cluster (e.g., cluster I, II, or III), clones demonstrated their divergent nature by the wide range of sequence similarity (73.2 to 99.5%). Furthermore, except for the two surface oceanic samples, all of the other communities contained five or six distinct genetic groups, with up to 29 different genotypes found in a single sample. The cyanophage diversity reported here is much greater than that reported from previous studies (11, 31). Wilson et al. (30) used g20 gene-based cyanomyovirus-specific primers and DGGE to investigate cyanophage diversity in the Atlantic Ocean. Only between 2 and 10 different genotypes were found by DGGE analysis in each concentrated virus community. The difference between our results and those of Wilson et al. (31) may result from the much higher sensitivity of our cloning-sequencing analysis and the degeneracy in our primers that was not included in the DGGE-oriented PCR primers. Many of the cyanophages used in this study were isolated from Georgia estuaries and exhibited diverse transmission electron microscopy morphology and restriction fragment length polymorphism patterns (11). Our results not only revealed much greater diversity among the cyanophages in marine environments but also demonstrated their phylogenetic affiliation via nucleic acid sequence analysis.
High phage diversity in natural marine environments could be due to either phage-host genetic exchange through transduction or genetic exchange between coinfecting phages. The g20 gene is within a mobile genetic module that is exchangeable between viruses during coinfection of bacterial hosts (12). The diverse cyanophage assemblages and dynamic phage-host interactions observed suggest that marine phages may be important in shaping the genetic diversity and composition of marine microbial communities.
It is likely that clusters I, II, and III represent lytic phages that infect WH7803-like Synechococcus spp. The wide geographical origins of isolates in clusters I and III further suggest that cyanophages in these clusters may infect Synechococcus spp. that are adapted to widely differing marine environments, whereas cluster II includes cyanophages whose hosts are adapted solely to oligotrophic environments. The identities of clusters A through F were enigmatic due to the absence of any isolates in them. Possibly, some clusters may represent cyanophages infecting Synechococcus cells that are genetically distant from the WH7803-like strains, from which most of our cyanophages were isolated. It is known that considerable genetic diversity exists within natural marine Synechococcus assemblages (13, 25). Isolation of more cyanophages that infect different marine Synechococcus spp. will shed light on the origins of the unknown clusters. We also cannot rule out the possibility that some g20 clones in clusters A to F are other bacteriophages whose hosts are closely related to Synechococcus spp. For example, marine Prochlorococcus spp. are abundant (4, 16) and genetically similar to Synechococcus spp. in terms of 16S rRNA (26) and RNA polymerase gene phylogeny (15). It is possible that our primers may amplify Prochlorococcus phages. This hypothesis remains to be tested when the phages infecting Prochlorococcus spp. become available.
The genetic composition and structure of natural cyanophage communities in estuaries and the open ocean were different from each other. While 56% of the estuarine g20 clones studied were within clusters I, II, and III, the vast majority (87% average) of the oceanic clones fell into clusters B to F. Moreover, unique phylogenetic clusters were found in each environment. For example, clusters A, F, and I were amplified from the estuarine sample whereas clusters B, C, D, and II were only present in the oceanic samples (Table 4). We also observed consistently higher cyanophage diversity in the oceanic DCM samples than in the surface samples, even though different waters were studied (i.e., the Gulf Stream versus the Sargasso Sea). Finally, it is evident that the dominant g20 clusters shifted from cluster E at the oceanic surface to cluster B at the DCM (Table 4). The distinct cyanophage population structures in estuarine water versus offshore water and in surface water versus the DCM layer suggested the presence of different host populations in response to different light, nutrient, and other physical conditions and a dynamic interaction between cyanophage and cyanobacterial populations. Little is known about the influence of environmental conditions on the genetic structure of natural Synechococcus assemblages. Presumably, the increasing cyanophage diversity from the surface to the DCM layer is related to the more diverse host populations at depth in stratified water columns.
The gene for viral capsid assembly protein g20 was successfully used in this study as a marker molecule to investigate genetic diversity and phylogenetic affiliation among cyanophage isolates and natural cyanophage assemblages in estuarine and oligotrophic environments. Phylogenetic analysis of cyanophage g20 gene sequences revealed that natural cyanophage populations are strikingly diverse and that their genetic structures vary greatly in different marine environments. This study suggests that marine viruses are under significant selection pressure and are able to maintain their diversity in response to host resistance. Further studies on the covariation of the genetic diversity of virus and host populations with changing environmental variables would provide new insights into the ecological roles of marine viruses.
We thank Jingrang Lu at the University of Georgia for providing the cyanophage isolates from Georgia estuaries. We also thank Ju Sheng for testing PCR primers on other bacteriophages isolated from the Baltimore Inner Harbor. Curtis Suttle at the University of British Columbia, John Waterbury at Woods Hole Oceanographic Institution, and William Wilson at the Plymouth Marine Laboratory, Plymouth, United Kingdom, kindly provided cyanophage isolates for this study.
Contribution no. 568 from the Center of Marine Biotechnology, University of Maryland Biotechnology Institute. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»