Journal of Oceanology and Limnology   2021, Vol. 39 issue(3): 846-864     PDF       
http://dx.doi.org/10.1007/s00343-020-0106-6
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

CHEN Ye, LI Siqi, XU Xiaoqing, MA Manman, MI Tiezhu, ZHEN Yu, YU Zhigang
Characterization of microbial communities in sediments of the South Yellow Sea
Journal of Oceanology and Limnology, 39(3): 846-864
http://dx.doi.org/10.1007/s00343-020-0106-6

Article History

Received Mar. 3, 2020
accepted in principle Apr. 6, 2020
accepted for publication May. 18, 2021
Characterization of microbial communities in sediments of the South Yellow Sea
Ye CHEN1,2,3, Siqi LI1,2,3, Xiaoqing XU3,4, Manman MA2,3,4, Tiezhu MI2,3,4, Yu ZHEN2,3,4, Zhigang YU2,5     
1 College of Marine Life Science, Ocean University of China, Qingdao 266003, China;
2 Laboratory for Marine Ecology and Environmental Science, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266237, China;
3 Key Laboratory of Marine Environment and Ecology, Ministry of Education, Qingdao 266100, China;
4 College of Environmental Science and Engineering, Ocean University of China, Qingdao 266100, China;
5 Key Laboratory of Marine Chemistry Theory and Technology, Ministry of Education/Institute for Advanced Ocean Study, Ocean University of China, Qingdao 266100, China
Abstract: Illumina sequencing and quantitative PCR (qPCR) based on the 16S ribosomal RNA (rRNA) gene were conducted to characterize the vertical distribution of bacterial and archaeal communities in the sediments of two sites from the South Yellow Sea. Both bacterial and archaeal communities showed a clear stratified distribution with sediment depth. The microbial communities in the upper layers were distinct from those in the deeper layers; the relative abundances of sequences of Thaumarchaeota, Gammaproteobacteria, and Actinobacteria were higher in the upper than in the deeper sediments, whereas the sequences of Bathyarchaeia, Lokiarchaeota, Euryarchaeota, Chloroflexi, and Deltaproteobacteria were relatively more abundant in the deeper sediments. Sediment depth and total organic carbon (TOC) can significantly influence both the bacterial and archaeal communities. Furthermore, bacterial and archaeal groups potentially involved in nitrogen, sulfur, and methane metabolism were detected in both sites. In our study, both ammonia-oxidizing bacteria (Nitrospira) and ammonia-oxidizing archaea (Candidatus Nitrosopumilus) were responsible for ammonia oxidization. Additionally, sulfur-reducing bacteria SEEP-SRB1 forming consortia with anaerobic methane-oxidizing archaea ANME-2a-2b were capable of anaerobic methane oxidation (AOM) in the 3400-02 sediment samples.
Keywords: microbial community    16S rRNA gene    high-throughput sequencing    South Yellow Sea    sediment    
1 INTRODUCTION

Continental margins are the main interface between terrestrial sediment sources and deep-sea depositional systems (Covault and Fildani, 2014), although they occupy less than 20% of the ocean surface area (Walsh, 1991). Sediment receives microbes and organic matter exporting from the upper water layer, and provides complex nutrients and solid surfaces for the growth of microorganisms (Wang et al., 2012). The biomass and taxon richness of the microbial community in sediments are much higher than those in the corresponding water bodies (Zinger et al., 2011). Marine microbial communities are major engines of biogeochemical processes, including carbon, nitrogen and sulfur, and they probably play pivotal roles in maintaining marine ecosystems to prevent environmental changes such as ocean acidification and global warming (Fuhrman, 2009; Graham et al., 2012). Analysis of the microbial community structure in sediments can enhance our understanding of benthic ecosystem processes and the roles that microbes play in overall oceanic processes.

The Southern Yellow Sea (SYS) is a semi-enclosed shallow shelf sea surrounded by the Chinese mainland and Korean Peninsula. The sediments of the SYS are mostly terrigenous matter, which is derived mainly from the Chinese mainland and secondarily from the Korean Peninsula (Shi et al., 2003). Simultaneously, the SYS has a complex dynamic system including monsoonal winds, waves, circulation (the Yellow Sea warm current, etc.), and tidal currents (Shi et al., 2003). This complex sedimentary environment in the SYS, including sediment sources and hydrodynamic conditions, can significantly influence the distribution of microbial communities and related biogeochemical processes. However, most studies have primarily focused on the water column and the surface sediments in the SYS (Dong et al., 2014; Liu et al., 2015b), and there is little information on this marine subsurface. Only Qiao et al. (2018) reported the vertical distribution of the bacterial community in the eastern China marginal seas (the SYS and East China Sea (ECS)); however, the potential functional role of archaea in this marine subsurface is still poorly understood. Once thought to inhabit only extreme environments, archaea are now known to inhabit diverse marine environments (DeLong, 1992; Fuhrman and Davis, 1997), where they constitute a relevant fraction of the microbial community (Karner et al., 2001). Archaea members, belonging to Euryarchaeota in marine systems, are recognized to perform the production and anaerobic oxidation of methane (Valentine, 2002; Liu and Whitman, 2008). The oxidation of ammonia to nitrite can be performed by the phylum Thaumarchaeota, which possesses the ammonia monooxygenase subunit A (amoA) genes (Brochier-Armanet et al., 2008; Spang et al., 2010; Berg et al., 2014). Therefore, understanding the variation in archaeal communities in sedimentary environments is crucial for predicting biogeochemical cycles in these environments.

This study was conducted to characterize the structure of the bacterial and archaeal communities at various depths collected in sediment samples from two sites of the SYS. To accomplish this goal, samples were evaluated using 16S rRNA gene high-throughput sequencing and quantitative PCR (qPCR). The present study aims to determine the abundance, diversity, and composition of the microbial community in sediments and their changes with depth, to identify functional microorganisms that can contribute to geochemical cycling in the sediments, and to reveal potential links between the abundance, community composition of bacteria and archaea and environmental factors.

2 MATERIAL AND METHOD 2.1 Sample collection

Sediment samples were collected from the nearshore site 3400-02 (length: 30 cm) along the Jiangsu Coastal Shoals and the offshore site 3500-07 (50 cm) in the center area of the SYS during a research cruise of R/V Kexue San Hao (Science #3 in Chinese) in September 2018 (Fig. 1). At each site, three parallel PVC tubes were used to subsample the collected sediments. The first core was prepared for pore water extraction. The pore water was collected by a Rhizon soil moisture sampler (Eijkelkamp Soil & Water, Giesbeek, Netherlands) and used to measure the pH, dissolved inorganic nitrogen (NH4+, NO3-, and NO2-), ions (SO42- and Cl) and elements (Fe and Mn). The second PVC core was used for the collection for gas measurement. Sediment samples were collected using sterile 3-mL cut-end syringe barrels through predrilled holes in the core tube, and sealed in glass vials containing 6 mL of 10% NaOH for methane analysis. The last PVC core was sliced at 2-cm intervals for organic matter measurement and DNA extraction.

Fig.1 Sampling locations in the South Yellow Sea
2.2 Determination of environmental parameters

The environmental parameters (temperature, salinity, and dissolved oxygen concentration) of the bottom water samples were recorded with a CTD (SBE 11 plus, Sea-Bird Electronics, Inc.) in situ. Other parameters were measured in the laboratory. pH was determined using an S479-uMix (METTLER TOLEDO SevenExcellence, Switzerland). The concentrations of inorganic nitrogen compounds (NH4+, NO3-, and NO2-) were determined by a QuAAtro nutrient autoanalyzer (Seal Analytical Ltd., United Kingdom), and the concentrations of SO42- and Cl were analyzed using ICS-3000 ion chromatography (Dionex, USA). The elemental concentrations of Fe and Mn were measured by Thermo Scientific iCAP RQ ICP-MS. The content of total organic carbon (TOC) and total nitrogen (TN) in the sediments were determined by elemental analyzer (EURO EA3000, Italy). Methane concentrations were determined using a gas chromatograph (GC; Agilent 7890B) equipped with a Porapark-Q column and a flame ionization detector (FID) by injecting 5 mL of headspace gas.

2.3 DNA extraction and quantitative PCR

Genomic DNA was extracted using the PowerSoil DNA Isolation Kit (Mo Bio, USA) in accordance with the manufacturer's protocol. The extracted DNA was then stored at -80 ℃ for further analyses. A SYBR Green method was used to quantify the total abundance of bacterial and archaeal communities, and all qPCR assays were performed on an ABI PRISM® 7500 Sequence Detection System (Applied Biosystems, USA). The bacterial 16S rRNA and archaeal 16S rRNA gene copies were quantified using the universal primer pair 338F/806R (Peiffer et al., 2013; Liu et al., 2015a) and U519F/806R (Porat et al., 2010; Shehab et al., 2013), respectively. Each 20-μL mixture contained 10-μL FastStart Universal SYBR Green Master (Rox) (Roche Diagnostics, Germany), 0.2-μL bovine serum albumin (BSA), 0.6 μL of each primer and 2.0-μL sediment DNA. Both bacterial and archaeal qPCR amplification conditions were as follows: an initial denaturation step at 95 ℃ for 10 min and 40 cycles of 95 ℃ for 15 s and 58 ℃ for 2 min. All qPCRs were performed in triplicate for each sample. Amplification products were verified by melting curve analysis.

Standard curves for qPCR assays were obtained with serial 10-fold dilutions of plasmids carrying the target gene fragments. Amplification efficiencies were 90% and 99% for bacteria and archaea, respectively, and a strong linear inverse relation (R2 > 0.99) was considered credible.

2.4 High-throughput sequencing and data analysis

Five sediment layers (2–4, 8–10, 14–16, 20–22, and 26–28 cm) from site 3400-02 and eight sediment layers (2–4, 8–10, 14–16, 20–22, 26–28, 32–34, 38– 40, and 44–46 cm) from site 3500-07 were chosen for bacterial and archaeal 16S rRNA gene tag sequencing. PCR amplifications of the bacterial 16S rRNA gene targeting the V3–V4 region and the archaeal 16S rRNA gene targeting the V4–V5 regions were conducted using primers 338F (5′-ACTCCTACGGGAGGCAGCA) and 806R (5′-GGACTACHVGGGTWTCTAAT) (Lee et al., 2012) and Arch519F (5′-CAGCCGCCGCGGTAA) and Arch915R (5′-GTGCTCCCCCGCCAATTCCT) (Coolen et al., 2004), respectively. The PCR conditions for the 338F/806R primers were as follows: 98 ℃ for 5 min; 24 cycles of 98 ℃ for 30 s, 52 ℃ for 30 s and 72 ℃ for 45 s; and then 72 ℃ for 5 min. The PCR conditions for the Arch519F/ Arch915R primers consisted of 98 ℃ for 5 min; 35 cycles of 98 ℃ for 30 s, 65 ℃ for 30 s and 72 ℃ for 45 s; followed by 72 ℃ for 5 min. PCR was carried out on an ABI 2720 (Applied Biosystems, USA), and PCR products were electrophoresed on a 1.2% polyacrylamide gel. The target fragment was recovered with Agencourt AMPure Beads (Beckman Coulter, Indianapolis, IN) and quantified using the PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). High-throughput sequencing of the 16S rRNA genes were conducted on the Illumina MiSeq platform (Personal Biotechnology Co., Ltd., Shanghai, China).

