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
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

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.

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.
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).
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.
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).

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.

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).

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 
Published ahead of print on 6 June 2008. 

REFERENCES
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 - 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.
3 - Campbell, B. J., A. S. Engel, M. L. Porter, and K. Takai. 2006. The versatile
-proteobacteria: key players in the sulphur cycle. Nat. Rev. Microbiol. 4:458-468.[CrossRef][Medline] 4 - 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]
5 - 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]
6 - Edgar, R. C. 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32:1792-1797.[Abstract/Free Full Text]
7 - 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]
8 - 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]
9 - 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]
10 - Felsenstein, J. 1989. PHYLIP—Phylogeny Inference Package (v 3.2). Cladistics 5:164-166.
11 - 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]
12 - Green, J., and B. J. M. Bohannan. 2006. Spatial scaling of microbial biodiveristy. Trends Ecol. Evol. 21:501-507.[CrossRef][Medline]
13 - 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]
14 - 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]
15 - 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]
16 - 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]
17 - 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]
18 - 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]
19 - 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]
20 - 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]
21 - 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]
22 - 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]
23 - 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]
24 - 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]
25 - 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]
26 - Posada, D., and K. A. Crandall. 1998. Modeltest: testing the model of DNA substitution. Bioinformatics 14:817-818.[Abstract/Free Full Text]
27 - 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]
28 - 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]
29 - 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.
30 - 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.