Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2009, p. 3513-3521, Vol. 75, No. 11
0099-2240/09/$08.00+0 doi:10.1128/AEM.00002-09
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

KAUST Global Academic Partnership Program, Department of Biology, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong SAR, China
Received 1 January 2009/ Accepted 31 March 2009
|
|
|---|
|
|
|---|
There are generally two pathways by which sponges may acquire their bacterial associates. The first pathway is filter feeding and selective retention of bacteria (44). Bacteria in the surrounding seawater can be captured by sponges when sponges filter the food particles out of the water column. The bacteria that resist digestion by sponge choanocytes and archaeocytes can survive and live inside the sponges. Sponge-associated bacterial communities acquired via this pathway are therefore heavily influenced by the type of bacteria in the water column. The second pathway is vertical transmission of bacterial associates from adult sponges to their progeny (10, 40, 49). Bacteria acquired through this pathway may exhibit specificity for genera or species of sponges due to coevolution (for a review, see reference 44).
Sponge-bacterium symbioses are often mutualistic; while bacteria may benefit from the favorable nutritional conditions in sponges (18, 44), some bacterial associates may help their hosts eliminate metabolic waste (4), stabilize the sponge skeleton (36), and defend against pathogens, predators, or competitors via the production of bioactive secondary metabolites (5, 21, 39, 47). In addition, some cyanobacterial symbionts are a source of nutrients for their hosts because of their photosynthetic and nitrogen-fixing abilities (3, 57).
There is a large body of knowledge concerning the possible roles of bacterial associates in sponges, but whether sponges harbor specific bacterial communities deserves more detailed study. Many previous investigations demonstrated that associated bacterial communities in certain species of sponges were highly similar and consistently different from the bacterial communities in the ambient environment (6, 22, 46). For instance, Hentschel et al. (16) showed that there were uniform microbial communities in the marine sponges Aplysina aerophoba, Rhopaloeides odorabile, and Theonella swinhoei from different geographic regions that were distinct from those in the water column or in sediments. In addition, Friedrich et al. (12) demonstrated that the composition of sponge-associated bacterial communities was resistant to environmental perturbations resulting from transplantation to different habitats. These studies suggest that there is a stable, specific, and perhaps mutualistic relationship between the two types of organisms (44). However, some bacterium-sponge symbioses do not appear to be consistent. For instance, Wichels et al. (56) demonstrated that the bacterial communities associated with the North Sea sponge Halichondria panicea varied substantially over time. Qian et al. (34) and Lee et al. (25) showed that the congeneric Callyspongia and Mycale sponges from different biogeographic regions had different bacterial associates. Although the discrepancies may be attributed to the methods employed in different studies, whether consistent sponge-bacterium associations occur in different species or genera of sponges remains unclear. Furthermore, most of the previous studies focused on only one sponge species or involved a few sponges from geographically separated regions. So far, there has been no large-scale study comparing the bacterial communities associated with sponges of different genera and species.
In this study, we compared the bacterial communities associated with marine sponges around San Juan Island, Washington, in order to investigate the specificity of congeneric and conspecific sponge-associated bacterial communities. Most previous studies have relied on only one method to assess the bacterial communities, limiting the resolution of community assessments. In addition, different studies have used different approaches to address the same question, creating uncertainties and making generalizations difficult. Here, we employed both culture-independent and -dependent approaches to compensate for the limitations of different methods and to obtain a more reliable assessment of the associated bacterial communities. We used two DNA fingerprinting techniques, terminal restriction fragment length polymorphism (TRFLP) (27) and denaturing gradient gel electrophoresis (DGGE) (11). TRFLP is an effective, sensitive, high-throughput technique that differentiates bacterial community structures, while DGGE allows subsequent identification of bacteria of interest by excision and sequencing of specific bands. Both techniques have been successfully and widely used for characterization of bacterial communities in marine samples (13, 26, 28, 43, 56). In addition, bacteria associated with sponges were isolated using cultivation methods and identified by comparative analysis of 16S rRNA gene sequences. Phylogenetic affiliations of the isolates were determined and compared for different samples in order to study the specificity of the sponge-associated culturable bacterial communities.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Sponge samples collected near San Juan Island, WA, in April 2005a
|
|
View this table: [in a new window] |
TABLE 2. Universal primers used for PCR amplification in this study
|
DGGE analysis of bacterial communities.
PCR products were mixed with loading buffer and loaded onto a 6% acrylamide gel with a 35 to 70% denaturing gradient using urea-formamide as the denaturing agent. Electrophoresis was performed using the D-Code system (Bio-Rad, United States) in 1x TAE (20 mM Tris base, 10 mM sodium acetate, 0.5 mM EDTA) at a constant temperature of 60°C and a voltage of 125 V for 18 h. The gel was stained with 1x SYBR Gold (Molecular Probes, United States) for 15 min and photographed with the Alpha Imager 2200 gel documentation system (Alpha Innotech).
Excision and sequencing of DGGE bands.
Bands of interest were excised from the gel with sterile razor blades, washed twice with 100 µl of autoclaved ddH2O, and incubated in 50 µl of autoclaved ddH2O for 24 h at 4°C. DNA dissolved from the gel sections was used as the template for PCR amplification with the 341F and 907R primers. The PCR mixture and conditions were the same as those mentioned above. PCR products were purified using a PCR purification mini kit (Watson Biotechnologies, China) according to the manufacturer's manual. The 16S rRNA gene amplicon of each band was sequenced from both ends using primers 341F and 907R. Each cycle sequencing mixture contained 5.5 µl of purified DNA, 4 µl of DYEnamic ET dye terminator (Amersham Biosciences), and 0.05 µM of one of the primers. Cycle sequencing was performed with the thermal cycler using the following conditions: 24 cycles of denaturing at 96°C for 30 s, annealing at 50°C for 15 s, and extension at 60°C for 4 min. Cycle sequencing products were purified by ethanol precipitation and subjected to electrophoresis with a MegaBACE genetic analyzer (Amersham Biosciences) operated in the sequencing mode. Nucleotide sequences obtained with the two primers were assembled using the computer software Sequencher (Gene Codes Corporation) and then compared with the DNA sequences deposited in the NCBI (National Center for Biotechnology Information) web server using BLAST (Basic Local Alignment Search Tool).
Bacterial isolation and identification.
Samples used for bacterial isolation were diluted 1:10 and 1:100 with AFSW. Two hundred microliters of each dilution was spread on plates containing 1x nutrient agar (3 g of yeast, 5 g of peptone, and 12 g of bacteriological agar in 1 liter of AFSW), 0.1x nutrient agar (0.3 g of yeast, 0.5 g of peptone, and 12 g of bacteriological agar in 1 liter of AFSW), and tryptone soya agar (30 g of tryptone soya broth and 12 g of bacteriological agar in 1 liter of AFSW) and incubated at 22°C with a photoperiod consisting of 15 h of light and 9 h of darkness. Bacterial colonies were examined with a dissecting microscope (at a magnification of x10) to determine morphological characteristics, including color, shape, size, and surface texture. Colonies with morphological differences were isolated and purified. The same morphological differentiation standard was applied to all species and replicates. DNAs were extracted from pure bacterial colonies using the method of Valsecchi (52). The 16S rRNA genes in the crude extracts were amplified by PCR using primers 26F and 1055R (Table 2). PCR was carried out using a mixture (total volume, 25 µl) containing 2 µl of DNA extract, 1 U of Taq polymerase (Amersham Biosciences), 0.25 mM of deoxynucleoside triphosphates, 0.1 µM of each primer, and 1x PCR buffer under the following conditions: 95°C for 8 min, followed by 37 cycles of 95°C for 1 min, annealing at 56°C for 2 min, and extension at 72°C for 2 min and then a final extension at 72°C for 5 min. PCR products were purified and sequenced as described above for the sequencing of DGGE bands except that the primers used (primers 26F and 1055R) were different. The closest phylogenetic affiliation for each isolate was obtained by comparing the assembled nucleotide sequences with those deposited in the GenBank database using BLAST. The phylogenetic affiliation was verified by using RDP classifier (53).
Statistical analysis.
The average numbers of TRFs generated from TRFLP analysis of different samples (n = 3) were calculated and compared by using one-way analysis of variance, followed by a Dunnett test. The similarity of the TRF profiles for different samples was determined using cluster analysis based on the presence or absence of TRFs that were different sizes in each replicate. Although the relative signal intensities of TRFs may reflect the relative abundance of different PCR amplicons (i.e., different bacterial types), they were not used in the statistical analysis because of the likelihood of preferential PCR amplification of numerically dominant bacterial types in the samples (hence, the signal intensity may not be a true reflection of the relative bacterial abundance) (42). TRFs that accounted for <0.01% of the total fluorescence intensity, TRFs that were <35 bp long, and TRFs that were >500 bp long were excluded from the statistical analysis to screen out background noise, to avoid pseudo-TRFs from primers, and to avoid inaccurate size determinations, respectively. The Pearson coefficient was calculated based on the total number of TRFs observed in all samples and the presence of these TRFs in individual replicates to produce a similarity matrix. Agglomerative hierarchical clustering based on the matrix was performed using the PRIMER v3.1 program (Plymouth Marine Laboratory, United Kingdom) to generate a dendrogram showing the similarity among samples (7). For DGGE analysis, the DGGE band position and intensity were determined by the GELCOMPAR II software (Applied Maths), with manual modifications. Band matching was performed with 1% position tolerance and 1% optimization. A similarity matrix was calculated based on the Pearson coefficient and used to construct a dendrogram (7). Since DNAs extracted from replicate samples were pooled for DGGE analysis, no replicates were shown in the dendrogram. For isolation of bacteria, the occurrence of different bacteria differentiated by comparative analysis of the 16S rRNA gene sequences in replicates was recorded.
Nucleotide sequence accession numbers.
The 16S rRNA gene sequences obtained in this study have been deposited in the GenBank database under accession numbers FJ596288 to FJ596553.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 1. Numbers of TRFs (three replicates) and DGGE bands generated using 16S rRNA amplicons from the bacterial communities associated with different sponges and in the water column. Samples for the water column analysis were collected at site 1 (Roche Harbor). Data that are significantly different at = 0.05 in a one-way analysis of variance followed by a Dunnett test are indicated by different letters above the filled bars. Table 1 explains sample abbreviations.
|
![]() View larger version (17K): [in a new window] |
FIG. 2. Dendrogram derived from the TRFLP analysis, showing the relatedness of bacterial communities (three replicates each) associated with different sponges and in the water column. The dendrogram was constructed by using the Pearson correlation of similarity. Table 1 explains sample abbreviations.
|
![]() View larger version (15K): [in a new window] |
FIG. 3. Dendrogram (left panel) derived from DGGE profiles (right panel) of PCR-amplified 16S rRNA fragments from the bacterial communities associated with different sponges and in the water column at site 1. A solid line in the DGGE profiles indicates a band with a strong signal detected by the software, while a dashed line indicates a band with a weak signal. Arrows a to i indicate prominent and unique bands that were found in the sponge samples but not in the seawater sample. These bands were excised and sequenced for identification of bacteria (Fig. 4). The dendrogram was constructed by using the Pearson correlation of similarity. Table 1 explains sample abbreviations.
|
![]() View larger version (33K): [in a new window] |
FIG. 4. Neighbor-joining tree based on 16S rRNA gene sequence comparisons, showing the estimated phylogenetic relationships between sequences retrieved from the DGGE bands (Fig. 3) and their closest relatives. Asterisks indicate nodes that are also present in a maximum-parsimony tree; the bold lines indicate branches that are also present in a tree constructed by the unweighted-pair group method using average linkages. Bootstrap values greater than 50%, based on 1,000 replicates, are indicated at the nodes. GenBank accession numbers are in parentheses. Sequences that fell into the sponge-specific clusters defined by Taylor et al. (44) are indicated by pound signs. Sequences retrieved from sponge-associated bacteria are indicated by bold type, while sequences retrieved from bacteria in seawater or sediments are indicated by light type. Bar, 2 nucleotide substitutions per 100 nucleotides.
|
1,000 bp) using BLAST and the RDP classifier placed these isolates into 45 different genera in four phyla, the Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes (Fig. 5). Of the four phyla, Proteobacteria and Firmicutes were the most dominant groups, accounting for 71% (in MyI-2) to 100% (in XH-4) of the culturable bacterial communities, but the relative percentages varied substantially among samples (Fig. 5). For instance, >80% of bacteria isolated from the seawater samples belonged to Proteobacteria, but only 17% of them belonged to Firmicutes. In contrast, Firmicutes, especially the genus Bacillus, was the most dominant group in the sponge samples and could account for up to 82% of the culturable bacterial communities (in XH-4), and the percentage of Proteobacteria could be as low as 2% (in MA-1). In general, the number and type of bacteria isolated from conspecific or congeneric sponges varied substantially.
![]() View larger version (51K): [in a new window] |
FIG. 5. Percentages of different genera recovered from the sponge and water samples. The phylogenetic affiliations of bacteria were determined by comparing the 16S rRNA gene sequences of the isolates with those deposited in the GenBank database. The letters in parentheses and on the right indicate phyla, as follows: P and P, Proteobacteria; F and F, Firmicutes; A and A, Actinobacteria; and B and B, Bacteroidetes. Dashed lines separate different phyla. Table 1 explains sample abbreviations.
|
|
|
|---|
The number of bacterial types (or ribotypes) in the sponge samples, as evaluated by the number of TRFs and DGGE bands, varied among species (Fig. 2); H. rufescens and M. incrustans had the highest and lowest diversity of ribotypes, respectively, and H. panicea had various numbers of TRFs and DGGE bands depending on the site. At site 1, the sponge-associated bacterial communities were different in different species and were distinct from those in the water column (Fig. 3 and 4). Similarly, distinct bacterial communities in sponges and the ambient water have been found in previous investigations (18, 19, 37, 55). It is generally believed that sponge-associated microbes benefit from an enriched nutrient environment due to the release of secondary metabolites from their sponge hosts (8, 18). However, this benefit may result in an increase or decrease in bacterial diversity depending on how the microbes utilize the metabolites and how they interact with each other.
Close associations of a variety of microorganisms, especially bacteria, are found in many sponges. So far, 14 bacterial phyla from sponges have been recognized, and the most frequently recovered sequences belonged to the phyla Acidobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, Cyanobacteria, Firmicutes, and Proteobacteria (44). In our study, sequences belonging to the phyla Proteobacteria and Firmicutes and isolates belonging to the phyla Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria were retrieved from the sponge samples (Fig. 5), indicating that these phyla were present in our sponge-associated bacterial communities. Five of the nine bands selected from our sponge samples were affiliated with bacteria whose sequences were retrieved from other sponges. Phylogenetic analysis showed that three of these five bands (bands g, h, and i) fell into the sponge-specific sequence clusters defined by Taylor et al. (44). Our results further support the argument that some sponges from geographically separated regions share some specific bacteria and sponge-specific clusters do exist (16, 44). In addition, the low level of similarity of some of our sequences with sequences in the GenBank database indicated that there is an unexplored bacterial group in sponges and the potential of identifying novel bacteria. In fact, we have already identified two novel bacteria isolated from M. incrustans (23) and H. panicea (24).
When the bacterial communities in conspecific sponges were compared, highly similar bacterial communities were observed for certain sponge species. We demonstrated for the first time that both H. rufescens and M. incrustans from different locations harbored highly similar bacterial communities, as shown by both TRFLP and DGGE analyses (Fig. 3 and 4). This is the first study of the bacterial communities in these two species, while only one previous study examined the fungal community in M. incrustans (31). In contrast to the great consistency in the bacterial communities associated with these two species at different sites, the bacterial communities associated with individuals of H. panicea from different sites varied substantially (Fig. 3 and 4). Bacterial communities in H. panicea have been investigated in a number of studies (1, 20, 29, 41, 56). A group of bacteria with 16S rRNA sequences closely related to sequences of the bacterial group Roseobacter was consistently recovered from geographically separated H. panicea sponges, yet different bacterial copopulations also occurred seasonally or on a small geographical scale with occasional dominance in this species. The two sequences retrieved from H. panicea in our study did not resemble any of the sequences previously recovered from this species. These findings indicated that there are considerable variations in the bacterial communities associated with H. panicea from different locations, and the previous results are in good agreement with our observations for the H. panicea samples.
It is not known why consistent bacterial communities were found in the sponges M. incrustans and H. rufescens but not in the sponge H. panicea. The type of bacteria present in a sponge depends on both environmental and sponge-related factors, such as structural organization and life history strategy (47). M. incrustans and H. rufescens are encrusting sponges commonly found on the surface of scallops that live in a defined geographical range and have a stable environment, while H. panicea occurs in a wide range of forms and locations and can be difficult to identify morphologically (R. von Soest, personal communication; our observations). Evolutionary history may lead to different degrees of specificity in the types of bacterial associates that a sponge can host. It is generally believed that there may have been three different evolutionary scenarios for the sponge-bacterium association, as summarized by Taylor et al. (44): vertical transmission, horizontal transmission, and environmental acquisition. At site 1, the sponge-associated bacterial communities varied among sponge species and were all different from the indigenous planktonic bacterial community; thus, neither purely environmental acquisition nor horizontal transfer between sponges can explain this observation perfectly. Vertical transmission of bacterial symbionts may explain our observations for M. incrustans and H. rufescens but not for other sponges. However, since the reproductive biology of these sponges is poorly understood and the host sponge phylogeny is still not clearly resolved, information on whether the sponge hosts and bacterial associates have coevolved is necessary to obtain further evidence of the symbiont transmission modes of these sponges.
As determined by both TRFLP and DGGE analyses, the bacterial communities in congeneric sponges (M. adhaerens and M. loveni; H. rufescens and H. foraminosa) were in different clusters (Fig. 3 and 4), indicating that there are substantial variations in the bacterial communities in the two genera. This is in good agreement with other studies which demonstrated that there are distinct and diverse bacterial communities in congeneric Mycale or Haliclona spp. (9, 25), providing further evidence that bacterial communities do not remain consistent in these genera.
In this study, both culture-dependent and culture-independent methods were employed to study the sponge-associated bacterial communities. TRFLP and DGGE allow easy and quick comparison of community profiles for microbial assemblages. Although these two methods resulted in slightly different clustering patterns (Fig. 3 and 4), the different patterns did not affect our broad overall conclusion that members of certain sponge species from different locations shared highly similar and consistent bacterial communities, suggesting that some sponges might exhibit specificity in their associated bacterial communities. There were even larger discrepancies when the results obtained with culture-dependent and -independent methods were compared (Fig. 3 to 5). None of the sequences retrieved from excised DGGE bands resembled the sequences from our bacterial isolates. In addition, the sequences of only two isolates from our study most closely matched other sponge-derived sequences from previous studies. The limitations of the culture-dependent methodology (2, 20) may explain the observed differences for the methods used. However, using the same medium selective pressure, we showed that the culturable bacterial communities in the sponges were dominated by Firmicutes, especially Bacillus, and that they were very different from the culturable bacterial communities in the water column, in which Proteobacteria, especially Pseudoalteromonas, was the major group (Fig. 5). This supported our observation resulting from TRFLP and DGGE analyses that the bacterial community in the water column was very distinct from the bacterial communities associated with the sponges. One advantage of the culture method is that it allows isolation of pure strains for further experiments. It should be noted that, although clone libraries may help to identify sponge-specific bacteria, considering the scope of this study, the enormous work involved, and the complex communities analyzed, such libraries were not used in this study.
This study provided information on variations in bacterial communities within and between a variety of sponge species. Clearly, generalizations concerning the specificity of sponge-associated bacterial communities are risky since the degrees of specificity may be different for different sponges. Our results indicated that the bacterial communities in certain species of sponges, including M. incrustans and H. rufescens, from different locations were highly similar, suggesting that these sponges consistently host specific bacterial communities. The mechanism that the sponges employ to maintain this consistent association and over what scales the associations are maintained are still not known. Further studies of whether these associations are consistent for larger geographic scales, as well as studies of the mode of symbiont transmission, are currently under way.
This study was supported by a grant from the China Ocean Mineral Resources Research and Development Association (COMRRDA06/07.SC02) and the KAUST Global Academic Partnership Program (KAUST005-CML.07/08) to P.-Y. Qian.
Published ahead of print on 10 April 2009. ![]()
|
|
|---|
-proteobacterium, "Candidatus Entheonella palauensis." Mar. Biol. 136:969-977.[CrossRef]
-proteobacterium. Mar. Biol. 138:843-851.[CrossRef]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»