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

 Previous Article  |  Next Article 

Applied and Environmental Microbiology, August 2008, p. 4973-4977, Vol. 74, No. 15
0099-2240/08/$08.00+0     doi:10.1128/AEM.02915-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Diversity of Uncultured Epsilonproteobacteria from Terrestrial Sulfidic Caves and Springs{triangledown}

Megan L. Porter1 and Annette Summers Engel2*

University of Maryland Baltimore County, Department of Biological Sciences, Baltimore, Maryland 21250,1 Louisiana State University, Department of Geology and Geophysics, Baton Rouge, Louisiana 708032

Received 24 December 2007/ Accepted 28 May 2008


arrow
ABSTRACT
 
This study expands the phylogenetic diversity of Epsilonproteobacteria using the 16S rRNA gene framework. Of the 73 lineages defined by sequence similarities at or greater than 99%, most were found at only one site. In contrast, eight lineages were retrieved from sites spanning geographic distances from 1,000 to >10,000 km.


arrow
INTRODUCTION
 
A recent review of the Epsilonproteobacteria taxonomic class revealed that 16S rRNA gene sequences retrieved from extreme environments comprised a large clade of mostly uncultivated microbes, provisionally designated the family "Thiovulgaceae" (3). While the marine representatives of the "Thiovulgaceae" represent a diverse assemblage (see, e.g., references 11, 13, 15, 18, 22, 23, and 30), nearly 80% of the terrestrial sequences in this family originated from just one site: Lower Kane Cave (LKC), Wyoming (3, 7, 8). LKC represented the first nonmarine, natural system demonstrably influenced by the activity of Epsilonproteobacteria (7). The Epsilonproteobacteria have been shown to contribute to ecosystem function via chemolithoautotrophic metabolisms (8) and play a significant role in rock dissolution and porosity formation (9). Although six Epsilonproteobacteria operational taxonomic units (OTUs) from LKC have been described (7, 8), this diversity was dominated by the two lineages: LKC group I (LKC-GI) and LKC group II (LKC-GII). Here, we extend the diversity of terrestrial Epsilonproteobacteria by investigating the 16S rRNA sequence diversity in 15 additional sulfidic cave and spring sites (Table 1) and demonstrate that these lineages are found in novel geographic locations.


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

 
TABLE 1. Epsilonproteobacterial diversity associated with sampling sites described in this and previous studiesa

16S rRNA gene sequences were amplified from DNA extracted from microbial mats collected at each site. Universal bacterial primers were used for cave sites (7) and newly developed lineage-specific primers for each LKC group (7) for the spring site DNA. Primer eps59f (forward, 5'-AGTCGAACGATGAGAGGA-3') was used for LKC-GI and eps174f (forward, 5'-CCCCATACTCCTTCTCAT-3') for LKC-GII; both of these primers were paired with a universal bacterial reverse primer. From these amplification products, 15 clone libraries were constructed and screened, and 2 to 57 clones were sequenced from each library using previously described methods (7, 8). Sequences were submitted to the RDP CHECK- CHIMERA program (20).

For phylogenetic analyses, 207 new epsilonproteobacterial sequences were combined with 139 previously published epsilonproteobacterial sequences (62 sulfidic cave and spring sequences; 77 from other habitats) and 16 outgroup taxa. Following sequence alignment (29) and refinement (6), GBlocks v0.91b (4) was used to eliminate alignment regions where there were sequence data for only a small proportion of the taxa. The full alignment file can be downloaded at http://geol.lsu.edu/Faculty/Engel/epsilon.htm. A maximum likelihood phylogeny with branch support from 100 bootstrap replicates was constructed (14) using the best-fit model (25, 26). Delineation of sequence clusters at the 99% sequence similarity predicted for species level relationships (see, e.g., references 16, 17, and 28) was done in DOTUR using the furthest-neighbor clustering algorithm (27). Species level OTUs were designated if DOTUR similarity groups were concordant with a monophyletic clade.

