Spatial heterogeneity of crenarchaeal assemblages within mesophilic soil ecosystems as revealed by PCR-single-stranded conformation polymorphism profiling

Microbial ecologists have discovered novel rRNA genes (rDNA) in mesophilic soil habitats worldwide, including sequences that affiliate phylogenetically within the division Crenarchaeota (domain Archaea). To characterize the spatial distribution of crenarchaeal assemblages in mesophilic soil habitats, we profiled amplified crenarchaeal 16S rDNA sequences from diverse soil ecosystems by using PCR-single-stranded-conformation polymorphism (PCR-SSCP) analysis. PCR-SSCP profiles provide a measure of relative microbial diversity in terms of richness (number of different phylotypes as estimated from the number of unique PCR-SSCP peaks) and evenness (abundance of each phylotype as estimated from the relative area under a peak). Crenarchaeal assemblages inhabiting prairie, forest, turf, and agricultural soils were characterized at six sampling locations in southern and central Wisconsin. Phylotype richness was found to be more stable than evenness among triplicate samples collected within 30 cm at each sampling location. Transformation of the PCR-SSCP data by principal-component analysis, followed by statistical testing (analysis of variance [P < 0.0001] and least-significant-difference analysis [alpha = 0.5]), supported the conclusion that each location exhibited a unique profile. To further characterize the spatial distribution of crenarchaeal assemblages at one location, additional soil samples (a total of 30) were collected from agricultural field plots at the Hancock Agricultural Research Station. PCR-SSCP revealed a patchy spatial distribution of crenarchaeal assemblages within and between these plots. This mosaic of crenarchaeal assemblages was characterized by differences in phylotype evenness that could not be correlated with horizontal distance (15 to 30 m) or with depth (0 to 20 cm below the surface). Crenarchaeal 16S rDNA clone libraries were produced and screened for unique SSCP peaks. Clones representing the dominant phylotypes at each location were identified, sequenced, and found to group phylogenetically with sequences in crenarchaeal clade C1b.     

Comparison of Crenarchaeal Consortia Inhabiting the Rhizosphere of Diverse Terrestrial Plants with Those in Bulk Soil in Native Environments

To explore whether the crenarchaeal consortium found in the rhizosphere is distinct from the assemblage of crenarchaeotes inhabiting bulk soil, PCR-single-stranded-conformation polymorphism (PCR-SSCP) profiles were generated for 76 plant samples collected from native environments. Divergent terrestrial plant groups including bryophytes (mosses), lycopods (club mosses), pteridophytes (ferns), gymnosperms (conifers), and angiosperms (seed plants) were collected for this study. Statistical analysis revealed significant differences between rhizosphere and bulk soil PCR-SSCP profiles (Hotelling paired T² test, P < 0.0001), suggesting that a distinct crenarchaeal consortium is associated with plants. In general, phylotype richness increased in the rhizosphere compared to the corresponding bulk soil, although the range of this increase was variable. Examples of a major change in rhizosphere (versus bulk soil) PCR-SSCP profiles were detected for all plant groups, suggesting that crenarchaeotes form associations with phylogenetically diverse plants in native environments. In addition, examples of minor to no detectable difference were found for all terrestrial plant groups, suggesting that crenarchaeal associations with plants are mediated by environmental conditions.     

Comparison of crenarchaeal consortia inhabiting the rhizosphere of diverse terrestrial plants with those in bulk soil in native environments

