- Open Access
Proliferation of hydrocarbon-degrading microbes at the bottom of the Mariana Trench
- Jiwen Liu†1, 2,
- Yanfen Zheng†1,
- Heyu Lin1,
- Xuchen Wang3,
- Meng Li4,
- Yang Liu4,
- Meng Yu3,
- Meixun Zhao2, 3,
- Nikolai Pedentchouk5,
- David J. Lea-Smith6,
- Jonathan D. Todd6,
- Clayton R. Magill7,
- Wei-Jia Zhang8,
- Shun Zhou1,
- Delei Song1,
- Haohui Zhong1,
- Yu Xin2, 3,
- Min Yu1, 2,
- Jiwei Tian9, 10Email author and
- Xiao-Hua Zhang1, 2Email author
© The Author(s). 2019
- Received: 25 November 2018
- Accepted: 22 February 2019
- Published: 12 April 2019
The Mariana Trench is the deepest known site in the Earth’s oceans, reaching a depth of ~ 11,000 m at the Challenger Deep. Recent studies reveal that hadal waters harbor distinctive microbial planktonic communities. However, the genetic potential of microbial communities within the hadal zone is poorly understood.
Here, implementing both culture-dependent and culture-independent methods, we perform extensive analysis of microbial populations and their genetic potential at different depths in the Mariana Trench. Unexpectedly, we observed an abrupt increase in the abundance of hydrocarbon-degrading bacteria at depths > 10,400 m in the Challenger Deep. Indeed, the proportion of hydrocarbon-degrading bacteria at > 10,400 m is the highest observed in any natural environment on Earth. These bacteria were mainly Oleibacter, Thalassolituus, and Alcanivorax genera, all of which include species known to consume aliphatic hydrocarbons. This community shift towards hydrocarbon degraders was accompanied by increased abundance and transcription of genes involved in alkane degradation. Correspondingly, three Alcanivorax species that were isolated from 10,400 m water supplemented with hexadecane were able to efficiently degrade n-alkanes under conditions simulating the deep sea, as did a reference Oleibacter strain cultured at atmospheric pressure. Abundant n-alkanes were observed in sinking particles at 2000, 4000, and 6000 m (averaged 23.5 μg/gdw) and hadal surface sediments at depths of 10,908, 10,909, and 10,911 m (averaged 2.3 μg/gdw). The δ2H values of n-C16/18 alkanes that dominated surface sediments at near 11,000-m depths ranged from − 79 to − 93‰, suggesting that these sedimentary alkanes may have been derived from an unknown heterotrophic source.
These results reveal that hydrocarbon-degrading microorganisms are present in great abundance in the deepest seawater on Earth and shed a new light on potential biological processes in this extreme environment.
- Challenger Deep
- Mariana Trench
- Hadal water
- Microbial community
- Hydrocarbon degradation
- Hydrocarbon biosynthesis
The hadal zone (6000–11,000 m) is the deepest part of the marine environment and is composed almost exclusively of trenches. Hadal zone planktonic microbial communities have recently been investigated in the Mariana [1–3], Puerto Rico , Japan , and Kermadec Trenches  due to rapid technological progress in deep-sea sampling. The hadal and abyssal (4000–6000 m) waters differ significantly in their micro- and macro-community compositions [1, 5, 7]. Specifically, heterotrophic Gammaproteobacteria and Bacteroidetes are dominant in the hadal ocean, compared to Thaumarchaeota in the abyssal ocean [1, 5, 6]. There are only a few depth profile studies of microbial communities in the hadal zone, and thus, little is known of the variability in community composition in this unique environment. In the Challenger Deep of the Mariana Trench, Nunoura et al. demonstrated via 16S rRNA gene sequencing that Gammaproteobacteria were the dominant class at ~ 10,257 m, comprising 62.7% of the total microbial population, compared to 21.6% and 1.7% at 9000 and 7998 m, respectively . Moreover, Oceanospirillales comprised 78.1% of the total Gammaproteobacteria population at 10,257 m, unlike at 9000 m where Pseudomonadales dominated. Tarn et al. also confirmed the dominance of Gammaproteobacteria in both the Challenger (~ 10,918 m) and Sirena (~ 10,667 m) Deeps of the Mariana Trench, using 16S rRNA gene sequencing . The drivers for this major microbial community shift at depths > 10,000 m are undetermined, since studies analyzing the genetic potential of microbes  and their abundance at the genus and species level at these depths are limited.
To address these knowledge gaps, we conducted the most detailed analysis to date of microbial populations and functional capacities in the Mariana Trench water samples from different depths by metagenome sequencing, 16S rRNA gene amplicon, and gene probe analyses. An unexpected enrichment for hydrocarbon-degrading bacteria, along with an increase in the abundance and transcription of hydrocarbon-degrading genes, was found at depths near the trench bottom (10,400 m and 10,500 m). Based on this information, we experimentally tested and confirmed the hydrocarbon degradation activity of bacterial strains isolated from 10,400 m water supplemented with hexadecane under conditions simulating the deep ocean. Stable hydrogen isotope analysis of surface sediments at the bottom of the trench suggests that n-C16/18 alkanes, the dominant hydrocarbons in these samples, are derived from an unknown biological source, and may support hydrocarbon-degrading microbial populations.
Microbial communities vary throughout the Mariana Trench water column
To determine the vertical distribution of microbial populations in the Challenger Deep, shotgun metagenome sequencing and 16S rRNA gene amplicon analyses were performed on free-living (0.22–3 μm) (FL) and particle-associated (> 3 μm) (PA) communities collected from surface waters and from depths of 4000 m, 9600 m, 10,400 m, and 10,500 m (Additional file 1: Table S1 and Figure S1). Also, environmental analysis was performed on samples collected from 4 to 8320 m, which demonstrated that hadal and abyssal waters shared many comparable physiochemical factors, including temperature, salinity, oxygen, and nutrients (Additional file 1: Table S2).
Quantitative PCR assays showed that the 16S rRNA gene copy numbers decreased significantly from 0 (7.73 × 107 copies/L) to 4000 m (1.21 × 107 copies/L), but then increased with depth to 7.43 × 107, 7.58 × 107, and 7.82 × 107 copies/L at 9600, 10,400, and 10,500 m, respectively (Additional file 1: Figure S2). Similarly, Nunoura et al. showed that the abundance of 16S rRNA gene decreased along with depth before 6000 m but then increased slightly from 6000 m to 10,257 m except for 9000 m . In this study, the increase of 16S rRNA gene copy number in trench bottom water compared to 4000 m may be due to increased Gammaproteobacteria populations (see below), which harbor ~5.7 copies of 16S rRNA gene on average, far higher than other Proteobacteria , suggesting that 16S rRNA gene copy number is not always a good proxy of cell number.
Next, we analyzed the genus-level composition to identify the genera that are significantly enriched in NBW samples. According to the metagenomic data annotated against the NCBI-nr database, seven genera were significantly more abundant in NBW than at other depths (P < 0.05) (Fig. 1b). All belonged to the order Oceanospirillales and included Oleibacter (19.5 ± 6.5%), Thalassolituus (2.5 ± 0.7%), Alcanivorax (1.0 ± 0.4%), Oleiphilus (0.9 ± 0.7%), Oceanobacter (0.6 ± 0.2%), Hahella (0.5 ± 0.4%), and Bermanella (0.3 ± 0.1%). 16S rRNA gene amplicon analysis also confirmed the predominance of Oleibacter in NBW (37.8 ± 8.2%). This is the first environment where Oleibacter has been reported as the most abundant genus.
Oceanospirillales species encode genes involved in alkane degradation
To the best of our knowledge, the proportion of Oceanospirillales-derived alkane degraders in NBW (~ 20% in the metagenome and ~ 38% in the 16S rRNA gene amplicon; Additional file 1: Figure S4a and c) is higher than in any other natural environments observed on Earth. However, the predominance of Oceanospirillales, along with MCP, alkB, and type IV pilus assembly genes in NBW, resembles the situation (dominant microbial communities and functional genes) in the Deepwater Horizon spill (DWH) [15–19], suggesting that the microorganisms in these two distinct environments most likely share similar hydrocarbon-degrading pathways. However, the predominant DWH Oceanospirillales showed high sequence similarity to Candidatus Bermanella macondoprimitus and/or Candidatus Oceanospirillales desum [20, 21]. Although O. marinus was detected in n-hexadecane enrichments using water collected from oil plume at 1000 m after the DWH, it only represented 3% of the total sequences .
Abundance and expression of alkane degradation genes increase in near bottom waters
To determine whether increases in the transcription of alkB and almA corelate with their gene abundance, we performed RT-qPCR on 0, 4000, 9600, 10,400, and 10,500 m samples. AlkB and almA genes were highly transcribed in NBW compared to 0, 4000, and 9600 m (up to 42, 61, 26 and 86, 109, 34-fold increased, respectively) (Fig. 4c), and 81.5–100% of transcripts were from Oleibacter. This further supports the hypothesis that Oceanospirillales species are degrading both medium- and long-chain alkanes in this environment.
Near bottom water isolates degrade alkanes under physiologically relevant conditions
O. marinus DSM 24913 was tested for hydrocarbon degradation activity, since no Oleibacter species was isolated from NBW (despite repeated attempts with alkane-dependent media and Marine Broth under 4/16 °C and 0.1/60 MPa) and O. marinus is the only cultured species available in this genus. Oleibacter MAG 01, the most abundant clade in NBW, demonstrated 99.0% 16S rRNA gene sequence identity to O. marinus DSM 24913. Moreover, MAG 01 contained a putative alkB and two almA genes, which encoded enzymes with > 99.5% amino acid identity to those of O. marinus DSM 24913 , and up to 34.0% and 49.7% identity to functionally characterized AlkB and AlmA from A. dieselolei, respectively [14, 23]. MAG 01 also harbored 11 putative MCP genes (O. marinus DSM 24913 possessed 14 MCP genes), with up to 43.9% amino acid identity to the functionally confirmed A. dieselolei MCP involved in sensing alkanes .
Unlike NBW Alcanivorax species, O. marinus DSM 24913 could not degrade hydrocarbons at 4 °C and 60 MPa, but it displayed a high rate of hydrocarbon degradation in cultures fed on n-eicosane at 16 °C and atmospheric pressure (Fig. 5b). This implies that NBW Oleibacter species, such as MAG 01, have evolved to withstand high pressure and low temperature typical of hadal conditions. Thus, we propose that the NBW Oleibacter species may represent novel and hadal-specific ecotypes despite their high 16S rRNA gene (i.e., 99.0% identity) and genome sequence (98.6% average nucleotide identity) identity to O. marinus DSM 24913. A similar phenomenon was shown for Alteromonas macleodii species, which are phylogenetically divided into “surface” or “deep” ecotypes (recently a new species A. mediterranea was separated ), according to seawater depth, despite sharing > 99% 16S rRNA gene similarity . In support of bioinformatic predictions, n-alkanes can be utilized as a carbon and/or energy source by all tested NBW Oceanospirillales isolates under high pressure and low temperature, and O. marinus, representative of the dominant Oleibacter genus in NBW under atmospheric pressure at 16 °C.
Potential alkane accumulation at deeper depths and production in sediments in the Mariana Trench
Given the NBW Alcanivorax isolates were more efficient at degrading C18–C20 compared to > C22 alkanes (Fig. 5a), it is possible that these bacteria have evolved to take advantage of these specific hydrocarbons. Thus, this unknown biological source of n-C16, n-C18, and n-C20 hydrocarbons in the surface sediments at the bottom of the trench may contribute, with sinking particle hydrocarbons, to support the observed NBW hydrocarbon-degrading microbial population. Further work, involving turnover rate measurements, is likely required to establish this hypothesis. A similar hydrocarbon cycle is likely present in surface waters, where hydrocarbon-degrading bacteria are supported by alkanes synthesized by cyanobacteria and algae . Currently, it is uncertain if these hydrocarbons are produced by microbes only found in the sediments or whether they are widespread in NBW, since we were unable to study their prevalence in NBW samples and the enzymes synthesizing even-chain alkanes have not been identified.
In summary, this study demonstrates the highly stratified metabolic potential of the microbial community in the hadal zone of the Mariana Trench. The deepest waters (> 10,000 m) are enriched with alkane-degrading bacteria of Oceanospirillales at a proportion higher than observed anywhere else on Earth, which are transcribing alkane-degrading genes. The study raises important questions that warrant further investigation in this unique environment. Which sources of alkanes support microorganisms in the deepest waters on Earth: novel biological alkanes identified in deep hadal sediments or alkanes derived from other sources (e.g., surface water photoautotrophs and/or petroleum), or both? Further, what causes the abrupt and unexpected shift in metabolic potential towards hydrocarbon degradation at depths below 10,000 m and are such shifts common in other deep trench environments?
Seawater at 0 (cast #101, 200 L), 4000 (cast #102, 200 L), 9600 (cast #103, 96 L), 10,400 (cast #104, 96 L), and 10,500 m (cast #105, 96 L) was collected using Niskin bottles from the Challenger Deep of the Mariana Trench (11° 20.605′ N, 142° 19.557′ E, 10,800 m) aboard the R/V Dong Fang Hong 2 in September 2016. Collection points are shown in Additional file 1: Figure S1 . Water salinity and temperature (Additional file 1: Table S2) were recorded by CTD. For each depth, water samples from three or four Niskin bottles were pooled and filtered onto one membrane. It took 4–5 h to retrieve seawater from the deepest depth (10,500 m). Upon resurfacing, biomass collection was completed within 1 h. Seawater was filtered serially through 3 μm (TSTP, 142 mm, Millipore) and 0.22 μm (GTTP, 142 mm, Millipore) polycarbonate membranes. The communities collected on the 3-μm and 0.22-μm filters were designated as particle-associated and free-living fractions, respectively. All filters were stored in liquid nitrogen onboard and at − 80 °C in the laboratory. Sinking particles were collected at the same location as seawater using time series sediment traps (McLane Research Laboratories, Mark-7) deployed at three depths (2000 m, 4000 m, and 6000 m) from January to September 2016. The surface sediment samples were collected from the Mariana Trench (11° 19.5132′ N, 142° 11.2906′ E) at a depth of ~ 10,910 m during research cruise TS-03 in January 2017 using the lander “TianYaHao” aboard the R/V TAN SUO YI HAO (see details in Additional file 2: Supplementary method).
DNA and RNA extraction
DNA was extracted from water samples of both fractions at all five depths as previously described . Briefly, biomass was washed from half of each filter using 15-ml extraction buffer (100 mM Tris-HCl [pH 8.0], 100 mM sodium EDTA [pH 8.0], 100 mM sodium phosphate [pH 8.0], 1.5 M NaCl, 1% CTAB) and then centrifuged at 5000×g for 20 min at room temperature. The concentrated biomass was ground within liquid nitrogen. Proteinase K and SDS were added in sequence and incubated at 37 °C and 65 °C, respectively, followed by phenol-chloroform extraction. DNA was precipitated with 0.6 volume of isopropanol, washed with 70% ethanol, air dried, and dissolved in TE buffer. The concentration and integrity of genomic DNA were analyzed by Qubit fluorometer and agarose gel electrophoresis, respectively. RNA was extracted using the Trizol reagent and cleaned using the Direct-zol RNA MiniPrep kit. DNA removal was performed using the TURBO DNA-free kit. Successful DNA removal was confirmed by 16S rRNA gene PCR. cDNA was generated using the Promega M-MLV reverse transcriptase.
Metagenomic sequencing, 16S rRNA gene high-throughput sequencing, and clone library construction
Nine DNA samples (each greater than 1.5 μg), except for the particle-associated fraction of 9600 m due to an insufficient quantity of DNA, were sent to BGI (BGI, Shenzhen, China) for metagenomic sequencing. Libraries were prepared without any amplification step for each sample. Metagenomic shotgun sequencing was performed on the Illumina HiSeq X-Ten platform, with 2 × 150-bp paired-end reads. Primers 515F and 806R (Additional file 1: Table S8) were used to amplify the V4 region of the 16S rRNA gene, and DNA extracted from three membranes were treated as triplicate. Sequencing was performed on the Illumina Miseq PE300 platform (MiSeq Reagent Kit v3) at Majorbio BiFo-Pharm Technology Co., Ltd., Shanghai, China. Sequences were processed with the pipeline of UPARSE . Operational taxonomic units (OTUs) were clustered at a 97% similarity level and taxonomically assigned against the SILVA SSU reference database (release 132). To examine the phylogenetic relationship between bacterial isolates and environmental sequences, a clone library was constructed using DNA from the FL10500 sample with primers B8F and B1510R (Additional file 1: Table S8). PCR products were inserted into pUCm-T vectors and cloned into Escherichia coli JM109. Positive clones were checked with PCR using primers M13F/R. Correct clones were sequenced at BGI (Qingdao, China) and taxonomically assigned using the EzTaxon-e server (https://www.ezbiocloud.net).
Quantitative PCR and reverse transcription qPCR
The abundance of bacterial 16S rRNA genes was quantified using quantitative PCR (qPCR) with primers Eub338F and Eub518R (Additional file 1: Table S8). Genomic DNA extracted from 10,500 m was used as a template to amplify 16S rRNA genes. PCR products were purified using a Gel Extraction Kit (Omega), ligated with pUCm-T vector, and then transformed into competent E. coli JM109 cells for sequencing. The sequence of target genes was submitted to the EzTaxon-e server (https://www.ezbiocloud.net) and identified as Ulvibacter sp. (Flavobacteria). The confirmed plasmid was extracted using a Plasmid Mini Kit (Omega) and digested with restriction enzyme XhoI. The linear plasmid was purified using a Gel Extraction Kit (Omega) and was serially diluted to construct standard curves (concentration ranged from 2.91 × 102 to 2.91 × 107 copies/L). The reverse transcription qPCR (RT-qPCR) standards for the alkB and almA gene (Additional file 1: Table S8) were prepared as detailed in the above description for the 16S rRNA gene. qPCR was performed on the StepOnePlus Real-Time PCR System based on SYBR Green I.
Assembly and binning
The reads that contain more than 10% of undefined bases and more than 40% of low-quality bases, and had a length of more than 15 bases matching to the adapters were removed. After filtering, clean data ranging in size from 13.67 to 16.54 Gb was obtained from each sample, and the average Q30 of each sample was above 93%, indicating that the trimmed dataset was of high quality. High-quality reads were assembled using SOAPdenovo2  with a kmer size of 49. Contigs shorter than 500 bp were discarded. Gene prediction was performed using MetaGeneMark  with default parameters (http://exon.gatech.edu/GeneMark/meta_gmhmmp.cgi).
Metagenomic binning was performed on four bottom water samples (FL10400, PA10400, FL10500, and PA10500). The initial paired-end reads from these four samples were co-assembled using IDBA-UD . Subsequently, these reads were mapped to the final contigs with BWA . Contigs were grouped into MAGs based on empirical probabilistic distances of genome abundance and tetranucleotide frequency by using program MetaBAT2 . Completeness and contamination of MAGs were assessed using CheckM , and MAGs with a completeness above 80% and contamination lower than 10% were considered for further analysis. PhyloPhlAn  was used to determine the phylogeny of MAGs, following by the average nucleotide identity (ANI) calculated to further confirm the taxonomy. The relative abundance of each MAG was calculated by mapping each of them to the raw reads. Gene calling of the MAGs was carried out using the RAST server . All the predicted genes were searched against the SEED subsystems, COG (ftp://ftp.ncbi.nlm.nih.gov/pub/COG/), and KEGG (http://www.kegg.jp/kegg/) databases.
Taxonomic assignment and functional characterization
Genes were clustered to remove redundant sequences using CD-Hit  at the 95% identity and 90% coverage. The gene annotation was done using BLASTp with an e-value cutoff of 10−5 against KEGG, COG, SEED, and NCBI-nr (released August 2016) databases for functional and taxonomic analysis, and only the best hits were retained. Taxonomy assignment was performed in MEGAN software  (version 4.6) based on the BLAST results of nr database using the lowest common ancestor (LCA) algorithm. To verify the taxonomic assignment results, metagenomic reads were mapped onto sequences of the SILVA SSU Ref NR99 database (release 132) using the BBMap package (http://sourceforge.net/projects/bbmap/) with 97% identity. The sequence coverage, which is the total length of mapped reads divided by the reference sequence length, was determined. The relative abundance of each taxon was calculated from the sequence coverage of each taxon (E) and the total sequence coverage (T) as E/T × 100%. Gene abundance was normalized by the number of the mapped reads of the target gene to the gene length. The taxa that were significantly enriched in the four NBW samples were described. Averaged abundances across the four NBW samples were provided. Standard deviations were calculated to account for within-group differences.
The marker genes (Additional file 3) for the major known types of aliphatic hydrocarbon-degradation enzymes, i.e., methane monooxygenase (MMO), gaseous short-chain alkane (C2–C4) monooxygenase, medium-chain alkane degradation enzymes (AlkB and cytochrome P450 monooxygenases, degrading C5–C26), long-chain alkane monooxygenase (represented by AlmA, degrading C14–C36), and other relevant genes potentially involved in hydrocarbon degradation were searched against the metagenomic data using BLASTp with an e-value cutoff of 10−5. An identity value of 40% was used for AlkB, cytochrome P450, and AlmA to increase the possibility of capturing the real functional proteins. All MMO sequences retrieved from the BLASTp results were further annotated in the NCBI-nr database, since they all showed low identity (< 42.9%) compared to marker genes. Only those sequences with top hits to methane monooxygenase in the NCBI-nr database were finally kept. Sequences of alkB and cytochrome P450 were retrieved from a functional gene pipeline (http://fungene.cme.msu.edu/). Short-chain (C2-C4) hydrocarbon monooxygenase and long-chain alkane monooxygenase were retrieved from NCBI according to Li et al.  and Wang and Shao , respectively. Reference sequences of methane monooxygenase were downloaded from the RefSeq database at NCBI. Genes involved in three intermediates (catechol, protocatechuate, and gentisate) of the aromatic degradation pathway were retrieved from the GeoChip 4.6 (http://www.ou.edu/ieg/tools) and RefSeq database at NCBI (Additional file 3). Statistics where applied were examined with the Wilcox ranking test.
Bacterial cultivation and n-alkane degradation
Water samples (10,400 m) were enriched using MMC medium (24 g NaCl, 7.0 g MgSO4·7H2O, 1 g NH4NO3, 0.7 g KCl, 2.0 g KH2PO4, 3.0 g Na2HPO4, and 10 ml of a trace element solution, pH 7.5) with hexadecane as a sole carbon source and cultivated at 10 °C, atmospheric pressure on a rotary shaker (170 rpm). A total of 38 bacterial strains were obtained. Taxonomic assignment was determined by using the EzBioCloud server (http://www.ezbiocloud.net/). Three Alcanivorax strains (A. jadensis ZYF844, A. venustensis ZYF848, and A. dieselolei ZYF854) were chosen for n-alkane degradation tests. Incubations were performed in 50 ml of MMC medium supplemented with 10 μl of alkane mixture at 4 °C, atmospheric pressure for 30 days. After confirmation that these species can utilize alkanes at low temperature, we further tested whether they could degrade n-alkanes at high pressure. The strains were incubated in 5 ml of MMC medium supplemented with 5 mg of n-eicosane at 4 °C, 60 MPa for 20 days. High pressure incubations were conducted in stainless steel reactors (380 ml, maximum pressure 60 MPa; Nantong Feiyu Oil Science and Technology Exploitation, China). Pressure was delivered by water using a manual pump. Once the incubation was finished, the remaining alkanes were extracted immediately using dichloromethane and analyzed using an Agilent 7890B Gas Chromatography (GC) with an FID detector. Oleibacter marinus DSM 24913 (= NBRC 105760) bought from the Marine Culture Collection of China was also tested for n-alkane degradation at 4 °C, 60 MPa, but demonstrated poor growth. Thus, this strain was incubated at 16 °C, atmospheric pressure with n-eicosane as a sole carbon source (see details in Additional file 2: Supplementary method).
Analysis of n-alkanes
The remaining alkanes in 50- or 5-ml bacterial cultures were extracted using 10 or 2 ml of dichloromethane and analyzed using an Agilent 7890B Gas Chromatograph (GC) with an FID detector. Two further rounds of extraction were performed. A deuterated alkane (C24D50) was used as an internal recovery standard (see details in Additional file 2: Supplementary method). The n-alkane analysis procedure for sinking particle and surface sediment samples was the same as that of the bacterial cultures except for extraction. Sinking particles and surface sediment samples were first freeze-dried and ground using a mortar and pestle. In order to obtain sufficient n-alkane yields from sinking particles at each depth, we combined the samples collected at nine time points (from January to September). Alkanes were extracted from 700 mg of sinking particles (for each depth) and 5.0 g surface sediments using 30 ml dichloromethane/methanol (9:1, v/v) for 10 min. This procedure was repeated four times and the extracts were combined. Dichloromethane and hexane without sample were used as a blank. All the extraction process was conducted in a clean room, and all the glassware were baked at 550 °C for 12 h to remove any organic matter. The alkanes were not analyzed in the original water samples since we were not expecting to find such a large proportion of hydrocarbon degraders below 10,000 m.
Analysis of carbon and hydrogen isotopic compositions of n-alkanes
Once isolated, biomarker n-alkanes were analyzed first by gas chromatography (Thermo Scientific Trace 1310) with a flame ionization detector (GC-FID) fitted to a 30-m TG-5SilMS (0.25 mm × 0.25 μm) fused silica column. Sample aliquots of 1 μl were injected by programmed temperature vaporizing (PTV) injector at 320 °C through a split/splitless silcosteel inlet liner. The chromatograph oven was ramped from 60 °C (1 min) to 320 °C (20 min) at 6 °C min−1. Helium was used as a carrier at 2 ml min−1, and the detector was set to a constant temperature of 320 °C. Stable isotope measurements of biomarker n-alkanes were measured by gas chromatography combustion (δ13C) or pyrolysis (δ2H) isotope ratio monitoring mass spectrometry (GC-C-irMS or GC-pyr-irMS, respectively) with a Thermo Scientific Trace 1310 GC with split/splitless injector and TriPlus RSH liquid autosampler coupled to a Delta V Plus MS (see details in Additional file 2: Supplementary method).
We would like to thank all of the scientists and crews on the R/V Dong Fang Hong 2 and R/V TAN SUO YI HAO for their assistance with sampling during the cruise. We also thank Prof. Curtis Suttle from the University of British Columbia, Canada and Dr. Boyd McKew from the University of Essex for critical reading and discussion. We thank Dr. James V. Gardner from the Center for Coastal & Ocean Mapping, University of New Hampshire, for providing Additional file 1: Figure S1.
This work was funded by the National Natural Science Foundation of China (91751202, 41730530, 41476112 and 41506154), Fundamental Research Funds for the Central Universities (201762009), Scientific and Technological Innovation Project of Qingdao National Laboratory for Marine Science and Technology (2016ASKJ14 and 2016ASKJ15), National Program on Global Change and Air-Sea Interaction (GASI-02-PAC-ST-MSwin) and the Strategic Priority Research Program from the Chinese Academy of Sciences (XDB06010203). JDT was funded by Natural Environmental Research (NERC) grants NERC standard grants NE/S001352, NE/P012671 and NE/N002385.
Availability of data and materials
Sequence data for metagenome and 16S rRNA reads have been deposited in the NCBI Sequence Read Archive under accession numbers SRP119520 and SRP119680, respectively. The Oceanospirillales MAGs have been deposited at GenBank under the accession numbers QAGL00000000 (MAG 01), PDUK00000000 (MAG 78), PDUL00000000 (MAG 5.60), PDUM00000000 (MAG 12) and PZPK00000000 (MAG 156).
X-HZ designed the experiments and analyzed the data. JT designed the cruise and the large-volume water sampler. JL and YZ collected the water samples and analyzed the environmental sequence data. YZ extracted the community DNA, isolated the strains from seawater, and detected the n-alkane degradation capacity. HL, ML, YL, MY, and HZ performed the metagenomic binning and pathway reconstruction. XW, MY, MZ, and YZ measured the n-alkanes quantity. YX analyzed the environmental factors. WZ collected the sediment samples. NP and CM performed the compound-specific stable isotope measurements and interpreted δ13C and δ2H data. SZ and DS extracted the RNA from seawater and did the RT-qPCR. DJL-S, NP, and JDT provided critical ideas for interpreting the genomic and n-alkane distribution data. X-HZ, JL, YZ, DJL-S, JDT, and NP wrote the manuscript. All authors edited and approved the final manuscript.
Ethics and consent for publication
To carry out sampling, the authors received a ‘Research and Training Permit for the Exclusive Economic Zone of the Federated States of Micronesia’ issued by the National Oceanic Resource Management Authority - Federated States of Micronesia.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Nunoura T, Takaki Y, Hirai M, Shimamura S, Makabe A, Koide O, et al. Hadal biosphere: insight into the microbial ecosystem in the deepest ocean on Earth. Proc Natl Acad Sci U S A. 2015;112:E1230–6.Google Scholar
- Tarn J, Peoples LM, Hardy K, Cameron J, Bartlett DH. Identification of free-living and particle-associated microbial communities present in hadal regions of the Mariana Trench. Front Microbiol. 2015;7:665.Google Scholar
- Tian J, Fan L, Liu H, Liu J, Li Y, Qin Q, et al. A nearly uniform distributional pattern of heterotrophic bacteria in the Mariana Trench interior. Deep Sea Res Part 1 Oceanogr Res Pap. 2018;142:116–26.Google Scholar
- Eloe EA, Shulse CN, Fadrosh DW, Williamson SJ, Allen EE, Bartlett DH. Compositional differences in particle-associated and free-living microbial assemblages from an extreme deep-ocean environment. Environ Microbiol Rep. 2011;3:449–58.View ArticleGoogle Scholar
- Nunoura T, Hirai M, Yoshida-Takashima Y, Nishizawa M, Kawagucci S, Yokokawa T, et al. Distribution and niche separation of planktonic microbial communities in the water columns from the surface to the hadal waters of the Japan Trench under the Eutrophic Ocean. Front Microbiol. 2016;7:1261.View ArticleGoogle Scholar
- Peoples LM, Donaldson S, Osuntokun O, Xia Q, Nelson A, Blanton J, et al. Vertically distinct microbial communities in the Mariana and Kermadec trenches. PLoS One. 2018;13:e0195102.View ArticleGoogle Scholar
- Lacey NC, Rowden AA, Clark MR, Kilgallen NM, Linley T, Mayor DJ, et al. Community structure and diversity of scavenging amphipods from bathyal to hadal depths in three South Pacific Trenches. Deep Sea Res Part 1 Oceanogr Res Pap. 2016;111:121–137.Google Scholar
- León-Zayas R, Novotny M, Podell S, Shepard CM, Berkenpas E, Nikolenko S, et al. Single cells within the Puerto Rico trench suggest hadal adaptation of microbial lineages. Appl Environ Microbiol. 2015;81:8265–76.View ArticleGoogle Scholar
- Sun DL, Jiang X, Wu QL, Zhou NY. Intragenomic heterogeneity in 16S rRNA genes causes overestimation of prokaryotic diversity. Appl Environ Microb. 2013;79:5962–9.View ArticleGoogle Scholar
- Uyaguari-Diaz MI, Chan M, Chaban BL, Croxen MA, Finke JF, Hill JE, et al. A comprehensive method for amplicon-based and metagenomic characterization of viruses, bacteria, and eukaryotes in freshwater samples. Microbiome. 2016;4:20.View ArticleGoogle Scholar
- Mattick JS. Type IV pili and twitching motility. Annu Rev Microbiol. 2002;56:289–314.View ArticleGoogle Scholar
- Teramoto M, Ohuchi M, Hatmanti A, Darmayati Y, Widyastuti Y, Harayama S, et al. Oleibacter marinus gen. nov., sp. nov., a bacterium that degrades petroleum aliphatic hydrocarbons in a tropical marine environment. Int J Syst Evol Microbiol. 2011;61:375–80.View ArticleGoogle Scholar
- Liu J, Bacosa HP, Liu Z. Potential environmental factors affecting oil-degrading bacterial populations in deep and surface waters of the northern Gulf of Mexico. Front Microbiol. 2017;7:2131.View ArticleGoogle Scholar
- Wang WP, Shao ZZ. The long-chain alkane metabolism network of Alcanivorax dieselolei. Nat Commun. 2014;5:5755.View ArticleGoogle Scholar
- Rivers AR, Sharma S, Tringe SG, Martin J, Joye SB, Moran MA. Transcriptional response of bathypelagic marine bacterioplankton to the Deepwater Horizon oil spill. ISME J. 2013;7:2315–29.View ArticleGoogle Scholar
- Mason OU, Hazen TC, Borglin S, Chain PS, Dubinsky EA, Fortney JL, et al. Metagenome, metatranscriptome and single-cell sequencing reveal microbial response to Deepwater Horizon oil spill. ISME J. 2012;6:1715–27.View ArticleGoogle Scholar
- Hazen TC, Dubinsky EA, DeSantis TZ, Andersen GL, Piceno YM, Singh N, et al. Deep-sea oil plume enriches indigenous oil-degrading bacteria. Science. 2012;330:204–8.View ArticleGoogle Scholar
- Redmond MC, Valentine DL. Natural gas and temperature structured a microbial community response to the Deepwater Horizon oil spill. Proc Natl Acad Sci U S A. 2012;109:20292–7.View ArticleGoogle Scholar
- Valentine DL, Kessler JD, Redmond MC, Mendes SD, Heintz MB, Farwell C, et al. Propane respiration jump-starts microbial response to a deep oil spill. Science. 2010;330:208–11.View ArticleGoogle Scholar
- Hu P, Dubinsky EA, Probst AJ, Wang J, Sieber CMK, Tom LM, et al. Simulation of Deepwater Horizon oil plume reveals substrate specialization within a complex community of hydrocarbon degraders. Proc Natl Acad Sci U S A. 2017;114:7432–7.View ArticleGoogle Scholar
- Delmont TO, Eren AM. Simulations predict microbial responses in the environment? This environment disagrees retrospectively. Proc Natl Acad Sci U S A. 2017;114:E8947–9.View ArticleGoogle Scholar
- Gutierrez T, Singleton DR, Berry D, Yang T, Aitken MD, Teske A. Hydrocarbon-degrading bacteria enriched by the Deepwater Horizon oil spill identified by cultivation and DNA-SIP. ISME J. 2013;7:2091–104.View ArticleGoogle Scholar
- Liu C, Wang W, Wu Y, Zhou Z, Lai Q, Shao Z. Multiple alkane hydroxylase systems in a marine alkane degrader, Alcanivorax dieselolei B-5. Environ Microbiol. 2011;13:1168–78.View ArticleGoogle Scholar
- Ivanova EP, López-Pérez M, Zabalos M, Nguyen SH, Webb HK, Ryan J, et al. Ecophysiological diversity of a novel member of the genus Alteromonas, and description of Alteromonas mediterranea sp. nov. Antonie Van Leeuwenhoek. 2015;107:119–32.View ArticleGoogle Scholar
- Ivars-Martinez E, D'Auria G, Rodríguez-Valera F, Sânchez-Porro C, Ventosa A, Joint I, et al. Biogeography of the ubiquitous marine bacterium Alteromonas macleodii determined by multilocus sequence analysis. Mol Ecol. 2008;17:4092–106.View ArticleGoogle Scholar
- Simoneit BRT, Grimait JO, Fischer K, Dymond J. Upward and downward flux of particulate organic material in abyssal waters of the Pacific Ocean. Naturwissenschaften. 1986;73:322–5.View ArticleGoogle Scholar
- Luo M, Gieskes J, Chen LY, Shi XF, Chen DF. Provenances, distribution, and accumulation of organic matter in the southern Mariana Trench rim and slope: implication for carbon cycle and burial in hadal trenches. Mar Geol. 2017;386:98–106.View ArticleGoogle Scholar
- Heeszel DS, Wiens DA, Shore PJ, Shiobara H, Sugioka H. Earthquake evidence for along-arc extension in the Mariana Islands. Geochem Geophy Geosy. 2008;9:12.View ArticleGoogle Scholar
- Grimalt J, Albaigés J. Sources and occurrence of C12-C22 n-alkane distributions with even carbon-number preference in sedimentary environments. Geochim Cosmochim Acta. 1987;51:1379–84.View ArticleGoogle Scholar
- Guan HX, Chen LY, Luo M, Liu LH, Mao SY, Ge HM, et al. Composition and origin of lipid biomarkers in the surface sediments from the southern Challenger Deep, Mariana Trench. Geosci Front. 2019;10:351–60.Google Scholar
- Pedentchouk N, Turich C. Carbon and hydrogen isotopic compositions of n-alkanes as a tool in petroleum exploration. Geol Soc London Spec Publ. 2018;468:105–25.View ArticleGoogle Scholar
- Li C, Sessions AL, Kinnaman FS, Valentine DL. Hydrogen-isotopic variability in lipids from Santa Barbara Basin sediments. Geochim Cosmochim Acta. 2009;73:4803–23.View ArticleGoogle Scholar
- Zhang X, Gillespie AL, Sessions AL. Large D/H variations in bacterial lipids reflect central metabolic pathways. Proc Natl Acad Sci U S A. 2009;106:12580–6.View ArticleGoogle Scholar
- Lea-Smith DJ, Biller SJ, Davey MP, Cotton CA, Perez Sepulveda BM, Turchyn AV, et al. Contribution of cyanobacterial alkane production to the ocean hydrocarbon cycle. Proc Natl Acad Sci U S A. 2015;112:13591–6.View ArticleGoogle Scholar
- Schirmer A, Rude MA, Li X, Popova E, Del Cardayre SB. Microbial biosynthesis of alkanes. Science. 2010;329:559–62.View ArticleGoogle Scholar
- Sorigué D, Légeret B, Cuiné S, Blangy S, Moulin S, Billon E, et al. An algal photoenzyme converts fatty acids to hydrocarbons. Science. 2017;357:903–7.View ArticleGoogle Scholar
- Gardner JV, Armstrong AA, Calder BR, Beaudoin J. So, how deep is the Mariana Trench? Mar Geod. 2014;37:1–13.View ArticleGoogle Scholar
- Zhou J, Bruns MA, Tiedje JM. DNA recovery from soils of diverse composition. Appl Environ Microbiol. 1996;62:316–22.PubMedPubMed CentralGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.View ArticleGoogle Scholar
- Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2, an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012;1:18.View ArticleGoogle Scholar
- Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucleic Acids Res. 2010;38:e132.View ArticleGoogle Scholar
- Peng Y, Leung HCM, Yiu SM, Chin FYL. IDBA-UD, a de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics. 2012;28:1420–8.View ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticleGoogle Scholar
- Kang DWD, Froula J, Egan R, Wang Z. MetaBAT, an efficient tool for accurately reconstructing single genomes from complex microbial communities. Peerj. 2015;3:e1165.View ArticleGoogle Scholar
- Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM, assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.View ArticleGoogle Scholar
- Segata N, Bornigen D, Morgan XC, Huttenhower C. PhyloPhlAn is a new method for improved phylogenetic and taxonomic placement of microbes. Nat Commun. 2013;4:2304.View ArticleGoogle Scholar
- Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42:D206–14.View ArticleGoogle Scholar
- Li W, Godzik A. Cd-hit, a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.View ArticleGoogle Scholar
- Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17:377–86.View ArticleGoogle Scholar
- Li M, Jain S, Baker BJ, Taylor C, Dick GJ. Novel hydrocarbon monooxygenase genes in the metatranscriptome of a natural deep-sea hydrocarbon plume. Environ Microbiol. 2014;16:60–71.View ArticleGoogle Scholar