Microbial community changes along the active seepage site of one cold seep in the Red Sea

The active seepage of the marine cold seeps could be a critical process for the exchange of energy between the submerged geosphere and the sea floor environment through organic-rich fluids, potentially even affecting surrounding microbial habitats. However, few studies have investigated the associated microbial community changes. In the present study, 16S rRNA genes were pyrosequenced to decipher changes in the microbial communities from the Thuwal seepage point in the Red Sea to nearby marine sediments in the brine pool, normal marine sediments and water, and benthic microbial mats. An unexpected number of reads from unclassified groups were detected in these habitats; however, the ecological functions of these groups remain unresolved. Furthermore, ammonia-oxidizing archaeal community structures were investigated using the ammonia monooxygenase subunit A (amoA) gene. Analysis of amoA showed that planktonic marine habitats, including seeps and marine water, hosted archaeal ammonia oxidizers that differed from those in microbial mats and marine sediments, suggesting modifications of the ammonia oxidizing archaeal (AOA) communities along the environmental gradient from active seepage sites to peripheral areas. Changes in the microbial community structure of AOA in different habitats (water vs. sediment) potentially correlated with changes in salinity and oxygen concentrations. Overall, the present results revealed for the first time unanticipated novel microbial groups and changes in the ammonia-oxidizing archaea in response to environmental gradients near the active seepages of a cold seep.

Introduction

Cold seeps mainly occur in geologically active and passive continental margins, and they transport dissolved and gaseous phase compounds to the ocean to sustain significant chemosynthetic biomass by providing bioactive reductants, sulfides, methane and hydrogen (Levin, 2005; Suess, 2010). Many early studies have focused mainly on carbon and sulfur cycling in these specialized ecosystems, with a specific focus on the anaerobic oxidation of methane (AOM) coupled to sulfate reduction (SR) in microbes in hypersaline cold seep sediments (Maignien et al., 2013). In particular, AOM in a cool seep environment has been extensively studied using metagenomic and metatranscriptomic methods (Stokke et al., 2012), while changes in the microbial composition along the environmental gradient near the seepage sites have only been marginally assessed.

In cold seep ecosystems, microbes must rapidly transform carbon, sulfur and nitrogen compounds. Nitrogen fixation plays important roles in cold seeps (Pernthaler et al., 2008; Dang et al., 2009; Dekas et al., 2009; Miyazaki et al., 2009), and high rates of nitrogen removal due to denitrification in cold seep sediments have been proposed (Bowles and Joye, 2011). The diversity and abundance of anaerobic ammonium oxidizing (anammox) bacteria in cold seep hydrocarbon-rich fluids have been reported (Russ et al., 2013; Shao et al., 2014). A reduced diversity and abundance of the ammonia oxidizing archaea (AOA) thaumarchaea were found in cold seep sediments in the Okhotsk Sea (Dang et al., 2010) and northeastern Japan Sea (Nakagawa et al., 2007). All of these studies have shed some light on the possible contributions of various types of microorganisms to nitrogenous nutrient recycling; however, the ecological functions of the thaumarchaea in cold seeps remain largely unresolved.

The oxidation of ammonia is the first and rate-limiting step in the process of nitrification (Kowalchuk and Stephen, 2001) performed by bacterial and archaeal groups. Ammonia monooxygenase (amo) catalyzes the oxidation of ammonia and consists of several subunits, among which the amoA gene encoding subunit A has been widely used as a reliable genetic marker to explore the diversity and abundance of AOA (AOA) in diverse ecosystems (Junier et al., 2010; Cao et al., 2011b). Various environmental parameters, such as pH, depth, nutrients, and dissolved oxygen, have been identified as potential factors determining the dominant ammonia oxidizer phylotypes and their diversity in ecosystems (Erguder et al., 2009; Cao et al., 2013). However, the ecological role of the AOA in cold seeps remains unexplored.

The Thuwal Seeps is a cold brine seep system located at a depth of ∼850 m and was first discovered on May 7, 2010 by a remotely operated vehicle (ROV) on the Saudi continental margin of the central Red Sea during a survey built into the framework of “KAUST Red Sea Expedition Spring 2010” (Batang et al., 2012). The seep is located at the base of a steep rocky wall that is closer to the shore (20 km) than to the axial trough (120 km). In fact, active brine ventings have been observed at two seep sites, which are named Thuwal Seeps I and II (I is located at 22° 17.3′ N–38° 53.8′ E; II is located at 22° 16.9′ N–38° 53.9′ E). A shallow brine pool was formed by fluids from seeps with a low temperature (21.7°C) and salinity (74‰) compared with the other brine pools in the Red Sea (Batang et al., 2012). Although the hypersaline brine pools at the Thuwal Seeps are harsh to organisms, high biomass production was observed (Batang et al., 2012). Brine waters likely originate from evaporitic deposits of submarine geological formations that flow from the faulting system at the base of the rocky scarp where the Thuwal brine pool formed. Extended chemosynthetic bacterial mats and dense aggregations of live and dead organisms have been observed (spatangoid urchins, anemones, serpulid tubeworms, sponges, clams, fishes, crabs and shrimps) (Batang et al., 2012). Thus, the Thawal seeps provides a good opportunity to understand the responses of microbes that reside near the seepage sites of cold seeps.