To explore whether the crenarchaeal consortium found in the rhizosphere is distinct from the assemblage of crenarchaeotes inhabiting bulk soil, PCR-single-stranded-conformation polymorphism (PCR-SSCP) profiles were generated for 76 plant samples collected from native environments. Divergent terrestrial plant groups including bryophytes (mosses), lycopods (club mosses), pteridophytes (ferns), gymnosperms (conifers), and angiosperms (seed plants) were collected for this study. Statistical analysis revealed significant differences between rhizosphere and bulk soil PCR-SSCP profiles (Hotelling paired T(2) test, P < 0.0001), suggesting that a distinct crenarchaeal consortium is associated with plants. In general, phylotype richness increased in the rhizosphere compared to the corresponding bulk soil, although the range of this increase was variable. Examples of a major change in rhizosphere (versus bulk soil) PCR-SSCP profiles were detected for all plant groups, suggesting that crenarchaeotes form associations with phylogenetically diverse plants in native environments. In addition, examples of minor to no detectable difference were found for all terrestrial plant groups, suggesting that crenarchaeal associations with plants are mediated by environmental conditions.     

Bacterial Diversity and Community Structure in an Aerated Lagoon Revealed by Ribosomal Intergenic Spacer Analyses and 16S Ribosomal DNA Sequencing

We investigated the bacterial community structure in an aerated plug-flow lagoon treating pulp and paper mill effluent. For this investigation, we developed a composite method based on analyses of PCR amplicons containing the ribosomal intergenic spacer (RIS) and its flanking partial 16S rRNA gene. Community percent similarity was determined on the basis of RIS length polymorphism. A community succession was evident in the lagoon, indicated by a progressive community transition through seven sample locations. The most abrupt changes in community structure were associated with a temperature change from 39 to 35°C and with increases in dissolved oxygen. The temporal differences in community structure, based on summer and winter samplings, were greater than the spatial differences during either season. Clone libraries of rDNA-RIS amplicons were constructed from each of three summer samples. Among 90 clones analyzed (30 clones from each sample), 56 phylotypes were distinguished by restriction fragment length polymorphism. Indices of phylotype richness, evenness, and diversity all increased in clone libraries from the beginning to the end of the lagoon. A representative clone of each phylotype was phylogenetically analyzed on the basis of its partial 16S rRNA gene sequence (ca. 450 bp). Phylogenetic analysis confirmed the increase in diversity and further indicated increasing richness of bacterial divisions. Pioneers in the community spatial succession appeared to include thermotolerant, microaerophilic methanol-oxidizing bacteria related to the genus Methylobacillus, as well as thermotolerant, microaerophilic nitrogen-fixing bacteria related to the genus Azospirillum.     

Prokaryotic diversity in continuous cropping and rotational cropping soybean soil

In spite of the techniques based on the amplification of 16S rRNA genes (16S rDNA) to compare bacterial communities that are now widely in use in microbial ecology, little is known about the composition of the soybean continuous cropping (CC) and rotational cropping (RC) soil microbial community. To address this, we compared the levels of bacterial community diversity in RC and 5-year CC rhizosphere soil samples. We selected 407 clones in RC and 490 clones in CC for restriction fragment length polymorphism analysis. A total of 123 phylotypes were identified among the 16S rDNA clones, while 78 unique and 21 common phylotypes were identified among the CC soil isolates. Analysis of sequences from a subset of the phylotypes showed that at least 11 bacterial divisions were represented in the clone libraries. The phylotype richness, frequency distribution (evenness), and composition of the two clone libraries were investigated using a variety of diversity indices. Although the analysis of diversity indices and LIBSHUFF comparisons revealed that the compared libraries were not significantly different ( P =0.05) between the RC vs. CC soils, some differences could be observed in terms of specific phyla and groups. We concluded that the group variance was not determined immediately by the cropping system's induction, but was a long-term and slow process.     

Biodiversity, Community Structural Shifts, and Biogeography of Prokaryotes within Antarctic Continental Shelf Sediment

