This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lim-Fong, G. E.
Right arrow Articles by Haygood, M. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lim-Fong, G. E.
Right arrow Articles by Haygood, M. G.
Agricola
Right arrow Articles by Lim-Fong, G. E.
Right arrow Articles by Haygood, M. G.

 Previous Article  |  Next Article 

Applied and Environmental Microbiology, June 2008, p. 3605-3609, Vol. 74, No. 11
0099-2240/08/$08.00+0     doi:10.1128/AEM.02798-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Evolutionary Relationships of "Candidatus Endobugula" Bacterial Symbionts and Their Bugula Bryozoan Hosts{triangledown}

Grace E. Lim-Fong,1,3 Lindsay A. Regali,3 and Margo G. Haygood2*

Scripps Institution of Oceanography, Marine Biology Research Division, University of California, San Diego, La Jolla, California 92093-0202,1 Department of Environmental and Biomolecular Systems, OGI School of Science & Engineering, Oregon Health & Science University, Beaverton, Oregon 97006-8921,2 Department of Biology, Randolph-Macon College, Ashland, Virginia 230053

Received 11 December 2007/ Accepted 24 March 2008


arrow
ABSTRACT
 
Ribosomal gene sequences were obtained from bryozoans in the genus Bugula and their bacterial symbionts; analyses of host and symbiont phylogenetic trees did not support a history of strict cospeciation. Symbiont-derived compounds known to defend host larvae from predation were only detected in two out of four symbiotic Bugula species.


arrow
INTRODUCTION
 
Like several other marine invertebrates (7, 8), some bryozoans harbor proteobacterial symbionts (12, 15, 28, 29). Woollacott (28) observed bacteria in the pallial sinus of larvae of the bryozoans Bugula neritina, Bugula simplex, and Bugula pacifica under the electron microscope, whereas bacteria were absent in Bugula turrita and Bugula stolonifera larvae. B. neritina was subsequently found to be a complex of at least three sibling species, Shallow (S), Deep (D) and North Atlantic (N) (4, 18). The D and S sibling species of B. neritina contain a bacterial symbiont called "Candidatus Endobugula sertula," the producer of bryostatins (3, 17), a class of defensive compounds which was shown to deter fish from eating B. neritina larvae in the S sibling species (16, 17). Bryostatins also have anticancer activity (20). Thus, "Ca. Endobugula sertula" confers chemical defense as a benefit to the bryozoan host (17). Neither "Ca. Endobugula sertula" nor bryostatins were detected in the N sibling species of B. neritina (17, 18). Bugula simplex also possesses a bacterial symbiont, "Candidatus Endobugula glebosa," that is closely related to "Ca. Endobugula sertula." This symbiont, which resides in the same location of the larva, is also associated with bryostatin activity (15).

The presence of the bacterial symbionts in larval tissue and their absence in the surrounding seawater are indicative of a vertical mode of transmission (26). Moreover, sympatric Bugula species do not appear to exchange symbionts (28). In general, vertically transmitted symbionts tend to cospeciate with their hosts, as evidenced by the aphid Buchnera (2) and ant "Candidatus Blochmannia" (23) examples. Topological congruence between the host and symbiont phylogenetic trees is usually used as evidence for cospeciation (9, 11). This study aimed to determine whether cospeciation could account for the pattern of symbiont distribution in the genus Bugula and whether the symbionts produce bryostatin-like molecules for host defense throughout the genus. If so, the symbiotic Bugula species should form a monophyletic group whose topology is congruent with the tree of "Ca. Endobugula" symbionts, and bryostatin activity should be present only in symbiotic Bugula species.


arrow
Sample collection and DNA extraction.
 