Campbell et al. (3) systematically evaluated 16S rRNA gene sequences derived from environmental sources and found that marine and terrestrial evolutionary lineages were phylogenetically distinct from each other, although there were far fewer terrestrially derived sequences than sequences from marine settings. The sequences recovered from the cave and spring locations represent extensive diversity within the Epsilonproteobacteria, with representatives from all of the major families except the Hydrogenimonaceae (order Campylobacterales) and Nautiliaceae (order Nautiliales) (Fig. 1). We recovered the same major clades designated by Campbell et al. (3), with the majority of new sequences belonging to the provisional family "Thiovulgaceae" (Fig. 1). These results expand the diversity and geographic ranges of the clades and add terrestrially derived sequences to marine clades I and II (3) (Fig. 1). From the current data set, there were no apparent relationships among phylogeny, habitat type (e.g., light versus dark), and geochemistry (data not shown) (Fig. 1; Table 1).


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

 
FIG. 1. Maximum likelihood phylogeny of the newly characterized Epsilonproteobacteria terrestrial sequences, rooted to basal epsilonproteobacterial sequences (e.g., Caminibacter, AY691429 and AY691430; Hydrogenimonas, AB105048; Lebetimonas, AB167820; Nautilia, AF357197 and AJ575809; Nitratiruptor, AB175499; Thioreductor, AB175498), other proteobacterial sequences (e.g., Alphaproteobacteria, AY180103 and D84529; Deltaproteobacteria, Y11560 and X95181; Gammaproteobacteria, AJ543327 and AY776275), and Aquifex aeolicus (AJ309733) and Thermotoga maritima (AJ401021) sequences (outgroups not shown). Accession numbers for previously acquired sequences are included in brackets. Nodal supports from 100 bootstrap replicates are indicated on branches; bootstrap proportions less than 70% are not shown. Species level OTUs are based on monophyletic clades identified by DOTUR clustering at the 99% similarity level and are indicated as follows: circles represent OTUs consisting of sequences from a single location, squares represent OTUs consisting of sequences from two locations found within the same region, and stars represent OTUs consisting of sequences from two or more sites located on different continents. For clarity, monophyletic clades (i.e., a group of sequences containing an ancestor and all of its descendants) consisting of sequences from a single location are collapsed to a single branch; some of these sequence groups consist of more than one OTU (as defined by sequences with ≥99% similarity), which is indicated by more than one symbol per branch. For collapsed clades, the number of sequences is indicated in parentheses behind each location name. Abbreviations correspond to sampling sites listed in Table 1. The four main "Thiovulgaceae" clades from Campbell et al. (3) are indicated.

Based on OTU designations, ~50% of the lineages recovered were novel and 70% of the OTUs were represented by singleton and doubleton sequence groups. Based on this high number of novel OTUs, it is highly likely that additional terrestrial epsilonproteobacterial diversity remains to be characterized. Future research should continue to explore these sites, as the goal of this initial study was not to exhaustively sample all sites using all possible primer sets. Nevertheless, by using the lineage-specific approach, we may have uncovered epsilonproteobacterial diversity that would have otherwise gone undetected by using domain-specific primers (see, e.g., reference 24).

It was expected that greater geographic separation between the sampled sites would correspond to increased genetic distance among the sequences from each site (see, e.g., references 12, 16, 17, and 24). However, this hypothesis was not supported. Only eight OTUs comprised sequences retrieved from more than one site (Fig. 1). One of these OTUs comprised sequences from more than two sites, specifically, nine locations spanning >10,000 km. This large OTU was also phylogenetically affiliated with LKC-GII and was represented by ≤1% sequence divergence (Fig. 2). A recent study of the Frasassi Caves in Italy also identified sequences belonging to this OTU (e.g., clone WM35) (19), further supporting the cosmopolitan distribution of this group (21). Interestingly, sequences affiliated with LKC-GI were not found in any of the newly investigated habitats in this study.


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

 
FIG. 2. Average pairwise sequence distances for each of the eight OTUs containing representatives from more than one site plotted against the average geographic distance among all of the sites within the OTU. The phylogenetic distances were calculated using the PHYLIP v3.6 DNADIST module (10), which was used as DOTUR input to determine OTUs. Error bars represent ± 1 standard deviation.

It is clear that the Epsilonproteobacteria are key players in biogeochemical cycling in marine habitats (see, e.g., references 3 and 23). Although it is possible that these organisms may also be significant players in terrestrial systems, to our knowledge, this is the first study to evaluate whether or not the Epsilonproteobacteria are widely distributed in terrestrial habitats. Our findings revealed that at least the two LKC lineages have contrasting patterns of geographic distribution, even within a limited 16S rRNA framework. Because we found no correlation between phylogeny and geochemical conditions or habitat type (data not shown), the disparate spatial distributions may be due to differing biogeographic, geologic, and/or hydrostratigraphic mechanisms (see, e.g., reference 3). At this time, however, any further discussion would be highly speculative because we currently do not understand the mechanisms for microbial dispersal and colonization in the terrestrial subsurface. Moreover, because the current range of sample locations is restricted to the Northern Hemisphere, future work from both noncave or karst systems and Southern Hemisphere sites will undoubtedly provide new evidence to test hypotheses related to diversity driven by geographic isolation (see, e.g., references 5, 12, and 24).


arrow
Nucleotide sequence accession numbers.
 