To explore potential changes in the microbial communities from the seepage site to nearby areas in this cold seep ecosystem, habitats including marine sediments, marine water, seep water, and microbial mats were sampled for 16S rRNA gene analysis. Because high concentrations of inorganic nitrogen species were detected in the environment, the AOA community structures were also examined based on the amoA gene. The results showed that the habitat type dictated the community structure, while the environmental gradient shaped the changes in the AOA community from the active seepage site to peripheral areas.

Materials and Methods

Sampling and Environmental Parameter Measurements

Field sampling was conducted in November 2011 in the Thuwal cold Seep II (22^o16’N–38^o53’E) via the ROV Max Rover developed by Deep Sea Systems International (DSSI), USA, during the KAUST Red Sea exploration cruise (Figure 1). The venting site of the seepage was only about 1.5 m wide and 1.0 m deep, while the brine pool was very shallow with a depth of approximately 1.0 m in most places. Four types of habitats around the brine pool approximating the seeping vents were sampled. The samples from the seep vent (Seep4), normal marine water (TS06W, normal marine water overlaying the pool), sediments (TS03S and TS06S from outside the brine pool; TS08S from inside the brine pool), and a microbial mat on the bank of the pool (top of the microbial mat) have been described in a previous study (Wang et al., 2014b). Two replicate each were sampled from the seep and marine water for 16S rRNA gene analysis.

FIGURE 1

FIGURE 1. Diagram of the sampling sites in the Red Sea modified from Batang et al. (2012). (A) Locations of brine pools in the Red Sea axial zone (blue shade); (B) Location of Thuwal Seeps (dark blue patch) near coral reefs (pink patches) along the Saudi coast; (C) Outline of Thuwal Seeps brine pool (thick line) superimposed on depth contours (Batang et al., 2012).

The in situ physicochemical parameters (temperature, salinity, and concentration of dissolved oxygen) were measured as described previously (Lee et al., 2014; Wang et al., 2014b). The concentrations of dissolved organic carbon (DOC) and total nitrogen (TN) in the water sample and pore-water from the sediment were measured using the combustion method (Trichet et al., 2001), while the concentrations of ammonium and nitrate plus nitrite were determined separately using a TNM-I analyzer (Simadzu, Kyoto, Japan). Detailed descriptions of the sampling sites and their physicochemical characteristics are provided in Figure 1 and Table 1.

TABLE 1

TABLE 1. Physicochemical parameters of the profiles from the cold seep to the far sediments in the Red Sea (some samples and parameters were adopted from Wang et al., 2014b).

Molecular Experiments

Genomic DNA was extracted from all of the sediment samples (collected in triplicate) using the PowerMax Soil DNA Isolation Kit (Epicentre Biotechnologies, Madison, WI, USA) according to the manufacturer’s instructions. For microbial mat and water filters, the modified SDS-based method described by Lee et al. (2009) was employed. The DNA quality and quantity were checked using a PicoGreen dsDNA quantitation kit (Life Technologies, Carlsbad, CA, USA) and gel electrophoresis.

The amoA gene was cloned into a plasmid, and 16S rRNA gene sequencing was conducted by pyrosequencing. PCR amplifications using the isolated genomic DNA as template were performed in 20-μl reactions consisting of 1.25 U of Taq DNA polymerase (New England Biolabs, England), 2 μl of 10× PCR buffer (15 mM Mg²⁺), 1 μl of deoxynucleoside triphosphates (dNTPs) (2.5 mM), 2 μl (∼30 ng) of DNA template, and 0.5 μl of each primer (10 μM) (Arch-amoAF: 5′-STAATGGTCTGGCTTAGACG-3′ and Arch-amoAR: 5′-GCGGCCATCCATCTGTATGT-3′) (Francis et al., 2005). The amplicons resulting from PCR performed in triplicate for each sample were pooled and used to construct clone libraries. The amplified 635-bp bands were excised from 1.0% (wt/vol) agarose gels and gel-purified using the Universal DNA Purification Kit (Tiangen Biotech, China). The purified PCR products were cloned into the pCR2.1-TOPO vector using the TOPO TA cloning kit for sequencing with One Shot TOP10 competent cells according to the manufacturer’s instructions (Invitrogen, USA). Each clone was randomly selected and screened by PCR with primers M13F and M13R to select positive clones. The positive clones were sequenced with the vector-specific primers M13F and M13R using the Applied Biosystems 3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA).

Regarding the 16S rRNA gene pyrosequencing, the hyper variable regions V4 to V8 of the bacterial and archaeal 16S rRNA gene were amplified for each sample using the universal forward primer U515F (5′- GTGYCAGCMGCCGCGGTAA -3′) and the reverse primer R1492 (5′- GACGGGCGGTGTGTRCAA -3′) (Wang et al., 2014a; Zhang et al., 2014) with unique 6-nucleotide (nt) barcodes. Five units of Pfu Turbo DNA polymerase (Stratagene, La Jolla, CA, USA), 1x Pfu reaction buffer, 0.2 mM dNTPs (TaKaRa, Dalian, China), 0.1 μM each barcoded primer, and 20 ng of genomic DNA template were mixed in a 100-μl PCR volume. The PCR procedure included an initial denaturation at 94°C for 5 min, 26 cycles at 94°C for 30 s, 53°C for 30 s, and 72°C for 45 s, and final extension at 72°C for 6 min using a Bio-Rad thermal cycler (MJ Research Inc., Bio-Rad). The PCR products were purified using a TaKaRa Agarose Gel DNA Purification Kit (TaKaRa, Dalian, China) and quantified with a NanoDrop device. Two hundred nanograms of each purified 16S amplicon were mixed and then subjected to pyrosequencing using the Roche 454 FLX Titanium platform at the Chinese National Human Genome Centre in Shanghai, China.

