Original Research ARTICLE
Proteomic Stable Isotope Probing Reveals Biosynthesis Dynamics of Slow Growing Methane Based Microbial Communities
- 1Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA, USA
- 2Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, TN, USA
Marine methane seep habitats represent an important control on the global flux of methane. Nucleotide-based meta-omics studies outline community-wide metabolic potential, but expression patterns of environmentally relevant proteins are poorly characterized. Proteomic stable isotope probing (proteomic SIP) provides additional information by characterizing phylogenetically specific, functionally relevant activity in mixed microbial communities, offering enhanced detection through system-wide product integration. Here we applied proteomic SIP to 15 and CH4 amended seep sediment microcosms in an attempt to track protein synthesis of slow-growing, low-energy microbial systems. Across all samples, 3495 unique proteins were identified, 11% of which were 15N-labeled. Consistent with the dominant anaerobic oxidation of methane (AOM) activity commonly observed in anoxic seep sediments, proteins associated with sulfate reduction and reverse methanogenesis—including the ANME-2 associated methylenetetrahydromethanopterin reductase (Mer)—were all observed to be actively synthesized (15N-enriched). Conversely, proteins affiliated with putative aerobic sulfur-oxidizing epsilon- and gammaproteobacteria showed a marked decrease over time in our anoxic sediment incubations. The abundance and phylogenetic range of 15N-enriched methyl-coenzyme M reductase (Mcr) orthologs, many of which exhibited novel post-translational modifications, suggests that seep sediments provide niches for multiple organisms performing analogous metabolisms. In addition, 26 proteins of unknown function were consistently detected and actively expressed under conditions supporting AOM, suggesting that they play important roles in methane seep ecosystems. Stable isotope probing in environmental proteomics experiments provides a mechanism to determine protein durability and evaluate lineage-specific responses in complex microbial communities placed under environmentally relevant conditions. Our work here demonstrates the active synthesis of a metabolically specific minority of enzymes, revealing the surprising longevity of most proteins over the course of an extended incubation experiment in an established, slow-growing, methane-impacted environmental system.
Marine methane seep sediments harbor complex microbial communities that play significant roles in the global carbon and sulfur biogeochemical cycles (Jørgensen and Kasten, 2006; Reeburgh, 2007). Studies of these features have predominately focused on biodiversity (e.g., Bidle et al., 1999; Knittel et al., 2005; Pernthaler et al., 2008), the energetic and biochemical basis of syntrophic partnerships (Alperin and Hoehler, 2009; Stams and Plugge, 2009), and ecosystem-wide contributions to global methane processing (Reeburgh, 2007). Such studies revealed that one of the dominant metabolisms at seep sites—the anaerobic oxidation of methane (AOM)—is enacted by consortia of anaerobic methanotrophs (ANME) most closely related to methanogens and sulfate-reducing deltaproteobacteria. Phylogenetically-linked measures of anabolic activity and growth of AOM consortia using stable isotope probing have identified intricately coupled microbial metabolisms (Krüger et al., 2008a; Dekas et al., 2009; Orphan et al., 2009), while metagenomics and metatranscriptomics studies have pointed to metabolic potential and associated biochemical pathways (Hallam et al., 2004; Pernthaler et al., 2008; Meyerdierks et al., 2010; Stokke et al., 2012; Wang et al., 2014). Proteomic stable isotope probing (SIP) couples these experimental objectives by offering a functionally- and phylogenetically-constrained enzymatic profile of constituent organisms as well as a temporal reporter of protein synthesis and metabolic response to distinct conditions (Pan et al., 2011; Seifert et al., 2012; Justice et al., 2014; Mohr et al., 2014).
In environments such as anoxic methane seep sediment and an array of subsurface habitats, where energy can become limiting, microbes frequently exhibit extremely slow growth rates and are particularly recalcitrant to activity-based analyses. Radiotracers are sensitive probes to rates on the order of 10−19–10−16 mol cell−1 day−1 (Parkes et al., 1990), but the range of discoverable metabolic reactions is severely constrained. Fluorescence in situ hybridization (FISH) coupled with nanoscale secondary ion mass spectrometry (nanoSIMS; e.g., Wagner, 2009) can detect assimilation rates as low as 10−17 mol cell−1 day−1 (Morono et al., 2014) and visualize phylogenetically constrained spatial associations, though experimental throughput is low and only assimilatory processes can be queried. Whole-cell bioorthogonal non-canonical amino acid tagging (BONCAT) coupled with FISH can be used to fluorescently visualize microorganisms active in protein synthesis (Hatzenpichler et al., 2014), but identification of specific newly synthesized proteins has only been successfully applied in a few cases (Babin et al., 2016) and requires further development for complex environmental systems. Proteomic SIP represents an important entry in the analysis of metabolic activity in low-energy microbial systems, due to its spatially broad, yet functionally- and phylogenetically-specific search space. The procedure is able to identify particular metabolic pathways or enzyme-mediated responses (Bozinovski et al., 2012; Justice et al., 2014) that can be integrated across constituents of a particular lineage, offering an opportunity to access a segment of the low-activity biosphere that might go undetected by other methods due to low levels of anabolism by individual organisms. Although challenges remain—particularly surrounding protein extraction, peptide quantification, database collation, and the interpretation of non-detections—proteomic SIP offers a promising method for assessing the in vivo activity and catalytic function of microbial communities.
Culture-independent studies of energy-limited methane seep settings have included meta-omics investigations focused largely on the pathway responsible for AOM. These investigations have revealed that AOM likely utilizes the same enzymes as methanogenesis, operating in the reverse direction (Hallam et al., 2004; Meyerdierks et al., 2010). ANME-1 draft genomes and fosmids, however, lack the mer gene (Meyerdierks et al., 2010; Stokke et al., 2012), prompting the proposition of a reverse-methanogenesis bypass (Meyerdierks et al., 2010). ANME-2 lineages, including ANME-2a (Wang et al., 2014), and nitrate-reducing ANME-2d (Haroon et al., 2013) genomes, as well as a magnetic enrichment of ANME-2c consortia (Pernthaler et al., 2008), contained the mer gene. Corresponding gene expression profiles revealed ANME-1 methylenetetrahydromethanopterin dehydrogenase (mtd), heterodisulfide reductase subunits A and B (hdrAB), and methyl-coenzyme M reductase subunit A (mcrA) transcripts (Meyerdierks et al., 2010). ANME-2a and ANME-2d transcriptomic datasets exhibited highly expressed methyl-tetrahydromethanopterin coenzyme M methyltransferase (mtr), mcr, mer, and methenyltetrahydromethanopterin cyclohydrolase (mch) and substantially lower levels of mtd and several formylmethanofuran dehydrogenase (fmd) subunits (Haroon et al., 2013; Wang et al., 2014).
Proteomic SIP in environmental samples enables the identification of isotopically enriched proteins synthesized by a microbial assemblage after incubation with an isotopically labeled substrate (e.g., 15N-ammonium). The application of proteomic SIP to slow growing methane seep microbial communities offers an opportunity to examine production patterns of reverse methanogenesis enzymes central to ANME-facilitated methane-oxidation, as well as to identify the proteins expressed by syntrophic bacterial partners and other organisms within the seep ecosystem. This approach can highlight metabolic shifts resulting from incubation conditions, reveal the fraction of proteins that are newly synthesized, and potentially facilitate the construction of a framework for interspecies metabolic relationships (Jehmlich et al., 2008; Hawley et al., 2014). The greatest insights from environmental metaproteomic studies to date have been associated with time-resolved analyses of low-complexity microbial systems (i.e., acid mine drainage microbial biofilms; Denef et al., 2010; Pan et al., 2011), while similar efforts with sponge symbionts (Liu et al., 2012), sinking marine particulates (Moore et al., 2012), and the human gut (Verberkmoes et al., 2009; Xiong et al., 2015) have clarified interspecies relationships and dominant biogeochemical functionalities. There have been comparatively few studies applying proteomics to deep-sea methane seep ecosystems; these have focused on individual enzymes (Krüger et al., 2003), explored specific biochemical questions (Glass et al., 2014), or achieved relatively low coverage of the full ecosystem (Stokke et al., 2012). Here we use proteomics and proteomic SIP to examine multiple methane seep sediment incubations containing diverse ANME and bacterial lineages and exhibiting a range of methanotrophic activity levels.
Sample Collection from Methane Seeps
Seafloor sediment was collected from active methane seeps at Hydrate Ridge, Oregon, an area that has served as a natural laboratory for the study of AOM for over a decade (Suess et al., 1999; Treude et al., 2003; Boetius and Suess, 2004; Marlow et al., 2014a). At Hydrate Ridge, the sulfate-methane transition zone occurs in the top several centimeters of seafloor sediment; within this horizon, microbial aggregate abundance, methane oxidation rates, and sulfate reduction rates are highest (Boetius et al., 2000; Boetius and Suess, 2004). Methane concentrations at active seep sediments have been measured and modeled at values up to 70 mM (Boetius and Suess, 2004) and 50 mM (Tryon et al., 2002), respectively; sulfate concentrations decrease with sediment depth, concomitant with augmented sulfide levels up to 26 mM (Sahling et al., 2002). Samples for this study were collected during two field campaigns, including R/V Atlantis leg AT-15-68 using the DSV Alvin (August 2010) and Atlantis leg AT-18-10 using the ROV Jason II (September 2011). Samples are referenced according to unique four-digit serial numbers. Push core (PC) 16 from Alvin dive 4635 at Hydrate Ridge South (44° 34.20′ N, 125° 8.87′ W, 775 m depth) was collected from a microbial mat indicative of active methane seepage. The recovered 12 cm of sediment from PC-16 was divided into the upper 0-6 cm (sample #3730, used for metagenomic sequencing) and the lower 6-12 cm (#3731, used for metaproteomic and nanoSIMS-based analyses). #5133 and #5579 refer to elevator 3A push cores 47 and 41, respectively, collected at Hydrate Ridge North during Jason II dive J2-593 from an area of active seepage marked by bubble ebullition and surface expression of white microbial mat (44° 40.17′ N, 125° 5.89′ W, 600 m depth). #5133 recovered 9 cm of sediment, with a dark sulfidic zone at ~4-5 cm depth; #5579 contained 12 cm of sediment.
Shipboard, push cores were immediately transferred to a 4°C walk-in cold room and processed within several hours. To prepare sediment for future experimentation, compacted sediment was stored in anoxic, Ar-flushed mylar bags at 4°C until use. Incubations prepared for metaproteomic analysis were amended with 1 mM of 14 or 15. Inhibition of microbial activity due to 15 was unlikely, as methanogen and sulfate reducing bacterial pure cultures, as well as methane seep sediment, showed no adverse effects when presented with 15N ammonia ratios up to 100% (Krüger et al., 2008b).
Laboratory-Based 15N-SIP Incubations
In advance of experimental set-up, separate sediment slurries from #3731 and #5133 were reconstituted under anoxic conditions using 0.22 μm filtered, anoxic N2-sparged Hydrate Ridge bottom water (at a 1:2 sediment:bottom water ratio by volume) and these mesocosm incubations (1L bottle) were maintained with a 2 bar CH4 headspace for 1 month. The concentration of the filtered bottom water was 496 μM, as measured with a Dionex DX-500 Ion Chromatograph using the protocol specified in Green-Saxena et al. (2014). Incubations for proteomic SIP were collected from this mesocosm experiment; four incubations were set up, corresponding to two sediment samples (#3731 and #5133), each receiving two experimental treatments (Figure 1). 60 mL of the #3731 slurry was transferred into each of two sterile 120 mL serum bottles in an anaerobic chamber and supplemented with 1 mM final concentration of either 14NH4Cl or 15NH4Cl (99% enrichment, Cambridge Isotope Laboratories, Inc., Andover, MA). The vials were sealed with butyl rubber stoppers, flushed for 5 min each with Ar and then CH4, and overpressured with CH4 to a final pressure of 2 bars. #5133 14 and experiments were set up in a similar manner, though the size of the serum bottles was 35 mL and initial slurry volumes were 20 mL. All samples were kept at 4°C in the dark during the incubation period.
Figure 1. A schematic diagram of the proteomic SIP sampling approach used in this study. Sediment was collected from active seep areas at Hydrate Ridge North for incubation #5133 and Hydrate Ridge South for #3731 (see text for full sampling details). Sediment from both sample sites was used as inoculum in parallel methane-infused incubations, with one bottle from each set receiving 1 mM 14NH4Cl and the other receiving 1 mM 15NH4Cl. Both #3731 incubations were sampled for proteomic SIP after 17 and 326 days; both #5133 incubations were sampled for proteomic SIP after 160 days.
Aggregate counts were conducted at incubation start and end points, and sulfide concentrations were determined at four time points. To quantify the number of cell aggregates, incubation sediment was subjected to a percoll density separation performed according to Orphan et al. (2009). Sediment slurry (~60 μL) was mixed with 1290 μL of TE (pH = 9) and heated for 3 min at 60°C. The sample tubes were put on ice for 5 min; 4.5 μL of 30% H2O2 was then added (to deactivate peroxidase) and incubated at room temperature for 10 min. Tubes were put back on ice, and 150 μL sodium pyrophosphate was introduced. Sample mixtures were sonicated on ice 3 times at 8 W (Branson sonifier 150) for 10 s each time and applied to the density gradient vials. The vials were then centrifuged at 4800 rpm for 15 min at 4°C. The supernatant was filtered through both 3 μm and 0.22 μm polycarbonate filters, which were then rinsed with 2 mL PBS and dehydrated with 2 mL 1:1 ethanol:PBS. Aggregate counts were performed under epifluorescence microscopy (60x objective lens, Olympus BX51 microscope, Olympus, Melville, NY) using DAPI; the corresponding morphologies have been shown via FISH analysis to consist dominantly of ANME-SRB aggregates active in anaerobic methane oxidation (Orphan et al., 2002, 2009; Marlow et al., 2014b). Sulfide concentrations were determined using the Cline Assay (Cline, 1969) as described in Dekas et al. (2009).
16S rRNA Gene Tag Sequencing
DNA from 14N incubation samples #3731 T0d, #3731 T326d, and #5133 T170d was extracted from ~0.5 g of sediment using the UltraClean Soil DNA isolation kit (Mo Bio Laboratories, Carlsbad, CA) and microbial community composition was profiled using the 515f (5′-GTGYCAGCMGCCGCGGTAA-3′) and 806r (5′-GGACTACNVGGGTWTCTAAT-3′) primer pair (Caporaso et al., 2011). PCR barcoding reactions were conducted in house at Caltech and amplicons were outsourced to Laragen Inc. (Laragen Inc., Culver City, CA) for sequencing on the Illumina MiSeq platform. Resulting sequence data consisted of 250 bp paired end reads, which were assembled into a single contig and analyzed using Qiime 1.8 (Caporaso et al., 2010). Sequences were trimmed to a minimum quality value of 30, and taxonomy was assigned to each OTU representative based on the SILVA 115 rRNA database uclust (Edgar et al., 2011) using a sequence similarity threshold of 99%. The full protocol is contained in a GNU makefile (Supplemental Data Sheet 1); the makefile must be in the same directory as the raw sequencing reads, and qiime 1.8 and usearch must be installed. The complete protocol can be run by typing “./orphanlab_itag_protocol.mk all” into the command prompt.
Metagenome Database Construction and Annotation
The metagenomic database against which tandem mass spectra were referenced to identify proteins was constructed from several sources (Table 1). DNA was extracted using the UltraClean Soil DNA isolation kit (Mo Bio Laboratories, Carlsbad, CA) from ~0.5 g of bulk methane seep sediment (before distribution to the proteomic SIP incubations) and from magneto-FISH capture experiments from Hydrate Ridge (Pernthaler et al., 2008; Dekas et al., 2009; Glass et al., 2014; Table 1). Bulk sediment metagenomes from samples #3730, #5133, and #5579 were sequenced using the Illumina HiSeq platform; four magneto-FISH metagenomes (targeting ANME-2c) were sequenced using 454 GS-20 (BC02, Pernthaler et al., 2008), 454 GS-FLX (BC03), and 454 GS-Titanium (BC04 and BC05). Metagenomes derived from Illumina HiSeq data were processed as follows. Raw sequencing reads were quality trimmed (using a threshold of 0.05) and assembled using CLCBio genomics workbench 6.0. Sequencing reads from each metagenome were assembled individually and two assemblies were performed by combining sequencing reads from the three datasets. For each assembly, the raw sequencing reads from all datasets were aligned to contigs >500 bp with BWA 0.7.5a (Li, 2013) using the “mem” algorithm and converted into bam files using samtools 0.1.19 (Li et al., 2009). These bam files were used as input for GroopM 0.2 (Imelfort et al., 2014), which generated genome bins. The “refine” command in GroopM was used to improve the genome bins manually by removing contigs with unusual coverage or tetramer frequencies. The completeness and contamination of genome bins was determined using CheckM 0.9.4 (Parks et al., 2014).
Bead capture metagenomes derived from 454 sequencing platforms generated ~25–145 Mb of sequence data per sample. Sequencing reads from all bead capture metagenomes were combined and assembled using Newbler. Open reading frames (ORFs) were called using prodigal 2.5 (Hyatt et al., 2010) using the metagenome mode (command line ‘prodigal –i < FILE> -p meta –a < translation.faa> -o < output.gbk>’). Only ORFs that were considered complete (contained both a start and stop codon) were considered for inclusion to the metagenomic database. This data set was supplemented with fosmid sequence from seep sediment and physiologically relevant genomes from cultured organisms (Supplemental Data Sheet 2). All ORFs were clustered at 100% identity using uclust (Edgar et al., 2011) to remove duplicates.
Functional annotations were derived from previously determined NCBI annotations (for cultured organism genomes) or the KEGG Automatic Annotation Server (KAAS; Moriya et al., 2007) for genes from the metagenomes. For cultured organisms, phylogenic affiliations of gene products were incorporated from NCBI annotations. For the metagenomes, the phylogeny of genome bins with low (< 10%) contamination was assessed using the CheckM parsimonious genome tree, and the proteins derived from that genome bin were all given the same taxonomy. For all other proteins derived from the metagenomes, individual genes were compared to the NCBI nr database using blastp 2.2.29+ (Camacho et al., 2009) and taxonomy was determined using the lowest common ancestor algorithm implemented by Blast2lca 0.600 (https://github.com/emepyc/Blast2lca). #5133 15N enriched proteins that were annotated as “hypothetical proteins” by KAAS (n = 110) were subjected to a secondary search process. Corresponding ORFs were blasted against UniRef 90 (Suzek et al., 2007) using the UniProt webserver (UniProt Consortium, 2008) and an e-value cutoff of 10−30. Annotations derived from this process were incorporated into the data analysis, while ORFs that remained “hypothetical” or “uncharacterized” were searched by hand against the Interpro database (Mitchell et al., 2015) for putative functional domains. The subset of these proteins that was also present in all six proteomes is discussed below.
Proteins involved in the reverse methanogenesis pathway were resolved to finer phylogenetic resolution to determine their ANME subtype affiliation. Proteins from the metagenomes were combined with proteins downloaded from the KEGG database and aligned using muscle 3.8.31 (Edgar, 2004). A phylogenetic tree was built from the alignment using FastTree 2.1.7 (Price et al., 2010) and visualized with Archaeopteryx 0.9901 (https://sites.google.com/site/cmzmasek/home/software/archaeopteryx). ANME subtypes were determined by identifying monophyletic groups with known ANME-1, ANME-2a, ANME-2b, and ANME-2c sequences and cross validating between genes from the same contigs to resolve inconsistencies in phylogenetic placement. Protein alignments were visualized in Jalview 2.8.2 (Waterhouse et al., 2009).
Metaproteomic Sample Processing
The proteomic SIP sampling approach is shown in Figure 1. Sample #3731 incubations were sampled after 17 days (T17d) and 326 days (T326d) of incubation with methane; #5133 was sampled after 160 days (T160d). These sampling times were determined based on the multi-month doubling times of ANME-SRB aggregates (Orphan et al., 2009), and the lack of significant 15N enrichment in the T17d protein sample that necessitated the extended incubation timescale. At the designated sampling times, 20 mL (#3731) or 5 mL (#5133, small volume due to sample limitation) of sediment slurry was removed and immediately flash frozen with liquid nitrogen to stop biological reactions and preserve proteins representative of the incubation conditions. Proteins were extracted and digested using an SDS-TCA-based protocol as outlined in Chourey et al. (2010). To summarize, frozen sediment slurry was thawed, and 10 g (#3731 T17d and T326d) or 5 g (#5133) of wet sediment was suspended in 20 mL of lysis buffer [5% SDS, 50 mM Tris-HCl, pH 8.5, 0.15 M NaCl, 0.1 mM EDTA, 1 mM MgCl2, 75 mM dithiothretiol (DTT)]. The mixture was vortexed and immersed in a boiling water bath for 20 min, allowed to cool for 5 min, and transferred to a sterile 50 mL Falcon tube. A 10-min centrifugation at 2500 × g separated soil particles and cellular debris from the whole cell lysate, which was retained in the supernatant. The supernatant was pipetted into sterile 2-mL Eppendorf tubes, amended with 100% trichloroacetic acid (TCA) to a final concentration of 25%, vortexed, and stored overnight at –20°C to precipitate the extracted proteins. Protein pellet was collected by centrifugation at 21,000 g for 20 min, and the supernatant was discarded. The protein pellet was washed with chilled acetone followed by centrifugation at 21,000 g for 5 min. After three acetone washes, the protein pellet was air dried and solubilized in Guanidine buffer [6 M guanidine HCl, 10 mM DTT in Tris buffer (50 mM Tris with 10 mM CaCl2, pH 7.8)] followed by an overnight incubation at 37°C. The mixture was diluted six-fold with Tris CaCl2 buffer to dilute the guanidine concentration to 1 M and avoid trypsin deactivation. Trypsin was added at a predetermined concentration of 40 μg trypsin/1–3 mg total protein in two separate doses of 20 μg each, separated by several hours of incubation at 37°C; the second addition was allowed to incubate overnight. DTT was again added to a final concentration of 10 mM and allowed to incubate for 45 min at room temperature to linearize the peptides. The peptide solution was de-salted using a C18+ Sep-Pak solid phase extraction cartridges (Waters, Milford, MA) solvent exchanged into acidified water (HPLC grade water, 0.1% formic acid) using a Savant SpeedVac (ThermoFisher Scientific, Waltham, MA), filtered through an Ultrafree-MC centrifugal filter device (Millipore, Billerica, MA), and stored at −80°C until two-dimensional LC−MS/MS analysis.
Duplicate runs of LC-MS/MS analysis were conducted for each sample as described previously (Pan et al., 2011; Justice et al., 2014). Peptides (~75 μg) were pressure-loaded onto a split-phase Reverse Phase (RP)-SCX column, which was aligned with an in house RP-packed PicoFrit tip (New Objective, Woburn, MA). Peptides were eluted and chromatographically separated via a 24-h Multi-Dimensional Protein Identification Technology (MuDPIT), using an 12-step salt gradient using liquid chromatography (LC) as described earlier (Sharma et al., 2012). Peptide separation, fragmentation, and measurements were carried out using a range of Orbitrap instruments in keeping with the rapid development of the field's best available technology: #3731 14N T17d and #3731 15N T17d were analyzed on an LTQ Orbitrap Velos, #3731 14N T326d on an LTQ Velos, #3731 15N T326d on an LTQ Orbitrap Velos Pro, and #5133 14N T160d and #5133 15N T160d on an LTQ Orbitrap Elite. The mass spectrometers (Thermo Fisher Scientific, Germany) were coupled to the Ultimate 3000 HPLC system (Dionex, USA) and operated in data dependent mode using the Thermo Xcalibur software V2.1.0. Each full scan was followed by individual MS/MS scans of the 20 most abundant parent ions, which were fragmented via collision-activated dissociation (CID) using 35% collision energy, a mass exclusion width of 0.2 m/z, and a dynamic exclusion duration of 60 s.
Protein Search Parameters
MS/MS data was searched using the Sipros/ProRata 3.0 program (Wang et al., 2013); precursor mass tolerance and fragment ion mass tolerance were 0.05 Da and 0.02 Da for the data generated on Orbitrap instruments and 3 Da and 0.5 Da for those generated on LTQ Velos, respectively. The metaproteomic database was predicted from the metagenomic database (Table 1) and concatenated with its reverse sequence to estimate its false positive detection rate (FDR; Elias and Gygi, 2007). The search results were filtered to achieve a 1% FDR at the peptide level, and proteins were inferred from the identified peptides using parsimony rules (Nesvizhskii and Aebersold, 2005). The most stringent settings, which have been used in several proteomics studies (Denef et al., 2009; Verberkmoes et al., 2009; Liu et al., 2012), are a 1% FDR, 2 total peptides, and 1 unique peptide (expressed with the notation [1%, 2TP, 1UP]). However, the [1%, 1TP, 1UP] approach has been promoted on statistical grounds (Gupta and Pevzner, 2009) and used in several environmental and culture-based proteomics studies (Ram et al., 2005; Dong et al., 2010; Webb et al., 2013; Hebert et al., 2014). With the acknowledgement that less stringent protein identification parameters result in lower confidence of a protein's presence, we present both [1%, 1TP, 1UP] and [1%, 2TP, 1UP] results, clearly noting the provenance of each identification.
The following PTMs were dynamically searched in McrA analysis: hydroxylation of proline and lysine; monomethylation of arginine, lysine, aspartic acid, and glutamic acid; dimethylation of arginine and lysine; trimethylation of lysine; phosphorylation of serine, threonine, tyrosine, histidine and aspartic acid; acetylation of lysine; S-nitrosylation of cysteine; nitration of tyrosine; and methylthiolation of aspartic acid. The raw search results from both the PTM search and the broader search against the full database were filtered together to achieve a 1% FDR at the peptide level, allowing the top peptide of a MS/MS spectrum from the PTM search to compete with the top peptide of the same spectrum from the regular search. PTMs were localized based on the DeltaP score as described in Li et al. (2014).
SIP Enrichment Threshold and Analysis
A smaller database constructed of all ORFs identified during the initial non-SIP 1TP, 1UP analyses (22,312 ORFs) and their reverse decoy sequences was used for SIP searches and filtering due to computational limitations. Because protein identification is predicated on the detection of constituent peptides, and because 1TP, 1UP identifications were considered in the analysis, designations of “enriched” proteins (defined as proteins that contain 15N with statistical confidence) must be made at the peptide level.
Procedural steps for determining 15N atom% are fully described elsewhere (Pan et al., 2011; Fischer et al., 2013; Justice et al., 2014). Briefly, for each peptide sequence from in silico digestion of target/reverse proteins in the database, the most abundant mass and isotopic distribution of b and y ions were predicted at a given 15N atom% ranging from 0 to 100% (1% intervals), assuming uniform 15N enrichment. For each experimental MS/MS spectrum, candidate peptides were found if their most abundant masses at a given 15N atom% were within seven mass windows centered at the parent mass of this experimental spectrum. Each candidate peptide was scored based on the mass accuracy and goodness of fit between calculated b and y ions and measured b and y ions. Those candidate peptides with the highest score from the entire MS run were filtered using the target/decoy approach (Elias and Gygi, 2007) to achieve 1% FDR at the peptide level.
To set the 15N value above which a peptide would be declared “enriched”, multiple factors were considered. First, peptide 15N distributions from each sample's 14 incubation were determined, and the value below which 99% of peptides fell set a 1% FDR for the declaration of enriched peptides. This value was 3% 15N in #3731 T17d, 4% in #3731 T326d, and 5% in #5133 T160d. Alternatively, we can consider the probability of unenriched peptides containing a range of 15N values using the binomial distribution. Using the mean enrichment value of the corresponding 14 incubation as the baseline 15N frequency (0.45% for #3731 14N T17d, 0.35% for #3731 14N T326d, and 0.49% for #5133 14N T160d), expected frequencies of X-length peptides containing Y% 15N can be calculated. When the probability of different 15N atom % values was plotted based on the binomial distribution and a best fit curve was calculated (R2 = 0.99924), the 15N enrichment value corresponding to a 1% likelihood of natural occurrence (i.e., the 15N value corresponding to a 99% confidence that it in fact was not enriched), was 0.1% for the shortest peptide (six amino acids). In a 53-residue peptide (the longest detected in this study), 70 N atoms are predicted; the likelihood of one 15N atom incorporation is much higher (24.6%) but the atom % 15N this represents is substantially lower (1.43%). There is a 4.2% chance of a 2.86% 15N enrichment, and a 0.47% chance of a 4.29% enrichment.
Of these statistical analyses, the most conservative 99% confidence interval surrounding an “enriched” designation is 5% 15N (derived from the 15N distribution in peptides recovered from the #5133 14N sample). This value therefore demarcates enriched proteins in this study's analysis; the same enrichment level has been used as a threshold in other SIP studies (Pan et al., 2011; Justice et al., 2012).
Bulk protein extract from sample #3731 and cell aggregates from sample #5133 were analyzed with the Cameca nanoSIMS 50 L. Targets were mapped on the indium tin oxide-coated glass squares using both epifluorescence and transmitted light microscopy with a DeltaVision RT (Applied Precision Inc., WA) according to Dekas and Orphan (2011). For our analyses, a 2.5 pA cesium beam with a spot size of 100–200 nm was used. The beam rastered over areas dependent on the size of the aggregates selected (generally between 5 and 20 μm) with a 256 × 256 pixel resolution, a dwell time of 5–10 ms per pixel, and a resolving power of ~5000. 15N12C− and 14N12C− ions were collected for multiple frames, each lasting ~30 min. NanoSIMS images were processed with L'image (developed by L. Nittler, Carnegie Institution of Washington, Washington D.C.). Each set of frames was corrected for drift and detector dead time. Clostridia spores with known 13C/12C and 15N/14N ratios were used as a standard and measured before and after sample analysis to correct for instrumental mass fractionation and drift.
Sequence data can be found at the NCBI bioproject database, accession number PRJNA290197.
Results and Discussion
Activity and Community Composition in Seep Sediment Microcosms
Samples of methane seep sediment from Hydrate Ridge (hereafter referred to as #3731 and #5133) were each allocated into two anaerobic microcosm incubations overpressured with methane and amended with either 14 (unlabeled treatment) or 15 (labeled treatment; Figure 1). Cell aggregate abundance and sulfide concentrations, both of which correlate positively with sulfate-based AOM activity (Iversen and Jorgensen, 1985), were monitored throughout the incubation period. Both #3731 and #5133 were collected from actively seeping areas, yet the #3731 incubation, recovered from the lower 6–12 cm sediment horizon, demonstrated indications of lower metabolic activity (Table 2) and 15N incorporation (see below). Potential reasons for this discrepancy include differing in situ methane concentrations despite general markers of activity, localized heterogeneity of methanotrophic constituents (Table 2), and sediment horizon depth-related variability in AOM activity.
Microbial communities were analyzed using 16S rDNA iTAG sequencing to compare the changes in diversity between the two time points of #3731 and between sediment #3731 and #5133 (Supplemental Data Sheet 3a). Across all three samples, the most abundant representatives were Deltaproteobacteria from the uncultivated SEEP-SRB1 (accounting for 7.2% of sequences on average, ± 1.3% SD) and SEEP-SRB2 (5.1 ± 1.7% SD) clades (Schreiber et al., 2010), Sulfurovum (6.7 ± 2.8% SD), and ANME-1a (6 ± 2.6% SD) and ANME-1b (5 ± 1.3% SD). ANME-2 subclades accounted for an average of 3.2% of sequences (±0.56% SD); however, the primer set used has been shown to preferentially detect ANME-1 over ANME-2 (Case et al., 2015), and independent FISH-based microscopy indicates a prevalence of ANME-2 consortia in the initial inoculum of both samples (data not shown). In all three samples, ANME sequences represented a substantial proportion of the archaeal relative abundance (67.8% on average, ± 6.3% SD). AOM-related enzymes of unspecified lineage putatively come from ANME organisms given that methanogenic lineages comprised < 0.1% of the archaeal community across all samples at the order level, as determined by 16S rRNA gene tag sequencing (Supplemental Data Sheet 3b).
The metaproteomic dataset produced here, using optimized extraction techniques and Orbitrap mass spectrometers, represents the deepest proteomic measurement of methane seep sediments to date and rivals most other environmental proteomics studies. Analyses at the 1% FDR level yielded 5664 unique proteins across all samples under the 1TP, 1UP condition, and 3495 unique proteins under the 2TP, 1UP condition (see Table 3 for quantification details). Compared with previous environmental studies (Supplemental Presentation 1; Table S1), the large number of proteins and selective incorporation of stable isotope label demonstrated here indicates that increased microbial diversity, low overall activity levels, and challenging physicochemical conditions are not insurmountable challenges in SIP-resolved metaproteomic studies. Nonetheless, the quantitative recovery, detection, and characterization of all proteins is not technically feasible because of protein adsorption to silicate minerals (Ding and Henrichs, 2002), complexation (Nguyen and Harvey, 2001) as well as the difficulty of solubilizing membrane-bound constituents (Tan et al., 2008).
A core set of 456 (1TP, 1UP) and 283 (2TP, 1UP) proteins were detected in all six samples; of these, the five most common categories were “proteins of unknown function” with best matches to hypothetical genes (comprising 31.1% of protein identifications averaged between the search conditions), adenylyl sulfate reductase subunit A (AprA, 7.4%), McrB (5.7%), McrG (4.9%), and ATPase subunit A (4.6%). Phylogenetic assignments of the cumulative metaproteome demonstrated a higher relative abundance of Bacteria-classified proteins (56.9% 1TP, 1UP; 58.7% 2TP, 1UP) than those assigned to Archaea (37.2% 1TP, 1UP; 35.2% 2TP, 1UP).
15N Enrichment of Proteins Expressed in AOM Sediment Incubations
To provide additional context for 15N incorporation-based activity levels, the isotopic composition of bulk protein extract and intact ANME-SRB consortia was measured by nanoscale Secondary Ion Mass Spectrometry (nanoSIMS). NanoSIMS analysis of #3731 sediment bulk protein extracts from the 15 amended incubation revealed a slow-growing community with a low but measurable enrichment at the early (0.45 atom % enrichment for 15N T17d) and late (0.53 atom % at 15N T326d) time points compared to the 14 treated controls (0.38 atom % at both time points; Table 4). In silico identification of enriched peptides from sample #3731 also exhibited a small but detectable isotopic enrichment of the protein pool during the 326-day incubation. The 15N atom % enrichment of the entire #3731 15N peptide pool based on Sipros analysis was 0.4% at T17d, increasing to 0.8% at T326d (while corresponding values for 14N control incubations were 0.5 and 0.4%, respectively).
Peptides recovered from sample #5133 15N after 160 days of incubation with 15 demonstrated significantly more uptake of the 15N label, with Sipros-calculated isotopic enrichment of 9.5 atom % 15N, while peptides from the corresponding 14 control exhibited a mean enrichment of 0.5%. NanoSIMS analysis of ANME-SRB consortia from the same sample showed substantially higher 15N incorporation levels, reaching 23.8% after 64 days (Supplemental Presentation 1; Figure S1). The broader peptide-based enrichment value may provide a more representative constraint on community-wide anabolic rates because ANME-SRB consortia constitute a particularly active subset of the community (Knittel and Boetius, 2009). Peptides that uniquely mapped to ANME archaea had a mean 15N value of 13.8%, while those attributed to the dominant deltaproteobacterial sulfate-reducing group (Desulfobacterales) had a mean enrichment value of 10.1%. Proteins previously linked to the reverse methanogenesis pathway had a mean enrichment of 18.3%, whereas proteins associated with dissimilatory sulfate reduction exhibited a mean enrichment of 9.3%. These values suggest that ANME lineages responded to the methane incubation conditions with heightened biosynthetic activity compared to SRB, and that a higher proportion of reverse methanogenesis enzymes were newly synthesized compared to other identified ANME proteins.
The observation of enhanced 15N incorporation in the nanoSIMS analysis of individual intact AOM consortia relative to the calculated 15N values from bulk ANME and SRB-associated proteins recovered from the same incubation is not fully understood. Potential explanations include the incomplete detection of ANME and deltaproteobacterial proteins, enrichment discrepancies among organisms or different cellular nitrogen sources, or variation in the turnover of proteins and other cellular nitrogen stores (Reitzer, 2003). Membrane proteins are believed to comprise 20–30% of prokaryotic proteomes based on genomic analysis (Stevens and Arkin, 2000), yet solubilizing membranes and accessing these proteins is a technical challenge (Trötschel and Poetsch, 2015); preferential incorporation of 15N into membrane proteins, such as receptors, transporters, or oxidoreductases may be partially responsible for the 15N unaccounted for in our SIP metaproteomic analysis.
Peptide enrichment values exhibited a marked bimodal distribution, with an unenriched peak and a 15N enriched peak at 49% for the #5133 15N T160d sample and 64% for the #3731 15N T326d sample (Figure 2). The presence of two distinct peptide pools in the data presented here implies that most of the detected proteins (made up of peptides with 15N values below 5%) were likely synthesized before the addition of 15 and were not fully degraded, while the remainder was newly synthesized. The lack of peptides with intermediate enrichment values argues against multiple protein precursor pools (Justice et al., 2014) or substantial heterotrophic turnover of the system's primary producers during the course of the 160-day and 326-day incubations.
Figure 2. 15N enrichment distributions for all peptides identified in (A) #5133 15N T160d and (B) #3731 15N T326d incubations. In both cases, the precursor pool of was 496 μM and 1 mM of 15 was added at the beginning of the incubation period.
Functional and Phylogenetic Distributions of Enriched Proteins
Of the proteins identified from the #5133 15N sample, 32% (n = 377) were isotopically enriched using the conservative (2TP, 1UP) criterion. To determine which functional traits or phylogenetic affiliations may be disproportionately represented in either the enriched or unenriched protein pools, annotated proteins were partitioned by enzyme type or phylogenetic assignment of the corresponding ORF. The relative abundance of these categories was calculated for both the enriched and unenriched proteins, and the variance contributed by each bin was determined (calculated by the square of the difference in relative abundances).
Eighteen functional protein categories accounted for 79.8% of the variance between enriched and unenriched proteins, with each contributing at least 1% of the variance (Figure 3A). The largest contributor, whose distribution produced 24.8% of the observed variance, was RNA polymerase (RNAP), which was more prominent in the unenriched protein pool. By examining the phylogenetic assignments of RNAP orthologs, the enzyme serves as a representative signal of organism-wide regulation. The detection of remnant, unenriched RNAP proteins reveals their resistance to degradation over the course of the 160-day experiment and offers the possibility that all transcription needs were met by pre-existing, long-lived copies. The relative lack of newly produced copies could indicate that the incubation conditions were deleterious to active transcription and, perhaps, survival of some microorganisms; this scenario is consistent with the detection of few enriched proteins from similar phylogenetic groupings, suggesting a broader shutdown of transcriptional machinery. In this context, the lineages responding most negatively to the anaerobic incubation conditions include Methylococcales, Sulfurovum, Chromatiales, Alteromonadales, and Pseudomonas (Figure 3A), indicating a shift away from putative oxygen- or nitrate-based metabolisms common within the oxidizing conditions of upper seep sediment layers.
Figure 3. (A) Protein functional groups accounting for at least 1% of variance between the enriched and unenriched #5133 15N T160d protein pools. Positive values reveal protein types that are more abundant in the enriched fraction; negative values indicate those more prevalent in the unenriched fraction. Inset table shows the phylogenetic affiliations of unenriched RNAP and enriched GroEL proteins. (B) Phylogenetic associations accounting for at least 1% of variance between the enriched and unenriched #5133 15N T160d protein pools. Positive values reveal phylogenetic assignments whose protein products are more abundant in the enriched fraction; negative values indicate those more prevalent in the unenriched fraction. Phylogenetic assignments were made at the family level; higher level assignments are provided if no family-level specificity was available, and genus-level assignments are provided if no other genera were observed in the same family. Archaeal Gzfosmids are ANME-affiliated sequences reported by Hallam et al. (2004).
Of the predominantly enriched protein types in #5133, chaperonin GroEL accounts for the most variance (10.9%) between the enriched and unenriched pools, followed by F-type ATPase subunits B (8.4%), and A (6.1%). GroEL is frequently detected in environmental metaproteomic studies (Benndorf et al., 2007; Sowell et al., 2009; Verberkmoes et al., 2009) and has been demonstrated to prevent amorphous protein aggregations (Hartl et al., 2011), facilitate proper folding of a wide variety of enzymes (Kerner et al., 2005), and mediate stress response (Kim et al., 2013); changing pressure or energetic regimes in our incubations may have prompted such a response. F-type ATPase and pyrophosphatase PpaC (1.4%) are both involved in ATP generation, concordant with previous observations of heightened energy production under conditions of physiological stress (Alexandre et al., 2001) such as those potentially experienced by slow-growing methanotrophic consortia (Valentine, 2007). Several high-variance proteins are involved in reverse methanogenesis and sulfate reduction pathways. AprA (5.1%), McrB (3.0%), McrG (2.8%), sulfite reductase subunit B (DsrB; 2.4%), McrA (2.1%), MtrA (1.7%), AprB (1.4%), and sulfate adenylyltransferase (Sat, 1.3%) were all more prevalent in the 15N enriched fraction. This observation is consistent with active sulfate-coupled AOM as demonstrated by sulfide production during the incubation as well as quantitative increases in and growth of ANME-SRB consortia (Table 2, Supplemental Presentation 1, Figure S1).
Of the high-variance proteins that were detected at higher levels within the unenriched fraction, methyl-accepting chemotaxis proteins (accounting for 1.6% of the variance) may be used in biofilm-linked cell signaling or orientation toward favorable environmental conditions (Morgan et al., 2006). Most of the observed chemotaxis proteins exhibit greatest homology with Gammaproteobacteria, and their near-absence among 15N-enriched proteins suggests decreased need for mobility-linked sensing and/or less communication among these constituents. Glutamine synthetase (3.0%) catalyzes the formation of glutamine from glutamate and ammonia, and its minimal presence among enriched proteins could indicate the accessibility of free glutamine (and likely other amino acids) from biomass degraded during the incubation, or a heightened reaction rate—facilitated by increased concentrations—that obviated new enzyme production. Predominantly unenriched proteins may reveal functions of lower importance under incubation conditions, while reflecting remnant capabilities that highlight the importance of proteins' ability to outlive upstream markers of activity (e.g., transcripts (Moran et al., 2013).
Fifteen phylogenetic bins each contributed at least 1% of the variance between the relative abundances of enriched and unenriched proteins recovered from #5133 15N, representing 92.5% of the total variance (Figure 3B). Among the more notable differences was the prevalence of enriched proteins linked to methane-metabolizing lineages (accounting for 32.3% of the variance) and Deltaproteobacteria (27.7%), particularly orders with known sulfate reducers (18.3%) expressing proteins from carbon fixation and sulfate reduction pathways, GroEL, and ribosomal proteins. Proteins attributed to Epsilonproteobacteria (20.8%) and especially Sulfurovum RNAP, Apr, and ATPase (15.3%) were markedly less prevalent in the enriched pool. The dearth of enriched epsilonproteobacterial protein products is consistent with these organisms' proposed susceptibility to changing conditions (Alain et al., 2004; Sievert et al., 2007; Toner et al., 2012) as well the potentially unfavorable conditions for supporting active oxygen- or nitrate-coupled sulfide oxidation. The high relative abundance of Sulfurovum in the #5133 tag sequencing data (Supplemental Data Sheet 3a) highlights the utility of SIP proteomics: the abundance and persistence of 16S rRNA genes and proteins may not always be reflective of biosynthetically active participants in an environment.
Finally, of the proteins that remained unannotated (i.e., proteins of unknown function using the 2TP, 1UP criterion), 26 were enriched in the #5133 15N T160d sample and were present in all six proteomes, suggesting they are important in seep sediment habitats and are actively synthesized under conditions supporting AOM. Approximately half of the unknown proteins were attributed to uncultured Archaea while most of the remainder was associated with sulfate-reducing Deltaproteobacteria (Table 5). Up to 40% of expressed bacterial proteins currently lack functional classification, and there is a growing need for further categorization and characterization of these proteins (Galperin and Koonin, 2004). These 26 consistently detected, uncharacterized proteins represent notable targets for functional biochemical study (Guengerich and Cheng, 2011; Goodacre et al., 2014).
Table 5. UniProt results for the 26 proteins that were detected in all six proteomes, enriched in #5133 15N T160d, and were unannotated via KAAS and UniProt.
The anaerobic oxidation of methane (AOM) is a dominant metabolism at Hydrate Ridge methane seeps. Many distinct orthologs of all enzymes involved in the reverse methanogenesis pathway were recovered from incubated seep sediment (Figure 4; phylogenetic affiliations provided in Supplemental Data Sheet 4a), signifying unprecedented completeness and diversity of pathway protein detection. Methyl-coenzyme M reductase (Mcr) is one of the most abundant enzymes in ANME, accounting for 7% of total extractable proteins from the ANME-1 rich Black Sea mats (Krüger et al., 2003). Across all six of our proteomic data sets, 10.4% of peptides were attributable to Mcr; 34.6% of 15N enriched peptides in sample #5133 were associated with this critical enzyme. Considerable diversity of detected Mcr proteins was observed, with approximately 20 orthologs of each Mcr subunit identified across all six samples, representing 47% of the mcr protein sequences within our seep metagenome database. The #5133 microcosm thus appeared to stably maintain a genetically diverse population of active ANME archaea over the course of the extended incubation. Of the expressed alpha, beta, and gamma subunits of MCR, one-third was affiliated with ANME-1 (31% of detected orthologs, see Supplemental Data Sheet 4a), with smaller contribution of ANME-2c (19%), Methanosarcinaceae (14%), and ANME-2a derived (12%) orthologs. The three non-catalytic Mcr components (McrA2, McrC, McrD) were also identified, but these putative (methanogenic) activation domains (Prakash et al., 2014) were notably less prevalent in the metaproteome. This result may indicate McrA2, McrC, and McrD re-use as activating components, or their decreased relevance under methanotrophic conditions, potentially reflective of a methane-oxidizing enzymatic mechanism that is distinct from a strict reversal of methane formation.
Figure 4. Metaproteomic data for enzymes involved in the reverse methanogenesis pathway. Filled boxes indicate an enzyme constituent's presence in the associated sample's proteome. Due to space constraints, Hdr detections are not shown, but relevant information is provided as Protein orthologsdetected(enrichedorthologs)∕orthologsinmetagenomicdatabase.
Several previous ANME-targeted meta-omics efforts have sought to reconstruct the reverse methanogenesis pathway (Table 6; e.g., Hallam et al., 2004; Meyerdierks et al., 2010; Stokke et al., 2012; Haroon et al., 2013; Wang et al., 2014), but the Mer protein in particular has proven challenging to identify (Hallam et al., 2004; Meyerdierks et al., 2010; Stokke et al., 2012). In our proteomic experiments, seven Mer orthologs were identified (six of which were 15N-enriched), exhibiting best homology to ANME-2a, M. burtonii, and other Methanosarcinales lineages. Our detection of 15N-enriched Mer orthologs under incubation conditions favoring AOM is consistent with a reverse methanogenesis metabolism. It is possible that ANME-1 methane oxidation proceeds via the mer bypass (Meyerdierks et al., 2010), whereby activated methane is converted first to methyl-containing intermediates and then to formaldehyde, which is processed to methylene-H4MPT by a formaldehyde-activating enzyme (Fae) and hexulose-6-phosphate (Hps). The viability of the methanogenic version of this pathway has been demonstrated with an M. barkeri mer deletion mutant (Welander and Metcalf, 2008). An ANME-1 mer copy has not yet been identified, and 12 orthologs of the Fae/Hps fusion protein were detected across all samples (1TP, 1UP).
Table 6. Central reverse methanogenesis enzymes identified in metagenomic, metatranscriptomic, and metaproteomic analyses.
Green, identified; orange, not addressed; red, not identified. ANME proteins newly detected in seep environments during this study are highlighted in bold text. (McrABG subunits associated with ANME-2c derived fosmids were also detected; no other ANME-2c genes related to reverse methanogenesis were in the metagenomic database.)
Building upon the 356 proteins detected from seep sediments in the Nyegga area (Norwegian Sea; Stokke et al., 2012), additional subunits of enzymes putatively involved in the reverse methanogenesis pathway were newly detected in our study, including McrC and McrD (Table 6). Of the components not found in the proteome, some of those present in the co-located metagenomes (MtrBCDEF, HdrDE) are membrane-bound proteins (Thauer et al., 2008) and thus are expected to be more difficult to quantitatively recover (Trötschel and Poetsch, 2015). Given the relatively poor genomic coverage of the ANME-2 lineages, combined with the observed dominance of ANME sequences relative to methanogens in our iTAG diversity surveys (Supplemental Data Sheet 3b), it is likely that most of the reverse methanogenesis pathway orthologs of poor homology derive from ANME representatives.
McrA Post Translational Modifications
Mcr is a critical enzyme in the reverse methanogenesis pathway, activating methane with a disulfide bond and initiating the multi-step anaerobic oxidation of methane (Krüger et al., 2003; Shima and Thauer, 2005). To better understand the range of variability found between McrA orthologs and potentially distinguish between methanogenic and methanotrophic forms of the enzyme, we conducted a computational search for selected post-translational modifications (PTMs; Camacho et al., 2009). MS-based PTM analysis is a developing field, and it is difficult to draw functional conclusions from detected modifications; nonetheless, the approach has been used to corroborate modifications (Selmer et al., 2000) and as a discovery tool supplemented by subsequent crystallographic confirmation (Moellering and Cravatt, 2013). Spectra from sample #5133 14N T160d were used due to its enhanced metabolic activity relative to #3731 (Table 2) and unlabeled peptides, which ensured the accuracy of PTM mass windows. PTMs were observed among 10 of the 18 McrA orthologs, 13 of 19 McrB orthologs, and 7 of 17 McrG orthologs; none were identified on five McrA2, one McrC, and three McrD orthologs. Twenty-eight distinct modification events of McrA were observed, including mono-, di-, and tri-methylations, hydroxylation, beta-methylthiolation, and acetylation (Figure 5). PTMs were detected with high frequency (occurring an average of once every 84 amino acids among ANME Mcr subunits), but only 11% were affiliated with multiple orthologs. This observation, in addition to their positioning away from active site residues, suggests a dominance of PTMs responsive to environmental variation (Olsen et al., 2010) that may be involved in enzyme regulation (Olsen and Mann, 2013).
Figure 5. A phylogenetic and PTM-based analysis of the McrA subunit. The left panel shows a Muscle-aligned phylogenetic tree of all McrA sequences in the metagenomic database, as well as additional cultured organisms to account for all seven methanogenic orders and other ANME representatives. The metaproteomic detection panel shows the presence/absence bars from Figure 4, positioned next to their corresponding ORF label. The active site/PTM analysis panel provides alignments generated by Jalview for selected amino acid positions of note that together account for all PTMs observed in this and previous studies, as well as residues implicated in the enzyme's active site. Amino acids detected in proteomics analyses are highlighted in gray; colored residues exhibited PTMs specified in the figure key. PTMs observed in previous studies (Elias and Gygi, 2007; Sharma et al., 2012; Wang et al., 2013) are portrayed by color-coding of amino acid position numbers at the top of the alignment; boxed numbers correspond with active site residues described in Sharma et al. (2012). ORF coverage was calculated by dividing the number of amino acids in an ORF's detected peptides by that ORF's full length. All PTM analysis only includes data from #5133 14N T160d.
The frequent occurrence and differential expression of PTMs among McrA orthologs is a discovery that suggests modifications are common and are phylogenetically and/or environmentally heterogeneous. Functional traits can be substantially affected by the expression of such modifications (Moellering and Cravatt, 2013) and they represent potentially fruitful targets for future work aimed at distinguishing enzymatic modulation or differentiating methanogenic from methanotrophic Mcr in vivo.
The dissimilatory sulfate reduction pathway was well-represented in the metaproteomes from all six seep sediment incubation samples (Figure 6). In total, 17 Sat, 37 AprA, 13 AprB, 11 DsrA, and 18 DsrB proteins were recovered, accounting for 28% of the potential metagenomic search space among these enzyme subunits. Previous protein-based examinations of sulfate reduction in methane seep environments have reported a complete sulfate reduction pathway within the ANME-1 dominated “reefs” in the Black Sea (Basen et al., 2011) and have localized AprB (Wrede et al., 2012), Sat, and Dsr enzymes (Milucka et al., 2013) to ANME-affiliated sulfate SRB via antibody detection. Sulfate reduction pathway proteins were also recovered from Nyegga seep sediments, with the exception of DsrB (Stokke et al., 2012).
Figure 6. Metaproteomic data for enzymes involved in the sulfate reduction pathway. For key, see Figure 4.
Phylogenetic affiliations of sulfate reduction pathway orthologs (Sat, AprA, AprB, DsrA, and DsrB) identified across all six proteomic experiments are shown in Supplemental Data Sheet 4b,c. Of those classified to class level, deltaproteobacterial enzymes (84%) dominated, with Desulfobacterales as the most common family (39% of all class-level partitioned orthologs). Based on the 16S rDNA community data, the majority of the deltaproteobacteria were members of the SEEP-SRB1 and SEEP-SRB2 clades within the Desulfobacterales (Supplemental Data Sheet 3a), an observation supported by FISH visualization (Trembath-Reichert et al., 2016) and the frequent SEEP-SRB association with ANME (Orphan et al., 2002; Knittel et al., 2005).
Sulfate-coupled oxidation of non-methane hydrocarbons can also play a substantial role in sulfur metabolic processing at hydrocarbon seeps (Bose et al., 2013). This process is typically attributed to the DSS clade, with additional contributions from non-DSS Desulfobacteraceae, Desulfobulbaceae, Syntrophobacteraceae, and Desulfurellales (Widdel et al., 2010; Kleindienst et al., 2014). It is difficult to attribute the presence of these organisms exclusively to higher-order alkane oxidation, but Sat enzymes linked to Syntrophobacteraceae were detected, and all of the above lineages were represented in the 16S rRNA gene tag data (Supplemental Data Sheet 3a). 90 (1TP, 1UP) proteins—89 of which did not correspond with previously annotated proteins—demonstrated closest homology to the Desulfosarcina strain BuS5, the only pure culture SRB known to degrade short chain alkanes (Kniemeyer et al., 2007).
Seven thiosulfate reductase orthologs attributed to Delta- and Epsilonproteobacteria were detected (2TP, 1UP), suggesting that thiosulfate reduction is a relatively common metabolic capability within seep sediments. This observation is consistent with previous evidence of thiosulfate reduction from ANME-1-rich seep sediment incubations from Eckernförde Bay (Jagersma et al., 2012) and thiosulfate disproportionation in ANME-1- and ANME-2-dominant samples (Nauhaus et al., 2005). In both cases, the process was seemingly disconnected from AOM, and we observe no evidence of thiosulfate reductase synthesized during our incubations. The connection between AOM-linked sulfate reduction and thiosulfate metabolism remains uncertain, but may impact sulfur cycling in important ways.
Enzymes facilitating sulfur oxidizing pathways were present in the co-localized metagenomes, indicating that the sediments possess a widespread potential for sulfur compound oxidation. Tag sequencing data revealed a substantial proportion of the Sulfurovum genus, accounting for between 4 and 9.6% of total 16S rRNA gene sequences (Supplemental Data Sheet 3a). Sulfurovum species are frequently invoked as sulfur-oxidizing constituents in anoxic zones (Sylvan et al., 2012; Schunck et al., 2013; Urich et al., 2014); their high 16S rRNA gene relative abundance, combined with the low recovery of associated sulfur-oxidizing proteins, suggests that DNA signatures of this lineage are retained over the experimental timescale despite the absence of substantial protein-based biosynthesis.
This study represents one of the most comprehensive SIP proteomics investigations of a complex environmental milieu carried out to date, revealing protein synthesis in slow-growing methane seep sediment communities with a high degree of functional and phylogenetic detail. The detection of multiple functionally relevant orthologs provides a broad sense of in situ active metabolic pathways, while the incorporation of stable isotope probing methods reveals the subset of proteins actively produced under laboratory-based methanotrophic conditions. Components of all of the proteins involved in the reverse methanogenesis pathway were identified, and Mcr was shown to account for 10.4% of all detected peptides. The prevalence of 15N enriched orthologs involved in reverse methanogenesis and sulfate reducing pathways bolsters our understanding that AOM is the dominant biogeochemical process in seep-simulating incubation conditions, indicating that metabolic activity need not scale with community composition as determined by 16S rRNA gene relative abundances. MS-based PTM identification is well-suited to reveal a wider functional diversity than that encoded by nucleic acid sequences alone, a potentially common reality with significant metabolic repercussions. We were able to detect several enriched, pervasive proteins that lack functional annotations, highlighting the utility of proteomic SIP as a discovery platform for ecologically important proteins that should be prioritized for subsequent biochemical characterization. Finally, proteomic SIP can inform more focused proteomic investigations to quantify enzymes of pathways of interest; such values could then be integrated into flux balance models of carbon, sulfur, and nitrogen cycles to better constrain the dynamics and rates of reactions at both the consortium and ecosystem scales. Similar analyses of methanogenic habitats or cultures would provide additional information on the differences between methanogenic and methanotrophic pathways, potentially revealing activity-based controls on methane-linked metabolism.
Tracking the metabolic activity of energy-limited microbial communities has long challenged analytical techniques that address specific microscale anabolic reactions, system-wide metabolic potential, or broad biogeochemical transformations. Proteomic SIP enables the functionally-, phylogenetically-, and temporally-constrained understanding of protein synthesis, opening access to a previously untapped array of distributed low-growth microorganisms active in selected catabolic pathways. Developing a more nuanced appreciation of metabolite flow and interspecies interaction in this way is crucial for the improved understanding and management of microbial ecosystems and their global impacts.
JM, CS, RH, and VO designed the study; JM performed incubation set-up, microbiological and geochemical analyses; CS constructed the metagenomic databases and developed the computational architecture for protein analysis; KC developed protein extraction procedures; ZL and CP performed proteomic SIP computational analysis; JM analyzed proteomic findings from an environmental microbiological perspective and wrote the manuscript with input from all authors.
This work was supported by the US Department of Energy, Office of Science, Office of Biological Environmental Research, under award numbers DE-SC0004949 and DE-SC0010574, and the Life Underground NASA Astrobiology Institute (NNA13AA92A) (to VO). JM was supported by a National Energy Technology Laboratory Methane Hydrate Research Fellowship funded by the National Research Council of the National Academies. The funders had no role in study design, data collection and interpretation, or the manuscript preparation and submission process.
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 thank the Captains, Crew, Alvin group, Jason group, and Science party members from RV Atlantis legs AT-15-68, and AT-18-10. David Case and Stephanie Connon assisted with iTAG data collection and interpretation; Shawn McGlynn performed the nanoSIMS analysis, and Joshua Steele, Rachel Poretsky, and Shulei Sun collected and processed earlier iterations of the metagenomic database. We thank Peter Girguis and the ARPA-E REMOTE team for helpful comments during the preparation of this manuscript. This research used resources of the Oak Ridge Leadership Computing Facility. Oak Ridge National Laboratory is supported by the Office of Science of the U.S. Department of Energy.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.00563
Alain, K., Zbinden, M., Le Bris, N., Lesongeur, F., Quérellou, J., Gaill, F., et al. (2004). Early steps in microbial colonization processes at deep-sea hydrothermal vents. Environ. Microbiol. 6, 227–241. doi: 10.1111/j.1462-2920.2003.00557.x
Alexandre, H., Ansanay-Galeote, V., Dequin, S., and Blondin, B. (2001). Global gene expression during short-term ethanol stress in Saccharomyces cerevisiae. FEBS Lett. 498, 98–103. doi: 10.1016/S0014-5793(01)02503-0
Alperin, M. J., and Hoehler, T. M. (2009). Anaerobic methane oxidation by archaea/sulfate-reducing bacteria aggregates: 1. Thermodynamic and physical constraints. Am. J. Sci. 309, 869–957. doi: 10.2475/10.2009.01
Babin, B. M., Bergkessel, M., Sweredoski, M. J., Moradian, A., Hess, S., Newman, D. K., et al. (2016). SutA is a bacterial transcription factor expressed during slow growth in Pseudomonas aeruginosa. Proc. Natl. Acad. Sci. U.S.A. 113, E597–E605. doi: 10.1073/pnas.1514412113
Basen, M., Krüger, M., Milucka, J., Kuever, J., Kahnt, J., Grundmann, O., et al. (2011). Bacterial enzymes for dissimilatory sulfate reduction in a marine microbial mat (Black Sea) mediating anaerobic oxidation of methane. Environ. Microbiol. 13, 1370–1379. doi: 10.1111/j.1462-2920.2011.02443.x
Benndorf, D., Balcke, G. U., Harms, H., and Von Bergen, M. (2007). Functional metaproteome analysis of protein extracts from contaminated soil and groundwater. ISME J. 1, 224–234. doi: 10.1038/ismej.2007.39
Bidle, K. A., Kastner, M., and Bartlett, D. H. (1999). A phylogenetic analysis of microbial communities associated with methane hydrate containing marine fluids and sediments in the Cascadia margin (ODP site 892B). FEMS Microbiol. Lett. 177, 101–108. doi: 10.1111/j.1574-6968.1999.tb13719.x
Boetius, A., and Suess, E. (2004). Hydrate Ridge: a natural laboratory for the study of microbial life fueled by methane from near-surface gas hydrates. Chem. Geol. 205, 291–310. doi: 10.1016/j.chemgeo.2003.12.034
Boetius, A., Ravenschlag, K., Schubert, C. J., Rickert, D., Widdel, F., Gleseke, 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
Bozinovski, D., Herrmann, S., Richnow, H.H., Von Bergen, M., Seifert, J., and Vogt, C. (2012). Functional analysis of an anaerobic m-xylene-degrading enrichment culture using protein-based stable isotope probing. FEMS Microbiol. Ecol. 81, 134–144. doi: 10.1111/j.1574-6941.2012.01334.x
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Caporaso, J. G., Lauber, C. L., Walters, W. A., Berg-Lyons, D., Lozupone, C. A., Turnbaugh, P. J., et al. (2011). Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. U.S.A. 108, 4516–4522. doi: 10.1073/pnas.1000080107
Case, D., Pasulka, A., Marlow, J., Grupe, B., Levin, L., and Orphan, V. (2015). Methane seep carbonates host distinct, diverse, and dynamic microbial assemblages. mBio 6, e01348-15. doi: 10.1128/mBio.01348-15
Chourey, K., Jansson, J., VerBerkmoes, N., Shah, M., Chavarria, K. L., Tom, L. M., et al. (2010). Direct cellular lysis/protein extraction protocol for soil metaproteomics. J. Proteome Res. 9, 6615–6622. doi: 10.1021/pr100787q
Dekas, A. E., and Orphan, V. J. (2011). Identification of diazotrophic microorganisms in marine sediment via fluorescence in situ hybridization coupled to nanoscale secondary ion mass spectrometry (FISH–NanoSIMS). Meth. Enzymol. 486, 281–305. doi: 10.1016/B978-0-12-381294-0.00012-2
Denef, V. J., Mueller, R. S., and Banfield, J. F. (2010). AMD biofilms: using model communities to study microbial evolution and ecological complexity in nature. ISME J. 4, 599–610. doi: 10.1038/ismej.2009.158
Denef, V. J., VerBerkmoes, N. C., Shah, M. B., Abraham, P., Lefsrud, M., Hettich, R. L., et al. (2009). Proteomics-inferred genome typing (PIGT) demonstrates inter-population recombination as a strategy for environmental adaptation. Environ. Microbiol. 11, 313–325. doi: 10.1111/j.1462-2920.2008.01769.x
Dong, H.-P., Wang, D.-Z., Dai, M., and Hong, H.-S. (2010). Characterization of particulate organic matter in the water column of the South China Sea using a shotgun proteomic approach. Limnol. Oceanogr. 55, 1565–1578. doi: 10.4319/lo.2010.55.4.1565
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
Elias, J. E., and Gygi, S. P. (2007). Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat. Methods 4, 207–214. doi: 10.1038/nmeth1019
Fischer, C. R., Bowen, B. P., Pan, C., Northen, T. R., and Banfield, J. F. (2013). Stable-isotope probing reveals that hydrogen isotope fractionation in proteins and lipids in a microbial community are different and species-specific. ACS Chem. Biol. 8, 1755–1763. doi: 10.1021/cb400210q
Glass, J. B., Yu, H., Steele, J. A., Dawson, K. S., Sun, S., Chourey, K., et al. (2014). Geochemical, metagenomic and metaproteomic insights into trace metal utilization by methane-oxidizing microbial consortia in sulphidic marine sediments. Environ. Microbiol. 16, 1592–1611. doi: 10.1111/1462-2920.12314
Green-Saxena, A., Dekas, A. E., Dalleska, N. F., and Orphan, V. J. (2014). Nitrate-based niche differentiation by distinct sulfate-reducing bacteria involved in the anaerobic oxidation of methane. ISME J. 8, 150–163. doi: 10.1038/ismej.2013.147
Guengerich, F. P., and Cheng, Q. (2011). Orphans in the human cytochrome P450 superfamily: approaches to discovering functions and relevance in pharmacology. Pharmacol. Rev. 63, 684–699. doi: 10.1124/pr.110.003525
Hallam, S. J., Putnam, N., Preston, C. M., Detter, J. C., Rokhsar, D., Richardson, P. M., et al. (2004). Reverse methanogenesis: testing the hypothesis with environmental genomics. Science 305, 1457–1462. doi: 10.1126/science.1100025
Haroon, M. F., Hu, S., Shi, Y., Imelfort, M., Keller, J., Hugenholtz, P., et al. (2013). Anaerobic oxidation of methane coupled to nitrate reduction in a novel archaeal lineage. Nature 500, 567–570. doi: 10.1038/nature12375
Hatzenpichler, R., Scheller, S., Tavormina, P. L., Babin, B. M., Tirrell, D. A., and Orphan, V. J. (2014). In situ visualization of newly synthesized proteins in environmental microbes using amino acid tagging and click chemistry. Environ. Microbiol. 16, 2568–2590. doi: 10.1111/1462-2920.12436
Hawley, A. K., Brewer, H. M., Norbeck, A. D., Paša-Tolić, L., and Hallam, S. J. (2014). Metaproteomics reveals differential modes of metabolic coupling among ubiquitous oxygen minimum zone microbes. Proc. Natl. Acad. Sci. U.S.A. 111, 11395–11400. doi: 10.1073/pnas.1322132111
Hebert, A. S., Richards, A. L., Bailey, D. J., Ulbrich, A., Coughlin, E. E., Westphall, M. S., et al. (2014). The one hour yeast proteome. Mol. Cell. Proteomics 13, 339–347. doi: 10.1074/mcp.M113.034769
Hyatt, D., Chen, G.-L., LoCascio, P. F., Land, M. L., Larimer, F. W., and Hauser, L. J. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 11:119. doi: 10.1186/1471-2105-11-119
Imelfort, M., Parks, D., Woodcroft, B. J., Dennis, P., Hugenholtz, P., and Tyson, G. W. (2014). GroopM: an automated tool for the recovery of population genomes from related metagenomes. PeerJ 2:e603. doi: 10.7717/peerj.603
Iversen, N., and Jorgensen, B. B. (1985). Anaerobic methane oxidation rates at the sulfate-methane transition in marine sediments from Kattegat and Skagerrak (Denmark)1. Limnol. Oceanogr. 30, 944–955. doi: 10.4319/lo.1985.30.5.0944
Jagersma, C. G., Meulepas, R. J. W., Timmers, P. H. A., Szperl, A., Lens, P. N. L., and Stams, A. J. M. (2012). Enrichment of ANME-1 from Eckernförde Bay sediment on thiosulfate, methane and short-chain fatty acids. J. Biotechnol. 157, 482–489. doi: 10.1016/j.jbiotec.2011.10.012
Jehmlich, N., Schmidt, F., Von Bergen, M., Richnow, H.-H., and Vogt, C. (2008). Protein-based stable isotope probing (Protein-SIP) reveals active species within anoxic mixed cultures. ISME J. 2, 1122–1133. doi: 10.1038/ismej.2008.64
Justice, N. B., Li, Z., Wang, Y., Spaudling, S. E., Mosier, A. C., Hettich, R. L., et al. (2014). 15N- and 2H proteomic stable isotope probing links nitrogen flow to archaeal heterotrophic activity. Environ. Microbiol. 16, 3224–3237. doi: 10.1111/1462-2920.12488
Justice, N. B., Pan, C., Mueller, R., Spaulding, S. E., Shah, V., Sun, C. L., et al. (2012). Heterotrophic archaea contribute to carbon cycling in low-ph, suboxic biofilm communities. Appl. Environ. Microbiol. 78, 8321–8330. doi: 10.1128/AEM.01938-12
Kerner, M. J., Naylor, D. J., Ishihama, Y., Maier, T., Chang, H.-C., Stines, A. P., et al. (2005). Proteome-wide analysis of chaperonin-dependent protein folding in Escherichia coli. Cell 122, 209–220. doi: 10.1016/j.cell.2005.05.028
Kim, Y. E., Hipp, M. S., Bracher, A., Hayer-Hartl, M., and Ulrich Hartl, F. (2013). Molecular chaperone functions in protein folding and proteostasis. Annu. Rev. Biochem. 82, 323–355. doi: 10.1146/annurev-biochem-060208-092442
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
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–901. doi: 10.1038/nature06200
Knittel, K., Lösekann, T., Boetius, A., Kort, R., and Amann, R. (2005). Diversity and distribution of methanotrophic archaea at cold seeps. Appl. Environ. Microbiol. 71, 467–479. doi: 10.1128/AEM.71.1.467-479.2005
Krüger, M., Blumenberg, M., Kasten, S., Wieland, A., Känel, L., Klock, J.-H., et al. (2008a). A novel, multi-layered methanotrophic microbial mat system growing on the sediment of the Black Sea. Environ. Microbiol. 10, 1934–1947. doi: 10.1111/j.1462-2920.2008.01607.x
Krüger, M., Meyerdierks, A., Glöckner, F. O., Amann, R., Widdel, F., Kube, M., et al. (2003). A conspicuous nickel protein in microbial mats that oxidize methane anaerobically. Nature 426, 878–881. doi: 10.1038/nature02207
Krüger, M., Wolters, H., Gehre, M., Joye, S. B., and Richnow, H. (2008b). Tracing the slow growth of anaerobic methane-oxidizing communities by 15N-labelling techniques. FEMS Microbiol. Ecol. 63, 401–411. doi: 10.1111/j.1574-6941.2007.00431.x
Li, Z., Wang, Y., Yao, Q., Justice, N. B., Ahn, T.-H., Xu, D., et al. (2014). Diverse and divergent protein post-translational modifications in two growth stages of a natural microbial community. Nat. Commun. 5. doi: 10.1038/ncomms5405
Marlow, J. J., Steele, J. A., Case, D. H., Connon, S. A., Levin, L. A., and Orphan, V. J. (2014b). Microbial abundance and diversity patterns associated with sediments and carbonates from the methane seep environments of Hydrate Ridge, OR. Front. Mar. Sci 1:44. doi: 10.3389/fmars.2014.00044
Marlow, J. J., Steele, J. A., Ziebis, W., Thurber, A. R., Levin, L. A., and Orphan, V. J. (2014a). Carbonate-hosted methanotrophy represents an unrecognized methane sink in the deep sea. Nat. Commun. 5, doi: 10.1038/ncomms6094
Meyerdierks, A., Kube, M., Kostadinov, I., Teeling, H., Glöckner, F. O., Reinhardt, R., et al. (2010). Metagenome and mRNA expression analyses of anaerobic methanotrophic archaea of the ANME-1 group. Environ. Microbiol. 12, 422–439. doi: 10.1111/j.1462-2920.2009.02083.x
Milucka, J., Widdel, F., and Shima, S. (2013). Immunological detection of enzymes for sulfate reduction in anaerobic methane-oxidizing consortia. Environ. Microbiol. 15, 1561–1571. doi: 10.1111/1462-2920.12003
Mitchell, A., Chang, H.-Y., Daugherty, L., Fraser, M., Hunter, S., Lopez, R., et al. (2015). The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Res. 43, D213–D221. doi: 10.1093/nar/gku1243
Mohr, W., Tang, T., Sattin, S. R., Bovee, R. J., and Pearson, A. (2014). Protein stable isotope fingerprinting: multidimensional protein chromatography coupled to stable isotope-ratio mass spectrometry. Anal. Chem. 86, 8514–8520. doi: 10.1021/ac502494b
Moore, E. K., Nunn, B. L., Goodlett, D. R., and Harvey, H. R. (2012). Identifying and tracking proteins through the marine water column: Insights into the inputs and preservation mechanisms of protein in sediments. Geochim. Cosmochim. Acta 83, 324–359. doi: 10.1016/j.gca.2012.01.002
Morgan, R., Kohn, S., Hwang, S.-H., Hassett, D. J., and Sauer, K. (2006). BdlA, a Chemotaxis Regulator Essential for Biofilm Dispersion in Pseudomonas aeruginosa. J. Bacteriol. 188, 7335–7343. doi: 10.1128/JB.00599-06
Moriya, Y., Itoh, M., Okuda, S., Yoshizawa, A. C., and Kanehisa, M. (2007). KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 35, W182–W185. doi: 10.1093/nar/gkm321
Morono, Y., Ito, M., and Inagaki, F. (2014). “Detecting slow metabolism in the subseafloor: analysis of single cells using NanoSIMS, in Microbial Life of the Deep Biosphere, eds J. Kallmeyer and D. Wagner (Berlin: De Gruyter), 101–120.
Nauhaus, K., Treude, T., Boetius, A., and Krüger, M. (2005). Environmental regulation of the anaerobic oxidation of methane: a comparison of ANME-I and ANME-II communities. Environ. Microbiol. 7, 98–106. doi: 10.1111/j.1462-2920.2004.00669.x
Nguyen, R. T., and Harvey, H. R. (2001). Preservation of protein in marine systems: hydrophobic and other noncovalent associations as major stabilizing forces. Geochim. Cosmochim. Acta 65, 1467–1480. doi: 10.1016/S0016-7037(00)00621-9
Olsen, J. V., Vermeulen, M., Santamaria, A., Kumar, C., Miller, M. L., Jensen, L. J., et al. (2010). Quantitative phosphoproteomics reveals widespread full phosphorylation site occupancy during mitosis. Sci. Signal. 3, ra3. doi: 10.1126/scisignal.2000475
Orphan, V. J., House, C. H., Hinrichs, K.-U., McKeegan, K. D., and DeLong, E. F. (2002). Multiple archaeal groups mediate methane oxidation in anoxic cold seep sediments. Proc. Natl. Acad. Sci. U.S.A. 99, 7663–7668. doi: 10.1073/pnas.072210299
Orphan, V. J., Turk, K. A., Green, A. M., and House, C. H. (2009). Patterns of 15N assimilation and growth of methanotrophic ANME-2 archaea and sulfate-reducing bacteria within structured syntrophic consortia revealed by FISH-SIMS. Environ. Microbiol. 11, 1777–1791. doi: 10.1111/j.1462-2920.2009.01903.x
Pan, C., Fischer, C. R., Hyatt, D., Bowen, B. P., Hettich, R. L., and Banfield, J. F. (2011). Quantitative tracking of isotope flows in proteomes of microbial communities. Mol. Cell. Proteomics 10:M110.006049 doi: 10.1074/mcp.M110.006049
Parkes, R., Cragg, B., Fry, J., Herbert, R., Wimpenny, J., Allen, J., et al. (1990). Bacterial biomass and activity in deep sediment layers from the Peru Margin [and Discussion]. Philos. Trans. R. Soc. Lond. A Math. Phys. Sci. 331, 139–153. doi: 10.1098/rsta.1990.0061
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2014). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114
Pernthaler, A., Dekas, A. E., Brown, C. T., Goffredi, S. K., Embaye, T., and Orphan, V. J. (2008). Diverse syntrophic partnerships from deep-sea methane vents revealed by direct cell capture and metagenomics. Proc. Natl. Acad. Sci. U.S.A. 105, 7052–7057. doi: 10.1073/pnas.0711303105
Ram, R. J., VerBerkmoes, N. C., Thelen, M. P., Tyson, G. W., Baker, B. J., Blake, R. C., et al. (2005). Community proteomics of a natural microbial biofilm. Science 308, 1915–1920. doi: 10.1126/science. 1109070
Sahling, H., Rickert, D., Lee, R. W., Linke, P., and Suess, E. (2002). Macrofaunal community structure and sulfide flux at gas hydrate deposits from the Cascadia convergent margin, NE Pacific. Mar. Ecol. Prog. Ser. 231, 121–138. doi: 10.3354/meps231121
Schreiber, L., Holler, T., Knittel, K., Meyerdierks, A., and Amann, R. (2010). Identification of the dominant sulfate-reducing bacterial partner of anaerobic methanotrophs of the ANME-2 clade. Environ. Microbiol. 12, 2327–2340. doi: 10.1111/j.1462-2920.2010.02275.x
Schunck, H., Lavik, G., Desai, D. K., Großkopf, T., Kalvelage, T., Löscher, C. R., et al. (2013). Giant hydrogen sulfide plume in the oxygen minimum zone off Peru supports chemolithoautotrophy. PLoS ONE 8:e68661. doi: 10.1371/journal.pone.0068661
Seifert, J., Taubert, M., Jehmlich, N., Schmidt, F., Völker, U., Vogt, C., et al. (2012). Protein-based stable isotope probing (protein-SIP) in functional metaproteomics. Mass Spectrom. Rev. 31, 683–697. doi: 10.1002/mas.21346
Selmer, T., Kahnt, J., Goubeaud, M., Shima, S., Grabarse, W., Ermler, U., et al. (2000). The biosynthesis of methylated amino acids in the active site region of Methyl-coenzyme M Reductase. J. Biol. Chem. 275, 3755–3760. doi: 10.1074/jbc.275.6.3755
Sharma, R., Dill, B. D., Chourey, K., Shah, M., VerBerkmoes, N. C., and Hettich, R. L. (2012). Coupling a detergent lysis/cleanup methodology with intact protein fractionation for enhanced proteome characterization. J. Proteome Res. 11, 6008–6018. doi: 10.1021/pr300709k
Sievert, S. M., Wieringa, E. B., Wirsen, C. O., and Taylor, C. D. (2007). Growth and mechanism of filamentous-sulfur formation by Candidatus Arcobacter sulfidicus in opposing oxygen-sulfide gradients. Environ. Microbiol. 9, 271–276. doi: 10.1111/j.1462-2920.2006.01156.x
Sowell, S. M., Wilhelm, L. J., Norbeck, A. D., Lipton, M. S., Nicora, C. D., Barofsky, D. F., et al. (2009). Transport functions dominate the SAR11 metaproteome at low-nutrient extremes in the Sargasso Sea. ISME J. 3, 93–105. doi: 10.1038/ismej.2008.83
Stevens, T. J., and Arkin, I. T. (2000). Do more complex organisms have a greater proportion of membrane proteins in their genomes? Proteins 39, 417–420. doi: 10.1002/(SICI)1097-0134(20000601)39:4 < 417::AID-PROT140>3.0.CO;2-Y
Stokke, R., Roalkvam, I., Lanzen, A., Haflidason, H., and Steen, I. H. (2012). Integrated metagenomic and metaproteomic analyses of an ANME-1-dominated community in marine cold seep sediments. Environ. Microbiol. 14, 1333–1346. doi: 10.1111/j.1462-2920.2012.02716.x
Suess, E., Torres, M., Bohrmann, G., Collier, R., 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
Suzek, B. E., Huang, H., McGarvey, P., Mazumder, R., and Wu, C. H. (2007). UniRef: comprehensive and non-redundant UniProt reference clusters. Bioinformatics 23, 1282–1288. doi: 10.1093/bioinformatics/btm098
Sylvan, J. B., Toner, B. M., and Edwards, K. J. (2012). Life and death of deep-sea vents: bacterial diversity and ecosystem succession on inactive hydrothermal sulfides. MBio 3:e00279-11. doi: 10.1128/mBio.00279-11
Thauer, R. K., Kaster, A.-K., Seedorf, H., Buckel, W., and Hedderich, R. (2008). Methanogenic archaea: ecologically relevant differences in energy conservation. Nat. Rev. Microbiol. 6, 579–591. doi: 10.1038/nrmicro1931
Toner, B. M., Lesniewski, R. A., Marlow, J. J., Briscoe, L. J., Santelli, C. M., Bach, W., et al. (2012). Mineralogy drives bacterial biogeography of hydrothermally inactive seafloor sulfide deposits. Geomicrobiol. J. 30, 313–326. doi: 10.1080/01490451.2012.688925
Trembath-Reichert, E., Case, D. H., and Orphan, V. J. (2016). Characterization of microbial associations with methanotrophic archaea and sulfate-reducing bacteria through statistical comparison of nested Magneto-FISH enrichments. Peer J. 4:e1913. doi: 10.7717/peerj.1913
Treude, T., Boetius, A., Knittel, K., Wallmann, K., and Jorgensen, 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
Tryon, M., Brown, K., and Torres, M. (2002). Fluid and chemical flux in and out of sediments hosting methane hydrate deposits on Hydrate Ridge, OR, II: Hydrological processes. Earth Planet. Sci. Lett. 201, 541–557. doi: 10.1016/S0012-821X(02)00732-X
Urich, T., Lanzén, A., Stokke, R., Pedersen, R. B., Bayer, C., Thorseth, I. H., et al. (2014). Microbial community structure and functioning in marine sediments associated with diffuse hydrothermal venting assessed by integrated meta-omics. Environ. Microbiol. 16, 2699–2710. doi: 10.1111/1462-2920.12283
Verberkmoes, N. C., Russell, A. L., Shah, M., Godzik, A., Rosenquist, M., Halfvarson, J., et al. (2009). Shotgun metaproteomics of the human distal gut microbiota. ISME J. 3, 179–189. doi: 10.1038/ismej.2008.108
Wagner, M. (2009). Single-cell ecophysiology of microbes as revealed by Raman microspectroscopy or secondary ion mass spectrometry imaging. Annu. Rev. Microbiol. 63, 411–429. doi: 10.1146/annurev.micro.091208.073233
Wang, F.-P., Zhang, Y., Chen, Y., He, Y., Qi, J., Hinrichs, K.-U., et al. (2014). Methanotrophic archaea possessing diverging methane-oxidizing and electron-transporting pathways. ISME J. 8, 1069–1078. doi: 10.1038/ismej.2013.212
Wang, Y., Ahn, T.-H., Li, Z., and Pan, C. (2013). Sipros/ProRata: a versatile informatics system for quantitative community proteomics. Bioinformatics 29, 2064–2065. doi: 10.1093/bioinformatics/btt329
Waterhouse, A. M., Procter, J. B., Martin, D. M., Clamp, M., and Barton, G. J. (2009). Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191. doi: 10.1093/bioinformatics/btp033
Webb, K. J., Xu, T., Park, S. K., and Yates, J. R. III. (2013). Modified MuDPIT separation identified 4488 proteins in a system-wide analysis of quiescence in yeast. J. Proteome Res. 12, 2177–2184. doi: 10.1021/pr400027m
Welander, P. V., and Metcalf, W. W. (2008). Mutagenesis of the C1 oxidation pathway in Methanosarcina barkeri: new insights into the Mtr/Mer Bypass Pathway. J. Bacteriol. 190, 1928–1936. doi: 10.1128/JB.01424-07
Widdel, F., Knittel, K., and Galushko, A. (2010). “Anaerobic hydrocarbon-degrading microorganisms: an overview, in Handbook of Hydrocarbon and Lipid Microbiology, ed K. Timmis (Berlin; Heidelberg: Springer), 1997–2021.
Wrede, C., Krukenberg, V., Dreier, A., Reitner, J., Heller, C., and Hoppert, M. (2012). Detection of metabolic key enzymes of methane turnover processes in cold seep microbial biofilms. Geomicrobiol. J. 30, 214–227. doi: 10.1080/01490451.2012.665150
Keywords: proteomics, stable isotope probing, methane seeps, anaerobic oxidation of methane, environmental microbiology
Citation: Marlow JJ, Skennerton CT, Li Z, Chourey K, Hettich RL, Pan C and Orphan VJ (2016) Proteomic Stable Isotope Probing Reveals Biosynthesis Dynamics of Slow Growing Methane Based Microbial Communities. Front. Microbiol. 7:563. doi: 10.3389/fmicb.2016.00563
Received: 21 February 2016; Accepted: 04 April 2016;
Published: 29 April 2016.
Edited by:Mark Alexander Lever, ETH Zürich, Switzerland
Reviewed by:Jennifer F. Biddle, University of Delaware, USA
Ronald Oremland, United States Geological Survey, USA
Copyright © 2016 Marlow, Skennerton, Li, Chourey, Hettich, Pan and Orphan. 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.
†Present Address: Jeffrey J. Marlow, Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA