Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2002, p. 6292-6299, Vol. 68, No. 12
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.12.6292-6299.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
and Wendy Morrill
Department of Biochemistry, Microbiology, and Molecular Biology, University of Maine, Orono, Maine
Received 29 May 2002/ Accepted 24 September 2002
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The mechanism by which shipworms digest wood is thought to involve gram-negative bacterial endosymbionts found within the gills. These bacteria were first observed by transmission electron microscopy in the shipworm Bankia australis (32) and were proposed to synthesize essential amino acids lacking in the shipworm's wood-based diet (31). However, the axenic cultivation of a cellulolytic, dinitrogen-fixing bacterium from the gills of several shipworm species suggests that the symbionts play a more direct role in wood feeding (41). This symbiotic bacterium, recently named Teredinibacter turnerae (10), was shown to reside within specialized epithelial cells (bacteriocytes) in the shipworm L. pedicellatus (9) in a gill region previously referred to as the Gland of Deshayes (36). Enzymes produced by this symbiont were proposed to facilitate the digestion of wood and to provide an internal source of combined nitrogen to supplement the host's protein-deficient diet (16, 41). Phylogenetic analysis, based on 16S rRNA sequences, showed T. turnerae to be a member of the gamma subdivision of the proteobacteria (9, 10).
As in most intracellular bacterial endosymbioses, the symbiont population L. pedicellatus has been presumed to be a monoculture (9, 37, 41). However, in a few cases, multiple intracellular endosymbiont species have been observed within a single host. For example, the gill bacteriocytes of certain mussel species from hydrothermal vents and cold seeps (8, 12) contain two coexisting symbiont species. Multiple endosymbioses have also been reported in other invertebrate taxa; however, in these cases, the symbionts either were extracellular (29, 30), were localized to separate tissues (28), or were pathogens that infect only a portion of the host population (34). The purpose of this investigation was to explore the diversity and phylogenetic composition of symbiont populations in the gills of L. pedicellatus. The results provide evidence for an intracellular bacterial consortium potentially more complex than others previously described.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Isolation of T. turnerae.
T. turnerae was isolated from the gill tissue of lab-reared L. pedicellatus as described previously (41) by serial dilution (101 to 1010) of gill homogenate in shipworm basal medium (SBM) containing 0.2% agar with cellulose (Sigmacell, 2 g/liter; Sigma Chemical) added as a carbon source and with no added source of combined nitrogen. Axenic stock cultures were maintained in the same medium. SBM liquid cultures, used for nucleic acid extractions, contained SBM supplemented with Sigmacell (2 g/liter) and ammonium chloride (5 mM).
Nucleic acid preparation, PCR amplification, and cloning.
DNA was extracted from
100 mg of fresh gill tissue from each adult specimen of L. pedicellatus as described previously (7). DNA extraction from T. turnerae cultures was done as described previously (10). Bacterial 16S rRNA genes were amplified by PCR (35) from bulk nucleic acids extracted from L. pedicellatus gill tissue (two adult specimens) and cultivated T. turnerae cells. Nearly complete 16S rRNA gene sequences were amplified from each DNA preparation independently by using bacterial-domain-specific primers 27f and 1492r (23). PCRs were performed in 50-µl volumes as described previously (7) with 20 or 35 cycles and either 1 ng of gill DNA/µl or 0.2 ng of bacterial genomic DNA/µl.
PCR amplification products from L. pedicellatus gill tissue were cloned into Escherichia coli by using the pGEM-T vector system (Promega) following the manufacturer's recommended protocol. Plasmid DNAs were recovered by using the Wizard Plus SV Miniprep DNA purification system (Promega) as directed by the manufacturer and were screened for the presence of insert DNA by agarose gel electrophoresis. A combined total of 43 clones containing inserts of the appropriate size (
1.5 kb as determined by agarose gel electrophoresis) were selected and screened by sequencing both strands of a region corresponding to E. coli positions 357 to 907 (
500 bp) (23).
Sequencing.
Sequencing reactions were performed with an ABI PRISM Dye Terminator cycle sequencing ready reaction kit by using AmpliTaq DNA polymerase FS on an ABI 373A DNA sequencer. Approximately 500 bp of the double-strand sequence was determined for each clone by using primers 357f and 907r (23). Nearly complete double-strand sequences of 16S rRNA genes (1,390 to 1,441 bp) were determined for selected clones and PCR products by using primers 27f, 357f, 704f, 690r, 907r, 1101r, and 1492r (23).
Fluorescent in situ hybridization. (i) Tissue and cell preparation.
Gill tissues were rinsed in sterile seawater and fixed in Davidson's fixative (10% glycerol, 8% formaldehyde, 28.5% ethanol, 10% acetic acid, 30% seawater) (20) overnight at room temperature. Tissues were then dehydrated through a series of graded ethanol wash steps (50% once, 70% once, 80% once, 95% twice, and 100% thrice for 30 min each), followed by two washes in xylene at room temperature (45 min each). Tissues were infiltrated with paraffin (three changes, 30 min, 60°C). Consecutive tissue sections (5-µm thickness) were cut and mounted on APTS (3-aminopropyltriethoxysilane; Sigma)-treated slides. Slides were then baked at 57°C for 1 h. Paraffin was removed by three washes in xylene (15 min) and rehydrated through a graded ethanol series (inverse of that described above). Slides were washed in 0.2 M HCl and rinsed in 0.02 M Tris-HCl, pH 8.2. Following a 5-min treatment with proteinase K (0.5 µl/ml), slides were rinsed as described above, postfixed in 3.7% formaldehyde in 0.02 M Tris-HCl, and rinsed again. After drying, slides were treated with triethanolamine-acetic anhydride (24). T. turnerae cells grown in vitro and symbiont cells from gill tissue homogenates were fixed for 1 h to overnight in 3.7% formaldehyde in 90% methanol. Cells were dotted on APTS-coated slides and air dried. Slides were postfixed in 90% methanol-3.7% formaldehyde for 20 min at room temperature (5) and then rinsed in distilled water (dH2O) for 1 min. Slides were then treated in triethanolamine-acetic anhydride, and hybridizations were performed as previously described.
(ii) Hybridization conditions.
In situ hybridizations were performed on gill tissue sections and gill tissue homogenates from four adult shipworms sacrificed on 31 October 1996, 6 November 1996, 17 August 1997, and 14 October 1997. Hybridization buffer (750 mM NaCl, 5 mM EDTA, 100 mM Tris-HCl [pH 7.8], 0.025% sodium dodecyl sulfate) was added to each section along with 50 ng of each probe. Sections were covered with coverslips cut from Parafilm (American National Can) to prevent evaporation and were incubated in a humid chamber at 50°C for 3 to 24 h. Slides were washed three times (5 min) in 0.2x SET (30 mM NaCl, 200 nM EDTA, 4 mM Tris-HCl [pH 7.8]) at 50°C (
11 to 14°C below melting temperature). Slides were then air dried and mounted with Vectashield mounting medium (Vector Laboratories).
(iii) Oligonucleotide design, synthesis, and labeling.
Oligonucleotide (DNA) test probes were designed as described previously (9) to hybridize specifically with a variable region (E. coli positions 637 to 662) of the RNA-like strand of the consensus sequence of clone groups LP1 through LP6. These clone-specific probes were labeled with either a red or green fluorochrome (Texas Red or fluorescein isothiocyanate [FITC], respectively) and were designated lpx-r or lpx-g, respectively, where x is a number designating the clone consensus group (1 through 6) for which the probe was designed, r is red fluorochrome, and g is green fluorochrome. Negative-control probes were labeled similarly and designated lpxn-r or lpxn-g. Each negative-control probe contained two differences (G
C) from its target sequence at E. coli positions 648 and 656. A bacterial-domain-specific probe corresponding to E. coli positions 338 to 355 (2) was labeled with FITC and designated bac-g. Probe sequences are listed in Table 1.
|
Optimal hybridization and wash conditions for each probe were determined empirically by in situ hybridization to fixed and mounted cells under successively increasing stringencies until fluorescence due to probe hybridization was undetectable. Hybridization specificity was also confirmed for each probe by competitive hybridization with unlabeled exact-match probes. Specimens were observed by using a Nikon Labophot 2 epifluoresence microscope with Texas Red, fluorescein, and Texas Red, fluorescein, dual-wavelength filter sets (Chroma). Photomicrographs were taken with a Nikon FX-35DX camera with Kodak Gold (200 ASA) color print film or Kodak Elite II (200 ASA) slide film. Exposure times ranged from 30 to 524 s. Exposure times for negative-control probes were doubled compared to the corresponding exposure times for experimental probes to rule out the possibility of faint, nonspecific staining.
Phylogenetic analysis.
Sequences of 16S rRNA gene clone groups LP1 (1,390 bp, three pooled clones), LP2 (1,397 bp, two pooled clones), LP3 (1,441 bp, one clone), and LP4(Tt) (1,352 bp, one clone) and the direct-PCR sequence of the T. turnerae 16S rRNA gene were determined bidirectionally and aligned with 16S rRNA genes from selected gamma proteobacteria and environmental samples contained in the GenBank and Ribosomal Database Project databases. Database and culture collection accession numbers are reported in Fig. 1. Reference taxa were selected from among members of the gamma subdivision of the proteobacteria most closely related to T. turnerae, including taxa with known cellulolytic activity and the gill symbionts of other bivalves, as described previously (10).
|
Phylogenetic analyses were performed using algorithms contained in PAUP* (Phylogenetic Analysis Using Parsimony, version 4.0b10) (D. L. Swofford, PAUP* 4.0. Sinauer, Sunderland, Mass., 1997). Maximum likelihood trees were inferred by using the HKY85 substitution model (18) with empirical base frequencies, equal rates for all sites, and a transition/transversion ratio of 2. Maximum parsimony analyses were performed with character states optimized by accelerated transformation (ACCTRAN). In all analyses, branch swapping was by tree bisection-reconnection, with characters weighted equally (weight = 1) and gaps treated as missing data. Bootstrap analyses were performed by using the full heuristic search option with 1,000 replicates. Bootstrap analyses under maximum likelihood were limited to a reduced taxon set (20 taxa) to accommodate the greater computational demands of this technique.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
99% sequence identity to the consensus sequence for that group over the
500 bp of sequence determined. Of the 43 clones, 31 were included in group LP1 (72%), 4 in group LP5 (9%), 2 in group LP2 (5%), and 1 each in groups LP3, LP4, and LP6 (2%). Consensus sequence LP4 was identical to the sequence determined for pure-cultured cells of T. turnerae and is hereafter referred to as LP4(Tt). Inserts contained in three remaining clones (3W16, 3W49, and 3W59) were determined to be probable chimeras and so were not assigned to consensus groups. All three putative chimeras were recovered from PCR 1 (35 cycles), consistent with the observation that chimera formation may increase with an increasing number of cycles during PCR amplification (26, 38, 39).
Sequence variability within consensus groups.
Sequence variability was observed within consensus groups; however, the frequency and pattern of most nucleotide substitutions were consistent with that expected due to polymerase error. For example, the 31 clones corresponding to the largest consensus group (LP1) contained a total of 22 divergences from the group consensus sequence, distributed over 18 positions (an average of <1 difference per 700 bp). At 15 (83%) of these variable positions, a difference from the consensus was observed in only 1 of the 31 clones. Only three positions (0.6% of all positions) contained differences from the consensus in two or more clones. In only two of these positions did differences covary (two clones shared the same difference [A
G] at E. coli position 498, and four clones shared the same difference [T
G] at E. coli position 610). The total number of differences observed among the 31 clones could be explained by a nucleotide misincorporation rate of about 5 x 10-5 incorporation errors/bp/cycle. This is within the range of values (2 x 10-4 to 7.2 x 10-5 errors/bp/cycle) reported for Taq polymerase under similar PCR conditions (11, 21, 27). Furthermore, none of the differences that occurred in stems were matched by compensating changes in the complementary strand, and five of these differences resulted in mismatches that could destabilize conserved helices. These results indicate that sequence polymorphism in the symbiont population represented by this consensus group is low by comparison to between-group variation (
1.5 to 8.1%). Therefore, in this investigation, each consensus group sequence was treated as a unique rRNA type, or ribotype. Each ribotype was considered to represent a proxy for a putative symbiont phylogenetic type, or phylotype.
Phylogenetic analyses of ribotypes.
The phylogenetic analyses presented here include only those ribotypes (LP1 through LP4) that were detected in tissue and tissue homogenates by in situ hybridization (see below). Pairwise sequence identity (calculated as 1-Kimura 2-parameter distance) between these ribotypes ranged from 91.9 to 93.4%, with the exception of LP1 versus LP3, which shared a 98.5% sequence identity.
Comparison to published sequences (Ribosomal Database Project, GenBank) demonstrated that the L. pedicellatus putative symbiont ribotypes are derived from gamma proteobacteria. The sequence most similar to ribotypes LP1 through LP4 (93.3 to 95.3%) is that reported for an uncultivated symbiont from the shipworm Bankia setacea (37). Ribotypes LP5 and LP6 were distantly related to those detected in tissue, the former showing the greatest similarity (
90%) to epsilon proteobacteria, while the latter was weakly associated with alpha proteobacteria.
A heuristic search by maximum likelihood analysis retained a single best tree with an ln likelihood of 11,037.2 (Fig. 1). Maximum parsimony retained the four best trees of 1,732 steps (consistency index = 0.4330, homoplasy index = 0.5670, retention index = 0.5290, and rescaled consistency index = 0.2291). Trees inferred by both methods were topologically identical with respect to all nodes receiving significant (>70%) (19) bootstrap support.
In these analyses, the T. turnerae sequence (identical to ribotype LP4) and ribotypes LP1, LP2, and LP3 were nested within a single clade that also included the 16S rRNA genes of an uncultivated symbiont from the shipworm B. setacea (37); an uncultivated bacterium, environmental sample BD2-13 (GenBank accession no. AB015541), from deep-sea sediments (1,521 m) in Sugura Bay, Japan (1,521 m; 34°55'N, 138°39'E) (25); and an undescribed marine heterotrophic bacterium cultivated from seawater collected along the coast of La Jolla, Calif. (SCB11, GenBank accession no. Z31658) (33). This clade receives significant bootstrap support (19) in both maximum parsimony (89%) and maximum likelihood (77%) analyses and is hereafter referred to informally as the Teredinibacter clade.
No strong support was observed for any specific relationship between the Teredinibacter clade and other named genera. Analyses described here suggest a peripheral relationship with either Pseudomonas spp. (sensu stricto) or Oceanospirillum spp., although with weak to insignificant bootstrap support. Consistent with previous reports, a possible relationship is suggested between the Teredinibacter clade and strain 2-40, an unnamed bacterium isolated from a marine sea grass bed in the Chesapeake Bay (15).
These results indicate that ribotypes LP1 through LP3 are derived from bacterial species that are candidates for classification within the genus Teredinibacter, the only formally described taxon in this lineage (10). Though no convention has been established to correlate degree of 16S rRNA sequence identity with taxonomic standing, most taxa that are considered to be separate species within the same genus differ by at least 1.5% (13). By this criterion, ribotypes LP1, LP2, and LP3, which differ from that of T. turnerae by >6.5%, may represent one or more new species within this genus. However, taxonomic classification must await more complete characterization of the putative symbionts and related free-living bacteria.
Fluorescent in situ hybridization.(i) Single-probe hybridizations.
When applied to tissue, probes lp1-r through lp4(Tt)-r [specific for ribotypes LP1 through LP4(Tt)] resulted in fluorescent staining of the symbiont-containing vesicles in sectioned bacteriocytes and individual bacterial cells in gill homogenates (single-probe data not shown). Staining and appearance of symbionts within bacteriocytes was similar to that observed in a previous study (9). Fluorescent staining by the bacterial-domain-specific positive-control probe bac-g was also observed in the same regions. No staining of gill tissue was detected by using ribotype-specific probes lp5-r or lp6-r or negative-control probes lp1n-r through lp4(Tt)n-r. These results indicate that ribotypes LP1 through LP4(Tt) derive from symbionts that reside within gill bacteriocytes. Only these symbiont ribotypes were examined further in this investigation.
(ii) Dual-probe hybridizations in individual symbiont cells.
Dual-probe hybridizations, performed on individual symbiont cells released from bacteriocytes by gentle homogenization, demonstrated that no single ribotype-specific probe stained all symbiont cells visible in a given homogenate. When any red-labeled ribotype-specific probe [lp1-r through lp4(Tt)-r] was applied in combination with the green-labeled bacterial-domain-specific probe, both green-stained cells and yellow- to orange-stained cells were observed (Fig. 2d and e). These results indicate that, although each ribotype-specific target sequence could be detected in some symbiont cells, no single ribotype could account for all bacteria present in the gill homogenates. This is consistent with the possibility that the four ribotypes correspond to distinct symbiont phylotypes that coexist in the gill tissue.
|
Although these results cannot rule out the possibility that some symbiont cells may express multiple ribotypes, hybridization results indicate that the cells of T. turnerae express only one of the four ribotypes when grown in pure culture. Only one of the four ribotype-specific probes, probe lp4(Tt)-r, was observed to stain pure-cultured log-phase cells of T. turnerae. Furthermore, in dual-probe experiments with probes lp4-r and bac-g, all cells were stained by both probes (Fig. 2a-c), indicating a homogeneous bacterial population expressing ribotype LP4.
(iii) Dual-probe hybridizations in bacteriocytes.
Dual-probe hybridizations performed in sectioned gill tissue produced distinct and specific staining patterns for each ribotype. This is evident in Fig. 3, panels 1 to 6, which shows a series of tissue sections (panels 1 to 5 are consecutive) from the gill of a single specimen of L. pedicellatus. For example, in section 1, treated with probes lp4(Tt)-r and bac-g, approximately one-third of the bacteriocytes were stained orange to yellow, indicating hybridization of both probes while the remaining bacteriocytes were stained green, indicating hybridization of only the bacterial-domain-specific probe bac-g. However, a different and nearly inverse staining pattern is shown for probe lp1-r in section 2. That is, cells that stained green in panel 1 appear orange to yellow in panel 2 and vice versa. This can be seen clearly in the boxed regions in Fig. 3, which indicate bacteriocytes that can be identified by shape, size, and position in several sections. Bacteriocyte a was stained by probe lp1-r (panel 2) but not by probe lp4(Tt)-r (panel 1), while the inverse is shown for bacteriocyte c. These results indicate that the symbiont ribotypes LP1 and LP4(Tt) occur in separate and mutually exclusive subsets of bacteriocytes. This conclusion is further supported by the section shown in panel 5, in which ribotype-specific probes lp1-r and lp4(Tt)-g show patterns of staining consistent with those found in panels 1 and 2 but stained only with the primary colors red and green.
|
Distribution of ribotypes.
Staining by ribotype-specific probes also showed consistent patterns of distribution with respect to gill structure. Bacteriocytes labeled by the lp4(Tt) and lp2 probes were distributed symmetrically in the more dorsal and medial regions of the gill, closely adjoining the central axis dividing the right and left gill demibranchs. This axis runs approximately along the horizontal midline of panels 1 to 6 in Fig. 3. Probes lp1 and lp3 were widely distributed in the more lateral and ventral regions of the gill filaments.
Cooccurrence of ribotypes.
There are six possible pairwise staining combinations that can be generated with a group of four probes (1:2, 1:3, 1:4, 2:3, 2:4, and 3:4). To date, simultaneous staining by only two of these six possible pairwise combinations (1:3 and 2:4) has been observed in single bacteriocytes (Fig. 3, panels 8 and 9). This pattern of colocalization suggests that these ribotypes correspond to four distinct symbiont phylotypes that may be paired in specific combinations within bacteriocytes. This might indicate that symbiotic associations also occur between bacteria within bacteriocytes, i.e., bacteria-bacteria symbioses may exist within the larger animal-bacteria symbiosis.
Again, the observed pattern of cooccurring ribotypes could be also be explained by the differential expression of multiple divergent 16S rRNA gene copies present in the genome(s) of fewer than four symbiont phylotypes. However, pairwise evolutionary distances among the four ribotypes are larger than those typically observed among operons in a single genome. Multiple rRNA operons are common in bacteria, and sequence variability between 16S rRNA gene copies as high as 6.4% (1, 40, 42) has been reported. However, such large differences appear to be unusual (3). A survey of 14 complete bacterial-genome sequences (representing 55 rRNA operons) revealed a range of inter-operon variation in 16S rRNA genes from 0 to 1.23%, with over half of all taxa showing no inter-operon differences and nearly 80% (11 of 14) showing a less than 0.2% sequence difference among 16S rRNA gene copies (22). In comparison, the sequence divergences observed among consensus groups LP1 through LP4(Tt) are greater than 6.5% for all pairwise comparisons, except that of LP1 versus LP3, which differs by 1.5%. Thus, these differences are larger than those typically observed between rRNA operons in a single genome and, in fact, with the exception of the LP1 versus LP3 comparison, exceed any inter-operon differences yet reported.
PCR amplification and cloning results are also inconsistent with multiple divergent operons, at least for ribotypes LP1 and LP4. In amplifications performed by using bacterial-domain-specific primers with genomic DNA extracted from pure cultures of T. turnerae, a single sequence identical to the LP4 consensus sequence was detected. In contrast, the LP1 consensus sequence predominates in similar amplifications from L. pedicellatus gill DNA. It is difficult to explain how such differences could arise if both sample types contain a monoculture of the same bacterial phylotype.
Although the data presented here support the existence of multiple symbionts in L. pedicellatus, the functional significance of this multiple symbiosis remains unknown. In nature, wood degradation is a complex process that is rarely accomplished by a single microbial species (4). Thus, it is possible that each symbiont type in L. pedicellatus contributes different enzymes to a complex cellulolytic system. Alternatively, multiple symbiont phylotypes could provide flexibility to the host in adapting to a variety of wood types or environmental conditions. In either case, it may be revealing to determine whether multiple endosymbiosis is unique to this species or is a more widespread phenomenon. Unfortunately, only a small number of adult specimens have been examined in this study and these have been collected and reared in similar locations and under similar conditions. Therefore, although uniform results were observed in these specimens, it remains to be determined whether the composition of symbiont populations is consistent among all members of this species or whether variation may occur with the differing condition, environment, and developmental stage of the hosts. Finally, although T. turnerae (LP4) has been isolated from the gills of several shipworm species (10), it remains uncertain whether the other putative symbionts detected here are distributed widely among shipworm species. Further investigations are presently under way in this laboratory with a broader range of host specimens and taxa and more quantitative molecular and histological methods to better describe the diversity and complexity of symbiont populations in teredinid bivalves.
In summary, (i) four distinct but related symbiont ribotypes have been detected in the gills of L. pedicellatus, (ii) no single ribotype was detected in all symbiont cells in gill homogenates or in all bacteriocytes in gill tissue, (iii) the 16S rRNA gene sequence of one ribotype (LP4) is identical to that determined for T. turnerae grown in pure culture, (iv) ribotypes LP1 through LP3 have not been detected in cultivated cells of T. turnerae, (v) two ribotypes were detected in some but not all bacteriocytes, and (vi) within individual bacteriocytes, the ribotype combinations LP1 and LP3 and LP2 and LP4 have been observed but other possible pairwise combinations have not. Taken together, these observations provide strong evidence that at least two symbiont phylotypes (represented by ribotypes LP1 and LP3 and LP2 and LP4, respectively) and as many as four distinct symbiont phylotypes coexist in the gills of L. pedicellatus.
| ACKNOWLEDGMENTS |
|---|
We express our thanks to Scott Gallager (Woods Hole Oceanographic Institution) for providing the breeding stock of L. pedicellatus and Martin Polz (MIT) for thoughtful discussion of the manuscript.
| FOOTNOTES |
|---|
Present address: Woods Hole Oceanographic Institution, Woods Hole, MA 02543. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|