Sequencing Analysis for 16S rRNA Gene Pyrosequencing Data

The pyrosequencing data were deposited in the NCBI Sequence Read Archive (SRA) database. The downstream bioinformatics analysis was performed using QIIME 1.7.4 (Caporaso et al., 2010b). The criteria used for quality control of all data have been described in a previous study (Lee et al., 2012).

Reads were assigned to their respective samples according to their barcodes and then subjected to a second round of quality control using Denoiser (Reeder and Knight, 2010). Qualified reads were clustered using uclust (Edgar, 2010) and then assigned to operational taxonomic units (OTUs) at a similarity of 97%. Representatives of the most abundant reads were selected from each OTU for subsequent analysis. Representative OTUs were aligned de novo using MUSCLE (Edgar, 2004), and a phylogenetic tree was produced using FastTree (Price et al., 2009). Representative OTUs were also aligned using PyNAST (Caporaso et al., 2010a) with the Silva108 database as a reference. Successfully aligned reads were submitted to ChimeraSlayer (Haas et al., 2011) to identify and discard chimeric reads. Species diversity, richness, and rarefaction curves were computed at a similarity of 97% as part of QIIME’s alpha diversity pipeline. Beta diversity analysis was conducted after rarefying the samples in the smallest library using QIIME. A step size of 100 was used with 100 repetitions at each step. Taxonomic assignment was conducted using the Ribosomal Database Project (RDP) classifier version 2.2 (Wang et al., 2007) against Silva108 (Pruesse et al., 2007) with a bootstrap confidence level of 50%. The number of reads assigned to the different genera were converted into percentages.

Diversity Indice Calculations for the amoA Gene

The DNA sequences of the amoA genes were transformed into MEGA 6.0 software (Tamura et al., 2013) and aligned using CLUSTALW in MEGA. OTUs were identified at a similarity level of 0.97 using the MOTHUR program (Schloss et al., 2009). Based on the OTU assignments, the nonparametric abundance estimators, Chao1 and ACE, and the diversity estimates (ACE) were calculated using summary.single in MOTHUR. Rarefaction curves were calculated for all samples at a distance cutoff of 0.03 using rarefaction.single (Schloss et al., 2009).

Phylogenetic Analysis of the amoA Gene

The selected OTUs from all the samples were imported into MEGA and translated into amino acids. Reference sequences from GenBank were downloaded to construct a comprehensive database and then clustered at an amino acid similarity level of 0.97. After merging with the reference sequences, including those from enrichments, the initial phylogenetic tree was constructed based on the neighbor-joining (NJ) algorithm with 1,000 bootstrap replicates using MEGA software (Tamura et al., 2013). The NJ phylogenetic tree was employed initially to construct the maximum likelihood (ML) tree in MEGA.

Microbial Community Analyses for the amoA Gene

To compare the phylogenetic diversity between different habitats, a genetic distance matrix of the sequences from each habitat was developed. Principal coordinate analysis (PCoA) and jackknife environment cluster analysis were conducted using the online software Fast UniFrac (Hamady et al., 2010), which utilizes the genetic distance matrix based on the gene sequence data. In addition, hierarchical clustering analysis (UPGMA algorithm with jackknife supporting values) was used. The environmental cluster tree was projected in MEGA 6.0 (Tamura et al., 2013).

Non-metric multidimensional scaling (NMDS) was conducted in MOTHUR based on a genetic cutoff distance of 0.03 for amoA to determine the similarity between samples. The command was modeled based on the function using the majorization algorithm (Borg and Groenen, 1997).

One-way ANOSIM methods with 999 permutations were performed in MOTHUR to test for the significance of differences in community composition between the clone libraries. Simultaneously, LIBSHUFF statistical comparisons were conducted using LIBSHUFF in MOTHUR.

Nucleotide Sequence Accession Numbers

All of the archaeal amoA gene sequences retrieved in the present study have been deposited in the GenBank database at NCBI under accession numbers KM109433–KM109966, and the 16S rRNA gene pyrosequencing data have been submitted to the NCBI SRA database under accession number SRX501833.

Results

Sampling Site Descriptions and Physiochemical Parameters

The features of our sampling sites have been documented in a previous study (Wang et al., 2014b), and the sampling map is also shown in Figure 1. Higher concentrations of inorganic nitrogen compounds have been observed in the seep and site TS08S compared with sites TS03S and TS06 (Table 1). High concentrations of hydrogen sulfide (>200 μmol/L) were also recorded in the nearly anoxic brine pool (Wang et al., 2014b).

16S rRNA Gene Diversity and Composition

A total of 47,353 reads (∼4304 reads per sample) were passed through the quality check. The number of OTUs and estimated species richness at the 3% dissimilarity level are listed (Supplementary Figure S1). The highest number of OTUs was found on top of the microbial mat followed by the marine sediments (Supplementary Figure S1).

