Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2000, p. 3454-3463, Vol. 66, No. 8
Oklahoma State University, Stillwater,
Oklahoma 740781; Florida State
University, Tallahassee, Florida 323082; and
Pacific Northwest Laboratory, Richland, Washington
993523
Received 27 January 2000/Accepted 30 May 2000
This study was undertaken in an effort to understand how the
population structure of bacteria within terrestrial
deep-subsurface environments correlates with the physical and chemical
structure of their environment. Phylogenetic analysis was performed
on strains of Arthrobacter that were collected from various
depths, which included a number of different sedimentary units from the
Yakima Barricade borehole at the U.S. Department of Energy's
Hanford site, Washington, in August 1992. At the same time that
bacteria were isolated, detailed information on the physical, chemical, and microbiological characteristics of the sediments was
collected. Phylogenetic trees were prepared from the 39 deep-subsurface
Arthrobacter isolates (as well as 17 related type strains)
based on 16S rRNA and recA gene sequences. Analyses based
on each gene independently were in general agreement. These analyses
showed that, for all but one of the strata (sedimentary layers
characterized by their own unifying lithologic composition), the
deep-subsurface isolates from the same stratum are largely
monophyletic. Notably, the layers for which this is true were composed
of impermeable sediments. This suggests that the populations
within each of these strata have remained isolated under constant,
uniform conditions, which have selected for a particular dominant
genotype in each stratum. Conversely, the few strains isolated from a
gravel-rich layer appeared along several lineages. This suggests that
the higher-permeability gravel decreases the degree of isolation of
this population (through greater groundwater flow), creating
fluctuations in environmental conditions or allowing migration, such
that a dominant population has not been established. No correlation was
seen between the relationship of the strains and any particular
chemical or physical characteristics of the sediments. Thus,
this work suggests that within sedimentary deep-subsurface
environments, permeability of the deposits plays a major role in
determining the genetic structure of resident bacterial populations.
One of the goals of microbial
ecology is to determine how environmental factors have shaped the
genetic structure of natural microbial populations, creating patterns
of genetic variability and affecting evolutionary change. The
physical and chemical conditions of the environment act directly on
microbes and therefore play an important role in determining
the genetic structure of their populations (although horizontal
gene transfer also plays a role [20, 31, 47]). Noting
the importance of environmental selection on bacteria, Maynard Smith
(30) included in his model of bacterial population genetics
the concept of ecotypic structure (different clones adapting to their
surrounding environment).
A handful of studies examining free-living bacterial populations have
explored the extent to which the environment has shaped the population
structure. McArthur et al. (32) have shown that the amount
of genetic diversity found in Burkholderia (formerly Pseudomonas) cepacia populations across a
landscape gradient, as measured by multilocus enzyme
electrophoresis, was directly proportional to levels of
local environmental diversity. Conversely, a later study examined
(by multilocus enzyme electrophoresis) the diversity within B. cepacia populations collected along a stream continuum and found
that the population diversity could not be correlated to habitat
diversity, nor was there a correlation between genetic and geographic
patterns (50). The authors suggested that this was the case
because fast mixing may predominate within the stream environment, in
contrast to the more isolated soil environment of the previous study.
Comparable analyses have been performed on various marine bacterial
populations. A study of marine Vibrio species, isolated from
across the water column at locations throughout the world, found that
most of the large taxonomic groups identified inhabited particular
depth ranges, and the substrates available at these depths could be
correlated with the substrates best utilized by the resident group
(41). Similarly, studies on marine cyanobacterial (Prochlorococcus) populations found a correlation between
depths that the populations inhabited and the light adaptation
characteristics of the organisms, dividing the population into
high-light- and low-light-adapted clades (11, 49). However,
a study on populations of marine cyanobacteria of the genus
Synechococcus found no such correlation (48).
Also, Field et al. (12) found that analysis of 16S rRNA
genes from members of the SAR11 cluster (an uncultured group of the
bacterioplankton), cloned from across the water column in locations in
both the Atlantic and Pacific Oceans, showed evidence of niche
partitioning, with certain lineages showing highly depth-specific distributions.
Studies performed over the past decade have firmly established that
microorganisms exist within deep-subsurface environments (100- to
1,000-m depth) (1-3, 5, 14, 15, 17, 23, 24, 33, 39, 42).
The types of microbes present and their relative abundances can be
quite variable (from below detection to 107 cells per g
[dry weight] of sediment) (18). Because the
deep-subsurface environments themselves represent a wide variety of
geologic, hydrologic, and geochemical conditions, it is difficult to
generalize the factors that will be important to the microbial ecology
of all subsurface environments. However, it is recognized that factors such as sediment porosity and texture, nutrient availability, oxygen
conditions, and degree of water saturation and rate of water flow will
be important in determining the types of microbes present and their
distribution, persistence, and activities at depth (18).
We wished to examine the genetic diversity within populations of
terrestrial deep-subsurface bacteria, both to determine the population
structure and to see if this structure could be correlated with the
physical and chemical structure of the environment. We chose to study
the bacterial population diversity in the deep subsurface of the U.S.
Department of Energy (DOE) Hanford site, Washington State, because the
distinct, varying lithofacies present over a narrow depth interval
represent a variety of geochemically and geophysically distinct
habitats. This site had been targeted by the DOE Deep-Subsurface
Science Program for detailed geomicrobiological investigation. Sediment
core samples were collected during drilling of the Yakima Barricade
borehole (YBB) at the Hanford site (well no. 699-48-96) in August 1992. We chose, as our focus group for this study, isolates from the YBB
identified as belonging to the genus Arthrobacter.
Arthrobacter species are common in soils and are aerobic,
high-G+C, chemoorganotrophic, gram-positive bacteria characterized by a
rod-to-coccus morphology change as they enter stationary phase.
Arthrobacter isolates were found to be present in fairly
large numbers in this, as well as other (4, 39, 46),
deep-subsurface environments.
We used phylogenetic sequence analysis of two genes, the 16S rRNA and
recA genes, to help determine the relationships among the
Arthrobacter species. The utility of 16S rRNA as a
phylogenetic molecular tool is well recognized (13, 36, 37,
51). RecA is a relatively conserved protein from bacteria and is
involved in DNA repair and homologous recombination. RecA protein and
gene sequences have been used in a number of studies to determine
bacterial phylogenies (9, 21, 28, 34, 52). These analyses
show phylogenies that are largely congruent with those based on 16S rRNA gene sequences (9, 21, 28). Phylogenies based on the two gene sequences were used to confirm one another in our study. Also,
since recA, as a protein-encoding sequence, is likely to show more variability at the nucleotide level than 16S rRNA gene sequences (since more fluctuation is allowed within protein-coding genes because of codon degeneracy), higher resolution of the more closely related species in trees based on recA might be
expected. Thus, it was hoped that the use of recA might help
to further distinguish differences within the group under study.
Therefore, using phylogenies based on these two molecules, we sought to
explore the genetic diversity of Arthrobacter species
collected from various depths at the Hanford site in order to determine
how the varying environmental conditions created by the different
sedimentary facies may have played a role in shaping the population structure.
Description of the deep-subsurface site.
The Hanford site
lies in the semiarid south-central portion of Washington State. The YBB
is located in the western portion of the Hanford site in an area of low
meteoric water recharge. The depth to the aquifer in this area is
100 m; groundwater recharge occurs laterally from uplands to the
west, making this site hydrologically upstream of any contamination
associated with Hanford operations. The subsurface samples used in this
study were collected at the YBB in August 1992 as part of a project to
explore in detail the microbiology, geohydrology, and geochemistry of
the subsurface (6, 15, 16, 24, 33, 35). The sediment
samples, taken as a set of cores at between 172.9 and 223.0 m, were
well within the saturated zone. This interval is composed of a series
of mostly fine-grained fluvial (river)-lacustrine (lake)-paleosol
(ancient soil) deposits, approximately 6 million to 8 million years
old, stratigraphically within the Ringold Formation, which overlies the
Columbia River Basalt Group (7). A stratigraphic column summarizing the sedimentary lithofacies is shown in Fig.
1, which also indicates the areas from
which core samples were taken and the lithologic descriptions of the
sediments. The top of the interval sampled is composed of fine-grained
lacustrine sediments (173- to 185-m depth) that are underlain by an
approximately 5-cm layer of volcanic tuff (tephra). Between the tephra
and the top of the basalt are represented two complete fluvial
sequences, beginning with gravel at the base and grading up to sand and
finally a paleosol. Thus, there is a paleosol layer (upper paleosol,
185 to 193 m), followed by an interval of fluvial sands (193 to
197 m), a coarse-grained fluvial sandy gravel (197 to 210.8 m),
and another paleosol sequence (lower paleosol, 210.8 to 222.5 m) that
is underlain by another fluvial sand and sandy gravel sequence. The
lowermost fluvial sand and gravel units were not evaluated in this
study.
0099-2240/00/$04.00+0
Copyright © 2000, American Society for Microbiology. All rights reserved.
Genetic Diversity among Arthrobacter
Species Collected across a Heterogeneous Series of Terrestrial
Deep-Subsurface Sediments as Determined on the Basis of 16S
rRNA and recA Gene Sequences