The 207 new epsilonproteobacterial sequences obtained in this study have been assigned GenBank accession numbers DQ295539 to DQ295541, DQ295543 to DQ295720, and DQ295722 to DQ295747.


arrow
ACKNOWLEDGMENTS
 
We thank the many field assistants who helped to sample the caves and springs over the years. B. J. Campbell, D. E. Northup, C. J. Schulz, and P. C. Bennett contributed to discussions regarding the implications of the results.

This research was initiated with support from the National Science Foundation LExEn grant (EAR-0085576) to P. C. Bennett at the University of Texas at Austin. Continuing work was partly supported by the College of Basic Sciences at Louisiana State University, the Louisiana Board of Regents (NSF/LEQSF 2005-Pfund-04) to sample in the Frasassi Caves, Italy, and the National Science Foundation (DEB-0640835).


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Louisiana State University, Department of Geology and Geophysics, E235 Howe-Russell Geoscience Complex, Baton Rouge, LA 70803. Phone: (225) 578-2469. Fax: (225) 578-2302. E-mail: aengel{at}lsu.edu Back

{triangledown} Published ahead of print on 6 June 2008. Back


arrow
REFERENCES
 
    1
  1. Angert, E. R., D. E. Northup, A. L. Reysenbach, A. S. Peek, B. M. Goebel, and N. R. Pace. 1998. Molecular phylogenetic analysis of a bacterial community in Sulphur River, Parker Cave, Kentucky. Am. Mineral. 83:1583-1592.[Abstract]
  2. 2
  3. Barton, H. A., and F. Luiszer. 2005. Microbial metabolic structure in a sulfidic cave hot spring: potential mechanisms of biospeleogenesis. J. Cave Karst Stud. 67:28-38.
  4. 3
  5. Campbell, B. J., A. S. Engel, M. L. Porter, and K. Takai. 2006. The versatile {varepsilon}-proteobacteria: key players in the sulphur cycle. Nat. Rev. Microbiol. 4:458-468.[CrossRef][Medline]
  6. 4
  7. Castresana, J. 2000. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17:540-552.[Abstract/Free Full Text]
  8. 5
  9. Cho, J.-C., and J. M. Tiedje. 2000. Biogeography and degree of endemicity of fluorescent Pseudomonas strains in soil. Appl. Environ. Microbiol. 66:5448-5456.[Abstract/Free Full Text]
  10. 6
  11. Edgar, R. C. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32:1792-1797.[Abstract/Free Full Text]
  12. 7
  13. Engel, A. S., N. Lee, M. L. Porter, L. A. Stern, P. C. Bennett, and M. Wagner. 2003. Filamentous "Epsilonproteobacteria" dominate microbial mats from sulfidic cave springs. Appl. Environ. Microbiol. 69:5503-5511.[Abstract/Free Full Text]
  14. 8
  15. Engel, A. S., M. L. Porter, L. A. Stern, S. Quinlan, and P. C. Bennett. 2004. Bacterial diversity and ecosystem function of filamentous microbial mats from aphotic (cave) sulfidic springs dominated by chemolithoautotrophic "Epsilonproteobacteria." FEMS Microbiol. Ecol. 51:31-53.[CrossRef][Medline]
  16. 9
  17. Engel, A. S., L. A. Stern, and P. C. Bennett. 2004. Microbial contributions to cave formation: new insights into sulfuric acid speleogenesis. Geology 32:369-372.[Abstract/Free Full Text]
  18. 10
  19. Felsenstein, J. 1989. PHYLIP—Phylogeny Inference Package (v 3.2). Cladistics 5:164-166.
  20. 11
  21. Goffredi, S., S. Johnson, and R. Vrijenhoek. 2007. Genetic diversity and potential function of microbial symbionts associated with newly discovered species of Osedax polychaete worms. Appl. Environ. Microbiol. 73:2314-2323.[Abstract/Free Full Text]
  22. 12
  23. Green, J., and B. J. M. Bohannan. 2006. Spatial scaling of microbial biodiveristy. Trends Ecol. Evol. 21:501-507.[CrossRef][Medline]
  24. 13
  25. Grote, J., M. Labrenz, D. Pfeiffer, G. Jost, and K. J.ürgens. 2007. Quantitative distributions of Epsilonproteobacteria and a Sulfurimonas subgroup in pelagic redoxclines of the central Baltic Sea. Appl. Environ. Microbiol. 73:7155-7161.[Abstract/Free Full Text]
  26. 14
  27. Guindon, S., and O. Gascuel. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52:696-704.[Abstract/Free Full Text]
  28. 15
  29. Hirayama, H., M. Sunamura, K. Takai, T. Nunoura, T. Noguchi, H. Oida, Y. Furushima, H. Yamamoto, T. Oomori, and K. Horikoshi. 2007. Culture-dependent and -independent characterization of microbial communities associated with a shallow submarine hydrothermal system occurring within a coral reef off Taketomi Island, Japan. Appl. Environ. Microbiol. 73:7642-7656.[Abstract/Free Full Text]
  30. 16
  31. Horner-Devine, M. C., K. M. Carney, and B. J. M. Bohannan. 2004. An ecological perspective on bacterial biodiversity. Proc. R. Soc. Lond. Ser. B Biol. Sci. 271:113-122.[Medline]
  32. 17
  33. Horner-Devine, M. C., M. Lage, J. B. Hughes, and B. J. M. Bohannan. 2004. A taxa-area relationship for bacteria. Nature 432:750-753.[CrossRef][Medline]
  34. 18
  35. Huber, J., D. M. Welch, H. Morrison, S. Huse, P. Neal, D. Butterfield, and M. Sogin. 2007. Microbial population structures in the deep marine biosphere. Science 318:97-100.[Abstract/Free Full Text]
  36. 19
  37. Macalady, J. L., E. H. Lyon, B. Koffman, L. K. Albertson, K. Meyer, S. Galdenzi, and S. Mariani. 2006. Dominant microbial populations in limestone-corroding stream biofilms, Frasassi cave system, Italy. Appl. Environ. Microbiol. 72:5596-5609.[Abstract/Free Full Text]
  38. 20
  39. Maidak, B. L., J. R. Cole, T. G. Lilburn, C. T. Parker, Jr., P. R. Saxman, R. J. Farris, et al. 2001. The RDP-II (Ribosomal Database Project). Nucleic Acids Res. 29:173-174.[Abstract/Free Full Text]
  40. 21
  41. Massana, R., E. F. DeLong, and C. Pedros-Alio. 2000. A few cosmopolitan phylotypes dominate planktonic archaeal assemblages in widely different oceanic provinces. Appl. Environ. Microbiol. 66:1777-1787.[Abstract/Free Full Text]
  42. 22
  43. Moussard, H., E. Corre, M. A. Cambon-Bonavita, Y. Fouquet, and C. Jeanthon. 2006. Novel uncultured Epsilonproteobacteria dominate a filamentous sulphur mat from the 13 degrees N hydrothermal vent field, East Pacific Rise. FEMS Microbiol. Ecol. 58:449-463.[CrossRef][Medline]
  44. 23
  45. Nakagawa, S., Y. Takaki, S. Shimamura, A.-L. Reysenbach, K. Takai, and K. Horikoshi. 2007. Deep-sea vent epsilon-proteobacterial genomes provide insights into emergence of pathogens. Proc. Natl. Acad. Sci. USA 104:12146-12150.[Abstract/Free Full Text]
  46. 24
  47. Papke, R. T., N. B. Ramsing, M. M. Bateson, and D. M. Ward. 2003. Geographical isolation in hot spring cyanobacteria. Environ. Microbiol. 5:650-659.[CrossRef][Medline]
  48. 25
  49. Posada, D., and T. R. Buckley. 2004. Model selection and model averaging in phylogenetics: advantages of the AIC and Bayesian approaches over likelihood ratio tests. Syst. Biol. 53:793-808.[Abstract/Free Full Text]
  50. 26
  51. Posada, D., and K. A. Crandall. 1998. Modeltest: testing the model of DNA substitution. Bioinformatics 14:817-818.[Abstract/Free Full Text]
  52. 27
  53. Schloss, P. D., and J. Handelsman. 2005. Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl. Environ. Microbiol. 71:1501-1506.[Abstract/Free Full Text]
  54. 28
  55. Stackebrandt, E., and B. M. Goebel. 1994. A place for DNA-DNA reassociation and 16S ribosomal-RNA sequence-analysis in the present species definition in bacteriology. Int. J. Syst. Evol. Microbiol. 44:846-849.[Abstract/Free Full Text]
  56. 29
  57. Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The Clustal_X Windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25:4876-4882.
  58. 30
  59. Yakimov, M., V. La Cono, R. Denaro, G. D'Auria, F. Decembrini, K. Timmis, P. Golyshin, and L. Giuliano. 2007. Primary producing prokaryotic communities of brine, interface and seawater above the halocline of deep anoxic Lake L'Atalante, eastern Mediterranean Sea. ISME J. 1:743-755.[CrossRef][Medline]


Applied and Environmental Microbiology, August 2008, p. 4973-4977, Vol. 74, No. 15
0099-2240/08/$08.00+0     doi:10.1128/AEM.02915-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.





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