Bryozoan colonies were either collected by hand off floating docks or via scuba dive and were rinsed in sterile seawater. DNA was extracted from the samples immediately after collection, or the samples were preserved in 100% ethanol for the bryostatin assays. Genomic DNA, which comprises host and symbiont DNA, was extracted from the upper portion of each colony using the DNeasy DNA extraction kit (Qiagen, Inc., Valencia, CA) by following the manufacturer's directions. DNA from B. neritina was separated from low-molecular-weight PCR inhibitors on agarose gels and was purified using a rapid gel extraction system kit (Marligen Biosciences, Inc., Ijamsville, MD) prior to PCR.


arrow
Symbiont 16S SSU rRNA gene sequencing and analysis.
 
We used a combination of universal and specific primers (Table 1) to amplify the symbiont 16S small-subunit (SSU) rRNA gene from Bugula dentata, B. pacifica, B. stolonifera, Bugula turbinata, and B. turrita. Only the B. pacifica and B. turbinata samples resulted in a positive amplification (Table 2), indicating that of the five Bugula species tested, only those two contained "Ca. Endobugula"-related symbionts. These results are consistent with Woollacott's electron microscopy results for the species he examined.


View this table:
[in this window]
[in a new window]

  		TABLE 1. List of primers used in this studya


View this table:
[in this window]
[in a new window]

  		TABLE 2. Gene sequences obtained in this study and deposited in GenBank

These symbiont 16S SSU rRNA sequences were edited in Sequencher 4.7 (Genecodes Corp., Ann Arbor, MI) and aligned by eye with sequences from "Ca. Endobugula sertula" and "Ca. Endobugula glebosa" (symbionts of B. neritina and B. simplex, respectively) and other gammaproteobacteria using secondary-structure information (1), and a phylogenetic tree was constructed using 1,376 characters with the alphaproteobacterium Rhodospirillum rubrum as an out-group.

Maximum parsimony (MP), neighbor-joining (NJ), and maximum likelihood (ML) trees were constructed in PAUP* version 4.0b (D. L. Swofford, Sinauer Associates, Sunderland, MA). Transversions (Tv) were weighted twice the transitions (Ti) in parsimony analyses based on the Ti:Tv ratio estimated by ML, and the evolution model used in NJ and ML analyses was chosen by restricting the general time-reversible (GTR) gamma I model using the likelihood ratios test. Nodal support of the tree topologies was calculated by bootstrapping the MP and NJ trees with 1,000 replicates and by Bayesian analysis (22). Bayesian posterior probabilities were obtained by MrBayes in two ways: (i) by applying a 4 x 4 model of DNA substitution to all characters and (ii) by applying a doublet model to stem regions of the SSU 16S rRNA molecule. Default settings were used except for the rate prior, which was set to variable. One million generations were run with sampling at every 100th generation; the first 1,000 Bayesian trees were discarded as burn-in. All analyses, regardless of the phylogenetic method, demonstrate that the "Ca. Endobugula" symbionts form a well-supported clade among the gammaproteobacteria and are closely related to the cellulolytic symbionts of the shipworm mollusk Lyrodus pedicellatus (Fig. 1). We used RRTree (21) to compare the G+C contents and the evolutionary rates of the 16S rRNA genes of "Ca. Endobugula" symbionts to those of closely related free-living bacteria and found no significant difference for either parameter (P = 0.44 and P = 0.12, respectively).


Figure 1
View larger version (16K):
[in this window]
[in a new window]

  		FIG. 1. Phylogenetic tree of "Ca. Endobugula" symbionts and closely related bacteria based on the 16S SSU rRNA gene with R. rubrum as an out-group estimated using MP analysis. Support for internal nodes is shown as follows (top to bottom): MP bootstrap, NJ bootstrap, Bayesian posterior probability for a 4 x 4 model, Bayesian posterior probability, and a 4 x 4 model for loops/doublet model for stems. Only values above 50% or 0.5 are shown. The scale bar represents 10 nucleotide changes.


arrow
Bryozoan host 16S LSU rRNA and COI gene sequencing and analysis.
 