At the 50% confidence threshold, qualified reads could be assigned to eighteen phyla based on the analyses using the QIIME pipeline (Supplementary Figure S2). The proportions of these phyla varied among different habitats. For example, 12.70% of the Proteobacteria was observed in TS06W-1. In contrast, Proteobacteria demonstrated a low abundance in marine sediments but were exceptionally enriched in the microbial mat. Marine sediments contained a higher abundance of Euryarchaeota than the seeps and marine water (<0.25%). The Thaumarchaea was an over-represented phylum in seeps and marine water compared with marine sediments and the microbial mat (Supplementary Figure S2). Nitrospirae were highly abundant in two marine sediment samples, TS08S and TS09S, whereas a high abundance of Deferribacteres was detected in the seeps and marine water (Supplementary Figure S2).

Substantial variations in the microbial communities associated with different habitats were detected down to the genus level. For example, the microbial communities from marine water and seeps were dominated by unclassified Thaumarchaea, ranging from 19.27% in Seep4-2 to 62.50% in TS06W-2 (Figure 2). In contrast, marine sediment TS09S and the top of the microbial mat contained proportions of unclassified Thaumarchaeota of 0.34% and 2.17%, respectively (Figure 2). In cold seeps and marine water, Thermoplasmata, Sar406 clade (Deferribacteres), Nitrospira, Sar324 clade and the E01-9c-26 marine group also exhibited a high abundance. In particular, the abundance of Sar406 was rather high in the seeps (Figure 2). Group_c3 in the Crenarchaeota, a miscellaneous crenarchaeotic group, marine benthic group d (MBGD) and dhveg-1 (currently designated as Thermoplasmata) were found in microbial mats and marine sediments. Opb95 (Nitrospira) was the dominant species group in two of the sediment samples from the brine pool (TS08S and TS09S) (Figure 2). The miscellaneous crenarchaeotic group was dominant in normal marine sediment samples (TS03S and TS06S) (Figure 2). In addition, TS03S was dominated by MBGD and the dhveg-1 (Thermoplasmata) group. At the genus level, another two groups, Nkb17 (Holophagae) and Rb25 (Acidobacteria), were also observed with a higher abundance in microbial mats than in the other samples.

FIGURE 2

FIGURE 2. Taxonomic classification at the genus level of microbial reads retrieved from different habitats in the Red Sea based on 16S rRNA gene pyrosequencing data presented as the relative abundance (MMT and MMB denote the top and bottom of the microbial mat, respectively).

amoA Gene Diversity

Rarefaction curves (Figure 3) and diversity indices (Table 2) were determined for each clone library. The results revealed that the AOA diversity was far from exhaustively sampled, in particular for the microbial mat sample (Figure 3). The highest diversity indices were found on the top of the microbial mat and in the marine sediments, which indicated that the sediment may represent the largest reservoir of AOA diversity (Table 2). In contrast, the diversity indices were low in water samples from seeps, brine, and normal bottom seawater (Table 2).

FIGURE 3

FIGURE 3. Rarefaction cures for amoA genes sequences from each sample based on a cutoff of 0.03 generated by MOTHUR.

TABLE 2

TABLE 2. Diversity indices for the amoA genes in all samples from the Red Sea.

Phylogenetic Tree Based on the amoA Gene

The phylogenetic analysis of the amoA gene revealed three major monophyletic clusters (i.e., Nitrosopumilus, Nitrosotalea, Nitrosocaldus) and a non-monophyletic cluster (Nitrososphaera) that comprised mostly genes from the soil and sediment environments (Pester et al., 2012; Cao et al., 2013). In addition, one subcluster, Nitrosoarchaeum within the Nitrosopumilus cluster, was characterized by organisms that survive in low salinity habitats. In the present study, only two sequences (TS03S-85 and TS06-52) were classified in the Nitrososphaera cluster, while both the Nitrosotalea and Nitrosocaldus clusters were absent (Figure 4). In contrast, clusters of Nitrosopumilus and Nitrosoarchaeum were dominant in this environment (Figure 4). Although the statistical supports for most of the nodes were weak, based on two criteria, ML and NJ, the topologies were remarkably consistent with one another. This result supported the stability of the topology, and subsequently, only ML phylogeny was presented (Figure 4).

FIGURE 4

FIGURE 4. Phylogenetic tree reconstructed from the deduced AmoA protein sequences using the maximum likelihood (ML) criterion. AmoA sequences from the present study are shown in bold, and sample names of representative sequences within the same operational taxonomic units (OTUs) are bracketed.

Most of the amoA gene sequences determined in the present study was distributed into the Nitrosopumilus cluster, with the highest contribution from marine sediments and the top of the microbial mat. Within the Nitrosopumilus cluster, most of the subclusters had low support values (Figure 4). We downloaded all of the archaeal amoA gene sequences from GenBank to create one database (updated to November, 2013) and would assign the gene sequences from the present study to the closest relatives. However, the closest relatives to most of the sequences from the present study affiliated with the Nitrosopumilus cluster could not be identified (Figure 4). In particular, the sequences from the top of the microbial mat did not cluster with the sequences from GenBank (Figure 4). This over-dispersion of sequences from marine sediments and the microbial mat in the Nitrosopumilus cluster might be explained by the influence of the chemocline in the sampling area. Sequences from the microbial mat exhibited a higher diversity than those from other habitats and were dispersed throughout the whole phylogenetic tree, although a low abundance of Thaumarchaea was observed in these two samples based on 16S rRNA gene analysis (Figure 2). In addition, one lineage with a long branch included a partial sequence of TS06W from normal marine water, and most sequences from the seep clustered with those from marine water in the Gulf of Mexico (Tolar et al., 2013) and East China Sea (Hu et al., 2011) (Figure 4).

