Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2007, p. 7664-7679, Vol. 73, No. 23
0099-2240/07/$08.00+0 doi:10.1128/AEM.01272-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Max Planck Institute for Marine Microbiology, Celsiusstrasse 1, D-28359 Bremen, Germany
Received 8 June 2007/ Accepted 26 September 2007
|
|
|---|
|
|
|---|
Indeed, PCR assays have been developed for the amplification of dsrAB genes, but up to now these approaches are restricted to diversity analyses of the SRP community (2, 3, 15, 42, 43, 57). Although multiple events of lateral gene transfer (LGT) have affected the Dsr phylogeny of certain SRP lineages, the usefulness of dsrAB as functional gene markers for environmental analysis was confirmed (32, 62, 67). Recently, new PCR assays have been developed for the amplification of the aprBA gene locus from SRP and sulfur-oxidizing bacteria (SOB). The ubiquitous presence of aprBA genes in SRP was confirmed (38); however, the PCR-based screening among SOB reference strains revealed its restricted distribution to photo- and chemotrophs with strict anaerobic or at least facultative anaerobic lifestyles, e.g., several Chlorobiaceae and most Chromatiaceae species, Thiobacillus, Thiothrix, as well as invertebrate symbionts and their free-living relatives (37). The functional role of a reverse-acting APS reductase in apr-possessing, aerobic members of the alphaproteobacterial SAR11 and SAR116 cluster (19) is unresolved. Phylogenetic studies revealed the occurrence of LGT events with the results showing some anoxygenic photo- and chemolithotrophic representatives (e.g., Thiocapsa roseopersicina, Thiobacillus spp.) harboring a second, SRP-related apr gene locus in addition to a LGT-unaffected gene locus. In some SOB species (e.g., Thiocystis violascens), only the LGT-derived gene copy is present due to loss of the authentic gene (37). The previous analyses allowed the establishment of a comprehensive AprBA database covering most of the genus level diversity of SRP and apr-harboring SOB which provides a high-resolution framework for the assignment of environmentally retrieved gene sequences in molecular ecological studies. Comparative sequence analyses confirmed the reductive- and oxidative-type APS reductases to be highly conserved among SRP and SOB. Although affected by several LGT events involving representatives of SRP and SOB lineages, the aprBA genes have in general been transmitted vertically during evolution, which supports their usage as functional gene markers in microbial diversity surveys (37, 38). Astonishingly, the genes have only been employed by two medical studies to assess the SRP biodiversity in human tissue (13, 65).
In this study, we present a new assay for the functional aprA gene analysis using a combination of PCR and double-gradient denaturing gradient gel electrophoresis (DG-DGGE). The species diversity coverage of the new aprA gene-targeting primer set is evaluated in comparison with a previously published one (13) by phylogenetic investigation of SRP and SOB reference strains as well as several environmental samples derived from distinct habitats (sediment, filtered seawater, microbial mats, and invertebrate tissue). The aprA gene analysis allowed the simultaneous determination of the environmental SRP and SOB diversity, in contrast to the published dsrAB-based assays (32, 62). Finally, the assay was applied to investigate the sulfate-reducing and sulfur-oxidizing microbial communities inhabiting deep-sea sediment, seawater, and manganese crust samples of the volcanic arc of the Lesser Antilles, Caribbean Sea (Caribflux project); the work is complemented with 16S rRNA gene analyses and SRP and SOB cultivation studies.
|
|
|---|
Study sites and sampling from the Caribbean Sea.
Sediment, seawater, and manganese crusts samples were collected from the volcanic arc of the Lesser Antilles, Caribbean Sea, during the RV Sonne cruise SO-154 (15 January to 8 February 2001) (Table 1). The deep-sea sediment cores and manganese crusts of the Kahouanne Basin (KB) and Montserrat Ridge (MR) were obtained by a multiple corer and a chain bag dredge, respectively. The cores were immediately subsampled by slicing the extruded sediment at 2-cm intervals. Fifteen liters of surface seawater was sampled from the photic water zone near St. Lucia Bay (SLB) by a rosette water sampler and subsequently filtered with Sterivex (pore size, 0.2 µm; Millipore) to collect the pelagic microbial community. For further molecular analyses, the sample material was immediately frozen and stored at –20°C until DNA extraction.
|
View this table: [in a new window] |
TABLE 1. Location and characteristics of environmental samples from habitats investigated in this study by aprA and 16S rRNA gene fragment analyses
|
DNA isolation.
Genomic DNAs from the SRP and SOB reference strains and the enrichment cultures of the MPN series were obtained by using a DNAeasy kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. Total environmental DNAs from (i) filtered seawater (SLB), (ii) the 0- to 2- and 4- to 6-cm horizons of the sediments (KB and MR), and (iii) the manganese crust (MR) were extracted in duplicate using a modified version of the original protocol of Zhou et al. (64) (five cycles of freezing [liquid nitrogen] and thawing [at 30°C] of either 2 g of sediment/manganese crust or one filter [in 5.4 ml of the extraction buffer] were included before starting the extraction procedure). Environmental DNAs from (i) the surface sediments of an intertidal sand flat from "Janssand" on Spiekeroog Island, German Wadden Sea, and of a coastal (reed-covered) area from Vilm Island, Baltic Sea, (ii) microbial mats collected from the thermophilic subsurface environment of the Toyoha Mine, Japan, (iii) filtered water of a hydrothermal spring from Vail, CO, (iv) sulfido- and methanogenic enrichment cultures of anoxic bottom waters of Mono Lake, CA, (v) microbial mats of cold-seep Beggiatoa fields at the Hydrate Ridge, Cascadia margin, OR, and (vi) tissue material of Echinocardium cordatum (Table 1) were kindly provided by M. Mussmann, J. Kuever, M. Fukui, J. C. M. Scholten, T. Treude, and S. Gomes da Silva, respectively (see references 20, 28, 34, 42, 48, and 58 for descriptions of sample collection and the DNA extraction procedures). The quality of the isolated DNAs was verified by PCR amplification of the 16S rRNA gene using the primer set GM3F/GM4R and PCR protocol as described by Muyzer et al. (41). The extracted DNAs (dissolved in 10 mM Tris-HCl, 1 mM EDTA, pH 7.5) were stored at –20°C until further molecular analysis.
PCR amplification of partial aprA and 16S rRNA genes for DGGE analysis.
A modified version of the primer set AprA-3-FW/APS-RV, published by Deplancke et al. (13), and a new primer set, AprA-1-FW/AprA-5-RV, were used to amplify an approximately 0.4-kb aprA gene fragment for subsequent DGGE analysis (GC clamp added to the reverse primers) (Table 2). PCR assays were performed with a REDTaq DNA polymerase kit from Sigma-Aldrich (St. Louis, MO). Reaction mixtures (50 µl total volume) contained 5 µl 10x REDTaq PCR buffer (with 11 mM MgCl2), 5 µl 10x bovine serum albumin solution (3 mg/ml), 200 µM deoxynucleoside triphosphates mixture, 1 µM each primer, and 100 ng genomic DNA as template (10 and 100 ng environmental DNA were used in the assays). "Hot start" PCR was performed in a thermocycler (Eppendorf, Hamburg, Germany) followed by a "touchdown" PCR protocol as described in Table 2.
|
View this table: [in a new window] |
TABLE 2. aprA gene-targeting primers and thermocycling conditionsa used in the PCR assays
|
The PCR products were visually analyzed by electrophoresis of aliquots (10% of the reaction volume) on a 2% (wt/vol) agarose gels run in 1x Tris-borate-EDTA buffer stained with ethidium bromide (0.5 mg liter–1) to verify correct amplicon size. Prior to further analysis, the amplicons of the expected gene fragment size were purified using either a QIAquick gel extraction kit (Qiagen, Hilden, Germany) or a Perfectprep gel cleanup sample kit (Eppendorf, Hamburg, Germany) following the supplier's recommendations.
DG-DGGE analysis of PCR-amplified aprA and 16S rRNA gene fragments.
DG-DGGE analyses (10) of the mixed-template aprA and 16S rRNA gene PCR products were performed using the D-GENE and D-CODE systems (Bio-Rad, München, Germany). For aprA gene fragment analysis, DG-DGGE gels (1.0-mm thick) were poured with an acrylamide gradient from 6 to 8% acrylamide/bis-acrylamide stock solution, 37.5:1 (vol/vol) (Bio-Rad), superimposed over a colinear denaturant gradient from 30 to 60% of denaturant (100% denaturant corresponds to 7 M urea and 40% formamide [vol/vol], deionized with AG501-X8 mixed bed resin [Bio-Rad]). For 16S rRNA gene analysis, DG-DGGE gels (1.0 mm thick) consisted of a 6 to 8% acrylamide gradient superimposed over a colinear 20 to 70% denaturant gradient. Gradients were formed using a Bio-Rad gradient former model 385. Twenty microliters of PCR product was mixed with 6 µl of dye solution (0.1% bromphenol blue [wt/vol], 70% [vol/vol] glycerol) and applied to the gels. The DG-DGGE electrophoresis runs were performed (i) for 2 h at 150 V and subsequently for 2 h at 200 V (aprA gene fragment) and (ii) for 3.5 h at a constant voltage of 200 V (16S rRNA gene fragment) using 1x Tris-acetate-EDTA as buffer and a constant temperature of 60°C. After gel electrophoresis, the gels were stained for 15 min with ethidium bromide (0.5 mg liter–1) and rinsed for 10 min in Milli Q water. The DNA bands were visualized on a UV transillumination table (Biometra, Göttingen, Germany); persisting and dominant bands were excised from the polyacrylamide gel by a scalpel, eluted in 50 µl Tris-HCl, pH 8.0, and reamplified using 1 µl of the eluate as template and PCR conditions as described above. The purity and migration behavior of the reamplification products of the DGGE bands were checked by DG-DGGE. The reamplification products were purified from free PCR primers using either a QIAquick gel extraction kit (Qiagen) or a Perfectprep gel cleanup sample kit (Eppendorf) following the suppliers' recommendations.
Nucleotide sequencing.
The reamplification products of the DG-DGGE bands were sequenced directly in both directions using the respective aprA and 16S rRNA gene amplification primers and an ABI Prism BigDye terminator cycle sequencing ready reaction kit (Applied Biosystems, Foster City, CA) according to the manufacturer's instructions. Sequencing reactions were run on an ABI PRISM 3100 genetic analyzer (Applied Biosystems).
Comparative analysis of the environmentally derived, partial aprA and 16S rRNA sequences.
The partial 16S rRNA sequences obtained from DGGE analysis were checked using the program CHECK_CHIMERA from the Ribosomal Database Project for detection of chimeric artifacts. The sequence data of each DGGE band were then added into the 16S rRNA sequence database of the Technical University Munich, Germany, using the software program package ARB (http://www.arb-home.de), and the data were assembled and first aligned automatically with the ARB_Align tool; the alignment was then manually corrected. The operational taxonomic units (OTUs) were compared to the GenBank database by the BLAST analysis tools of the National Centre for Biotechnology Information (www.ncbi.nlm.nih.gov/BLAST/) in order to determine their phylogenetic affiliation and to identify their closest phylogenetic relatives.
The aprA sequence data obtained from DGGE analysis were assembled and manually corrected using the BioEdit sequence alignment editor (version 7.0.5) (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). The deduced, partial AprA sequences of the environmental samples were implemented into the persisting AprBA alignment of SRP and SOB references strains (see the introduction), including all full-length Apr sequences available from the public databases. The AprA data sets were phylogenetically analyzed with the tree inference methods included in the ARB software package and PHYLIP. Regions of ambiguous homology as well as insertions and deletions (indels) were excluded, yielding a final data set of 110 compared amino acid positions. Unrooted phylogenetic trees were calculated based on 411 sequences by performing distance matrix (neighbor-joining), maximum-parsimony, and maximum-likelihood analyses. The global rearrangement and randomized species input order options as well as the JTT matrix as an amino acid replacement model were used in the phylogenetic analysis.
GenBank accession numbers.
The nucleotide sequence data are available under the GenBank accession numbers EF614307 to EF614329 (archaeal 16S rRNA), EF614394 to EF614449 (bacterial 16S rRNA), and EF618618 to EF618673 and EF633037 to EF633105 (aprA).
|
|
|---|
Species coverage: PCR amplification results of the new aprA gene-targeting primer pair from reference strains.
By using DNA templates from single SRP and SOB reference strains in the PCR assays, no 0.4-kb aprA amplicon was obtained with AprA-3-FW/APS-RV (with a GC clamp) from the thermophilic sulfate reducer lineages, from several members of the Desulfobulbaceae, Desulfobacteraceae, Syntrophobacteraceae, and the Desulfobacterium anilini-related groups, nor from some Desulfotomaculum and Desulfosporosinus spp., Desulfobacca acetoxidans, Caldivirga maquilingensis, and some aprBA gene-containing, sulfur oxidizer lineages of the Gammaproteobacteria. In contrast, the new primer pair allowed the amplification of correct-sized PCR products from SRP and SOB of a wider phylogenetic range (amplification results are summarized in Table S1 and Fig. S1 in the supplemental material), reflecting its improved species coverage compared to the previous set. Nevertheless, we were still unable to yield amplicons from the Desulfosporosinus spp. tested, which might have been caused by inhibited primer annealing (caused by differential codon usage and/or unusual amino acid substitutions in the target sequence). Indeed, the comparison of the primer sequences to the novel apr database revealed the presence of single mismatches of AprA-1-FW and AprA-5-RV with apr sequences of, e.g., Desulfotomaculum spp. (subcluster Ia), Desulfobacteraceae, and some representatives of sulfur-oxidizing Alphaproteobacteria and Gammaproteobacteria. However, the presence of single mismatches at internal or 5'-end sequence positions in the new PCR primers did not affect their PCR efficiency, as it has been demonstrated for 16S rRNA gene-targeting primers (35, 50). In contrast, the presence of several mismatches at the 3' ends of AprA-3-FW and APS-RV inhibited an efficient annealing and elongation in the PCR assays and thus caused the limited species coverage of this primer set.
Preferential amplification: PCR amplification results of the new aprA gene-targeting primer from defined DNA template ratio mixtures by DG-DGGE analysis.
The aprA gene-based approach may suffer from the same methodical constraints as reported for the 16S rRNA gene amplification of mixed-template PCR assays (46, 61), which bias the coverage of phylotypes relative to the original microbial communities they were derived from. To evaluate the template-to-PCR product bias introduced by preferential amplification of certain target sequences, assays were performed using defined template mixtures of genomic DNA from four different SRP reference strains (Desulfosarcina variabilis, Desulfobacter sp. strain DSM 2035, Desulfovibrio profundus, and Desulfotomaculum thermobenzoicum). The DNA concentration of each SRP species was varied to yield relative template mixing ratios of 1:1:1:1, 1:1:1:0.1, or 1:1:1:0.01 in the PCR assays. The denaturing gradient analysis of the obtained mixed-template amplicons indicated that although the AprA-1-FW/AprA-5-RV(GC) primer set appeared to preferentially amplify D. thermobenzoicum (see band intensities on the DGGE gels, Fig. 1), the aprA gene fragments of the remaining SRP species were amplified and detectable in the gels even when the DNA template ratio of one species is only 1:10 relative to the other ones (Fig. 1A). In contrast, the primer pair AprA-3-FW/APS-RV(GC) preferentially amplified D. profundus, leaving only the amplicon of D. thermobenzoicum detectable when the DNA template ratio declined to 1:10 (relative to D. profundus) in the defined mixtures (Fig. 1B). The aprA gene fragment of D. variabilis was only amplified when its template ratio was above 100:1 relative to D. profundus. Thus, the bias introduced by preferential amplification of certain aprA phylotypes appears to be less pronounced with the new primer pair of this study, which is demonstrated by the improved congruence between PCR-detected diversity coverage and actual genetic diversity of SRP in the defined multitemplate mixtures as an example for environmentally derived DNA.
![]() View larger version (64K): [in a new window] |
FIG. 1. Preferential amplification of aprA gene fragments from defined genomic DNA ratio mixtures of the four SRP species Desulfosarcina variabilis (DSM 2060), Desulfobacter sp. (DSM 2035), Desulfovibrio profundus (DSM 11384), and Desulfotomaculum thermobenzoicum (6193) by the primer sets used for DGGE analysis. (A) Primer set AprA-1-FW/AprA-5-RV (GC clamp). (B) Primer set AprA-3-FW/APS-RV (GC clamp). In the PCR assays, the DNA concentration of each SRB species was varied to mixing ratios of 1:1, 1:0.1, or 1:0.01.
|
![]() ![]() View larger version (76K): [in a new window] |
FIG. 2. Phylogenetic tree based on 411 SRP/SOB reference strain and environmentally derived partial AprA sequences. The tree was inferred using the maximum-likelihood method. Sequences of the SOB Apr lineage I and Pyrobaculum aerophilum were used as outgroup references. The 16S rRNA gene-based taxonomic classification of SOB and SRP reference strains is indicated. The environmental AprA sequences are highlighted in boldface type and colored according to the habitat they were retrieved from (Spiekeroog Island, orange [Spkg]; Toyoha Mine, red [MatA/B/D]; Hydrate Ridge, green [HR28Begg]; Echinodermata cordatum, blue [Ec2, -4, -9]; Mono Lake, yellow [HS2, -3, -7, -8]; Vilm Island, brown; Vail, violet). The environmental AprA sequences received with primer set AprA-3-FW/APS-RV are termed with a "D" while those obtained with the new primer set AprA-1-FW/AprA-5-RV are termed with a "N" in the sequence name. For comparison of the microbial diversity coverage of both aprA gene-targeting primer pairs, identical phylotypes detected with both primer pairs are marked with open boxes, whereas those phylotypes that were obtained with only one primer pair are highlighted by dark gray (AprA-3-FW/APS-RV) and light gray (AprA-1-FW/AprA-5-RV) shaded boxes. The scale bar corresponds to 10% estimated sequence divergence.
|
16S rRNA and aprA gene-based analysis: microbial communities of the sulfur cycle from sediment, manganese crust, and seawater samples of the Caribbean Sea.
The microbial communities of the sulfur cycle at four study sites in the area of the Lesser Antilles volcanic arc (Table 1) was explored by molecular and cultivation approaches. Because of its improved coverage of environmental species diversity, only the new aprA gene-targeting primer pair was used to investigate the sediment, manganese crust, and filtered seawater samples (Fig. 3) in comparison with the 16S rRNA gene-based DG-DGGE analysis. The 16S rRNA analysis (summarized in Table 3) revealed a homogenous microbial community structure within the oxic sediment surface layers of the KB and MR (temperature, 4.5 to 4.8°C), comprising members of the Crenarchaeota, Gammaproteobacteria and Deltaproteobacteria, Bacteroidetes, low G+C gram-positive organisms, and Acidobacteria. The OTUs were affiliated with environmental sequences derived from various deep-sea sediments (1, 4, 23, 27, 45, 55) but only distantly with validated SOB and SRP. According to the MPN counting technique, the number of SRP increased in the surface sediments of KB and MR with depth from 4.3 x 105 to 2.4 x 106 cells/g, whereas the number of aerobic SOB remained constant between 4.3 x 106 to 7.5 x 106 cells/g. A single thiosulfate-oxidizing isolate revealed distantly related free-living relatives of thiotrophic mussel symbionts (1), whereas most strains belonged to the commonly culturable alphaproteobacterial lineages, e.g., Sphingomonadaceae and Roseobacter (14, 18). Sulfate-reducing relatives of Desulfovibrio acrylicus were isolated from both study sites. Similar to the 16S rRNA analysis, the aprA gene-based results (Fig. 4) demonstrated a homogenous microbial diversity in the sediments of KB and MR; however, all identified environmental AprA phylotypes belonged to novel alphaproteobacterial SOB clusters (AprA lineage I) that lack reference sequences of cultured SOB species. AprA sequences related to further SOB lineages or SRP were not obtained from the sediment samples. Since aprA genes were absent from the thiosulfate-oxidizing isolates, the aprA gene-detected Alphaproteobacteria are not identical to the cultivated ones but represent yet-unrecognized alphaproteobacterial SOB. Consistent with the 16S rRNA analysis, the AprA sequences derived from the sulfidogenic cultures are closely related to the sequence of Desulfovibrio acrylicus.
![]() View larger version (94K): [in a new window] |
FIG. 3. DGGE analysis of aprA gene fragments (amplified with primer pair AprA-1-FW/AprA-5-RV) using DNA samples from sediment (sed.), manganese crust, and seawater of the Caribbean Sea as templates. Lanes 1 to 8, DNA samples from Kahouanne Basin sediment (horizon 0 to 2 cm, lanes 1 to 4; horizon 4 to 6 cm, lanes 5 to 8); lanes 9 to 12, DNA samples from Montserrat Ridge sediment (horizon 0 to 2 cm); lanes 13 to 16, DNA samples from Montserrat Ridge manganese crust; lanes 17 to 19, DNA samples from aerobic sulfur oxidizer enrichment (enrichm.) cultures inoculated with manganese crust; lanes 20 and 21, DNA samples from St. Lucia Bay filtered (filt.) seawater. Bands indicated with numbers were excised from the gels and sequences.
|
|
View this table: [in a new window] |
TABLE 3. List of phylogenetic affiliations of 16S rRNA partial sequences obtained from DGGE bands of environmental samples and enrichment cultures from the Caribbean Sea (collected in January 2001)a
|
![]() View larger version (44K): [in a new window] |
FIG. 4. Phylogenetic tree based on 411 SRP/SOB reference strain and environmentally derived partial AprA sequences. The tree was inferred using the maximum-likelihood method. Sequences of the Apr lineage I and Pyrobaculum aerophilum were used as outgroup references. The environmental AprA sequences of sediment, seawater, and manganese crust samples of the Caribbean Sea are highlighted in boldface type (for abbreviations used for study sites and enrichment cultures, see Table 1 and Table 3). The 16S rRNA gene-based taxonomic classification of SOB and SRP reference strains is indicated. The scale bar corresponds to 10% estimated sequence divergence.
|
The 16S rRNA-based microbial diversity from seawater of SLB (collected from the thermocline zone; temperature, 22°C) was characteristic for bacterioplankton of photic zones in temperate oceanic waters (12, 18, 30) with OTUs related to marine group II euryarchaeota (12, 17) as well as Prochlorococcus and Synechococcus spp. (66), dominating the archaeal and bacterial community; no phylotypes affiliated with validated SOB and SRP were detected. The MPN counting technique, however, pointed to the presence of 2.4 x 109 aerobic SOB cells/ml and 3.9 x 105 SRP cells/ml in the oxic SLB water column. Thiosulfate-oxidizing members of the commonly cultured Alphaproteobacteria, e.g., Erythrobacter, Roseobacter lineage, and Rhodobacter spp. (7, 14, 18), were predominantly cultivated, whereas no gammaproteobacterial SOB were isolated. Besides the Desulfovibrio-, Desulfobulbus- and Desulfocapsa-related SRB strains isolated, three marine green sulfur bacteria (related to Chlorobium vibrioforme f. sp. thios NCBI 8346, Chlorobium bathyomarinum, and Prosthecochloris sp. strain Vk) could be enriched in coculture with a Desulfovibrio species. Consistent with the 16S rRNA, the aprA gene analysis demonstrated the absence of SAR11 or SAR116 cluster-like AprA phylotypes in the SLB water column. The three environmentally derived AprA phylotypes form a separate planktonic sequence cluster within the Alphaproteobacteria (SOB-Apr lineage I) and a novel Apr lineage of yet-uncertain phylogenetic affiliation (distant relationship to the gamma/betaproteobacterial SOB reference strains). Since aprA genes were absent from the cultivated Alphaproteobacteria, the environmental AprA sequences originated from a yet-undiscovered group of planktonic SOB. No aprA gene fragment could be obtained from the C. vibrioforme f. sp. thios relative, whereas the AprA sequences of the C. bathyomarinum and Prosthecochloris sp. cultures represented a novel Apr lineage of marine Chlorobiaceae. No AprA sequences related to SRP were identified from the oxygenated water column, although the results of the MPN cultures pointed to their presence. However, the phylogenetic placements of these cultured sulfate-reducing strains by AprA and 16S rRNA analysis were consistent.
|
|
|---|
aprA as a functional marker gene for SRP and SOB diversity surveys in comparison with 16S rRNA and dsrAB gene-based studies.
A direct comparison of the diversity analysis results for the aprA and 16S rRNA (and dsrAB) gene approaches is possible in the cases of the subsurface microbial mats of Toyoha Mine (42) and digestive tissue of E. cordatum specimens (20). The aprA gene-based analysis of the microbial mats provided a more complete picture of the dissimilatory microbial sulfur cycle: the Thermodesulforhabdus, Thermodesulfobacterium, and Desulfotomaculum, as well as the (sulfur-oxidizing) Thiobacillus plumbophilus, and free-living invertebrate symbiont relatives remained unidentified by the 16S rRNA gene approaches when applying universal primers (42). According to aprA gene analysis, the symbiotic SRP community of the nodules of E. cordatum is less diverse as it has been suggested by the 16S rRNA data (20). The aprA gene-undetected, putative sulfate-reducing strains were affiliated with sequence clusters of uncultured Deltaproteobacteria, which lack a closely related validated SRP species (only moderate relationships to Desulfobacterium catecholicum, Desulfobacula toluolica, and Desulfotomaculum acetoxidans). This reflects the above-mentioned limitations of the rRNA gene analysis in unambiguously linking the genetic identity of an unknown microorganism to its physiology.
Based on the dsrAB gene analysis of the Toyoha Mine microbial mats, the existence of two novel gram-positive SRP lineages was proposed by Nakagawa and coworkers (42). Indeed, newly discovered Dsr lineages of environmental sequences that lack sequences of cultured reference strains have frequently been found in various ecological studies and used to propose the presence of yet-unrecognized sulfate-respiring representatives (2, 3, 15, 42, 43, 57). However, many sulfite-reducing strains that are incapable of sulfate respiration are closely related to validated SRP, e.g., members of Desulfitobacterium (affiliated with Desulfosporosinus) (53), Sporotomaculum and Pelotomaculum (affiliated with Desulfotomaculum subclusters Ib and Ih) (26), Bilophila and Lawsonia (affiliated with Desulfovibrio), as well as the gram-positive genera Moorella, Carboxydothermus (24, 26), and archaeal Pyrobaculum (25). Since the dsrAB genes of these sulfite reducers have target sites that are fully complementary to the published primer sets (26, 67), the usage of the latter does not allow a differentiation between physiological groups of sulfite and sulfate reducers in PCR amplification. Therefore, the proposal of novel SRP lineages based on Dsr sequences that belong to environmental clusters is ambiguous. In contrast, an identified Apr lineage that lacks a close relationship to SRP reference sequences is a straightforward indication for the presence of an unrecognized phylogenetic lineage of SRP. Thus, our failure to detect additional Desulfotomaculum-related AprA phylotypes in MatB and MatD of Toyoha Mine, besides those affiliated to the thermophilic Desulfotomaculum (the dsrAB-consistent result), might be explained by the presence of gram-positive sulfite respirers.
All published dsr gene-targeting PCR assays are currently restricted to the investigation of SRP diversity because of the low sequence similarity between Dsr in SRP and SOB (32, 62); no universal or SOB-specific amplification primers and no comprehensive DsrAB database of sulfur-oxidizing reference strains exist. The aprA gene-targeting primers of this study, however, allowed the identification of the divergent SOB lineages and even the concomitant investigation of the SRP and SOB species diversity, as demonstrated for some environmental samples e.g., in hydrothermal water (Vail, CO) and microbial mats of Beggiatoa fields (Hydrate Ridge, OR) (Fig. 2). In the latter, the existence of a basal-branching Desulfobulbaceae representative was demonstrated within a microbial community that was otherwise predominated by different chemotrophic SOB species. (Note: a Desulfobulbus relative was identified in these sediment samples by 16S rRNA cloning analysis [T. Losekann, unpublished data]). Nevertheless, the aprA gene-based analysis will underestimate the microbial diversity involved in the oxidative processes of inorganic sulfur compounds in nature as a result of the restricted distribution of the apr genes among the SOB species (see the introduction). As a drawback, this molecular approach will allow their selected detection in environmental samples. However, the low bootstrap values in the SOB lineages show that the partial AprA sequences have only limited phylogenetic information which renders the direct assignment of the oxidative-type environmental AprA sequences to distinct SOB genera/species often difficult. This is also complicated by the LGT-affected evolutionary history of these genes in SOB with the result of some SOB species harboring a second, SRP-related apr gene locus besides an authentic one or only the LGT-derived gene copy (37). Concerning the phylogenetic complexity of SOB in the Beggiatoa field microbial mats, the identified AprA phylotypes pointed to the presence of (at least) one chemolithotrophic sulfur-oxidizing gammaproteobacterium (free-living relative of invertebrate symbionts) and one chemolithotrophic SOB forming a novel lineage of currently uncertain affiliation.
Comparison of 16S rRNA and aprA gene-based analyses: the sulfur-cycling microbial community in the area of the Lesser Antilles volcanic arc (Caribbean Sea).
The results of the 16S rRNA analyses revealed that the microbial community compositions from the aerated deep-sea sediment surfaces of KB and MR and the SLB photic water zone are similar to those of other nonhydrothermal, subtropical study sites in the Mediterranean and East China Seas as well as in the Western Pacific (12, 18, 23, 27, 45, 55). Most identified sediment-associated bacteria might have been involved in the decomposition of detritus-derived organic matter, e.g., the Bacteroidetes comprising chemoorganoheterotrophic degraders of various biopolymers (31), while the archaeal community appeared to be dominated by potential chemolithoautotrophic, nitrifying crenarchaeota (33, 63). In contrast, the SLB bacterioplankton consisted predominantly of phototrophic microorganisms, e.g., (putative) proteorhodopsin-containing euryarchaeota (17) and diverse unicellular cyanobacteria (66). However, no environmental 16S rRNA sequence was closely related to a validated SRP or SOB species; thus, the analyses provided no insights into the microbial community of the sulfur cycle at the investigated study sites. Their presence could only be demonstrated by selective cultivation techniques, with the highest taxa diversity apparent in the SLB photic water zone. The collected geochemical data collected during sampling indicated the absence of a direct, external input of reduced inorganic sulfur compounds at these study sites (Caribflux project final report). Interestingly, the aprA gene analysis revealed the bacterioplankton assemblages and the deep-sea sediments of the Caribbean Sea to be dominated by chemotrophic sulfur oxidizers that belong to four distinct and habitat-specific phylotype clusters within the alphaproteobacterial lineage (closest relatives are the SAR116 group members) and at least one separately branching lineage of yet-uncertain taxonomic affiliation. No phylotypes corresponding to SRP and beta/gammaproteobacterial chemolithoautotrophic SOB were detected at the study sites. The wide distribution of apr genes among members of the Alphaproteobacteria was unexpected because previous studies (37) demonstrated their absence in reference strains of this lineage, with the exceptions of Pelagibacter ubique of the SAR11 group (19) and an uncultivated representative of the SAR116 group. The alphaproteobacterial SAR11 cluster and Roseobacter clade members are the numerically dominant bacterioplankton groups in the decomposition of abundant organic sulfur compounds, e.g., dimethylsulfoniopropionate (DMSP) (21, 22, 36), and thus are implicated in sulfur cycling. Several representatives of the Roseobacter clade have also been reported to harbor the ability to oxidize inorganic sulfur intermediates, e.g., sulfite or thiosulfate to facilitate sulfur-based lithoheterotrophy (6, 21, 39), since reduced inorganic sulfur compounds are produced during the degradation of DMSP (21) or organosulfonates (8, 9). Indeed, thiosulfate-oxidizing, sulfate-producing members of this marine heterotrophic SOB group (e.g., Sulfitobacter spp.) have frequently been identified and isolated from various nonhydrothermal (deep-sea) sediments (51, 56, 59). Based on the results of the aprA gene analysis, we postulate that these chemolithoheterotrophic, sulfur-oxidizing members of the Roseobacter clade are also present in the investigated deep-sea sediment and photic water zone samples from the Caribbean Sea. Indeed, thiotrophic Pelagibaca bermudensis relatives were isolated from both habitats. However, no aprA partial sequence was obtained from the latter; this indicates that the culturable members do not represent the abundant species in these habitats. The functional role of the enzyme "dissimilatory APS reductase" in the alphaproteobacterial strains is unclear: a reverse-acting APS reductase might primarily be used for intracellular detoxification of sulfite as an accumulating product of the organic sulfur compound degradation; the oxidation product APS will be incorporated into the assimilatory sulfur metabolism or converted to sulfate (by the activity of a dissimilatory ATP sulfurylase), allowing additional energy conservation as previously proposed for species of SAR11 and SAR116 (37). The dominant role of chemolithoheterotrophic alphaproteobacterial instead of chemolithoautotrophic gammaproteobacterial SOB in the microbial sulfur transformation processes might be realistic for the investigated nonhydrothermal sediment and water column samples that lack a direct input of inorganic sulfur compounds. Nevertheless, the recognized limitations of PCR amplification (46, 61) and the DGGE method (41) might probably have biased the analysis towards the abundant groups in the microbial community, the SOB, with the consequence of the nondetection of minor present sulfate reducers. The 16S rRNA and aprA gene analyses of the sulfidogenic enrichment cultures point consistently to the presence of Desulfovibrio acrylicus-related strains. Interestingly, this SRB species is capable of DMSP cleavage for subsequent reduction of the acrylate moiety to propionate (60).
In contrast to the previous habitats, the geochemical data pointed to a (at least) subrecent, external (hydrothermal) input of inorganic sulfur compounds at the MR manganese crust collection site, reflected in the sulfur-based, resident microbial community structure. Indeed, OTUs affiliated with validated chemolithoautotrophic SOB species (52) were detected directly from the environmental samples by 16S rRNA gene analysis, while close relatives of SOB derived from hydrothermal systems in the Aegean Sea and Mid-Atlantic Ridge (34, 44) could even be isolated. Nevertheless, sulfur-oxidizing members of phylogenetic lineages recognized to be generally dominant at active deep-sea hydrothermal vents, e.g., Epsilonproteobacteria (44, 54), were not identified. The aprA gene analysis revealed the high phylogenetic complexity of SOB in the crust surface which was not discovered by 16S rRNA analysis. The sulfur-cycling microbial community is suggested to consist of (i) putative chemolithoautotrophic SOB strains of the Betaproteobacteria and Gammaproteobacteria, (ii) sulfur-oxidizing representatives of an uncertain-affiliated, novel SOB lineage, and (iii) members of the previously mentioned chemolithoheterotrophic alphaproteobacterial SOB. The betaproteobacterium appeared to be a marine representative of the genus Thiobacillus; however, a clear assignment of the other phylotypes to certain SOB genera/species is not possible with the amount of sequence information of the partial AprAs and the current database because Apr reference sequences of Beggiatoa and chemolithoautotrophic Epsilonproteobacteria (54) are missing. Previous work on Thiomicrospira spp. indicated the general absence of apr genes in this genus (see also Table S1 in the supplemental material), whereas the APS pathway seemed to be restricted to a few chemolithoautotrophic, marine members of Beggiatoa.
Published ahead of print on 5 October 2007. ![]()
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
-subclass of the class Proteobacteria. Appl. Environ. Microbiol. 65:3810-3819.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»