The bryozoan mitochondrial 16S large-subunit (LSU) rRNA and cytochrome oxidase I (COI) genes were amplified with universal primers (Table 2) as previously described (5, 10) for all the collected bryozoans. The N sibling species was not tested, as specimens were not available. B. dentata did not yield a clean COI PCR product, so a B. dentata COI sequence from GenBank was used instead. The 16S LSU rRNA (409 bp) and the COI genes (507 bp) were used to determine the phylogeny of the host bryozoans, and these genes were chosen as they are routinely used to resolve bryozoan phylogenies at the species level (6, 14). The incongruent-length difference test (9) demonstrated that the two genes were not incongruent (P = 0.56) and can thus be concatenated for data analysis. The likelihood criterion estimated Ti:Tv ratios of 2:1 and 3:1 for the 16S LSU rRNA and COI genes, respectively, and these estimates were used in MP reconstructions of the concatenated data set. Bayesian posterior probabilities were obtained in a manner similar to that described for the symbiont phylogeny with the following exceptions: a 4 x 4 DNA substitution model was applied to the 16S LSU rRNA partition, and a codon model was applied to the COI partition. Both MP and Bayesian analyses indicate that the symbiotic Bugula species do not form a monophyletic clade (Fig. 2). This observation was also supported by separate analyses of the two genes (data not shown). Using the Shimodaira-Hasegawa test (25) in PAUP*, we found that constraining the symbiotic Bugula into a monophyletic group could not explain the observed data as well as the consensus tree (P = 0.003). Even though we can rule out the possibility of strict cospeciation (without host switching or symbiont loss) explaining the distribution of bacterial symbionts in Bugula, other evolutionary scenarios should be explored.


Figure 2
View larger version (14K):
[in this window]
[in a new window]

  		FIG. 2. Incongruent Bayesian host and symbiont phylogenetic trees based on mitochondrial COI and LSU 16S rRNA genes and bacterial 16S rRNA genes, respectively. Support for internal nodes is shown on the host tree as follows (top to bottom): MP bootstrap, Bayesian posterior probability for the combined data set, Bayesian posterior probability for the COI data set, and Bayesian posterior probability for the16S LSU rRNA data set. Support for internal nodes is shown on the symbiont tree as follows (top to bottom): MP bootstrap, NJ bootstrap, and Bayesian posterior probability for a 4 x 4 model. Only values above 50% or 0.5 are shown. Bryostatin activity was detected in bryozoans whose names are in a bold font. The scale bar represents 0.1 substitutions per site.


arrow
Cospeciation analyses.
 
One possible scenario is the early acquisition of a symbiont followed by a subsequent loss of symbionts in some taxa. In this case, the symbiotic Bugula species tree should be congruent with the "Ca. Endobugula" tree. This does not appear to be the case (Fig. 2). We performed parametric bootstrapping (13) and found that the difference in host and symbiont topologies is not due to sampling error (P < 0.0001), meaning that the host and symbiont truly have different evolutionary histories. We then reconstructed a possible cospeciation scenario in TreeMap 2.0 (19) by comparing the symbiotic host topology to that of its symbionts. The most parsimonious reconstruction (weighting cospeciations:host switching:symbiont loss as 0:1:1) required three cospeciation events and one host switch. One thousand random symbiont trees were generated, and we found that it was possible to arrive at the same reconstruction by chance (P = 0.11). The TreeMap analysis assumes that all nodes are strongly supported, and this is not the case for the node leading to B. pacifica. We thus swapped the positions of the B. neritina cluster and B. pacifica on the tree and performed the Shimodaira-Hasegawa test (25) to determine if the observed data could be explained equally well by this rearranged tree. The swapped tree was equally supported by the data (P = 0.10), indicating that the ambiguous placement of B. pacifica will make it hard to determine if the similarity between the host and symbiont topologies arose by chance.


arrow
Bryostatin activity.
 
The bryostatin activities of various Bugula species were tested using the phorbol dibutyrate (PdBU) displacement assay, which tests for the ability of a compound, like a bryostatin, to displace PdBU from protein kinase C (24). Previous research has demonstrated activity in extracts of B. neritina and B. simplex (15); these extracts were used as positive controls. The assay was performed on adult extracts of B. pacifica, B. turbinata, B. turrita, B. stolonifera, and B. dentata, using the protocol described previously (15), and none of the tested extracts showed any bryostatin activity. One drawback of the PdBU displacement assay is that polyketides structurally and biosynthetically related to bryostatins but which possess altered functional groups will not be detected, and these polyketides may still be able to function as defensive metabolites. These data suggest that B. pacifica, B. turbinata, B. turrita, B. stolonifera, and B. dentata may lack symbiotic chemical defense; however, direct feeding deterrence assays are required to confirm this inference.

In this study, we obtained symbiont 16S SSU rRNA gene sequences from B. pacifica and B. turbinata, and these sequences form a monophyletic group with the corresponding gene in "Ca. Endobugula sertula" (B. neritina symbiont) and "Ca. Endobugula glebosa" (B. simplex symbiont). Though the symbionts are very closely related, we do not currently have sufficient evidence supporting cospeciation. Alternative hypotheses include one in which "Ca. Endobugula" colonized their respective hosts independently or one in which rampant host switching masks a few cospeciation events. Some nodes of the host tree are poorly supported, and perhaps the use of additional genetic markers will clarify the position of B. pacifica on the tree; these additional data should make cospeciation analyses more robust. The absence of bryostatin activity detectable by the PdBU displacement assay in symbiotic B. pacifica and B. turbinata is surprising. One possible scenario is that the bryostatin genes were present in the ancestral "Ca. Endobugula" species and were only retained in B. neritina and B. simplex, as these bryozoans have large and brightly colored larvae that need chemical protection (14). We postulate that "Ca. Endobugula" species might enhance the fitness of their hosts by perhaps playing either nutritional (26) or developmental (27) roles, which have been documented in other marine invertebrates. The addition of more species to the host tree, and the investigation of the presence of bacterial symbionts of these host species, will shed more light on the evolutionary history of the Bugula-"Ca. Endobugula" symbiosis and on the alternative roles that these symbionts might play in their host bryozoans.


arrow
Nucleotide sequence accession numbers.
 
The nucleotide sequence accession numbers for the gene sequences obtained in this study and deposited in GenBank are listed in Table 2.


arrow
ACKNOWLEDGMENTS
 
G.E.L.-F. was a Howard Hughes Medical Institute predoctoral fellow, and L.A.R. was a Randolph-Macon College Schapiro undergraduate research fellow.

This study was supported by a California Sea Grant (grant R/MP-88) and by the National Institutes of Health (grant 5R01CA079678-03).

We are grateful to Ray Andersen and Scottie Henderson for help with sample collection, to John Huelsenbeck for advice on phylogenetic analyses, and to Matthew Dick for his feedback on the manuscript.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Department of Environmental and Biomolecular Systems, OGI School of Science & Engineering, Oregon Health & Science University, Beaverton, OR 97006-8921. Phone: (503) 748-1993. Fax: (503) 748-1464. E-mail: haygood{at}ebs.ogi.edu Back

{triangledown} Published ahead of print on 4 April 2008. Back