The subcluster (designated as the Nitrosoarchaeum subcluster) in the Nitrosopumilus cluster has been previously proposed to reside in a Low Salinity Environment Cluster that includes Candidatus Nitrosoarchaeum limnia (Mosier et al., 2012) and Candidatus Nitrosoarchaeum koreensis (Kim et al., 2011), which were isolated from low salinity environments (Cao et al., 2013). Surprisingly, in the present study, the sequences from all of the high salinity samples were embedded in this subcluster (Figure 4).

Phylogenetic Ecology of the amoA Gene

Five types of habitats of archaeal ammonia oxidizers were examined herein (Figure 5). AOA communities were more similar within habitats than among habitats, as deduced from NMDS, UniFrac, Anosim, and Libshuff analyses (Figures 5 and 6; Table 3). For example, the clone libraries of the amoA gene from marine water shared more similarities with those from the cold seep and clustered together, as supported by UniFrac analyses (Figure 6). Based on the OTUs defined in each sample, NMDS plotted almost all of the clones from the seeps and marine water together. These samples were separated from those retrieved from marine sediments and the microbial mat (Figure 5). Although the water from the cold seeps had a higher salinity than normal marine water, similarities were observed among the communities of archaeal ammonia oxidizers between the cold seep and normal marine water (Figures 5 and 6).

FIGURE 5

FIGURE 5. Non-metric multidimensional scaling (NMDS) function analysis of all amoA genes in the present study. Each point represents one OTU of the amoA gene.

FIGURE 6

FIGURE 6. Principal component analysis (PCA) and jackknife environment cluster analysis of all amoA genes from various microbial habitats in the Red Sea.

TABLE 3

TABLE 3. Anosim and Libshuff analyses to compare the similarity of archaeal ammonia oxidizer communities between two clone libraries based on amoA gene sequences.

Discussion

An interaction between the fluid composition and microbes in the cold seep environment has been observed in some active cold seeps (Levin, 2005; Suess, 2010). The responses of microbes, such as microbial nitrogen utilizers, to the organic-rich fluids, however, has not been well investigated. In the present study, we confirmed the previously described presence of groups of archaeal ammonia oxidizers in a cold seep (Dang et al., 2010) based on the amoA gene and revealed clear changes in microbial communities along the environmental gradient from the seep vent to the cold seep brine pool and its surroundings (Figures 2 and 4).

Based on the 16S rRNA gene analysis, most of the reads from the microbial habitats could be sorted into known phyla. At the genus level, however, a variety of uncultured groups were assigned as endemic in different habitats in the Thuwal seeps, indicating a unique repertoire of novel microbial lineages. The potential ecological function of these groups is unclear, and further investigations using cultures or single cell genomics methods are necessary. For example, the OPB95 group affiliated with the candidate division OP8 was identified as the dominant group in marine sediment from the brine pool (Hugenholtz et al., 1998). MBGD and dhveg-1 (Thermoplasmata) were relatively abundant in microbial mats and in normal marine sediments in the present study. A previous study has already described the role of MBGD in protein remineralization in anoxic marine sediments (Lloyd et al., 2013). Altogether, these uncultured microbial groups support the need for additional functional studies.

The transition of archaeal ammonia oxidizers from the seep vent to the surrounding microbial mats, marine sediments, and water was clearly observed. Although relatively high concentrations of ammonia were recorded in the brine sediments, a high abundance of ammonia oxidizers was not observed. This result could be attributed to the low oxygen concentration; archaeal ammonia oxidizers still require oxygen as an electron acceptor. Another potential explanation could be a strong effect of salinity in determining AOA community structure patterns. Previous studies of prokaryotic phylogenies have revealed a clear separation between freshwater and marine lineages (Lozupone and Knight, 2007; Auguet et al., 2010; Cao et al., 2013), suggesting that salinity is one of the most important evolutionary barriers preventing frequent environmental transitions (Logares et al., 2009). In the present study, the various levels of salinity in different areas of the brine pools and cold seeps (Wang et al., 2014b) may have had a great influence on the evolution of AOA. In addition to the salinity, the concentrations of inorganic nitrogen-related compounds could represent another influential environmental parameter because the diversity indices changed along the concentration gradient of ammonia from the seep vent to the surrounding areas.

A lower abundance but a high diversity of unclassified Thaumarchaea was detected in the marine sediments and microbial mats, consistent with the findings of previous studies (Dang et al., 2010; Cao et al., 2011a, 2012). However, the marine water and seep vent were largely composed of uncultured Thaumarchaea with a low diversity but a high abundance, which is in agreement with the low diversity indices observed for the amoA sequences of marine plankton in a recent global-scale study (Cao et al., 2013). It is possible that habitats with low concentrations of substrate contain more diverse archaeal ammonia oxidizers, while habitats with high concentrations of substrate have a low diversity but a high abundance of this group. One interesting observation was that the amoA gene sequences from the normal marine sediments and the top of the microbial mat in the Red Sea shared a high similarity with sequences from the South China Sea but with limited sequences from other marine sediments (Cao et al., 2011a, 2012; Dang et al., 2013) (Figure 4). The ecological significance and allopatric distribution of this group of archaeal ammonia oxidizers in the marine sediments of Red Sea and South China Sea remains enigmatic and should be assessed in future studies.

