Genomic and Transcriptomic Evidence for Carbohydrate Consumption among Microorganisms in a Cold Seep Brine Pool

Abstract

The detailed lifestyle of microorganisms in deep-sea brine environments remains largely unexplored. Using a carefully calibrated genome binning approach, we reconstructed partial to nearly-complete genomes of 51 microorganisms in biofilms from the Thuwal cold seep brine pool of the Red Sea. The recovered metagenome-assembled genomes (MAGs) belong to six different phyla: Actinobacteria, Proteobacteria, Candidatus Cloacimonetes, Candidatus Marinimicrobia, Bathyarchaeota, and Thaumarchaeota. By comparison with close relatives of these microorganisms, we identified a number of unique genes associated with organic carbon metabolism and energy generation. These genes included various glycoside hydrolases, nitrate and sulfate reductases, putative bacterial microcompartment biosynthetic clusters (BMC), and F₄₂₀H₂ dehydrogenases. Phylogenetic analysis suggested that the acquisition of these genes probably occurred through horizontal gene transfer (HGT). Metatranscriptomics illustrated that glycoside hydrolases are among the most highly expressed genes. Our results suggest that the microbial inhabitants are well adapted to this brine environment, and anaerobic carbohydrate consumption mediated by glycoside hydrolases and electron transport systems (ETSs) is a dominant process performed by microorganisms from various phyla within this ecosystem.

Introduction

Microorganisms play important roles in the biogeochemical cycles of deep-sea environments, such as deep-sea hydrothermal vents and brine pools. Studies of microbial lifestyles in brine pools are less common than those of hydrothermal vents. The brine pools studied so far include the Shaban Deep (Ferrer et al., 2012), Discovery Deep (Wang et al., 2013), Atlantis II Deep (Ngugi et al., 2015), and Kebrit Deep (Guan et al., 2015) of the Red Sea, which contains more than 25 deep hypersaline anoxic pools (Antunes et al., 2011). Previous investigations of other deep-sea brine pools included those in the Mediterranean Sea (van der Wielen et al., 2005) and the Gulf of Mexico (Formolo and Lyons, 2013). Brine pools probably share a common formation process of tectonically induced brine formation, resulting in frequently observed high methane and sulfite fluxes (Faber et al., 1998). As a consequence, the prevalent processes in the microbial inhabitants are anaerobic oxidation of methane and sulfite and sulfate reduction (Antunes et al., 2011). Surveys of brine pools identified a few widespread groups including Proteobacteria, Actinobacteria, Cyanobacteria, and Deferribacteres, and showed a higher diversity of bacteria over the archaea (Ferrer et al., 2012; Guan et al., 2015; Ngugi et al., 2015). Despite the common origin and shared features of the investigated brine pools, regional variation may lead to distinct combinations of physicochemical parameters and sustains microorganisms with unique adaptive strategies.