16S ribosomal DNA (rDNA) clone library analysis was conducted to assess prokaryotic diversity and community structural changes within a surficial sediment core obtained from an Antarctic continental shelf area (depth, 761 m) within the Mertz Glacier Polynya (MGP) region. Libraries were created from three separate horizons of the core (0- to 0.4-cm, 1.5- to 2.5-cm, and 20- to 21-cm depth positions). The results indicated that at the oxic sediment surface (depth, 0 to 0.4 cm) the microbial community appeared to be dominated by a small subset of potentially r-strategist (fast-growing, opportunistic) species, resulting in a lower-than-expected species richness of 442 operational taxonomic units (OTUs). At a depth of 1.5 to 2.5 cm, the species richness (1,128 OTUs) was much higher, with the community dominated by numerous gamma and delta proteobacterial phylotypes. At a depth of 20 to 21 cm, a clear decline in species richness (541 OTUs) occurred, accompanied by a larger number of more phylogenetically divergent phylotypes and a decline in the predominance of Proteobacteria. Based on rRNA and clonal abundance as well as sequence comparisons, syntrophic cycling of oxidized and reduced sulfur compounds appeared to be the dominant process in surficial MGP sediment, as phylotype groups putatively linked to these processes made up a large proportion of clones throughout the core. Between 18 and 65% of 16S rDNA phylotypes detected in a wide range of coastal and open ocean sediments possessed high levels of sequence similarity (>95%) with the MGP sediment phylotypes, indicating that many sediment prokaryote phylotype groups defined in this study are ubiquitous in marine sediment.     

Analysis of the Rumen Bacterial Diversity under two Different Diet Conditions using Denaturing Gradient Gel Electrophoresis,Random Sequencing,and Statistical Ecology Approaches

Bacterial community structure and diversity in the rumen of steers in conditions of hay and corn diets was assessed by in vitro retrieval and analysis of the variable region (V3) of 16S rDNA. Two types of libraries were generated in this study: DGGE libraries, which further were analysed by excising, reamplification, and sequencing, and random shotgun sequence libraries. Phylogenetic and sequence similarity analyses of the resultant 68 clone sequences in DGGE libraries revealed the presence of 42 operational taxonomic units (OTUs) or phylotypes defined as having more than 97% of sequence similarity. One hundred and thirty four clone sequences in shotgun libraries were clustered into 72 phylotypes. The phylotype similarity, diversity, richness, and evenness in these libraries were estimated using a variety of diversity indices. In relation to diet, the corn-fed animals displayed more diverse and rich bacterial populations, which were mostly contributed by CFB-related phylotypes. Proteobacteria were also numerically prevalent on this diet (27%) but were represented by a few phylotypes thus diminishing the overall diversity and species richness values. On hay diet, the principal contributors to general diversity and species richness appeared to be low-G + C gram-positives. Although the ruminal Treponemaes were encountered only in hay-fed animals, their impact on species diversity on hay diet was low because of the limited number of phylotypes.     

Characterization of the Dominant and Rare Members of a Young Hawaiian Soil Bacterial Community with Small-Subunit Ribosomal DNA Amplified from DNA Fractionated on the Basis of Its Guanine and Cytosine Composition

The small-subunit ribosomal DNA (rDNA) diversity was found to be very high in a Hawaiian soil community that might be expected to have lower diversity than the communities in continental soils because the Hawaiian soil is geographically isolated and only 200 years old, is subjected to a constant climate, and harbors low plant diversity. Since an underlying community structure could not be revealed by analyzing the total eubacterial rDNA, we first fractionated the DNA on the basis of guanine-plus-cytosine (G+C) content by using bis-benzimidazole and equilibrium centrifugation and then analyzed the bacterial rDNA amplified from a fraction with a high biomass (63% G+C fraction) and a fraction with a low biomass (35% G+C fraction). The rDNA clone libraries were screened by amplified rDNA restriction analysis to determine phylotype distribution. The dominant biomass reflected by the 63% G+C fraction contained several dominant phylotypes, while the community members that were less successful (35% G+C fraction) did not show dominance but there was a very high diversity of phylotypes. Nucleotide sequence analysis revealed taxa belonging to the groups expected for the G+C contents used. The dominant phylotypes in the 63% G+C fraction were members of the Pseudomonas, Rhizobium-Agrobacterium, and Rhodospirillum assemblages, while all of the clones sequenced from the 35% G+C fraction were affiliated with several Clostridium assemblages. The two-step rDNA analysis used here uncovered more diversity than can be detected by direct rDNA analysis of total community DNA. The G+C separation step is also a way to detect some of the less dominant organisms in a community.     