arrow
REFERENCES
 
    1
  1. Cannone, J. J., S. Subramanian, M. N. Schnare, J. R. Collett, L. M. D'Souza, Y. Du, B. Feng, N. Lin, L. V. Madabusi, K. M. Müller, N. Pande, Z. Shang, N. Yu, and R. R. Gutell. 2002. The comparative RNA web (CRW) site: an online database of comparative sequence and structure information for ribosomal, intron, and other RNAs. BMC Bioinformatics 3:2.[CrossRef][Medline]
    2. 2
  2. Clark, M. A., N. A. Moran, P. Baumann, and J. J. Wernegreen. 2000. Cospeciation between bacterial endosymbionts (Buchnera) and a recent radiation of aphids (Uroleucon) and pitfalls of testing for phylogenetic congruence. Evolution Int. J. Org. Evolution 54:517-525.[CrossRef][Medline]
    3. 3
  3. Davidson, S. K., S. W. Allen, G. E. Lim, C. M. Anderson, and M. G. Haygood. 2001. Evidence for the biosynthesis of bryostatins by the bacterial symbiont "Candidatus Endobugula sertula" of the bryozoan Bugula neritina. Appl. Environ. Microbiol. 67:4531-4537.[Abstract/Free Full Text]
    4. 4
  4. Davidson, S. K., and M. G. Haygood. 1999. Identification of sibling species of the bryozoan Bugula neritina that produce different anticancer bryostatins and harbor distinct strains of the bacterial symbiont "Candidatus Endobugula sertula." Biol. Bull. 196:273-280.[Abstract]
    5. 5
  5. Dick, M. H., J. R. Freeland, L. P. Williams, and M. Coggeshall-Burr. 2000. Use of 16S mitochondrial ribosomal DNA sequences to investigate sister-group relationships among gymnolaemate bryozoans, p. 197-210. In A. Herrera Cubilla and J. B. C. Jackson (ed.), Proceedings of the 11th International Bryozoology Association Conference. International Bryozoological Association, Balboa, Republic of Panama.
    6. 6
  6. Dick, M. H., A. Herrera-Cubilla, and J. B. Jackson. 2003. Molecular phylogeny and phylogeography of free-living Bryozoa (Cupuladriidae) from both sides of the Isthmus of Panama. Mol. Phylogenet. Evol. 27:355-371.[CrossRef][Medline]
    7. 7
  7. Distel, D. L., E. F. DeLong, and J. B. Waterbury. 1991. Phylogenetic characterization and in situ localization of the bacterial symbiont of shipworms (Teredinidae: Bivalvia) by using 16S rRNA sequence analysis and oligodeoxynucleotide probe hybridization. Appl. Environ. Microbiol. 57:2376-2382.[Abstract/Free Full Text]
    8. 8
  8. Eisen, J. A., S. W. Smith, and C. M. Cavanaugh. 1992. Phylogenetic relationships of chemoautotrophic bacterial symbionts of Solemya velum Say (Mollusca: Bivalvia) determined by 16S rRNA gene sequence analysis. J. Bacteriol. 174:3416-3421.[Abstract/Free Full Text]
    9. 9
  9. Farris, J. D., M. Kallersjo, A. G. Kluge, and C. Bult. 1995. Constructing a significance test for incongruence. Syst. Biol. 44:570-572.[Free Full Text]
    10. 10
  10. Folmer, O., M. Black, W. Hoeh, R. A. Lutz, and R. C. Vrijenhoek. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 3:294-299.[Medline]
    11. 11
  11. Hafner, M. S., and S. A. Nadler. 1988. Phylogenetic trees support the coevolution of parasites and their hosts. Nature 332:258-259.[CrossRef][Medline]
    12. 12
  12. Haygood, M. G., and S. K. Davidson. 1997. Small subunit ribosomal RNA genes and in situ hybridization of the bacterial symbionts in the larvae of the bryozoan Bugula neritina and proposal of "Candidatus Endobugula sertula." Appl. Environ. Microbiol. 63:4612-4616.[Abstract]
    13. 13
  13. Huelsenbeck, J. P., B. Rannala, and Z. Yang. 1997. Statistical tests of host-parasite cospeciation. Evolution Int. J. Org. Evolution 51:410-419.[CrossRef]
    14. 14
  14. Hughes, R. N., A. Gómez, P. J. Wright, H. I. Moyano, J. M. Cancino, G. R. Carvalho, and D. H. Lunt. 2008. Molecular phylogeny supports division of the ‘cosmopolitan’ taxon Celleporella (Bryozoa; Cheilostomata) into four major clades. Mol. Phylogenet. Evol. 46:369-374.[CrossRef][Medline]
    15. 15
  15. Lim, G. E., and M. G. Haygood. 2004. "Candidatus Endobugula glebosa," a specific bacterial symbiont of the marine bryozoan Bugula simplex. Appl. Environ. Microbiol. 70:4921-4929.[Abstract/Free Full Text]
    16. 16
  16. Lindquist, N., and M. E. Hay. 1996. Palatability and chemical defense of marine invertebrate larvae. Ecol. Monogr. 66:431-450.[CrossRef]
    17. 17
  17. Lopanik, N. B., N. Lindquist, and N. Targett. 2004. Potent cytotoxins produced by a microbial symbiont protect host larvae from predation. Oecologia 139:131-139.[CrossRef][Medline]
    18. 18
  18. McGovern, T. M., and M. E. Hellberg. 2003. Cryptic species, cryptic endosymbionts, and geographical variation in chemical defences in the bryozoan Bugula neritina. Mol. Ecol. 12:1207-1215.[CrossRef][Medline]
    19. 19
  19. Page, R. D. M., and M. A. Charleston. 1998. Trees within trees: phylogeny and historical associations. Trends Ecol. Evol. 13:356-359.[CrossRef]
    20. 20
  20. Pettit, G. R. 1991. The bryostatins. Fortschr. Chem. Org. Naturst. 57:153-195.[Medline]
    21. 21
  21. Robinson-Rechavi, M., and D. Huchon. 2000. RRTree: relative-rate tests between groups of sequences on a phylogenetic tree. Bioinformatics 16:296-297.[Abstract/Free Full Text]
    22. 22
  22. Ronquist, F., and J. P. Huelsenbeck. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19:1572-1574.[Abstract/Free Full Text]
    23. 23
  23. Sauer, C., E. Stackebrandt, J. Gadau, B. Holldobler, and R. Gross. 2000. Systematic relationships and cospeciation of bacterial endosymbionts and their carpenter ant host species: proposal of the new taxon Candidatus Blochmannia gen. nov. Int. J. Syst. Evol. Microbiol. 50:1877-1886.[Abstract]
    24. 24
  24. Schaufelberger, D. E., B. Alvarado, P. Andrews, and J. A. Beutler. 1990. Detection and quantitation of bryostatin 1 and 2 in Bugula neritina by combined high-performance liquid chromatography and 3H-phorbol dibutyrate displacement. J. Liq. Chromatogr. 13:583-598.[CrossRef]
    25. 25
  25. Shimodaira, H., and M. Hasegawa. 1999. Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol. Biol. Evol. 16:1114-1116.
    26. 26
  26. Sipe, A. R., A. E. Wilbur, and S. C. Cary. 2000. Bacterial symbiont transmission in the wood-boring shipworm Bankia setacea (Bivalvia: Teredinidae). Appl. Environ. Microbiol. 66:1685-1691.[Abstract/Free Full Text]
    27. 27
  27. Visick, K. L., and M. J. McFall-Ngai. 2000. An exclusive contract: specificity in the Vibrio fischeri-Euprymna scolopes partnership. J. Bacteriol. 182:1779-1787.[Free Full Text]
    28. 28
  28. Woollacott, R. M. 1981. Association of bacteria with bryozoan larvae. Mar. Biol. 65:155-158.[CrossRef]
    29. 29
  29. Zimmer, R., and R. M. Woollacott. 1983. Mycoplasma-like organisms: occurrence with the larvae and adults of a marine bryozoan. Science 220:208-210.[Abstract/Free Full Text]

Applied and Environmental Microbiology, June 2008, p. 3605-3609, Vol. 74, No. 11
0099-2240/08/$08.00+0     doi:10.1128/AEM.02798-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lim-Fong, G. E.
Right arrow Articles by Haygood, M. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lim-Fong, G. E.
Right arrow Articles by Haygood, M. G.
Agricola
Right arrow Articles by Lim-Fong, G. E.
Right arrow Articles by Haygood, M. G.

Copyright © 2009 by the American Society for Microbiology.   For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»