In general, the advective supply of methane leads to dense microbial communities with high metabolic rates and anaerobic methane oxidation that is presumably coupled to SR and thus facilitates the formation of carbonates and generates extremely high concentrations of hydrogen sulfide in pore waters (Stokke et al., 2012; Maignien et al., 2013; Vigneron et al., 2013). Nonetheless, methane oxidizers were less abundant in the present study, which was in stark contrast to the typical cold seeps. However, fairly low concentrations of methane in the Thuwal seeping water (unpublished results) might explain the absence of microbes that oxidize methane in this cold seep ecosystem. The reduced abundance of methane oxidizers and scarcity of methane support the ubiquity of the cold seep in the Red Sea.

Overall, our analysis provides further insight regarding the changes in the microbial communities that reside in different habitats of the cold seep ecosystem in the Red Sea. We found (1) an unexpected number of unclassified microbial groups that contained several groups from candidate divisions with unknown ecological functions, (2) community structures that tended to be grouped by habitat types and salinity and oxygen concentrations appeared to be the driving force in shaping community structures, and (3) archaeal ammonia oxidizers in seeps and marine water differed from those in microbial mats and two types of marine sediments, consistent with many previous studies.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was supported by the National Basic Research Program of China (973 Program, No: 2012CB417304), the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB06010100 and XDB06010200), Exploration and Evaluation of Potential Resources from Deep-sea Microorganisms (DY125-15-R-02) from China Ocean Mineral Resources R&D Association (COMRRDA12SC02), and an award (SA-C0040/UK-C0016) from the King Abdullah University of Science and Technology in Saudi Arabia to P-YQ. The authors are grateful to Alex Shek for his technical assistance. We thank Dr. Abdulaziz M Al-Suwailem and his team from the King Abdullah University of Science and Technology for organizing research cruise and assistance in collecting samples.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2015.00739

References

Auguet, J. C., Barberan, A., and Casamayor, E. O. (2010). Global ecological patterns in uncultured Archaea. ISME J. 4, 182–190. doi: 10.1038/ismej.2009.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Batang, Z. B., Papathanassiou, E., Al-Suwailem, A., Smith, C., Salomidi, M., Petihakis, G., et al. (2012). First discovery of a cold seep on the continental margin of the central Red Sea. J. Marine Syst. 94, 247–253. doi: 10.1016/j.jmarsys.2011.12.004

CrossRef Full Text | Google Scholar

Borg, I., and Groenen, P. J. F. (1997). Modern Multidimensional Scaling: Theory and Applications. (New York: Springer). doi: 10.1007/978-1-4757-2711-1

CrossRef Full Text | Google Scholar

Bowles, M., and Joye, S. (2011). High rates of denitrification and nitrate removal in cold seep sediments. ISME J. 5, 565–567. doi: 10.1038/ismej.2010.134

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, H., Auguet, J. C., and Gu, J. D. (2013). Global ecological pattern of ammonia-oxidizing archaea. PLoS ONE 8:e52853. doi: 10.1371/journal.pone.0052853

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, H., Hong, Y., Li, M., and Gu, J. D. (2012). Community shift of ammonia-oxidizing bacteria along an anthropogenic pollution gradient from the Pearl River Delta to the South China Sea. Appl. Microbiol. Biotechnol. 94, 247–259. doi: 10.1007/s00253-011-3636-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, H., Hong, Y., Li, M., and Gu, J. D. (2011a). Phylogenetic diversity and ecological pattern of ammonia-oxidizing archaea in the surface sediments of the western Pacific. Microb. Ecol. 62, 813–823. doi: 10.1007/s00248-011-9901-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, H., Li, M., Dang, H., and Gu, J. D. (2011b). Responses of aerobic and anaerobic ammonia/ammonium-oxidizing microorganisms to anthropogenic pollution in coastal marine environments. Methods Enzymol. 496, 35–62. doi: 10.1016/B978-0-12-386489-5.00002-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Bittinger, K., Bushman, F. D., DeSantis, T. Z., Andersen, G. L., and Knight, R. (2010a). PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26, 266–267. doi: 10.1093/bioinformatics/btp636

PubMed Abstract | CrossRef Full Text | Google Scholar

Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010b). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, H., Luan, X., Chen, R., Zhang, X., Guo, L., and Klotz, M. G. (2010). Diversity, abundance and distribution of amoA-encoding archaea in deep-sea methane seep sediments of the Okhotsk Sea. FEMS Microbiol. Ecol. 72, 370–385. doi: 10.1111/j.1574-6941.2010.00870.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, H., Luan, X., Zhao, J., and Li, J. (2009). Diverse and novel nifH and nifH-like gene sequences in the deep-sea methane seep sediments of the Okhotsk Sea. Appl. Environ. Microbiol. 75, 2238–2245. doi: 10.1128/AEM.02556-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, H., Zhou, H., Yang, J., Ge, H., Jiao, N., Luan, X., et al. (2013). Thaumarchaeotal signature gene distribution in sediments of the northern South China Sea: an indicator of the metabolic intersection of the marine carbon, nitrogen, and phosphorus cycles? Appl. Environ. Microbiol. 79, 2137–2147. doi: 10.1128/AEM.03204-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Dekas, A. E., Poretsky, R. S., and Orphan, V. J. (2009). Deep-sea archaea fix and share nitrogen in methane-consuming microbial consortia. Science 326, 422–426. doi: 10.1126/science.1178223

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461