DNA-based and culture-based characterization of a hydrocarbon-degrading consortium enriched from Arctic soil

A hydrocarbon-degrading consortium was enriched from fuel-contaminated soil from the northeastern tip of Ellesmere Island (82 degrees 30'N, 62 degrees 19'W). The enrichment culture was grown on Jet A-1 fuel at 7 degrees C. Bacterial 16S RNA gene (rDNA) fragments were amplified by polymerase chain reaction (PCR) from members of the above consortium and cloned into a plasmid vector. Partial sequences (approximately 500 bp) were determined for 29 randomly selected rDNA clones. The majority of sequences were most similar to the corresponding rDNA sequences of Rhodococcus erythropolis (15 sequences), Sphingomonas spp. (six sequences), and Pseudomonas synxantha (four sequences). Amplified ribosomal DNA restriction analysis confirmed that a larger set of 50 clones had frequencies of the three phylotypes similar to those above. Phylotype-specific PCR assays were developed and validated for the above three phylotypes. The consortium was plated and grown on Jet A-1 fuel vapors, and randomly selected isolated colonies were screened with the above PCR assays. Of 17 colonies, six matched the Rhodococcus phylotype, and three matched the Pseudomonas phylotype. A representative strain of each phylotype was physiologically characterized. Both isolates grew on alkanes at low temperature and had general characteristics consistent with their respective phylotypes. During growth of the consortium, the three phylotype populations were monitored by a most probable number PCR assay. All three phylotypes were detected, but their relative abundance was not consistent with that of the phylotypes in the clone library. The relative abundance of all three phylotypes changed substantially during long-term incubation of the consortium. The DNA-based approach used identified phylotypes consistently present in the consortium, but it failed to predict the relative abundance of their populations.     

Bacterial diversity and community structure in an aerated lagoon revealed by ribosomal intergenic spacer analyses and 16S ribosomal DNA sequencing

We investigated the bacterial community structure in an aerated plug-flow lagoon treating pulp and paper mill effluent. For this investigation, we developed a composite method based on analyses of PCR amplicons containing the ribosomal intergenic spacer (RIS) and its flanking partial 16S rRNA gene. Community percent similarity was determined on the basis of RIS length polymorphism. A community succession was evident in the lagoon, indicated by a progressive community transition through seven sample locations. The most abrupt changes in community structure were associated with a temperature change from 39 to 35 degrees C and with increases in dissolved oxygen. The temporal differences in community structure, based on summer and winter samplings, were greater than the spatial differences during either season. Clone libraries of rDNA-RIS amplicons were constructed from each of three summer samples. Among 90 clones analyzed (30 clones from each sample), 56 phylotypes were distinguished by restriction fragment length polymorphism. Indices of phylotype richness, evenness, and diversity all increased in clone libraries from the beginning to the end of the lagoon. A representative clone of each phylotype was phylogenetically analyzed on the basis of its partial 16S rRNA gene sequence (ca. 450 bp). Phylogenetic analysis confirmed the increase in diversity and further indicated increasing richness of bacterial divisions. Pioneers in the community spatial succession appeared to include thermotolerant, microaerophilic methanol-oxidizing bacteria related to the genus Methylobacillus, as well as thermotolerant, microaerophilic nitrogen-fixing bacteria related to the genus Azospirillum.     

Niche Partition of Bacteriovorax Operational Taxonomic Units Along Salinity and Temporal Gradients in the Chesapeake Bay Reveals Distinct Estuarine Strains

