Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2006, p. 412-417, Vol. 72, No. 1
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.1.412-417.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Ocean Genome Legacy Foundation, Center for Marine Genomic Research and Conservation, 240 County Road, Ipswich, Massachusetts 01938,1 Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139,2 Department of Biochemistry, Microbiology, and Molecular Biology, University of Maine, 5735 Hitchner Hall, Orono, Maine 04469-57353
Received 13 July 2005/ Accepted 12 October 2005
|
|
|---|
subgroup of proteobacteria which contained at least five well-supported internal clades (phylotypes). By competitive quantitative PCR and constant denaturant capillary electrophoresis, we estimated the number and abundance of symbiont phylotypes in gill samples of 13 individual shipworm specimens. Phylotype composition varied greatly; however, in all specimens the numerically dominant symbiont belonged to one of two nearly mutually exclusive phylotypes, each of which was detected with similar frequencies among specimens. A third phylotype, containing the culturable symbiont Teredinibacter turnerae, was identified in nearly all specimens, and two additional phylotypes were observed more sporadically. Such extensive variation in ribotype and phylotype composition among host specimens adds to a growing body of evidence that microbial endosymbiont populations may be both complex and dynamic and suggests that such genetic variation should be evaluated with regard to physiological and ecological differentiation. |
|
|---|
The previous detection of multiple symbiont ribotypes in L. pedicellatus raises questions regarding the extent to which symbiont communities differ within this single host species. Specifically, does symbiont community structure vary among individuals or are these communities under tight selection or regulation? To address this question, we have quantitatively examined the diversity, abundance, and frequency of occurrence of ribotypes present in specimens drawn from a population of L. pedicellatus reared under controlled conditions in a long-term laboratory culture.
|
|
|---|
T. turnerae was isolated from gill tissue of laboratory-reared L. pedicellatus by the method of Waterbury et al. (31) as modified by Distel et al. (6), and DNA was extracted as described in reference 6.
PCR clone libraries.
Bacterial 16S rRNA genes were amplified by PCR (26) with Bacteria domain-specific primers 27f and 1492r (22) and the high-fidelity polymerase PfuTurbo (Stratagene, La Jolla, CA). Each PCR mixture (50-µl total volume) contained 50 ng of bulk nucleic acid extracted from gill tissue, 1x PFU buffer, 250 µM each deoxynucleoside triphosphate, 1.25 U of Pfu polymerase, and each primer at 200 nM. Cycling parameters were 2 min at 94°C, followed by 25 cycles of 1 min at 94°C, 1 min at 55°C, and 2 min at 72°C, followed by 8 min at 72°C. To eliminate potential heteroduplex artifacts, 5 µl of the completed PCR mixture was diluted 1:10 into a fresh PCR mixture and three cycles of reconditioning PCR were performed (28).
Amplified 16S rRNA gene fragments were cloned into plasmid vector pCR-Blunt II-TOPO (Invitrogen, Carlsbad, CA) and transformed into Escherichia coli (TOP10) cells according to the manufacturers instructions. Two 16S rRNA clone libraries were created and screened, one constructed from the whole gill of a single specimen (single-specimen library [SSL]) and the second from the pooled gills of four specimens (multiple-specimen library [MSL]) collected from a single lab-reared population of this host species. Plasmids were purified and prepared for sequencing with a RevPrep Orbit automated plasmid purification workstation (GeneMachines, San Carlos, CA).
Sequence analysis.
Cloned 16S rRNA gene fragments were sequenced bidirectionally with primers 27f, 1392r, 926f, 530f (22), and 519r2 (5'-GTGCCAGCMGCCGCGGTAAT-3') with the ABI PRISM Big Dye Terminator cycle sequencing kit, Amplitaq DNA polymerase, and an ABI 3730X1 DNA analyzer. Raw ABI sequence data were edited manually with Sequencher software v.3.1.2 (Gene Codes Corp., Ann Arbor, MI). For phylogenetic analyses, sequences were aligned manually with SeqLab (GCG v.10.1; Accelrys Inc.), taking into consideration secondary-structure information (3). Gaps were inserted to compensate for length variation in inferred loops and helices. Sites within regions of uncertain alignment were identified and eliminated from further analyses. In addition, nucleotide substitutions at positions conserved in greater than 98% of the known sequences and substitutions that result in mispairings that interrupt conserved helices (3) were treated as putative nucleotide misincorporation errors and were disregarded in all further analyses, leaving a final character set of 1,234 nucleotide positions. To identify putative chimeras, phylogenetic analyses were performed independently on 5' and 3' portions of the sequence alignments. Putative chimeras were considered confirmed if (i) different parent ribotypes were identified for the 3' and 5' ends, (ii) these identifications were supported by multiple sequence differences in both ends, (iii) all differences could be explained by a single recombination event, and (iv) both parent ribotypes occurred multiple times in the clone libraries.
Phylogenetic analyses were performed with algorithms contained in PAUP* (version 4.0b10) (27). Evolutionary-distance trees were inferred by using the minimum evolution optimality criterion with distance correction by HKY85. 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 the full heuristic search option with 1,000 replicates.
Oligonucleotide design, synthesis, and labeling.
Based on the sequence diversity encountered in clone libraries in this study and a previous study (6), a region of the 16S rRNA gene 115 bp in length (nucleotide positions 948 to 1063, E. coli numbering [3]) was identified that (i) contained sufficient variability to differentiate symbiont phylotypes (clades containing several ribotypes) and (ii) was bounded by priming sites conserved among all members of the symbiont clade. For constant denaturant capillary electrophoresis (CDCE), primers 1063r (E. coli numbering) and 948f were synthesized with a 54-bp GC-rich region and a fluorescein isothiocyanate label, respectively, at their 5' ends (Table 1). Primers were synthesized and polyacrylamide gel electrophoresis purified by Integrated DNA Technologies (Coralville, IA).
|
View this table: [in a new window] |
TABLE 1. CDCE primer sequences
|
Internal competitor DNA was prepared from a linearized plasmid containing ribotype RT5 (clone A01) modified by PCR to contain mutations at positions 968 (A
T) and 972 (C
T). The mutations were introduced by using an extended version of the forward amplification primer (948fmut; Table 1) containing mismatches adjacent to the 948f primer used in QPCR-CDCE. The resulting PCR product was cloned as described above, purified, and linearized by overnight digestion with BamHI (Invitrogen, Carlsbad, CA). Concentration was determined by spectrophotometry after purification with a QIAquick PCR purification kit (QIAGEN, Valencia, CA).
For CDCE analysis, PCR mixtures were diluted 50-fold into distilled and deionized water (Millipore, Billerica, MA) and electroinjected for 45 s at 2 µA into a fused silica capillary (75-µm inside diameter) filled with a replaceable 5% linear polyacrylamide gel matrix (Scientific Polymers, Ontario, NY). Electrophoretic separation was performed at constant current (10 µA) and temperature (75.5°C). A 488-nm argon ion laser (ILT, Salt Lake City, UT) served as the light source for fluorescence detection of the fluorescein isothiocyanate label on individual amplicons, and a photomultiplier tube (Oriel, Stratford, CT) was used for fluorescence detection. The output signal was collected with the Workbench Data Acquisition Program (Strawberry Tree, Inc., Sunnyvale, CA).
To determine whether PCR amplification bias occurred during coamplification of multiple templates, linearized plasmids containing target sequences representative of the five phylotypes examined and the internal competitor were mixed in equal proportions (105 copies of each) and a 25-cycle amplification was performed with primers 948f and 1063r-GC. The PCR product was diluted 1:1,000 and reamplified for 15 cycles for a total of three times. For each target-competitor pair, peak area ratios were compared after each round of amplification. After the final round of amplification (a total of >106-fold), no significant difference in amplification efficiency was observed among target sequences relative to the internal competitor, confirming that QPCR-CDCE could be used to estimate the copy numbers of the identified phylotypes by coamplification from gill samples.
Data analysis.
All QPCRs were performed in triplicate, and replicates were analyzed individually by CDCE. Peak areas were measured with AcqKnowledge 2.1 software (Biopac Systems, Santa Barbara, CA) and were used to estimate target sequence copy number with the formula Pi = (CStd)(APi/AStd)/Wg(Ns/Nt), where Pi is the number of small-subunit rRNA gene copies per gram (wet weight) of gill tissue corresponding to the ith phylotype in the initial symbiont population, CStd is the number of copies of internal competitor added to the reaction mixture, APi is the area of the CDCE peak corresponding to an individual bacterial phylotype, AStd is the area of the peak corresponding to the internal competitor, Wg is the wet weight of the gill tissue, Nt is the total nucleic acid extracted from that tissue, and Ns is the nucleic acid used for amplification. Small-subunit rRNA gene copy numbers are reported as means ± standard deviations (n = 3).
Nucleotide sequence accession numbers.
The nucleotide sequences determined in this work have been submitted to GenBank and assigned accession no. DQ272300 to DQ272317.
|
|
|---|
Phylogenetic analyses show that the observed ribotypes form a unique clade within the
subgroup of proteobacteria whose closest relative among cultivated and characterized bacteria is "Saccharophagus degradans" 2-40, a free-living marine heterotrophic bacterium that degrades a broad spectrum of plant polysaccharides (19). This shipworm symbiont clade, in turn, contains several well-supported internal clades. Hereafter, we refer to these distinct clades or phylogenetic types as phylotypes (Fig. 1). Five of these phylotypes, designated P1 to P5, could be distinguished by CDCE-based analysis with a single variable region of the 16S rRNA gene (see below). Sequence variation between groups, among the five phylotypes, ranged from 0.2 to 8.2% (Phylip 3.6b, DNADIST, uncorrected distance) (Table 2), and ribotype diversity observed within phylotypes ranged from 0 to 1.6%.
![]() View larger version (20K): [in a new window] |
FIG. 1. Phylogenetic relationships among 18 unique symbiont ribotypes (RTs) identified in the combined SSL and MSL. Brackets (P1 to P5) indicate clades (phylotypes) that can be distinguished by CDCE with the primers 948f and 1063r-GC. Bold values in parentheses indicate the numbers of clones in the combined libraries belonging to a specific ribotype. The tree was inferred by an evolutionary-distance method (PAUP 4.0, with minimum evolution optimality criterion and HKY85 correction; 1,234 bp were considered in the analysis). Bootstrap proportions are expressed as a percentage of 1,000 replicates and are indicated on the tree for nodes with values >50%. The reference sequences were T. turnerae T7902T AY028398 and "S. degradans" 2-40 AF055269.
|
|
View this table: [in a new window] |
TABLE 2. Distance matrixa comparing sequence divergences among the 18 unique symbiont ribotypes observed in this study and the "S. degradans" reference sequence
|
With a combination of CDCE and QPCR (23, 29), symbiont phylotype abundance was estimated in the gills of 13 individual shipworm specimens. Summing peak intensities from all five phylotypes for 11 specimens (2 were omitted from this calculation because gill weights were not determined) yielded estimates of total symbiont rRNA gene copy abundance ranging from 0.81 x 1010 to 8.99 x 1010 g1 (wet weight) of gill tissue, with a mean of 3.55 x 1010 ± 0.36 x 1010 g1. These numbers cannot be readily translated into symbiont cell abundance estimates because the average numbers of genome copies per symbiont cell and 16S rRNA operons per symbiont genome are not known. Nonetheless, the16S rRNA gene copy abundances observed here fall within the range of cell abundances (109 to 1011 cells g1 [wet weight] of symbiont-containing tissue) estimated for the sulfur-oxidizing symbionts of the marine bivalves Solemya reidi, Bathymodiolus thermophilus, and Calyptogena magnifica and the vestimentiferan tubeworm Riftia pachyptila (4, 25).
Peak patterns in CDCE spectra varied significantly among the gills of the 13 individual specimens examined (Fig. 2B to H), revealing considerable variation in the composition of symbiont populations (Fig. 3A). The number of phylotypes detected per specimen ranged from two to five and included all phylotypes recovered in the combined libraries plus an additional unidentified target sequence (U1) that appeared as a minor component in one specimen (LpS10) (Fig. 3A). In the majority of specimens (11 of 13), two or three symbiont phylotypes were detected (Fig. 3A). The numerically dominant phylotype among all specimens was P3 (Fig. 3C), and the symbiont phylotype most frequently observed was P4 (Fig. 3A). No correlation was observed between rearing location and pattern of symbiont composition.
![]() View larger version (13K): [in a new window] |
FIG. 2. CDCE profiles obtained from gills of the shipworm L. pedicellatus after PCR amplification with primers 948f and 1063r-GC. (A) Combined gills from four specimens used for the MSL but not included among the 13 individual specimens examined. (B to H) CDCE profiles from gills of seven representative specimens of L. pedicellatus, each showing a different combination of identified peaks. The identities of peaks P1 to P5 were confirmed by comigration with PCR products from sequenced clones spiked into samples after PCR and prior to CDCE. U1 in panel E represents an unidentified peak observed in one specimen only.
|
![]() View larger version (52K): [in a new window] |
FIG. 3. Phylotype abundance in the gills of 13 individual specimens of L. pedicellatus as determined by QPCR-CDCE expressed as a percentage of the total CDCE peak area for each specimen. (A) Phylotypes are as follows: P1, diagonal stripes; P2, vertical stripes; P3, horizontal stripes; P4, diamonds; P5, checkering; unidentified peak U1, black. (B) Number of phylotypes detected per specimen ranked in order of decreasing frequency. (C) Average abundance of phylotypes expressed as percentages of the total combined CDCE peak area for gills of 13 specimens of L. pedicellatus.
|
1% of the total population would not be detected. Therefore, it is possible that P4, or any of the other identified symbiont phylotypes, was present in some samples at levels below this detection limit. Thus, in general, the symbiont populations observed in L. pedicellatus are composed of a single dominant phylotype (P1 or P3), a consistent minor phylotype (P4), and in some cases one or more additional minor phylotypes that occur more sporadically. Within each phylotype, there may be considerable ribotype variation and indeed internal clades are present that were not distinguished by the CDCE target sequence used here. The observation that phylotypes P1 and P3 show inverse covariation with one another but not with P4 indicates that different factors determine the abundance and tissue distribution of these three phylotypes. Thus, at least to this extent, they appear to be ecologically and physiologically differentiated. Indeed, the fact that these two phylotypes are sometimes codetected within the same bacteriocytes (6) suggests that their inverse proportions could be explained by competition for the same intracellular niche. It is interesting, however, that either phylotype may be dominant in a given host and that this pattern holds true even when both phylotypes are detected in the same specimen. Moreover, at least in this small sampling of the host population, dominance by either phylotype occurs with similar frequencies. Thus, a simple dominance hierarchy of the phylotypes does not explain the patterns observed. Instead, other currently unknown factors, such as the timing and mode of symbiont transmission, the physiological state and growth dynamics of the host and symbionts, or the physicochemical parameters of the environment, may influence the composition of symbiont populations.
The question remains to what extent intraspecies variation occurs in other symbiont populations. In recent years, multiple coexisting symbionts have been recognized in a small but growing number of marine animal species. For example, two divergent symbiont ribotypes, one a sulfur oxidizer and the other a methane oxidizer, coexist in the bacteriocytes of certain mussels from deep-sea hydrothermal vents and cold-water seeps (7, 16, 32). Two phylogenetically distinct ribotypes, both presumed to derive from sulfur-oxidizing bacteria, have been reported in separate regions of the gill filaments of a thyasirid clam (17). And in reef-building corals, multiple ribotypes of zooxanthellae (eukaryotic algal endosymbionts) have been detected within individual hosts (2). Multiple coexisting symbiont ribotypes have also been observed in naturally occurring and experimentally manipulated extracellular endosymbioses, including light organ symbioses of sepiolid squid (14, 24) and the subcuticular sulfur-oxidizing, sulfate-reducing symbiont communities in an oligochaete worm (11). In each of these cases, the potential for intraspecies variation in the composition of symbiont populations exists and indeed such variation could have adaptive significance for host survival.
The structure of symbiont communities in L. pedicellatus closely resembles that recently observed within free-living marine bacterioplankton (1) and sediment (21) communities in which clusters of closely related co-occurring ribotypes predominate. In these studies, it has been proposed that these microdiverse ribotype clusters (equivalent to the phylotypes observed here) represent important units of ecological differentiation in these communities and that much of the diversity within such clusters may be selectionally neutral (29). It will be important to determine to what extent the phylotypes observed here represent physiologically and ecologically differentiated populations whose distribution correlates with factors such as the developmental state, condition, or environmental fitness of the host.
We thank Luisa A. Marcelino and Aoy Tomita-Mitchell for advice and assistance in developing CDCE protocols.
|
|
|---|
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»