PubMed Abstract | CrossRef Full Text | Google Scholar

Erguder, T. H., Boon, N., Wittebolle, L., Marzorati, M., and Verstraete, W. (2009). Environmental factors shaping the ecological niches of ammonia-oxidizing archaea. FEMS Microbiol. Rev. 33, 855–869. doi: 10.1111/j.1574-6976.2009.00179.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Francis, C. A., Roberts, K. J., Beman, J. M., Santoro, A. E., and Oakley, B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. U.S.A. 102, 14683–14688. doi: 10.1073/pnas.0506625102

PubMed Abstract | CrossRef Full Text | Google Scholar

Haas, B. J., Gevers, D., Earl, A. M., Feldgarden, M., Ward, D. V., Giannoukos, G., et al. (2011). Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21, 494–504. doi: 10.1101/gr.112730.110

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamady, M., Lozupone, C., and Knight, R. (2010). Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 4, 17–27. doi: 10.1038/ismej.2009.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, A., Jiao, N., Zhang, R., and Yang, Z. (2011). Niche partitioning of marine group I Crenarchaeota in the euphotic and upper mesopelagic zones of the East China Sea. Appl. Environ. Microbiol. 77, 7469–7478. doi: 10.1128/AEM.00294-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Hugenholtz, P., Pitulle, C., Hershberger, K. L., and Pace, N. R. (1998). Novel division level bacterial diversity in a Yellowstone hot spring. J. Bacteriol. 180, 366–376.

PubMed Abstract | Google Scholar

Junier, P., Molina, V., Dorador, C., Hadas, O., Kim, O. S., Junier, T., et al. (2010). Phylogenetic and functional marker genes to study ammonia-oxidizing microorganisms (AOM) in the environment. Appl. Microbiol. Biotechnol. 85, 425–440. doi: 10.1007/s00253-009-2228-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, B. K., Jung, M. Y., Yu, D. S., Park, S. J., Oh, T. K., Rhee, S. K., et al. (2011). Genome sequence of an ammonia-oxidizing soil archaeon, “Candidatus Nitrosoarchaeum koreensis” MY1. J. Bacteriol. 193, 5539–5540. doi: 10.1128/JB.05717-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowalchuk, G. A., and Stephen, J. R. (2001). Ammonia-oxidizing bacteria: a model for molecular microbial ecology. Annu. Rev. Microbiol. 55, 485–529. doi: 10.1146/annurev.micro.55.1.485

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, O. O., Wang, Y., Tian, R., Zhang, W., Shek, C. S., Bougouffa, S., et al. (2014). In situ environment rather than substrate type dictates microbial community structure of biofilms in a cold seep system. Sci. Rep. 4, 3587. doi: 10.1038/srep03587

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, O. O., Wong, Y. H., and Qian, P. Y. (2009). Inter- and intraspecific variations of bacterial communities associated with marine sponges from san juan island, washington. Appl. Environ. Microbiol. 75, 3513–3521. doi: 10.1128/AEM.00002-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, O. O., Yang, J. K., Bougouffa, S., Wang, Y., Batang, Z., Tian, R. M., et al. (2012). Spatial and Species Variations in bacterial communities associated with corals from the Red Sea as revealed by pyrosequencing. Appl. Environ. Microbiol. 78, 7173–7184. doi: 10.1128/AEM.01111-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Levin, L. A. (2005). Ecology of cold seep sediments: interactions of fauna with flow, chemistry and microbes. Oceanogr. Mar. Biol. 43, 1–46. doi: 10.1201/9781420037449.ch1

CrossRef Full Text | Google Scholar

Lloyd, K. G., Schreiber, L., Petersen, D. G., Kjeldsen, K. U., Lever, M. A., Steen, A. D., et al. (2013). Predominant archaea in marine sediments degrade detrital proteins. Nature 496, 215–218. doi: 10.1038/nature12033

PubMed Abstract | CrossRef Full Text | Google Scholar

Logares, R., Brate, J., Bertilsson, S., Clasen, J. L., Shalchian-Tabrizi, K., and Rengefors, K. (2009). Infrequent marine-freshwater transitions in the microbial world. Trends Microbiol. 17, 414–422. doi: 10.1016/j.tim.2009.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozupone, C. A., and Knight, R. (2007). Global patterns in bacterial diversity. Proc. Natl. Acad. Sci. U.S.A. 104, 11436–11440. doi: 10.1073/pnas.0611525104

PubMed Abstract | CrossRef Full Text | Google Scholar

Maignien, L., Parkes, R. J., Cragg, B., Niemann, H., Knittel, K., Coulon, S., et al. (2013). Anaerobic oxidation of methane in hypersaline cold seep sediments. FEMS Microbiol. Ecol. 83, 214–231. doi: 10.1111/j.1574-6941.2012.01466.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyazaki, J., Higa, R., Toki, T., Ashi, J., Tsunogai, U., Nunoura, T., et al. (2009). Molecular characterization of potential nitrogen fixation by anaerobic methane-oxidizing archaea in the methane seep sediments at the number 8 Kumano Knoll in the Kumano Basin, offshore of Japan. Appl. Environ. Microbiol. 75, 7153–7162. doi: 10.1128/AEM.01184-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Mosier, A. C., Allen, E. E., Kim, M., Ferriera, S., and Francis, C. A. (2012). Genome sequence of “Candidatus Nitrosoarchaeum limnia” BG20, a low-salinity ammonia-oxidizing archaeon from the San Francisco Bay Estuary. J. Bacteriol. 194, 2119–2120. doi: 10.1128/JB.00007-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakagawa, T., Mori, K., Kato, C., Takahashi, R., and Tokuyama, T. (2007). Distribution of cold-adapted ammonia-oxidizing microorganisms in the deep-ocean of the northeastern Japan Sea. Microbes Environ. 22, 365–372. doi: 10.1264/jsme2.22.365

