Transcriptional response of bathypelagic marine bacterioplankton to the Deepwater Horizon oil spill

ISME Journal, Aug 2013

The Deepwater Horizon blowout released a massive amount of oil and gas into the deep ocean between April and July 2010, stimulating microbial blooms of petroleum-degrading bacteria. To understand the metabolic response of marine microorganisms, we sequenced ∼66 million community transcripts that revealed the identity of metabolically active microbes and their roles in petroleum consumption. Reads were assigned to reference genes from ∼2700 bacterial and archaeal taxa, but most assignments (39%) were to just six genomes representing predominantly methane- and petroleum-degrading Gammaproteobacteria. Specific pathways for the degradation of alkanes, aromatic compounds and methane emerged from the metatranscriptomes, with some transcripts assigned to methane monooxygenases representing highly divergent homologs that may degrade either methane or short alkanes. The microbial community in the plume was less taxonomically and functionally diverse than the unexposed community below the plume; this was due primarily to decreased species evenness resulting from Gammaproteobacteria blooms. Surprisingly, a number of taxa (related to SAR11, Nitrosopumilus and Bacteroides, among others) contributed equal numbers of transcripts per liter in both the unexposed and plume samples, suggesting that some groups were unaffected by the petroleum inputs and blooms of degrader taxa, and may be important for re-establishing the pre-spill microbial community structure.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

Transcriptional response of bathypelagic marine bacterioplankton to the Deepwater Horizon oil spill

The ISME Journal Transcriptional response of bathypelagic marine bacterioplankton to the Deepwater Horizon oil spill Adam R Rivers 1 2 Shalabh Sharma 1 2 Susannah G Tringe 0 1 Jeffrey Martin 0 1 Samantha B Joye 1 2 Mary Ann Moran 1 2 0 DOE Joint Genome Institute , Walnut Creek, CA , USA 1 University of Georgia , Athens, GA 30602-3636 , USA 2 Department of Marine Sciences, University of Georgia , Athens, GA , USA The Deepwater Horizon blowout released a massive amount of oil and gas into the deep ocean between April and July 2010, stimulating microbial blooms of petroleum-degrading bacteria. To understand the metabolic response of marine microorganisms, we sequenced B66 million community transcripts that revealed the identity of metabolically active microbes and their roles in petroleum consumption. Reads were assigned to reference genes from B2700 bacterial and archaeal taxa, but most assignments (39%) were to just six genomes representing predominantly methane- and petroleum-degrading Gammaproteobacteria. Specific pathways for the degradation of alkanes, aromatic compounds and methane emerged from the metatranscriptomes, with some transcripts assigned to methane monooxygenases representing highly divergent homologs that may degrade either methane or short alkanes. The microbial community in the plume was less taxonomically and functionally diverse than the unexposed community below the plume; this was due primarily to decreased species evenness resulting from Gammaproteobacteria blooms. Surprisingly, a number of taxa (related to SAR11, Nitrosopumilus and Bacteroides, among others) contributed equal numbers of transcripts per liter in both the unexposed and plume samples, suggesting that some groups were unaffected by the petroleum inputs and blooms of degrader taxa, and may be important for re-establishing the pre-spill microbial community structure. The ISME Journal (2013) 7, 2315-2329; doi:10.1038/ismej.2013.129; published online 1 August 2013 Subject Category: Integrated genomics and post-genomics approaches in microbial ecology Introduction The Deepwater Horizon blowout released 6.9 1011 g of oil and gas into Gulf of Mexico (GOM) waters (Reddy et al., 2012) . Approximately one-third of the hydrocarbons initially concentrated in a lens near 1100 m depth (Camilli et al., 2010; Ryerson et al., 2012) , inducing a shift in the composition of the microbial community and the proliferation of Gammaproteobacteria related to Oceanospirillales, Colwellia and Cycloclasticus (Hazen et al., 2010) . Previous work has examined community composition and the abundance of functional genes and transcripts (Hazen et al., 2010; Lu et al,. 2011; Mason et al., 2012) , as well as microbial consumption of hydrocarbons (Redmond and Valentine, 2012) . This paper extends that work by analyzing B66 million microbial messenger RNAs (mRNAs) and 4350 000 16S rRNA genes for a comprehensive metabolic and ecological reconstruction of microbial processes in response to the oil and gas inputs. The petroleum released from the Macondo Well contained short- to long chain alkanes, aromatic compounds, polycyclic aromatic hydrocarbons, asphaltenes and polar compounds. Methane accounted for 15% of the total petroleum and gas leaving the riser pipe, short alkanes excluding methane (C2–C5) comprised 9% and medium chain alkanes (n-C6 to n-C20) comprised 13% (Reddy et al., 2012) . High rate constants were reported for the oxidation of ethane and propane (Valentine et al., 2010) , while bacteria enriched from the subsurface plume had the potential to degrade long chain n-alkanes and branched alkanes (Hazen et al., 2010) . The role of bacteria in methane metabolism is less clear, however, as methane may have been removed by microbial degradation or by venting to the atmosphere. Methane oxidation rates measured in May (3 weeks after the uncontrolled blowout began) were high (Joye, unpublished data), but rate constants for methane oxidation measured a few weeks later and again in August and September (after the well had been capped) were low (Valentine et al., 2010; Kessler et al., 2011) and functional gene microarrays showed no enrichment of methane monooxygenase genes in the plume communities (Lu et al., 2011). Counts of 16S rRNA genes diagnostic for methanotrophs varied depending on the timing of sampling and methodology (Hazen et al., 2010; Valentine et al., 2010; Kessler et al., 2011; Redmond and Valentine, 2012; Mason et al., 2012) . Estimates of atmospheric venting, the alternative route for methane loss, did not resolve the gas fate question as flux estimates were high when based on surface water measurements (Joye, unpublished data), but low when based on remote-sensing data (Yvon-Lewis et al., 2011). This study builds on the previous meta-omic analyses conducted following the Deepwater Horizon accident (Hazen et al., 2010; Lu et al., 2011; Mason et al., 2012) by generating a sizeable transcript data set with low levels of rRNA contamination and robust coverage of more than 561 000 functional genes from exposed and unexposed communities. The high coverage of expressed oxygenase sequences allowed application of evolutionary placement methods to address the relative importance of methane, ethane and propane metabolism, whereas the use of internal standards generated information on changes in transcript numbers per liter. These approaches created a detailed picture of the metabolic and ecological responses of the bathypelagic microbial community. Materials and methods Microbial mRNA was sequenced from four samples collected at two different stations located B6–8 km from the wellhead. Two samples were collected within the oil/gas plume (samples P16 at 1116 m and P52 at 1198 m) from areas with elevated colored dissolved organic matter (CDOM, a proxy for oil) fluorescence and oxygen depletion. Two samples were collected below the plume, and although intended to represent two non-impacted communities, the sample IP16 from 1240 m was adjacent to an oxygen anomaly and subsequent data analysis suggested that it represented an intermediate condition. The sample NP52 from 1286 m was the least affected by the oil/gas release, and although it did not have detectable oil or oxygen anomalies by instruments on the conductivity, temperature and depth (CTD) rosette instrumentation, gas chromatographic analysis revealed that methane was elevated above background levels (Table 1). This paired-in/ below-plume sampling scheme minimized the influence of geographic location. In addition to the transcript data, seventeen 16S rRNA gene libraries were created from seven stations at depths above, in and below the plume. An oceanographic station map and depth profiles (Supplementary Figure S1) were plotted using the R package OCE (Kelley, 2013) . Sample collection and sequencing Samples were collected between 26 May 2010 and 3 June 2010 aboard the RV Walton Smith using a rinsed CTD rosette. For RNA analysis, 9 l of seawater was immediately filtered through a 0.45-mm-pore Physical Depth (m) Hydrocarbon plume Oxygen anomaly Chemicalw CDOM (relative fluorescence units) Methane (nmol l 1) Nitrate (mmol l 1) Nitrite (mmol l 1) Phosphate (mmol l 1) Biological Cell density (cells ml 1) Total RNA (mg per l seawater) Sequencing Median read length (bases) Total joined reads rRNA reads (%) Internal standard reads (%) Possible proteins (%) Archaeal and bacterial hits RefSeq genes RefSeq taxa KEGG assignments P16 1116 Yes Yes size, 142-mm-diameter polyethersulfone membrane filter (Pall Corporation, Port Washington, NY, USA). The larger pore size was chosen to minimize clogging by oil droplets and to keep average filtration times to o10 min to better preserve mRNA integrity, although it may have biased cell collection toward larger cells or those associated with oil droplets. All tubing, filters and containers were soaked in 1 M HCl and rinsed in distilled water before sampling. Filters were immediately preserved in RNAlater (Life Technologies, Carlsbad, CA, USA) and kept at 4 1C. A smaller volume of water (0.5–1 l) was filtered through a 0.45-mm pore size, 47-mm diameter polyethersulfone filter and stored on dry ice for subsequent DNA extraction. RNA was extracted according to Poretsky et al. (2009) . Six internal RNA standards (artificial mRNAs synthesized by in vitro transcription of DNA templates) were added as internal controls (1.4 1010 transcripts per sample) (Satinsky et al., 2013) . The two smallest standards (B200 nt) were not appreciably recovered on the solid phase extraction matrix used for RNA purification and were not considered further. RNA extracts were depleted of rRNA by creating sample-specific removal probes targeting the large and small subunits of archaeal, bacterial and eukaryotic rRNA genes using the DNA of each sample (Stewart et al., 2010) . The protocol was modified by pooling probes from multiple DNA samples (see below) and including an additional round of probe removal. rRNA-depleted mRNA was amplified using the MessageAmp II bacterial kit (Life Technologies) and complementary DNA was synthesized using the universal Riboclone complementary DNA synthesis system (Promega, Madison, WI, USA). Complementary DNA was sheared ultrasonically to B225 bp, and each sample was run on one lane of an Illumina Genome Analyzer IIx (Illumina Inc., San Diego, CA, USA) using 150-bp paired-end chemistry. Sample processing and bioinformatic analysis is summarized in Supplementary Figure S2. DNA was extracted from samples using the DNA Power Water Kit (Mo Bio Laboratories Inc., Carlsbad, CA, USA) and further purified by ethanol precipitation. The V6–V8 region of the 16S rRNA gene was amplified using the primers 926F (50-cctatcccctg tgtgccttggcagtctcag AAA CTY AAA KGA ATT GRC GG-30) and 1392R (50-ccatctcatccctgcgtgtctccgactcag[barcode] -ACGGGCGGTGTGTRC-30), with 454A or B adapter sequences indicated in lower case. PCR amplicons were purified using Ampure SPRI Beads (Beckman Coulter, Brea, CA, USA). Emulsion PCR and sequencing of the PCR amplicons were performed following the Roche/454 (Branford, CT, USA) GS FLX Titanium technology instructions. Bioinformatics 16S rRNA gene sequence data were analyzed through the QIIME pipeline, version 1.3 (Caporaso et al., 2010) . Sequencing artifacts were minimized using the Denoiser package of QIIME (Reeder and Knight, 2010). Default settings were used except that the Greengenes database (DeSantis et al., 2006) was used to train the Ribosomal Database Project classifier (Wang et al., 2007) . Technical replicates showed good agreement and the technical replicate with the highest number of reads was used for analysis. The Unifrac distance matrix (Lozupone et al., 2011) generated by QIIME was used for singlelinkage clustering of the sequences and for nonmetric multidimensional scaling and environmental variable fitting in the R package Vegan (Jari et al., 2013) . Transcriptomic data were checked for quality using FastX-Toolkit (Gordon and Hannon, 2011, unpublished data). Paired-end reads were joined using the program SHERA with default parameters (Rodrigue et al., 2010) and a quality metric of 0.5. Joined reads were trimmed using Seqtrim (Falgueras and Lara, 2010) with default values. To remove remaining rRNA and the internal standard sequences, a Blastn search was performed against a small database containing B1000 sequences that included the internal standards and representative ribosomal and tRNA sequences; matching reads with a bit score 450 were removed (Gifford et al., 2011) . The remaining paired and trimmed reads were queried against RefSeq protein database version 45 using Blastx for the best match with a bit score 440. Reads were mapped to KEGG ortholog groups (Kanehisa and Goto, 2000) and COG categories (Tatusov et al., 2003) by querying the matching RefSeq gene against the KEGG gene database or the COG gene database (bit scores 440). Blast results were stored for analysis in a MySQL database. Internal standards were tallied and used to measure sequencing accuracy, as well as to estimate absolute numbers of transcripts (transcripts l 1; Gifford et al., 2011; Satinsky et al., 2013) , except for sample P16 that had an error in the amount of standard added. The plot of log-fold change and COG category assignment was generated with Circos (Krzywinski et al., 2009). The analysis of differential abundance was carried out using the R package DESeq (Anders and Huber, 2010) by comparing P16 and P52 (plume samples) with NP52 (non-plume sample); IP16 was not included because of its intermediate physical and biological characteristics. Significance was calculated using dispersion values estimated from the two plume samples, as no replicate was available for NP52, and was defined as a P-value of o0.05 in a negative binomial test following correction for false discovery rate (Benjamini and Hochberg, 1995) . For pathway analysis, statistical significance was calculated based on the abundance of reads in KEGG orthology groups (K numbers). For organism-specific analyses, significance was calculated for each RefSeq gene present for the reference genome. The pathways inferred to be involved in hydrocarbon degradation were generated by modifying existing pathway maps (Anthony, 1982; Caspi et al., 2012) with manual curation. The de novo assembly of operons was carried out using the Rnnotater pipeline (Martin et al., 2010) ; the mean assembly length was 886 bp (Supplementary Figure S3). The candidate particulate methane monooxygenase (pmoCAB) operons in the assembled data were identified by querying known pmoA and pxmA genes from Methylobacter tundripaludum against the database of transcript assemblies using Tblastn with an expect value cutoff of 10 5. Open reading frames were identified using MetaGeneMark (Zhu et al., 2010) . Full-length monooxygenase subunit A ORFs were extracted from the operons with Blastx and aligned to known pmoA, pxmA and hmoA genes using MAFFT (Katoh et al., 2002) (L-INS-I algorithm, Blosum62 matrix) and a high-quality alignment region consisting of 232 amino-acid sites was selected using TrimAl (Capella-Gutierrez et al., 2009) . PhyML mixture was used to generate a maximum likelihood mixture model tree using the UL3 substitution matrices with a gamma rate model, assuming four rate categories and estimating the invariant sites (Le et al., 2008) . Hmmalign was run with the pmoA/amoA hidden Markov model (PF02461) to align 33 257 short reads, 5 short reference sequences and 28 full-length reference sequences (Eddy, 2011) . After alignment, the reads were placed on the tree with the maximum likelihood placement algorithm in RAxML (Stamatakis, 2006) . The results were converted to phyloXML format (Han and Zmasek, 2009) by pplacer/guppy (Matsen et al., 2010) and visualized in Archaeopteryx (Han and Zmasek, 2009). The tree of pyrroloquinoline quinone (PQQ) alcohol and methanol dehydrogenases was constructed using the same method. The site variability of assembled pmoA genes was assessed by realigning reads from evolutionary tree placement to the reference sequence using Geneious version 5.3.6 (Biomatters, Auckland, NZ, USA). The number of positions with variation (defined as 425% of reads different from the assembly) ranged from 3 (DH4550) to 120 (DH311), and variation co-occurred on reads with other variations, suggesting that reads within pmoA assemblies originated from just a few closely related bacteria. To assess the fate of the pre-spill microbial community, transcripts per liter were estimated across samples for 50 taxa that were most abundant in sample NP52. The ratios of transcripts per liter for NP52:P52 and NP52:IP16 were used to create a hierarchical clustering of the taxa (complete linkage, Manhattan distance). Two major clusters of the tree were tested to see if there were significant differences in the number of reads per sample; testing was done with a Welch corrected one-way analysis of variance on the log-transformed read counts. Sequences from metatranscriptomic and 16S rRNA sequencing are available from the Community Cyberinfrastructure for Advanced Microbial Ecology Research and Analysis database under project number CAM_P_0001088. The pmoA homolog assemblies are available at the National Center for Biotechnology Information database under accession numbers KF296332-KF296332. Results The depth of the Deepwater Horizon hydrocarbon plume was identified at seven stations in the vicinity of the well by CDOM fluorescence and low-oxygen anomalies (Supplementary Figure S1 B-H, Table 1) (Joye et al., 2011) . Seventeen seawater samples were collected at depths from directly in, above or below the plume for analysis of bacterial community composition based on amplified 16S rRNA genes. Four of these samples were selected for further metatranscriptome analysis as follows: cast 16 at 1116 and 1198 m, and cast 52 at 1240 and 1286 m. Chemical data indicated that two of the metatranscriptome samples represented plume conditions (characterized by both a hydrocarbon signal and an oxygen anomaly; these were designated as P16 and P52) and one represented largely non-impacted seawater (below both the hydrocarbon and oxygen signals; NP52). The fourth sample, initially selected as a second below-plume metatranscriptome replicate because of a low hydrocarbon signal, was later found to have a small oxygen anomaly and was designated as an intermediate plume sample (IP16); this sample might represent a region receiving low hydrocarbon inputs or an aged portion of the main plume. The two in-plume samples had greater biological activity than the non-plume or intermediate plume samples, with cell concentrations that were almost 10-fold higher (geometric mean in plume: 2.6 105, intermediate/non-plume: 2.4 104 cells ml 1) and total RNA yields that were B200-fold higher (Table 1, Supplementary Figure S4). 16S rRNA composition The 17 microbial communities were characterized based on 4350 000 16S rRNA amplicon sequences (Figure 1b), and composition was compared using nonmetric multidimensional scaling of a Unifrac taxonomic distance matrix (Lozupone et al., 2011) . The communities clustered in patterns that were generally consistent with CDOM and oxygen anomaly features (Figure 1d). Among the plume samples, four bacterial families accounted for 95% of the amplicon sequences. These had closest 16S rRNA sequence matches to Colwelliaceae (genus Colwellia), Oceanospirillaceae (genera Neptuniibacter and Bermanella), Piscirickettsiaceae (genus Cycloclasticus) and Methylococcaceae (genera Methylobacter and Methylococcus) families (Figure 1b). The importance of these taxonomic groups in the samples was significantly positively correlated with increased concentrations of methane and oil, and significantly negatively correlated with concentrations of dissolved oxygen (Po0.04 for each factor, Figure 1d). In the non-plume sample NP52, the top four taxonomic groups were members of the family Oceanospirillaceae (28%), the class Alphaproteobacteria (11%), the SAR406 group (14%) and the class Deltaproteobacteria (11%). The composition of the samples selected for metatranscriptome sequencing was representative of the broader microbial community. 16S rRNA amplicons from the two plume samples, P16 and P52, and clustered with other plume samples (Figure 1). Sample NP52 was taxonomically and chemically distinct from plume samples, and clustered with non-plume samples, and sample IP16 had an intermediate taxonomic profile. The four dominant families indicative of the plume (see above) constituted 99% of 16S rRNA amplicons in plume samples P16 and P52, 75% in intermediate sample IP16 and 36% in non-plume sample NP52. Transcript composition The 104 million joined (that is, with overlapping paired ends), quality-filtered reads had an average length of 217 bp (2.2 1010 bp total). Internal RNA standards added at the initiation of cell lysis (467, 576, 917 and 970 bp) were used to measure the error rate due to sequencing or mis-joining of the paired-end reads, and the combined error rate was estimated at 1.40% per bp, or three errors in an average 217-bp sequence. After removal of rRNAs, tRNAs and internal standards, 66 million possible protein-encoding sequences remained. Of those, B23 million sequences were classified as one of 591 000 bacterial or archaeal genes using Blastx searches against the NCBI RefSeq database (Table 1). The rest did not have a hit in the database that met our cutoff criteria or were assigned to sequences from eukaryotes or viruses. RefSeq hits were matched to KEGG ortholog groups by Blastp, with 15 million sequences assigned to 5511 KEGG groups. In the plume transcriptomes (P16 and P52), most reads binned to six reference genomes as follows: Bermanella marisrubri RED65 (19%), Colwellia psychrerythraea 34H (16%), M. tundripaludum SV96 (5.6%), Methylococcus capsulatus str. Bath (1.8%), Neptuniibacter caesariensis MED92 (1.2%) and Methylophaga thiooxidans DMS010 (1.1%) (Supplementary Table S1). Because sequences binning to a genome were quite divergent in some cases, these estimates are only approximate. For example, the Cycloclasticus-like bacterium that was abundant in 16S rRNA amplicon libraries from the plume had no reference genome available at the time of our analysis or that of Mason et al. (2012) ; transcripts from this taxon likely binned to the genome of a Gammaproteobacterial relative, with M. tundripaludum SV96 being the closest based on 16S rRNA gene similarity. Of the dominant plume taxa, the Colwellia-like reads most closely resembled their reference genome based on protein identity (86%) and consistency of genome coverage (Supplementary Figure S5). The transcripts that recruited to other reference genomes were more distantly related (75% identity for B. marisrubri, 72% for M. tundripaludum, 69% for M. capsulatus, 67% for N. caesariensis and 68% for M. thiooxidans). The non-plume transcriptome was more diverse. Reference genomes accumulating the most transcripts were B. marisrubri (5.5%), the archaeal ammonia oxidizer Nitrosopumilus maritimus (2.9%), M. tundripaludum (2.6%) and C. psychrerythraea (1.4%), with the remaining reads broadly distributed across multiple taxa. Community metabolic reconstruction Relative abundance of transcripts matching the KEGG database pathways for the oxidation of hydrocarbons was determined (Figure 2), and statistical comparisons were made between plume and non-plume samples (R package DESeq) (Anders and Huber, 2010) . With only a single nonplume sample (NP52), the dispersion calculations for the statistic were based only on the plume samples (P16 and P52; see Materials and methods) and significance values should be interpreted as approximate. Alkanes were among the most abundant hydrocarbons measured in the Deepwater Horizon plume (Reddy et al., 2012) , and correspondingly transcripts for the enzyme alkane-1-monooxygenase, which degrades C5–C13 n-alkanes, were abundant in the transcriptomes (n ¼ 7304 total reads) (Figure 2). Monooxygenases and dioxygenases involved in aromatic compound degradation were significantly enriched in the plume transcriptomes (Po0.05), including phenol-2 monooxygenase (involved in both phenol and toluene metabolism, phn genes), cyclohexanone monooxygenase, aryl alcohol dehydrogenase, benzoate 1,2-dioxygenase and biphenyl2,3-diol 1,2, monooxygenase (Figure 2). Consistent with aromatic compounds comprising o2% of the hydrocarbons from the Deepwater Horizon reservoir fluid (Reddy et al. 2012) , the abundance of transcripts related to aromatic degradation was considerably lower than for alkane degradation. Mason et al. (2012) also reported high levels of alkane-1-monooxygenase and cyclohexanone monooxygenase in metatranscriptomes, and obtained a single-cell genome of an Oceanospirillales member with the capacity to degrade cyclohexanone. Following initial oxidation, the hydrocarbons appeared to be incorporated into biomass primarily by the glyoxylate shunt of the tricarboxylic acid (TCA) cycle (Figure 2). In this pathway, acetyl-CoA (likely originating from the oxidation of alkanes via fatty acid beta oxidation or from aromatic compounds) are processed by isocitrate lyase and malate synthase, and both of these enzymes were significantly overrepresented in the plume transcriptomes. There was also evidence that bacterioplankton were using the methylcitrate pathway to metabolize propionyl-CoA generated from the oxidation of odd chain fatty acids, shorter fatty acids or possibly aromatic compounds (Figure 2). Methane oxidation is initiated by the methane monooxygenase-mediated conversion of methane to methanol. Transcripts with homology to genes encoding methane monooxygenases were highly abundant in the transcriptomes from all samples (all subunits, n ¼ 74 405) and significantly enriched in the plume samples (Po0.05) (Figure 2). Most of these were homologous to copper-containing pmo rather than soluble methane monooxygenases (Figure 2), as was also reported by Mason et al. (2012) . The hits to pmo genes typically had low amino-acid identity with experimentally verified sequences; however, (Supplementary Figure S6), prompting construction of a phylogenetic tree with full-length monooxygenase assemblies from the GOM transcriptomes and pmoA homologs demonstrated experimentally to oxidize methane, ammonia, ethane and propane (Figure 3a). The GOM metatranscriptomic assemblies fell into four clades as follows: Group 1 contained the canonical Gammaproteobacterial methanotroph clade, Group 2 contained the Group Z sequence with unknown function (Tavormina et al., 2010) , Group 3 contained the ethane-oxidizing cluster as defined by strains ET-HIRO and ET-SHO (Suzuki et al., 2012), along with environmental ethane-assimilating clone pmo-7 (Redmond et al., 2010) , and Group X clones (Tavormina et al., 2010) , and Group 4 contained a cluster of highly divergent assemblies, most closely related to the ethane-oxidizing Group 3. Mapping of the 33 285 pmoA-like short reads from the GOM metatranscriptome indicated that 50% of reads clustered with Group 2 (39% locus 180 þ 11% locus 311), 42% with Group 3, 4% with Group 1 and 2% with Group 4. Thus reads assigned to M. tundripaludum and M. capsulatus by Blast analysis against existing genome sequences were distributed widely across the phylogenetic tree, with most resembling the Group Z-like pmoA genes of unknown function or an ethane dehydrogenase function (Group 3) (Figure 3b). The subsequent oxidation of methanol to formaldehyde is mediated by a PQQ-containing methanol dehydrogenase (Anthony, 1982) , and homologs of this gene were significantly enriched in the plume transcriptomes (Po0.05) (Figure 2). The taxonomic distribution again suggests that some genes in this category were used to oxidize methanol, but most to oxidize ethanol or larger alcohols (Supplementary Figure S7). In the case of the former, known pathways for the resultant formaldehyde include oxidation for energy or assimilation via the serine pathway (typical of alphaproteobacterial ‘type II’ and Gammaproteobacterial ‘type X’ methanotrophs) or the ribulose monophosphate pathway (RuMP; typical of Gammaproteobacteria ‘type I’ and ‘type X’ methanotrophs). In the GOM transcriptome, RuMP sequences were found in all samples, whereas three key enzymes in the serine pathway (glycerate2-kinase, hydroxypyruvate reductase and malate thiokinase) were either not found or were present in low abundance. The best-supported dissimilatory pathways for the generation of energy from methane oxidation were the dissimilatory branch of the RuMP pathway (used by most obligate type I methanotrophs) and the tetrahydrofolate pathway (used by obligate and facultative methylotrophs). Taxon-specific metabolic reconstruction The six bacterial genomes recruiting the most reads from the plume transcriptomes included three taxa specializing on non-C1 hydrocarbons (C. psychrerythraea, B. marisrubri and N. caesariensis) and three specializing on C1 compounds (M. tundripaludum, M. capsulatus and M. thiooxidans). The genome of B. marisrubri, a bacterium related to known oil degraders, recruited the largest number of transcripts (3 409 771; Supplementary Table S1). These populations expressed several genes for alkane degradation, including alkane hydroxylase (alkB), rubredoxin (alkG) and rubredoxin reductase (alkT); of those mentioned genes, alkG was significantly enriched in the plume (Supplementary Tables S2, S4). An iron-containing alcohol dehydrogenase was also overrepresented in the plume and may be responsible for subsequent oxidation steps. Twenty flagellar proteins from this genome bin were enriched in plume samples and eight were depleted, suggesting that motility was regulated in response to the hydrocarbon inputs; several chemotaxis-related proteins were also enriched (Supplementary Table S4), and this was seen by Mason et al. (2012) , as well. The genome of Gammaproteobacterium C. psychrerythraea, a psychrophilic hydrocarbon degrader (Methe et al., 2005) , recruited 3 310 791 transcripts (Supplementary Table S1). The GOM population represented by this genome bin appeared to be consuming alkanes, as putative monooxygenases and dioxygenases identified in the Colwellia genome were enriched in samples collected within the plume (Supplementary Table S4). Homologs of known hydrocarbon-degrading genes included the ringdegrading genes phnAB and their related transcriptional regulators phnS and phnR (Supplementary Table S2). The PQQ alcohol dehydrogenase of C. psychrerythraea was relatively more abundant in the plume and a number of genes encoding its cofactor were significantly overrepresented (Supplementary Table S4). The genes isocitrate lyase and malate synthase were abundant in plume samples, whereas isocitrate dehydrogenase was underrepresented, indicating that the glyoxylate shunt of the TCA cycle was highly active in the Colwellia-like populations. Genes for motility, including three flagellar and five type IV pilus genes, were significantly overrepresented in the plume, as was the ferrous iron transporter feoB. Multiple subunits of cytochrome cbb(3) were enriched in the plume, a cytochrome with high affinity for oxygen that is typically expressed in suboxic conditions (Pitcher and Watmough, 2004) . In the N. caesariensis transcriptome bin (2 98 158 reads), homologs of the ring-hydroxylating dioxygenase phnAb, phenol-1-monooxygenase and the permease component of an ABC transporter for toluene suggested populations involved in the consumption of aromatic hydrocarbons (Supplementary Table S4). Quinoprotein-containing alcohol dehydrogenase transcripts were well expressed, as were glyoxylate shunt genes isocitrate lyase and malate synthase. One of the most abundant transcripts was for poly(R)-hydroxyalkanoic acid synthase, implying that Neptuniibacterlike populations were storing carbon derived from the hydrocarbons. Among the C1-utilizing taxa, the largest number of transcripts was recruited to M. tundripaludum (1 408 081 reads; Supplementary Table S4). At least one of the populations recruiting to this genome was living at the expense of methane, based on the presence of transcripts clustering with the verified pmo genes (pmoCAB), although the dominance of transcripts for ethane monooxygenases and divergent particulate monooxygenases (see above) suggests that this reference genome captured a number of populations with varying substrate specialization. Methanol dehydrogenase sequences were evident (mxaFI), although read counts were two orders of magnitude lower than for pmo, consistent with the possibility that most pmo-like genes were degrading other short chain alkanes. All three pathways for oxidizing formaldehyde predicted in this genome (Svenning et al., 2011) were represented in the transcriptome (the pterin-dependent pathways involving H4folate and H4MPT, as well as formaldehyde dehydrogenase), along with a formaldehydeactivating enzyme that accelerates the condensation of formaldehyde with H4MPT (Trotsenko and Murrell, 2008) . Transcripts for carbon assimilation via the RuMP pathway were evident, as well as for all the steps in the TCA cycle. Although the M. tundripaludum reference genome contains genes for nitrogen fixation and denitrification, these were not represented in the transcripts of the GOM populations. Genes for nitrite transport and assimilation were represented, including nitrite reductase (nirBD). Overall, the plume populations contained relatively more transcripts for flagellum synthesis, chemotaxis, pilus assembly, and cobalamin and pterin synthesis than the non-plume sample (Supplementary Table S3). The C1-utilizing genome of M. capsulatus recruited 4 53 969 reads from the GOM transcriptomes. The reads assigned to the two pmo operons (pmoCAB; 53 947 reads) also had greater similarity to non-methane oxidizing enzymes or those with unknown function (see above and Supplementary Table S4). There was lower but measurable expression of the soluble methane monooxygenase operon (378 reads; mmoX,Y,Z,B,C,D). Similar to M. tundripaludum, the M. capsulatus genome bin contained transcripts for all three types of formaldehyde oxidation in approximately equal numbers. M. capsulatus is categorized as a type X methanotroph based on evidence for low levels of enzymes from the serine pathway and Calvin cycle that may supplement carbon assimilation by the RuMP pathway (Hanson and Hanson, 1996) , and transcript patterns in the GOM populations are consistent with this assignment (Supplementary Table S3). Whether or not M. capsulatus has a functional TCA cycle has been a point of controversy (Ward et al., 2004; Kelly et al., 2005) , but the GOM populations were indeed operating the TCA cycle based on significant expression of the 2-oxoglutarate dehydrogenase E1 and E2 subunits. Nitrite reductase transcripts suggest that nitrite was a major source of N; the nitrate transporter was not expressed and the ammonium transporter had low relative expression. Plume populations had greater relative expression of a fatty acid degradation gene (long chain fatty acid CoA ligase) (Supplementary Table S3). The C1-utilizing M. thiooxidans (recruiting 266 672 reads) is capable of oxidizing methanol but not methane (Boden et al., 2011) . Transcripts for methanol dehydrogenase were present in this genome bin and the two formaldehyde oxidation pathways involving pterins (H4folate and H4MPT) were being transcribed. A surprisingly large number of genes encoding extracellular processes were among the most highly expressed genes in the M. thiooxidans genome bin, including type I, II and III secretion systems, a type IV pilus system, signal transduction and second messenger molecules, chemotaxis genes and LPS and lipid A synthesis. Transcripts for transport of both nitrite and ammonia were abundant, as were those for nitrite reductase. Ecology of the plume community The transcription shifts induced in the GOM bathypelagic bacterial community by the Deepwater Horizon hydrocarbons could represent either of the two ecological scenarios: (1) blooming hydrocarbonutilizing taxa that competed directly (for example, for nutrients or substrates) with the pre-existing bacterioplankton community, and whose rapid growth displaced typical GOM bathypelagic bacteria; or (2) hydrocarbon-utilizing taxa that responded to the new substrates without displacing the pre-existing species, with the latter representing a smaller proportion of the total transcriptome yet not adversely affected in numbers or activity. To address this issue, we calculated the numbers of transcripts l 1 for each genome bin in each sample, based on recovery of the internal standards added during sample processing (see Materials and Methods). The 50 best-represented taxa in the NP transcriptome fell into four qualitative ecological categories based on abundance patterns (Figure 4). Approximately 1/3 of the taxa had significantly higher numbers of transcripts in the plume by up to two orders of magnitude (one-way analysis of variance, Po0.001); this category consisted entirely of Gammaproteobacteria and included all six dominant bloom taxa (Bermanella, Colwellia, Neptuniibacter, Methylobacter, Methylococcus and Methylophaga) along with several members of the Oligotrophic Marine Gammaproteobacteria group (HTCC2207, HTCC2143, HTCC2148 and NOR5-3) and taxa affiliated with the Chromatiales (Figure 4). An additional B1/4 of the taxa similarly had increased contribution to the plume transcriptome, but only by two- to fivefold; these included Actinobacteria, marine Rhodospirallales relatives, Planctomyces and groups related to sulfur-oxidizing Gammaproteobacterial symbionts. Transcript numbers for another 1/3 of the taxa, including SAR11 and marine Thaumarchaeota, were not detectibly affected by the presence of the fast growing Gammaproteobacteria in the plume (analysis of variance, P ¼ 0.53). Finally, a surprisingly low 5% of the taxa, including two binning to Firmicutes genomes, had lower numbers of transcripts in the plume (Figure 4), apparently adversely affected by either competition or hydrocarbon toxicity. Discussion Microbial community composition and gene expression were profoundly altered by petroleum contamination, even in samples where plume signals were not detectible with in situ instrumentation. In NP52, CDOM fluorescence was below the 1 p.p.m. limit of detection for crude oil, yet 36% of 16S rRNA sequences and 17% of transcripts were from the four families that dominated plumeassociated communities. These same organisms made up o4% of a nearby bathypelagic community before the spill (King et al., 2012) . Methane was indeed elevated in the control sample compared with historical background concentrations (106 nM in NP52 versus B2 nM background; Lamontagne et al., 1973) . Yet with concentrations still B2000-fold less than the in-plume samples, the microbial community response was quite strong. Because sample NP52 contained low-level hydrocarbon contamination, differences in relative expression detected between impacted and non-impacted samples are likely to underestimate the effects of hydrocarbon inputs on bacterial gene expression. mRNAs from the four metatranscriptomes mapped to over 2700 bacterial and archaeal taxa, but the similarity of reads to their reference sequences varied considerably. Of the six highestrecruiting genome bins, only the Colwellia-like sequences closely resembled the reference genome (87% average protein identity; Supplementary Figure S5), whereas fragment recruitment plots for the other five taxa show genome-wide protein identities of B70%. The low sequence similarity to reference genomes affects interpretation of the metatranscriptomic data, as physiologies of the captured populations may vary considerably among each other, as well as compared with the sequenced representative. One relevant example is that reads mapping to the pmoCAB genes of M. tundripaludum and M. capsulatus represent multiple clusters with varying percent identities and likely originating from several populations (Supplementary Figure S6). The evolution of bacterial monooxygenases is complex: pmoA genes are paraphyletic; the superfamily encodes enzymes targeting methane, ammonia, ethane and propane; and several newly discovered homologs have unknown function, including Group Z pmoA and pxmA (Tavormina et al., 2011) . Of the 33 000 reads annotated as pmoA based on RefSeq Blast analysis, only 4% clustered with experimentally verified methane oxidation genes, 50% clustered with a Group 2 (Group Z-like) pmoA of unknown function and 42% clustered with ET-SHO/ ET-HIRO-like short chain hydrocarbon monooxygenases (Figure 3). Approximately 2% of the pmoA homologs were from a new, highly divergent group that has not been described previously (Group 4; Figure 3). Thus, the metatranscriptomic data cannot definitively address the question of the importance of methane oxidation in the Deepwater Horizon plume until the diverse pmo-like sequences are better characterized functionally. Likewise, only 3% of the 36 000 PQQ and PQQheme dehydrogenase sequences that catalyze the next step in methane oxidation grouped with verified methanol dehydrogenases. Transcripts representing the alkane-1-monooxygenase (alkB), rubredoxin and rubredoxin reductase system responsible for degradation of C6 to C13 n-alkanes were also abundant and enriched in the plume, whereas transcripts for oxidation of medium chain (C10–C30) alkanes (alkM) were particularly abundant in the C. psychrerythraea and B. marisrubri genome bins. With the exception of C1 metabolism, transcripts counts for alkane degradation were negatively related to hydrocarbon length and positively related to hydrocarbon solubility and concentration. Nitrate and nitrite transporters were abundant in the plume, and nitrite concentrations were elevated (Table 1). Concentrations of nitrate were generally lower in plume samples P16 and P52 and other hydrocarbon plume samples collected from the spill site (data not shown) compared with non-impacted samples, suggesting that dissimilatory nitrate reduction was occurring. Partial nitrification by ammonia oxidizers (such as N. maritimus) could also have been producing nitrite. The most abundant organism bins had enriched nitrite transporter transcripts and may have been taking advantage of the nitrite as a nutrient source. Nitrogen and phosphate did not appear to be limiting growth at the time of sampling based on measured concentrations (Table 1). The natural seepage of hydrocarbons from the GOM seafloor (4–10 1010 g per year; National Research Council, 2003) may set the context for the observed ecological response of the bathypelagic microbial community by sustaining a native community of hydrocarbon-degrading bacteria that quickly responds to periodic fluxes of oil and gas. In support of this idea, most of the pmoA homologs represented in the plume and non-plume transcriptomes were similar to the monooxygenases observed at natural hydrocarbon seep sites and submarine canyons (Redmond et al., 2010; Tavormina et al., 2010; Lesniewski et al., 2012; Suzuki et al., 2012; Li et al., 2013) . Further, populations related to several of the bloom-forming taxa (Cycloclasticus and Methylococcus) were present in bathypelagic GOM seawater 1 month before the Deepwater Horizon event, although at low abundance (King et al., 2012) . On the other hand, the presence of taxa whose transcript abundance did not vary between plume and unexposed samples suggests the presence of a community of non-responding bacterioplankton with a potential role in community structure recovery. Expression patterns of representatives of these non-responding microbes showed fewer differences in gene transcription patterns between control and plume samples than bloom taxa (Figure 5), and may have continued to occupy prespill niches in the petroleum-impacted seawater. This metatranscriptomic analysis has revealed the primary pathways used for hydrocarbon degradation following the Deepwater Horizon accident, based on highly abundant transcripts for degradation of C2–C4 gaseous alkanes and C5–C12 alkanes, and less abundant but still considerable transcripts for degradation of C12–C18 alkanes, aromatic compounds and methane. Many of the monooxygenase sequences binned to methane monooxygenases from sequenced bacteria, but are more likely to be involved in short chain alkane oxidation based on recent studies from deep-sea seeps. The hydrocarbons were metabolized by a few key pathways, predominantly beta oxidation, glyoxylate, 2-methylcitrate and RuMP. A few groups of Gammaproteobacteria accounted for most of the bacterial response. Despite the enormous bloom of hydrocarbon-degrading Gammaproteobacteria that increased bacterial cell counts by two orders of magnitude, members of the natural microbial community persisted at their pre-bloom activity levels and may be important in the re-establishment of the original microbial community. Conflict of Interest The authors declare no conflict of interest. Acknowledgements We are grateful to B Satinsky for preparation of the internal standards, H Luo for phylogenetic support, Y Sun for help with Circos, and S Malfatti and E Lindquist for sequencing assistance. This work was supported by the Gordon and Betty Moore Foundation, US National Science Foundation grant OCE-1043225, and resources and technical expertise from the University of Georgia’s Georgia Advanced Computing Resource Center. The work conducted by the US Department of Energy (DOE) Joint Genome Institute is supported by the Office of Science, Biological and Environmental Research Program of the DOE under contract number DE-AC02-05CH11231. The views and opinions of the authors expressed herein do not necessarily state or reflect those of the United States government, or any agency thereof, or the Regents of the University of California. Supplementary Information accompanies this paper on The ISME Journal website ( Anders S , Huber W. ( 2010 ). Differential expression analysis for sequence count data . Genome Biol 11 : R106 . Anthony C. ( 1982 ). The biochemistry of methylotrophs . Academic Press: London, UK. Benjamini Y , Hochberg Y. ( 1995 ). Controlling the false discovery rate - a practical and powerful approach to multiple testing . J R Stat Soc Series B Stat Methodol 57 : 289 - 300 . Boden R , Ferriera S , Johnson J , Kelly DP , Murrell JC , Scha¨fer H. ( 2011 ). Draft genome sequence of the chemolithoheterotrophic, halophilic methylotroph Methylophaga thiooxydans DMS010 . J Bacteriol 193 : 3154 - 3155 . Camilli R , Reddy CM , Yoerger DR , Van Mooy BAS , Jakuba MV , Kinsey JC et al. ( 2010 ). Tracking hydrocarbon plume transport and biodegradation at Deepwater Horizon . Science 330 : 201 - 204 . Capella-Gutierrez S , Silla-Martinez JM , Gabaldon T. ( 2009 ). trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses Bioinformatics 25 : 1972 - 1973 . Caporaso JG , Kuczynski J , Stombaugh J , Bittinger K , Bushman FD , Costello EK et al. ( 2010 ). QIIME allows analysis of high-throughput community sequencing data . Nat Methods 7 : 335 - 336 . Caspi R , Altman T , Dreher K , Fulcher CA , Subhraveti P , Keseler IM et al. ( 2012 ). The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases . Nucleic Acids Res 40 : D742 - D753 . DeSantis TZ , Hugenholtz P , Larsen N , Rojas M , Brodie EL , Keller K et al. ( 2006 ). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB . Appl Environ Microbiol 72 : 5069 - 5072 . Eddy SR. ( 2011 ). Accelerated Profile HMM Searches . PLoS Comput Biol 7 : e1002195 . Falgueras J , Lara A. ( 2010 ). SeqTrim: a high-throughput pipeline for pre-processing any type of sequence read . BMC Bioinformatics 11 : 38 . Gifford SM , Sharma S , Rinta-Kanto JM , Moran MA . ( 2011 ). Quantitative analysis of a deeply sequenced marine microbial metatranscriptome . ISME J 5 : 461 - 472 . Han M , Zmasek C. ( 2009 ). phyloXML: XML for evolutionary biology and comparative genomics . BMC Bioinformatics 10 : 356 . Hanson RS , Hanson TE . ( 1996 ). Methanotrophic bacteria . Microbiol Rev 60 : 439 - 471 . Hazen TC , Dubinsky EA , DeSantis TZ , Andersen GL , Piceno YM , Singh N et al. ( 2010 ). Deep-sea oil plume enriches indigenous oil-degrading bacteria . Science 330 : 204 - 208 . Jari O , GB F , Roeland K , Pierre L , Peter RM , BO R et al. ( 2013 ). Vegan: Community Ecology Package. Joye SB , MacDonald IR , Leifer I , Asper V. ( 2011 ). Magnitude and oxidation potential of hydrocarbon gases released from the BP oil well blowout . Nat Geosci 4 : 160 - 164 . Kanehisa M , Goto S. ( 2000 ). KEGG: Kyoto encyclopedia of genes and genomes . Nucleic Acids Res 28 : 27 - 30 . Katoh K , Misawa K , Kuma K , Miyata T. ( 2002 ). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform . Nucleic Acids Res 30 : 3059 - 3066 . Kelley D. ( 2013 ). Oce: Analysis of Oceanographic data . Kelly DP , Anthony C , Murrell JC. ( 2005 ). Insights into the obligate methanotroph Methylococcus capsulatus . Trends Microbiol 13 : 195 - 198 . Kessler JD , Valentine DL , Redmond MC , Du M , Chan EW , Mendes SD et al. ( 2011 ). A persistent oxygen anomaly reveals the fate of spilled methane in the deep Gulf of Mexico . Science 331 : 312 - 315 . King GM , Smith CB , Tolar B , Hollibaugh JT ( 2012 ). Analysis of composition and structure of coastal to mesopelagic bacterioplankton communities in the northern Gulf of Mexico . Front Microbio 3 : 438 . Krzywinski M , Schein J , Birol I , Connors J , Gascoyne R , Horsman D et al. ( 2009 ). Circos: an information aesthetic for comparative genomics . Genome Res 19 : 1639 - 1645 . Lamontagne RA , Swinnerton JW , Linnenbom VJ , Smith WD . ( 1973 ). Methane concentrations in various marine environments . J Geophys Res 78 : 5317 - 5324 . Le SQ , Lartillot N , Gascuel O. ( 2008 ). Phylogenetic mixture models for proteins . Philos Trans R Soc Lond B Biol Sci 363 : 3965 - 3976 . Lesniewski RA , Jain S , Anantharaman K , Schloss PD , Dick GJ . ( 2012 ). The metatranscriptome of a deep-sea hydrothermal plume is dominated by water column methanotrophs and lithotrophs . ISME J 6 : 2257 - 2268 . Li M , Jain S , Baker BJ , Taylor C , Dick GJ . ( 2013 ). Novel hydrocarbon monooxygenase genes in the metatranscriptome of a natural deep-sea hydrocarbon plume . Environ Microbiol e-pub ahead of print Doi: 10 .1111/ 1462 - 2920 . 12182 . Lozupone C , Lladser ME , Knights D , Stombaugh J , Knight R. ( 2011 ). UniFrac: an effective distance metric for microbial community comparison . ISME J 5 : 169 - 172 . Lu Z , Deng Y , Van Nostrand JD , He Z , Voordeckers J , Zhou A et al. ( 2011 ). Microbial gene functions enriched in the Deepwater Horizon deep-sea oil plume . ISME J 6 : 451 - 460 . Martin J , Bruno VM , Fang Z , Meng X , Blow M , Zhang T et al. ( 2010 ). Rnnotator: an automated de novo transcriptome assembly pipeline from stranded RNA-Seq reads . BMC Genomics 11 : 663 . Mason OU , Hazen TC , Borglin S , Chain PSG , Dubinsky EA , Fortney JL et al. ( 2012 ). Metagenome, metatranscriptome and single-cell sequencing reveal microbial response to Deepwater Horizon oil spill . ISME J 6 : 1715 - 1727 . Matsen FA , Kodner RB , Armbrust EV . ( 2010 ). Pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree . BMC Bioinformatics 11 : 538 . Methe BA , Nelson KE , Deming JW , Momen B , Melamud E , Zhang XJ et al. ( 2005 ). The psychrophilic lifestyle as revealed by the genome sequence of Colwellia psychrerythraea 34H through genomic and proteomic analyses . Proc Natl Acad Sci USA 102 : 10913 - 10918 . National Research Council. ( 2003 ). Oil in the Sea III: Inputs, Fates, and Effects . National Academies Press: Washington, DC. Pitcher RS , Watmough NJ . ( 2004 ). The bacterial cytochrome cbb3 oxidases . Biochim Biophys Acta 1655 : 388 - 399 . Poretsky RS , Gifford S , Rinta-Kanto J , Vila-Costa M , Moran MA . ( 2009 ). Analyzing gene expression from marine microbial communities using environmental transcriptomics . J Vis Exp 3 : 1 - 6 . Reddy CM , Arey JS , Seewald JS , Sylva SP , Lemkau KL , Nelson RK et al. ( 2012 ). Composition and fate of gas and oil released to the water column during the Deepwater Horizon oil spill . Proc Natl Acad Sci USA 109 : 20229 - 20234 . Redmond MC , Valentine DL . ( 2012 ). Natural gas and temperature structured a microbial community response to the Deepwater Horizon oil spill . Proc Natl Acad Sci USA 109 : 20292 - 20297 . Redmond MC , Valentine DL , Sessions AL . ( 2010 ). Identification of novel methane- , ethane- , and propane-oxidizing bacteria at marine hydrocarbon seeps by stable isotope probing . Appl Environ Microbiol 76 : 6412 - 6422 . Reeder J , Knight R. ( 2010 ). Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions . Nat Methods 7 : 668 - 669 . Rodrigue S , Materna AC , Timberlake SC , Blackburn MC , Malmstrom RR , Alm EJ et al. ( 2010 ). Unlocking short read sequencing for metagenomics . PLoS ONE 5 : e11840 . Ryerson TB , Camilli R , Kessler JD , Kujawinski EB , Reddy CM , Valentine DL et al. ( 2012 ). Chemical data quantify Deepwater Horizon hydrocarbon flow rate and environmental distribution . Proc Natl Acad Sci USA 109 : 20246 - 20253 . Satinsky BM , Gifford SM , Crump BC , Moran MA . ( 2013 ). Use of internal standards for quantitative metatranscriptome and metagenome analysis . Meth Enzym (in press). Stamatakis A. ( 2006 ). RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models . Bioinformatics 22 : 2688 - 2690 . Stewart FJ , Ottesen EA , DeLong EF . ( 2010 ). Development and quantitative analyses of a universal rRNAsubtraction protocol for microbial metatranscriptomics . ISME J 4 : 896 - 907 . Suzuki T , Nakamura T , Fuse H. ( 2012 ). Isolation of two novel marine ethylene-assimilating bacteria, Haliea species ETY-M and ETY-NAG, containing particulate methane monooxygenase-like genes . Microbes Environ 27 : 54 - 60 . Svenning MM , Hestnes AG , Wartiainen I , Stein LY , Klotz MG , Kalyuzhnaya MG et al. ( 2011 ). Genome sequence of the Arctic methanotroph Methylobacter tundripaludum SV96 . J Bacteriol 193 : 6418 - 6419 . Tatusov RL , Fedorova ND , Jackson JD , Jacobs AR , Kiryutin B , Koonin EV et al. ( 2003 ). The COG database: an updated version includes eukaryotes . BMC Bioinformatics 4 : 41 . Tavormina PL , Orphan VJ , Kalyuzhnaya MG , Jetten MSM , Klotz MG . ( 2011 ). A novel family of functional operons encoding methane/ammonia monooxygenase-related proteins in gammaproteobacterial methanotrophs . Environ Microbiol Rep 3 : 91 - 100 . Tavormina PL , Ussler W , Joye SB , Harrison BK , Orphan VJ . ( 2010 ). Distributions of putative aerobic methanotrophs in diverse pelagic marine environments ISME J 4 : 700 - 710 . Trotsenko YA , Murrell JC. ( 2008 ). Metabolic aspects of aerobic obligate methanotrophy . Adv Appl Microbiol 63 : 183 - 229 . Valentine DL , Kessler JD , Redmond MC , Mendes SD , Heintz MB , Farwell C et al. ( 2010 ). Propane respiration jump-starts microbial response to a deep oil spill . Science 330 : 208 - 211 . Wang Q , Garrity GM , Tiedje JM , Cole JR . ( 2007 ). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy . Appl Environ Microbiol 73 : 5261 - 5267 . Ward N , Larsen O , Sakwa J , Bruseth L , Khouri H , Durkin AS et al. ( 2004 ). Genomic insights into methanotrophy: the complete genome sequence of Methylococcus capsulatus (Bath) . PLoS Biol 2 : e303 . Yvon-Lewis SA , Hu L , Kessler J . ( 2011 ). Methane flux to the atmosphere from the Deepwater Horizon oil disaster . Geophys Res Lett 38 : L01602 . Zhu W , Lomsadze A , Borodovsky M. ( 2010 ). Ab initio gene identification in metagenomic sequences . Nucleic Acids Res 38 : e132 .

This is a preview of a remote PDF:

Adam R Rivers, Shalabh Sharma, Susannah G Tringe, Jeffrey Martin, Samantha B Joye, Mary Ann Moran. Transcriptional response of bathypelagic marine bacterioplankton to the Deepwater Horizon oil spill, ISME Journal, 2013, 2315-2329, DOI: 10.1038/ismej.2013.129