Abstract
Ecological evidence suggests that heterotrophic diazotrophs fueled by organic carbon respiration in sediments play an important role in marine nitrogen fixation. However, fundamental knowledge about the identities, abundance, diversity, biogeography, and controlling environmental factors of nitrogen-fixing communities in open ocean sediments is still elusive. Surprisingly, little is known also about nitrogen-fixing communities in sediments of the more research-accessible marginal seas. Here we report on an investigation of the environmental geochemistry and putative diazotrophic microbiota in the sediments of Bohai Sea, an eutrophic marginal sea of the western Pacific Ocean. Diverse and abundant nifH gene sequences were identified and sulfate-reducing bacteria (SRB) were found to be the dominant putative nitrogen-fixing microbes. Community statistical analyses suggested bottom water temperature, bottom water chlorophyll a content (or the covarying turbidity) and sediment porewater Eh (or the covarying pH) as the most significant environmental factors controlling the structure and spatial distribution of the putative diazotrophic communities, while sediment Hg content, sulfide content, and porewater -Si content were identified as the key environmental factors correlated positively with the nifH gene abundance in Bohai Sea sediments. Comparative analyses between the Bohai Sea and the northern South China Sea (nSCS) identified a significant composition difference of the putative diazotrophic communities in sediments between the shallow-water (estuarine and nearshore) and deep-water (offshore and deep-sea) environments, and sediment porewater dissolved oxygen content, water depth and in situ temperature as the key environmental factors tentatively controlling the species composition, community structure, and spatial distribution of the marginal sea sediment nifH-harboring microbiota. This confirms the ecophysiological specialization and niche differentiation between the shallow-water and deep-water sediment diazotrophic communities and suggests that the in situ physical and geochemical conditions play a more important role than geographical contiguity in determining the community similarity of the diazotrophic microbiota in marginal sea sediments.
Introduction
N2-fixing prokaryotes play a key role in marine nitrogen (N) cycling and ecosystem functioning such as carbon sequestration by providing newly fixed nitrogenous nutrients, particularly important in oligotrophic environments of the ocean (Karl et al., ; Dore et al., ; Montoya et al., ; Steppe and Paerl, ; Karl and Letelier, ; Sohm et al., ). Marine N2 fixation studies have been conducted for more than half a century, yet many diazotrophic microbes, their ecophysiology and environmental response have been revealed only in recent decades (reviewed by Zehr and Kudela, ; Voss et al., ). The marine N cycle appears to be a conundrum because the estimated N input by N2 fixation is significantly lower than the estimated N loss via denitrification and anaerobic ammonium oxidation (anammox) (Mahaffey et al., ). This suggested that the inventory of reactive nitrogen in the oceans is unbalanced and dwindling (Moisander et al., ). However, this conundrum may be an artifact caused by uncertainties in previous research results and underestimations of N2-fixation contribution to the marine N budget. These include insufficient and inaccurate measurements of the N2 fixation rate (Deutsch et al., ; Großkopf et al., ), undiscovered N2-fixing microbes (Zehr et al., , ; Pernthaler et al., ; Dekas et al., ; Zehr, ; Voss et al., ), overlooked robustness of N2-fixing physiology of environmental diazotrophs (Knapp, ), and overlooked and undiscovered N2-fixing environments such as nutrient-rich estuarine and coastal seas and marine sediments (Mehta et al., ; Pernthaler et al., ; Bonnet et al., ; Voss et al., ). Historically, marine N2 fixation was thought to be carried out mainly by aggregate-forming cyanobacteria such as Trichodesmium and important only in surface and near-surface waters of the ocean (Paerl, ; Dore et al., ; Karl et al., ; Bergman et al., ). However, the water column below the photic layer has been found to harbor significant N2 fixation activities in many marine waters, which were mainly carried out by heterotrophic diazotrophs either in hypoxic and anoxic environments or in marine particle-associated microenvironments (Farnelid et al., , ; Fernandez et al., ; Jayakumar et al., ; Bird and Wyman, ; Rahav et al., ; Dang and Lovell, ). Furthermore, the recent detection of numerous diazotrophic bacteria and archaea in marine sediments led to the hypothesis that marine sediments might constitute an important environment for N2 fixation in the oceans (Pernthaler et al., ; Dang et al., , ; Dekas et al., , , ; Fulweiler, ; Miyazaki et al., ).
Hypoxic and anoxic environments can be formed and maintains more easily in marine sediments than in the water column, particularly in eutrophic coastal seas. Nitrogen gas is abundant in marine sediments, partially provided by microbial denitrification and anammox processes (Dang et al., ; Trimmer and Nicholls, ; Shao et al., ). Therefore, heterotrophic rather than cyanobacterial diazotrophy may play a dominant role of N2 fixation in marine sediments (Bertics et al., ; Dekaezemacker et al., ). N2 fixation by heterotrophic diazotrophs requires a high amount of cellular energy that is mainly provided by respiration of large amounts of organic carbon (Shanmugam et al., ; Brill, ; Dang and Jiao, ). It has been speculated that the low N2-fixing rates of heterotrophic diazotrophs in the open ocean is caused by the lack of sufficient metabolic energy due to the scarcity of bioavailable organic carbon (Moisander et al., ). On the other hand, estuarine and coastal sediments may support high rates of N2 fixation by heterotrophic diazotrophs due to enhanced phytoplankton labile organic matter production under the impact of anthropogenic eutrophication and terrigenous nutrient inputs (Herbert, ; Boesch, ; Smith, ). Moreover, iron, phosphorus or both have been found to be the key factors limiting microbial N2 fixation and export production in many waters of the oligotrophic open oceans (Sañudo-Wilhelmy et al., ; Mills et al., ; Moore et al., , ; Boyd and Ellwood, ; Sohm et al., ; Jacq et al., ). In estuarine and coastal sediments, iron and phosphorus may be sufficiently abundant to support high rates of microbial N2 fixation (Street and Paytan, ; Homoky et al., ; Karl, ). Furthermore, it is evident that the sediment diazotrophic communities are highly resistant to the inhibition of high environmental and concentrations (McGlathery et al., ; Knapp, ; Bertics et al., ). It was thus reasonably hypothesized that the marginal sea sediments may prove to be key environments of N2 fixation. However, fundamental knowledge about the identities, diversity, biogeography, and controlling environmental factors of sediment N2-fixing microbes is still lacking.
It has been reported that the sediments of the northern South China Sea (nSCS), a large and relatively oligotrophic marginal sea of the western Pacific Ocean, might harbor the highest diazotroph diversity among all the marine environments ever investigated using molecular ecology approaches targeting the nifH gene (encoding the nitrogenase reductase subunit) (Dang et al., ). The Bohai Sea, another marginal sea of the western Pacific Ocean, is conversely characterized by its eutrophic status caused by high degree of river inputs and anthropogenic pollutions and by its low water exchange with the outer ocean due to its semi-enclosed topography (Figure 1; Dang et al., ). To answer the question how diazotrophic communities vary and respond to distinct environmental conditions, sediment nifH-harboring microbial assemblages in the Bohai Sea were investigated in an environmental geochemistry context and comparatively analyzed against the putative diazotrophic assemblages of the nSCS.
Figure 1
Materials and methods
Site description, sample collection, and environment factor measurements
Bohai Sea is a large shallow water basin with an area of 77 × 103 km2 and an average water depth of only 18 m. It is nearly enclosed and thus its water exchange capacity with the outer Yellow Sea of the western Pacific Ocean is highly limited (Figure 1). Moreover, more than 40 rivers discharge into the Bohai Sea. Thus, the Bohai Sea receives intense and extensive terrigenous and anthropogenic impacts, especially in its estuarine and coastal bay areas (PEMSEA and BSEMP, ).
Samples of surface sediments in the top 0–5 cm layer were collected from 14 stations of the Bohai Sea (Figure 1) in the August of 2008 as described previously (Dang et al., ). At each station, seawater physicochemical and biological parameters (Supplementary Table S1) were measured in situ at various water depths with a Compact-CTD equipped with a TCDKU sensor (Alec Electronics, Japan) as reported in a previous publication (Dang et al., ), which also contained most of the sediment and sediment porewater environmental parameters, except the sediment V and S contents that were measured in the current study via a VARIAN 725-ES inductively coupled plasma-optical emission spectrometer (Varian, Palo Alto, CA, USA) (Supplementary Table S1).
DNA extraction and nifH gene clone library construction and analyses
Community genomic DNA of sediment microbes was extracted by using a FastPrep DNA Extraction Kit for Soil and a FastPrep-24 Cell Disrupter (MP Biomedicals, Solon, OH, USA) as previously described (Dang et al., ,). DNA concentrations were measured with dye PicoGreen (Molecular Probes, Eugene, OR, USA) and a Modulus Single Tube Multimode Reader fluorometer (Turner BioSystems, Sunnyvale, CA, USA). Partial nifH gene sequences were amplified with PCR primers nifHfw and nifHrv, which have a very broad coverage for both bacterial and archaeal nifH genes (Mehta et al., ). In order to test reproducibility of our experimental procedure and to identify any potential within-site variability of the sediment diazotrophic community, two separate nifH gene clone libraries (B22-1 and B22-2) were constructed for the sediments of sampling station B22, each from a distinct subcore DNA sample, respectively. PCR product cloning followed previous procedures (Dang et al., , ). A miniprep method was used for recombinant plasmid extractions (Dang and Lovell, ). Cloned gene fragments were reamplified to check the correct size of the DNA inserts using cloning vector PCR primers M13-D and RV-M (Dang et al., ), which were also used for sequencing with an ABI 3770 sequencer (Applied BioSystems, Foster City, CA, USA). The obtained nifH gene sequences have been submitted to GenBank with accession numbers KM524369 to KM525656.
The nifH sequences were translated into conceptual protein sequences and the BLASTp program was used for retrieval of the top-hit NifH sequences from GenBank (last accessed 25 January 2014) (Altschul et al., ). The NifH protein sequences were grouped into operational taxonomic units (OTUs) using a 0.05 sequence distance cutoff calculated by using the DOTUR program (Schloss and Handelsman, ). These OTU sequences, along with reference sequences retrieved from GenBank, were aligned with program CLUSTAL X (version 2.1; Larkin et al., ), and used for phylogeny inference with the distance and neighbor-joining method implemented within the PHYLIP software package (version 3.69; Felsenstein, ), following a previously reported procedure (Dang et al., , ).
Quantification of sediment nifH gene copy numbers
The technique of real-time fluorescent quantitative PCR (qPCR) was employed to measure the abundance of the sediment nifH genes using primers nifHfw and nifHrv (Mehta et al., ; Dang et al., ). Triplicate sediment DNA samples from each station were assayed with an ABI Prism 7500 Sequence Detection System (Applied Biosystems, Foster City, CA, USA), following a previously published SYBR Green qPCR protocol (Dang et al., ). Plasmids containing the nifH gene inserts were extracted from Escherichia coli hosts using a Mini Plasmid Kit (Qiagen, Valencia, CA, USA), linearized with an endonuclease specific in the cloning vector region, and quantified using PicoGreen and a Modulus Single Tube Multimode Reader fluorometer (Dang et al., ). A qPCR standard curve was then generated with serially diluted linearized plasmids containing the target nifH gene fragment.
The qPCR condition for nifH gene quantification was optimized based on a previous study (Dang et al., ). The efficiency and sensitivity of the qPCR were shown in Supplementary Table S2, along with the range of the plasmid copy numbers used for the qPCR standard curve construction. In all experiments, negative controls lacking template DNA were run with the same qPCR procedure to prevent contamination or carryover. Agarose gel electrophoresis and melting curve analysis were performed to confirm the specificity of the qPCR reactions.
Statistical analyses
Clone library coverage (C) was calculated as C = [1 − (n1/N)] × 100, where n1 is the number of unique OTUs and N the total number of clones in a library (Mullins et al., ). Shannon-Wiener (H), Simpson (D), and evenness (J) indices were calculated with the OTU data (Dang et al., ). The DOTUR program (Schloss and Handelsman, ) was used for clone library rarefaction analysis and to calculate the abundance-based coverage estimator (SACE) and the bias-corrected richness estimator Chao1 (SChao1). The Fast UniFrac program (Hamady et al., ) was used for nifH-harboring microbial community clustering and principal coordinates analyses (PCoA). The Canoco software (version 4.5, Microcomputer Power, Ithaca, NY, USA) was used for the canonical correspondence analysis (CCA) (ter Braak and Šmilauer, ; Lepš and Šmilauer, ), to investigate any correlations between the nifH-harboring microbial communities and environmental factors by following a previously described procedure (Dang et al., ). The SPSS software (version 17.0) was used for Pearson correlation analyses (significance level α = 0.05) of the sediment nifH gene abundance with environmental factors (Dang et al., ).
The Fast UniFrac program (Hamady et al., ) was also used for clustering analyses of the sediment nifH-harboring microbial communities in both the Bohai Sea and the previously studied nSCS (Dang et al., ), to identify larger geographical-scale characteristics of the putative diazotrophic microbiota. CCA analysis (ter Braak and Šmilauer, ; Lepš and Šmilauer, ) was also performed to tentatively identify the key environmental factors that may control the biogeographical distribution of the sediment nifH-harboring microbial communities in both marginal seas.
Results
Diversity of nifH gene sequences from bohai sea sediments
Community clustering (Supplementary Figure S1) and PCoA (Supplementary Figure S2) analyses using Fast UniFrac software revealed that the two nifH gene clone libraries (B22-1 and B22-2) constructed from separate sediment subcore samples of station B22 were highly similar. These statistical results demonstrated the reproducibility of our experimental procedures and the negligibility of within-site variability of the putative diazotrophic community detected by using the nifH gene biomarker. Therefore, in subsequent analyses, these two nifH gene clone libraries (B22-1 and B22-2) were pooled as a single B22 clone library.
Of all the 14 nifH gene clone libraries constructed, a total of 1309 clones were found to contain a valid nifH gene fragment, which were further identified as 1288 unique nifH DNA sequences yielding 619 deduced unique NifH protein sequences and 246 OTUs. Based on the obtained biodiversity index values (Table 1), the sediment nifH genes were highly diverse and their richness was heterogeneously distributed in the different sampling stations of the Bohai Sea. The high diversity (Supplementary Figure S3, Table 1) indicated that the sequences in the constructed nifH gene clone libraries might only represent the dominant putative diazotrophs in the Bohai Sea sediments. In general, the coastal site B15 had the lowest nifH gene diversity and the sites B5, B16, B20, and B22 might have the highest nifH gene diversity (Supplementary Figure S3, Table 1).
Table 1
| Station | No. of clones | No. of unique sequences | No. of OTUs | C (%) | H | 1/D | J | SACE | SChao1 |
|---|---|---|---|---|---|---|---|---|---|
| B1 | 82 | 63 | 40 | 69.51 | 4.86 | 28.38 | 0.91 | 94.49 | 77.50 |
| B3 | 80 | 58 | 44 | 63.75 | 5.05 | 31.60 | 0.92 | 96.33 | 102.00 |
| B5 | 85 | 59 | 48 | 61.18 | 5.17 | 36.06 | 0.93 | 141.89 | 114.00 |
| B7 | 87 | 58 | 35 | 73.56 | 4.39 | 15.72 | 0.86 | 101.42 | 77.17 |
| B8 | 92 | 63 | 44 | 71.74 | 5.08 | 36.40 | 0.93 | 86.16 | 90.43 |
| B10 | 88 | 58 | 41 | 71.59 | 4.82 | 23.92 | 0.90 | 96.60 | 83.86 |
| B11 | 83 | 54 | 34 | 73.49 | 4.35 | 14.80 | 0.86 | 97.45 | 72.50 |
| B12 | 89 | 61 | 37 | 74.16 | 4.51 | 15.92 | 0.87 | 75.28 | 79.17 |
| B14 | 85 | 55 | 36 | 75.29 | 4.53 | 16.76 | 0.88 | 73.26 | 62.25 |
| B15 | 91 | 57 | 27 | 84.62 | 3.89 | 10.06 | 0.82 | 54.73 | 45.20 |
| B16 | 93 | 71 | 47 | 68.82 | 5.14 | 36.56 | 0.93 | 104.87 | 87.60 |
| B19 | 86 | 63 | 37 | 73.26 | 4.60 | 19.86 | 0.88 | 91.17 | 73.14 |
| B20 | 83 | 61 | 48 | 59.04 | 5.17 | 35.82 | 0.93 | 154.44 | 118.13 |
| B22 | 185 | 111 | 66 | 78.38 | 5.22 | 24.38 | 0.86 | 143.74 | 152.67 |
| B22-1 | 96 | 70 | 50 | 66.67 | 5.16 | 32.11 | 0.91 | 126.55 | 91.30 |
| B22-2 | 89 | 65 | 40 | 70.79 | 4.73 | 21.52 | 0.89 | 98.8 | 105.00 |
Biodiversity and predicted richness of the sediment nifH gene sequences obtained from the sampling stations of the Bohai Sea.
Phylogeny of deduced NifH protein sequences from bohai sea sediments
The obtained 1288 unique nifH gene sequences shared 37.5–99.7% sequence identity with one another and 51–99% sequence identity with the top-match sequences obtained from GenBank. Interestingly, the majority (87.7%) of the queried Bohai Sea sediment nifH gene sequences resulted in nifH top-match sequences derived from samples obtained from nSCS marine sediments (Dang et al., ). The deduced 619 unique NifH protein sequences shared 33.0–99.2% sequence identity with one another and shared 53–100% sequence identity with the top-match GenBank NifH sequences. More than half (57.4%) of the 619 unique NifH sequences shared >90% sequence identity with known bacteria, such as Alphaproteobacteria, Gammaproteobacteria, Deltaproteobacteria, Chlorobi, Firmicutes, and Verrucomicrobia (Supplementary Table S3). Importantly, most of these bacteria with top-match NifH sequences are known diazotrophs. This result indicates that diverse obligately anaerobic sulfate-reducing bacteria (SRB) in the Deltaproteobacteria class, which accounted for 18.3% of the total OTUs and 48.2% of the total clones, probably constituted the majority of taxa of diazotrophic communities in Bohai Sea sediments (Supplementary Table S3, Supplementary Figure S4).
The deduced NifH sequences were affiliated with six major groups in the constructed NifH phylogenetic tree (Figure 2, Supplementary Figure S4, Table 2), according to a recently proposed NifH sequence phylogenetic classification (Dang et al., , ). Class I contained 141 OTUs and could be further divided into four clusters (Clusters I, II, III, and Cluster IV subcluster A), in which Cluster III was the largest containing 102 OTUs (Supplementary Figure S4). Clusters I and III were the only clusters that contained the NifH sequences obtained from all 14 sampling stations; therefore, these sequences tentatively represent the most prevalent diazotrophic microbes in the Bohai Sea sediments. Sequence Classes II, III, V, VI, and VIII represented 19, 32, 10, 43, and 1 OTUs (Supplementary Figure S4), respectively, and the pertinent putative diazotrophic microbes might thus be specifically distributed only in certain sediment environments of the Bohai Sea.
Figure 2
Table 2
| Station | Class I | Class II | Class III | Class V | Class VI | Class VIII | |||
|---|---|---|---|---|---|---|---|---|---|
| Cluster I | Cluster II | Cluster III | Cluster IV(A) | Cluster IV(B) | Cluster V | Cluster VI(A) | Cluster VI(B) | ||
| B1 | 28 (7) | 41 (22) | 1 (1) | 2 (2) | 6 (4) | 2 (2) | 2 (2) | ||
| B3 | 24 (7) | 37 (19) | 6 (6) | 4 (4) | 6 (5) | 1 (1) | 2 (2) | ||
| B5 | 23 (10) | 33 (22) | 3 (1) | 11 (5) | 2 (2) | 3 (2) | 10 (6) | ||
| B7 | 38 (9) | 41 (19) | 4 (4) | 4 (3) | |||||
| B8 | 17 (5) | 1 (1) | 37 (18) | 2 (1) | 4 (2) | 21 (10) | 3 (1) | 7 (6) | |
| B10 | 26 (5) | 45 (23) | 4 (2) | 5 (3) | 1 (1) | 6 (6) | 1 (1) | ||
| B11 | 30 (6) | 40 (16) | 2 (2) | 6 (6) | 5 (4) | ||||
| B12 | 21 (4) | 49 (22) | 1 (1) | 14 (8) | 3 (1) | 1 (1) | |||
| B14 | 29 (7) | 47 (22) | 1 (1) | 4 (2) | 4 (4) | ||||
| B15 | 14 (5) | 43 (14) | 4 (2) | 2 (2) | 25 (2) | 3 (2) | |||
| B16 | 21 (7) | 59 (28) | 1 (1) | 3 (3) | 1 (1) | 8 (7) | |||
| B19 | 19 (8) | 57 (21) | 2 (2) | 1 (1) | 6 (4) | 1 (1) | |||
| B20 | 28 (10) | 39 (25) | 1 (1) | 3 (2) | 9 (7) | 3 (3) | |||
| B22 | 35 (10) | 128 (44) | 2 (2) | 1 (1) | 16 (6) | 1 (1) | 2 (2) | ||
The detail composition of the sequences (and OTUs) in the NifH phylogenetic tree obtained from the Bohai Sea sediments.
Key factors controlling the nifH-harboring microbial communities
Community classification based on Fast UniFrac clustering analysis identified three distinct clusters of the nifH-harboring microbial assemblages in Bohai Sea sediments (Figure 3). The nifH-harboring microbial assemblages of stations B8 and B15 were grouped together and represented a distinct cluster, the assemblage of station B5 represented another distinct cluster formed by a singleton member, and the assemblages of the remaining stations were grouped together and formed a third distinct cluster. This community classification pattern of nifH-harboring microbial assemblages in the Bohai Sea sediments was further supported by the PCoA analysis (Figure 4).
Figure 3
Figure 4
CCA analysis was performed to decode the putative diazotroph-environment relationship in the Bohai Sea sediments (Supplementary Figure S5). Bottom water temperature (p = 0.001; 1000 Monte Carlo permutations), bottom water chlorophyll a content (or the covarying bottom water turbidity) (p = 0.015; 1000 Monte Carlo permutations) and sediment porewater Eh (or the covarying porewater pH; p = 0.032; 1000 Monte Carlo permutations) were identified as the most significant environmental factors that might control the structure and spatial distribution of the nifH-harboring microbial communities in Bohai Sea sediments.
Abundance of the nifH-harboring microbes in bohai sea sediments
The qPCR results showed that the abundance of the nifH genes ranged from 2.95 × 107 copies g−1 sediment (station B3) to 2.63 × 109 copies g−1 sediment (station B19) in the Bohai Sea (Table 3). The highest nifH gene abundance occurred at the sampling station B1, B19, and B14, respectively, in Liaodong Bay, Laizhou Bay, and Bohai Bay sediments, respectively. Our previous study showed that the total bacterial abundance was also heterogeneously distributed in the sediments of the Bohai Sea, with the bacterial 16S rRNA genes ranging from 3.25 × 109 copies g−1 sediment (station B16) to 2.10 × 1010 copies g−1 sediment (station B20) (Dang et al., ). The ratios of the nifH gene abundance to the bacterial 16S rRNA gene abundance ranged from 0.36% (station B1) to 14.95% (station B19).
Table 3
| Sampling station | Mean copy no. of target genes g−1sediment (SE) | |
|---|---|---|
| nifH | Bacterial 16S rRNA* | |
| B1 | 7.21 × 107 (2.70 × 106) | 1.99 × 1010 (1.30 × 109) |
| B3 | 2.95 × 107 (2.65 × 106) | 5.04 × 109 (2.09 × 108) |
| B5 | 5.99 × 107 (5.75 × 106) | 5.17 × 109 (1.71 × 108) |
| B7 | 1.45 × 108 (3.07 × 106) | 8.00 × 109 (6.19 × 108) |
| B8 | 1.22 × 108 (9.57 × 106) | 1.86 × 1010 (1.60 × 109) |
| B10 | 2.77 × 108 (1.01 × 107) | 1.45 × 1010 (1.33 × 109) |
| B11 | 7.42 × 108 (4.26 × 107) | 1.66 × 1010 (5.43 × 108) |
| B12 | 8.70 × 108 (7.46 × 107) | 1.33 × 1010 (2.82 × 108) |
| B14 | 9.28 × 108 (5.40 × 107) | 1.16 × 1010 (1.16 × 109) |
| B15 | 4.34 × 108 (3.88 × 107) | 1.06 × 1010 (8.64 × 108) |
| B16 | 2.42 × 108 (1.90 × 107) | 3.25 × 109 (4.08 × 108) |
| B19 | 2.63 × 109 (2.22 × 108) | 1.76 × 1010 (7.52 × 108) |
| B20 | 9.45 × 108 (3.87 × 107) | 2.10 × 1010 (6.29 × 108) |
| B22 | 5.54 × 108 (4.14 × 107) | 1.83 × 1010 (6.02 × 108) |
The abundance of the nifH and bacterial 16S rRNA genes in the sediments of the Bohai Sea.
The bacterial 16S rRNA gene abundance data were obtained from a previous study (Dang et al., ).
Pearson correlation analyses indicated that the sediment Hg content (p = 0.018), sulfide content (p = 0.046), and porewater -Si content (p = 0.026) were positively correlated with the nifH gene abundance in all investigated Bohai Sea sediments.
Comparison of nifH-harboring microbiota in bohai sea and nSCS sediments
In order to identify the difference of the diazotrophic microbial communities in sediments from distinct marine environments and to detect any general ecological characteristics of the diazotrophic microbiota in sediments of different marginal seas of the western Pacific Ocean, comparative analyses were made between sediment samples from the Bohai Sea and the previously studied nSCS (Dang et al., ). Community classification analyses indicated that the nifH-harboring microbial communities could be classified into two groups based on both unweighted (Figure 5A) and weighted and normalized (Figure 5B) Fast UniFrac clustering analyses, with Group I containing exclusively the nifH-harboring microbial assemblages at the nSCS offshore and deep-sea sampling stations (water depth > 130 m) and Group II containing the nifH-harboring microbial assemblages at all the Bohai Sea sampling stations as well as the shallow-water sampling stations A3 and E501 of the nSCS (water depth < 70 m). The estuarine and nearshore sediments probably harbor N2-fixing microbial assemblages distinctly different from those of the offshore and deep-sea sediments, regardless of the large geographical distance between the Bohai Sea and the nSCS. The CCA analysis further identified the sediment porewater dissolved oxygen content (DO), water depth and in situ sediment temperature (the bottom water temperature was used to approximate the sediment temperature of the Bohai Sea) as the key environmental factors (p = 0.001, p = 0.002, and p = 0.014, respectively; 1000 Monte Carlo permutations) influencing significantly the community structure and biogeographical distribution of the putative sediment diazotrophic microbiota in Bohai Sea and nSCS sediments (Figure 6).
Figure 5
Figure 6
Discussion
Nitrogenous nutrient-rich environments such as estuarine and coastal seawaters and marine sediments have long been regarded as environments lacking significant diazotrophic activities, inferred previously from bacterial culture-based physiological studies. This inferred popular opinion resulted in a long-time negligence of the N2-fixing microorganisms and activities in these environments whereas newer work led to the hypothesis that marginal sea sediments may instead harbor diverse and abundant N2-fixing microorganisms (Knapp, ; Bertics et al., ; Dang et al., ; Dekaezemacker et al., ; Voss et al., ). To test this hypothesis, we probed the genetic potential of N2 fixation in marine sediment microbial communities: Indeed, highly diverse and abundant nifH gene sequences were obtained from the sediments of the Bohai Sea (Supplementary Figure S4, Table 2), consistent with the results obtained recently from marine sediments of the nSCS (Dang et al., ). Many of the Bohai Sea nifH gene sequences are related to known bacterial diazotrophs (Supplementary Figure S4, Supplementary Table S3). While these results revealed the N2-fixation potential of microbial communities in Bohai Sea sediments, they also suggested microbial diazotrophy as a common ecological property of marginal sea sediments.
Besides being majorly related to environmental nifH gene sequences obtained mainly from the nSCS sediments (Dang et al., ) and occasionally from salt marsh sediments (Lovell et al., ; Lovell and Davis, ), many of the nifH gene sequences from Bohai Sea sediments are related to diverse known bacterial species and a small number of methanogenic archaea (Supplementary Figure S4). This is consistent with previous findings that diazotrophic bacteria may be the predominant N2-fixing microbes in shallow water sediments (Burns et al., ). A recently defined cluster of archaeal nifH gene sequences, called Cluster III-x or Methane Seep Group, was associated with microbes frequently found in deep-sea methane seep, gas hydrate, and mud volcano sediments of the Pacific Ocean marginal seas (Pernthaler et al., ; Dang et al., , ; Dekas et al., , , ; Miyazaki et al., ). However, Cluster III-x nifH gene sequences were not detected in samples obtained from Bohai Sea sediments. Because Cluster III-x nifH gene sequences are affiliated with anaerobic methane-oxidizing archaea (ANME), they may occur only in methane-rich environments. None of our Bohai Sea sampling stations were located in methane seeps, gas hydrates or mud volcanoes. The lack of the chemoautotrophic ANME nifH gene sequences in samples from shallow coastal sediments is consistent with the in situ geochemical condition of the Bohai Sea.
More than half of the deduced NifH sequences associated with Bohai Sea sediment samples shared high (>90%) identity with NifH sequences from known bacteria, of which the majority are SRB. This result indicates that SRB may be the dominant and prevalent N2-fixing microbes in the sediments of the Bohai Sea. Many SRB harbor functional N2 fixation genetic inventories (Barton and Fauque, ) and SRB were previously found to be important N2-fixing bacteria in marine intertidal microbial mats (Zehr et al., ; Steppe and Paerl, , ; Stal et al., ). It was recently reported that the benthic microbial N2 fixation rate is coupled to sulfate reduction activities in lagoons and coastal bays supporting the notion that SRB may play a key role in sediment N2 fixation and new nitrogen production (Bertics et al., , ). Furthermore, N2 fixation activities of sediment SRB are tolerant of high (e.g., up to 0.8 or 1.2 mM) environmental concentrations (McGlathery et al., ; Bertics et al., ) suggesting that benthic N2 fixation is robust and widespread in marine sediments. In our current study, sediment sulfide, the product of microbial sulfate reduction, was found to correlate positively with nifH gene abundance in Bohai Sea sediments. This suggests that SRB may contribute substantially to N2 fixation and sulfate reduction in sediments and thus determine the abundance, spatial distribution, structure, and activity of diazotroph communities in marginal sea sediments.
SRB have been found to be the principal contributors to the accumulation and persistence of environmental organic Hg in coastal marine sediments (Sunderland et al., ; Schaefer et al., ; Parks et al., ). It is also known that sulfide produced by SRB facilitates precipitation and thus accumulation and persistence of Hg in anoxic marine sediments (Baeyens et al., ). These ecophysiological properties of SRB may explain the correlation between nifH gene abundance and Hg content in the studied Bohai Sea sediments. Hg is one of the key environmental contaminants in many estuaries and coastal bays of the Bohai Sea (Wang et al., ; Gao et al., ), where eutrophication is common, especially in summer times (Wang et al., ). Our data suggest that the accumulation and persistence of Hg in the sediments of the Bohai Sea may be attributable, at least partially, to the activity of nifH-harboring SRB.
The in situ bottom water temperature was putatively identified as a key environmental factor controlling the community structure and spatial distribution of the sediment putative diazotrophs in the Bohai Sea (Supplementary Figure S5). This finding agrees with the previously reported relationship of sediment in situ temperature with the community structure and spatial distribution of nifH-harboring microbial assemblages in the nSCS sediments (Dang et al., ). Our current study supports the importance of temperature as a key environmental factor controlling universally the distribution of the diazotroph communities in marine water columns and sediments; also because temperature has previously been implicated as the most important environmental factor in the spatial distribution of seawater N2-fixing cyanobacteria in the oceans (Stal, ). Temperature influences not only nitrogenase activity but also O2 solubility in seawater and sediment porewater. Enhanced respiration is usually employed by diazotrophs under elevated O2 concentrations (such as under low temperature conditions) for providing energy to replace O2-damaged nitrogenase by de novo synthesis and/or for O2 consumption to maintain anaerobic N2-fixing activity (Stal, ; Großkopf and Laroche, ; Bandyopadhyay et al., ; Brauer et al., ). These biochemical and physiological mechanisms explain well the importance of environmental temperature on the community structure, distribution, and activity of marine diazotrophic microbiota.
The bottom water chlorophyll a content (or the covarying bottom water turbidity) was putatively identified as another key environmental factor controlling the community structure and spatial distribution of the putative diazotrophs in Bohai Sea sediments (Supplementary Figure S5). Chlorophyll a is related to the biomass and primary production of phytoplankton and both chlorophyll a and turbidity may indicate the potential of organic matter export from water column to marine sediments (Sobczak et al., ; Volkman and Tanoue, ). The positive correlation between nifH gene abundance with -Si content in Bohai Sea sediments suggests that diatom productivity affects size and activity of the diazotroph community in sediments (Yool and Tyrrell, ; Wei et al., ). The phylogenetic analysis of deduced NifH sequences (Supplementary Figure S4) suggests that most of the nifH-harboring microbes in Bohai Sea sediments are heterotrophic bacteria. Respiration of bioavailable organic carbon is the major process for heterotrophic diazotrophs to acquire the metabolic energy to fix N2 (Riemann et al., ; Moisander et al., ), which is highly energy demanding (Shanmugam et al., ; Brill, ). Therefore, organic carbon supplies, especially those from metabolizable phytoplankton and benthic diatom products, may constitute the major energy sources to fuel the obligately or facultatively anaerobic diazotrophic microbiota in sediments of coastal seas. In addition, obligately anaerobic bacteria such as SRB and, in particular, the process of N2 fixation, are highly sensitive to oxygen stress (Dixon and Kahn, ; Riemann et al., ; Zhou et al., ). Benthic respiration on organic matter consumes oxygen, which may thus help maintaining the hypoxic to anoxic state of the sediment environment and facilitate N2 fixation by anaerobic diazotrophs. In line with this, sediment porewater Eh, an indicator of the environmental redox status, was found to correlate significantly with the community structure and biogeographical distribution of nifH-harboring microbes in Bohai Sea sediments (Supplementary Figure S5). To maintain the hypoxic to anoxic condition may also be necessary for N2 supply in marine sediments: for instance, the heterotrophic N2 fixation process was linked directly to the N2 production processes by anammox and denitrifying bacteria in marine OMZ waters off the Peru coast (Loescher et al., ). A similar coupling of N2 fixation and N2 production microbial processes may also exist in marine sediments. Indeed, the Bohai Sea sediments harbor diverse anammox bacteria including a new species, Candidatus Scalindua pacifica (Dang et al., ). Therefore, the in situ physical and geochemical condition, rather than the geographical location (Figures 3, 4, Supplementary Figure S5), may play a key role in determining the abundance, community structure, niche availability, biogeographical distribution, and activity of the diazotrophic microbiota in the sediment environments of the Bohai Sea.
In order to examine the general biogeographical and ecological characteristics of diazotrophic communities in sediments of the western Pacific Ocean, we performed a comparative analysis of nifH-harboring microbial communities in Bohai Sea and nSCS sediments (Dang et al., ). The results of community classification analyses indicate that the estuarine and nearshore benthic environments harbor distinctly different putative diazotrophic assemblages from those of the offshore and deep-sea benthic environments (Figure 5). Furthermore, the consistency of the unweighted (only taking into account the NifH OTU composition for the analysis; Figure 5A) and weighted (taking into account both the composition of the NifH OTUs and their relative abundances for the analysis; Figure 5B) community classification results suggests that the difference in composition (Figure 5A) may be the key contributor to the difference between sediment diazotroph communities in the western Pacific Ocean. Although the A3 and E501 sampling stations in the nSCS are very distant and segregated from the Bohai Sea, the sediment diazotrophic assemblages of these two stations are separated from the remaining nSCS stations and clustered with the sampling stations of the Bohai Sea (Figure 5), indicating that the in situ physical and geochemical condition may play a more important role than geographical contiguity in determining the community similarity of diazotrophic microbiota in sediments. In line with this, water depth, sediment porewater DO and in situ temperature were identified as the key environmental factors tentatively controlling the composition, community structure and biogeographical distribution of the N2-fixing microbiota in the western Pacific Ocean sediments (Figure 6). The identification of sediment porewater DO as a controlling environmental factor further verifies the importance of the anaerobic heterotrophs in N2 fixation in marginal sea sediments, while the importance of water depth and in situ temperature suggests niche specialization and segregation between the shallow-water (estuarine and nearshore) and the deep-water (offshore and deep-sea) diazotrophic sediment communities. Our current investigation tentatively identified important general ecological, biogeochemical, and biogeographical characteristics and key environmental factors that influence the sediment diazotrophic microbiota and their potential activity in marginal seas.
Statements
Author contributions
HD conceived and designed the experiments; HZ performed the experiments and analyzed the data; HD, HZ, and MK wrote the paper.
Acknowledgments
The authors thank Prof. Zhinan Zhang from the Ocean University of China for providing the sediment samples and certain environmental parameter measurements, the crew and scientists onboard R/V “Dong Fang Hong 2” for assisting with sample collection and sharing in situ environmental parameter measurements, and Chao Cheng, Guo Chen, Lei Yang, Qiang Li, and Zhiguo Chen for their assistance in the project. This work was supported by the SOA grant GASI-03-01-02-05, China MOST 973 program grant 2013CB955700, China National Key Research and Development Plan grant 2016YFA0601303, NSFC grants 91328209, 91428308, 41076091, and 91028011, CNOOC grants CNOOC-KJ125FZDXM00TJ001-2014 and CNOOC-KJ125FZDXM00ZJ001-2014, and US-NSF grant MCB-1202648.
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.01111
References
1
AltschulS. F.MaddenT. L.SchäfferA. A.ZhangJ.ZhangZ.MillerW.et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res.25, 3389–3402. 10.1093/nar/25.17.3389
2
BaeyensW.MeulemanC.MuhayaB.LeermakersM. (1998). Behaviour and speciation of mercury in the Scheldt estuary (water, sediments and benthic organisms). Hydrobiologia366, 63–79. 10.1023/A:1003124310848
3
BandyopadhyayA.ElvitigalaT.LibertonM.PakrasiH. B. (2013). Variations in the rhythms of respiration and nitrogen fixation in members of the unicellular diazotrophic cyanobacterial genus Cyanothece. Plant Physiol.161, 1334–1346. 10.1104/pp.112.208231
4
BartonL. L.FauqueG. D. (2009). Biochemistry, physiology and biotechnology of sulfate-reducing bacteria. Adv. Appl. Microbiol.68, 41–98. 10.1016/S0065-2164(09)01202-7
5
BergmanB.SandhG.LinS.LarssonJ.CarpenterE. J. (2013). Trichodesmium - a widespread marine cyanobacterium with unusual nitrogen fixation properties. FEMS Microbiol. Rev.37, 286–302. 10.1111/j.1574-6976.2012.00352.x
6
BerticsV. J.LöscherC. R.SalonenI.DaleA. W.GierJ.SchmitzR. A.et al. (2013). Occurrence of benthic microbial nitrogen fixation coupled to sulfate reduction in the seasonally hypoxic Eckernförde Bay, Baltic Sea. Biogeosciences10, 1243–1258. 10.5194/bg-10-1243-2013
7
BerticsV. J.SohmJ. A.TreudeT.ChowC. T.CaponeD. G.FuhrmanJ. A.et al. (2010). Burrowing deeper into benthic nitrogen cycling: the impact of bioturbation on nitrogen fixation coupled to sulfate reduction. Mar. Ecol. Prog. Ser.409, 1–15. 10.3354/meps08639
8
BirdC.WymanM. (2013). Transcriptionally active heterotrophic diazotrophs are widespread in the upper water column of the Arabian Sea. FEMS Microbiol. Ecol.84, 189–200. 10.1111/1574-6941.12049
9
BoeschD. F. (2002). Challenges and opportunities for science in reducing nutrient over-enrichment of coastal ecosystems. Estuaries25, 886–900. 10.1007/BF02804914
10
BonnetS.DekaezemackerJ.Turk-KuboK. A.MoutinT.HamersleyR. M.GrossoO.et al. (2013). Aphotic N2 fixation in the eastern tropical south Pacific Ocean. PLoS ONE8:e81265. 10.1371/journal.pone.0081265
11
BoydP. W.EllwoodM. J. (2010). The biogeochemical cycle of iron in the ocean. Nat. Geosci.3, 675–682. 10.1038/ngeo964
12
BrauerV. S.StompM.RossoC.van BeusekomS. A.EmmerichB.StalL. J.et al. (2013). Low temperature delays timing and enhances the cost of nitrogen fixation in the unicellular cyanobacterium Cyanothece. ISME J.7, 2105–2115. 10.1038/ismej.2013.103
13
BrillW. J. (1980). Biochemical genetics of nitrogen fixation. Microbiol. Rev.44, 449–467.
14
BurnsJ. A.ZehrJ. P.CaponeD. G. (2002). Nitrogen-fixing phylotypes of Chesapeake Bay and Neuse River estuary sediments. Microb. Ecol.44, 336–343. 10.1007/s00248-002-1000-9
15
DangH.JiaoN. (2014). Perspectives on the microbial carbon pump with special reference to microbial respiration and ecosystem efficiency in large estuarine systems. Biogeosciences11, 3887–3898. 10.5194/bg-11-3887-2014
16
DangH.LovellC. R. (2000). Bacterial primary colonization and early succession on surfaces in marine waters as determined by amplified rRNA gene restriction analysis and sequence analysis of 16S rRNA genes. Appl. Environ. Microbiol.66, 467–475. 10.1128/AEM.66.2.467-475.2000
17
DangH.LovellC. R. (2016). Microbial surface colonization and biofilm development in marine environments. Microbiol. Mol. Biol. Rev.80, 91–138. 10.1128/MMBR.00037-15
18
DangH. Y.LiT. G.ChenM. N.HuangG. Q. (2008). Cross-ocean distribution of Rhodobacterales bacteria as primary surface colonizers in temperate coastal marine waters. Appl. Environ. Microbiol.74, 52–60. 10.1128/AEM.01400-07
19
DangH. Y.LuanX. W.ZhaoJ. Y.LiJ. (2009a). 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. 10.1128/AEM.02556-08
20
DangH. Y.WangC. Y.LiJ.LiT. G.TianF.JinW.et al. (2009b). Diversity and distribution of sediment nirS-encoding bacterial assemblages in response to environmental gradients in the eutrophied Jiaozhou Bay, China. Microb. Ecol.58, 161–169. 10.1007/s00248-008-9469-5
21
DangH. Y.YangJ. Y.LiJ.LuanX. W.ZhangY. B.GuG. Z.et al. (2013a). Environment-dependent distribution of the sediment nifH-harboring microbiota in the Northern South China Sea. Appl. Environ. Microbiol.79, 121–132. 10.1128/AEM.01889-12
22
DangH. Y.ZhouH. X.ZhangZ. N.YuZ. S.HuaE.LiuX. S.et al. (2013b). Molecular detection of Candidatus Scalindua pacifica and environmental responses of sediment anammox bacterial community in the Bohai Sea, China. PLoS ONE8:e61330. 10.1371/journal.pone.0061330
23
DekaezemackerJ.BonnetS.GrossoO.MoutinT.BressacM.CaponeD. G. (2013). Evidence of active dinitrogen fixation in surface waters of the eastern tropical South Pacific during El Niño and La Niña events and evaluation of its potential nutrient controls. Global Biogeochem. Cycles27, 768–779. 10.1002/gbc.20063
24
DekasA. E.ChadwickG. L.BowlesM. W.JoyeS. B.OrphanV. J. (2014). Spatial distribution of nitrogen fixation in methane seep sediment and the role of the ANME archaea. Environ. Microbiol.16, 3012–3029. 10.1111/1462-2920.12247
25
DekasA. E.ConnonS. A.ChadwickG. L.Trembath-ReichertE.OrphanV. J. (2016). Activity and interactions of methane seep microorganisms assessed by parallel transcription and FISH-NanoSIMS analyses. ISME J.10, 678–692. 10.1038/ismej.2015.145
26
DekasA. E.PoretskyR. S.OrphanV. J. (2009). Deep-sea archaea fix and share nitrogen in methane-consuming microbial consortia. Science326, 422–426. 10.1126/science.1178223
27
DeutschC.SarmientoJ. L.SigmanD. M.GruberN.DunneJ. P. (2007). Spatial coupling of nitrogen inputs and losses in the ocean. Nature445, 163–167. 10.1038/nature05392
28
DixonR.KahnD. (2004). Genetic regulation of biological nitrogen fixation. Nat. Rev. Microbiol.2, 621–631. 10.1038/nrmicro954
29
DoreJ. E.BrumJ. R.TupasL. M.KarlD. M. (2002). Seasonal and interannual variability in sources of nitrogen supporting export in the oligotrophic subtropical North Pacific Ocean. Limnol. Oceanogr.47, 1595–1607. 10.4319/lo.2002.47.6.1595
30
FarnelidH.AnderssonA. F.BertilssonS.Al-SoudW. A.HansenL. H.SørensenS.et al. (2011). Nitrogenase gene amplicons from global marine surface waters are dominated by genes of non-cyanobacteria. PLoS ONE6:e19223. 10.1371/journal.pone.0019223
31
FarnelidH.Bentzon-TiliaM.AnderssonA. F.BertilssonS.JostG.LabrenzM.et al. (2013). Active nitrogen-fixing heterotrophic bacteria at and below the chemocline of the central Baltic Sea. ISME J.7, 1413–1423. 10.1038/ismej.2013.26
32
FelsensteinJ. (1989). PHYLIP - phylogeny inference package (Version 3.2). Cladistics5, 164–166.
33
FernandezC.FaríasL.UlloaO. (2011). Nitrogen fixation in denitrified marine waters. PLoS ONE6:e20539. 10.1371/journal.pone.0020539
34
FulweilerR. W. (2009). Fantastic fixers. Science326, 377–378. 10.1126/science.1181129
35
GaoX.ZhouF.ChenC. T. A. (2014). Pollution status of the Bohai Sea: an overview of the environmental quality assessment related trace metals. Environ. Int.62, 12–30. 10.1016/j.envint.2013.09.019
36
GroßkopfT.LarocheJ. (2012). Direct and indirect costs of dinitrogen fixation in Crocosphaera watsonii WH8501 and possible implications for the nitrogen cycle. Front. Microbiol.3:236. 10.3389/fmicb.2012.00236
37
GroßkopfT.MohrW.BaustianT.SchunckH.GillD.KuypersM. M.et al. (2012). Doubling of marine dinitrogen-fixation rates based on direct measurements. Nature488, 361–364. 10.1038/nature11338
38
HamadyM.LozuponeC.KnightR. (2010). Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J.4, 17–27. 10.1038/ismej.2009.97
39
HerbertR. A. (1999). Nitrogen cycling in coastal marine ecosystems. FEMS Microbiol. Rev.23, 563–590. 10.1111/j.1574-6976.1999.tb00414.x
40
HomokyW. B.JohnS. G.ConwayT. M.MillsR. A. (2013). Distinct iron isotopic signatures and supply from marine sediment dissolution. Nat. Commun.4, 2143. 10.1038/ncomms3143
41
JacqV.RidameC.L'HelguenS.KaczmarF.SaliotA. (2014). Response of the unicellular diazotrophic cyanobacterium Crocosphaera watsonii to iron limitation. PLoS ONE9:e92655. 10.1371/journal.pone.0086749
42
JayakumarA.Al-RshaidatM. M.WardB. B.MulhollandM. R. (2012). Diversity, distribution, and expression of diazotroph nifH genes in oxygen-deficient waters of the Arabian Sea. FEMS Microbiol. Ecol.82, 597–606. 10.1111/j.1574-6941.2012.01430.x
43
KarlD.LetelierR.TupasL.DioreJ.ChristianJ.HemelD. (1997). The role of nitrogen fixation in biogeochemical cycling in the subtropical North Pacific Ocean. Nature388, 533–538. 10.1038/41474
44
KarlD.MichaelsA.BergmanB.CaponeD.CarpenterE.LetelierR.et al. (2002). Dinitrogen fixation in the world's oceans. Biogeochemistry57, 47–98. 10.1023/A:1015798105851
45
KarlD. M. (2014). Microbially mediated transformations of phosphorus in the sea: new views of an old cycle. Annu. Rev. Mar. Sci.6, 279–337. 10.1146/annurev-marine-010213-135046
46
KarlD. M.LetelierR. M. (2008). Nitrogen fixation-enhanced carbon sequestration in low nitrate, low chlorophyll seascapes. Mar. Ecol. Prog. Ser.364, 257–268. 10.3354/meps07547
47
KnappA. N. (2012). The sensitivity of marine N2 fixation to dissolved inorganic nitrogen. Front. Microbiol.3:374. 10.3389/fmicb.2012.00374
48
LarkinM. A.BlackshieldsG.BrownN. P.ChennaR.McGettiganP. A.McWilliamH.et al. (2007). Clustal W and Clustal X version 2.0. Bioinformatics23, 2947–2948. 10.1093/bioinformatics/btm404
49
LepšJ.ŠmilauerP. (2003). Multivariate Analysis of Ecological Data Using CANOCO. New York, NY: Cambridge University Press.
50
LoescherC. R.GroßkopfT.DesaiF. D.GillD.SchunckH.CrootP. L.et al. (2014). Facets of diazotrophy in the oxygen minimum zone waters off Peru. ISME J.8, 2180–2192. 10.1038/ismej.2014.71
51
LovellC. R.DavisD. A. (2012). Specificity of salt marsh diazotrophs for vegetation zones and plant hosts: results from a North American marsh. Front. Microbiol.3:84. 10.3389/fmicb.2012.00084
52
LovellC. R.FriezM. J.LongshoreJ. W.BagwellC. E. (2001). Recovery and phylogenetic analysis of nifH sequences from diazotrophic bacteria associated with dead aboveground biomass of Spartina alterniflora. Appl. Environ. Microbiol.67, 5308–5314. 10.1128/AEM.67.11.5308-5314.2001
53
MahaffeyC.MichaelsA. F.CaponeD. G. (2005). The conundrum of marine N2 fixation. Am. J. Sci.305, 546–595. 10.2475/ajs.305.6-8.546
54
McGlatheryK. J.Risgaard-PetersenN.ChristensenP. B. (1998). Temporal and spatial variation in nitrogen fixation activity in the eelgrass Zostera marina rhizosphere. Mar. Ecol. Prog. Ser.168, 245–258. 10.3354/meps168245
55
MehtaM. P.ButterfieldD. A.BarossJ. A. (2003). Phylogenetic diversity of nitrogenase (nifH) genes in deep-sea and hydrothermal vent environments of the Juan de Fuca Ridge. Appl. Environ. Microbiol.69, 960–970. 10.1128/AEM.69.2.960-970.2003
56
MillsM. M.RidameC.DaveyM.La RocheJ.GeiderR. J. (2004). Iron and phosphorus co-limit nitrogen fixation in the eastern tropical North Atlantic. Nature429, 292–294. 10.1038/nature02550
57
MiyazakiJ.HigaR.TokiT.AshiJ.TsunogaiU.NunouraT.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. 10.1128/AEM.01184-09
58
MoisanderP. H.SerrosT.PaerlR. W.BeinartR. A.ZehrJ. P. (2014). Gammaproteobacterial diazotrophs and nifH gene expression in surface waters of the South Pacific Ocean. ISME J.8, 1962–1973. 10.1038/ismej.2014.49
59
MontoyaJ. P.HollC. M.ZehrJ. P.HansenA.VillarealT. A.CaponeD. G. (2004). High rates of N2 fixation by unicellular diazotrophs in the oligotrophic Pacific Ocean. Nature430, 1027–1032. 10.1038/nature02824
60
MooreC. M.MillsM. M.AchterbergE. P.GeiderR. J.LaRocheJ.LucasM. I.et al. (2009). Large-scale distribution of Atlantic nitrogen fixation controlled by iron availability. Nat. Geosci.2, 867–871. 10.1038/ngeo667
61
MooreC. M.MillsM. M.ArrigoK. R.Berman-FrankI.BoppL.BoydP. W.et al. (2013). Processes and patterns of oceanic nutrient limitation. Nat. Geosci.6, 701–710. 10.1038/ngeo1765
62
MullinsT. D.BritschgiT. B.KrestR. L.GiovannoniS. J. (1995). Genetic comparisons reveal the same unknown bacterial lineages in Atlantic and Pacific bacterioplankton communities. Limnol. Oceanogr.40, 148–158. 10.4319/lo.1995.40.1.0148
63
PEMSEA and BSEMP (2005). Bohai Sea Environmental Risk Assessment. PEMSEA Technical Report No. 12. 114. Partnerships in Environmental Management for the Seas of East Asia, Bohai Sea Environmental Management Project of the People's Republic of China.
64
PaerlH. W. (1990). Physiological ecology and regulation of N2 fixation in natural waters. Adv. Microb. Ecol.11, 305–344. 10.1007/978-1-4684-7612-5_8
65
ParksJ. M.JohsA.PodarM.BridouR.HurtR. A.Jr.SmithS. D.et al. (2013). The genetic basis for bacterial mercury methylation. Science339, 1332–1335. 10.1126/science.1230667
66
PernthalerA.DekasA. E.BrownC. T.GoffrediS. K.EmbayeT.OrphanV. 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. 10.1073/pnas.0711303105
67
RahavE.Bar-ZeevE.OhayonS.ElifantzH.BelkinN.HerutB.et al. (2013). Dinitrogen fixation in aphotic oxygenated marine environments. Front. Microbiol.4:227. 10.3389/fmicb.2013.00227
68
RiemannL.FarnelidH.StewardG. F. (2010). Nitrogenase genes in non- cyanobacterial plankton: prevalence, diversity and regulation in marine waters. Aquat. Microb. Ecol.61, 235–247. 10.3354/ame01431
69
Sañudo-WilhelmyS. A.KustkaA. B.GoblerC. J.HutchinsD. A.YangM.LwizaK.et al. (2001). Phosphorus limitation of nitrogen fixation by Trichodesmium in the central Atlantic Ocean. Nature411, 66–69. 10.1038/35075041
70
SchaeferJ. K.RocksS. S.ZhengW.LiangL.GuB.MorelF. M. M. (2011). Active transport, substrate specificity, and methylation of Hg(II) in anaerobic bacteria. Proc. Natl. Acad. Sci. U.S.A.108, 8714–8719. 10.1073/pnas.1105781108
71
SchlossP. D.HandelsmanJ. (2005). Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl. Environ. Microbiol.71, 1501–1506. 10.1128/AEM.71.3.1501-1506.2005
72
ShanmugamK. T.O'GaraF.AndersenK.ValentineR. C. (1978). Biological nitrogen fixation. Annu. Rev. Plant Physiol.29, 263–276. 10.1146/annurev.pp.29.060178.001403
73
ShaoS. D.LuanX. W.DangH. Y.ZhouH. X.ZhaoY. K.LiuH. T.et al. (2014). Deep-sea methane seep sediments in the Okhotsk Sea sustain diverse and abundant anammox bacteria. FEMS Microbiol. Ecol.87, 503–516. 10.1111/1574-6941.12241
74
SmithV. H. (2003). Eutrophication of freshwater and coastal marine ecosystems: a global problem. Environ. Sci. Pollut. Res. Int.10, 126–139. 10.1065/espr2002.12.142
75
SobczakW. V.CloernJ. E.JassbyA. D.Müller-SolgerA. B. (2002). Bioavailability of organic matter in a highly disturbed estuary: the role of detrital and algal resources. Proc. Natl. Acad. Sci. U.S.A.99, 8101–8105. 10.1073/pnas.122614399
76
SohmJ. A.WebbE. A.CaponeD. G. (2011). Emerging patterns of marine nitrogen fixation. Nat. Rev. Microbiol.9, 499–508. 10.1038/nrmicro2594
77
StalL. J. (2009). Is the distribution of nitrogen-fixing cyanobacteria in the oceans related to temperature?Environ. Microbiol. 11, 1632–1645. 10.1111/j.1758-2229.2009.00016.x
78
StalL. J.SeverinI.BolhuisH. (2010). The ecology of nitrogen fixation in cyanobacterial mats. Adv. Exp. Med. Biol.675, 31–45. 10.1007/978-1-4419-1528-3_3
79
SteppeT. F.PaerlH. W. (2002). Potential N2 fixation by sulfate-reducing bacteria in a marine intertidal microbial mat. Aquat. Microb. Ecol.28, 1–12. 10.3354/ame028001
80
SteppeT. F.PaerlH. W. (2005). Nitrogenase activity and nifH expression in a marine intertidal microbial mat. Microb. Ecol.49, 315–324. 10.1007/s00248-004-0245-x
81
StreetJ. H.PaytanA. (2005). Iron, phytoplankton growth, and the carbon cycle. Met. Ions. Biol. Syst. 43, 153–193. 10.1201/9780824751999.ch7
82
SunderlandE. M.GobasF. A. P. C.HeyesA.BranfireunB. A.BayerA. K.CranstonR. E.et al. (2004). Speciation and bioavailability of mercury in well-mixed estuarine sediments. Mar. Chem.90, 91–105. 10.1016/j.marchem.2004.02.021
83
ter BraakC. J. F.ŠmilauerP. (2002). CANOCO Reference Manual and CanoDraw for Windows User's Guide: Software for Canonical Community Ordination (Version 4.5).Ithaca, NY: Microcomputer Power.
84
TrimmerM.NichollsJ. C. (2009). Production of nitrogen gas via anammox and denitrification in intact sediment cores along a continental shelf to slope transect in the North Atlantic. Limnol. Oceanogr.54, 577–589. 10.4319/lo.2009.54.2.0577
85
VolkmanJ. K.TanoueE. (2002). Chemical and biological studies of particulate organic matter in the ocean. J. Oceanogr.58, 265–279. 10.1023/A:1015809708632
86
VossM.BangeH. W.DippnerJ. W.MiddelburgJ. J.MontoyaJ. P.WardB. (2013). The marine nitrogen cycle: recent discoveries, uncertainties and the potential relevance of climate change. Philos. Trans. R. Soc. B.368:20130121. 10.1098/rstb.2013.0121
87
WangJ.ChenS.XiaT. (2010). Environmental risk assessment of heavy metals in Bohai Sea, North China. Procedia Environ. Sci.2, 1632–1642. 10.1016/j.proenv.2010.10.174
88
WangX. L.CuiZ. G.GuoQ.HanX. R.WangJ. T. (2009). Distribution of nutrients and eutrophication assessment in the Bohai Sea of China. Chin. J. Oceanol. Limnol.27, 177–183. 10.1007/s00343-009-0177-x
89
WeiH.SunJ.MollA.ZhaoL. (2004). Phytoplankton dynamics in the Bohai Sea-observations and modelling. J. Mar. Syst.44, 233–251. 10.1016/j.jmarsys.2003.09.012
90
YoolA.TyrrellT. (2003). Role of diatoms in regulating the ocean's silicon cycle. Global Biogeochem. Cycles17, 1103–1121. 10.1029/2002GB002018
91
ZehrJ. P. (2011). Nitrogen fixation by marine cyanobacteria. Trends Microbiol.19, 162–173. 10.1016/j.tim.2010.12.004
92
ZehrJ. P.BenchS. R.CarterB. J.HewsonI.NiaziF.ShiT.et al. (2008). Globally distributed uncultivated oceanic N2-fixing cyanobacteria lack oxygenic photosystem II. Science322, 1110–1112. 10.1126/science.1165340
93
ZehrJ. P.KudelaR. M. (2011). Nitrogen cycle of the open ocean: from genes to ecosystems. Annu. Rev. Mar. Sci.3, 197–225. 10.1146/annurev-marine-120709-142819
94
ZehrJ. P.MellonM.BraunS.LitakerW. (1995). Diversity of heterotrophic nitrogen fixation genes in a marine cyanobacterial mat. Appl. Environ. Microbiol.61, 2527–2532.
95
ZehrJ. P.WaterburyJ. B.TurnerP. J.MontoyaJ. P.OmoregieE.StewardG. F.et al. (2001). Unicellular cyanobacteria fix N2 in the subtropical Pacific Ocean. Nature412, 635–638. 10.1038/news010809-11
96
ZhouJ.HeQ.HemmeC. L.MukhopadhyayA.HilleslandK.ZhouA.et al. (2011). How sulphate-reducing microorganisms cope with stress: lessons from systems biology. Nat. Rev. Microbiol.9, 452–466. 10.1038/nrmicro2575
Summary
Keywords
coastal sediments, ecophysiological specialization, estuary, marginal sea, niche differentiation, nifH, nitrogen fixation, sulfate-reducing bacteria
Citation
Zhou H, Dang H and Klotz MG (2016) Environmental Conditions Outweigh Geographical Contiguity in Determining the Similarity of nifH-Harboring Microbial Communities in Sediments of Two Disconnected Marginal Seas. Front. Microbiol. 7:1111. doi: 10.3389/fmicb.2016.01111
Received
09 March 2016
Accepted
04 July 2016
Published
20 July 2016
Volume
7 - 2016
Edited by
Dennis A. Bazylinski, University of Nevada at Las Vegas, USA
Reviewed by
Zhe-Xue Quan, Fudan University, China; Zhanfei Liu, The University of Texas at Austin, USA
Updates
Copyright
© 2016 Zhou, Dang and Klotz.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongyue Dang danghy@xmu.edu.cn
This article was submitted to Aquatic Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.