Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2008, p. 7265-7271, Vol. 74, No. 23
0099-2240/08/$08.00+0 doi:10.1128/AEM.00540-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

University of Colorado—Boulder, Department of Geological Sciences, UCB 399, Boulder, Colorado 80309,1 University of Colorado—Boulder, Department of Ecology and Evolutionary Biology, Campus Box 334, Boulder, Colorado 803092
Received 5 March 2008/ Accepted 30 September 2008
|
|
|---|
|
|
|---|
The microbial oxidation of manganese (Mn) influences a number of biogeochemical processes, since Mn(III/IV) oxides are chemically reactive and can (i) sequester numerous trace metals through sorption processes, (ii) serve as terminal electron acceptors for microbial respiration, and (iii) directly oxidize complex organic matter (7, 8, 30, 33, 35). Thus, to unravel several linked chemical cycles, it is important to understand the evolution of microbial Mn(II) oxidation capabilities and the microorganisms and genes involved.
In Bacillus species, mnxG, a putative multicopper oxidase, acts to catalyze the oxidation first of Mn(II) to Mn(III) and then of Mn(III) to Mn(IV) (6, 39). mnxG is localized in the exosporium, a layer that occurs outside of the spore coat in some Bacillus species (11, 12, 38). Mn-oxidizing Bacillus strains have been obtained from coastal surface sediments near San Diego, CA; in sediments and particulate plumes at deep-sea hydrothermal vents in the Guaymas Basin in the Gulf of California (6, 13, 24); and now at the Henderson Molybdenum Mine in Empire, CO.
The Henderson Mine was proposed as a site for a Deep Underground Science and Engineering Laboratory (DUSEL) in 2003 and has only recently been investigated for its potential to harbor unique microbial life in the deep subsurface (29). Anoxic, CO2- and metal-rich brines are released into aerated tunnels through boreholes at depths of >1 mile below the surface. Millimolar concentrations of Mn were detected in the fluids, and Mn oxide minerals were precipitating around the mouths of the boreholes. Eight Bacillus sp. isolates capable of Mn(II) oxidation were isolated from the mineral precipitates, while the mnxG gene was amplified from only one of these Mn-oxidizing isolates.
The inability to amplify mnxG from the other isolates led us to question if these isolates were more phylogenetically distinct from the mnxG-positive isolate than they were from each other and, if so, if 16S rRNA sequences could be used to predict the presence or absence of an amplifiable mnxG gene within the genomes of the Henderson Mine isolates. We hypothesized that the isolates would be dispersed throughout the Bacillus clade and predicted that the mnxG-positive isolate, HM06-02, would cluster with other mnxG-positive Bacillus isolates. We also predicted that their placement within the 16S rRNA phylogeny would relate to the presence or absence of mnxG within the genome of the isolates. We discovered a strong relationship between 16S rRNA sequence similarity and the presence or absence of an amplifiable mnxG gene in the Henderson Mine isolates. Additionally, the data revealed a complex phylogenetic distribution of the mnxG gene that we can use as a tool to aid in unraveling the history and dispersal of the mnxG gene within the Bacillus clade.
|
|
|---|
3,000 ft below the surface, where a series of boreholes were drilled into granite in March 2005 to release deep-seated fluids. At the time of sampling, the draining fluids were
45°C and the ambient air temperature was
24°C. Secondary minerals, including Fe oxides, Mn oxides, and Ca sulfates, precipitated around the mouth of each borehole. Sterile Falcon tubes were filled with mineral precipitates and stored on ice (
8 h) until return to the laboratory.
Isolation and identification of Mn(II)-oxidizing bacteria.
Samples were diluted 100-fold with water from the site that had been filtered through a sterile 0.2-µm polycarbonate filter. Samples were plated onto solid medium containing organic carbon and amended with 1 mM soluble Mn(II) and incubated in the dark at 25°C or 50°C. Growth occurred on Leptothrix medium, a chemically defined freshwater medium for the isolation of Mn-oxidizing bacteria (25), and K medium, a seawater-based rich medium containing yeast extract and peptone as carbon sources (34) (Table 1). Single colonies were transferred until they were morphologically uniform (
6 generations), and then isolates were grown in liquid medium for DNA extraction (DNeasy blood and tissue kit; Qiagen).
|
View this table: [in a new window] |
TABLE 1. Organisms isolated from the Henderson molybdenum mine
|
Determination of Mn(II)-oxidizing capabilities and amplification of mnxG sequences.
To evaluate their Mn(II) oxidation capability, isolates were grown on plates amended with soluble Mn(II). After 1 week, the experimentally determined amount of time needed for sporulation to occur, dark colonies, potentially encrusted with brown Mn(IV) oxides, were tested for the presence of oxidized Mn by using a colorimetric leucoberbelin blue (LBB) assay (19). All Bacillus sp. isolates were initially tested for the presence of the mnxG gene by amplification with degenerate mnxG primers previously designed by Dick et al. (6), but only unrelated products were generated. More-specific primers, mnxGIF, 5'-ACGCATGTCTTTCACTATCATGTTCAT-3', and mnxGIR, 5'-AAATAAGTGGTCATGGAAGAACCATGC-3', previously designed by Francis et al. (13), yielded more-appropriate products. The PCR protocol was the following: 30 cycles of 94°C for 30 s, 45°C for 30 s, and 60°C for 1 min, followed by 1 cycle of 72°C for 15 min (13). PCR products were either directly sequenced by SeqWright DNA Technology Services, Houston, TX, or were cloned into chemically competent Escherichia coli cells by using a TOPO-TA cloning kit (Invitrogen). Clones containing the appropriate
800-bp PCR product were purified and sequenced by using the primer M13F. A BLAST search (1) of each sequence was conducted to determine if the sequences represented an mnxG gene.
Phylogenetic analysis. (i) Sequence attainment.
Of the 11 organisms isolated, eight were Bacillus, two were Paenibacillus sp., and one was a Georgenia sp. A total of 57 16S rRNA sequences were used in this study and include 24 Bacillus species not known to possess an mnxG gene, 23 sequences from mnxG-positive Bacillus organisms, the 7 sequences of mnxG-negative Henderson Mine Bacillus isolates, 1 sequence from the mnxG-positive Henderson Mine isolate, and the 2 sequences of the Henderson Mine isolates most similar to Paenibacillus species. 16S rRNA sequences were obtained by conducting nucleotide-nucleotide BLAST searches of Bacillus sp. strain SG-1 (1). Bacillus sp. strain SG-1 is a model organism used for the study of biotic manganese oxidation and possesses the mnxG gene (12, 24, 33, 35). The possession of an mnxG gene was determined by conducting a protein BLAST search using the protein sequence of the mnxG gene of Bacillus sp. strain SG-1. All Bacillus organisms identified in the mnxG search were included in the analysis to represent Mn-oxidizing members of the genus Bacillus. The protein product and nucleotide sequences of the mnxG gene in these Bacillus species were acquired from GenBank. BLAST searches were conducted on the 16S rRNA sequences of the isolates, and sequences of known Bacillus organisms with high identity to the Henderson Mine isolates (HM06-01 to -10) were included in the analysis. Many of the sequences were also used in previous Mn oxidation studies (6, 13).
(ii) Sequence alignment.
16S rRNA sequences were aligned using ClustalX (37) with default settings as follows: gap opening, 10.00; gap extension, 0.20; delay divergent sequences, 30%; DNA transition weight, 0.50; and negative matrix, OFF. An alignment of the mnxG protein sequence was constructed in ClustalX by using the following parameters: gap opening, 10.00; gap extension, 0.20; delay divergent sequences, 30%; protein weight matrix, Gonnet series; protein matrix, IUB; residue-specific penalties, ON; hydrophilic penalties, ON; hydrophilic residues, GPSNDQEKR; gap separation distance, 4; and end gap separation, OFF. GB02-29 was not included in the mnxG gene alignment because a protein sequence was not available for this organism. The nucleotide sequence of the mnxG gene of Henderson Mine isolate HM06-02 was translated into a protein sequence by using the translation tool ExPASy (http://www.expasy.ch/tools/dna.html) and included in the alignment.
(iii) Construction of trees.
A Bayesian analysis of the aligned 16S rRNA genes from all organisms was conducted using MrBayes (17). One million generations of an MCMC run were conducted with the following parameters: nst, 6; rates, gamma; sample frequency, 100; and nchains, 4. A total of 10,000 trees were constructed; the first 1,000 trees (equal to the first 10% or 100,000 generations of the analysis) were attributed to the burn-in period and excluded from further analyses. A 50% majority consensus tree was built in PAUP (version 4.0b10) (31) by using the results of the Bayesian analysis.
PAUP (version 4.0b10) (31) was used to construct a maximum parsimony and 50% majority rule bootstrap tree by using the protein sequence of the mnxG gene.
(iv) Character mapping.
Stochastic character mapping implemented with SIMMAP (2) allowed the estimation of the number of gains and losses of mnxG genes from the presence or absence data. The analysis was performed on 1,000 Bayesian 16S rRNA trees. We tested the hypothesis of no association between phylogeny and the presence or absence of mnxG by comparing the number of gains and losses for the observed data with the number of gains and losses after randomizing the presence or absence data across taxa 30 times, yielding a total of 30,000 mappings of gains and losses. The probability that this hypothesis is true was determined by following the method of Jones and Martin (18).
Nucleotide sequence accession numbers.
The partial 16S rRNA sequences of strains HM06-01 to HM06-11 were submitted to GenBank under accession numbers EU004563 to EU004572 and EU525854 (Table 1). The single mnxG gene identified was submitted to GenBank under accession number EU004573.
|
|
|---|
mnxG from a Henderson Mine Bacillus isolate is deeply nested within mnxG gene phylogeny.
We successfully amplified and sequenced a single mnxG gene from isolate HM06-02. The gene, labeled HM06-02 mnxG, is most closely related to sequences recovered from marine sediments from Point Loma, CA, and is deeply nested within a well-supported clade (bootstrap value = 99) of mnxG genes sampled from four widely different localities: Point Loma, Guaymas Basin, Mission Bay, and Henderson Mine (6, 11, 13; this work) (Fig. 1).
![]() View larger version (12K): [in a new window] |
FIG. 1. Maximum parsimony bootstrap 50% majority rule consensus tree based on mnxG amino acid sequences. Bootstrap values are indicated at selected nodes and are based on 100 bootstrap replicates. The SG-1 cluster is labeled and is consistent across both the 16S rRNA and mnxG gene trees. The HM06-02 cluster is also labeled.
|
![]() View larger version (21K): [in a new window] |
FIG. 2. Fifty-percent majority rule consensus phylogenetic tree based on 16S rRNA sequences of Bacillus species. Henderson Mine isolates are in boldface type, and the single asterisk indicates the mnxG-positive Bacillus organism isolated in this study. Consensus values are indicated at select nodes and are based on 9,000 Bayesian trees. The SG-1 and the HM06-02 clusters are labeled.
|
![]() View larger version (13K): [in a new window] |
FIG. 3. Histogram depicting the probabilities of numbers of gains and losses of the mnxG gene needed to explain the observed data.
|
![]() View larger version (16K): [in a new window] |
FIG. 4. Histogram comparing the number of changes needed within the trees constructed by Bayesian analysis to explain the distribution of the mnxG gene in the observed data versus the randomized data. Note that the distribution of observed values falls to the left of the distribution of random values; this indicates that the null hypothesis can be rejected.
|
|
|
|---|
In the case of the Henderson Mine isolates, the results of 16S rRNA consensus tree analyses illustrate that the mnxG-positive isolate HM06-02 is most closely related to other mnxG-positive Bacillus organisms, while the mnxG-negative isolates were most closely related to other Bacillus organisms not known to possess an mnxG gene (Fig. 1). The evidence of phylogenetic conservation of the mnxG gene was further supported by the occurrence of an SG-1 cluster on both the 16S rRNA and mnxG trees (Fig. 1 and 2).
The phylogenetic inheritance of the mnxG gene within the Bacillus clade is further supported by stochastic character-mapping analyses. The most-probable number of gains and losses of the mnxG gene for the data set (13) is less than the number of changes (21) expected if the gene is randomly distributed on the phylogenetic trees (Fig. 4). These results indicate that phylogenetic inheritance plays a dominant role in the pattern of distribution of the gene. However, the multiple occurrences of gene loss suggest that there is a dynamic aspect to the evolution of the mnxG gene within the Bacillus clade.
The results of stochastic character-mapping analyses provided evidence for the multifaceted evolutionary process of mnxG within the Bacillus clade. Gene loss may explain relationships between taxa as seen on the 16S rRNA tree (Fig. 2), such as within the HM06-02 cluster or within the cluster of four isolates most closely related to GB02-2A, of which Bacillus novalis LMG21835 is the only isolate that lacks the gene (Fig. 2). The loss of genes is not uncommon and may be enhanced if the adaptation that the gene codes for is not used by the organism (4, 9, 16, 20, 22, 26). In Bacillus species, mnxG is only located in the exosporium of metabolically dormant spores (10). This may, in part, explain the frequency of loss of the gene, as the capability to oxidize Mn may not be necessary to the survival of metabolically dormant spores.
Despite the potential for the simultaneous occurrence of vertical inheritance and gene loss, it is still possible to use 16S rRNA-based phylogeny to predict to a first order the presence or absence of an amplifiable mnxG gene within the genome of a Bacillus organism.
The results of LBB tests for the presence of Mn(III/IV) oxides reveal that seven of eight Henderson Mine Bacillus isolates are capable of oxidizing manganese (Table 1), yet the mnxG gene was only amplified in one isolate, HM06-02. We present two lines of evidence supporting the observation that the mnxG gene is absent from the genomes of isolates HM06-01 and HM06-03 to -10: (i) the mnxG gene was not amplifiable, and (ii) the mnxG-negative organisms fall within other Bacillus clusters that are also not known to possess an mnxG gene. These data are in contrast to data from Guaymas Basin, where the authors were able to amplify an mnxG gene in all of their Mn(II)-oxidizing Bacillus species (6). There are two possible explanations for the lack of an amplifiable mnxG gene in Mn(II)-oxidizing Bacillus isolates from the Henderson Mine: (i) the mnxG-negative Bacillus sp. isolates may possess an entirely different gene that is responsible for the oxidation of manganese, or (ii) the primers designed to amplify mnxG may target only a subset of mnxG genes possessed by Bacillus sp. organisms. Henderson Mine Bacillus species are not the only examples of this paradox. Bacillus pumilus isolates capable of vigorously oxidizing manganese have been isolated from Mn crusts on basaltic rocks (A. S. Templeton and L. E. Mayhew et al., unpublished data), and yet mnxG could not be amplified from the genomic DNA in tests done in our and other labs (e.g., G. J. Dick, personal communication [regarding results of amplification using primers and conditions found in reference 6]). This was further corroborated by 16S rRNA phylogenies where B. pumilus does not fall within clusters of Mn(II)-oxidizing Bacillus species (6). Fortunately, the genome of Bacillus pumilus SAFR-032 has recently been completed (15). Protein sequences of all the putative manganese multicopper oxidases, mnxG (Bacillus sp. strain SG-1), cumA (Pseudomonas putida MT-2), and mofA (Leptothrix discophora SS-1), were searched against the Bacillus pumilus SAFR-032 genome by using BLAST. No genes homologous to mnxG were evident from the BLAST search; thus, it appears that mnxG was not amplified because it is not present in the genome. However, both cumA and mofA were most similar to cotA (24% and 27% identity, respectively), an outer spore coat protein of B. pumilus, supporting the idea that perhaps Bacillus sp. organisms possess multiple genes whose products have the functional capability of Mn(II) oxidation.
Conclusions.
16S rRNA phylogenies were determined to be a tool for predicting the presence or absence of a functional gene. The results of Bayesian analyses and character mapping indicate that 16S rRNA relationships between Bacillus organisms can be used to predict the presence or absence of the mnxG gene within the genome of a Bacillus organism. The one mnxG-positive Henderson Mine isolate falls in a cluster consisting predominantly of other mnxG-positive organisms. The mnxG-negative Henderson Mine isolates fall within clusters built solely of other mnxG-negative Bacillus organisms. These observations are supported by the high consensus values obtained from the Bayesian analyses. These conclusions agree with those found by other researchers investigating the relationship of 16S rRNA phylogenies to functional gene phylogenies (21, 23, 27, 28).
Gene loss was found to influence the occurrence of functional genes within the genomes of a subset of Bacillus species. A greater number of losses than gains have occurred (12 versus 1). In previous research, gene loss has been shown to influence the appearance of microbial genomes (4, 9). This research indicates that to fully understand the distribution of the mnxG gene within the genus Bacillus, the processes of phylogenetic inheritance and genomic remodeling through gene loss need to be considered.
Funding for this work was provided by National Science Foundation grant MCB-0623815 (A. S. Templeton) and graduate fellowships awarded to L. E. Mayhew and E. D. Swanner by the Department of Geological Sciences, University of Colorado at Boulder.
Published ahead of print on 10 October 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»