The predatory Bacteriovorax are Gram-negative bacteria ubiquitous in saltwater systems that prey upon other Gram-negative bacteria in a similar manner to the related genus Bdellovibrio. Among the phylogenetically defined clusters of Bacteriovorax, cluster V has only been isolated from estuaries suggesting that it may be a distinct estuarine phylotype. To assess this hypothesis, the spatial and temporal distribution of cluster V and other Bacteriovorax phylogenetic assemblages along the salinity gradient of Chesapeake Bay were determined. Cluster V was expected to be found in significantly greater numbers in low to moderate salinity waters compared to high salinity areas. The analyses of water and sediment samples from sites in the bay revealed cluster V to be present at the lower salinity and not high salinity sites, consistent with it being an estuarine phylotype. Cluster IV had a similar distribution pattern and may also be specifically adapted to estuaries. While the distribution of clusters V and IV were similar for salinity, they were distinct on temperature gradients, being found in cooler and in warmer temperatures, respectively. The differentiation of phylotype populations along the salinity and temporal gradients in Chesapeake Bay revealed distinct niches inhabited by different phylotypes of Bacteriovorax and unique estuarine phylotypes.     

Cultivation of Mesophilic Soil Crenarchaeotes in Enrichment Cultures from Plant Roots

Because archaea are generally associated with extreme environments, detection of nonthermophilic members belonging to the archaeal division Crenarchaeota over the last decade was unexpected; they are surprisingly ubiquitous and abundant in nonextreme marine and terrestrial habitats. Metabolic characterization of these nonthermophilic crenarchaeotes has been impeded by their intractability toward isolation and growth in culture. From studies employing a combination of cultivation and molecular phylogenetic techniques (PCR-single-strand conformation polymorphism, sequence analysis of 16S rRNA genes, fluorescence in situ hybridization, and real-time PCR), we present evidence here that one of the two dominant phylotypes of Crenarchaeota that colonizes the roots of tomato plants grown in soil from a Wisconsin field is selectively enriched in mixed cultures amended with root extract. Clones recovered from enrichment cultures were found to group phylogenetically with sequences from clade C1b.A1. This work corroborates and extends our recent findings, indicating that the diversity of the crenarchaeal soil assemblage is influenced by the rhizosphere and that mesophilic soil crenarchaeotes are found associated with plant roots, and provides the first evidence for growth of nonthermophilic crenarchaeotes in culture.     

Soil Bacterial Community Shift Correlated with Change from Forest to Pasture Vegetation in a Tropical Soil

The change in vegetative cover of a Hawaiian soil from forest to pasture led to significant changes in the composition of the soil bacterial community. DNAs were extracted from both soil habitats and compared for the abundance of guanine-plus-cytosine (G+C) content, by analysis of abundance of phylotypes of small-subunit ribosomal DNA (SSU rDNA) amplified from fractions with 63 and 35% G+C contents, and by phylogenetic analysis of the dominant rDNA clones in the 63% G+C content fraction. All three methods showed differences between the forest and pasture habitats, providing evidence that vegetation had a strong influence on microbial community composition at three levels of taxon resolution. The forest soil DNA had a peak in G+C content of 61%, while the DNA of the pasture soil had a peak in G+C content of 67%. None of the dominant phylotypes found in the forest soil were detected in the pasture soil. For the 63% G+C fraction SSU rDNA sequence analysis of the three most dominant members revealed that their phyla changed from Fibrobacter and Syntrophomonas assemblages in the forest soil to Burkholderia and Rhizobium–Agrobacterium assemblages in the pasture soil.     

Assessing the diversity of arbuscular mycorrhizal fungi in semiarid shrublands dominated by Artemisia tridentata ssp. wyomingensis

