Original Research ARTICLE
High Diversity of Anaerobic Alkane-Degrading Microbial Communities in Marine Seep Sediments Based on (1-methylalkyl)succinate Synthase Genes
- 1Department of Molecular Ecology, Max Planck Institute for Marine Microbiology, Bremen, Germany
- 2HGF MPG Joint Research Group for Deep-Sea Ecology and Technology, Max Planck Institute for Marine Microbiology, Bremen, Germany
Alkanes comprise a substantial fraction of crude oil and are prevalent at marine seeps. These environments are typically anoxic and host diverse microbial communities that grow on alkanes. The most widely distributed mechanism of anaerobic alkane activation is the addition of alkanes to fumarate by (1-methylalkyl)succinate synthase (Mas). Here we studied the diversity of MasD, the catalytic subunit of the enzyme, in 12 marine sediments sampled at seven seeps. We aimed to identify cosmopolitan species as well as to identify factors structuring the alkane-degrading community. Using next generation sequencing we obtained a total of 420 MasD species-level operational taxonomic units (OTU0.96) at 96% amino acid identity. Diversity analysis shows a high richness and evenness of alkane-degrading bacteria. Sites with similar hydrocarbon composition harbored similar alkane-degrading communities based on MasD genes; the MasD community structure is clearly driven by the hydrocarbon source available at the various seeps. Two of the detected OTU0.96 were cosmopolitan and abundant while 75% were locally restricted, suggesting the presence of few abundant and globally distributed alkane degraders as well as specialized variants that have developed under specific conditions at the diverse seep environments. Of the three MasD clades identified, the most diverse was affiliated with Deltaproteobacteria. A second clade was affiliated with both Deltaproteobacteria and Firmicutes likely indicating lateral gene transfer events. The third clade was only distantly related to known alkane-degrading organisms and comprises new divergent lineages of MasD homologs, which might belong to an overlooked phylum of alkane-degrading bacteria. In addition, masD geneFISH allowed for the in situ identification and quantification of the target guild in alkane-degrading enrichment cultures. Altogether, these findings suggest an unexpectedly high number of yet unknown groups of anaerobic alkane degraders and underline the need for comprehensive surveys of microbial diversity based on metabolic genes in addition to ribosomal genes.
Alkanes are found throughout nature and belong to the most abundant organic compounds in the biogeosphere (Wilkes et al., 2002). They are main components of crude oil and natural gas and are either formed by living organisms (Tissot and Welte, 1984) or through geological transformation of biomass (Claypool and Kvenvolden, 1983). In the marine environment, alkanes occur in gas hydrates, petroleum-rich hydrothermal sediments, seeps or areas anthropogenically contaminated such as occurred during the Exxon Valdez oil spill in 1989 or the Deep Water Horizon blowout in 2010. Marine seeps can differ remarkably in their hydrocarbon composition and concentration. Methane seeps mainly emit methane in micromolar to millimolar ranges m-2 d-1 (e.g., Hydrate Ridge in the NE Pacific, Suess et al., 1999), gas seeps emit a substantial amount of C2-C5 gasses in addition to methane (e.g., Mediterranean Amon mud volcano, Mastalerz et al., 2009) and hydrocarbon seeps emit a broad range of alkanes, alkenes and aromatics (e.g., at Guaymas Basin in the Gulf of California, Byrne and Emery, 1960; Simoneit and Lonsdale, 1982). Many microorganisms are able to utilize these hydrocarbons as their carbon source (Widdel et al., 2010). In seep sediments, degradation mainly takes place under anoxic conditions as oxygen is generally depleted within the first few millimeters of the sediment. A large fraction of sulfate reduction (SR) at gas and hydrocarbon seeps is fueled by the anaerobic oxidation of methane (AOM, Reeburgh, 2007). However, as indicated by a global median ratio of SR to AOM of 10.7, a major part of total SR is fueled by the oxidation of non-methane hydrocarbons, in particular the oxidation of alkanes (Bowles et al., 2011).
Microbial anaerobic oxidation of alkanes has been described for a large range of alkanes: gaseous (C2–C5) alkanes (Kniemeyer et al., 2007; Savage et al., 2010; Jaekel et al., 2012; Adams et al., 2013; Bose et al., 2013), mid-chain (C6–C12) alkanes (Ehrenreich et al., 2000; Davidova et al., 2006) and long-chain (C13–C20) alkanes (Aeckersberg et al., 1991; So and Young, 1999; Zengler et al., 1999; Cravo-Laureau et al., 2004a). Cultivated anaerobic alkane degraders use sulfate (sulfate-reducing bacteria, SRB), nitrate, manganese or ferric iron Fe(III) as electron acceptors (Weelink et al., 2009; Widdel et al., 2010). Furthermore, hydrocarbon-degrading enrichment cultures have been established under methanogenic conditions (Zengler et al., 1999; Chang et al., 2006; Berdugo-Clavijo and Gieg, 2014; Embree et al., 2014). Isolated or enriched anaerobic alkane degraders belong to two phyla: Proteobacteria and Firmicutes (Figure 1). Only recently an archaeon, Archaeoglobus fulgidus (Euryarchaeota), has been shown to degrade long-chain alkanes (C10–C21) with thiosulfate or sulfate (Khelifi et al., 2014). In marine environments alkane degradation is predominantly performed by SRB within the class Deltaproteobacteria and in particular with members of the family Desulfobacteraceae (Figure 1). Members of the Desulfosarcina/Desulfococcus (DSS) clade have been shown to be key players in seep sediments (Knittel et al., 2003; Acosta-González et al., 2013; Kleindienst et al., 2014). Stable isotope probing identified four specialized DSS clades that are active in the oxidation of short- and long-chain alkanes (Kleindienst et al., 2014). Cultivation-independent studies using clone libraries or next generation sequencing techniques (Teske et al., 2002; Orcutt et al., 2010; Ruff et al., 2015) and alkane incubation studies (Savage et al., 2010; Adams et al., 2013; Bose et al., 2013) showed a high diversity of SRB at hydrocarbon seeps suggesting potential for a much larger diversity of alkane degraders as is currently known from available isolates.
FIGURE 1. Phylogenetic tree showing the affiliation of 16S rRNA gene sequences from isolated or enriched anaerobic hydrocarbon degraders to selected reference sequences of the domain Bacteria. Nitrate-reducing bacteria are printed in blue, sulfate-reducing bacteria in red, iron-reducing bacteria in purple, phototrophic bacteria in green, fermentative bacteria in orange and syntrophic bacteria degrading hydrocarbons in a consortium under methanogenic conditions in light blue. Georgfuchsia toluolica, printed in light green, has been shown to use Fe(III), Mn(IV), and nitrate as terminal electron acceptor for growth on aromatic compounds. Cd. Methylomirabilis oxyfera, printed in brown, can oxidize methane anaerobically by utilizing oxygen produced internally from dismutation of nitric oxide into nitrogen and oxygen gas. Substrate usage is given within parenthesis. The bar represents 10% estimated sequence divergence.
Several biochemical reactions have been described for alkane activation under anoxic conditions (Callaghan, 2013; Musat, 2015). The most well-described and particularly dominant pathway is the addition of alkanes to fumarate yielding alkylsuccinates. The enzyme involved in this initial activation step is a glycyl radical enzyme of the pyruvate formate lyase family, the (1-methylalkyl)succinate synthase, Mas (Grundmann et al., 2008), which has also been referred to as alkylsuccinate synthase, Ass (Callaghan et al., 2010). The putative catalytic subunit of Mas is subunit D (MasD), which is equivalent to Ass subunit A (AssA). Alkane activation by MasD has been described for a range of n-alkanes and cycloalkanes in sulfate-reducing and nitrate-reducing isolates and enrichment cultures (Kropp et al., 2000; Callaghan et al., 2006; Kniemeyer et al., 2007; Grundmann et al., 2008; Musat et al., 2010; Jaekel et al., 2015). Furthermore, fumarate addition has also been suggested for methanogenic alkane degradation based on the detection of assA/masD genes in enrichments (Davidova et al., 2011; Mbadinga et al., 2011; Zhou et al., 2012; Aitken et al., 2013; Cheng et al., 2013). Recently, an assABC operon was detected on a Smithella single cell genome (Tan et al., 2014). Considering the wide distribution of alkane activation via fumarate addition, MasD/AssA serves as valid biomarker for anaerobic alkane degradation. For alkane degradation by achaeon Archaeoglobus fulgidus, an alkylsuccinate synthase activity was hypothesized for pyruvate formate lyase (Pfl) D based on a higher similarity of PflD with AssA compared with reference Pfls (Khelifi et al., 2014). Yet, there are only few environmental studies on AssA/MasD diversity in contaminated soils and groundwater (Wang et al., 2012; Zhou et al., 2012; von Netzer et al., 2013), contaminated river sediments and aquifers (Callaghan et al., 2010) or marine hydrocarbon-impacted sediments (Acosta-González et al., 2013; von Netzer et al., 2013; Johnson et al., 2015). Considering the high number of long branches in AssA/MasD phylogenetic trees these studies point to the existence of a broad diversity of microorganisms involved in the degradation of alkanes.
In this study we addressed the diversity of the alkane-degrading microbial community by massive parallel 454-tag sequencing of masD genes retrieved from 12 globally distributed marine seep sediments and correlate their characteristics with environmental parameters like hydrocarbon composition, water depth, temperature and SR rates. We hypothesize that in marine seep sediments the diversity of the anaerobic alkane-degrading microbial community is much higher than previously known, including new deeply branching taxonomic lineages. Due to the narrow range of substrate used by isolated alkane degraders (Musat, 2015) we hypothesize that the MasD-microbial community at methane and gas seeps differ significantly from that at hydrocarbon seeps. As an integral part of the study, we optimized the geneFISH protocol to identify alkane-degrading communities in situ.
Materials and Methods
Sediments were sampled from seven globally distributed marine seeps differing in their geographical, biological, chemical, and geological features (Table 1, Figure 2). Investigated methane seeps predominantly releasing methane were located on the Cascadian Margin at Hydrate Ridge (Pacific Ocean, station HR19), on the Hikurangi continental margin (New Zealand, Wairarapa, station NZ315) and in the central North Sea (Tommeliten, station Tomm). Gas seeps at the Mediterranean Amon mud volcano (AMV, stations AMV760, AMV825) release significant amounts of other gasses (C2–C5 alkanes) in addition to methane. Hydrocarbon seeps, in the northern (stations GoM4463 and GoM156), and southern Gulf of Mexico (Chapopote Asphalt Volcano: stations GoM140, GoM13, GoM17) and in the Guaymas Basin (stations GB4573 and GB4484) in the Gulf of California are characterized by seepage of complex hydrocarbons (Table 1).
FIGURE 2. Marine seep sites investigated in this study. Green: hydrocarbon seeps, orange: methane seeps; purple: gas seeps.
Nucleic Acid Extraction
DNA was extracted by mechanical, chemical and enzymatic cell lysis steps according to the protocol of Zhou et al. (1996) from sediments frozen immediately after sampling at -20°C. Extracted nucleic acids were washed with ice-cold ethanol (80% v/v), centrifuged at 14,000 × g for 10 min at RT, dried for 15–30 min at RT and gently resuspended without pipetting in TE buffer (10 mM Tris-HCl, pH 8.0, 1 mM EDTA, pH 8.0) for 1 h at 4°C.
MasD Amplification and Pyrosequencing
The masD/assA gene was amplified using one of two primer pairs: 7757f-1,f-2 (TCG GAC GCG TGC AAC GMY CTG A; MasD amino acid position 395 in strain HxN1; accession number CAO03074)/8543r (TCG TCR TTG CCC CAY TTN GG; position 657 in HxN1) or primer pair 7766f (TGT AAC GGC ATG ACC ATT GCG CT; position 398 in HxN1)/8543r (TCG TCR TTG CCC CAY TTN GG) (von Netzer et al., 2013). Primers were barcoded and extended with an SfiI restriction site at the 5′ end for ligation with the 454-adapters. For each sample, eight replicate PCRs (20 μl volume) per primer pair were carried out containing, 0.5 μM primer each, 250 μM dNTPs, 0.3 μg μl-1 BSA, 1 × PCR buffer, 0.25 U Taq polymerase (5Prime, Germany) under the following conditions: initial denaturation at 95°C for 5 min, followed by 35 cycles of denaturation (96°C, 1 min), annealing (58°C, 1 min), elongation (72°C, 2 min), and a final elongation step (72°C, 10 min). Replicate PCR reactions of both primer pairs were pooled and the 800 bp-amplicons were then extracted from an agarose gel (1.5% w/v) and purified using the MiniElute PCR Purification Kit (Qiagen) according to the manufacturer’s recommendations. Massive parallel tag sequencing of the amplicons was carried out on a 454 Life Sciences GS FLX sequencer (Roche, Basel, Switzerland) at the Max Planck-Genome-Center, Cologne, Germany.
For pairwise comparison of 16S rRNA genes and MasD (see below) we amplified and Sanger sequenced masD from the following alkane-degrading strains and enrichments (Table 2): PF2803 (DSM16219, Cravo-Laureau et al., 2004b), LM2801 (DSM18843, Cravo-Laureau et al., 2007), and Propane60-GuB (Kniemeyer et al., 2007).
Raw reads were submitted to a rigorous quality control using a mothur version 1.29.1 routine (Schloss et al., 2009) including stringent quality filtering and chimera check using UCHIME (Edgar et al., 2011). Sequences were removed from the data set if they had ≥1 mismatch to the forward primer, were <200 bp in length, had >8 homopolymers, had an average quality score <20 (qthreshold = 20) or contained any ambiguities. Sequences were then translated in silico (RevTrans 1.4 Server) and screened for MasD based on the presence of the amino acid motive FECIR, FECIK, FECQR, FECVR, FDCIR or FDNIA. Ninety percent of MasD/AssA proteins in our database possess one of these motives at position 435 (HxN1). Sequences with stop codons were removed from the dataset. Sequences were further checked for the presence of the catalytically active cysteine (Grundmann et al., 2008) at position 477 of HxN1. We could not confirm the presence of the conserved motif (RVXG) that harbors the radical-storing glycine (Becker et al., 1999) and is characteristic for all glycyl radical enzymes including MasD and BssA as it lies outside the region amplified by our primers (position 811; HxN1).
Establishment of a Protein Database for Glycyl Radical Enzymes
We built a comprehensive protein database containing more than 10,000 sequences for MasD, BssA, NmsA and pyruvate formate lyase (Pfl). Sequences were retrieved from this study or from public databases GenBank, NCBI and DDBJ or from publically available complete or draft genomes, single cell genomes, and metagenomes accessible via the Integrated Microbial Genome (IMG) with microbiome samples (IMG/M) system at JGI1 (Markowitz et al., 2014). The mas operon sequence from strain HxN1 (AM748709) was used for a BlastX search versus the IMG/M database resulting in >1400 potential glycyl radical enzyme sequences. We then imported all amino acid sequences retrieved from these databases along with our own sequences retrieved from seep sediments into the software package ARB (Ludwig et al., 2004) and aligned them using MAFFT v7 (Katoh and Standley, 2013). Manual correction of frame shifts based on BLASTX was necessary for about 15 of our seep sequences due to insertions or deletions caused by 454-pyrosequencing. The amino acid alignment is provided in Supplementary information (Supplementary Table S1).
Definition of OTU Cut-Offs for MasD Amino Acid Sequences
Taxonomic units for MasD were defined based on 15 alkane-degrading pure cultures by pairwise comparison of their amino acid sequences with the corresponding 16S rRNA gene sequences. A list of strains used for this analysis is provided in Table 2. The sequence difference D and similarities S (S = 1 – D) were calculated. Finally, the similarity of MasD amino acid sequence pairs were plotted versus the similarity of the 16S rRNA gene sequence pairs of the same strains.
Phylogenetic Tree Reconstruction
The phylogenetic tree based on small subunit ribosomal RNA genes was calculated with nearly full-length sequences (>1350 bp) available in Arb Silva database release 111 (Pruesse et al., 20072) by neighbor-joining analysis in combination with filters which consider only 50% conserved regions of the 16S rRNA. Partial sequences were subsequently inserted into the reconstructed consensus tree by parsimony criteria, without allowing changes in the overall tree topology. MasD-based phylogenetic tree was constructed by maximum likelihood analysis (PhyML algorithm, Blosum 62 substitution model) considering 95 amino acid positions (position 436 to 506, strain HxN1; CAO03074) using 441 deduced amino acid sequences. Only one representative sequence per MasD family level OTU0.72 is shown in the final tree.
Community Diversity Analysis
Sequence abundance tables were generated by clustering the retrieved MasD seep sequences at 96% amino acid identity based on 120 amino acid positions (Pos. 398-500, HxN1) using a distance matrix in Mothur (Schloss et al., 2009) and used to calculate inverse Simpson diversity indices and species rarefaction. Bray–Curtis dissimilarities (Bray and Curtis, 1957) between all samples were calculated and used for two-dimensional non-metric multidimensional scaling (NMDS) ordinations with 20 random starts (Kruskal, 1964). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. Hierarchical clustering of all samples was performed using Ward’s method (Ward, 1963), which minimizes the total within-cluster variance. A network was built based on a presence absence matrix. The network vertices (nodes) were plotted using a Fruchterman and Reingold (1991) force-directed algorithm, which causes an increase in the nodes attraction to each other with increasing similarity between them; the more OTU0.96 shared between two samples, the closer they are in the network. Analyses were carried out with the R statistical environment and the packages vegan (Oksanen et al., 2012), ggplot2 (Wickham, 20093) and network (Butts, 2008) in addition to custom R scripts. In silico coverage (C) of MasD sequences was calculated per station according to the following equation
where n is the number of singletons (SSOabs + SSOrel.) and N the total number of sequences analyzed.
Design of Probes for masD Detection
To set-up a masD geneFISH assay, an enrichment culture with n-butane was established under sulfate-reducing conditions from hydrocarbon seep sediments of site GB4573. 5 ml sediment slurry was made using a 1:1 mix of sediment from the upper 10 cm with artificial anoxic seawater (Widdel and Bak, 1992) and incubated in hungate tubes containing 5 ml anoxic media for anaerobic sulfate-reducing microorganisms (Widdel and Pfennig, 1982). As substrate, 1 bar n-butane gas (Messer Griesheim GmbH, Krefeld, Germany) was added to the headspace. All tubes were kept under N2/CO2 (90/10, v/v) atmosphere and were horizontally incubated at 28°C. When sulfide production reached 15 mM, 10% of the enrichment culture was subsequently transferred as inoculum into fresh n-butane supplemented media resulting in sediment-free enrichment cultures. DNA was extracted from the enrichment GB4573_14 and masD genes were amplified using primers ass/bss F and ass/bss R as described previously (Callaghan et al., 2010). Cloning and sequencing of 689 bp-amplicons was performed as described previously (Kleindienst et al., 2012). Based on 28 sequenced clones used to represent the breadth of diversity of masD in our sample, we designed a probe mix consisting of 9 masD-targeting dsDNA polynucleotide probes (316 bp; GC content 65%) using the PolyPro software (Moraru et al., 2010, 2011, Supplementary Table S2). The individual probes showed >75% nucleotide sequence identity to the target regions of all members of the masD gene clusters retrieved from the GB4573_14 enrichment and >72% nucleotide sequence identity to all retrieved MasD OTUs0.72. For probe synthesis, plasmid DNA was extracted from 9 selected clones using the SpinMiniprep Kit (Qiagen, Hilden, Germany) and dsDNA probes were synthesized by PCR with enzymatic incorporation of Dig-labeled nucleotides (dUTP) using the PCR Dig Probe Synthesis Kit (Roche, Diagnostics, Mannheim, Germany). As negative control, probe NonPolyPr350 was synthesized. Synthesized dsDNA probes were purified with the GeneClean Turbo kit (Q-Biogene).
Fluorescence In Situ Hybridization of masD (geneFISH)
One ml was subsampled from the enrichment culture GB4573_14 at three different time points, fixed with formaldehyde (2% final concentration in 1x PBS (pH 7.4) for 1 h at RT) and an aliquot was filtered on polycarbonate filters (GTTP, pore size: 0.2 μm). Filters were incubated in 0.01 M HCl for 10 min at room temperature to inactivate endogenous peroxidases, followed by incubation in 10 mg/ml lysozyme for 1 h at 37°C to permeabilize cell walls. CARD-FISH targeting 16S rRNA with probe DSS658 (Manz et al., 1998) specific for the deltaproteobacterial clade of Desulfosarcina/Desulfococcus was performed as described earlier (Ishii et al., 2004) using Alexa488-labeled tyramide for signal amplification. Following hybridization of 16S rRNA, inactivation of the probe-coupled horseradish peroxidase enzymes (HRP) was achieved by incubations of the filters in 3% H2O2 in 1xPBS for 30 min at RT followed by incubation in 0.1 M HCl for 10 min. RNAs in the cells were digested by incubating the filters in RNase solution (0.5 U μl-1 RNase I, Ambion), 30 μg ml-1 RNase A (Sigma), 0.1 M Tris-HCl pH 8 for 4 – 5 h at 37°C. Filter sections, with either the masD probe mix or negative control probe NonPolyPr350 were incubated in hybridization buffer containing 45% formamide as calculated by the PolyPro software. After initial denaturation at 85°C for 25 min, hybridization lasted for 18 – 22 h at 50°C followed by binding of the anti-Dig HRP-conjugated antibody (Fab fragments) and signal amplification with a Alexa 594-labeled tyramide. Filter sections were embedded in SlowFadeGold antifade reagent (Invitrogen), containing 1 μg ml-1 4′,6-diamidino-2-phenylindole (DAPI). Microscopy was performed on an epifluorescence microscope (Axioplan, Carl Zeiss). To correct for false positive signals, i.e., extracellular signals and non-specific probe binding to particles or filter matrix, we subtracted negative control counts from the total gene counts.
Nucleotide Sequence Accession Numbers
MasD clone sequences from enrichment culture GB4573_14 have been deposited in the EMBL, GenBank and DDBJ nucleotide sequence database under accession numbers LN610408 to LN610424 and HG764719 to HG764728 and MasD sequences from strains PF2803 (DSM 16219), LM2801 (DSM 18843) and the phylotype Propane60GuB under numbers LN879420- LN879422. Raw masD pyrotag sequences have been stored in the sequence read archive under SRA bioproject number 278019.
Definition of OTU Cut-Offs for masD
Thresholds for OTU clustering were calculated by linear correlation of 15 cultured alkane-degrading strains with both 16S rRNA genes and MasD sequence information (Table 2). The taxonomic threshold for microbial species is 98.7% for 16S rRNA genes (Stackebrandt and Ebers, 2006; Yarza et al., 2014), which corresponded to a threshold of 96% for MasD (OTU0.96; Figure 3). The 16S rRNA gene genus-level (94.5%) and family level thresholds (86.5%; Yarza et al., 2014) corresponded to cut-off values for MasD of 86% (OTU0.86) and 72% (OTU0.72), respectively.
FIGURE 3. Pairwise comparison of MasD amino acid and 16S rRNA gene sequence identity for 15 strains capable of anaerobic alkane degradation. Accession numbers of used strains are given in Table 2. The intersection of vertical dashed lines and the regression line gives the cutoff values for defining a species (96%), genus (86%), and family (72%) based on MasD identity corresponding to the proposed values based on 16S rRNA genes (Yarza et al., 2014).
Diversity of MasD Community in Seep Sediments
The diversity of anaerobic alkane-degrading bacteria was studied by pyrosequencing of masD present in sediments from 12 globally distributed stations at methane, gas or hydrocarbon seep sites. We retrieved a total of 12,745 raw sequences. Strict quality filtering of raw sequences was an important part of our data analysis as sequencing errors and chimeras affect OTU clustering of protein sequences more severely than clustering of nucleic acid sequences and would result in an overestimation of environmental MasD diversity. After filtering, 40% of raw reads (5,131 masD sequences, 652 nmsA sequences and 79 bssA sequences) were included in further analyses. After translation and clustering the protein sequences on species-level (OTU0.96) we obtained 420 MasD OTU0.96, 1 NmsA OTU0.96, and 1 BssA OTU0.96 (Table 3). The retrieval of nmsA and bssA sequences from our seep sites was a result of unspecific binding of our masD primers. It is likely then that diversity of bssA and nmsA is not fully covered by these primers and is omitted from further analysis.
MasD Richness and Evenness
Observed richness of alkane degraders ranged between 13 and 93 MasD species-level OTU0.96 (Table 3). Coverage ranged between 83 and 94% indicating sufficient sampling effort (Table 3, Supplementary Figure S1) except for the site with the highest inverse Simpson index, GoM4463, for which the coverage was only 74%. OTU abundance for other taxonomic levels is provided in Supplementary Table S3. After subsampling to standardize sequencing effort, OTU0.96 diversity differed by a maximum factor of 2.25 between sites. Diversity was highest at GoM4463 (27 OTU0.96) and GoM156 (21 OTU0.96) but clearly lower at GB4484 and all methane seep sites (13–18 OTU0.96). Inverse Simpson diversity index (1/D), which takes into account both OTU richness and evenness, ranged between 3 and 9.
Similarity of MasD Communities at Different Seep Sites
The MasD community structure at different seep sites was visualized by NMDS based on amino acid sequences clustered at species level (OTU0.96; Figure 4A). Sites were grouped according to hydrocarbon type at the individual sites, i.e., methane (=methane seeps), gaseous hydrocarbons (=gas seeps), or diverse hydrocarbons (=hydrocarbon seeps) and connected to the weighted centroid of the within-group distances. MasD diversity of samples within these three groups showed a high similarity with a high shared proportion of taxa (31%). Dissimilarity was most pronounced between methane seeps and all other sites. Methane seep sites HR19, NZ315 and Tomm revealed quite similar community structures (43% shared taxa) but were clearly different from those at gas and diverse hydrocarbon seeps (<6% shared taxa). AMV gas seeps host MasD communities partially overlapping with those of the hydrocarbon seeps GB and GoM. The frequency of co-occurrence was highest between the hydrocarbon seeps GoM13 and GoM156, between hydrocarbon and gas seeps GB4484, AMV825 and GoM17 and between methane seeps Tomm and HR19.
FIGURE 4. Similarity of MasD community structure in investigated seep sediments. (A) Similarity of MasD community structure visualized by non-metric multidimensional scaling. Ordination based on species level OTU0.96 dissimilarities (Bray–Curtis). Each sample (dot) is connected to the weighted averaged mean of the within group distances. Ellipses represent one standard deviation of the weighted averaged mean. Dot size reflects Inverse Simpson index (1/D), hence the larger a dot the higher the diversity. The mean percentage of shared MasD OTU0.96 between any two sites was 9%, lowest was 0% and highest was 43%. (B) Abundance and occurrence of MasD OTU0.96. Given abundance is based on number of reads per OTU0.96. OTU0.96 that were present in only 1 of the 12 investigated sediments were defined as locally restricted while those OTU0.96 that were present in at least 10 of the 12 sites were defined as cosmopolitan OTU0.96. (C) Co-occurrence of OTU0.96 among the investigated seep sites. The darker the color the more frequent the co-occurrence.
Hierarchical clustering of co-occurrence frequencies identified two main clusters with cluster A consisting of only gas and hydrocarbon seeps (GoM13, GB4573, GB4484, AMV825, GoM17) and cluster B including all seep types (AMV760, GoM140, NZ315, GoM4463, Tomm, HR19; Figure 4C).
Cosmopolitan and Locally Restricted OTUs
Two out of the 420 species-level OTU0.96 were cosmopolitan as defined by their presence in at least 10 of the 12 stations investigated (Figures 4B and 5). These two OTU comprised the majority of MasD sequences. Fifteen percent of total MasD sequences were assigned to OTU#1 that dominated at almost all sites with 5–55% of quality sequences retrieved from the individual sites. OTU#1 was rare only at methane seep sites NZ315 and Tomm with <2%. Conversely, OTU#2 dominated these two methane seep sites with 53% of total quality reads at NZ315 and 22% at Tomm, in addition to 52% at HR19. The third most abundant OTU#3, which was present at six sites (Supplementary Table S4), strongly dominated in GoM13 (44% of the reads), GoM156 (26%) and AMV825 (14%) sediments. In general, the next ten abundant OTU#4 to #13 were only occasionally abundant and found to be dominant at a single seep site (Supplementary Table S4). For example, OTU#11 comprised only 2.8% of total sequences and were found at five sites but dominated GoM17 (22%) and AMV760 (65%).
FIGURE 5. Network graph displaying the connectivity among sampling sites based on presence-absence of MasD OTU0.96. The different colors indicate the number of investigated seeps that contained an individual OTU. Black dots are locally restricted OTU0.96, the stars are cosmopolitan OTU0.96. OTU are connected to the sampling sites (black polygons) they occur at.
Rare OTU0.96 appearing only once in the whole data set (i.e., absolute single sequence OTUs; SSOabs) are referred to as ‘locally restricted.’ A majority of OTU0.96 (61%) were assigned to SSOabs. Conditionally rare organisms are those appearing only once in a given sample but more often in one or more of the other samples (i.e., relative single sequence OTU0.96, SSOrel, Gobet et al., 2012). The conditionally rare microbial MasD community comprised 16% of total OTU0.96 of which >50% were represented by a low read number (<10 reads per OTU).
Phylogenetic Affiliation of MasD
For phylogenetic analysis, MasD was clustered on a proposed family level of 72% similarity (OTU0.72). In total, 83 family level OTUs were detected. Representative sequences formed three clusters that shared <60% sequence similarity (Figure 6). The three clusters might comprise organisms from different phyla because the inter-cluster identity of <60% is close to the estimated MasD phylum level OTU threshold at 52% (according to the phylum threshold of 75% for 16S rRNA genes; Yarza et al., 2014). The proposed three phyla are mostly consistent with the 16S rRNA taxonomy of the cultured organisms. Strains within cluster I all belong to Proteobacteria, cluster II does not contain any sequences from cultured strains and cluster III includes sequences from Proteobacteria (strain BuS5, Deltaproteobacteria) but also from Firmicutes.
FIGURE 6. Phylogenetic tree showing the affiliations of MasD amino acid sequences retrieved from seep sediments with selected reference sequences. The phylogenetic tree was calculated using the maximum likelihood algorithm (PhyML) with 1000 bootstrap replicates and blosum62 correction considering 95 amino acid positions. The tree was rooted using pyruvate formate lyase (Pfl) as outgroup. Sequences from cultivated alkane degraders or metagenomic analysis were included as additional reference sequences; substrate usage is given within parenthesis. Strains PF2803, LM2801 and enriched phylotype Propane60GuB were added to the tree using parsimony criteria. Nitrate-reducing bacteria are in blue, sulfate-reducing bacteria in red and syntrophic bacteria degrading hydrocarbons under methanogenic conditions in light blue. Sequences from this study are in boldface type. Only one representative amino acid sequence of the individual family level OTUs0.72 is shown. MasD sequences for geneFISH probe design, which were retrieved from a clone library constructed from Guaymas Basin enrichment with butane, are underlined. The scale bar gives 10% estimated sequence divergence. Pie charts represent the portion of sequences from a certain seep type assigned to the cluster. Orange = methane seep, purple = gas seep, green = hydrocarbon seep. Abbreviations: Ass, alkylsuccinate synthase, Mas, 1-methyl alkyl succinate synthase, Bss, benzylsuccinate synthase, Nms, naphthyl-2-methyl-succinate synthase.
OTU richness was highest in cluster I: 372 species-level OTU (OTU0.96) were identified comprising 4099 sequences (80% of total) with 72 total family level OTU (OTU0.72). In cluster II, 47 species-level OTU0.96 were identified comprising 1022 sequences (20% of total) with 10 family level OTU0.72. And finally cluster III contained only 10 sequences clustered into a single OTU on species, genus and family level (Figure 6). For cluster I we defined seven monophyletic subclades: cluster Ia to Ig. Clusters Ia, Ie, and Ig only contain MasD from uncultivated organisms. Clusters Ia and Ib are comprised of a particularly high number of family level OTU0.72 (18 and 16, respectively). Cluster Ib also included the AssA found in a draft genome from Smithella SCADC, a syntrophic deltaproteobacterium derived from different methanogenic alkane-degrading enrichment cultures (Embree et al., 2014; Tan et al., 2014), and Candidatus Magnetoglobus multicellularis (ATBP010000388 on IMG/M). The only cultivated relative in cluster Id is the long-chain alkane- and alkene-degrading Desulfatibacillum aliphaticivorans. This cluster includes the two most abundant and cosmopolitan OTU0.96 #1 and #2. Cluster If is quite diverse and comprised of MasD from betaproteobacterial nitrate reducers and deltaproteobacterial sulfate reducers but also from medium- or short-chain alkane degraders (Figure 6).
In Situ Identification of Alkane-Degrading Bacteria by masD geneFISH
A geneFISH assay for masD was established to identify the alkane-degrading bacterial community in the environment. A sediment-free culture from GB4573 sediments grown under sulfate-reducing conditions with n-butane as sole carbon source was used to optimize the protocol. The enrichment is dominated by alkane-degrading SRB of the SCA2 clade belonging to the deltaproteobacterial Desulfococcus/Desulfosarcina (DSS) branch (Kleindienst et al., 2014). As a basis for probe design, a masD gene library was constructed from this enrichment culture. The obtained sequences fall into all three defined clusters (Figure 6). Because similarity between clusters I, II, and III was too low on the DNA level to allow the design of a single dsDNA polynucleotide probe covering the whole diversity of the gene (Figure 6, Supplementary Table S1), clusters were targeted independently. For cluster II, it was possible to cover all sequences with one dsDNA polynucleotide probe. For cluster I and III, we designed probe mixes of 4 dsDNA polynucleotide probes each for complete coverage. The DNA-based sequence similarity between MasD OTU0.72 from the 12 seeps and the individual probes was above the threshold of 72%, which has been reported to be detected by a 350 bp probe (Moraru et al., 2010). The dsDNA masD probes were synthesized separately and applied as mix on the GB enrichment culture. Hybridization of the enrichment culture showed masD gene signals co-localized with signals of probe DSS658 (Figure 7A), which targets the Desulfosarcina/Desulfococcus branch of Deltaproteobacteria, indicating a successful in situ identification of alkane-degrading bacteria. Members of the DSS made up between 73 and 90% of total DAPI-stained cells. Of these cells, 49–58% showed positive masD gene signals with the masD probe mix (Figure 7B).
FIGURE 7. GeneFISH on alkane-degrading enrichment cultures from Guaymas Basin sediments grown under sulfate-reducing conditions. (A) Fluorescence micrograph showing CARD-FISH stained Desulfosarcina/Desulfococcus cells in green (16S rRNA-targeted probe DSS658) and geneFISH-stained masD-carrying cells in red (masD probe mix). White arrows indicate DSS cells with masD-geneFISH signals; yellow arrow indicates non-specific signal. Scale bar = 5 μm. (B) Relative abundance of DSS and masD-carrying cells in enrichments.
Diversity of Alkane-Degrading Bacteria
To date, only about 20 bacterial strains or phylotypes in enrichments have been found to degrade alkanes anaerobically. They belong to a total of ten families within the phyla Proteobacteria and Firmicutes (Figure 1). In contrast, a quite high bacterial 16S rRNA gene diversity has been reported for hydrocarbon seeps in the Guaymas Basin and Gulf of Mexico, in particular for deltaproteobacterial SRB (Teske et al., 2002; Dhillon et al., 2003; Orcutt et al., 2010; Wegener and Knittel, unpublished data), suggesting that there are likely many more SRB and other bacteria than currently known that thrive anaerobically on alkanes. Our data clearly support this hypothesis as we found a total of 420 MasD species-level OTU0.96 and 83 family level OTU0.72 from the 12 investigated hydrocarbon-impacted sites. Based on rarefaction curves, diversity was not fully covered; therefore increased sequencing efforts might reveal even more diversity. Furthermore, we also might have missed new lineages which did not have the sequence motifs we were searching for. For example, we excluded about 70% of reads retrieved from GoM4463 and HR19 from the analysis only due to the absence of these motifs. The high overall MasD diversity can be explained either by the presence of a microbial seed bank, which comprises dormant alkane-degrading organisms that are resuscitated following environmental changes (Lennon and Jones, 2011) or by the presence of many niches as alkanes are abundant in nature.
The oily sediment site from the northern Gulf of Mexico (GoM4463) had the highest MasD α-diversity (Inverse Simpson Index), followed by site GB4573 in the Guaymas Basin, indicating that these habitats are species-rich because of the presence of very diverse hydrocarbons. This high diversity might also enable the microbes to withstand certain environmental changes, like the strong temperature gradients found at Guaymas Basin. In contrast, the asphaltic sample GoM17 was the least diverse, supporting the notion that in an extreme environment with a rather limited buffet of degradable hydrocarbons, small changes would have serious impact on the microbial communities where few new species accumulate.
Cosmopolitan and Rare Alkane-Degrading Bacteria
Two abundant cosmopolitan species-level MasD OTU0.96 were detected in this study indicating that only a small number of dominant alkane degraders are globally distributed. These OTU were assigned to cluster Id, which included Desulfatibacillum aliphaticivorans str. CV2803, a sulfate-reducing long-chain alkane- (C13–C18) and alkene- (C7–C23) degrading bacterium isolated from hydrocarbon-polluted sediments in the Gulf of Fos (France, Cravo-Laureau et al., 2004a). Their closest relatives were MasD sequences (e.g., Guaymas clone LCA_537) retrieved from an enrichment with dodecane (Kleindienst et al., 2014) suggesting long chain alkanes as substrates for the organisms representing OTU#1 and #2. Their cosmopolitan presence suggests an as of yet unknown environmental importance of members in this group for the degradation of long-chain alkanes and alkenes in marine seep sediments. Cosmopolitan distribution was also shown for other seep-associated organisms, such as certain ANME that have been found in seep sediments worldwide (Ruff et al., 2015).
A tremendous amount of species-level single sequence OTU0.96 (SSOabs and SSOrel) was retrieved from the 12 different sediments. SSOabs are locally restricted, permanently rare species (Gobet et al., 2012). In contrast, SSOrel are organisms that are rare in one ecosystem, but very common or even dominant in another. These organisms may start growth when the conditions change, which was shown for pelagic communities where rare organisms became abundant after disturbance (Sjöstedt et al., 2012) or showed seasonal patterns (Hugoni et al., 2013). Guaymas Basin site 4573 had the most SSOrel and SSOabs supporting the uniqueness of this habitat with strong geochemical gradients, in particular temperature gradients, previously reported (Teske et al., 2002; Teske and Sørensen, 2008). These minor MasD variants have developed under specialized circumstances that may be linked to the diversity and the structure of the alkanes in that environment.
Factors Driving the Alkane-Degrading Community Structure
The MasD community structure was so clearly driven by the hydrocarbon source at the various seeps that the pressure of other factors such as water depth has little to no effect. MasD OTU from methane, gas and hydrocarbon seep sites was most similar within the habitat type as shown by non-metric multidimensional scaling. Thus, the range of available alkanes seems to strongly influence the diversity of MasD-carrying microbes. This is consistent with the narrow substrate range of cultured strains. For example, strain BuS5 oxidizes only propane and butane (Kniemeyer et al., 2007), Azoarcus sp. str. HxN1 oxidizes C6–C8 (Ehrenreich et al., 2000) and Desulfatibacillum aliphaticivorans str. CV2803 oxidizes C13–C18 (Cravo-Laureau et al., 2004a). It has been shown that seep-associated anaerobic methanotrophs (ANME clades) are structured by sediment depth and sediment temperature, (Ruff et al., 2015) while sulfate reducers (SEEP-SRB clades) are structured more by faunal activity and thus biogeochemistry (Felden et al., 2014). In all studies, water depth was not identified as an influencing factor.
Phylogeny of Anaerobic Alkane-Degrading Bacteria
For 16S rRNA genes there are widely accepted thresholds for the definition of a species, genus or family (Rosselló-Móra and Amann, 2015). With the clustering of protein-coding genes, however, we must consider the high variation in sequence conservation. To establish taxonomic levels based on MasD we used cultivated alkane-degrading species and correlated their 16S rRNA nucleotide sequences with their MasD amino acid sequences. This approach has also been implemented for methane-oxidizing bacteria (PmoA, Degelmann et al., 2010) and nitrogen-fixing bacteria (NifH, Bowen et al., 2013). We determined the following threshold values for MasD: 96% amino acid identity for species, 87% for genus, and 72% for families. The species threshold of 96% that we defined for MasD was higher than those determined for PmoA (93%) and NifH (88%), which can in part be explained by the recent re-evaluation and implementation of a new 16S rRNA gene species threshold of 98.7% (Yarza et al., 2014) used in our study compared to the 97% identity threshold used in the above previous studies.
Phylogenetic analysis of MasD resulted in three clusters of which the first was dominant, the second regularly present and the third very rare. The 60% amino acid identity between these clusters is likely in the range of a phylum level threshold for MasD. The only available non-deltaproteobacterial MasD is from phylotype Propane60GuB (cluster III) and showed a maximum of 56% amino acid identity to any sequence of cluster I and cluster II.
Based on the MasD protein sequences alone we cannot exclude that the respective enzymes have similar substrate spectra or affinities. In the MasD phylogenetic tree (Figure 6) there was also no obvious separation between sequences retrieved from freshwater (Callaghan et al., 2010) and those marine habitats (Kleindienst et al., 2012; Acosta-González et al., 2013). Thus, cultivation attempts as well as enzyme characterizations are necessary to evaluate phylum-level functional diversity of MasD.
The dominant cluster I was most diverse including several isolated deltaproteobacterial and betaproteobacterial sulfate- and nitrate-reducing bacteria. It was present at all sites suggesting that these microbes are able to adapt to the local environmental features which differ remarkably between seeps. Within subcluster Ib, eight OTU0.72 cluster tightly with MasD from methanogenic enrichments and Smithella sp. We hypothesize that these MasD belong to syntrophic methanogenic hydrocarbon-degrading communities that are relevant for alkane degradation at marine seeps. However, so far studies on methanogenic alkane-degradation are limited to enrichment cultures (Zengler et al., 1999; Chang et al., 2006; Berdugo-Clavijo and Gieg, 2014; Embree et al., 2014). In situ studies that aim to quantify the responsible syntrophic bacteria and archaea directly in their environment have not been conducted, although the importance of crude oil biodegradation via methanogenesis has been proposed for subsurface oil reservoirs (Jones et al., 2008). The anaerobic biodegradation of octacosane (C28), a solid paraffin, under methanogenic conditions (Davidova et al., 2011; Callaghan, 2013), could be one process of interest for future in situ studies.
MasD of the only cultured short-chain alkane-degrading strain BuS5 (Musat, 2015) grouped into cluster III. Although this strain is found in numbers of ca. 2% of total cell counts at AMV and Guaymas Basin (Kleindienst and Knittel, unpublished data), we did not retrieve any OTU related to BuS5-MasD. Their absence is most likely due to a substantial primer mispairing; the forward primers had 11 (7757f1-f2, 22mer) and 13 (7766f, 23mer) mismatches, respectively, to the BuS5-masD sequence retrieved from the isolate genome (JGI gene ID 2513990058). In addition to deltaproteobacterial BuS5, cluster III is also comprised of two Firmicutes, i.e., a single cell from family Peptococcaceae and a phylotype from the enrichment Propane60GuB, indicative of a lateral gene transfer event. We also report a second potential instance of lateral gene transfer in MasD cluster I, which contains both Beta- and Deltaproteobacteria in subclusters.
Despite the broad diversity of cluster II (10 family level OTU0.72 from this study), this cluster is comprised of only two other sequences from a metagenomic study from Santa Barbara oil seep sediments (Hess, JGI metagenome, IMG submission ID 26744). To date few environmental studies have been conducted and most have investigated fresh water environments (Callaghan et al., 2010; Cheng et al., 2013) and few marine environments (Acosta-González et al., 2013; Kimes et al., 2013; von Netzer et al., 2013; Johnson et al., 2015). All of them implement clone libraries with very limited sequencing depth.
Considering the high overall diversity found within MasD, the range of substrates currently known to be activated by MasD/AssA might expand with further studies. For example, ethane is the most abundant short-chain non-methane alkane at several of our sites (Bazylinski et al., 1988; Boetius et al., 2000; Orcutt et al., 2008; Felden et al., 2013). However, isolation of anaerobic ethane-degrading organisms is lacking. Recently, sequences retrieved from a Gulf of Mexico batch reactor with ethane and sulfate have been shown to be closely related to Desulfosarcina variabilis and strain BuS5 (Bose et al., 2013). Nevertheless, cultivation or enrichment is necessary in order to get detailed insights into the metabolic pathways and substrate spectrum. Large metagenomic studies with systematic screening for masD will help to further unravel MasD diversity and improve current primer sets for future studies. Another approach that could prove useful in future masD marker gene surveys is oligotyping, which uses only the most informative nucleotide positions to define phylogenetically distinct oligotypes (Eren et al., 2013).
In Situ Identification of Alkane-Degrading Communities
In this study we also used a comparative sequence analysis of masD for the development of a polynucleotide probe set. We demonstrated that these probes could be used for the in situ identification of bacteria catalyzing the anaerobic degradation of alkanes by geneFISH. The individual geneFISH probes designed in this study could also be used to target subpopulations, i.e., MasD cluster I, cluster II, and cluster III (Figure 6). geneFISH was originally invented for pmoA and applied on Escherichia coli clones in order to link function with phylogeny (Moraru et al., 2010). Later, the method was applied on thin sections of the hydrothermal vent mussel Bathymodiolus to detect uptake hydrogenases (hupL) in its symbionts (Petersen et al., 2011). Here, we showed that direct application of geneFISH to sediment samples from Guaymas Basin was also possible but challenging due to high background fluorescence.
The in situ quantification of masD will help to identify the size of alkane-degrading communities in the environment independent of cultivation and enrichment and independent from a priori knowledge regarding which taxa are capable of alkane degradation. This method is and will be a major step forward to begin to describe the global relevance and diversity of alkane degradation at natural hydrocarbon seeps and in oil-polluted environments.
As alkane is the predominant component of crude oil, our findings point to a so far overlooked high potential of marine benthic microbes to react to natural changes in hydrocarbon seepage or to massive hydrocarbon input as encountered during anthropogenic oil spills. The number of MasD OTU0.72 that we detected in this study indicated that there may be eight times more family level bacterial clades that are capable of hydrocarbon degradation than previously assumed. This observed high diversity of MasD might allow for higher rates of survival and adaptability within a microbial community exposed to such unstable and changing environmental conditions. Exploring the prevalence and diversity of MasD can help to identify novel lineages of alkane degraders as well as to differentiate closely related phylotypes. Several MasD clusters were discovered that do not yet have cultured representatives. Therefore, one major objective for future studies should be enrichment and isolation of these lineages for further characterization of substrate specialization. Further methodological approaches could include metagenomics and single cell genomics of FACS sorted cells carrying masD genes.
Furthermore, this study represents a first step toward establishing comparative MasD sequencing and masD geneFISH as powerful tools for targeted investigation of the diversity and abundance of alkane-degrading bacteria in anoxic environments.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We greatly acknowledge Antje Boetius for providing samples from expeditions with the research vessels SONNE, METEOR, ALKOR, and SEWARD JOHNSON II. Samples were taken in the framework of the GEOTECHNOLOGIEN programs MUMM I and II (grants 03G0554A and 03G0608A) funded by the German Ministry of Education and Research (BMBF) and the German Research foundation and by the EU 5th FP METROL and 6th FP HERMES. Andreas Teske and Marc Mussmann are acknowledged for providing samples from an expedition with ATLANTIS. Special thanks go to Sara Kleindienst who initiated the enrichment cultures, Florin Musat for fruitful discussions and sharing unpublished MasD sequences, Jimena Barrero Canosa for introduction and technical support on geneFISH, Hanno Teeling for help with 454 tag sequencing analysis as well as Manabu Fukui and Hisaya Kojima for providing the MasD sequence of strain PL12. We are grateful to Meghan Chafee for critical reading, valuable suggestions and language editing. SER was supported by the Leibniz program awarded to Antje Boetius. Further support was provided by the Max Planck Society, Germany.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2015.01511
Acosta-González, A., Rosselló-Móra, R., and Marqués, S. (2013). Characterization of the anaerobic microbial community in oil-polluted subtidal sediments: aromatic biodegradation potential after the Prestige oil spill. Environ. Microbiol. 15, 77–92. doi: 10.1111/j.1462-2920.2012.02782.x
Adams, M. M., Hoarfrost, A. L., Bose, A., Joye, S. B., and Girguis, P. R. (2013). Anaerobic oxidation of short-chain alkanes in hydrothermal sediments: potential influences on sulfur cycling and microbial diversity. Front. Microbiol. 4:110. doi: 10.3389/fmicb.2013.00110
Aitken, C. M., Jones, D. M., Maguire, M. J., Gray, N. D., Sherry, A., Bowler, B. F. J., et al. (2013). Evidence that crude oil alkane activation proceeds by different mechanisms under sulfate-reducing and methanogenic conditions. Geochim. Cosmochim. Acta 109, 162–174. doi: 10.1016/j.gca.2013.01.031
Bazylinski, D. A., Farrington, J. W., and Jannasch, H. W. (1988). Hydrocarbons in surface sediments from a Guaymas Basin hydrothermal vent site. Org. Geochem. 12, 547–558. doi: 10.1016/0146-6380(88)90146-5
Becker, A., Fritz-Wolf, K., Kabsch, W., Knappe, J., Schultz, S., and Wagner, A. F. V. (1999). Structure and mechanism of the glycyl radical enzyme pyruvate formate-lyase. Nat. Struct. Biol. 6, 969–975. doi: 10.1038/13341
Berdugo-Clavijo, C., and Gieg, L. M. (2014). Conversion of crude oil to methane by a microbial consortium enriched from oil reservoir production waters. Front. Microbiol. 5:197. doi: 10.3389/fmicb.2014.00197
Boetius, A., Ravenschlag, K., Schubert, C. J., Rickert, D., Widdel, F., Gieseke, A., et al. (2000). A marine microbial consortium apparently mediating anaerobic oxidation of methane. Nature 407, 623–626. doi: 10.1038/35036572
Bose, A., Rogers, D. R., Adams, M. M., Joye, S. B., and Girguis, P. R. (2013). Geomicrobiological linkages between short-chain alkane consumption and sulfate reduction rates in seep sediments. Front. Microbiol. 4:386. doi: 10.3389/fmicb.2013.00386
Bowen, J. L., Byrnes, J. E. K., Weisman, D., and Colaneri, C. (2013). Functional gene pyrosequencing and network analysis: an approach to examine the response of denitrifying bacteria to increased nitrogen supply in salt marsh sediments. Front. Microbiol. 4:342. doi: 10.3389/fmicb.2013.00342
Bowles, M. W., Samarkin, V. A., Bowles, K. M., and Joye, S. B. (2011). Weak coupling between sulfate reduction and the anaerobic oxidation of methane in methane-rich seafloor sediments during ex situ incubation. Geochim. Cosmochim. Acta 75, 500–519. doi: 10.1016/j.gca.2010.09.043
Callaghan, A. V., Davidova, I. A., Savage-Ashlock, K., Parisi, V. A., Gieg, L. M., Suflita, J. M., et al. (2010). Diversity of benzyl- and alkylsuccinate synthase genes in hydrocarbon-impacted environments and enrichment cultures. Environ. Sci. Technol. 44, 7287–7294. doi: 10.1021/es1002023
Callaghan, A. V., Gieg, L. M., Kropp, K. G., Suflita, J. M., and Young, L. Y. (2006). Comparison of mechanisms of alkane metabolism under sulfate-reducing conditions among two bacterial isolates and a bacterial consortium. Appl. Environ. Microbiol. 72, 4274–4282. doi: 10.1128/AEM.02896-05
Chang, Y. T., Lee, J. F., Chao, H. P., and Liao, W. L. (2006). Bacterial community changes with N′-N′ dimethylforamide (DMF) additives during polycyclic aromatic hydrocarbons (PAH) biodegardation. Environ. Technol. 27, 1–14. doi: 10.1080/09593332708618613
Cheng, L., Rui, J., Li, Q., Zhang, H., and Lu, Y. (2013). Enrichment and dynamics of novel syntrophs in a methanogenic hexadecane-degrading culture from a Chinese oilfield. FEMS Microbiol. Ecol. 83, 757–766. doi: 10.1111/1574-6941.12031
Cravo-Laureau, C., Labat, C., Joulian, C., Matheron, R., and Hirschler-Réa, A. (2007). Desulfatiferula olefinivorans gen. nov., sp. nov., a long-chain n-alkene-degrading, sulfate-reducing bacterium. Int. J. Syst. Evol. Microbiol. 57, 2699–2702. doi: 10.1099/ijs.0.65240-0
Cravo-Laureau, C., Matheron, R., Cayol, J. L., Joulian, C., and Hirschler-Réa, A. (2004a). Desulfatibacillum aliphaticivorans gen. nov., sp nov., an n-alkane- and n-alkene-degrading, sulfate-reducing bacterium. Int. J. Syst. Evol. Microbiol. 54, 77–83. doi: 10.1099/ijs.0.63104-0
Cravo-Laureau, C., Matheron, R., Joulian, C., Cayol, J. L., and Hirschler-Rea, A. (2004b). Desulfatibacillum alkenivorans sp nov., a novel n-alkene-degrading, sulfate-reducing bacterium, and emended description of the genus Desulfatibacillum. Int. J. Syst. Evol. Microbiol. 54, 1639–1642. doi: 10.1099/ijs.0.63104-0
Davidova, I., Callaghan, A., Duncan, K., Sunner, J., Biri, B., Wawrik, B., et al. (2011). “Long-chain paraffin metabolism by a methanogenic bacterial consortium enriched from marine sediments (poster presentation),” in Proceedings of the 8th International Symposium of Subsurface Microbiology, Garmisch-Partenkirchen.
Davidova, I. A., Duncan, K. E., Choi, O. K., and Suflita, J. M. (2006). Desulfoglaeba alkanexedens gen. nov., sp nov., an n-alkane-degrading, sulfate-reducing bacterium. Int. J. Syst. Evol. Microbiol. 56, 2737–2742. doi: 10.1099/ijs.0.64398-0
Degelmann, D. M., Borken, W., Drake, H. L., and Kolb, S. (2010). Different atmospheric methane-oxidizing communities in European beech and Norway spruce soils. Appl. Environ. Microbiol. 76, 3228–3235. doi: 10.1128/AEM.02730-09
Dhillon, A., Teske, A., Dillon, J., Stahl, D. A., and Sogin, M. L. (2003). Molecular characterization of sulfate-reducing bacteria in the Guaymas Basin. Appl. Environ. Microbiol. 69, 2765–2772. doi: 10.1128/AEM.69.5.2765-2772.2003
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Embree, M., Nagarajan, H., Movahedi, N., Chitsaz, H., and Zengler, K. (2014). Single-cell genome and metatranscriptome sequencing reveal metabolic interactions of an alkane-degrading methanogenic community. ISME J. 8, 757–767. doi: 10.1038/ismej.2013.187
Eren, A. M., Maignien, L., Sul, W. J., Murphy, L. G., Grim, S. L., Morrison, H. G., et al. (2013). Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol. Evol. 4, 1111–1119. doi: 10.1111/2041-210X.12114
Felden, J., Lichtschlag, A., Wenzhöfer, F., De Beer, D., Feseker, T., Pop Ristova, P., et al. (2013). Limitations of microbial hydrocarbon degradation at the Amon mud volcano (Nile deep-sea fan). Biogeosciences 10, 3269–3283. doi: 10.5194/bg-10-3269-2013
Felden, J., Ruff, S. E., Ertefai, T., Inagaki, F., Hinrichs, K. U., and Wenzhöfer, F. (2014). Anaerobic methanotrophic community of a 5346-m-deep vesicomyid clam colony in the Japan Trench. Geobiology 12, 183–199. doi: 10.1111/gbi.12078
Gobet, A., Boer, S. I., Huse, S. M., Van Beusekom, J. E. E., Quince, C., Sogin, M. L., et al. (2012). Diversity and dynamics of rare and of resident bacterial populations in coastal sands. ISME J. 6, 542–553. doi: 10.1038/ismej.2011.132
Grundmann, O., Behrends, A., Rabus, R., Amann, J., Halder, T., Heider, J., et al. (2008). Genes encoding the candidate enzyme for anaerobic activation of n-alkanes in the denitrifying bacterium, strain HxN1. Environ. Microbiol. 10, 376–385. doi: 10.1111/j.1462-2920.2007.01458.x
Grünke, S., Felden, J., Lichtschlag, A., Girnth, A. C., De Beer, D., Wenzhöfer, F., et al. (2011). Niche differentiation among mat-forming, sulfide-oxidizing bacteria at cold seeps of the Nile Deep Sea Fan (Eastern Mediterranean Sea). Geobiology 9, 330–348. doi: 10.1111/j.1472-4669.2011.00281.x
Hugoni, M., Taib, N., Debroas, D., Domaizon, I., Jouan Dufournel, I., Bronner, G., et al. (2013). Structure of the rare archaeal biosphere and seasonal dynamics of active ecotypes in surface coastal waters. Proc. Natl. Acad. Sci. U.S.A. 110, 6004–6009. doi: 10.1073/pnas.1216863110
Ishii, K., Mußmann, M., Mac Gregor, B. J., and Amann, R. (2004). An improved fluorescence in situ hybridization protocol for the identification of bacteria and archaea in marine sediments. FEMS Microbiol. Ecol. 50, 203–212. doi: 10.1016/j.femsec.2004.06.015
Jaekel, U., Musat, N., Adam, B., Kuypers, M., Grundmann, O., and Musat, F. (2012). Anaerobic degradation of propane and butane by sulfate-reducing bacteria enriched from marine hydrocarbon cold seeps. ISME J. 7, 885–895. doi: 10.1038/ismej.2012.159
Jaekel, U., Zedelius, J., Wilkes, H., and Musat, F. (2015). Anaerobic degradation of cyclohexane by sulfate-reducing bacteria from hydrocarbon-contaminated marine sediments. Front. Microbiol. 6:116. doi: 10.3389/fmicb.2015.00116
Johnson, J. M., Wawrik, B., Isom, C., Boling, W. B., and Callaghan, A. V. (2015). Interrogation of Chesapeake Bay sediment microbial communities for intrinsic alkane-utilizing potential under anaerobic conditions. FEMS Microbiol. Ecol. 91, 1–14. doi: 10.1093/femsec/fiu035
Jones, D. M., Head, I. M., Gray, N. D., Adams, J. J., Rowan, A. K., Aitken, C. M., et al. (2008). Crude-oil biodegradation via methanogenesis in subsurface petroleum reservoirs. Nature 451, 176–180. doi: 10.1038/nature06484
Khelifi, N., Amin Ali, O., Roche, P., Grossi, V., Brochier-Armanet, C., Valette, O., et al. (2014). Anaerobic oxidation of long-chain n-alkanes by the hyperthermophilic sulfate-reducing archaeon, Archaeoglobus fulgidus. ISME J. 8, 2153–2166. doi: 10.1038/ismej.2014.58
Kimes, N. E., Callaghan, A. V., Aktas, D. F., Smith, W. L., Sunner, J., Golding, B. T., et al. (2013). Metagenomic analysis and metabolite profiling of deep-sea sediments from the Gulf of Mexico following the Deepwater Horizon oil spill. Front. Microbiol. 4:50. doi: 10.3389/fmicb.2013.00050
Kleindienst, S., Herbst, F.-A., Stagars, M., Von Netzer, F., Von Bergen, M., Seifert, J., et al. (2014). Diverse sulfate-reducing bacteria of the Desulfosarcina/Desulfococcus clade are the key alkane degraders at marine seeps. ISME J. 8, 2029–2044. doi: 10.1038/ismej.2014.51
Kleindienst, S., Ramette, A., Amann, R., and Knittel, K. (2012). Distribution and in situ abundance of sulfate-reducing bacteria in diverse marine hydrocarbon seep sediments. Environ. Microbiol. 14, 2689–2710. doi: 10.1111/j.1462-2920.2012.02832.x
Kniemeyer, O., Musat, F., Sievert, S. M., Knittel, K., Wilkes, H., Blumenberg, M., et al. (2007). Anaerobic oxidation of short-chain hydrocarbons by marine sulphate-reducing bacteria. Nature 449, 898–910. doi: 10.1038/nature06200
Knittel, K., Boetius, A., Lemke, A., Eilers, H., Lochte, K., Pfannkuche, O., et al. (2003). Activity, distribution, and diversity of sulfate reducers and other bacteria in sediments above gas hydrate (Cascadia margin, Oregon). Geomicrobiol. J. 20, 269–294. doi: 10.1080/01490450303896
Kropp, K. G., Davidova, I. A., and Suflita, J. M. (2000). Anaerobic oxidation of n-dodecane by an addition reaction in a sulfate-reducing bacterial enrichment culture. Appl. Environ. Microbiol. 66, 5393–5398. doi: 10.1128/AEM.66.12.5393-5398.2000
Manz, W., Eisenbrecher, M., Neu, T. R., and Szewzyk, U. (1998). Abundance and spatial organization of gram-negative sulfate-reducing bacteria in activated sludge investigated by in situ probing with specific 16S rRNA targeted oligonucleotides. FEMS Microbiol. Ecol. 25, 43–61. doi: 10.1111/j.1574-6941.1998.tb00459.x
Markowitz, V. M., Chen, I.-M. A., Chu, K., Szeto, E., Palaniappan, K., Pillay, M., et al. (2014). IMG/M 4 version of the integrated metagenome comparative analysis system. Nucleic Acids Res. 42, D568–D573. doi: 10.1093/nar/gkt919
Mastalerz, V., De Lange, G. J., and Dählmann, A. (2009). Differential aerobic and anaerobic oxidation of hydrocarbon gases discharged at mud volcanoes in the Nile deep-sea fan. Geochim. Cosmochim. Acta 73, 3849–3863. doi: 10.1016/j.gca.2008.12.030
Mbadinga, S. M., Wang, L.-Y., Zhou, L., Liu, J.-F., Gu, J.-D., and Mu, B.-Z. (2011). Microbial communities involved in anaerobic degradation of alkanes. Int. Biodeter. Biodegradation 65, 1–13. doi: 10.1016/j.ibiod.2010.11.009
Moraru, C., Lam, P., Fuchs, B. M., Kuypers, M. M. M., and Amann, R. (2010). GeneFISH – an in situ technique for linking gene presence and cell identity in environmental microorganisms. Environ. Microbiol. 12, 3057–3073. doi: 10.1111/j.1462-2920.2010.02281.x
Moraru, C., Moraru, G., Fuchs, B. M., and Amann, R. (2011). Concepts and software for a rational design of polynucleotide probes. Environ. Microbiol. Rep. 3, 69–78. doi: 10.1111/j.1758-2229.2010.00189.x
Musat, F., Wilkes, H., Behrends, A., Woebken, D., and Widdel, F. (2010). Microbial nitrate-dependent cyclohexane degradation coupled with anaerobic ammonium oxidation. ISME J. 4, 1290–1301. doi: 10.1038/ismej.2010.50
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., O’Hara, R. B., Simpson, G. L., et al. (2012). vegan: Community Ecology Package. Available at: http://CRAN.R-project.org/package=vegan.
Orcutt, B., Samarkin, V., Boetius, A., and Joye, S. B. (2008). On the relationship between methane production and oxidation by anaerobic methanotrophic communities from cold seeps of the Gulf of Mexico. Environ. Microbiol. 10, 1108–1117. doi: 10.1111/j.1462-2920.2007.01526.x
Orcutt, B. N., Joye, S. B., Kleindienst, S., Knittel, K., Ramette, A., Reitz, A., et al. (2010). Impact of natural oil and higher hydrocarbons on microbial diversity, distribution, and activity in Gulf of Mexico cold-seep sediments. Deep Sea Res. Part II 57, 2008–2021. doi: 10.1016/j.dsr2.2010.05.014
Petersen, J. M., Zielinski, F. U., Pape, T., Seifert, R., Moraru, C., Amann, R., et al. (2011). Hydrogen is an energy source for hydrothermal vent symbioses. Nature 476, 176–180. doi: 10.1038/nature10325
Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., and Peplies, J. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864
Ruff, S. E., Arnds, J., Knittel, K., Amann, R., Wegener, G., Ramette, A., et al. (2013). Microbial communities of deep-sea methane seeps at Hikurangi continental margin (New Zealand). PLoS ONE 8:e72627. doi: 10.1371/journal.pone.0072627
Ruff, S. E., Biddle, J. F., Teske, A. P., Knittel, K., Boetius, A., and Ramette, A. (2015). Global dispersion and local diversification of the methane seep microbiome. Proc. Natl. Acad. Sci. U.S.A. 112, 4015–4020. doi: 10.1073/pnas.1421865112
Savage, K. N., Krumholz, L. R., Gieg, L. M., Parisi, V. A., Suflita, J. M., Allen, J., et al. (2010). Biodegradation of low-molecular-weight alkanes under mesophilic, sulfate-reducing conditions: metabolic intermediates and community patterns. FEMS Microbiol. Ecol. 72, 485–495. doi: 10.1111/j.1574-6941.2010.00866.x
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Sjöstedt, J., Koch-Schmidt, P., Pontarp, M., Canbäck, B., Tunlid, A., Lundberg, P., et al. (2012). Recruitment of members from the rare biosphere of marine bacterioplankton communities after an environmental disturbance. Appl. Environ. Microbiol. 78, 1361–1369. doi: 10.1128/AEM.05542-11
Suess, E., Torres, M., Bohrmann, G., Collier, R. W., Greinert, J., Linke, P., et al. (1999). Gas hydrate destabilization: enhanced dewatering, benthic material turnover and large methane plumes at the Cascadia convergent margin. Earth Planet. Sci. Lett. 170, 1–15. doi: 10.1016/S0012-821X(99)00092-8
Tan, B., Nesbo, C., and Foght, J. (2014). Re-analysis of omics data indicates Smithella may degrade alkanes by addition to fumarate under methanogenic conditions. ISME J. 8, 2353–2356. doi: 10.1038/ismej.2014.87
Teske, A., Hinrichs, K.-U., Edgcomb, V., De Vera Gomez, A., Kysela, D., Sylva, S. P., et al. (2002). Microbial diversity of hydrothermal sediments in the Guaymas Basin: evidence for anaerobic methanotrophic communities. Appl. Environ. Microbiol. 68, 1994–2007. doi: 10.1128/AEM.68.4.1994-2007.2002
Treude, T., Boetius, A., Knittel, K., Wallmann, K., and Jørgensen, B. B. (2003). Anaerobic oxidation of methane above gas hydrates at Hydrate Ridge, NE Pacific Ocean. Mar. Ecol. Prog. Ser. 264, 1–14. doi: 10.3354/meps264001
von Netzer, F., Pilloni, G., Kleindienst, S., Krüger, M., Knittel, K., Gründger, F., et al. (2013). Enhanced gene detection assays for fumarate-adding enzymes allow uncovering of anaerobic hydrocarbon degraders in terrestrial and marine systems. Appl. Environ. Microbiol. 79, 543–552. doi: 10.1128/AEM.02362-12
Wang, L.-Y., Li, W., Mbadinga, S. M., Liu, J.-F., Gu, J.-D., and Mu, B.-Z. (2012). Methanogenic microbial community composition of oily sludge and its enrichment amended with alkanes incubated for over 500 days. Geomicrobiol. J. 29, 716–726. doi: 10.1080/01490451.2011.619634
Weelink, S. A. B., Van Doesburg, W., Saia, F. T., Rijpstra, W. I. C., Röling, W. F. M., Smidt, H., et al. (2009). A strictly anaerobic betaproteobacterium Georgfuchsia toluolica gen. nov., sp. nov. degrades aromatic compounds with Fe(III), Mn(IV) or nitrate as an electron acceptor. FEMS Microbiol. Ecol. 70, 575–585. doi: 10.1111/j.1574-6941.2009.00778.x
Wegener, G., Shovitri, M., Knittel, K., Niemann, H., Hovland, M., and Boetius, A. (2008). Biogeochemical processes and microbial diversity of the Gullfaks and Tommeliten methane seeps (Northern North Sea). Biogeosciences 5, 1127–1144. doi: 10.5194/bg-5-1127-2008
Widdel, F., and Bak, F. (1992). “Gram-negative mesophilic sulfate-reducing bacteria,” in The Prokaryotes, 2nd Edn, eds A. Balows, H. G. Truper, M. Dworkin, W. Harder, and K. H. Schleifer (New York, NY: Springer-Verlag).
Widdel, F., Knittel, K., and Galushko, A. (2010). “Anaerobic hydrocarbon-degrading microorganisms: an overview,” in Handbook of Hydrocarbon and Lipid Microbiology, eds K. N. Timmis, T. McGenity, J. R. Van Der Meer, and V. De Lorenzo (Heidelberg: Springer), 1997–2021.
Widdel, F., and Pfennig, N. (1982). Studies on dissimilatory sulfate-reducing bacteria that decompose fatty acids II. Incomplete oxidation of propionate by Desulfobulbus propionicus gen. nov., sp. nov. Arch. Microbiol. 131, 360–365. doi: 10.1007/BF00411187
Wilkes, H., Rabus, R., Fischer, T., Armstroff, A., Behrends, A., and Widdel, F. (2002). Anaerobic degradation of n-hexane in a denitrifying bacterium: further degradation of the initial intermediate (1-methylpentyl)succinate via C-skeleton rearrangement. Arch. Microbiol. 177, 235–243. doi: 10.1007/s00203-001-0381-3
Yarza, P., Yilmaz, P., Pruesse, E., Glockner, F. O., Ludwig, W., Schleifer, K.-H., et al. (2014). Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat. Rev. Microbiol. 12, 635–645. doi: 10.1038/nrmicro3330
Zhou, L., Li, K.-P., Mbadinga, S., Yang, S.-Z., Gu, J.-D., and Mu, B.-Z. (2012). Analyses of n-alkanes degrading community dynamics of a high-temperature methanogenic consortium enriched from production water of a petroleum reservoir by a combination of molecular techniques. Ecotoxicology 21, 1680–1691. doi: 10.1007/s10646-012-0949-5
Keywords: alkyl succinate synthase, MasD, AssA, sulfate reduction, anaerobic alkane oxidation, microbial diversity, geneFISH
Citation: Stagars MH, Ruff SE, Amann R and Knittel K (2016) High Diversity of Anaerobic Alkane-Degrading Microbial Communities in Marine Seep Sediments Based on (1-methylalkyl)succinate Synthase Genes. Front. Microbiol. 6:1511. doi: 10.3389/fmicb.2015.01511
Received: 30 July 2015; Accepted: 14 December 2015;
Published: 07 January 2016.
Edited by:Hans H. Richnow, Helmholtz Centre for Environmental Research, Germany
Reviewed by:Beth Orcutt, Bigelow Laboratory for Ocean Sciences, USA
Zhidan Liu, China Agricultural University, China
Copyright © 2016 Stagars, Ruff, Amann and Knittel. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Katrin Knittel, email@example.com
†Present address: S. Emil Ruff, University of Calgary, 2500 University Drive NW, Calgary, AB T2N 1N4, Canada