CrossRef Full Text | Google Scholar

Pernthaler, A., Dekas, A. E., Brown, C. T., Goffredi, S. K., Embaye, T., and Orphan, V. J. (2008). Diverse syntrophic partnerships from-deep-sea methane vents revealed by direct cell capture and metagenomics. Proc. Natl. Acad. Sci. U.S.A. 105, 7052–7057. doi: 10.1073/pnas.0711303105

PubMed Abstract | CrossRef Full Text | Google Scholar

Pester, M., Rattei, T., Flechl, S., Grongroft, A., Richter, A., Overmann, J., et al. (2012). amoA-based consensus phylogeny of ammonia-oxidizing archaea and deep sequencing of amoA genes from soils of four different geographic regions. Environ. Microbiol. 14, 525–539. doi: 10.1111/j.1462-2920.2011.02666.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Price, M. N., Dehal, P. S., and Arkin, A. P. (2009). FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol. Biol. Evol. 26, 1641–1650. doi: 10.1093/molbev/msp077

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864

PubMed Abstract | CrossRef Full Text | Google Scholar

Reeder, J., and Knight, R. (2010). Rapidly denoising pyrosequencing amplicon reads by exploiting rank-abundance distributions. Nat. Methods 7, 668–669. doi: 10.1038/nmeth0910-668b

PubMed Abstract | CrossRef Full Text | Google Scholar

Russ, L., Kartal, B., Op den Camp, H. J., Sollai, M., Le Bruchec, J., Caprais, J. C., et al. (2013). Presence and diversity of anammox bacteria in cold hydrocarbon-rich seeps and hydrothermal vent sediments of the Guaymas Basin. Front. Microbiol. 4:219. doi: 10.3389/fmicb.2013.00219

PubMed Abstract | CrossRef Full Text | Google Scholar

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Shao, S., Luan, X., Dang, H., Zhou, H., Zhao, Y., Liu, H., et al. (2014). Deep-sea methane seep sediments in the Okhotsk Sea sustain diverse and abundant anammox bacteria. FEMS Microbiol. Ecol. 87, 503–516. doi: 10.1111/1574-6941.12241

PubMed Abstract | CrossRef Full Text | Google Scholar

Stokke, R., Roalkvam, I., Lanzen, A., Haflidason, H., and Steen, I. H. (2012). Integrated metagenomic and metaproteomic analyses of an ANME-1-dominated community in marine cold seep sediments. Environ. Microbiol. 14, 1333–1346. doi: 10.1111/j.1462-2920.2012.02716.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Suess, E. (2010). “Marine cold seeps,” in Handbook of Hydrocarbon and Lipid Microbiology, ed. K. N. Timmis (Berlin: Springer-Verlag), 187–203.

Google Scholar

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197

PubMed Abstract | CrossRef Full Text | Google Scholar

Tolar, B. B., King, G. M., and Hollibaugh, J. T. (2013). An analysis of thaumarchaeota populations from the northern gulf of Mexico. Front. Microbiol. 4:72. doi: 10.3389/fmicb.2013.00072

PubMed Abstract | CrossRef Full Text | Google Scholar

Trichet, J., Défarge, C., Tribble, J., Tribble, G. W., and Sansone, F. J. (2001). Christmas Island lagoonal lakes, models for the deposition of carbonate-evaporite-organic laminated sediments. Sedim. Geol. 140, 177–189. doi: 10.1016/S0037-0738(00)00177-9

CrossRef Full Text | Google Scholar

Vigneron, A., Cruaud, P., Pignet, P., Caprais, J. C., Gayet, N., Cambon-Bonavita, M. A., et al. (2013). Bacterial communities and syntrophic associations involved in AOM process of the Sonora Margin cold seeps, Guaymas basin. Environ. Microbiol. 16, 2777–2790. doi: 10.1111/1462-2920.12324

PubMed Abstract | CrossRef Full Text

Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Tian, R. M., Gao, Z. M., Bougouffa, S., and Qian, P.-Y. (2014a). Optimal eukaryotic 18S and universal 16S/18S ribosomal RNA primers and their application in a study of symbiosis. PLoS ONE 9:e90053. doi: 10.1371/journal.pone.0090053

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhang, W., Cao, H., Shek, C. S., Tian, R., Wong, Y. H., et al. (2014b). Diversity and dispersal of eukaryotic microbes around a brine pool adjacent to the Thuwal cold seeps in the Red Sea. Front. Microbiol. 5:37. doi: 10.3389/fmicb.2014.00037

CrossRef Full Text

Zhang, W. P., Wang, Y., Tian, R. M., Bougouffa, S., Yang, B., Cao, H. L., et al. (2014). Species sorting during biofilm assembly by artificial substrates deployed in a cold seep system. Sci. Rep. 4, 6647. doi: 10.1038/srep06647

PubMed Abstract | CrossRef Full Text | Google Scholar