Variation in the abiotic environment and host plant preferences can affect the composition of arbuscular mycorrhizal (AMF) assemblages. This study analyzed the AMF taxa present in soil and seedlings of Artemisia tridentata ssp. wyomingensis collected from sagebrush steppe communities in southwestern Idaho, USA. Our aims were to determine the AMF diversity within and among these communities and the extent to which preferential AMF–plant associations develop during seedling establishment. Mycorrhizae were identified using molecular methods following DNA extraction from field and pot culture samples. The extracted DNA was amplified using Glomeromycota specific primers, and identification of AMF was based on phylogenetic analysis of sequences from the large subunit-D2 rDNA region. The phylogenetic analyses revealed seven phylotypes, two within the Claroideoglomeraceae and five within the Glomeraceae. Four phylotypes clustered with known species including Claroideoglomus claroideum, Rhizophagus irregularis, Glomus microaggregatum, and Funneliformis mosseae. The other three phylotypes were similar to several published sequences not included in the phylogenetic analysis, but all of these were from uncultured and unnamed glomeromycetes. Pairwise distance analysis revealed some phylotypes with high genetic variation. The most diverse was the phylotype that included R. irregularis, which contained sequences showing pairwise differences up to 12 %. Most of the diversity in AMF sequences occurred within sites. The smaller genetic differentiation detected among sites was correlated with differences in soil texture. In addition, multiplication in pot cultures led to differentiation of AMF communities. Comparison of sequences obtained from the soil with those from A. tridentata roots revealed no significant differences between the AMF present in these samples. Overall, the sites sampled were dominated by cosmopolitan AMF taxa, and young seedlings of A. tridentata ssp. wyomingensis were colonized in relation to the abundance of these taxa in the soil.     

Biodiversity,community structural shifts,and biogeography of prokaryotes within Antarctic continental shelf sediment

16S ribosomal DNA (rDNA) clone library analysis was conducted to assess prokaryotic diversity and community structural changes within a surficial sediment core obtained from an Antarctic continental shelf area (depth, 761 m) within the Mertz Glacier Polynya (MGP) region. Libraries were created from three separate horizons of the core (0- to 0.4-cm, 1.5- to 2.5-cm, and 20- to 21-cm depth positions). The results indicated that at the oxic sediment surface (depth, 0 to 0.4 cm) the microbial community appeared to be dominated by a small subset of potentially r-strategist (fast-growing, opportunistic) species, resulting in a lower-than-expected species richness of 442 operational taxonomic units (OTUs). At a depth of 1.5 to 2.5 cm, the species richness (1,128 OTUs) was much higher, with the community dominated by numerous gamma and delta proteobacterial phylotypes. At a depth of 20 to 21 cm, a clear decline in species richness (541 OTUs) occurred, accompanied by a larger number of more phylogenetically divergent phylotypes and a decline in the predominance of Proteobacteria. Based on rRNA and clonal abundance as well as sequence comparisons, syntrophic cycling of oxidized and reduced sulfur compounds appeared to be the dominant process in surficial MGP sediment, as phylotype groups putatively linked to these processes made up a large proportion of clones throughout the core. Between 18 and 65% of 16S rDNA phylotypes detected in a wide range of coastal and open ocean sediments possessed high levels of sequence similarity (>95%) with the MGP sediment phylotypes, indicating that many sediment prokaryote phylotype groups defined in this study are ubiquitous in marine sediment.     

Community Dynamics of Arbuscular Mycorrhizal Fungi in High-Input and Intensively Irrigated Rice Cultivation Systems