As important components of marine ecosystems, biofilms play an important role in ocean carbon transformation (Bhaskar and Bhosle, 2005; D'Arcy et al., 2013). However, biofilms are often very complex in nature. Heterogeneity of natural substrates induces heterogeneous biofilm development that is difficult to measure and observe (Claret, 1998). Devices for sampling deep-sea biofilms are generally lacking, making their study challenging. To avoid these problems, one approach to studying the ecological roles of biofilms is to use artificial substrates to develop biofilms under steady-state conditions. This method has been used in previous studies of biofilm development in hydrothermal vent areas (Guezennec et al., 1998), estuary systems (Jones et al., 2007), and intertidal zones (Salta et al., 2013; Zhang et al., 2013).

Our previous investigation on the metagenomics of biofilms developed on artificial surfaces (Zhang et al., 2014, 2015) has illustrated the developmental processes of biofilms in a brine pool associated with the Thuwal cold seeps (Batang et al., 2012) on the floor of the Red Sea. This brine pool is located at a depth of approximately 850 m and is characterized by a high concentration of potential electron acceptors such as nitrate (18.5 mg L⁻¹) and sulfate (2.7 g L⁻¹). The particulate organic carbon (POC) concentration was 3.5 mg L⁻¹ which was 10 times greater than that of the adjacent normal deep-sea water (NDW) (Zhang et al., 2014, 2015). In contrast, methane and hydrogen sulfide were depleted with respect to NDW. Metagenomics were performed for biofilms developed on three different types of substrates: aluminum (Al), titanium (Ti) and polyvinyl chloride (PVC). Comparisons between brine and NDW biofilms indicated that polysaccharide metabolism related genes were highly enriched in the brine biofilms (Zhang et al., 2015). We also reconstructed two dominant bacterial genomes from the brine biofilm metagenomics datasets, a novel deltaproteobacterium and a novel epsilonproteobacterium. The former metagenome-assembled genome (MAG) possessed enhanced polysaccharide fermentation pathways, whereas the later possessed complete gene clusters for nitrate reduction and nitrogen fixation. We proposed that polysaccharide fermentation and proteolysis interacted with nitrogen cycling to form a complex chain for energy generation in these biofilms (Zhang et al., 2015).

Although, the biofilm development process in the Thuwal cold seep brine pool and the genomics of the two dominant microorganisms has been analyzed, the lifestyle of the other microorganisms in the brine biofilms remains unclear. In the present study, using the genome binning technique, we analyzed 51 partial to near-complete microbial genomes reconstructed from the Al and PVC biofilm metagenomes and the metatranscriptomes derived from these biofilms, toward a systematic understanding of microbial lifestyles in this ecosystem.

Materials and methods

Biofilm sampling, DNA extraction and metagenomic analysis

Biofilm sampling was conducted in May 2013 in the Thuwal seeps II (22°16N-38°53E) by the remotely operated underwater vehicle (ROV) Max Rover, DSSI, USA during the Red Sea exploration cruise. Biofilms were developed on six different types of materials in the brine pool for 3 and 6 days. The devices for biofilm development has been reported in Zhang et al. (2014) and the schematic diagrams are also shown in Figure S1 in the present study. The materials were deployed into the brine pool by the ROV. Environmental parameters, such as temperature, salinity, and dissolved oxygen (DO) were measured using conductivity-temperature-depth (CTD) and SBE43 DO sensors attached to the ROV.

The biofilms used in the present study were developed on aluminum and PVC in the brine pool for 6 days and were referred to as “biofilm_Al and biofilm_PVC.” Biofilms on these two materials were selected because biofilm_Al had the highest biomass and PVC represented a biofilm developed on plastic materials. Cells were recovered from the substrates using sterile cell scrapers and suspended in Tris-HCl buffer before extracting nucleic acid using the AllPrep DNA/RNA Mini Kit (Qiagen, Hilden, Germany). The extracted DNA was sequenced using an Illumina HiSeq 2000 platform. Metagenomic datasets of biofilm_Al and biofilm_PVC have been documented in Zhang et al. (2015), and the information of these two metagenomes are also listed in Table S2 in the present study. Assembly of the combined metagenomes of biofilm_Al and biofilm_PVC was performed using SPAdes Genome Assembler 3.6.1 (Bankevich et al., 2012). The specified K values for assembly included 21, 31, 41, 51, 61, 71, and 81, and the “—careful” and “—pe” options were used. The metagenome-based taxonomic structures were obtained by ribosomal database project (RDP) classification (Wang et al., 2007) of the extracted 16S rRNA sequences, and the coverage information.

Genome binning and validation

Draft genome binning was conducted according to previously described methods (Albertsen et al., 2013; Tian et al., 2014). The genomes were separated mainly on the genome coverage, GC content, and tetranucleotide frequency. The Illumina reads were mapped to the contigs using Bowtie2 (version 2.0.0; Langmead and Salzberg, 2012). The genome coverage was calculated using SAMtools on a local server. The contigs were searched against a set of single-copy protein-encoding genes (Albertsen et al., 2013) using default cut-off values. The single-copy protein-encoding genes identified were searched against the National Center for Biotechnology Information non-redundant (NCBInr) database with BLASTP (e-value < 1e-07) with “xml” output format. Taxonomic information for the contigs was obtained by importing the BLAST results into MEGAN 5.0 (Huson et al., 2011). The metagenome pairs biofilm_Al and biofilm_PVC were selected to obtain the most effective separation and extraction of contigs. It must be noted that the diversity of the biofilm communities was not high and that many contigs belonging to different genomes formed clear, separate clusters based on coverage, facilitating genome binning and elimination of contamination (Figure S2).

Potential contig contamination was examined using CheckM (Parks et al., 2015), and PhylopythiaS, which is a sequence composition-based classifier that utilizes the hierarchical relationships between clades (Patil et al., 2012). Moreover, the Thaumarchaeota MAGs in the present study were compared with the eight genomes obtained using single-cell genomics by Ngugi et al. (2015), based on BLASTP (e-value < 1e-07) and BLASTN (e-value < 1e-07) with predicted open reading frames (ORFs) as input. The results were visualized using the Artemis Comparison Tool (ACT; Carver et al., 2005).

Genome annotation

The MAGs were annotated by searching the Carbohydrate-Active EnZymes (CAZys) (Cantarel et al., 2009), Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa and Goto, 2000), and Clusters of Orthologous Groups (COGs) (Galperin et al., 2014) databases. The ORFs in the MAGs were predicted using Prodigal (Hyatt et al., 2010) using a local server. Genes were annotated by BLASTP searches against the above-listed databases using an e-value cut-off of 1e-07 on a local server. The online tool KAAS (KEGG Automatic Annotation Server) was also used to validate the KEGG annotation by local BLAST. Metabolic pathways were analyzed using online tools in KEGGMAPPER (http://www.genome.jp/kegg/mapper.html). Key enzymes revealed by KEGGMAPPER were rechecked against the NCBInr database based on the BLAST results.

Phylogenetic tree construction

In the phylogenetic analysis of the functional genes (e.g., polysaccharide deacetylase), the protein sequence alignment was determined by the molecular evolutionary genetics analysis using Muscle (MEGA 6.06; Tamura et al., 2013) with the following parameters: gap open -50, cluster method UPGMB and min diagonal length 24. Gblock software (Castresana, 2000) was applied to eliminate less informative sites in the alignments. The construction of maximum-likelihood trees was conducted using MEGA 6.06 with the Tamura-Nei model, the Nearest-Neighbor-Interchange (NNI) method and 500 bootstrap replicates.

RNA extraction, cDNA amplification and metatranscriptomic sequencing

Biofilms on different substrates were harvested (10 slides per sample) onboard immediately after recovery using sterile cotton tips and stored in RNAlater buffer. RNA extraction was performed using the AllPrep DNA/RNA Mini Kit (Qiagen, Hilden, Germany). Bacterial cells in RNAlater buffer were pelleted by centrifugation at 4000 g for 10 min and then lysed with lysozyme, proteinase K and RLT buffer provided in the kit. Total nucleic acid was extracted, and DNA was removed using the DNA column provided with the AllPrep DNA/RNA Mini Kit. The integrity and quantity of the RNA were checked by agarose gel electrophoresis and a Nanodrop device (ND-1000 spectrophotometer, DiaMed China Limited, Hong Kong), respectively. Prior to first-strand cDNA synthesis, genomic DNA contamination was further removed using the DNA-free™ DNA Removal Kit (Ambion, Austin, TX). The extracted RNA was then used as a template for whole-cDNA amplification using the Ovation Single-Cell RNA-Seq System, which enriches for non-rRNA during cDNA synthesis using proprietary whole transcriptome primers targeting non-rRNA sequences in the transcriptome. The synthesized cDNA was sequenced using an Illumina HiSeq 2000 platform.

Metatranscriptomic analysis

The information for metatranscriptomic sequences are summarized in Table S2. The sequenced cDNA was subjected to quality control using the next-generation sequencing (NGS) QC toolkit (Patel and Jain, 2012). Assembly of the metatranscriptomic reads was performed using SPAdes Genome Assembler 3.6.1 (Bankevich et al., 2012) on a local server. The specified K values 21, 31, 41, 51, 61, 71, and 81 were used under the “—careful” and “—pe” options. The ORFs were predicted following the steps mentioned above and annotated by BLASTP against the COG -database. To depict the gene expression profiling for individual microorganisms, the reads were aligned to contigs of the MAGs using Bowtie2 (version 2.0.0; Langmead and Salzberg, 2012). Gene expression abundance is the number of transcripts mapped to the CAZys of each MAG that normalized to the length of the genes.

Data accessibility

The MAGs of Cloacimonetes were deposited in the Whole Genome Shotgun (WGS) database under accession numbers of SAMN04423165, SAMN04423166, SAMN04423167, and SAMN04423168. The MAGs of Marinimicrobia were deposited in the WGS database under accession numbers of SAMN04423157, SAMN04423158, SAMN04423159, SAMN04423160, SAMN04423161, SAMN04423162, and SAMN04423164. The MAG of Bathyarchaeota was deposited in the WGS database under accession number of SAMN04423169.

Results

Genome information

Genome binning using two metagenomes, biofilm_Al and biofilm_PVC, generated 51 MAGs of the biofilm-inhabited microorganisms (TCS1-51; “TCS” served as the abbreviation for “Thuwal cold seep”), including members of Actinobacteria (two MAGs belonging to Candidatus Microthrix), Gammaproteobacteria, Deltaproteobacteria, Epsilonprtoteobacteria, Candidatus Marinimicrobia, Candidatus Cloacimonetes, Bathyarchaeota (previously referred to as Miscellaneous Crenarchaeotic Group, MCG) and Thaumarchaeota (Table 1). Based on the metagenome-derived 16S rRNA genes, these phyla were the dominant microbial groups in the biofilms (Figure S3). According to the results of CheckM, these MAGs had a maximum completeness value of 100% and a maximum potential contamination of 2.4%. Based on the results of PhylopythiaS, contigs of most of the MAGs were very consistent, as suggested by the percentage of contigs assigned to a defined taxonomy. For example, 97.2% of the Pseudomonas sp. nov (TCS13) contigs could be assigned to the genus Pseudomonas. Moreover, contigs of the newly defined bacterial phyla, were compared with their relatives with complete or single-cell derived genomes to assess the potential contamination. For example, 95.8% of the Thaumarchaeota sp. nov. 1 (TCS50) and 91.3% of the Thaumarchaeota sp. nov. 2 (TCS51) could be assigned to the eight single-cell Thaumarchaeota genomes, which were from another brine pool in the Red Sea (Ngugi et al., 2015). The assigned ORFs were distributed in all of the contigs of the TCS MAGs as exemplified by Thaumarchaeota sp. nov. 1 (Figure S4). These results indicated that the MAGs reconstructed in the present study were of high quality.

Central metabolic pathway

A large proportion of the genes identified in many of the MAGs were found to be involved in fermentative carbohydrate metabolism. The distribution of CAZys (Cantarel et al., 2009) responsible for the degradation of complex carbohydrates is listed in Table S1, such as β-glucosidase (GH3), α-amylase (GH13) and endo-β-1,4-galactanase (GH53). In total 17 unique CAZy categories corresponding to 1692 genes were identified in the 51 MAGs (on average 33 genes/genome). Genes encoding proteins playing a role in electron transport and energy production, including those encoding F₄₂₀H₂ dehydrogenase (fpo), cytochrome C (cytoC), the rhodobacter nitrogen fixation (rnf) complex, nitrate reductase and sulfate reductase, were identified in many of the MAGs. Annotation by COGs (Galperin et al., 2014) showed a high diversity in the gene arrangement of the F₄₂₀ biosynthesis gene clusters (Figure S5), which were present in eight of the MAGs belonging to different phyla. Moreover, putative bacterial microcompartment (BMC) clusters, types of protein-bound bacterial organelles that function in carbon catabolism and storage (Nobu et al., 2016), were present in five MAGs including those of Actinobacteria, Marinimicrobia and Desulfobacter (Figure S6). As reported by Nobu et al. (2016), genes related to aldehyde and sugar metabolism are present in the BMC clusters, including rnf-C like NADH dehydrogenase (COG4656), sugar isomerase (COG0363) and aldehyde dehydrogenase (COG1012). However, other genes were also identified in the present study, such as NTP pyrophosphohydrolases and oxidative damage repair enzymes (COG0494), making the BMC clusters different from those in previous study.

The prevalence and versatility of carbohydrate metabolism and energy generation was further supported by KEGG-based pathway reconstruction. The central metabolic strategies are summarized in Figure 1. Many of the MAGs encoded proteins associated with monosaccharide and polysaccharide transport, nitrate and sulfate reduction, electron transport systems (ETSs), glycolysis, and reoxidation of energy-conserving elements via acetyl-CoA reduction to ethanol or acetate. For example, monosaccharide and polysaccharide transport was indicated by the presence of a number of ABC transporters in all the bacterial genomes; dissimilatory nitrate reduction was indicated by both respiratory nitrate reductase (nar) and periplasmic nitrate reductases (nap); dissimilatory sulfate reduction was indicated by sulfur oxidation (SOX) gene cluster and adenylylsulfate reductase; and acetate generation was indicated by acetate kinase. Notably, the Marinimicrobia bacterium sp. nov. 2 (TCS39) and Bathyarchaeota sp. nov (TCS49) possessed both dissimilatory nitrogen and sulfur reduction pathways.

Figure 1

Unique genomic features

When the TCS MAGs were compared with reference genomes within the same phylum, unique features concerning carbohydrate metabolism and energy generation were observed. For example, genes encoding NADH dehydrogenase, such as nuoB, nuoL, ndufv2, NdhJ, and NoxF were only identified in four Cloacimonetes MAGs (TCS45-47) when compared with complete or single-cell Cloacimonetes genomes (Figure S7); similar results were observed for chitin deacetylase (EC3.5.1.41) and the α-D-glucose transformation pathway, which are important steps in glycolysis. Moreover, the complete pyruvate to acetate pathway was identified in Marinimicrobia bacterium sp. nov. 2 (TCS39) and Marinimicrobia bacterium sp. nov. 3 (TCS40), but the enzymes included in this pathway could not be found in the four reference genomes, the single-cell Marinimicrobia genomes from previous studies (Rinke et al., 2013; Nobu et al., 2015) (Figure S8). In addition, the Bathyarchaeota MAG (TCS49) harbored specific genes for chitin, cellulose, and β-D-glucose metabolism, such as a glycoside hydrolase family 3 gene (probably chitin deacetylase; EC3.5.1.41) and β-glucosidase (EC3.2.1.4) (Figure S9); the presence of aprAB and cysC genes, which are involved in dissimilatory sulfate reduction, differentiated TCS49 from the reference genomes.

Given that a number of unique genes were identified from these MAGs as compared with their relatives, we examined the phylogeny of these genes. The polysaccharide deacetylase (EC3.2.1.21) in Cloacimonetes bacterium sp. nov. 1 (TCS45) was phylogenetically close to the polysaccharide deacetylase from Lactobacillus ceti (Figure 2). Although, the bootstrap values in the tree were not high, BLAST analysis also revealed close relationship between polysaccharide deacetylase in Cloacimonetes bacterium sp. nov. 1 and those from firmicutes. The nitrate reductase α and β units in Marinimicrobia bacterium sp. nov. 3 (TCS40) were phylogenetically close to the nitrate reductases from Geothrix fermentans (Figures S10, S11). In addition, the glycoside hydrolase family 3 gene in Bathyarchaeota sp. nov (TCS49) was affiliated with Hungatella hathewayi and Clostridium hathewayi (Figure S12). The adjacent genes in the same contigs, which were affiliated with Cloacimonetes, Marinimicrobia or Bathyarchaeota, were displayed to confirm that the presence of these novel genes was probably not attributable to contig contamination.

Figure 2

Genome-wide expression profiling

To explore the major active functions, we obtained ~600 Mb of metatranscriptomic datasets for the biofilms developed on Al and PVC. The metatranscriptomic sequences were searched against the COG database, revealing an overall picture of the active functions. Transcripts of the genes involved in carbohydrate transport and metabolism were among the most abundant in the metatranscriptomes, particularly the glycosidases (COG0366), β-glucosidase-related glycosidases (COG1472) and β -galactosidase/beta-glucuronidase (COG3250) (Figure 3). In contrast, few peptidases that are essential for protein metabolism were identified, and rank analysis suggested relative lower expression levels for these genes (Figure S13).

Figure 3

The reads from these two metatranscriptomes were mapped to the 51 MAGs. Abundant active bacteria could be identified for nearly all of the phyla, whereas rare active bacteria were also included in each phylum. Key functional genes involved in nitrate and sulfate respiration were expressed, such as the napAB, narGH, nirK, norC, nosZ, and sulfur oxidation (SOX) gene cluster. The gene expression information was integrated in to pathways shown in Figure 1. Moreover, when the metatranscriptomic sequences were mapped to CAZys in each MAG, a heatmap showing expression pattern of CAZys was generated (Figure 4). The results confirmed the expression of diverse carbohydrate metabolism related genes, and indicated relative higher expression of the glycoside hydrolases (GHs), glycosyltransferase (GTs) than enzymes in other families.

Figure 4

Discussion

We reconstructed 51 partial to nearly complete MAGs with the aim of investigating the microbial lifestyles in brine pool biofilm communities. We validated the genome binning results using different approaches before the subsequent analyses. Genomic and transcriptomic analyses provided evidence that organic carbon metabolism utilizing nitrate, nitrite, and sulfate as potential electron acceptors is prevalent among these microorganisms. In particular, we identified several unique metabolic pathways for Cloacimonetes, Marinimicrobia, and Bathyarchaeota, which differ from previously reported lifestyles for these microbial groups.

It is important to validate the results from genome binning, although this approach has been repeatedly employed in various studies (Brown et al., 2015; Emerson et al., 2016). We suggest that validation using PhylopythiaS is suitable for genome bins of well-documented microbial taxa, such as Pseudomonas. But for novel phyla without many available genomes in public databases, comparing to single-cell amplified genomes (SAGs) from similar habitats is a more effective approach. The Thaumarchaeota genome binned in the present study possesses a gene inventory that is nearly identical to the SAGs from another Red Sea brine pool (Ngugi et al., 2015); this is strong evidence that the genome binning approach employed here is highly accurate.

The overall profiling of the genomics suggested the prevalence of organic carbon metabolism using diverse electron transfer elements. This notion is supported by the high diversity or specificity of particular pathways. Previous studies have discussed the distribution of F₄₂₀ biosynthesis genes in the genomes of microorganisms other than Archaea and Actinobacteria (Selengut and Haft, 2010; Li et al., 2014). In the present study, we showed a higher diversity of the F₄₂₀ biosynthesis gene cluster than previously assumed, suggesting a possible correlation between these genes and microbial adaptation to the brine pool. The presence of BMC compartments provides additional evidence that carbohydrate metabolism is particularly important for microbial diversity in this habitat, consistent with the concept that carbon storage is a mechanism conferring microbial resistance to hypersaline stresses (Werner et al., 2014). Notably, environmental adaptation of the Thuwal cold seep brine residents seems to have been facilitated by horizontal gene transfer (HGT) from various bacteria, as indicated by the phylogeny of several novel genes in the genome bins. The metatranscriptomic results suggested that anaerobic carbohydrate catabolism modulates not only genomic content but also post-genomic activity.

With respect to the detailed microbial lifestyles, this study is the first to report the presence of certain fermentative groups in a brine environment. To date, Microthrix have only been recorded in activated sludge ecosystems (Muller et al., 2012; McIlroy et al., 2013); Actinobacterium MedAcidi-G2B, for which little genomic information is available, have been reported only in the deep-sea habitat of the Mediterranean Sea (Mizuno et al., 2015); and Cloacimonetes are usually found in wastewater treatment plants or bioreactors (Pelletier et al., 2008; Nobu et al., 2015). Cloacimonetes and Marinimicrobia were found to be proteolytic amino-acid degraders and syntrophic propionate degraders, respectively (Nobu et al., 2015). However, the genomic and metatranscriptomic analyses of the present study showed their capacity to perform carbohydrate fermentation. The Bathyarchaeota may also be able to degrade polysaccharides, as suggested by the presence and expression of unique chitinase, cellulase and a rather complete acetate fermentation pathway in the MAG. The Bathyarchaeota MAG harbors a dissimilatory sulfate reduction pathway, differentiating it from other Bathyarchaeota members (Lloyd et al., 2013; Meng et al., 2014; Fillol et al., 2016) and suggesting the metabolic plasticity of this archaeal group. Nitrate or sulfate reducers, such as unclassified members of the Deltaproteobacteria and Epsilonproteobacteria, are also probable important players in polysaccharide fermentation, supporting the metabolic interconnection of carbohydrate metabolism with sulfur and nitrogen cycles in the brine pool.

Conclusion

Microbial consumption of organic carbons that sink from the surface ocean to deep-sea environments contributes to the remineralization of so called “dead organic matter” or “fecal material.” Although, organic carbon consumption by deep-sea bacteria and archaea has not been quantified, our findings demonstrated the participation of diverse and ubiquitous microbial taxa in this process. The niche-specific functions suggest a long history of microbial adaptation to the brine pool, where anaerobic organic carbon metabolism is likely to be an important driving force. Our findings shed new light on the roles of microorganisms in carbon cycling in deep-sea environments.

References

• 1

  AlbertsenM.HugenholtzP.SkarshewskA.NielsenK. L.TysonG. W. (2013). Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat. Biotechnol.31, 533–538. 10.1038/nbt.2579

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AntunesA.NgugiD. K.StinglU. (2011). Microbiology of the Red Sea (and other) deep-sea anoxic brine lakes. Environ. Microbiol. Rep.3, 416–433. 10.1111/j.1758-2229.2011.00264.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.19, 455–477. 10.1089/cmb.2012.0021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BatangZ. B.PapathanassioubE.Al-SuwailemaA.SmithbC.SalomidibM.PetihakisbG. M.et al. (2012). First discovery of a cold seep on the continental margin of the central Red Sea. J. Mar. Syst.94, 247–253. 10.1016/j.jmarsys.2011.12.004

  • CrossRef
  • Google Scholar

• 5

  BhaskarP. V.BhosleN. B. (2005). Microbial extracellular polymeric substances in marine biogeochemical processes. Curr. Sci.88, 45–53.

  • Google Scholar

• 6

  BrownC. T.HugL. A.ThomasB. C.SharonI.CastelleC. J.SinghA.et al. (2015). Unusual biology across a group comprising more than 15% of domain Bacteria. Nature523, 208–211. 10.1038/nature14486

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  CantarelB. L.CoutinhoP. M.RancurelC.BernardT.LombardV.HenrissatB. (2009). The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res.37, 233–238. 10.1093/nar/gkn663

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  CarverT. J.RutherfordK. M.BerrimanM.RajandreamM. A.BarrellB. G.ParkhillJ. (2005). ACT: the Artemis comparison tool. Bioinformatics21, 3422–3423. 10.1093/bioinformatics/bti553

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  CastresanaJ. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol.17, 540–552. 10.1093/oxfordjournals.molbev.a026334

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  ClaretC. (1998). A method based on artificial substrates to monitor hyporheic biofilm development. Int. Rev. Hydrobiol.83, 135–143. 10.1002/iroh.19980830204

  • CrossRef
  • Google Scholar

• 11

  D'ArcyR.AmendJ. P.OsburnM. R. (2013). Microbial diversity and potential for arsenic and iron biogeochemical cycling at an arsenic rich, shallow-sea hydrothermal vent (Tutum Bay, Papua New Guinea). Chem. Geol.348, 37–47. 10.1016/j.chemgeo.2012.02.024

  • CrossRef
  • Google Scholar

• 12

  EmersonJ. B.ThomasB. C.AlvarezW.BanfieldJ. F. (2016). Metagenomic analysis of a high carbon dioxide subsurface microbial community populated by chemolithoautotrophs and bacteria and archaea from candidate phyla. Environ. Microbiol.18, 1686–1703. 10.1111/1462-2920.12817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  FaberE.BotzR.PoggenburgJ.SchmidtM.StoffersP.HartmannM. (1998). Methane in Red Sea brines. Org. Geochem.29, 363–379. 10.1016/S0146-6380(98)00155-7

  • CrossRef
  • Google Scholar

• 14

  FerrerM.WernerJ.ChernikovaT. N.BargielaR.FernándezL.La ConoV.et al. (2012). Unveiling microbial life in the new deep-sea hypersaline Lake Thetis. Part II: a metagenomic study. Environ. Microbiol.14, 268–281. 10.1111/j.1462-2920.2011.02634.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  FillolM.AuguetJ. C.CasamayorE. O.BorregoC. M. (2016). Insights in the ecology and evolutionary history of the Miscellaneous Crenarchaeotic Group lineage. ISME J.10, 665–677. 10.1038/ismej.2015.143

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  FormoloM. J.LyonsT. W. (2013). Sulfur biogeochemistry of cold seeps in the Green Canyon region of the Gulf of Mexico. Geochim. Cosmochim. Acta119, 264–285. 10.1016/j.gca.2013.05.017

  • CrossRef
  • Google Scholar

• 17

  GalperinM. Y.MakarovaK. S.WolfY. I.KooninE. V. (2014). Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res.43, D261–D269. 10.1093/nar/gku1223

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  GuanY.HikmawanT.AntunesA.NgugiD.StinglU. (2015). Diversity of methanogens and sulfate-reducing bacteria in the interfaces of five deep-sea anoxic brines of the Red Sea. Res. Microbiol.166, 688–699. 10.1016/j.resmic.2015.07.002

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  GuezennecJ.Ortega-MoralesO.RaguenesG.GeeseG. (1998). Bacterial colonization of artificial substrate in the vicinity of deep-sea hydrothermal vents. FEMS Microbiol. Ecol.26, 89–99. 10.1111/j.1574-6941.1998.tb00495.x

  • CrossRef
  • Google Scholar

• 20

  HusonD. H.MitraS.RuscheweyhH. J.WeberN.SchusterS. C. (2011). Integrative analysis of environmental sequences using MEGAN4. Genome Res.21, 1552–1560. 10.1101/gr.120618.111

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  HyattD.ChenG. L.LoCascioP. F.LandM. L.LarimerF. W.HauserL. J.et al. (2010). Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics11:119. 10.1186/1471-2105-11-119

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  JonesP. R.CottrellM. T.KirchmanD. L.DexterS. C. (2007). Bacterial community structure of biofilms on artificial surfaces in an estuary. Microb. Ecol.53, 153–162. 10.1007/s00248-006-9154-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  KanehisaM.GotoS. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res.28, 27–30. 10.1093/nar/28.1.27

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods9, 357–359. 10.1038/nmeth.1923

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  LiX.FengF.ZengY. (2014). Genome of betaproteobacterium Caenimonas sp. strain SL110 contains a coenzyme F420 biosynthesis gene cluster. J. Microbiol. Biotechnol.24, 1490–1494. 10.4014/jmb.1405.05039

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  LloydK. G.SchreiberL.PetersenD. G.KjeldsenK. U.LeverM. A.SteenA. D.et al. (2013). Predominant archaea in marine sediments degrade detrital proteins. Nature496, 215–218. 10.1038/nature12033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  McIlroyS. J.KristiansenR.AlbertsenM.KarstS. M.RossettiS.NielsenJ. L.et al. (2013). Metabolic model for the filamentous ‘Candidatus Microthrix parvicella’ based on genomic and metagenomic analyses. ISME J.7, 1161–1172. 10.1038/ismej.2013.6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  MengJ.XuJ.QinD.HeY.XiaoX.WangF. (2014). Genetic and functional properties of uncultivated MCG archaea assessed by metagenome and gene expression analyses. ISME J.8, 650–659. 10.1038/ismej.2013.174

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  MizunoC. M.Rodriguez-ValeraF.GhaiR. (2015). Genomes of planktonic Acidimicrobiales: widening horizons for marine Actinobacteria by Metagenomics. mBio6, e02083–e02014. 10.1128/mBio.02083-14

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  MullerE. E.PinelN.GilleceJ. D.SchuppJ. M.PriceL. B.EngelthalerD. M.et al. (2012). Genome sequence of ‘Candidatus Microthrix parvicella’ Bio17-1, a long-chain-fatty-acid-accumulating filamentous actinobacterium from a biological wastewater treatment plant. J. Bacteriol.194, 6670–6671. 10.1128/JB.01765-12

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  NgugiD. K.BlomJ.AlamI.RashidM.Ba-AlawiW.ZhangG.et al. (2015). Comparative genomics reveals adaptations of a halotolerant thaumarchaeon in the interfaces of brine pools in the Red Sea. ISME J.9, 396–411. 10.1038/ismej.2014.137

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  NobuM. K.DodsworthJ. A.MurugapiranS. K.RinkeC.GiesE. A.WebsterG.et al. (2016). Phylogeny and physiology of candidate phylum ‘Atribacteria’ (OP9/JS1) inferred from cultivation-independent genomics. ISME J.10, 273–286. 10.1038/ismej.2015.97

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  NobuM. K.NarihiroT.RinkeC.KamagataY.TringeS. G.WoykeT.et al. (2015). Microbial dark matter ecogenomics reveals complex synergistic networks in a methanogenic bioreactor. ISME J.9, 1710–1722. 10.1038/ismej.2014.256

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  ParksD. H.ImelfortM.SkennertonC. T.HugenholtzP.TysonG. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res.25, 1043–1055. 10.1101/gr.186072.114

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  PatelR. K.JainM. (2012). NGS QC Toolkit: a toolkit for quality control of next generation sequencing data. PLoS ONE7:e30619. 10.1371/journal.pone.0030619

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  PatilK. R.RouneL.McHardyA. C. (2012). The PhyloPythiaS web server for taxonomic assignment of metagenome sequences. PLoS ONE7:e38581. 10.1371/journal.pone.0038581

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  PelletierE.KreimeyerA.BocsS.RouyZ.GyapayG.ChouariR.et al. (2008). Candidatus Cloacamonas acidaminovorans: genome sequence reconstruction provides a first glimpse of a new bacterial division. J. Bacteriol.190, 2572–2579. 10.1128/JB.01248-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  RinkeC.SchwientekP.SczyrbaA.IvanovaN. N.AndersonI. J.ChengJ. F.et al. (2013). Insights into the phylogeny and coding potential of microbial dark matter. Nature499, 431–437. 10.1038/nature12352

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  SaltaM.WhartonJ. A.BlacheY.StokesK. R.BriandJ. F. (2013). Marine biofilms on artificial surfaces: structure and dynamics. Environ. Microbiol.15, 2879–2893. 10.1111/1462-2920.12186

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  SelengutJ. D.HaftD. H. (2010). Unexpected abundance of coenzyme F-420-dependent enzymes in Mycobacterium tuberculosis and other Actinobacteria. J. Bacteriol.192, 5788–5798. 10.1128/JB.00425-10

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol.30, 2725–2729. 10.1093/molbev/mst197

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  TianR. M.WangY.BaugouffaS.GaoZ. M.CaiL.BajicV.et al. (2014). Genomic analysis reveals versatile heterotrophic capacity of a potentially symbiotic sulfur oxidizing bacterium in sponge. Environ. Microbiol.16, 3548–3561. 10.1111/1462-2920.12586

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  van der WielenP. W.BolhuisH.BorinS.DaffonchioD.CorselliC.GiulianoL.et al. (2005). The enigma of prokaryotic life in deep hypersaline anoxic basins. Science307, 121–123. 10.1126/science.1103569

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  WangQ.GarrityG. M.TiedjeJ. M.ColeJ. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol.73, 5261–5267. 10.1128/AEM.00062-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  WangY.CaoH.ZhangG.BougouffaS.LeeO. O.Al-SuwailemA.et al. (2013). Autotrophic microbe metagenomes and metabolic pathways differentiate adjacent Red Sea brine pools. Sci. Rep.3:1748. 10.1038/srep01748

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  WernerJ.FerrerM.MichelG.MannA. J.HuangS.JuarezS.et al. (2014). Halorhabdus tiamatea: proteogenomics and glycosidase activity measurements identify the first cultivated euryarchaeon from a deep-sea anoxic brine lake as potential polysaccharide degrader. Environ. Microbiol.16, 2525–2537. 10.1111/1462-2920.12393

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  ZhangW.WangY.BougouffaS.TianR.CaoH.LiY.et al. (2015). Synchronized dynamics of bacterial niche-specific functions during biofilm development in a cold seep brine pool. Environ. Microbiol.17, 4089–4104. 10.1111/1462-2920.12978

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  ZhangW.WangY.LeeO. O.TianR.CaoH.GaoZ.et al. (2013). Adaptation of intertidal biofilm communities is driven by metal ion and oxidative stresses. Sci. Rep.3:3180. 10.1038/srep03180

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  ZhangW. P.WangY.TianR. M.BougouffaS.YangB.CaoH. L.et al. (2014). Species sorting during biofilm assembly by artificial substrates deployed in a cold seep system. Sci. Rep.4:6647. 10.1038/srep06647

  • Pubmed Abstract
  • CrossRef
  • Google Scholar