Sequenced reads were de-multiplexed, qualityfiltered, and processed using the QIIME software package (http://qiime.org/) according to the following three criteria (Caporaso et al., 2012). The paired-end reads were joined with at least a 10-bp overlap and no mismatches using FLASH (Magoč and Salzberg, 2011). Sequences were clustered into operational taxonomic units (OTUs) based on their sequence similarity at a 97% dissimilarity level using UCLUST. A representative sequence for each OTU was assigned against the SILVA SSU Ref NR database (https://www.arb-silva.de/no_cache/download/archive/release_132ARB_files/) via QIIME.

2.5 Statistical analyses

Any OTU containing less than 0.001% of the total sequences in samples was discarded to reduce the effect of spurious OTUs on further analysis (Bokulich et al., 2013). To fairly compare all the samples at the same sequencing depth, sequences were randomly reduced based on the smallest read numbers in each sample. Then, the diversity indices, including Good's coverage, Shannon, Simpson, and Chao 1, were calculated using the QIIME software package. For beta diversity, the Bray-Curtis dissimilarity-based dendrogram was plotted to examine the differences in bacterial and archaeal communities between the sediment samples. Comparisons of parameters between the two sites were performed using SPSS by one-way ANOVA. The Spearman correlation test was used to evaluate correlations between microbial abundance, diversity, and the 20 most abundant bacterial and archaeal genera with environmental parameters. Then, heatmaps were constructed to visualize Spearman's correlation using an online tool (http://www.lc-bio.cn).

3 RESULT 3.1 Environmental parameter characterization

The environmental parameters of sediment samples from the South Yellow Sea were different between the two sites (Supplementary Figs.S1 & S2). The pH was 7.80–8.21 and 7.56–8.05 at site 3400-02 and 3500-07, respectively. The pH in the upper 22 cm of the site 3400-02 was higher than at site 3500-07 and then decreased in the deeper sediments. The concentrations of NH4+ (110.74–1 320.90 μmol/L, average 471.46 μmol/L) in the pore water at site 3400-02 were higher than at site 3500-07 (24.53–68.35 μmol/L, average 57.91 μmol/L) (one-way ANOVA, P < 0.05). Higher concentrations of NO3- (0.57–15.14 μmol/L, average 3.39 μmol/L) were also detected at site 3400-02 than at site 3500-07 (0.52–5.49 μmol/L, average: 1.15 μmol/L) (one-way ANOVA, P < 0.05). The concentrations of SO42- (average 20.34 mmol/L) at site 3400-02 were lower than at site 3500-07 (average 26.07 mmol/L) (ANOVA, P < 0.05). Lower contents of TOC and TN were also observed at site 3400-02. The CH4 concentration at site 3400-02 and 3500-07 ranged from 1.30 to 2.90 μmol/kg (average 1.96 μmol/ kg) and from 1.22 to 2.84 μmol/kg (2.09 μmol/kg), respectively (Supplementary Fig.S3). There were also distinct differences for some environmental parameters in the vertical distribution of the two sites (Supplementary Figs.S1–S2). The NO3- and NO2- concentration in pore water and TOC content in sediment had a tendency to decrease with depth, while NH4+ had an opposite tendency.

3.2 Abundance of total bacterial and archaeal 16S rRNA genes

Quantitative PCR was used to detect the abundance of the bacterial and archaeal communities based on the 16S rRNA genes (Fig. 2). The observed bacterial 16S rRNA gene abundance at 3400-02 and 3500-07 ranged from 8.50×107 to 3.71×108 copies/g and 1.25×108 to 1.03×109 copies/g (wet weight), respectively (Fig. 2a). The bacterial abundance showed a decreasing trend with sediment depth at both sites. The relative abundance of bacteria in the total microbial community ranged from 97.12% to 99.41% at site 3400-02 and 93.00% to 99.78% at site 3500-07, and the relative bacterial abundance at site 3500-07 decreased with sediment depth (Fig. 2b).

Fig.2 Vertical profiles of total bacterial and archaeal 16S rRNA gene copy numbers in sediment cores at the two sites (a) and the proportion of bacterial abundance or archaeal abundance relative to the total microbial abundance in sediment cores at the two sites (b)

The observed archaeal 16S rRNA gene abundance ranged from 5.38×105 to 7.93×106 copies/g and 1.55×106 to 2.92×107 copies/g at 3400-02 and 3500-07 (wet weight), respectively (Fig. 2a). In contrast to the vertical variation of bacteria abundance along the sediment depth, the archaea abundance decreased with the sediment depth in the upper 10 cm at site 3400-02 and then showed marked fluctuations. At site 3500-07 in contrast, the archaea abundance increased with sediment depth. For the two sites, the relative abundance of archaea within the total microbial community was 0.59%–2.88% at site 3400-02 and 0.22%–7.00% at site 3500-07, and the relative archaeal abundance at site 3500-07 increased with sediment depth (Fig. 2b).

The bacterial 16S rRNA gene abundance was significantly positively correlated with TOC, TN, SO42-, and Cl (P < 0.05) and negatively correlated with NH4+ and NO3- (P < 0.05), while the archaeal 16S rRNA gene abundance was positively correlated with sediment depth, TOC, TN and SO42- (P < 0.05) and negatively correlated with NH4+, NO3-, and C/N (P < 0.05) (Fig. 3).

Fig.3 Heatmap of Spearman's correlation coefficients between the abundance and diversity of bacterial and archaeal communities and environmental factors Correlation coefficient with P-value score less than 0.05 and 0.01 is labeled with "*" and "**", respectively.
3.3 Bacterial and archaeal diversity and richness

A total of 419 741 high-quality bacterial 16S sequences and 510 751 archaeal 16S sequences were recovered from Illumina-based analysis (Tables 1 & 2). The number of OTUs varied from 1 809 to 2 696 per sample and 804 to 3 005 per sample at a 97% sequence similarity level for the bacterial and archaeal communities, respectively. Good's coverage of bacterial and archaeal communities at site 3400-02 and 3500-07 ranged from 93.21%–100.00% and 95.70%–100.00%, respectively, which indicated that the sequencing depth was sufficient to cover most of bacterial and archaeal communities. At the corresponding sediment layer of the two sites, the bacterial and archaeal diversity (Shannon-Wiener and Simpson's diversity index) was higher at site 3500-07 than 3400-02, with the exception of the bacterial community in the 14–16-cm layer. At site 3400-02, the bacterial Shannon-Wiener index values were higher than the archaeal community in the corresponding sediment layer, while at site 3500-07, the bacterial Shannon-Wiener index values were higher than the archaeal community in the 2–22-cm layers but became lower in the deeper sediments (26– 46 cm).

Table 1 Vertical distribution of sediment bacterial richness and diversity estimates based on the 16S rRNA gene
Table 2 Vertical distribution of sediment archaeal richness and diversity estimates based on the 16S rRNA gene

With respect to bacterial diversity, Simpson's diversity index was negatively correlated with sediment depth (P < 0.05). With respect to archaeal diversity, both Simpson's diversity index and Shannon-Wiener index were positively correlated with sediment depth, TN, and SO42- (P < 0.05) but negatively correlated with NO3- (P < 0.05), and the archaeal Chao 1 index was positively correlated with sediment depth and TN (P < 0.05) but negatively correlated with NO3- (P < 0.05) (Fig. 3).

3.4 Bacterial 16S rDNA analysis

The composition of bacterial communities of sediment samples from the SYS is shown in Fig. 4. Most of the bacterial sequence reads were associated with Proteobacteria, accounting for 43.41% of all the sequences, followed by Chloroflexi (16.35%), Acidobacteria (9.31%) and Actinobacteria (8.18%). Among the Proteobacteria, Deltaproteobacteria (19.33%) and Gammaproteobacteria (16.94%) were the most abundant, followed by Alphaproteobacteria (7.14%). Chloroflexi mainly consisted of two classes: Anaerolineae (10.84%) and Dehalococcoidia (4.10%).

Fig.4 Relative abundance of the bacterial community (including classes of Proteobacteria and Chloroflexi) observed in sediment cores collected from the SYS at the phylum level with an abundance > 1%

At site 3400-02, Gammaproteobacteria represented the predominant group in each sediment layer, followed by Deltaproteobacteria and Anaerolineae. At site 3500-07, sequences from Gammaproteobacteria dominated in the 2–10-cm sections and decreased with sediment depth. The relative abundances of Deltaproteobacteria increased with sediment depth and dominated in the 14–46-cm sections. Additionally, Anaerolineae and Dehalococcoidia belonging to the Chloroflexi were enriched in the deeper sediments (26–46 cm) of site 3500-07.

3.5 Archaeal 16S rDNA analysis

Nine different archaeal phyla were identified in the sediment samples (Fig. 5a). Crenarchaeota (44.12%) and Thaumarchaeota (32.58%) were the most abundant phyla of archaeal sequences, followed by Asgardaeota (11.97%) and Euryarchaeota (5.67%). Within the Crenarchaeota, Bathyarchaeia (43.81%) (Formerly named MCG) had the highest relative abundance (Fig. 5b). Sequences belonging to Bathyarchaeia had lower abundances in the surface sediments of sites 3400-02 (18.93%) and 3500-07 (3.64%); however, they predominated in the subsurface sediment layers of the two sites, ranging from 36.58% to 60.83% (3400-02) and 36.81% to 57.28% (3500-07). Nitrososphaeria was the most abundant class within the Thaumarchaeota, which was found to predominate in the surface sediment of the two sites and to decrease with sediment depth at site 3500-07. Lokiarchaeia (11.48%) was the most dominant class within Asgardaeota, and it was mainly present in the subsurface sediments (8–46 cm) of site 3500-07. Among Euryarchaeota, the majority of sequences were assigned to Thermoplasmata, which was particularly enriched in deeper sediment layers (26–46 cm) at site 3500-07 (Fig. 5b).

Fig.5 Relative abundance of archaeal groups observed in sediment cores collected from the SYS a. archaeal community at the phylum level; b. archaeal community at the class level with an abundance > 1%.
3.6 Community comparison and relationships between community structure and environmental factors

The bacterial and archaeal community structures clearly differed between the two sites and among different sediment depths, as displayed by the clustering analysis (Fig. 6). Overall, the bacterial and archaeal community structures in the sediment layers from site 3400-02 formed distinct clusters, and those in the sediment layers from site 3500-07 tended to cluster together. The upper layer samples (2–22 cm) of site 3500-07 clustered, whereas the deep layer samples (26–46 cm) tended to cluster together for the bacterial community. For the archaeal community, the surface sediment layer of site 3500-07 was obviously different from the other layers, and the other layers clustered into two groups: the upper layer samples (8–22 cm) clustered into one group, whereas the deep layer samples (26–46 cm) clustered together.

Fig.6 UPGMA hierarchical clustering based on a Bray-Curtis distance matrix for the bacterial composition (a) and archaeal composition (b) at the OTU level

Biota-environment (BIOENV) identified that sediment depth, TOC, C/N, NH4+, and Fe were most strongly correlated with the bacterial community, while sediment depth, TOC, TN, C/N, and Fe were most strongly correlated with the archaeal community. These variables explained 52.87% and 53.14% of the variation in the bacterial and archaeal communities, respectively, as determined by CCA (Fig. 7). The composition of the bacterial community was significantly correlated with sediment depth, TOC and NH4+ and the archaeal community was significantly correlated with sediment depth, TOC and TN via Monte Carlo permutation significance tests.

Fig.7 Canonical correspondence analysis (CCA) biplot of environmental parameters and microbial community structure at the OTU level a. bacterial community; b. archaeal community. 3500-07: B1 (2-4 cm); B2 (8-10 cm); B3 (14-16 cm); B4 (20-22 cm); B5 (26-28 cm); B6 (32-34 cm); B7 (38-40 cm); B8 (44-46 cm); 3400-02: A1 (2-4 cm); A2 (8-10 cm); A3 (14-16 cm); A4 (20-22 cm); A5 (26-28 cm). Correlation coefficient with P-value score less than 0.05 and 0.01 is labeled with "*" and "**", respectively.
3.7 Microbes involved in biogeochemical cycles

The 16S rRNA genes showed that ammonia-oxidizing bacteria Nitrosospira and ammonia-oxidizing archaea Candidatus Nitrosopumilus were the dominant groups associated with ammonia oxidation at the two sites. The relative abundances of Nitrosospira (ammonia-oxidizing bacteria (AOB)) and Candidatus Nitrosopumilus (ammonia-oxidizing archaea (AOA)) were higher in the upper sediment layers at both sites, and their proportions decreased with increasing depth (Fig. 8a & b). Furthermore, the proportion of AOA in the archaeal community was much higher than AOB in the bacterial community at each site. The bacterial genus Nitrospira was the main group associated with nitrate oxidation and exhibited relative abundances of 0.14%–0.29% and 0.02%– 1.17% at 3400-02 and 3500-07, respectively (Fig. 8c). With respect to denitrifying bacteria, Bacillus was the main group in most sediment layers, except the 2–4 cm and 20–22 cm at site 3400-02, in which Vibrio or Microbacterium was the predominant group (Fig. 8d).

Fig.8 Relative abundances of microbes involved in the nitrogen process at the two sites a. ammonia-oxidizing bacteria (AOB); b. ammonia-oxidizing archaea (AOA); c. nitrate-oxidizing bacteria (NOB); d. denitrifying bacteria.

Thirty-four genera (Desulfatiglans, SEEP-SRB1, Desulfuromonas, etc.) of sulfate-reducing bacteria (SRBs) affiliated with Delta-proteobacteria and one genus (Thermodesulfovibrio) belonging to Nitrospirae were detected to participate in sulfate reduction. Desulfatiglans and SEEP-SRB1 dominated at the two sites, and their relative abundance increased with sediment depth (Fig. 9a). Sulfur-oxidizing bacteria (SOB) were affiliated with Alpha-, Delta-proteobacteria, Epsilon-proteobacteria, and Spirochaetae (Fig. 9b). The genus Spirochaeta within Spirochaetae represented the dominant SOB at site 3400-02 and deeper sediment layers (26–46 cm) at site 3500-07 and tended to increase with sediment depth at each site. The SOB bacteria AqS1 were dominant in the upper sediment layers (2–22 cm) of site 3500-07.

Fig.9 Relative abundances of microbes involved in the sulfur process a. sulfur-reducing bacteria (SRB); b. sulfur-oxidizing bacteria (SOB).

The methanogens had a higher relative abundance at site 3400-02, accounting for 0.15%–0.66% (0.35% on average), than at site 3500-07, with a relative abundance of 0.01%–0.17% (on average 0.05%) (Fig. 10a). The composition of methanogens varied significantly between the two sites. Fifteen genera (Methanosarcina, Methanolobus, Methanobacterium, etc.) were discovered at site 3400-02, whereas seven genera (Methanococcoides, Methanosarcina, Methanobrevibacter, etc.) were discovered at site 3500-07. Among the methanogens, Methanosarcina, with a relative abundance of 0.13% (on average), was the dominant methanogen at site 3400-02, whereas Methanococcoides, with a relative abundance of 0.04% (on average), was the dominant methanogen at site 3500-07, except at 14–16 cm, in which Methanosarcina was the dominant group. Sequence analysis of the 16S rRNA genes also detected a minor fraction of microbial communities involved in methane aerobic oxidation. Methylobacterium was the predominant methanotroph at site 3500-07, while the methanotrophs at site 3400-02 consisted of three genera: Methylobacterium, Methylomarinum, and Methyloprofundus (Fig. 10b). Anaerobic methane-oxidizing archaeas (ANMEs), including ANME-1b, ANME-2a-2b, and ANME-3, associated with CH4 anaerobic oxidation, were detected in the sediment samples of site 3400-02 (Fig. 10c). ANME-1b dominated in the 2–4-cm and 8–10-cm layers, while ANME-2a-2b dominated in the other sediment layers. At site 3500-07, only a minor fraction of ANME-3 was detected in the 32–34-cm layer.

Fig.10 Relative abundances of microbes involved in the methane process at the two sites a. methanogens; b. aerobic methanotrophic bacteria (MOB); c. anaerobic methanotrophic archaea (ANME).
4 DISCUSSION 4.1 Bacterial and archaeal abundance in sediments of the SYS

In this study, the composition and distribution of bacterial and archaeal communities were characterized in the sediments of two sites of the SYS. The abundance of the bacterial community was significantly higher than the abundance of the archaeal community at each site (one-way ANOVA, P < 0.01) (Fig. 2a). This result is in agreement with many previous studies, which demonstrated that bacteria rather than archaea dominate in sediments from various marine environments (Schippers and Neretin, 2006; Schippers et al., 2012; Liu et al., 2015b). Organic matter availability in marine sediments is recognized as an important factor that affects the abundance of the microbial community (Säwström et al., 2016; Qiao et al., 2018). In this study, the TOC and TN contents at site 3500-07 were 2.35–7.73 times and 1.05–14.75 times higher than at site 3400-02 in the corresponding layer. Accordingly, the bacterial and archaeal 16S rRNA gene abundances at site 3500-07 were 1.18–6.80 times and 0.10–12.65 times higher than at site 3400-02 in the corresponding layer, with the exception of the archaeal community in the upper 8-cm layers. It was indicated that fresh and marine organic matter (higher TOC and TN) can support higher microbial abundance.

4.2 Bacterial and archaeal community compositions in sediments of the SYS

The bacterial community composition at each site was dominated by Proteobacteria, which mainly consisted of the classes Alpha-, Delta- and Gammaproteobacteria (Fig. 4). The phylum Proteobacteria, as the most dominant group, was also observed in shallow sea and deep-sea sediments (Zhu et al., 2013; Liu et al., 2015b). Most members of Delta-proteobacteria are well known as SRBs (Varon-Lopez et al., 2014), while groups within the Gamma- and Alpha-proteobacteria have been shown to perform sulfur oxidation (Bagchi et al., 2005; Lenk et al., 2011). Previous studies have demonstrated their dominance in the sediments of the northern Chinese marginal seas, implying that sulfate and/or sulfur-related processes could be the most important biogeochemical reactions in these areas (Zhu et al., 2013; Liu et al., 2015b; Zhang et al., 2017; Qiao et al., 2018). Consistent with previous studies of marine sediments (Lentini et al., 2014; Mahmoudi et al., 2015; Qiao et al., 2018), Chloroflexi dominated in deep sediments (3500-07: 26–46 cm). Compared with other microorganisms, Chloroflexi can better utilize recalcitrant organic matter buried in deep sediments, which likely explains their higher abundance in such niches (Wilms et al., 2006; Zhou et al., 2017). In this study, the major component of Chloroflexi was Anaerolineae, which contains members of gram-negative chemoorganotrophic bacteria and can only grow under strict anaerobic conditions as previously described (Yamada et al., 2006). Interestingly, the phylum Chlamydiae accounted for 0.11%–0.22% and 0.08%–1.23% of the total bacterial sequences at site 3400-02 and 3500-07, respectively (Fig. 4). Traditionally, Chlamydiae are regarded as important pathogens and symbionts with eukaryotic hosts. However, a recent report showed that in Arctic Mid-Ocean, Chlamydiae dominated in deep anoxic marine, with relative abundances up to 43% of the bacterial community (Dharamshi et al., 2020). The features of host dependency in marine-sediment Chlamydiae suggested that these members play an important and overlooked ecological role in marine ecology and may have alternate lifestyle strategies (Dharamshi et al., 2020).

In our study, the archaeal community showed a clear stratified distribution as the aerobic group. Thaumarchaeota was abundant in the surface layer and decreased rapidly with depth, whereas Bathyarchaeia and Lokiarchaeia tended to be dominant in the lower layers and increased accordingly along the depth profiles (Fig. 5b). Sequences affiliated with Thaumarcheaota accounted for a high percentage of the total archaeal sequences in the surface sediments of the two sites. Thaumarchaeota is the most abundant archaea on earth and widely distributed in freshwater, soils, marine waters, and sediments (Offre et al., 2013). Species within the Nitrososphaeria class of Thaumarchaeota have been identified as aerobic ammonia-oxidizers (Stahl and de la Torre, 2012). The high relative abundance of Nitrososphaeria observed herein suggested that this group could play a major role in nitrification and nitrogen cycling in sediments. Bathyarchaeia (formerly named MCG) and Lokiarchaeia (formerly named DSAG or MBGB) have been proposed as major groups in heterotrophic metabolism of organic matter in anoxic sediments (Inagaki et al., 2006; Nunoura et al., 2012). Moreover, previous studies have shown that members of the Bathyarchaeia are anaerobic methanotrophs (Evans et al., 2015), methanogens (Evans et al., 2015), acetogens (He et al., 2016b; Lazar et al., 2016) or anaerobic heterotrophs (Lloyd et al., 2013, Lazar et al., 2016), indicating that Bathyarchaeia can play an important role in the organic carbon cycle. Lokiarchaeia or MBGB has been suggested to play a role in biogeochemical processes such as sulfate reduction and methane oxidation (Inagaki et al., 2006). The majority of sequences of Euryarchaeota are affiliated with Thermoplasmatales, which are regularly detected in marine sediments (Durbin and Teske, 2012). Microorganisms from the archaeal class Thermoplasmata have been found to be linked to methanogenic activities (Poulsen et al., 2013). Taken together, our results indicate that the complex archaeal communities in our given sediments might be involved in diverse metabolic processes, especially carbon and nitrogen cycling.

4.3 Archaea plays important roles in deep-water sediments

In our study, the bacterial community abundance and their ratio within the total microbial community showed a slightly decreasing tendency from the surface layers down to the subsurface layers of each site (Fig. 2). With respect to the archaeal community abundance, a different distribution pattern was identified at site 3500-07. The abundance of archaea and their ratio within the total microbial community increased with sediment depth (Fig. 2). The changing pattern for bacterial and archaeal abundance observed herein was consistent with many previous reports in the shallow sediments of marine benthic habitats (Urakawa et al., 2000; Jiang et al., 2011; Zhang et al., 2017; Pala et al., 2018; Qiao et al., 2018). Heterotrophic bacteria respiration consumes oxygen rapidly in upper layers of the sediment (Qiao et al., 2018). As a result, the microbial carbon oxidation rates declined in deeper sediments followed by a subsequent decrease in bacterial abundance (Qiao et al., 2018). However, archaea can unite by a universal ecological ability to adapt to energy stress in the subsurface environment where food and energy resources are poor (Gold, 1992). Moreover, the archaeal community had a higher richness and diversity in the deeper sediments (3500-07: 20–46 cm). Archaea Bathyarchaeia, Lokiarchaeia, and Thermoplasmata were enriched in these deeper water sediments, leading to a higher archaeal diversity. These members belong to anaerobic heterotrophs that can consume buried carbon (Biddle et al., 2006; Lin et al., 2015), indicating that they could contribute most to the anaerobic degradation of organic matter in such niches.

4.4 Microbes with potential ecological functions

Although accurate measures of metabolic functions are impossible based on 16S rRNA gene sequencing data alone (Tringe and Hugenholtz, 2008; Sun et al., 2018), we can use 16s rRNA taxonomies to predict the functional potential of microbes. Ecological function can be inferred from a genus level phylogeny, which provides better phylogenetic resolution than higher taxonomic categories (Ho et al., 2016; Sun et al., 2018). Nitrification is a two-step biological process that converts ammonia to nitrite and subsequently to nitrate. Ammonia oxidation is generally regarded as the rate-limiting step in nitrification and is carried out by both AOB and AOA. In this study, 16S rRNA gene sequencing analysis of microbial communities indicated that the majority of AOB fell within the Nitrosospira and that the majority of AOA fell within the Candidatus Nitrosopumilus at both sites (Fig. 8a & b). This result is consistent with findings showing that AOB Nitrosospira and AOA Nitrosopumilus usually dominate the ammoniaoxidizing community in marine systems (He et al., 2016a, 2018; Yu et al., 2016).

Sulfate reduction, mediated by SRB groups, is the dominant organic matter mineralization process in sulfate-rich shelf sediments. In this study, the genera Desulfatiglans and SEEP-SRB1 dominated at both sites (Fig. 9a). SEEP-SRB1 forming consortia with methane-oxidizing archaea ANME-2a-2b are capable of anaerobic methane oxidation (AOM) by using sulfate as the terminal electron acceptor (Timmers et al., 2016). ANME-2a-2b was also detected in the samples of site 3400-02 except for the surface layer (0–2 cm), suggesting that sulfate-dependent anaerobic methane oxidation (SAMO) might occur at this site. Robador et al. (2016) found that Desulfatiglans was one of the major groups of sulfate-reducing bacteria in polar, temperate, and tropical marine sediments. Groups within Desulfatiglans can completely oxidize organic matter such as phenol, benzoate, or aromatic compounds to carbon dioxide (Suzuki et al., 2014). In our study, AqS1 (of the Gammaproteobacteria class) and Spirochaeta 2 (of the Spirochaetae phylum) were the dominant sulfur-oxidizing bacteria. Previous studies have indicated that AqS1 can also take part in carbon monoxide oxidation and inorganic phosphate assimilation (Gauthier et al., 2016). Interestingly, the trends of the relative abundance of Spirochaeta 2 at the two sites were consistent with those of the relative abundance of SRB. This finding can be explained by the possible coexistence of Spirochaeta phylotypes with SRB by supplying them with carbon sources and electron donors (acetate, H2, and CO2) (Berlanga et al., 2008).

Methanogenesis is the final degradation process of organic matter in anoxic environments. SRB usually outcompete methanogens for substrates (H2, CO2, and acetate) due to their higher energy yield and higher substrate affinity in sulfate-rich environments (Schönheit et al., 1982; Lovley and Klug, 1983). In our study, a low abundance of methanogen groups was detected due to a high level of sulfate. Methanosarcina was predominant at site 3400-02 (Fig. 10a), which is made up of highly versatile species that are capable of methanogenesis using various substrates (acetate, C1 compounds, CO2/H2) (Liu and Whitman, 2008). The genus Methanococcoides predominated at site 3500-07 (Fig. 10a) and can only use C1 compounds to produce methane (Liu and Whitman, 2008). The dominance of methylotrophic methanogens in our study implies that methylotrophic methanogenesis might be the main methanogenic pathway at site 3500-07. Studies have shown that methylotrophic methanogenesis dominated CH4 production in shallow marine sediments with higher levels of sulfate in Aarhus Bay, Orca Basin, the Gulf of Mexico, and the western Mediterranean Sea (Xiao et al., 2018; Zhuang et al., 2018), because methylotrophic methanogens can be competitive when they compete with sulfate reducers for methylated substrates (e.g., methanol, methylamines or other C1 compounds). However, according to the correlation analysis, there was no significant correlation between the relative abundance of methanogens and the concentration of CH4 in our studied sediments (R=-0.223, P > 0.05), probably because the CH4 concentration was mediated by both methanogen and methanotrophs. Both aerobic methanotrophic bacteria and ANME were detected at the two sites, suggesting that aerobic and anaerobic methane oxidation are possible in such a niche. ANME often forms consortia with SRB to catalyze sulfate-dependent aerobic methane oxidation (S-DAMO). ANME-1 and ANME-2 are always associated with the Desulfosarcina/Desulfococcus cluster, while ANME-3 is associated with the Desulfobulbus branch (Boetius et al., 2000; Orphan et al., 2002; Niemann et al., 2006). No Desulfobulbus sequences or minor Desulfococcus and Desulfosarcina sequences were detected in our samples. However, ANME-2a-2b and SEEP-SRB1 coexisted in the 8–28-cm layer at site 3400-02 (Figs. 9a & 10c), further confirming that AOM can occur at site 3400-02 through the combination of these two groups.

4.5 Potential influence of environmental parameters on the microbial community

Marine microbial communities are affected extensively by environmental changes. Many previous studies have focused on the environmental parameters driving the distribution patterns of the microbial community on the surface of the northern Chinese marginal seas (Dong et al., 2014; Xiong et al., 2014; Liu et al., 2015b). Temperature and dissolved oxygen in bottom water have been detected to be vital factors in determining both bacterial and archaeal communities, whereas chlorophyll a in sediment was significant only in structuring archaeal communities in the Chinese marginal seas (Bohai Sea, Yellow Sea, and the northern East China Sea) (Liu et al., 2015b). Shifts in the bacterial community were observed in nitrogen-polluted sediments in the East China Sea (Xiong et al., 2014), indicating that the composition of nutrients is an important factor controlling the microbial community. Few studies have focused on the vertical profile of microbial communities in the northern Chinese marginal seas, with the exception of Qiao et al. (2018)'s study.

In this study, both bacterial and archaeal communities were found to be significantly influenced by sediment depth and TOC (Fig. 7). Sediment depth has been suggested a primary factor determining the difference in bacterial and archaeal communities in various marine ecosystems (Reed et al., 2006; Devereux et al., 2015; Qiao et al., 2018). We also confirmed in our study that a majority of the identified taxa abundance at the genus level increased or decreased along the sediment depth profile with significant statistical support (Supplementary Fig. S4). Sediment depth may function as a comprehensive proxy of multiple environmental parameters rather than directly affecting microbial diversity. The changes in the availability of oxygen or redox state with sediment depth could be the major factor contributing to the stratified distribution of the microbial community. The most abundant archaeal genus Candidatus Nitrosopumilus was abundant in the surface layer and rapidly decreased along the depth profile. This group of ammonia oxidizers has a metabolic trait for aerobic respiration. Most Deltaproteobacteria, including Desulfatiglans and SEEP-SRB1, belong to strictly anaerobic genera and are sulfate- and sulfur-reducing bacteria, and their relative abundance was very positively correlated with sediment depth. Additionally, microbes preferentially degrade the easily utilizable organic matter in surface sediments, whereas recalcitrant organic matter accumulates in deeper sediments (Cowie and Hedges, 1994; Wakeham et al., 1997). Consequently, microorganisms inhabiting deeper sediments have to meet their metabolic demands by relying on more recalcitrant organic matter (Zonneveld et al., 2010). In this study, the dominant bacterial phylum Chloroflexi and dominant archaeal classes Bathyarchaeia and Lokiarchaeia were believed to have an anaerobic and heterotrophic lifestyle depending on the buried recalcitrant organic matter, probably explaining their higher abundance in the deep sediments (Biddle et al., 2006; Wilms et al., 2006).

The organic matter content is another factor that shapes the bacterial and archaeal communities in marine sediment (Liu et al., 2015b; Oni et al., 2015). At the genus level of the bacterial community, OM60 (NOR5) andAqS1affiliated with Gammaproteobacteria and Lutimonas and Muriicola affiliated with Flavobacteriales preferred substrates with a higher content of TOC (Supplementary Fig.S3). Species within OM60 (NOR5) are aerobic anoxygenic phototrophic bacteria that can use both light and organic compounds as energy sources (Spring et al., 2013). Genomic evidence suggests that AqS1 has a versatile heterotrophic metabolism depending on the organic compounds (Gauthier et al., 2016). The class Flavobacteria is known to play an important role in the recycling of organic carbon in the marine environment (Lozada et al., 2014). They are also known to metabolize high-molecular-weight polymeric substances (Fernández-Gómez et al., 2013). With respect to the archaeal community, we identified that some methanogen genera, including Methanobacterium, Methanolobus, Methanosarcina, and Methanogenium, were highly negatively correlated with TOC and TN (Supplementary Fig.S3). This phenomenon can be explained by the observation that methanogens cannot use macromolecular organic matter and fail to compete with other microorganisms in sediments with a high organic carbon content.

5 CONCLUSION

Bacterial and archaeal diversity and abundance at site 3500-07 were higher than at 3400-02. Both bacterial and archaeal communities exhibited a clearly stratified distribution due to sediment changes. Archaea Thaumarchaeota and bacteria Gammaproteobacteria and Actinobacteria were higher in the upper sediment layers, whereas the archaea Bathyarchaeia, Lokiarchaeia, and Euryarchaeota, and bacteria Chloroflexi and Deltaproteobacteria were abundant in the deep sediments. Sediment depth and TOC were the vital factors in explaining both the bacterial and archaeal communities. Bacterial communities were abundant in microbes potentially involved in nitrogen and sulfur cycling, while archaeal communities were abundant in microbes potentially involved in ammonia oxidation and methanogenesis.

6 DATA AVAILABILITY STATEMENT

All data generated and/or analyzed during this study are available from the corresponding author upon request.

7 ACKNOWLEDGMENT

We thank all of the scientists and crew members on-board the R/V Kexue San Hao for the assistance provided in the collection of samples and geochemical data during the cruise.

References
Bagchi A, Roy D, Roy P. 2005. Homology modeling of a transcriptional regulator SoxR of the Lithotrophic sulfur oxidation (Sox) operon in α-Proteobacteria. Journal of Biomolecular Structure and Dynamics, 22(5): 571-577. DOI:10.1080/07391102.2005.10507027
Berg C, Listmann L, Vandieken V, Vogts A, Jürgens K. 2014. Chemoautotrophic growth of ammonia-oxidizing Thaumarchaeota enriched from a pelagic redox gradient in the Baltic Sea. Frontiers in Microbiology, 5: 786. DOI:10.3389/fmicb.2014.00786
Berlanga M, Aas J A, Paster B J, Boumenna T, Dewhirst F E, Guerrero R. 2008. Phylogenetic diversity and temporal variation in the Spirochaeta populations from two Mediterranean microbial mats. International Microbiology, 11(4): 267-274. DOI:10.2436/20.1501.01.71
Biddle J F, Lipp J S, Lever M A, Lloyd K G, Sorensen K B, Anderson R, Fredricks H F, Elvert M, Kelly T J, Schrag D P, Sogin M L, Brenchley J E, Teske A, House C H, Hinrichs K U. 2006. Heterotrophic Archaea dominate sedimentary subsurface ecosystems off Peru. Proceedings of the National Academy of Sciences of the United States of America, 103(10): 3846-3851. DOI:10.1073/pnas.0600035103
Boetius A, Ravenschlag K, Schubert C J, Rickert D, Widdel F, Gieseke A, Amann R, Jørgensen B B, Witte U, Pfannkuche O. 2000. A marine microbial consortium apparently mediating anaerobic oxidation of methane. Nature, 407(6804): 623-626. DOI:10.1038/35036572
Bokulich N A, Subramanian S, Faith J J, Gevers D, Gordon J I, Knight R, Mills D A, Caporaso J G. 2013. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nature Methods, 10(1): 57-59. DOI:10.1038/nmeth.2276
Brochier-Armanet C, Boussau B, Gribaldo S, Forterre P. 2008. Mesophilic crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota. Nature Reviews Microbiology, 6(3): 245-252. DOI:10.1038/nrmicro1852
Caporaso J G, Lauber C L, Walters W A, Berg-Lyons D, Huntley J, Fierer N, Owens S M, Betley J, Fraser L, Bauer M, Gormley N, Gilbert J A, Smith G, Knight R. 2012. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. The ISME Journal, 6(8): 1621-1624. DOI:10.1038/ismej.2012.8
Coolen M J L, Hopmans EC, Rijpstra WIC, Muyzer G, Schouten S, Volkman J K, Damsté JSS. 2004. Evolution of the methane cycle in Ace Lake (Antarctica) during the Holocene: response of methanogens and methanotrophs to environmental change. Organic Geochemistry, 35(10): 1151-1167. DOI:10.1016/j.orggeochem.2004.06.009
Covault J A, Fildani A. 2014. Continental shelves as sediment capacitors or conveyors: source-to-sink insights from the tectonically active Oceanside shelf, southern California, USA. Geological Society, London, Memoirs, 41(1): 315-326. DOI:10.1144/M41.23
Cowie G L, Hedges J I. 1994. Biochemical indicators of diagenetic alteration in natural organic matter mixtures. Nature, 369(6478): 304-307. DOI:10.1038/369304a0
DeLong E F. 1992. Archaea in coastal marine environments. Proceedings of the National Academy of Sciences of the United States of America, 89(12): 5685-5689. DOI:10.1073/pnas.89.12.5685
Devereux R, Mosher J J, Vishnivetskaya T A, Brown S D, Beddick D L Jr, Yates D F, Palumbo A V. 2015. Changes in northern Gulf of Mexico sediment bacterial and archaeal communities exposed to hypoxia. Geobiology, 13(5): 478-493. DOI:10.1111/gbi.12142
Dharamshi J E, Tamarit D, Eme L, Stairs C W, Martijn J, Homa F, Jørgensen S L, Spang A, Ettema T J G. 2020. Marine sediments illuminate Chlamydiae diversity and evolution. Current Biology, 30(6): 1032-1048.e7. DOI:10.1016/j.cub.2020.02.016
Dong Y, Zhao Y, Zhang W Y, Li Y, Zhou F, Liu C G, Wu Y, Liu S M, Zhang W C, Xiao T. 2014. Bacterial diversity and community structure in the East China Sea by 454 sequencing of the 16S rRNA gene. Chinese Journal of Oceanology and Limnology, 32(3): 527-541. DOI:10.1007/s00343-014-3215-2
Durbin A M, Teske A. 2012. Archaea in organic-lean and organic-rich marine subsurface sediments: an environmental gradient reflected in distinct phylogenetic lineages. Frontiers in Microbiology, 3: 168. DOI:10.3389/fmicb.2012.00168
Evans P N, Parks D H, Chadwick G L, Robbins S J, Orphan V J, Golding S D, Tyson G W. 2015. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science, 350(6259): 434-438. DOI:10.1126/science.aac7745
Fernández-Gómez B, Richter M, Schüler M, Pinhassi J, Acinas S G, González J M, Pedrós-Alió C. 2013. Ecology of marine Bacteroidetes: a comparative genomics approach. The ISME Journal, 7(5): 1026-1037. DOI:10.1038/ismej.2012.169
Fuhrman J A, Davis A A. 1997. Widespread archaea and novel bacteria from the deep sea as shown by 16S rRNA gene sequences. Marine Ecology Progress Series, 150(1): 275-285. DOI:10.3354/meps150275
Fuhrman J A. 2009. Microbial community structure and its functional implications. Nature, 459(7244): 193-199. DOI:10.1038/nature08058
Gauthier M E A, Watson J R, Degnan S M. 2016. Draft genomes shed light on the dual bacterial symbiosis that dominates the microbiome of the coral reef sponge Amphimedon queenslandica. Frontiers in Marine Science, 3: 196. DOI:10.3389/fmars.2016.00196
Gold T. 1992. The deep, hot biosphere. Proceedings of the National Academy of Sciences of the United States of America, 89(13): 6045-6049. DOI:10.1073/pnas.89.13.6045
Graham D E, Wallenstein M D, Vishnivetskaya T A, Waldrop M P, Phelps T J, Pfiffner S M, Onstott T C, Whyte L G, Rivkina E M, Gilichinsky D A, Elias D A, Mackelprang R, VerBerkmoes N C, Hettich R L, Wagner D, Wullschleger S D, Jansson J K. 2012. Microbes in thawing permafrost: the unknown variable in the climate change equation. The ISME Journal, 6(4): 709-712. DOI:10.1038/ismej.2011.163
He H, Zhen Y, Mi T Z, Fu L L, Yu Z G. 2018. Ammonia-oxidizing archaea and bacteria differentially contribute to ammonia oxidation in sediments from adjacent waters of Rushan Bay, China. Frontiers in Microbiology, 9: 116. DOI:10.3389/fmicb.2018.00116
He H, Zhen Y, Mi T Z, Yu Z G. 2016a. Community composition and abundance of ammonia-oxidizing archaea in sediments from the Changjiang Estuary and its adjacent area in the East China Sea. Geomicrobiology Journal, 33(5): 416-425. DOI:10.1080/01490451.2014.986695
He Y, Li M, Perumal V, Feng X, Fang J, Xie J, Sievert S M, Wang F. 2016b. Genomic and enzymatic evidence for acetogenesis among multiple lineages of the archaeal phylum Bathyarchaeota widespread in marine sediments. Nature Microbiology, 1(6): 16035. DOI:10.1038/nmicrobiol.2016.35
Ho A, Angel R, Veraart A J, Daebeler A, Jia Z J, Kim S Y, Kerckhof F M, Boon N, Bodelier P L E. 2016. Biotic interactions in microbial communities as modulators of biogeochemical processes: methanotrophy as a model system. Frontiers in Microbiology, 7: 1285. DOI:10.3389/fmicb.2016.01285
Inagaki F, Nunoura T, Nakagawa S, Teske A, Lever M, Lauer A, Suzuki M, Takai K, Delwiche M, Colwell F S, Nealson K H, Horikoshi K, D'Hondt S, Jørgensen B B. 2006. Biogeographical distribution and diversity of microbes in methane hydrate-bearing deep marine sediments on the Pacific Ocean Margin. Proceedings of the National Academy of Sciences of the United States of America, 103(8): 2815-2820. DOI:10.1073/pnas.0511033103
Jiang L J, Zheng Y P, Chen J Q, Xiao X, Wang F P. 2011. Stratification of archaeal communities in shallow sediments of the Pearl River Estuary, Southern China. Antonie Van Leeuwenhoek, 99(4): 739-751. DOI:10.1007/s10482-011-9548-3
Karner M B, DeLong E F, Karl D M. 2001. Archaeal dominance in the mesopelagic zone of the Pacific Ocean. Nature, 409(6819): 507-510. DOI:10.1038/35054051
Lazar C S, Baker B J, Seitz K, Hyde A S, Dick G J, Hinrichs K U, Teske A P. 2016. Genomic evidence for distinct carbon substrate preferences and ecological niches of Bathyarchaeota in estuarine sediments. Environmental Microbiology, 18(4): 1200-1211. DOI:10.1111/1462-2920.13142
Lee C K, Barbier B A, Bottos E M, McDonald I R, Cary S C. 2012. The inter-valley soil comparative survey: the ecology of Dry Valley edaphic microbial communities. The ISME Journal, 6(5): 1046-1057. DOI:10.1038/ismej.2011.170
Lenk S, Arnds J, Zerjatke K, Musat N, Amann R, Mußmann M. 2011. Novel groups of Gammaproteobacteria catalyse sulfur oxidation and carbon fixation in a coastal, intertidal sediment. Environmental Microbiology, 13(3): 758-774. DOI:10.1111/j.1462-2920.2010.02380.x
Lentini V, Gugliandolo C, Bunk B, Overmann J, Maugeri T L. 2014. Diversity of prokaryotic community at a shallow marine hydrothermal site elucidated by Illumina sequencing technology. Current Microbiology, 69(4): 457-466. DOI:10.1007/s00284-014-0609-5
Lin X J, Handley K M, Gilbert J A, Kostka J E. 2015. Metabolic potential of fatty acid oxidation and anaerobic respiration by abundant members of Thaumarchaeota and Thermoplasmata in deep anoxic peat. The ISME Journal, 9(12): 2740-2744. DOI:10.1038/ismej.2015.77
Liu C S, Zhao D F, Yan L H, Wang A J, Gu Y Y, Lee D J. 2015a. Elemental sulfur formation and nitrogen removal from wastewaters by autotrophic denitrifiers and anammox bacteria. Bioresource Technology, 191: 332-336. DOI:10.1016/j.biortech.2015.05.027
Liu J W, Liu X S, Wang M, Qiao Y L, Zheng Y F, Zhang X H. 2015b. Bacterial and archaeal communities in sediments of the North Chinese Marginal Seas. Microbial Ecology, 70(1): 105-117. DOI:10.1007/s00248-014-0553-8
Liu Y C, Whitman W B. 2008. Metabolic, phylogenetic, and ecological diversity of the methanogenic archaea. Annals of the New York Academy of Sciences, 1125(1): 171-189. DOI:10.1196/annals.1419.019
Lloyd K G, Schreiber L, Petersen D G, Kjeldsen K U, Lever M A, Steen A D, Stepanauskas R, Richter M, Kleindienst S, Lenk S, Schramm A, Jørgensen B B. 2013. Predominant archaea in marine sediments degrade detrital proteins. Nature, 496(7444): 215-218. DOI:10.1038/nature12033
Lovley D R, Klug M J. 1983. Sulfate reducers can outcompete methanogens at freshwater sulfate concentrations. Applied and Environmental Microbiology, 45(1): 187-192. DOI:10.1128/AEM.45.1.187-192.1983
Lozada M, Marcos M S, Commendatore M G, Gil M N, Dionisi H M. 2014. The bacterial community structure of hydrocarbon-polluted marine environments as the basis for the definition of an ecological index of hydrocarbon exposure. Microbes and Environments, 29(3): 269-276. DOI:10.1264/jsme2.ME14028
Magoč T, Salzberg S L. 2011. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics, 27(21): 2957-2963. DOI:10.1093/bioinformatics/btr507
Mahmoudi N, Robeson M S, Castro H F, Fortney J L, Techtmann S M, Joyner D C, Paradis C J, Pfiffner S M, Hazen T C. 2015. Microbial community composition and diversity in Caspian Sea sediments. FEMS Microbiology Ecology, 91(1): 1-11. DOI:10.1093/femsec/fiu013
Niemann H, Lösekann T, De Beer D, Elvert M, Nadalig T, Knittel K, Amann R, Sauter E J, Schlüter M, Klages M, Foucher J P, Boetius A. 2006. Novel microbial communities of the Haakon Mosby mud volcano and their role as a methane sink. Nature, 443(7113): 854-858. DOI:10.1038/nature05227
Nunoura T, Takaki Y, Kazama H, Hirai M, Ashi J, Imachi H, Takai K. 2012. Microbial diversity in deep-sea methane seep sediments presented by SSU rRNA gene tag sequencing. Microbes and Environments, 27(4): 382-390. DOI:10.1264/jsme2.ME12032
Offre P, Spang A, Schleper C. 2013. Archaea in biogeochemical cycles. Annual Review of Microbiology, 67(1): 437-457. DOI:10.1146/annurev-micro-092412-155614
Oni O E, Schmidt F, Miyatake T, Kasten S, Witt M, Hinrichs K U, Friedrich M W. 2015. Microbial communities and organic matter composition in surface and subsurface sediments of the Helgoland mud area, North Sea. Frontiers in Microbiology, 6: 1290. DOI:10.3389/fmicb.2015.01290
Orphan V J, House C H, Hinrichs K U, McKeegan K D, DeLong E F. 2002. Multiple archaeal groups mediate methane oxidation in anoxic cold seep sediments. Proceedings of the National Academy of Sciences of the United States of America, 99(11): 7663-7668. DOI:10.1073/pnas.072210299
Pala C, Molari M, Nizzoli D, Bartoli M, Viaroli P, Manini E. 2018. Environmental drivers controlling bacterial and archaeal abundance in the sediments of a Mediterranean lagoon ecosystem. Current Microbiology, 75(9): 1147-1155. DOI:10.1007/s00284-018-1503-3
Peiffer J A, Spor A, Koren O, Jin Z, Tringe S G, Dangl J L, Buckler E S, Ley R E. 2013. Diversity and heritability of the maize rhizosphere microbiome under field conditions. Proceedings of the National Academy of Sciences of the United States of America, 110(16): 6548-6553. DOI:10.1073/pnas.1302837110
Porat I, Vishnivetskaya T A, Mosher J J, Brandt C C, Yang Z K, Brooks S C, Liang L Y, Drake M M, Podar M, Brown S D, Palumbo A V. 2010. Characterization of archaeal community in contaminated and uncontaminated surface stream sediments. Microbial Ecology, 60(4): 784-795. DOI:10.1007/s00248-010-9734-2
Poulsen M, Schwab C, Jensen B B, Engberg R M, Spang A, Canibe N, Højberg O, Milinovich G, Fragner L, Schleper C, Weckwerth W, Lund P, Schramm A, Urich T. 2013. Methylotrophic methanogenic Thermoplasmata implicated in reduced methane emissions from bovine rumen. Nature Communications, 4: 1428. DOI:10.1038/ncomms2847
Qiao Y L, Liu J W, Zhao M X, Zhang X H. 2018. Sediment depth-dependent spatial variations of bacterial communities in mud deposits of the eastern China marginal seas. Frontiers in Microbiology, 9: 1128. DOI:10.3389/fmicb.2018.01128
Reed A J, Lutz R A, Vetriani C. 2006. Vertical distribution and diversity of bacteria and archaea in sulfide and methanerich cold seep sediments located at the base of the Florida Escarpment. Extremophiles, 10(3): 199-211. DOI:10.1007/s00792-005-0488-6
Robador A, Müller A L, Sawicka J E, Berry D, Hubert C R J, Loy A, Jørgensen B B, Brüchert V. 2016. Activity and community structures of sulfate-reducing microorganisms in polar, temperate and tropical marine sediments. The ISME Journal, 10(4): 796-809. DOI:10.1038/ismej.2015.157
Säwström C, Serrano O, Rozaimi M, Lavery P S. 2016. Utilization of carbon substrates by heterotrophic bacteria through vertical sediment profiles in coastal and estuarine seagrass meadows. Environmental Microbiology Reports, 8(5): 582-589. DOI:10.1111/1758-2229.12406
Schippers A, Kock D, Höft C, Köweker G, Siegert M. 2012. Quantification of microbial communities in subsurface marine sediments of the Black Sea and off Namibia. Frontiers in Microbiology, 3: 16. DOI:10.3389/fmicb.2012.00016
Schippers A, Neretin L N. 2006. Quantification of microbial communities in near-surface and deeply buried marine sediments on the Peru continental margin using real-time PCR. Environmental Microbiology, 8(7): 1251-1260. DOI:10.1111/j.1462-2920.2006.01019.x
Schönheit P, Kristjansson J K, Thauer R K. 1982. Kinetic mechanism for the ability of sulfate reducers to outcompete methanogens for acetate. Archives of Microbiology, 132(3): 285-288. DOI:10.1007/BF00407967
Shehab N, Li D, Amy G L, Logan B E, Saikaly P E. 2013. Characterization of bacterial and archaeal communities in air-cathode microbial fuel cells, open circuit and sealed-off reactors. Applied Microbiology and Biotechnology, 97(22): 9885-9895. DOI:10.1007/s00253-013-5025-4
Shi X F, Shen S X, Yi H I, Chen Z H, Meng Y. 2003. Modern sedimentary environments and dynamic depositional systems in the southern Yellow Sea. Chinese Science Bulletin, 48(S1): 1-7. DOI:10.1007/BF02900933
Spang A, Hatzenpichler R, Brochier-Armanet C, Rattei T, Tischler P, Spieck E, Streit W, Stahl D A, Wagner M, Schleper C. 2010. Distinct gene set in two different lineages of ammonia-oxidizing archaea supports the phylum Thaumarchaeota. Trends in Microbiology, 18(8): 331-340. DOI:10.1016/j.tim.2010.06.003
Spring S, Riedel T, Spröer C, Yan S, Harder J, Fuchs B M. 2013. Taxonomy and evolution of bacteriochlorophyll a-containing members of the OM60/NOR5 clade of marine gammaproteobacteria: description of Luminiphilus syltensis gen. nov., sp. nov., reclassification of Haliea rubra as Pseudohaliea rubra gen. nov., comb. nov., and emendation of Chromatocurvus halotolerans. BMC Microbiology, 13(1): 118. DOI:10.1186/1471-2180-13-118
Stahl D A, de la Torre J R. 2012. Physiology and diversity of ammonia-oxidizing archaea. Annual Review of Microbiology, 66(1): 83-101. DOI:10.1146/annurev-micro-092611-150128
Sun W M, Xiao E Z, Pu Z L, Krumins V, Dong Y R, Li B Q, Hu M. 2018. Paddy soil microbial communities driven by environment- and microbe-microbe interactions: a case study of elevation-resolved microbial communities in a rice terrace. Science of the Total Environment, 612: 884-893. DOI:10.1016/j.scitotenv.2017.08.275
Suzuki D, Li Z L, Cui X X, Zhang C F, Katayama A. 2014. Reclassification of Desulfobacterium anilini as Desulfatiglans anilini comb. nov. within Desulfatiglans gen. nov., and description of a 4-chlorophenol-degrading sulfate-reducing bacterium, Desulfatiglans parachlorophenolica sp. nov. International Journal of Systematic and Evolutionary Microbiology, 64(Pt 9): 3081-3086. DOI:10.1099/ijs.0.064360-0
Timmers P H A, Suarez-Zuluaga D A, van Rossem M, Diender M, Stams A J M, Plugge C M. 2016. Anaerobic oxidation of methane associated with sulfate reduction in a natural freshwater gas source. The ISME Journal, 10(6): 1400-1412. DOI:10.1038/ismej.2015.213
Tringe S G, Hugenholtz P. 2008. A renaissance for the pioneering 16S rRNA gene. Current Opinion in Microbiology, 11(5): 442-446. DOI:10.1016/j.mib.2008.09.011
Urakawa H, Yoshida T, Nishimura M, Ohwada K. 2000. Characterization of depth-related population variation in microbial communities of a coastal marine sediment using 16S rDNA-based approaches and quinone profiling. Environmental Microbiology, 2(5): 542-554. DOI:10.1046/j.1462-2920.2000.00137.x
Valentine D L. 2002. Biogeochemistry and microbial ecology of methane oxidation in anoxic environments: a review. Antonie van Leeuwenhoek, 81(1-4): 271-282. DOI:10.1023/A:1020587206351
Varon-Lopez M, Dias A C F, Fasanella C C, Durrer A, Melo I S, Kuramae E E, Andreote F D. 2014. Sulphur-oxidizing and Sulphate-reducing communities in Brazilian mangrove sediments. Environmental Microbiology, 16(3): 845-855. DOI:10.1111/1462-2920.12237
Wakeham S G, Lee C, Hedges J I, Hernes P J, Peterson M J. 1997. Molecular indicators of diagenetic status in marine organic matter. Geochimica et Cosmochimica Acta, 61(24): 5363-5369. DOI:10.1016/S0016-7037(97)00312-8
Walsh J J. 1991. Importance of continental margins in the marine biogeochemical cycling of carbon and nitrogen. Nature, 350(6313): 53-55. DOI:10.1038/350053a0
Wang Y, Sheng H F, He Y, Wu J Y, Jiang Y X, Tam N F Y, Zhou H W. 2012. Comparison of the levels of bacterial diversity in freshwater, intertidal wetland, and marine sediments by using millions of Illumina tags. Applied and Environmental Microbiology, 78(23): 8264-8271. DOI:10.1128/AEM.01821-12
Wilms R, Kopke B, Sass H, Chang T S, Cypionka H, Engelen B. 2006. Deep biosphere-related bacteria within the subsurface of tidal flat sediments. Environmental Microbiology, 8(4): 709-719. DOI:10.1111/j.1462-2920.2005.00949.x
Xiao K Q, Beulig F, Røy H, Jørgensen B B, Risgaard-Petersen N. 2018. Methylotrophic methanogenesis fuels cryptic methane cycling in marine surface sediment. Limnology and Oceanography, 63(4): 1519-1527. DOI:10.1002/lno.10788
Xiong J B, Ye X S, Wang K, Chen H P, Hu C J, Zhu J L, Zhang D M. 2014. Biogeography of the sediment bacterial community responds to a nitrogen pollution gradient in the East China Sea. Applied and Environmental Microbiology, 80(6): 1919-1925. DOI:10.1128/AEM.03731-13
Yamada T, Sekiguchi Y, Hanada S, Imachi H, Ohashi A, Harada H, Kamagata Y. 2006. Anaerolinea thermolimosa sp. nov., Levilinea saccharolytica gen. nov., sp. nov. and Leptolinea tardivitalis gen. nov., sp. nov., novel filamentous anaerobes, and description of the new classes Anaerolineae classis nov. and Caldilineae classis nov. in the bacterial phylum Chloroflexi. International Journal of Systematic and Evolutionary Microbiology, 56(6): 1331-1340. DOI:10.1099/ijs.0.64169-0
Yu S L, Yao P, Liu J W, Zhao B, Zhang G L, Zhao M X, Yu Z G, Zhang X H. 2016. Diversity, abundance, and niche differentiation of ammonia-oxidizing prokaryotes in mud deposits of the eastern China marginal Seas. Frontiers in Microbiology, 7: 137. DOI:10.3389/fmicb.2016.00137
Zhang Y, Wang X G, Zhen Y, Mi T Z, He H, Yu Z G. 2017. Microbial diversity and community structure of sulfate-reducing and sulfur-oxidizing bacteria in sediment cores from the East China Sea. Frontiers in Microbiology, 8: 2133. DOI:10.3389/fmicb.2017.02133
Zhou Z C, Meng H, Liu Y, Gu J D, Li M. 2017. Stratified bacterial and archaeal community in mangrove and intertidal wetland mudflats revealed by high throughput 16S rRNA gene sequencing. Frontiers in Microbiology, 8: 2148. DOI:10.3389/fmicb.2017.02148
Zhu D C, Tanabe S H, Yang C, Zhang W M, Sun J Z. 2013. Bacterial community composition of South China Sea sediments through pyrosequencing-based analysis of 16S rRNA genes. PLoS One, 8(10): e78501. DOI:10.1371/journal.pone.0078501
Zhuang G C, Heuer V B, Lazar C S, Goldhammer T, Wendt J, Samarkin V A, Elvert M, Teske A P, Joye S B, Hinrichs K U. 2018. Relative importance of methylotrophic methanogenesis in sediments of the western Mediterranean Sea. Geochimica et Cosmochimica Acta, 224: 171-186. DOI:10.1016/j.gca.2017.12.024
Zinger L, Amaral-Zettler L A, Fuhrman J A, Horner-Devine M C, Huse S M, Welch D B M, Martiny J B H, Sogin M, Boetius A, Ramette A. 2011. Global patterns of bacterial beta-diversity in seafloor and seawater ecosystems. PLoS One, 6(9): e24570,. DOI:10.1371/journal.pone.0024570
Zonneveld K A F, Versteegh G J M, Kasten S, Eglinton T I, Emeis K C, Huguet C, Koch B P, de Lange G J, de Leeuw J W, Middelburg J J, Mollenhauer G, Prahl F G, Rethemeyer J, Wakeham S G. 2010. Selective preservation of organic matter in marine environments; processes and impact on the sedimentary record. Biogeosciences, 7(2): 483-511. DOI:10.5194/bg-7-483-2010