Application of a mycorrhizal inoculum could be one way to increase the yield of rice plants and reduce the application of fertilizer. We therefore studied arbuscular mycorrhizal fungi (AMF) in the roots of wetland rice (Oryza sativa L.) collected at the seedling, tillering, heading, and ripening stages in four paddy wetlands that had been under a high-input and intensively irrigated rice cultivation system for more than 20 years. It was found that AMF colonization was mainly established in the heading and ripening stages. The AMF community structure was characterized in rhizosphere soils and roots from two of the studied paddy wetlands. A fragment covering the partial small subunit (SSU), the whole internal transcribed spacer (ITS), and the partial large subunit (LSU) rRNA operon regions of AMF was amplified, cloned, and sequenced from roots and soils. A total of 639 AMF sequences were obtained, and these were finally assigned to 16 phylotypes based on a phylogenetic analysis, including 12 phylotypes from Glomeraceae, one phylotype from Claroideoglomeraceae, two phylotypes from Paraglomeraceae, and one unidentified phylotype. The AMF phylotype compositions in the soils were similar between the two surveyed sites, but there was a clear discrepancy between the communities obtained from root and soil. The relatively high number of AMF phylotypes at the surveyed sites suggests that the conditions are suitable for some species of AMF and that they may have an important function in conventional rice cultivation systems. The species richness of root-colonizing AMF increased with the growth of rice, and future studies should consider the developmental stages of this crop in the exploration of AMF function in paddy wetlands.     

Cultivation of mesophilic soil crenarchaeotes in enrichment cultures from plant roots

Because archaea are generally associated with extreme environments, detection of nonthermophilic members belonging to the archaeal division Crenarchaeota over the last decade was unexpected; they are surprisingly ubiquitous and abundant in nonextreme marine and terrestrial habitats. Metabolic characterization of these nonthermophilic crenarchaeotes has been impeded by their intractability toward isolation and growth in culture. From studies employing a combination of cultivation and molecular phylogenetic techniques (PCR-single-strand conformation polymorphism, sequence analysis of 16S rRNA genes, fluorescence in situ hybridization, and real-time PCR), we present evidence here that one of the two dominant phylotypes of Crenarchaeota that colonizes the roots of tomato plants grown in soil from a Wisconsin field is selectively enriched in mixed cultures amended with root extract. Clones recovered from enrichment cultures were found to group phylogenetically with sequences from clade C1b.A1. This work corroborates and extends our recent findings, indicating that the diversity of the crenarchaeal soil assemblage is influenced by the rhizosphere and that mesophilic soil crenarchaeotes are found associated with plant roots, and provides the first evidence for growth of nonthermophilic crenarchaeotes in culture.     

Coexistence of Bacterial Sulfide Oxidizers, Sulfate Reducers, and Spirochetes in a Gutless Worm (Oligochaeta) from the Peru Margin

Olavius crassitunicatus is a small symbiont-bearing worm that occurs at high abundance in oxygen-deficient sediments in the East Pacific Ocean. Using comparative 16S rRNA sequence analysis and fluorescence in situ hybridization, we examined the diversity and phylogeny of bacterial symbionts in two geographically distant O. crassitunicatus populations (separated by 385 km) on the Peru margin (water depth, ~300 m). Five distinct bacterial phylotypes co-occurred in all specimens from both sites: two members of the γ-Proteobacteria (Gamma 1 and 2 symbionts), two members of the δ-Proteobacteria (Delta 1 and 2 symbionts), and one spirochete. A sixth phylotype belonging to the δ-Proteobacteria (Delta 3 symbiont) was found in only one of the two host populations. Three of the O. crassitunicatus bacterial phylotypes are closely related to symbionts of other gutless oligochaete species; the Gamma 1 phylotype is closely related to sulfide-oxidizing symbionts of Olavius algarvensis, Olavius loisae, and Inanidrilus leukodermatus, the Delta 1 phylotype is closely related to sulfate-reducing symbionts of O. algarvensis, and the spirochete is closely related to spirochetal symbionts of O. loisae. In contrast, the Gamma 2 phylotype and the Delta 2 and 3 phylotypes belong to novel lineages that are not related to other bacterial symbionts. Such a phylogenetically diverse yet highly specific and stable association in which multiple bacterial phylotypes coexist within a single host has not been described previously for marine invertebrates.