![]()
ABSTRACT
Top
Abstract
Introduction
Materials and Methods
Results and Discussion
References
![]()
INTRODUCTION
Top
Abstract
Introduction
Materials and Methods
Results and Discussion
References
![]()
MATERIALS AND METHODS
Top
Abstract
Introduction
Materials and Methods
Results and Discussion
References

View larger version (27K):
[in a new window]
FIG. 1.
Lithographic and stratigraphic column of the YBB. The
areas from which core samples were taken are indicated, along with
their respective sediment sample identifiers.
Drilling, tracer introduction, and soil sample collection. Coring was performed using cable-tool drilling (a drilling method that relies on percussion to advance the borehole without the use of drilling fluids, which can contaminate core samples). Samples were collected using a 2.5-ft-long, 5-in.-outside-diameter split spoon lined with clear, sterilized lexan liners fitted inside the split spoon. Between core runs, 8-in. steel casing was advanced to the next core interval and cleaned out with a core barrel.
Fluorescent microspheres and potassium bromide were added during coring as particle and solute tracers, respectively, to help detect any contamination that may have infiltrated the inner core that was used for analysis (24). Examination for the presence of the tracers in the inner core material showed that no significant contamination existed in the material used for study (24). Once cores were recovered, they were immediately capped and processed on site in an argon-filled glove bag using flame-sterilized instruments. Processing included cutting away the core liner and paring away approximately 1 cm of the outer surface of the core to remove potentially contaminated material. The topmost approximately 12 to 50 cm of inner core material was then homogenized to provide a common sample with which microbiological and geochemical assays were performed, while the bottom approximately 10 to 12 cm was used for permeability and other physical-property measurements. Physical, chemical, and microbiological tests were done by various laboratories on portions of core samples following overnight shipment on ice of samples packed in argon-filled jars.Physical and geochemical analyses of sediments. Detailed physical and chemical analyses were performed on sediment samples (16, 24, 33). Lithologic descriptions and estimation of grain size distribution were made on site. Hydraulic conductivity measurements were performed by Core Petrophysics, Inc. (Houston, Tex.). The pH, Eh, and dissolved ions (measured using an ion chromatograph) were analyzed on expressed porewaters (extracted by centrifugation). Total organic carbon was measured by Huffman Laboratories (Golden, Colo.) as the difference between total carbon (measured by combustion) and carbonate carbon (measured by acidification).
Bacterial strains.
For isolation of subsurface strains,
serial dilutions were made of homogenized sediment samples (stored at
4°C) in phosphate buffer and plated on a variety of media (2,
4). All colonies arising on plates were collected and preserved
in the DOE Subsurface Microbial Culture Collection (SMCC) (D. Balkwill,
Florida State University). The isolates, as part of the SMCC, were
subjected to a battery of tests, including morphological and
physiological characterizations as well as 16S rRNA similarity rank
analysis (29) and sequence analysis (4). Isolates
that were chosen for these analyses (Table
1) were any that were identified as being
most closely related to Arthrobacter globiformis from
sediment samples that were processed immediately or stored no longer
than 3 to 8 weeks. Because the samples were not all equally screened for isolates, the relative distribution of Arthrobacter
species used in this study should not be interpreted as the relative
abundance of Arthrobacter species along the sampled
interval. Deep-subsurface Arthrobacter isolates were grown
at 25 to 30°C in nutrient broth (Difco Laboratories, Detroit, Mich.).
|
Primer design and PCR amplification, cloning, and sequencing of 16S rRNA and recA genes. DNA was isolated from the Arthrobacter and Micrococcus strains by lysozyme digestion followed by extraction in phenol and chloroform in the presence of hexadecyltrimethylammonium bromide (38) or as described previously (4). PCR amplification and sequencing of 16S rRNA genes from subsurface Arthrobacter species were done as described by Balkwill et al. (4). 16S rRNA sequences for the type strains were obtained from GenBank. The GenBank accession numbers for these sequences are as follows: A. agilis, X80748; A. aurescens, X83405; A. citreus, X80737; A. globiformis, X80736; A. histidinolovorans, X83406; A. nicotianae, X80739; A. nicotinovorans, X80743; A. oxydans, X83408; A. pascens, X80740; A. polychromogenes, X80741; A. protophormiae, X80745; A. sulfureus, X83409; A. ureafaciens, X80744; A. uratoxydans, X83410; M. luteus, M38242; and M. lylae, X80750. These 16S rRNA sequences were determined for DSM type strains, except M. luteus 16S rRNA, for which the ATCC strain was used; we chose ATCC strains for recA sequence determination that matched these DSM strains (listed above). Also used in these analyses were sequences for the Streptomyces ambofaciens ATCC 23877 16S rRNA gene (GenBank accession number M27245) and the S. ambofaciens DSM 40697 recA gene (accession number Z30324).
Primers used in the PCR amplification and sequencing of the Arthrobacter recA genes are shown in Table 2. Since, prior to this study, no recA sequences had been determined for any Arthrobacter species, for initial analysis of recA, degenerate primers GPRA-FB and GPRA-R2, designed against highly conserved areas of the RecA protein, were used in the PCR to amplify an interior ~350-bp fragment from total DNA of A. globiformis. This fragment was gel purified (Qiagen QIAquick Gel Purification Kit), cloned into the PCR product-cloning vector pCR2 (Invitrogen TA Cloning Kit), and sequenced. In order to amplify a larger region of recA for phylogenetic studies, GPRA-UF2 and GPRA-UR2 were designed against other conserved areas based on an alignment of the amino acid sequences of seven previously characterized recA genes (three from other high-G+C gram-positive genera) that were most related to the translation of the A. globiformis 350-bp fragment. Using these primers, successful amplification of an ~830-bp band was obtained from a number of Arthrobacter species. These products were cloned as described above from A. globiformis and one of the deep-subsurface Arthrobacter species. Sequencing of these clones facilitated design of more specific primers, A19-F2 and A1-R, that were used to amplify and sequence many of the recA genes. The growing collection of sequences was used to design alternate primers (AU-F1 and AU-R1) that were used to amplify and sequence other Arthrobacter recA genes. AU-FM1 and AU-RM1 were used to complete the sequence in both directions. PCR products were gel purified and sequenced (200 fmol per reaction) on an Applied Biosystems model 373A automated DNA sequencer.
|
Phylogenetic data analysis. The 16S rRNA gene sequences were aligned, and after the ends of the sequences were trimmed to equal lengths, 1,160 bases (including gaps) of the sequences were included in the phylogenetic analyses (corresponding to nucleotide positions 197 to 1325 of the A. globiformis 16S rRNA gene). The recA sequences were aligned and trimmed to equal lengths, and 360 bases were included in the phylogenetic analyses (corresponding to nucleotide positions 312 to 669 of the Escherichia coli recA gene). Alignment of 16S rRNA sequences was performed by eye, using the secondary structure of A. globiformis 16S rRNA (Ribosomal Database Project) (29) as a guide. The recA nucleotide sequences were aligned using the program Clustal V (19), aided by the high conservation within the translations of the sequence. Analysis of sequence composition was done using the MEGA program (27). Phylogenetic analyses were done using the PHYLIP 3.572 package of programs (10), using both neighbor-joining/distance matrix (DNADIST and NEIGHBOR) and parsimony (DNAPARS) methods.
Nucleotide sequence accession numbers. The GenBank accession numbers for the 16S rRNA sequences generated for this report are as follows (SMCC strain, accession number): G915, AF197020; G919, AF197021; G950, AF197022; G954, AF197023; G959, AF197024; G960, AF197025; G961, AF197026; G962, AF197027; G963, AF197028; G964, AF197029; G965, AF197030; G966, AF197031; G968, AF197032; G969, AF197033; G970, AF197034; G979, AF197035; G980, AF197036; G982, AF197037; G984, AF197038; G986, AF197039; G991, AF197040; G993, AF197041; ZAT001, AF197042; ZAT004, AF197044; ZAT005, AF197045; ZAT012, AF197046; ZAT013, AF197047; ZAT014, AF197048; ZAT031, AF197049; ZAT054, AF197050; ZAT055, AF197051; ZAT056, AF197052; ZAT200, AF197053; ZAT255, AF197054; ZAT262, AF197055; ZAT263, AF197056; ZAT277, AF197057; ZAT351, AF196342; and ZAT352, AF196343. The GenBank accession numbers for the recA sequences generated for this report are as follows (SMCC or type strain, accession number): G915, AF214757; G919, AF214758; G950, AF214759; G954, AF214760; G959, AF214761; G960, AF214744; G961, AF214762; G962, AF214763; G963, AF214764; G964, AF214776; G965, AF214777; G966, AF214765; G968, AF214766; G969, AF214739; G970, AF214740; G979, AF214767; G980, AF214742; G982, AF214768; G984, AF214741; G986, AF214769; G991, AF214743; G993, AF214746; ZAT001, AF214745; ZAT004, AF214748; ZAT005, AF214747; ZAT012, AF214755; ZAT013, AF214770; ZAT014, AF214773; ZAT031, AF214756; ZAT054, AF214771; ZAT055, AF214772; ZAT056, AF214774; ZAT200, AF214775; ZAT255, AF214750; ZAT262, AF214754; ZAT263, AF214751; ZAT277, AF214753; ZAT351, AF214749; ZAT352, AF214752; A. agilis, AF214779; A. aurescens, AF214793; A. citreus, AF214781; A. globiformis, AF214780; A. histidinolovorans, AF214788; A. nicotianae, AF214792; A. nicovorans, AF214784; A. oxydans, AF214789; A. pascens, AF214786; A. polychromogenes, AF214785; A. protophormiae, AF214790; A. sulfureus, AF214787; A. ureafaciens, AF214782; A. uratoxydans, AF214791; M. luteus, AF214783; and M. lylae, AF214778.
| |
RESULTS AND DISCUSSION |
|---|
|
|
|---|
Sediment chemical, physical, and microbiological characteristics. In this study, we explored the phylogenetic diversity of deep-subsurface Arthrobacter populations within several sedimentary lithofacies of the YBB at the DOE Hanford site. The strains were isolated from sediment samples taken at depths of 172.9 to 217.7 m (Fig. 1) during the August 1992 coring of the YBB. With the same sediment samples, detailed microbiological, physical, and chemical properties were analyzed. Table 1 details physical and chemical information for the sediment samples and indicates from which samples Arthrobacter species used in these analyses were isolated.
The sampled interval consists of varying lithologies, each with distinct physical and chemical characteristics. In general, total organic carbon (TOC) was highest in the lacustrine layer, except for the uppermost portion, which includes sample G1, at 172.9-173.8 m. The top 2 m of the lacustrine unit is an oxidized layer, characterized by iron and manganese oxides and low organic carbon, associated with highly oxidizing groundwater in the gravels immediately above the lacustrine sediments. TOC was also relatively low in the upper paleosol and in the fluvial sands and gravels. In general, the pH was slightly alkaline throughout the lacustrine interval and was more neutral to acidic throughout the upper paleosol and lower layers (except for fluvial sand sample G31 [196.3 to 197.0 m], at pH 8.8). The Eh indicated, in general, moderately oxidizing conditions throughout the interval, being somewhat less oxidizing in the upper paleosol. Overall, nitrate concentrations varied little, and phosphate concentrations were extremely low. Sulfate concentrations were relatively low throughout the lacustrine unit but began increasing in the upper paleosol, reaching a maximum within sample G21 (188.7 to 189.6 m), after which the concentration dropped off dramatically and remained low throughout the rest of the interval. Ammonium concentrations were very low in all samples, mostly below detection (data not shown). The salinity (major dissolved cations) (data not shown) was relatively constant throughout the sampled interval, with the exception of three samples (G2 [173.8 to 174.3 m], G17 [185.2 to 186.0 m], and G27 [192.9 to 193.6 m]) which exhibited lower cation concentrations than other samples (none of these were samples in which Arthrobacter species used in these studies originated). The lacustrine, fluvial sand (which was well cemented), and paleosol units have particularly low porosity and hydraulic conductivities, which inhibit movement of groundwater, nutrients, and bacteria through these layers. Thus, the lacustrine, upper paleosol, and fluvial sand units, combined, form a low-permeability hydraulic barrier between two highly permeable gravel layers. Extremely little chemical, physical, or microbial data are available for the lower paleosol unit. Several groups have investigated the microbial characteristics of the sampled interval (16, 24, 33). In general, investigators have found microbial activities and numbers (both aerobic and anaerobic) to be very low across the interval. The numbers of cells and direct measurements of their metabolism were highest in the lacustrine layers, where TOC was highest, and were lowest in the fluvial sands, where TOC was also low (24). The persistence of organic matter in the lacustrine layer may be because of the lack of its availability to microbes, probably due to the extremely low permeability of this unit (24). The results of these investigations suggest that patchy aerobic and anaerobic microenvironments exist throughout these sediments, allowing low levels of various anaerobes as well as strict aerobes, including Arthrobacter species, to persist within the sediments.Sequencing of 16S rRNA and recA genes from Arthrobacter strains. The recA and 16S rRNA gene sequences were determined for 39 YBB deep-subsurface Arthrobacter species (Table 1). Seventeen Arthrobacter and related-genus type strains were also included in the analyses for comparison. The overall pairwise nucleotide sequence identity for deep-subsurface Arthrobacter strains for 16S rRNA genes was 95.3 to 100%, and that for recA was 78.1 to 100%. Occasionally, multiple deep-subsurface strains were found to have identical 16S rRNA gene sequences and identical recA gene sequences. However, these strains were not always isolates from the same sediment sample. For example, ZAT004 and ZAT013 (from the upper lacustrine facies, sediment sample G1) were found to have identical sequences for both genes, but their sequence identity is shared by ZAT054 (from the lower paleosol unit, sediment sample G37). Similarly, ZAT055 (from the lower paleosol unit, sediment sample G37) was found to have sequences identical to those from G984 and G993 (both from the fluvial gravel unit, sediment sample G33). Therefore, the fact that isolates exhibiting identical sequences may originate from different sediment samples indicates that the identities cannot simply be the result of isolating essentially identical organisms (originating from the breakup of microcolonies during homogenization of the sediment sample or from siblings that may have multiplied during the short-term storage of some sediment samples). The majority of nucleotide differences in recA sequences between deep-subsurface strains were synonymous substitutions: while there were 124 variable sites in the 360 nucleotide positions analyzed, there were only 17 variable sites in the 120 amino acid residues. Preliminary trees generated using the RecA protein sequences showed that there was not enough sequence variation in the amino acid sequence to be phylogenetically useful, and therefore the nucleic acid sequences of recA were used in the phylogenetic analyses.
Phylogenetic analyses.
Phylogenetic trees based on 16S rRNA
and recA gene sequences were generated by maximum-parsimony
analyses using S. ambofaciens as an outgroup to root the
trees (Fig. 2 and
3). Bootstrap analyses were performed,
and bootstrap values of 50% or greater are shown at the appropriate
nodes. Distance matrix analyses gave trees with topologies very similar
to those obtained by the parsimony method for the two genes. The
recA gene trees give relationships in general agreement with
those generated by 16S rRNA analysis. This suggests that true
phylogenetic relationships are being observed and that the
relationships have not been obscured by any horizontal gene transfer
that may have taken place. (Moreover, recent horizontal gene transfer
within these deep-subsurface sediments is unlikely given the low cell
numbers present and low permeability of sediments.) The differences
that do exist between the 16S rRNA and recA gene trees are
in the deeper (more ancestral) nodes of the trees. These deeper nodes
in the recA gene trees have bootstrap values of less than
50%, which indicates that they are not statistically significant. There is more variability within the recA nucleotide
sequences of the group than there is within the 16S rRNA gene
sequences. (For instance, for our data sets the average percent
sequence identities within clusters and between clusters based on 16S
rRNA sequence analysis are 99.2 to 99.9% and 96.9 to 98.1%,
respectively; the values based on recA sequence analysis are
89.4 to 100% and 85.6 to 89.1%, respectively.) This leads to the
lower bootstrap values at the deeper nodes but higher resolution at the
shallower nodes, as seen by the high bootstrap values on these nodes.
Thus, within recA sequences there has been so much change
between distantly related groups (between clusters) that these deeper
relationships cannot be properly distinguished using recA,
but they can be distinguished using the more conserved 16S rRNA
sequences. However, the less distant relationships (within clusters)
can be better resolved with recA than with 16S rRNA.
|
|
Conclusions. We have sampled the diversity within Arthrobacter populations in a diverse assemblage of deep-subsurface strata. The genetic structure of the populations at the YBB appears to be predominantly created and controlled by the degree of physical isolation and physical-chemical homogeneity of individual sedimentary lithofacies. For the majority of facies, a single population, most likely one that was present at the time of sediment deposition or one that migrated into the layer after deposition and was most suited to survival in that particular sediment type, has persisted and become the dominant culturable Arthrobacter type. It will be interesting to see if similar relationships between the structures of the environment and resident microbial populations hold true for other genera within the YBB or for other deep-subsurface sites and how this population structure may be altered by human activities in the subsurface. An understanding of the structure of deep-subsurface bacterial populations and how this is related to the structure of the environment would help direct deep-subsurface bioremediation efforts, as well as help direct the search within the deep subsurface for microbes that possess unique, not previously documented, characteristics for biotechnology applications.
| |
ACKNOWLEDGMENTS |
|---|
We thank Brendan Bohannan for helpful discussions in the preparation of the manuscript.
This research was supported by the Deep Microbiology Subprogram of the Subsurface Science Program, Office of Energy Research, U.S. Department of Energy, under grants DE-FG02-93ER61680 (to R.V.M.) and DE-FG02-96ER62210 and DE-FG05-91ER61159 (to D.L.B.).
| |
FOOTNOTES |
|---|
* Corresponding author. Mailing address: Microbiology and Molecular Genetics, 307 Life Science East, Oklahoma State University, Stillwater, OK 74078-3020. Phone: (405) 744-6243. Fax: (405) 744-6790. E-mail: rum67{at}okway.okstate.edu.
Present addresses: Carnegie Institution of Washington, Stanford, CA 94305.
Present address: Center for Microbial Ecology, Michigan State
University, East Lansing, MI 48824.
| |
REFERENCES |
|---|
|
|
|---|
| 1. |
Amy, P. S.,
D. L. Haldeman,
D. Ringelberg,
D. H. Hall, and C. Russell.
1992.
Comparison of identification systems for classification of bacteria isolated from water and endolithic habitats within the deep subsurface.
Appl. Environ. Microbiol.
58:3367-3373 |
| 2. | Balkwill, D. L. 1989. Numbers, diversity, and morphological characteristics of aerobic, chemoheterotrophic bacteria in deep subsurface sediments from a site in South Carolina. Geomicrobiol. J. 7:33-52. |
| 3. |
Balkwill, D. L.,
J. K. Fredrickson, and J. M. Thomas.
1989.
Vertical and horizontal variations in the physiological diversity of the aerobic chemoheterotrophic bacterial microflora in deep southeast coastal plain subsurface sediments.
Appl. Environ. Microbiol.
55:1058-1065 |
| 4. | Balkwill, D. L., R. H. Reeves, G. R. Drake, J. Y. Reeves, F. H. Crocker, M. B. King, and D. R. Boone. 1997. Phylogenetic characterization of bacteria in the subsurface microbial culture collection. FEMS Microbiol. Rev. 20:201-216[CrossRef][Medline]. |
| 5. | Boivin-Jahns, V., R. Ruimy, A. Bianchi, S. Daumas, and R. Christen. 1996. Bacterial diversity in a deep-subsurface clay environment. Appl. Environ. Microbiol. 62:3405-3412[Abstract]. |
| 6. | Chandler, D. P., F. J. Brockman, T. J. Bailey, and J. K. Fredrickson. 1998. Phylogenetic diversity of archaea and bacteria in a deep subsurface paleosol. Microb. Ecol. 36:37-50[CrossRef][Medline]. |
| 7. | Department of Energy. 1988. Consultation draft, site characterization plan, reference repository location, Hanford Site, Washington. DOE/RW-0164, vol. 1. U.S. Department of Energy, Washington, D.C. |
| 8. |
Duwat, P.,
S. D. Ehrlich, and A. Gruss.
1992.
A general method for cloning recA genes of gram-positive bacteria by polymerase chain reaction.
J. Bacteriol.
174:5171-5175 |
| 9. | Eisen, J. A. 1995. The RecA protein as a model molecule for molecular systematic studies of bacteria: comparison of trees of RecAs and 16S rRNAs from the same species. J. Mol. Evol. 41:1105-1123[Medline]. |
| 10. | Felsenstein, J. 1989. PHYLIP-phylogeny interference package (version 3.2). Cladistics 5:164-166. |
| 11. | Ferris, M. J., and B. Palenik. 1998. Niche adaptation in ocean cyanobacteria. Nature 396:226-240[CrossRef]. |
| 12. | Field, K. G., D. Gordon, T. Wright, M. Rappé, E. Urbach, K. Vergin, and S. J. Giovannoni. 1997. Diversity and depth-specific distribution of SAR11 cluster rRNA genes from marine planktonic bacteria. Appl. Environ. Microbiol. 63:63-70[Abstract]. |
| 13. |
Fox, G. E.,
E. Stackebrandt,
R. B. Hespell,
J. Gibson,
J. Maniloff,
T. A. Dyer,
R. S. Wolfe,
W. E. Balch,
R. S. Tanner,
L. J. Mgrum,
L. B. Zablen,
R. Blakemore,
R. Gupta,
L. Bonen,
B. J. Lewis,
D. A. Sahl,
K. H. Leuhrsen,
K. N. Chen, and C. R. Woese.
1980.
The phylogeny of the prokaryotes.
Science
209:457-463 |
| 14. | Fredrickson, J. K., T. R. Garland, R. J. Hicks, J. M. Thomas, S. W. Li, and K. M. McFadden. 1989. Lithotrophic and heterotrophic bacteria in deep subsurface sediments and their relation to sediment properties. Geomicrobiol. J. 7:53-66. |
| 15. | Fredrickson, J. K., S. W. Li, F. J. Brockman, D. L. Haldeman, P. S. Amy, and D. L. Balkwill. 1995. Time-dependent changes in viable numbers and activities of aerobic heterotrophic bacteria in subsurface samples. J. Microbiol. Methods 21:253-265[CrossRef]. |
| 16. | Fredrickson, J. K., J. P. McKinley, S. A. Nierzwicki-Bauer, D. C. White, D. B. Ringelberg, S. A. Rawson, S.-M. Li, F. J. Brockman, and B. N. Bjornstad. 1995. Microbial community structure and biogeochemistry of Miocene subsurface sediments: implications for long-term microbial survival. Mol. Ecol. 4:619-626. |
| 17. | Ghiorse, W. C. 1989. Special issue on deep subsurface microbiology. Geomicrobiol. J. 7:1-136. |
| 18. | Ghiorse, W. C., and J. T. Wilson. 1988. Microbial ecology of the terrestrial subsurface. Adv. Appl. Microbiol. 33:107-172[Medline]. |
| 19. |
Higgins, D. G.,
A. J. Bleasby, and R. Fuchs.
1991.
CLUSTAL V: improved software for multiple sequence alignment.
CABIOS
8:189-191 |
| 20. |
Jain, R.,
M. C. Rivera, and J. A. Lake.
1999.
Horizontal gene transfer among genomes: the complexity hypothesis.
Proc. Natl. Acad. Sci. USA
96:3801-3806 |
| 21. |
Karlin, S.,
G. M. Weinstock, and V. Brendel.
1995.
Bacterial classifications derived from RecA protein sequence comparisons.
J. Bacteriol.
177:6881-6893 |
| 22. | Keddie, R. M., M. D. Collins, and D. Jones. 1986. Genus Arthrobacter Conn and Dimmick 1947, 300AL, p. 1228-1301. In P. H. A. Sneath, N. S. Mair, M. E. Sharpe, and J. G. Holt (ed.), Bergey's manual of systematic bacteriology, vol. 2. Williams and Wilkins, Baltimore, Md. |
| 23. | Kieft, T. L., P. S. Amy, F. J. Brockman, J. K. Fredrickson, B. N. Bjornstad, and L. L. Rosacker. 1993. Microbial abundance and activities in relation to water potential in the vadose zones of arid and semiarid sites. Microb. Ecol. 26:59-78. |
| 24. | Kieft, T. L., J. K. Fredrickson, J. P. McKinley, B. N. Bjornstad, S. A. Rawson, T. J. Phelps, F. J. Brockman, and S. M. Pfiffner. 1995. Microbiological comparisons within and across contiguous lacustrine, paleosol, and fluvial subsurface sediments. Appl. Environ. Microbiol. 61:749-757[Abstract]. |
| 25. | Koch, C., F. A. Rainey, and E. Stackebrandt. 1994. 16S rDNA studies on members of Arthrobacter and Micrococcus: an aid for their future taxonomic restructuring. FEMS Microbiol. Lett. 123:167-172[CrossRef]. |
| 26. |
Koch, C.,
P. Schumann, and E. Stackebrandt.
1995.
Reclassification of Micrococcus agilis (Ali-Chohen 1889) to genus Arthrobacter agilis comb. nov. and emendation of the genus Arthrobacter.
Int. J. Syst. Bacteriol.
45:837-839 |
| 27. | Kumar, S., K. Tamura, and M. Nei. 1993. MEGA: molecular evolutionary genetics analysis, version 1.01. The Pennsylvania State University, University Park. |
| 28. | Lloyd, A. T., and P. M. Sharp. 1993. Evolution of the recA gene and the molecular phylogeny of bacteria. J. Mol. Evol. 37:399-407[Medline]. |
| 29. |
Maidak, B. L.,
J. R. Cole,
J. C. T. Parker,
G. M. Garrity,
N. Larsen,
B. Li,
T. G. Lilburn,
M. J. McCaughey,
G. J. Olsen,
R. Overbeek,
S. Pramanik,
T. M. Schmidt,
J. M. Tiedje, and C. R. Woese.
1999.
A new version of the RDP (Ribosomal Database Project).
Nucleic Acids Res.
27:171-173 |
| 30. | Maynard Smith, J. 1991. The population genetics of bacteria. Proc. R. Soc. Lond. B Biol. Sci. 245:37-41[CrossRef]. |
| 31. |
Maynard Smith, J.,
N. H. Smith,
M. O'Rourke, and B. G. Spratt.
1993.
How clonal are bacteria?
Proc. Natl. Acad. Sci. USA
90:4384-4388 |
| 32. |
McArthur, J. V.,
D. A. Kovacic, and M. H. Smith.
1988.
Genetic diversity in natural populations of a soil bacterium across a landscape gradient.
Proc. Natl. Acad. Sci. USA
85:9621-9624 |
| 33. | McKinley, J. P., T. O. Stevens, J. K. Fredrickson, J. M. Zachara, F. S. Colwell, K. B. Wagnon, S. C. Smith, S. A. Rawson, and B. N. Bjornstad. 1997. Biogeochemistry of anaerobic lacustrine and paleosol sediments within an aerobic unconfined aquifer. Geomicrobiol. J. 14:23-39. |
| 34. | Nowak, A., and J. Kur. 1995. Genomic species typing of acinetobacters by polymerase chain reaction amplification of the recA gene. FEMS Microbiol. Lett. 130:327-332[CrossRef][Medline]. |
| 35. | Ogram, A., W. Sun, F. J. Brockman, and J. K. Fredrickson. 1995. Isolation and characterization of RNA from low-biomass deep-subsurface sediments. Appl. Environ. Microbiol. 61:763-768[Abstract]. |
| 36. | Olsen, G. J. 1988. Phylogenetic analysis using ribosomal RNA. Methods Enzymol. 164:793-812[Medline]. |
| 37. |
Olsen, G. J.,
C. R. Woese, and R. Overbeek.
1994.
The winds of (evolutionary) change: breathing new life into microbiology.
J. Bacteriol.
176:1-6 |
| 38. | Peloquin, L., and C. W. Greer. 1993. Cloning and expression of the polychlorinated biphenyl-degradation gene cluster from Arthrobacter M5 and comparison to analogous genes from gram negative bacteria. Gene 125:35-40[CrossRef][Medline]. |
| 39. | Phelps, T. 1995. Special issue: deep subsurface microbiology. J. Microbiol. Methods 21:225-328. |
| 40. | Retallack, G. J. 1990. Soils of the past. Unwin Hyman, Boston, Mass. |
| 41. |
Simidu, U., and K. Tsukamoto.
1985.
Habitat segregation and biochemical activities of marine members of the family Vibrionaceae.
Appl. Environ. Microbiol.
50:781-790 |
| 42. | Sinclair, J. L., and W. C. Ghiorse. 1989. Distribution of aerobic bacteria, protozoa, algae, and fungi in deep subsurface sediments. Geomicrobiol. J. 7:15-31. |
| 43. |
Stackebrandt, E.,
C. Koch,
O. Gvozdiak, and P. Schumann.
1995.
Taxonomic dissection of the genus Micrococcus: Kocuria gen. nov., Nesterenkonia gen. nov., Kytococcus gen. nov., Dermacoccus gen. nov., and Micrococcus Cohn 1872 gen. emend.
Int. J. Syst. Bacteriol.
45:682-692 |
| 44. | Stanley, J. T. 1999. Bacterial biodiversity: a time for place. ASM News 65:681-687. |
| 45. | Stanley, J. T., and J. J. Gosink. 1999. Poles apart: biodiversity and biogeography of sea ice bacteria. Annu. Rev. Microbiol. 53:189-215[CrossRef][Medline]. |
| 46. | Stim, K. P. 1995. A phylogenetic analysis of micro-organisms isolated from subsurface environments. Mol. Ecol. 4:1-10[Medline]. |
| 47. | Syvanen, M. 1994. Horizontal gene transfer: evidence and possible consequences. Annu. Rev. Genet. 28:237-261[Medline]. |
| 48. | Toledo, G., and B. Palenik. 1997. Synechococcus diversity in the California Current as seen by RNA polymerase (rpoC1) gene sequences of isolated strains. Appl. Environ. Microbiol. 63:4298-4303[Abstract]. |
| 49. |
West, N. J., and D. J. Scanlan.
1999.
Niche-partitioning of Prochlorococcus populations in a stratified water column in the eastern North Atlantic Ocean.
Appl. Environ. Microbiol.
65:2585-2591 |
| 50. | Wise, M. G., L. J. Shimkets, and J. V. McArthur. 1995. Genetic structure of a lotic population of Burkholderia (Pseudomonas) cepacia. Appl. Environ. Microbiol. 61:1791-1798[Abstract]. |
| 51. |
Woese, C. R.
1987.
Bacterial evolution.
Microbiol. Rev.
51:221-271 |
| 52. | Zhou, J., and B. G. Spratt. 1992. Sequence diversity within the argF, fbp and recA genes of natural isolates of Neisseria meningitidis: interspecies recombination within the argF gene. Mol. Microbiol. 6:2135-2146[CrossRef][Medline]. |
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||