Skip to main content

Microbial river-to-sea continuum: gradients in benthic and planktonic diversity, osmoregulation and nutrient cycling



Coastal aquatic ecosystems include chemically distinct, but highly interconnected environments. Across a freshwater-to-marine transect, aquatic communities are exposed to large variations in salinity and nutrient availability as tidal cycles create periodic fluctuations in local conditions. These factors are predicted to strongly influence the resident microbial community structure and functioning, and alter the structure of aquatic food webs and biogeochemical cycles. Nevertheless, little is known about the spatial distribution of metabolic properties across salinity gradients, and no study has simultaneously surveyed the sediment and water environments. Here, we determined patterns and drivers of benthic and planktonic prokaryotic and microeukaryotic community assembly across a river and tidal lagoon system by collecting sediments and planktonic biomass at nine shallow subtidal sites in the summer. Genomic and transcriptomic analyses, alongside a suite of complementary geochemical data, were used to determine patterns in the distribution of taxa, mechanisms of salt tolerance, and nutrient cycling.


Taxonomic and metabolic profiles related to salt tolerance and nutrient cycling of the aquatic microbiome were found to decrease in similarity with increasing salinity, and distinct trends in diversity were observed between the water column and sediment. Non-saline and saline communities adopted divergent strategies for osmoregulation, with an increase in osmoregulation-related transcript expression as salinity increased in the water column due to lineage-specific adaptations to salt tolerance. Results indicated a transition from phosphate limitation in freshwater habitats to nutrient-rich conditions in the brackish zone, where distinct carbon, nitrogen and sulfur cycling processes dominated. Phosphorus acquisition-related activity was highest in the freshwater zone, along with dissimilatory nitrate reduction to ammonium in freshwater sediment. Activity associated with denitrification, sulfur metabolism and photosynthesis were instead highest in the brackish zone, where photosynthesis was dominated by distinct microeukaryotes in water (Cryptophyta) and sediment (diatoms). Despite microeukaryotes and archaea being rare relative to bacteria, results indicate that they contributed more to photosynthesis and ammonia oxidation, respectively.


Our study demonstrates clear freshwater–saline and sediment–water ecosystem boundaries in an interconnected coastal aquatic system and provides a framework for understanding the relative importance of salinity, planktonic-versus-benthic habitats and nutrient availability in shaping aquatic microbial metabolic processes, particularly in tidal lagoon systems.

Video abstract


Across connected riverine, estuarine and marine environments, freshwater flow and tidal influx helps to redistribute nutrients and establish a strong salinity gradient [1,2,3]. Such gradients directly influence the resident community composition and distribution, resulting in lower benthic macrofaunal species richness within the horohalinicum (where salinities are between 5 and 8 [4]) and higher microbial richness [5,6,7,8]. Global surveys of microbial diversity have demonstrated that both salinity and environment type (such as water or sediment) are major factors influencing bacterial [9] and archaeal community composition [10]. Besides nutrients, river discharge and phytoplankton community composition are also known to drive changes in planktonic estuarine communities [11, 12]. Previous studies of planktonic prokaryotic communities across freshwater-to-marine transects have highlighted significant differences in gene abundances related to glycolysis, respiration, catabolic pathways, osmolyte and metal transport, and the biosynthesis of quinones and isoprenoids between low and high salinities, in line with the shifts in taxonomic groups [13, 14]. However, little is known about the collective microbial ecology of bacteria, archaea and microeukaryotes in sediment and water across river-to-sea salinity gradients; how these factors (water, sediment, salinity) affect microbial gene expression across this environmental continuum also remains largely unexplored.

A key factor shaping microbial community distributions across salinity gradients is lineage-specific adaptations for salt tolerance [6, 8,9,10, 12, 15, 16]. Bacteria adopt two distinct strategies to cope with osmolarity changes in their environment: (1) a salt-in strategy which promotes K+ accumulation and (2) the synthesis or uptake of organic compatible solutes which function as osmoprotectants [17]. Various K+ uptake systems have been reported to date, for example, both trkA and kup systems in hyperosmotic environments, and the kdp system for low external K+ concentrations [18,19,20,21]. Bacteria can also employ various Na+/H+ antiporters and Na+ translocating dehydrogenases to export Na+ extracellularly in highly saline environments [13, 15]. Taxa and strategies employed are therefore predicted to change across coastal salinity gradients. Although, certain lineages, for instance, members of the order Pelagibacterales (SAR11), have been found at different salinities, suggesting ecological adaptations within microbial lineages to salt [16, 22].

In addition to salinity, dramatic differences in nutrient availability and types are typically observed across freshwater to marine systems and have been predicted to strongly influence community composition and function [1,2,3]. These systems mix nitrogen-rich freshwater with sulfate and phosphate-rich marine water [23, 24]. Sulfate-rich saline water promotes the formation of ferrosulfides, resulting in less phosphate sorption in sediments and higher phosphate concentrations in saline water [23, 25, 26]. Riverine systems, however, are often phosphate-limited due to greater phosphate sorption in non-saline sediments [27]; it is therefore likely that freshwater microbial communities rely on high-affinity phosphate uptake mechanisms. Estuarine conditions are comparatively rich in both nitrogen and phosphorus, which support the growth of primary producers [28], and in turn, heterotrophs [29, 30]. Benthic photosynthesis has been shown to promote nutrient-regulating biogeochemical processes in estuaries, such as coupled nitrification–denitrification [31]. However, the relative contribution of each of the taxa within the microbiome toward biogeochemical cycling across river–estuary–marine transects remains understudied.

To determine microbial taxonomic and functional diversity, and metabolic hotspots, across a river-to-sea transect, we sampled benthic sediment (hereafter referred to as “benthic” samples) and overlying water (“planktonic” samples) along a 5-km stretch connecting the Waiwera river, estuary and beach (Auckland, New Zealand). The Waiwera river predominantly drains rural pastoral catchments [32], and is therefore impacted by anthropogenic activities, as is typical for estuaries globally [33, 34]. The river and estuary have a water depth of approximately 0.2–1 m, and the estuary is classified as a permanently open tidal lagoon estuary [35] and is subjected to daily tidal influxes from the Hauraki Gulf (Pacific Ocean). Tidal lagoon systems, also known as bar-built coastal lagoons, or barrier enclosed lagoons, can be found throughout the world (e.g. Laguna Madre in Texas, USA [36]; Mosquito Lagoon in Florida, USA [37]; and Paravur estuary in Kerala, India [38]) and occupy about 13% of the global coastline [39]. They are also the most common estuary systems in New Zealand [40] and the UK [41]. A natural continuous salinity gradient occurs along the Waiwera river and estuary, potentially exerting a strong influence on the assembly of microbial communities. We hypothesised that non-saline, saline, water and sediment communities would be strongly demarcated by distinct prokaryotic and eukaryotic compositions, and utilise distinct osmoregulatory and nutrient acquisition mechanisms, reflecting, for example, nutrient limitations in the freshwater (P) and marine (N) end-member environments. Combining metagenomic, metatranscriptomic and geochemical data, we identified the impact of salinity and nutrient availability on prokaryotic and microeukaryotic community diversity and on critical physiological and ecological processes, such as osmoregulation, primary production and nutrient cycling.

Results and discussion

Community partitioning based on salinity and sediment–water habitat preferences

Contribution of salinity versus other environmental parameters in shaping microbial communities

Nutrient and salinity gradients are major features of estuarine systems, leading to higher estuarine productivity [1,2,3]. Accordingly, we observed a general decrease in nitrate concentrations along the Waiwera river–estuary–sea transect, higher organic carbon concentrations in brackish sediments, and an increase in dissolved reactive phosphorus and sulfate with salinity (Fig. 1a, b). We therefore tested the contribution of salinity relative to other geochemical variables in differentiating prokaryotic and microeukaryotic composition along this transect. To achieve this, near full-length small subunit (SSU) rRNA gene sequences were assembled and clustered into 6964 bacterial, 291 microeukaryal and 79 archaeal operational taxonomic units (OTUs, unique at 97% sequence identity). Results showed salinity best explains differences in sediment bacterial and microeukaryal community composition (Spearman’s correlation coefficient ρ = 0.43–0.91) and planktonic bacterial and microeukaryal composition (ρ = 0.78–0.91, Supplementary Table 1). Two subsets of parameters were maximally correlated with dissimilarities among planktonic (salinity, ammonium, and phosphate, ρ = 0.19) and benthic archaea (salinity and nitrate, ρ = 0.90, Supplementary Table 1), suggesting that besides salinity, nutrients play important roles in structuring these aquatic communities, in agreement with previous studies [11, 12, 22]. Overall, our results correlate strongly with research showing salinity is a major global driver of prokaryotic community composition [9, 10].

Fig. 1

Plots showing variation in nutrient concentrations and community composition. a The Waiwera river sampling locations (1 to 9) with water column salinity gradient (0–35) shown using kriging interpolation. b Nutrient concentrations in sediment, sediment pore water and water samples taken in triplicate across the sampling sites. Grey background represents brackish sites (2.0–4.5 km along the transect). Abbreviations: TRP total recoverable phosphorus, TS total sulfur, TN total nitrogen, TOC total organic carbon, DRP dissolved reactive phosphorus, DNPOC dissolved non-purgeable organic carbon. Error bars represent standard errors of means. Asterisks () indicate a significant difference between all water column and sediment samples (Wilcoxon p < 0.05). c Multivariate regression tree of microbial community abundance data associated with the environmental variables. Abbreviations: s salinity, NO2 nitrite, n sample size. Units: nitrite (g/m3). d NMDS ordination of small subunit rRNA gene sequences based on Bray–Curtis dissimilarities. Samples are coloured according to salinity. Environmental vectors were fitted onto the NMDS scores of the microbial community by the R-function envfit (p < 0.05, permutation = 999). e Distance–decay relationships of community dissimilarity (Bray–Curtis index) with increasing geographical distance. The regression lines were fitted to the data points using log generalised linear (GLM, solid line) and LOESS (dotted line) models. Each data point represents a pairwise comparison of samples from the water column (closed circles) or sediment (open circles)

Although salinity strongly differentiated each of the planktonic and benthic communities relative to other geochemical factors, water-vs-sediment environment type was found to be the single most important factor distinguishing microbial communities overall (Fig. 1c). This was due to consistently large differences in benthic and planktonic communities across the gradient, regardless of salinity (Fig. 1d; Supplementary Table 2), which resulted in significantly different benthic and planktonic microbial communities overall based on Bray–Curtis dissimilarities (Benjamini–Hochberg adjusted p < 0.05, Supplementary Table 3; Fig. 1d and Supplementary Figure 1). Results, therefore, support research highlighting the importance of these environment types, along with salinity, in structuring microbial communities [9, 10]. However, few studies, as here, have compared sediment and water communities directly within the same system [42, 43]. We found that only when comparing the effect of end-member salinities (non-saline vs marine) on benthic and planktonic communities were differences in microbial community composition comparable to those between water and sediment at the same salinity (Bray–Curtis dissimilarities 0.97 to 0.99 on average; Supplementary Table 2). Predictably, the combination of salinity and environment type contributed to almost entirely distinct communities (Bray–Curtis dissimilarities of 0.99 to 1.0, or 1.0 on average). Sediment and water environments are typically separated by large differences in redox processes, oxygen, nutrient and terminal electron acceptor availability [44, 45]. This is also evident in our study, where concentrations of ammonium, nitrate and dissolved non-purgeable organic carbon (DNPOC) were significantly higher in the sediment porewater, compared with the overlying water column (Wilcoxon, p < 0.05, Fig. 1b).

Salinity gradients and the non-saline–saline divide in sediment and water microbial communities

Consistent with salt as a major driver of aquatic community composition, reconstructed SSU rRNA and Braken2-estimated species abundance profiles indicated large differences in the composition of bacteria, archaea and microeukaryotes across the salinity gradient within sediment and water environments, although bacteria invariably dominated both the benthic and planktonic communities (Fig. 2a, trends exhibited by key taxa are described in Supplementary Information). An overall decrease in both planktonic and benthic community similarities occurred with increased difference in salinity and geographic distance (Fig. 1e and Supplementary Figure 2a), resulting in a clear separation between non-saline and saline community compositions (Fig. 1d). Results reflect the substantial physiological adaptations required for salt tolerance [14,15,16]. The rate of increasing dissimilarity was far greater among water column communities (analysis of covariance, p < 0.0001), owing to the relatively high similarity between neighbouring communities. Greater spatial dissimilarity overall was observed among communities at benthic sites, including among sites located only 5 m apart (Fig. 1e). This suggests sediment is a more heterogeneous environment compared to the overlying water, and is consistent with previous comparisons of planktonic and benthic microbial communities in river or coastal settings [42, 43]. The strong spatial heterogeneity of benthic communities is likely due to greater stability within sediments, which are less impacted by tidal mixing [1, 46].

Fig. 2

Community distribution and alpha diversity. a Bar plots indicate relative abundances of EMIRGE-reconstructed 16S and 18S rRNA gene sequences. Site salinities are shown in the upper left plots, and transect distances are given in the lower plots. Taxa are coloured at the order level. For each plot, samples from the water column (left) and sediment (right) are further categorised into freshwater, brackish or marine, as indicated in the top left plot. b Boxplots illustrate Shannon indices and relative OTU richness in both sediment and water column environments across the salinity gradient for archaeal, bacterial and microeukaryal communities (n = 36). Colours indicate environment type (freshwater, brackish and marine) based on the salinity gradient. Boxes represent the interquartile range (IQR) between the first and third quartiles and the horizontal line inside the box represents the median. Whiskers represent the lowest and highest values within 1.5 times the IQR from the first and third quartiles, respectively. Individual sample values are shown as dots. Asterisks () indicate a significant difference between all water column and sediment samples (Wilcoxon p < 0.05)

To further compare community differences across the transect, we categorised sediment and overlying water samples as each belonging to three distinct environments (non-saline, 0–0.5; brackish, 0.5–30; marine, 30–35), based on the salinity of porewater or water samples, respectively (Fig. 1a and Supplementary Table 4) [47]. The benthic community composition (based on OTU abundances) was significantly different among non-saline, brackish and marine environments (Benjamini–Hochberg adjusted p < 0.05, Supplementary Table 3). Fewer benthic and planktonic OTUs were shared between freshwater and brackish (6.4–10.9%) than between saline brackish and marine (13.5–17.9%) environments (Supplementary Figure 2b). This is reflected by the planktonic distance–decay curve, which shows greater dissimilarity at the freshwater–brackish transition (peak at ~1.5 km) and higher similarity at the brackish–marine transition (trough at ~2.5 km; Fig. 1e). These results support findings that salinity strongly partitions freshwater and marine planktonic prokaryotic communities [14, 22], and also demonstrate that the relatively small increase in salinity (of 8) spanning the freshwater–saltwater transition represents a powerful ecological barrier for both benthic and planktonic communities.

Influence of salinity on microbial alpha diversity

Microbial community alpha diversity, measured via the Shannon index and relative OTU richness, was significantly higher in sediments across the entire transect (Fig. 2a; Wilcoxon, p < 0.001), likely due to greater nutrient availability (Fig. 1b), availability of terminal electron acceptors, greater diversity in environmental niche space and greater community stability associated with more static environmental conditions [1, 46, 48]. Benthic and planktonic communities exhibited opposing trends in richness and Shannon diversity across the salinity gradient (Fig. 2b). In contrast to previous planktonic studies that showed no clear trend in diversity [6, 7] or a reduction in diversity with increasing salinity [11, 12], the median taxon richness increased from a minimum in freshwater to a maximum in marine water (123 to 166 OTUs), while the Shannon index increased from 2.7 in freshwater to 3.5–3.8 in brackish and marine water. Sediment community Shannon diversity and richness decreased along this gradient (Shannon indices from 6.0 to 5.9; OTUs from 596 to 502). While greater diversity or taxon richness in marine water versus freshwater may be attributed to OTUs being flushed out to sea, we observed no evidence for this, as very few OTUs (<3%) were shared between these systems (Supplementary Figures 2b and 3). The same patterns were observed when considering specific lineages (archaea, bacteria and microeukaryotes), except benthic microeukaryotes displayed no clear trend (Fig. 2b). Results showed no obvious reduction in prokaryotic or microeukaryote richness in the brackish zone, in contrast with Remane’s species minimum concept, which proposes minimum richness for benthic macrofauna in the brackish horohalinicum [4]. Our results on prokaryotes and microeukaryotes, taken together with other studies that have found a protistan species-maximum in brackish water [5, 49], suggest essential differences in taxa diversity associated with salinity gradients.

Osmoregulation, nutrient metabolism and primary production across environments

Differences in functional gene abundance and expression across sampled gradients and relationship to taxa composition

Overall, metabolic profiles (genes and genes expressed) decreased in similarity with increasing salinity divergence (Supplementary Figure 2c). However, neither taxonomic nor functional gene data were strong indicators of gene expression, as previously observed for freshwater-to-marine planktonic [14], ocean [50] and hypersaline desert environments [51]. Mantel tests showed that while taxonomic beta-diversity was strongly correlated with differences in functional potential (ρ = 0.90, p < 0.01), it was only weakly correlated with differences in gene expression (ρ = 0.54, p < 0.01). Similarly, we identified a moderate correlation between differences in gene and transcript abundances (ρ = 0.70, p < 0.01). Microbial composition is thought to be a poor predictor of ecological processes due to functional redundancy and variations in environmental response, leading to marked differences in transcript expression [50,51,52,53]. To ascertain which metabolic features, associated with osmoregulation, nutrient metabolism and primary production, were strongly differentiated among freshwater, brackish and marine environments, and between water and sediment, we applied linear discriminant analysis effect size (LEfSe) and weighted gene co-expression network analysis (WGCNA). LEfSe results revealed 25 discriminative features from the metagenomes and 42 from the metatranscriptomes (Fig. 3a; Kruskal–Wallis p < 0.05, Supplementary Table 5), indicating a large shift in metabolic potential and, most notably, activity across this transect (Supplementary Figure 4). In agreement with LEfSe results, WGCNA analysis showed 4 out of 6 co-expressed gene clusters were significantly associated with environment (Supplementary Figure 5; details are described in Supplementary Information), revealing a clear division in metabolism related to nutrient acquisition and osmoregulation between sediment and water and between saline and non-saline environments. Genes and transcript abundances related to phosphorus metabolism, photosynthesis and osmoregulation were significantly higher in fresh, brackish and marine water, respectively (Fig. 3 and Supplementary Figure 5). In contrast, brackish sediment was characterised by significantly higher sulfur metabolism (Fig. 3a and Supplementary Figure 5). Nutrients in tidal lagoon or bar-built estuaries, such as in this study, are generally derived from runoff from land or groundwater inputs, and promote primary productivity and nutrient regeneration in brackish systems [54]. Fewer genes and transcripts were identified as significantly and differentially abundant among benthic environments due to large variations in their relative abundance, reflecting the high spatial heterogeneity in benthic community composition and nutrient availability (Figs. 1b and 2a).

Fig. 3

Discriminative features and distribution of functional potential and metabolic activity across the salinity gradient. a Differentially abundant metabolic functions from metagenomic (top, n = 36) and metatranscriptomic (bottom, n = 33) data, as determined by LEfSe (specific pathways from Supplementary Table 5). Arrows indicate functions were significantly over-represented in their specific habitat. b Pie charts showing the relative fraction of functional gene categories across the freshwater-to-marine gradient based on metagenome and metatranscriptome data. Transcriptomic data were not obtained for marine sands due to low RNA quality. Colours in both panels indicate key nutrient acquisition and osmoregulation pathways

Distinct non-saline and saline osmoregulation strategies

In line with previous research indicating that microbial taxonomic distributions are primarily driven by salinity [14], largely due to lineage-specific adaptations to salt tolerance [15, 17, 55], our results show that dominant microbial taxa in non-saline (Betaproteobacteria) and saline (Alphaproteobacteria and Gammaproteobacteria) habitats adopted distinct osmoregulation strategies (Figs. 4 and 5 and Supplementary Figure 6). While K+ acquisition and Na+ export are considered important mechanisms for osmoregulation in brackish or marine environments [14, 15], results here illustrate that K+ transport is also crucial in freshwater communities. We observed a significantly greater abundance of the potassium (K+) transporter gene kdp, in non-saline sediment (Kruskal–Wallis, p < 0.05, Fig. 4 and Supplementary Figure 4). In contrast, trkA was more abundant in marine water and sediment (Fig. 4). The trkA system is prominent under hyperosmotic stress [19], characteristic of the marine environment. The repressible high-affinity K+ transporter, kdp, is known to be active at low external K+ concentrations [21], consistent with freshwater environments.

Fig. 4

Gene and gene expression profiles across the salinity gradient. Line plots showing log10-normalised gene and transcript abundances related to osmoregulation, nutrient acquisition and primary production pathways across the salinity gradient. Genes per category are given in Supplementary Table 9. Asterisks (*) indicate significant correlations between metagenomes and metatranscriptomes (Spearman ρ > 0.5, p < 0.05). Abbreviations: TPM genes/transcripts per million mapped reads

Fig. 5

Schematic illustrating metabolic genes and gene expression in benthic and planktonic environments across the freshwater-to-marine transect. Genes with relative abundances and expression levels at least twice higher than in those of other environments are shown. No transcriptomics data were available for marine sediments at sites 8 and 9. Key on right: rectangles with numbers represent metabolic pathways; ovals with letters indicate the main taxa expressing these genes

We found a significantly higher abundance of Na+/H+ antiporter genes (mnhABCDEFG) in marine sediment, and higher expression of Na+ translocating dehydrogenase (nqrF) and glycine betaine transporter (bcct and proXV) genes in marine water (Figs. 3a and 4, Supplementary Figure 4), primarily contributed by Alphaproteobacteria and Gammaproteobacteria, with an apparent increase in gene expression starting from site 7 (salinity = 26; Supplementary Figure 6). Glycine betaine transporter and Na+/H+ antiporter gene expression in the water was also significantly correlated with salinity (ρ = 0.90–0.95, p < 0.01, Supplementary Figure 6). In the marine environment, some bacteria rely on Na+-motive force to routinely export Na+ via various Na+/H+ antiporters and Na+ translocating dehydrogenases [56,57,58,59]. Increased abundance and expression of Na+ transporter genes with salinity suggests that brackish and marine microorganisms increasingly exported destabilising Na+ in exchange for greater accumulation of K+ within cells, which stabilise acidic salt-adapted proteins [17, 60, 61]. Glycine betaine transporters, in contrast, uptake organic solutes to achieve osmotic equilibrium [17, 62], or for metabolism and ATP generation [63].

Spatial distribution of light and dark primary production

Autotrophic microbial and microeukaryotic communities play crucial roles in aquatic primary production and form the base of the aquatic food chain. We observed taxonomically diverse organisms across the river–marine transect with capabilities for oxygenic or anoxygenic photosynthesis and carbon fixation via the Calvin–Benson, Wood–Ljungdahl and reverse tricarboxylic acid (TCA) pathways (Figs. 5 and 6). Gene expression related to primary production was the highest in water in the mid-brackish zone (salinity = 24; Figs. 4 and 6). This corresponded with increased phytoplankton relative abundance, particularly Cryptophyta (genus Hemiselmis; Figs. 2a and 6), and a peak in ammonium (0.04 g/m3, Fig. 1b)—the preferred nitrogen source for Cryptomonads [64] (Supplementary Information). Data indicate phytoplankton actively assimilated carbon via the Calvin–Benson cycle (RuBisCO) and conserved energy via oxygenic photosynthesis (photosystems I and II; Figs. 5 and 6). In contrast, diatoms can photosynthesise and fix carbon efficiently under low light [65]. Coastal lagoon systems, including the Waiwera river, are shallow water bodies (<5-m depth) [66], where light can penetrate to the bottom sediments and promote benthic primary productivity [54]. Consistent with this, Waiwera had a low median river turbidity of 7.0 Nephelometric Turbidity Units (NTU) reported for the sampling year (2018,, and results show that diatoms dominated sediment primary production across the brackish transect (Fig. 6). Increased phytoplankton primary production, and consequently, sinking biomass, may explain the higher concentrations of total sedimentary organic carbon, nitrogen and phosphorus measured at the mid-brackish site (Fig. 1b). High benthic primary production (Figs. 3b and 6) and mud content may also contribute to higher organic carbon in brackish sediment [67,68,69,70] (average 17.8% versus 4.3% elsewhere, Supplementary Table 4).

Fig. 6

Distribution of taxon-specific genes and transcripts involved in photosynthesis and carbon fixation across the study sites. NA indicates that no transcriptomic data are available for marine sediments. Abbreviations: WGS metagenome, metatranscriptome (WTS), genes/transcripts per million mapped reads (TPM)

There is limited information on carbon fixation pathways utilised across salinity gradients, but a few studies show that estuaries are global hotspots for carbon cycling and assimilation [28, 71, 72]. Here, the Calvin cycle, in conjunction with photosynthesis, was significantly more highly expressed in brackish water (significant, LDA score ≥ 2), whereas genes related to dark carbon fixation, such as the reductive TCA cycle and Wood–Ljungdahl pathway, were primarily expressed in sediment environments. While we did not measure dissolved oxygen in this study, microorganisms undertaking dark carbon fixation typically thrive across a redox gradient, such as the oxic–anoxic interface in aquatic sediments [73, 74]. Both the Wood–Ljungdahl pathway and reductive TCA cycle are strongly associated with anaerobes, although the latter is also found in microaerobes and some aerobes [75, 76]. Among organisms undertaking dark carbon fixation, we found that the reverse TCA cycle genes (aclA) were significantly more highly expressed in non-saline sediment (LDA score ≥2), and predominantly by Nitrospira (Fig. 6), which our results indicate also undertook complete ammonia oxidation (comammox) in this environment (Fig. 5). In contrast, Wood–Ljungdahl pathway codh genes were highly expressed in brackish sediments, primarily by Chloroflexi and Deltaproteobacteria (Figs. 5 and 6). Wood–Ljungdahl pathway gene expression in the brackish sites was at least 100 times greater than reverse TCA cycle gene expression, but 10 times less than the Calvin cycle gene expression (Figs. 5 and 6), highlighting the significance of the Calvin cycle and Wood–Ljungdahl pathway for carbon assimilation in this brackish system.

Nitrogen and phosphorus acquisition in freshwater and brackish environments

Ammonium was the main inorganic nitrogen source observed in benthic environments (0.1–2.7 g/m3), but was relatively scarce in the water column (<0.04 g/m3), where nitrate and nitrite concentrations were also lower (Fig. 1b and Supplementary Table 4). While ammonium, derived from anthropogenic activities, can be the dominant form of N entering shallow tidal lagoon settings [54], such riverine inputs can also be nitrate-dominated [77], as observed in this study (Fig. 1b). This suggests considerable benthic regeneration of ammonium [54] in the sampled sediment (top 2 cm). Higher ammonium concentrations in the benthic zones may also be due to sedimentation of phytoplankton blooms, boosting both nitrate uptake and ammonium release in benthic environments [78, 79]. Nitrate is therefore depleted during spring bloom progression [54]. To deal with fluctuating nitrogen availability in aquatic ecosystems, microbial communities are equipped with strategies for nitrogen assimilation (e.g. fixation) and storage (cyanophycin) [80, 81]. Genes related to nitrogen fixation (nifABDH) were detected and expressed in non-saline and brackish sediment (Figs. 4 and 5; no transcriptomic data are available for marine sediment), although ammonium levels were around 1 g/m3 in some of these sites (Fig. 1b and Supplementary Table 4). Excess ammonium (>50 g/m3) or oxygen (>2 %) is reported to inhibit nitrogenase activity [82,83,84]. However, as ammonium concentrations tested were orders of magnitude higher than those we detected, ammonium concentrations in the present study may not have affected nitrogenase activity. In contrast, there were few to no nifABDH transcripts detected in the overlying water, likely due to low predicted iron and molybdenum concentrations [85], which is characteristic of marine water [81]. High iron and molybdenum are required for the FeMo nitrogenase cofactor in all known diazotrophs [86]. In some bacteria, excess nitrogen assimilated is stored intracellularly as cyanophycin granules comprised of nitrogen and carbon storage polymers [80, 87]. Our results indicate that genes related to cyanophycin (cphA) were significantly more abundant and highly expressed by Betaproteobacteria in non-saline water (Figs. 3a and 5, and Supplementary Figure 7), suggesting active storage and utilisation of cyanophycin granules.

Dissolved reactive phosphorus (DRP) concentrations were the highest in marine sediment pore water (0.06 g/m3 ± 0.03), presumably owing to the naturally higher rate of phosphorus mobilisation in marine sediments [23], and were lowest in freshwater environments (<0.01 g/m3, Fig. 1b). Phosphorus retention/release from sediment is complex, but is mainly controlled by increased desorption under saline conditions [88, 89]. Previous studies have shown that the high-affinity inorganic phosphate transporter pst system and phosphate regulon phoB-phoR (which regulates inorganic phosphate uptake) are highly upregulated under phosphate-limited conditions [90, 91]. Accordingly, data here shows that genes for phosphate uptake (pstSABC) and regulation (phoU, phoB and phoR) were negatively correlated with phosphate concentration (ρ = −0.49, p < 0.05) and were significantly more highly expressed in phosphate-limited non-saline water (Fig. 3 and Supplementary Figures 5 and 7)—primarily by Betaproteobacteria (Burkholderiales), which were highly abundant in this environment (Fig. 2a). Depending on phosphate availability, phoR can act as a kinase (phosphate-limited conditions) or phosphatase (phosphate-sufficient) and thereby activate or interrupt phoB phosphate regulation activity, respectively [91]. In phosphate-limited conditions, the phosphorylated phoB activates the high-affinity pst operon [90, 91]. Our results indicate a reliance on efficient phosphate uptake and regulation mechanisms in phosphate-limited freshwater.

Nitrogen cycling dynamics across the salinity gradient

Genes associated with nitrification and denitrification were predominantly expressed in brackish sediments (Fig. 5), corresponding with higher total organic carbon (TOC), total nitrogen (TN) and total recoverable phosphorus (TRP) concentrations (Fig. 1b). Results indicate that metabolic coupling between distinct taxa, to oxidise ammonia fully to nitrate, was far more common than organisms undertaking comammox along the transect (Supplementary Figure 8). We detected relatively low expression of amoA genes involved in comammox by the class Nitrospira (the genus Nitrospira is well known for nitrification and comammox [92]) at a single benthic freshwater site (0.9 km along the transect, Supplementary Figure 8), where nitrate was highest and ammonium concentration was lowest (Fig. 1b). In contrast, genes involved in nitrification via ammonia oxidation to nitrite (amoA in archaea) and nitrite oxidation (nxrAB in Nitrospirae) were highly expressed in the brackish sediments (Fig. 5 and Supplementary Figure 8). The relative abundances of ammonia-oxidising bacteria and archaea can vary in coastal and estuarine sediments, depending on the environmental conditions (salinity, oxygen and ammonium concentration) [93,94,95]. Although archaea constituted only a minor fraction of the overall microbiome (Fig. 2a), gene expression data indicates they were the primary ammonia oxidiser in brackish sediments and are important contributors to estuarine biogeochemical cycling.

As for nitrification, genes related to nitrous oxide reduction (nosD, nosZ), the final step of the denitrification pathway, were predominantly expressed in brackish sediments (Fig. 4) and were primarily associated with Bacteroidetes (Fig. 5 and Supplementary Figure 8). Results therefore suggest a strong coupling of nitrification and denitrification in brackish sediments, leading to nitrogen loss from this estuary. Steep redox gradients promote the co-occurrence of aerobic and anaerobic metabolic pathways, including the aforementioned dark carbon fixation reactions and coupled nitrification–denitrification [31, 96]. Such a gradient may be expected in brackish sediments given the combination of high mud and organic carbon contents (Fig. 1b and Supplementary Table 4), which are associated with hypoxic/anoxic conditions due to lower diffusion [97] and rapid consumption of oxygen within sediments, combined with oxygen generated by photosynthesizing diatoms at the sediment–water interface (Figs. 5 and 6).

Unlike denitrification, dissimilatory nitrate reduction to ammonium (DNRA) retains bioavailable nitrogen [98]. The first two steps in denitrification and DNRA are nitrate and nitrite reduction, of which the second step involves distinct nitrite reductase enzymes (NirKS for denitrification; NirBD/NrfA for DNRA) [99]. Results show that DNRA genes (nrfA, nrfH) were highly expressed in non-saline sediment (Fig. 4 and Supplementary Figure 8). The ratio of denitrification-nosZ to DNRA-nrfA gene expression in freshwater sediment (ratio = 0.06) was on average 20 times lower than in brackish sediment (ratio = 1.56, Supplementary Table 6), suggesting a partitioning between DNRA in freshwater sediment and denitrification in brackish sediment. Previous studies have found that a higher ratio of organic carbon (electron donor) to nitrate and a greater supply of nitrate relative to nitrite favour DNRA over denitrification [78, 100,101,102]. However, we did not find any relationship between carbon/nitrate or nitrate/nitrite ratios and DNRA-related gene expression in sediment (Supplementary Table 6). Although nitrate/nitrite ratios decreased across the transect (primarily due to a corresponding decrease in nitrate), nitrate concentrations were on average over ten times higher for most brackish sediment (ratios 4–34, distances ≤ 2.6 km along the transect) (Fig. 1b and Supplementary Table 4). This suggests that other unmeasured factors, such as sulfide concentrations or microbial generation times [100, 103], may instead have influenced denitrifier and DNRA activity.

Sulfur metabolism was most active in the brackish–marine environment

Like nitrification–denitrification, genes for sulfur metabolism, including sulfur reduction (cysN, sat and aprA for sulfate; dsrABD for sulfite) and oxidation (rdsrAB for sulfur), were significantly more highly expressed in brackish sediments (Fig. 3a, Supplementary Figures 9-10). Expression was concomitant with highly abundant Gammaproteobacteria and Deltaproteobacteria capable of oxidising sulfur and reducing sulfate, respectively (Figs. 2a and 5). This is consistent with greater sulfide accumulation in brackish sediments owing to high mud and limited oxygen diffusion [30]. Taken together, these results reinforce the findings of studies indicating the brackish environment is a hotspot for biogeochemical cycling [104, 105]. However, the expression of thiosulfate oxidation genes (soxBCY) by Alphaproteobacteria was significantly higher in marine water (LDA score ≥ 2, Fig. 3a and Supplementary Figure 9) and positively correlated with high sulfate concentrations typical of marine water (ρ = 0.76, p < 0.05; Fig. 1b and Supplementary Figure 5). Thiosulfate can serve as an electron donor for anoxygenic photosynthesis [106] and be oxidised to sulfate via the Sox pathway [107]. Bacterial groups that perform aerobic anoxygenic photosynthesis (AAP), including alpha-, beta- and gamma-proteobacteria, are widely distributed in the surface ocean [108, 109]. Greater expression of genes related to thiosulfate oxidation (Supplementary Figure 9) and anoxygenic photosynthesis (pufLM) by planktonic marine Alphaproteobacteria (Figs. 5 and 6) suggests these bacteria may have coupled thiosulfate oxidation with photoheterotrophy in relatively carbon-limited marine water (DNPOC 1.7 versus 2.2 g/m3 in brackish water) [110].


Across the river–marine transect, in a tidal lagoon setting, the microbial community and diversity were dominated by bacteria. Contrasting diversity trends were observed between the water column and sediments, and between prokaryotes and microeukaryotes, reflecting the strong geochemical differences between these environments, and the distinct physiologies and ecosystem roles of these taxa groups. While archaea and microeukaryotes represented only a minor fraction of the overall communities, they played significant roles in ammonia oxidation and photosynthesis, respectively. Microorganisms employed distinct osmoregulation strategies between freshwater and saline habitats by actively expressing diverse osmoregulation genes encoding potassium, sodium and glycine betaine transporters across the transect. Results indicate a significant predominance of phosphate acquisition in non-saline water, along with microeukaryotic-driven primary production, and prokaryotic nitrification–denitrification and sulfur metabolism in the brackish sediments, corresponding with large differences in resource availability. Overall, our study demonstrates that differences in salinity and nutrients impose significant biological boundaries between sediment and water, and among non-saline, brackish and marine environments.

Materials and methods

Sampling strategy

Sampling was undertaken on 12 November 2018 within a low tide period. Nine sites (sites 1–9, Fig. 1a; site photos, Supplementary Figure 11) were selected covering a freshwater-to-marine salinity (S) gradient (0.2 < S < 35) across a 5-km-long transect of the Waiwera river and estuary, north of Auckland, New Zealand (36° 33′ 02.4″ S, 174° 39′ 08.1″ E). For molecular analyses, at least 10 L of water was collected between 0.2 and 1 m below the surface of all sites and sequentially filtered through sterile 1.2-μm and 0.22-μm mixed cellulose ester filters (Merck Millipore, MA, USA). Triplicate water samples (unfiltered) were similarly collected in polyethylene bottles for chemical analysis at Hill Laboratories (Hamilton, New Zealand). The top 2 cm of sediment was collected at three locations approximately 5 m apart at each site and wet sieved through a 1-mm mesh sieve to remove stones, plant material and macrofauna. All filter and sediment samples for molecular analyses were immediately preserved with LifeGuard Soil Preservation Solution (Qiagen, MD, USA) in sterile 50-mL Falcon tubes (1:3 wet sediment-to-LifeGuard ratio, individual water samples were filtered and preserved within 5–10 min of sample collection) and stored at −80 °C until extraction. Samples for mud content and chemical analyses were stored at 4 °C.

Mud content measurement

Sediment samples were oven dried at 105 °C for 12 h prior to sieving. Sediments were serially sieved using mesh sizes ranging from 1 mm to 63 μm with a sieve shaker. Fractions retained on each mesh were weighed to determine mud and sand contents. Mud (clay and silt) and sand are identified as sediment with sizes of <63 μm and >63 μm, respectively [111]. Mud content was calculated as follows: the weight of mud fraction/weight of sediments from all fractions × 100 %.

Physiochemical and nutrient analyses

Temperature, pH and salinity were measured on-site using a Pocket Pro+ Multi 2 Tester (HACH, CO, USA). All sediment and water samples were sent to Hill Laboratories (Hamilton, New Zealand) for processing. All nutrients were measured according to standard American Public Health Association (APHA) protocols [112]. TRP, TS, TN and TOC were only measured from sediment samples. All other nutrients were measured from water column and sediment pore water samples. Briefly, sediment TOC, following pre-treatment with acid to remove carbonates, was analysed using an Elementar Vario MAX Combustion Analyser (Langenselbold, Germany) with catalytic combustion (at 900 °C) and TN was measured via an Lachat Quikchem Series 2 Flow injection Analyser (Lachat Instruments, CO, USA). Total sulfur was measured using a LECO SC-32 Sulfur Determinator (MI, USA), and TRP was analysed using a Aquakem Konelab 600 Discrete Analyser (Thermo Fisher Scientific, MA, USA). Water samples were prefiltered through 0.45-μm membrane filters and saline samples diluted prior to analyses. Total ammonical-N, nitrate, nitrite and DRP were measured via a Lachat Quikchem Series 2 flow injection analyser (Lachat Instruments, CO, USA). Sulfate was analysed using a Dionex Ion Chromatography system (Sunnyvale, CA, USA). Samples for DNPOC were acidified to remove inorganic C before addition of persulfate, heating to a temperature of above 375 °C and analysis via a Sievers Innovox TOC analyser (SUEZ Analytical Instruments, CO, USA). River turbidity data, collected by Auckland Council at GPS coordinates 36° 33′ 01.4″ S 174° 39′ 59.3″ E, were obtained from the Land, Air, Water Aotearoa (

Nucleic acid extraction and sequencing

Both RNA and DNA were co-extracted using RNeasy PowerSoil Total RNA Kits and RNeasy PowerSoil DNA Elution Kits (Qiagen, MD, USA). Total RNA and DNA from each sample were isolated from 2–8 g of sediment or 1–5 g of filter (1/3-1 filter, 10L water per filter) with collected planktonic biomass (both 1.2 μm and 0.22 μm combined). Extracted nucleic acids were stored at −80 °C. RNA samples were treated with DNaseI. DNA removal was verified by PCR amplification of the 16S rRNA genes using 515F/806R primers [113, 114], JumpStart REDTaq Readymix (Sigma-Aldrich, MO, USA) and 2 μL of the template (sample or positive control), under the following conditions: 95 °C for 5 min; 55 cycles: 95 °C for 45 s, 50 °C for 60 s, 72 °C for 90 s; 72 °C for 10 min. Visualisation by gel electrophoresis was used to confirm no amplification in DNase-treated samples. RNA was purified using RNA Clean & Concentrator-5 Kits (Zymo Research, CA, USA). DNA extracts were further purified by mixing with nuclease-free 0.3 M sodium acetate (final concentration; Sigma-Aldrich, MO, USA) and 0.1 μg/μL glycogen (final concentration) and re-precipitating with 2:1 volumes of cold ethanol following incubation for 30 min at −20 °C. Pellets were washed with 70% ethanol and re-eluted in nuclease-free water. High molecular weight DNA was confirmed via gel electrophoresis. DNA and RNA concentrations and quality were further assessed using a Nanodrop 2000c Spectrophotometer (Thermo Fisher Scientific, MA, USA) and also, for RNA, an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA) using the Agilent RNA 6000 Nano Kit. Low-quality RNA extracts (without obvious 16S/23S peaks in bioanalyzer electropherogram) from sandy sediment at marine sites 8 and 9 were excluded from further analysis. The remaining 36 DNA and 30 RNA extracts were prepared for metagenomic and metatranscriptomic sequencing. Prior to RNA library preparation, ribosomal RNAs were depleted using the Ovation Universal Prokaryotic RNASeq, Prokaryotic AnyDeplete (TECAN, Zürich, Switzerland) system, according to the manufacturer’s instructions. Library preparation and sequencing were undertaken by the Otago Genomics Facility (University of Otago, Dunedin, New Zealand). Metagenome libraries were prepared using the Illumina TruSeq DNA Nano kit (with 550-bp insert sizes). Metatranscriptome libraries were prepared using the Ovation Complete Prokaryotic RNASeq kit. Paired-end 2 × 250 bp sequencing for DNA (over 14 lanes) and 2 × 125 bp sequencing for RNA (over 14 lanes) were performed using the Illumina HiSeq 2500 platform with V2 (DNA) and V4 (RNA) chemistry. All sequence data have been deposited with NCBI under BioProject PRJNA668816.

Sequence read processing

Ribosomal RNA reads were first removed from metatranscriptomes using SortMeRNA v2.1 [115] with bacterial, archaeal and eukaryotic rRNA databases, and results were passed to (BBMap [116]) to regenerate paired-end reads. Metagenomic and non-ribosomal metatranscriptomic read quality was inspected using FastQC [117]. Reads were trimmed using Trimmomatic v0.38 [118] with a minimum Phred score of 30 and adapter sequence removal via the ILLUMINACLIP option with Illumina universal adapter sequence (AGATCGGAAGAG). Reads were also trimmed to a maximum of 240 bp (DNA) and 115 bp (RNA) long to remove additional low-quality bases, and reads <80 bp (DNA) and <60 bp (RNA) long were discarded. A sediment sample from site 9, replicate 1 (S9R1), was subjected to additional trimming (headcrop 10) due to low base quality.

Small subunit (SSU) rRNA gene reconstruction

Near full-length 16S and 18S rRNA gene sequences were reconstructed from each metagenome, over 50 iterations and using a joining threshold of 97%, with EMIRGE [119] and the SILVA SSU Ref NR 99 138 database. We detected and removed 2670 chimeras with the uchime_ref command [120]. Sequences from across samples were concatenated and clustered into OTUs (97% threshold) with UCLUST [120], then rarefied to 3343 sequences per sample, resulting in a total of 7334 unique OTUs with a minimum count across samples ≥2 (Supplementary Table 7). SSU rRNA gene relative abundances were based on normalised priors from the final EMIRGE iteration.

Metagenome assembly, taxonomic and functional annotations

Trimmed reads from individual water samples or from spatial triplicate sediment samples from the same site were assembled or co-assembled, respectively, using metaSPAdes [121] with k-mer values: 41, 61, 81, 101 and 127, resulting in 143 million contigs (Supplementary Figure 12). A k-mer-based profiler, Kraken2 v2.1.2 [122], was used to assign taxonomy to the trimmed reads based on compositional similarity to publicly available genome databases from RefSeq [123]. Species relative abundances were further estimated using a Bayesian approach with Bracken v2.6.0 [124]. Data were rarefied to 2,543,737 reads per sample, and species belonged to Streptophyta and Chordata were removed, resulting in a total of 5136 unique species (Supplementary Table 8). Open reading frames (ORFs) were predicted and translated into protein sequences using Prodigal v2.6.3 [125]. Protein sequences were searched against a library of hidden Markov models (HMMs) consisting of TIGRFAMs [126], Pfam [127] and custom HMMs for metabolic pathways using HMMER3 [128]. Marker genes for the pathway of interest and the cut-off values for HMM scores were derived from Anantharaman et al. [129], Pfam [127], TIGRFAM [126] and FunGene [130] databases (Supplementary Table 9). Phylogenetic analyses were conducted on concatenated dsrAB genes to further classify them into dsrAB (in sulfate/sulfite-reducing bacteria) and rdsrAB (in sulfur-oxidising bacteria) using FastTree v2.1.10 [131]. Reference protein sequences for phylogenetic analyses were obtained by BLAST searches against the NCBI nr protein database [132]. Taxonomic affiliations of functionally annotated genes were assigned with DIAMOND BLASTP [133] with default settings against the NCBI nr protein database (downloaded on 16 February 2020), using a weighted lowest common ancestor (LCA) approach. We retained only hits with at least 70% query coverage, 80% identity and bit scores within 5% of the top hit for LCA identification.

Metagenome and metatranscriptome read mapping and counts

DNA and RNA (non-rRNA) reads were mapped back to metagenomic contigs using BBMap v37.93 [116] with default parameters. Read counts were calculated for each predicted open reading frame (ORF) using featureCounts v2.0.0 [134]. Gene and transcript counts were normalised using gene or transcript per million mapped reads (TPM) to measure the bulk contribution of the community [135, 136]. TPM = (number of reads mapped to the gene/gene length)/sum (number of reads mapped to the gene/gene length) × 1,000,000.

Statistical analyses

Unless otherwise stated, the following statistical analyses were carried out in R environment version 3.5.1 [137]. Shannon–Wiener indices, Simpson’s indices and relative OTU richness were determined using the vegan package [138] and the rarefied OTU table combining reconstructed 16S and 18S rRNA genes. A Wilcoxon signed-rank test was performed to determine significant differences in observed diversity between the water column and sediment. Distance–decay relationships for microbial communities and functional gene groups were constructed by first calculating pairwise Bray–Curtis dissimilarities among sites using the vegdist function in the vegan package [138], and fitting negative exponential functions via a log-linked generalised linear model (GLM) or local polynomial regression (LOESS) model. NMDS ordinations were plotted based on Bray–Curtis dissimilarities, and environmental variables were fitted with 999 Monte Carlo permutation tests using the vegan package [138]. Pairwise permutational multivariate analyses of variance (PERMANOVA) with Benjamini–Hochberg corrections were carried out to test for statistically significant variance among multivariate microbial community data using the pairwiseAdonis package [139]. Spearman’s correlation coefficients between the relative abundances of specific taxa and environmental parameters were calculated using the R cor function. Taxa with p < 0.05 were determined to be significantly correlated with the environmental parameters. The bioenv function in R’s vegan package [138] was used with Spearman’s rank correlations to search over subsets of the continuous environmental variables to determine the best explanatory variables for community composition. A multivariate regression tree was plotted using mvpart package [140], with 1000 cross-validations and tree size = 8 to identify the main predictors that correlate with microbial community structure. As both salinity and distance were highly correlated and exhibited multicollinearity, we used salinity (and excluded distance) as one of the predictor variables for the analysis. The decision tree identifies variance in microbial community composition (Bray–Curtis dissimilarity) caused by threshold values of key site physicochemical attributes. The values attached to each branch of the tree mark the physicochemical criteria used by the regression tree to group samples based on differences in community composition.

The correlation between gene and transcript relative abundances was measured using the R cor.test function. Mantel tests were run to assess Spearman’s correlations between functional and taxonomic community dissimilarity matrices based on Bray–Curtis dissimilarities. LEfSe was applied to metagenomic and metatranscriptomics data to identify environment-associated metabolic functions by utilising Kruskal–Wallis and Wilcoxon rank-sum tests to determine variations occurring between and within environments with a linear discriminant analysis (LDA) score of ≥2 [141]. To determine co-expressed genes associated with environmental parameters, weighted correlation network analysis was also performed using transcript counts normalised by TPM as input data to find clusters of highly correlated gene transcripts by using WGCNA package [142] with minimum module size = 5 and soft-thresholding power = 14 (R2 = 0.81) to relate the functions or modules to one another and to environmental parameters.

Availability of data and materials

The data generated in the current study are all publicly available. All sequence data have been deposited with NCBI under BioProject PRJNA668816. R scripts for plotting figures and statistical analyses are available at


  1. 1.

    Allan JD, Ibañez Castillo MM. Stream ecology: structure and function of running waters. 2nd ed. Dordrecht: Springer; 2007.

    Book  Google Scholar 

  2. 2.

    Heinrichs ME, Mori C, Dlugosch L. Complex interactions between aquatic organisms and their chemical environment elucidated from different perspectives. In: Jungblut S, Liebich V, Bode-Dalby M, editors. YOUMARES 9 - Oceans Our Res Our Future. Cham: Springer; 2020. p. 279–97.

    Chapter  Google Scholar 

  3. 3.

    Cloern JE, Jassby AD, Schraga TS, Nejad E, Martin C. Ecosystem variability along the estuarine salinity gradient: examples from long-term study of San Francisco Bay. Limnol Oceanogr. 2017;62(S1):272–91.

    CAS  Article  Google Scholar 

  4. 4.

    Remane A. Die brackwasserfauna. Verh Dt Zool Ges. 1934;36:34–74.

    Google Scholar 

  5. 5.

    Filker S, Kühner S, Heckwolf M, Dierking J, Stoeck T. A fundamental difference between macrobiota and microbial eukaryotes: protistan plankton has a species maximum in the freshwater-marine transition zone of the Baltic Sea. Environ Microbiol. 2019;21(2):603–17.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Herlemann DP, Labrenz M, Jürgens K, Bertilsson S, Waniek JJ, Andersson AF. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 2011;5(10):1571–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Herlemann DPR, Lundin D, Andersson AF, Labrenz M, Jürgens K. Phylogenetic signals of salinity and season in bacterial community composition across the salinity gradient of the Baltic Sea. Front Microbiol. Frontiers. 2016;7:1883.

    PubMed  PubMed Central  Article  Google Scholar 

  8. 8.

    Wang J, Yang D, Zhang Y, Shen J, van der Gast C, Hahn MW, et al. Do patterns of bacterial diversity along salinity gradients differ from those observed for macroorganisms? PLoS One. Public Library Sci. 2011;6(11):e27597.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Lozupone CA, Knight R. Global patterns in bacterial diversity. Proc Natl Acad Sci. 2007;104(27):11436–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  10. 10.

    Thompson LR, Sanders JG, McDonald D, Amir A, Ladau J, Locey KJ, et al. A communal catalogue reveals earth’s multiscale microbial diversity. Nature. Nature Publishing Group. 2017;551:457–63.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Mason OU, Canter EJ, Gillies LE, Paisie TK, Roberts BJ. Mississippi River plume enriches microbial diversity in the Northern Gulf of Mexico. Front Microbiol. Frontiers. 2016;7:1048.

    PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Doherty M, Yager PL, Moran MA, Coles VJ, Fortunato CS, Krusche AV, et al. Bacterial biogeography across the Amazon river-ocean continuum. Front Microbiol. Frontiers. 2017;8:882.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Dupont CL, Larsson J, Yooseph S, Ininbergs K, Goll J, Asplund-Samuelsson J, et al. Functional tradeoffs underpin salinity-driven divergence in microbial community composition. PLoS One. Public Library of Science. 2014;9:e89549.

    PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Fortunato CS, Crump BC. Microbial gene abundance and expression patterns across a river to ocean salinity gradient. PLoS One. 2015;10(11):e0140578.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Walsh DA, Lafontaine J, Grossart H-P. On the eco-evolutionary relationships of fresh and salt water bacteria and the role of gene transfer in their adaptation. In: Gophna U, editor. Lateral Gene Transf Evol. New York: Springer; 2013. p. 55–77.

    Chapter  Google Scholar 

  16. 16.

    Henson MW, Lanclos VC, Faircloth BC, Thrash JC. Cultivation and genomics of the first freshwater SAR11 (LD12) isolate. ISME J. Nature Publishing Group. 2018;12:1846–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Oren A. Life at high salt concentrations. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes. 4th ed. Berlin: Springer; 2013. p. 421–40.

    Chapter  Google Scholar 

  18. 18.

    Trchounian A, Kobayashi H. Kup is the major K+ uptake system in Escherichia coli upon hyper-osmotic stress at a low pH. FEBS Lett. 1999;447(2-3):144–8.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Trchounian A, Kobayashi H. K+ uptake by fermenting Escherichia coli cells: pH dependent mode of the TrkA system operating. Biosci Rep. 2000;20(4):277–88.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Bañuelos MA, Klein RD, Alexander-Bowman SJ, Rodríguez-Navarro A. A potassium transporter of the yeast Schwanniomyces occidentalis homologous to the Kup system of Escherichia coli has a high concentrative capacity. EMBO J. 1995;14(13):3021–7.

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Rhoads DB, Waters FB, Epstein W. Cation transport in Escherichia coli. VIII. Potassium transport mutants. J Gen Physiol. 1976;67(3):325–41.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Paver SF, Muratore D, Newton RJ, Coleman ML. Reevaluating the salty divide: phylogenetic specificity of transitions between marine and freshwater systems. Flynn TM, editor. mSystems. 2018;3:e00232–18.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Caraco N, Cole J, Likens GE. A comparison of phosphorus immobilization in sediments of freshwater and coastal marine systems. Biogeochemistry. 1990;9(3):277–90.

    CAS  Article  Google Scholar 

  24. 24.

    Xia X, Zhang S, Li S, Zhang L, Wang G, Zhang L, et al. The cycle of nitrogen in river systems: sources, transformation, and flux. Environ Sci Process Impacts. 2018;20(6):863–91.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Berner RA. Sedimentary pyrite formation. Am J Sci. 1970;268(1):1–23.

    CAS  Article  Google Scholar 

  26. 26.

    Emeis KC, Struck U, Leipe T, Pollehne F, Kunzendorf H, Christiansen C. Changes in the C, N, P burial rates in some Baltic Sea sediments over the last 150 years - relevance to P regeneration rates and the phosphorus cycle. Mar Geol. Elsevier. 2000;167(1-2):43–59.

    CAS  Article  Google Scholar 

  27. 27.

    Mcdowell RW, Larned ST, Houlbrooke DJ. Nitrogen and phosphorus in New Zealand streams and rivers: control and impact of eutrophication and the influence of land management. N Z J Mar Freshw Res. 2009;43(4):985–95.

    CAS  Article  Google Scholar 

  28. 28.

    Cloern JE, Foster SQ, Kleckner AE. Phytoplankton primary production in the world’s estuarine-coastal ecosystems. Biogeosciences. 2014;11(9):2477–501.

    Article  Google Scholar 

  29. 29.

    Andersson A, Brugel S, Paczkowska J, Rowe OF, Figueroa D, Kratzer S, et al. Influence of allochthonous dissolved organic matter on pelagic basal production in a northerly estuary. Estuar Coast Shelf Sci. 2018;204:225–35.

    CAS  Article  Google Scholar 

  30. 30.

    Hargrave BT, Holmer M, Newcombe CP. Towards a classification of organic enrichment in marine sediments based on biogeochemical indicators. Mar Pollut Bull. 2008;56(5):810–24.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    An S, Joye SB. Enhancement of coupled nitrification-denitrification by benthic photosynthesis in shallow estuarine sediments. Limnol Oceanogr. 2001;46(1):62–74.

    CAS  Article  Google Scholar 

  32. 32.

    Holland KR, Kleinmans A, Hussain E. State of the environment monitoring: river water quality annual report 2016. New Zealand: Auckland Council; 2018. p. 1–57. Report No.: TR2018/003

    Google Scholar 

  33. 33.

    Kennish MJ. Estuaries, anthropogenic impacts. In: Schwartz ML, editor. Encycl Coast Sci. Dordrecht: Springer Netherlands; 2005. p. 434–6.

    Google Scholar 

  34. 34.

    Kennish MJ. Environmental threats and environmental future of estuaries. Environ Conserv. Cambridge University Press. 2002;29(1):78–107.

    Article  Google Scholar 

  35. 35.

    Hume T, Gerbeaux P, Hart D, Kettles H, Neale D. A classification of New Zealand’s coastal hydrosystems. New Zealand: Ministry of the Environment; 2016. Report No.: HAM2016-062

    Google Scholar 

  36. 36.

    Onuf CP. Laguna Madre. In: Handley L, Altsman D, DeMay R, editors. Seagrass status trends North Gulf Mex 1940-2002. 1.0. Reston: U.S. Dept. of the Interior, U.S. Geological Survey; 2007. p. 29–40.

    Google Scholar 

  37. 37.

    Dame R, Alber M, Allen D, Mallin M, Montague C, Lewitus A, et al. Estuaries of the South Atlantic Coast of North America: their geographical signatures. Estuaries. 2000;23(6):793–819.

    CAS  Article  Google Scholar 

  38. 38.

    John SE, Rajimol TR, Mohan SV, Maya K, Padmalal D. Environmental degradation of a tropical estuary due to human interferences—a case study from southern Kerala, SW India. Arab J Geosci. 2017;10(16):352.

    CAS  Article  Google Scholar 

  39. 39.

    Cromwell JE. Barrier coast distribution: a world-wide survey. Abstr Vol 2nd Natl Coast Shallow Water Res Conf. 2nd ed., Los Angeles, C.A.: University Press, University of Southern California; 1971. p. 50.

  40. 40.

    Mclay C. An inventory of the status and origin of New Zealand estuarine systems. Proc N Z Ecol Soc. 1976;23:8–26.

    Google Scholar 

  41. 41.

    Davidson NC. Estuaries of Great Britain. In: Finlayson CM, Milton GR, Prentice RC, Davidson NC, editors. Wetl Book II Distrib Descr Conserv. Dordrecht: Springer Netherlands; 2018. p. 997–1009.

    Chapter  Google Scholar 

  42. 42.

    Liu T, Zhang AN, Wang J, Liu S, Jiang X, Dang C, et al. Integrated biogeography of planktonic and sedimentary bacterial communities in the Yangtze River. Microbiome. 2018;6(1):16.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Kong J, Wang Y, Warren A, Huang B, Sun P. Diversity distribution and assembly mechanisms of planktonic and benthic microeukaryote communities in intertidal zones of Southeast Fujian, China. Front Microbiol. Frontiers. 2019;10:2640.

    PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Griffiths JR, Kadin M, Nascimento FJA, Tamelander T, Törnroos A, Bonaglia S, et al. The importance of benthic–pelagic coupling for marine ecosystem functioning in a changing world. Glob Change Biol. 2017;23(6):2179–96.

    Article  Google Scholar 

  45. 45.

    van der Lee GH, Kraak MHS, Verdonschot RCM, Vonk JA, Verdonschot PFM. Oxygen drives benthic-pelagic decomposition pathways in shallow wetlands. Sci Rep. Nature Publishing Group. 2017;7:1–8.

    Article  Google Scholar 

  46. 46.

    Mansfeldt C, Achermann S, Men Y, Walser J-C, Villez K, Joss A, et al. Microbial residence time is a controlling parameter of the taxonomic composition and functional profile of microbial communities. ISME J. 2019;13(6):1589–601.

    PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Ohrel RL, Register KM. Chapter 14 Salinity. Volunt Estuary Monit Methods Man. In: The United States of America: Environmental Protection Agency & The Ocean Conservancy. 2nd ed; 2006. p. 219–28.

    Google Scholar 

  48. 48.

    Hood RR, Coles VJ, Gross TF. Marine models. In: Jørgensen SE, Fath BD, editors. Encycl Ecol. Oxford: Academic; 2008. p. 2217–25.

    Chapter  Google Scholar 

  49. 49.

    Telesh I, Schubert H, Skarlato S. Revisiting Remane’s concept: evidence for high plankton diversity and a protistan species maximum in the horohalinicum of the Baltic Sea. Mar Ecol Prog Ser. 2011;421:1–11.

    Article  Google Scholar 

  50. 50.

    Louca S, Parfrey LW, Doebeli M. Decoupling function and taxonomy in the global ocean microbiome. Science. American Association for the Advancement of Science. 2016;353(6305):1272–7.

    CAS  PubMed  Article  Google Scholar 

  51. 51.

    Uritskiy G, Tisza MJ, Gelsinger DR, Munn A, Taylor J, DiRuggiero J. Cellular life from the three domains and viruses are transcriptionally active in a hypersaline desert community. Environ Microbiol. 2020;0(7):1–17.

    CAS  Article  Google Scholar 

  52. 52.

    Louca S, Polz MF, Mazel F, Albright MBN, Huber JA, O’Connor MI, et al. Function and functional redundancy in microbial systems. Nat Ecol Evol. Nature Publishing Group. 2018;2(6):936–43.

    PubMed  Article  Google Scholar 

  53. 53.

    Martiny AC, Treseder K, Pusch G. Phylogenetic conservatism of functional traits in microorganisms. ISME J. Nature Publishing Group. 2013;7:830–8.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Glibert P, Boynton W, Flemer D, Heil C, Sharp J, Hagy J III, et al. Nutrients in estuaries: a summary report of the National Estuaries Experts Workgroup, 2005-2007. U.S.A: U.S. Environmental Protection Agency; 2010.

    Google Scholar 

  55. 55.

    Gunde-Cimerman N, Plemenitaš A, Oren A. Strategies of adaptation of microorganisms of the three domains of life to high salt concentrations. FEMS Microbiol Rev. Oxford Academic. 2018;42(3):353–75.

    CAS  PubMed  Article  Google Scholar 

  56. 56.

    Barquera B. The sodium pumping NADH:quinone oxidoreductase (Na+-NQR), a unique redox-driven ion pump. J Bioenerg Biomembr. 2014;46(4):289–98.

    CAS  PubMed  Article  Google Scholar 

  57. 57.

    Hiramatsu T, Kodama K, Kuroda T, Mizushima T, Tsuchiya T. A putative multisubunit Na+/H+ antiporter from Staphylococcus aureus. J Bacteriol. 1998;180(24):6642–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  58. 58.

    Vimont S, Berche P. NhaA, an Na+/H+ antiporter involved in environmental survival of Vibrio cholerae. J Bacteriol. 2000;182(10):2937–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Mulkidjanian AY, Dibrov P, Galperin MY. The past and present of sodium energetics: may the sodium-motive force be with you. Biochim Biophys Acta BBA Bioenerg. 1777;2008(7-8):985–92.

    CAS  Article  Google Scholar 

  60. 60.

    Lanyi JK. Salt-dependent properties of proteins from extremely halophilic bacteria. Bacteriol Rev. 1974;38(3):272–90.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  61. 61.

    Vrbka L, Vondrášek J, Jagoda-Cwiklik B, Vácha R, Jungwirth P. Quantification and rationalization of the higher affinity of sodium over potassium to protein surfaces. Proc Natl Acad Sci. National Academy of Sciences. 2006;103(42):15440–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  62. 62.

    Oren A. Bioenergetic aspects of halophilism. Microbiol Mol Biol Rev. American Society for Microbiology. 1999;63(2):334–48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  63. 63.

    Noell SE, Giovannoni SJ. SAR11 bacteria have a high affinity and multifunctional glycine betaine transporter. Environ Microbiol. 2019;21(7):2559–75.

    CAS  PubMed  Article  Google Scholar 

  64. 64.

    Antia NJ, Berland BR, Bonin DJ, Maestrini SY. Comparative evaluation of certain organic and inorganic sources of nitrogen for phototrophic growth of marine microalgae. J Mar Biol Assoc U K. 1975;55(3):519–39.

    CAS  Article  Google Scholar 

  65. 65.

    Taylor WR. Light and photosynthesis in intertidal benthic diatoms. Helgoländer Wiss Meeresunters. 1964;10(1-4):29–37.

    CAS  Article  Google Scholar 

  66. 66.

    Kjerfve B. Chapter 1 Coastal lagoons. In: Kjerfve B, editor. Coast Lagoon Process [Internet]: Elsevier; 1994. p. 1–8. Available from:

    Chapter  Google Scholar 

  67. 67.

    Keil RG, Tsamakis E, Fuh CB, Giddings JC, Hedges JI. Mineralogical and textural controls on the organic composition of coastal marine sediments: hydrodynamic separation using SPLITT-fractionation. Geochim Cosmochim Acta. 1994;58(2):879–93.

    CAS  Article  Google Scholar 

  68. 68.

    Medeiros PM, Seidel M, Ward ND, Carpenter EJ, Gomes HR, Niggemann J, et al. Fate of the Amazon River dissolved organic matter in the tropical Atlantic Ocean. Glob Biogeochem Cycles. 2015;29(5):677–90.

    CAS  Article  Google Scholar 

  69. 69.

    Prahl FG, Small LF, Eversmeyer B. Biogeochemical characterization of suspended particulate matter in the Columbia River estuary. Mar Ecol Prog Ser. 1997;160:173–84.

    CAS  Article  Google Scholar 

  70. 70.

    Seidel M, Yager PL, Ward ND, Carpenter EJ, Gomes HR, Krusche AV, et al. Molecular-level changes of dissolved organic matter along the Amazon River-to-ocean continuum. Mar Chem. 2015;177:218–31.

    CAS  Article  Google Scholar 

  71. 71.

    Eyre BD, Ferguson AJP, Webb A, Maher D, Oakes JM. Metabolism of different benthic habitats and their contribution to the carbon budget of a shallow oligotrophic sub-tropical coastal system (southern Moreton Bay, Australia). Biogeochemistry. 2011;102(1-3):87–110.

    Article  Google Scholar 

  72. 72.

    Maher DT, Eyre BD. Carbon budgets for three autotrophic Australian estuaries: implications for global estimates of the coastal air-water CO2 flux. Glob Biogeochem Cycles. 2012;26:1–18.

    Article  Google Scholar 

  73. 73.

    Dyksma S, Bischof K, Fuchs BM, Hoffmann K, Meier D, Meyerdierks A, et al. Ubiquitous Gammaproteobacteria dominate dark carbon fixation in coastal sediments. ISME J. Nature Publishing Group. 2016;10:1939–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  74. 74.

    Santoro AL, Bastviken D, Gudasz C, Tranvik L, Enrich-Prast A. Dark carbon fixation: an important process in lake sediments. PLoS One. Public Library of Science. 2013;8:e65813.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. 75.

    Berg IA. Ecological aspects of the distribution of different autotrophic CO2 fixation pathways. Appl Environ Microbiol. American Society for Microbiology. 2011;77(6):1925–36.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  76. 76.

    Lücker S, Wagner M, Maixner F, Pelletier E, Koch H, Vacherie B, et al. A Nitrospira metagenome illuminates the physiology and evolution of globally important nitrite-oxidizing bacteria. Proc Natl Acad Sci. National Academy of Sciences. 2010;107(30):13479–84.

    PubMed  PubMed Central  Article  Google Scholar 

  77. 77.

    Domingues RB, Barbosa AB, Sommer U, Galvão HM. Ammonium, nitrate and phytoplankton interactions in a freshwater tidal estuarine zone: potential effects of cultural eutrophication. Aquat Sci. 2011;73(3):331–43.

    CAS  Article  Google Scholar 

  78. 78.

    Herbert RA. Nitrogen cycling in coastal marine ecosystems. FEMS Microbiol Rev. 1999;23(5):563–90.

    CAS  PubMed  Article  Google Scholar 

  79. 79.

    Jensen MH, Lomstein E, Ssrensen J. Benthic NH4+ and NO3- flux following sedimentation of a spring phytoplankton bloom in Aarhus Bight, Denmark. Mar Ecol Prog Ser. 1990;61:87–96.

    CAS  Article  Google Scholar 

  80. 80.

    Ghylin TW, Garcia SL, Moya F, Oyserman BO, Schwientek P, Forest KT, et al. Comparative single-cell genomics reveals potential ecological niches for the freshwater acI Actinobacteria lineage. ISME J. Nature Publishing Group. 2014;8:2503–16.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  81. 81.

    Howarth RW, Marino R, Cole JJ. Nitrogen fixation in freshwater, estuarine, and marine ecosystems. 2. Biogeochemical controls1. Limnol Oceanogr. 1988;33:688–701.

    CAS  Google Scholar 

  82. 82.

    Daesch G, Mortenson LE. Effect of ammonia on the synthesis and function of the N2-fixing enzyme system in Clostridium pasteurianum. J Bacteriol. 1972;110(1):103–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  83. 83.

    Drozd JW, Tubb RS, Postgate JR. A chemostat study of the effect of fixed nitrogen sources on nitrogen fixation, membranes and free amino acids in Azotobacter chroococcum. J Gen Microbiol. 1972;73(2):221–32.

    CAS  PubMed  Article  Google Scholar 

  84. 84.

    Allen JF, Thake B, Martin WF. Nitrogenase inhibition limited oxygenation of earth’s proterozoic atmosphere. Trends Plant Sci. 2019;24(11):1022–31.

    CAS  PubMed  Article  Google Scholar 

  85. 85.

    Paerl HW, Crocker KM, Prufert LE. Limitation of N2 fixation in coastal marine waters: relative importance of molybdenum, iron, phosphorus, and organic matter availability. Limnol Oceanogr. 1987;32(3):525–36.

    CAS  Article  Google Scholar 

  86. 86.

    Raven JA. The iron and molybdenum use efficiencies of plant growth with different energy, carbon and nitrogen sources. New Phytol. 1988;109(3):279–87.

    CAS  Article  Google Scholar 

  87. 87.

    Watzer B, Forchhammer K. Cyanophycin synthesis optimizes nitrogen utilization in the unicellular cyanobacterium Synechocystis sp. PCC 6803. Appl Environ Microbiol. 2018;84(20):e01298–18.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. 88.

    Gardolinski PCFC, Worsfold PJ, McKelvie ID. Seawater induced release and transformation of organic and inorganic phosphorus from river sediments. Water Res. 2004;38(3):688–92.

    CAS  PubMed  Article  Google Scholar 

  89. 89.

    Howarth R, Chan F, Conley DJ, Garnier J, Doney SC, Marino R, et al. Coupled biogeochemical cycles: eutrophication and hypoxia in temperate estuaries and coastal marine ecosystems. Front Ecol Environ. 2011;9(1):18–26.

    Article  Google Scholar 

  90. 90.

    Lubin EA, Henry JT, Fiebig A, Crosson S, Laub MT. Identification of the PhoB regulon and role of PhoU in the phosphate starvation response of Caulobacter crescentus. de Boer P, editor. J Bacteriol. 2016;198(1):187–200.

    CAS  PubMed  Article  Google Scholar 

  91. 91.

    Santos-Beneit F. The Pho regulon: a huge regulatory network in bacteria. Front Microbiol. Frontiers. 2015;6:402.

    PubMed  PubMed Central  Article  Google Scholar 

  92. 92.

    Daims H, Lebedeva EV, Pjevac P, Han P, Herbold C, Albertsen M, et al. Complete nitrification by Nitrospira bacteria. Nature. Nature Publishing Group. 2015;528:504–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    Martens-Habbena W, Berube PM, Urakawa H, de la Torre JR, Stahl DA. Ammonia oxidation kinetics determine niche separation of nitrifying Archaea and Bacteria. Nature. Nature Publishing Group. 2009;461:976–9.

    CAS  PubMed  Article  Google Scholar 

  94. 94.

    Bouskill NJ, Eveillard D, Chien D, Jayakumar A, Ward BB. Environmental factors determining ammonia-oxidizing organism distribution and diversity in marine environments. Environ Microbiol. 2012;14(3):714–29.

    CAS  PubMed  Article  Google Scholar 

  95. 95.

    Mosier AC, Francis CA. Relative abundance and diversity of ammonia-oxidizing archaea and bacteria in the San Francisco Bay estuary. Environ Microbiol. 2008;10(11):3002–16.

    CAS  PubMed  Article  Google Scholar 

  96. 96.

    Paerl HW, Pinckney JL. A mini-review of microbial consortia: their roles in aquatic production and biogeochemical cycling. Microb Ecol. 1996;31(3):225–47.

    CAS  PubMed  Article  Google Scholar 

  97. 97.

    Hallett CS, Valesini FJ, Kilminster K, Wells NS, Eyre BD. A rapid protocol for assessing sediment condition in eutrophic estuaries. Environ Sci Process Impacts. 2019;21(6):1021–37.

    CAS  PubMed  Article  Google Scholar 

  98. 98.

    Kuypers MMM, Marchant HK, Kartal B. The microbial nitrogen-cycling network. Nat Rev Microbiol. Nature Publishing Group. 2018;16(5):263–76.

    CAS  PubMed  Article  Google Scholar 

  99. 99.

    Mohan SB, Cole JA. Chapter 7 - the dissimilatory reduction of nitrate to ammonia by anaerobic bacteria. In: Bothe H, Ferguson SJ, Newton WE, editors. Biol Nitrogen Cycle. Amsterdam: Elsevier; 2007. p. 93–106.

    Chapter  Google Scholar 

  100. 100.

    Kraft B, Tegetmeyer HE, Sharma R, Klotz MG, Ferdelman TG, Hettich RL, et al. The environmental controls that govern the end product of bacterial nitrate respiration. Science. 2014;345(6197):676–9.

    CAS  PubMed  Article  Google Scholar 

  101. 101.

    Tiedje JM, Sexstone AJ, Myrold DD, Robinson JA. Denitrification: ecological niches, competition and survival. Antonie Van Leeuwenhoek. 1983;48(6):569–83.

    Article  Google Scholar 

  102. 102.

    Wang S, Pi Y, Song Y, Jiang Y, Zhou L, Liu W, et al. Hotspot of dissimilatory nitrate reduction to ammonium (DNRA) process in freshwater sediments of riparian zones. Water Res. 2020;173:115539.

    CAS  PubMed  Article  Google Scholar 

  103. 103.

    Plummer P, Tobias C, Cady D. Nitrogen reduction pathways in estuarine sediments: influences of organic carbon and sulfide. J Geophys Res Biogeosciences. John Wiley & Sons, Ltd. 2015;120:1958–72.

    CAS  Article  Google Scholar 

  104. 104.

    Smith MW, Herfort L, Rivers AR, Simon HM. Genomic signatures for sedimentary microbial utilization of phytoplankton detritus in a fast-flowing estuary. Front Microbiol. 2019;10:2475.

    PubMed  PubMed Central  Article  Google Scholar 

  105. 105.

    Baker BJ, Lazar CS, Teske AP, Dick GJ. Genomic resolution of linkages in carbon, nitrogen, and sulfur cycling among widespread estuary sediment bacteria. Microbiome. 2015;3(1):14.

    PubMed  PubMed Central  Article  Google Scholar 

  106. 106.

    Yurkov VV, Beatty JT. Aerobic anoxygenic phototrophic bacteria. Microbiol Mol Biol Rev. 1998;62(3):695–724.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  107. 107.

    Grimm F, Franz B, Dahl C. Thiosulfate and sulfur oxidation in purple sulfur bacteria. In: Dahl C, Friedrich CG, editors. Microb Sulfur Metab. Berlin, Heidelberg: Springer; 2008. p. 101–16.

    Chapter  Google Scholar 

  108. 108.

    Béjà O, Suzuki MT, Heidelberg JF, Nelson WC, Preston CM, Hamada T, et al. Unsuspected diversity among marine aerobic anoxygenic phototrophs. Nature. 2002;415(6872):630–3.

    PubMed  Article  Google Scholar 

  109. 109.

    Koblížek M. Ecology of aerobic anoxygenic phototrophs in aquatic environments. FEMS Microbiol Rev. Oxford Academic. 2015;39(6):854–70.

    CAS  PubMed  Article  Google Scholar 

  110. 110.

    Ghosh W, Dam B. Biochemistry and molecular biology of lithotrophic sulfur oxidation by taxonomically and ecologically diverse bacteria and archaea. FEMS Microbiol Rev. 2009;33(6):999–1043.

    CAS  PubMed  Article  Google Scholar 

  111. 111.

    Sessitsch A, Weilharter A, Gerzabek MH, Kirchmann H, Kandeler E. Microbial population structures in soil particle size fractions of a long-term fertilizer field experiment. Appl Environ Microbiol. 2001;67(9):4215–24.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  112. 112.

    Rice EW, Bridgewater L, Association APH, Association AWW, Federation WE. Standard methods for the examination of water and wastewater. Washington D.C.: American Public Health Association; 2012.

    Google Scholar 

  113. 113.

    Hugerth LW, Wefer HA, Lundin S, Jakobsson HE, Lindberg M, Rodin S, et al. DegePrime, a program for degenerate primer design for broad-taxonomic-range PCR in microbial ecology studies. Appl Environ Microbiol. 2014;80(16):5116–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  114. 114.

    Walters W, Hyde ER, Berg-Lyons D, Ackermann G, Humphrey G, Parada A, et al. Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems. 2016;1:e00009–15.

    PubMed  Article  Google Scholar 

  115. 115.

    Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. Oxford Academic. 2012;28(24):3211–7.

    CAS  PubMed  Article  Google Scholar 

  116. 116.

    Bushnell B. BBMap: a fast, accurate, splice-aware aligner. 2014; Available from:

    Google Scholar 

  117. 117.

    Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010; Available from:

    Google Scholar 

  118. 118.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. Oxford Academic. 2014;30(15):2114–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  119. 119.

    Miller CS, Baker BJ, Thomas BC, Singer SW, Banfield JF. EMIRGE: reconstruction of full-length ribosomal genes from microbial community short read sequencing data. Genome Biol. 2011;12(5):R44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  120. 120.

    Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26(19):2460–1.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  121. 121.

    Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27(5):824–34.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  122. 122.

    Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20(1):257.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  123. 123.

    O’Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016;44(D1):D733–45.

    CAS  PubMed  Article  Google Scholar 

  124. 124.

    Lu J, Breitwieser FP, Thielen P, Salzberg SL. Bracken: estimating species abundance in metagenomics data. PeerJ Comput Sci. PeerJ Inc. 2017;3:e104.

    Article  Google Scholar 

  125. 125.

    Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11(1):119.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  126. 126.

    Haft DH, Selengut JD, White O. The TIGRFAMs database of protein families. Nucleic Acids Res. 2003;31(1):371–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  127. 127.

    Sonnhammer EL, Eddy SR, Birney E, Bateman A, Durbin R. Pfam: multiple sequence alignments and HMM-profiles of protein domains. Nucleic Acids Res. 1998;26(1):320–2.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  128. 128.

    Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7(10):e1002195.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  129. 129.

    Anantharaman K, Brown CT, Hug LA, Sharon I, Castelle CJ, Probst AJ, et al. Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system. Nat Commun. Nature Publishing Group. 2016;7:13219.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  130. 130.

    Fish JA, Chai B, Wang Q, Sun Y, Brown CT, Tiedje JM, et al. FunGene: the functional gene pipeline and repository. Front Microbiol. Frontiers. 2013;4:291.

    PubMed  PubMed Central  Article  Google Scholar 

  131. 131.

    Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. Public Library of Science. 2010;5:e9490.

    PubMed  PubMed Central  Article  Google Scholar 

  132. 132.

    Sayers EW, Beck J, Brister JR, Bolton EE, Canese K, Comeau DC, et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2020;48(D1):9–16.

    CAS  Article  Google Scholar 

  133. 133.

    Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. Nature Publishing Group. 2015;12(1):59–60.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  134. 134.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. Oxford Academic. 2014;30:923–30.

    CAS  Article  Google Scholar 

  135. 135.

    Wagner GP, Kin K, Lynch VJ. Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012;131(4):281–5.

    CAS  PubMed  Article  Google Scholar 

  136. 136.

    Li B, Ruotti V, Stewart RM, Thomson JA, Dewey CN. RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010;26(4):493–500.

    CAS  PubMed  Article  Google Scholar 

  137. 137.

    R Core Team. R: a language and environment for statistical computing. 2018; Available from:

    Google Scholar 

  138. 138.

    Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: community ecology package. 2019; Available from:

    Google Scholar 

  139. 139.

    Martinez Arbizu P. pairwiseAdonis: pairwise multilevel comparison using adonis. 2019; Available from:

    Google Scholar 

  140. 140.

    De’ath G, Therneau T. mvpart: multivariate partitioning. 2014; Available from:

    Google Scholar 

  141. 141.

    Segata N, Boernigen D, Tickle TL, Morgan XC, Garrett WS, Huttenhower C. Computational meta’omics for microbial community studies. Mol Syst Biol. 2013;9(1):666.

    PubMed  PubMed Central  Article  Google Scholar 

  142. 142.

    Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9(1):559.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references


Computational resources were provided by New Zealand eScience Infrastructure. We thank B Weir (Landcare Research) and C Astudillo-García, JS Boey, O Mosley and E Gios (University of Auckland) for assisting in sample collection and M Hoggard (University of Auckland) for manuscript edits.


This research was supported by Genomics Aotearoa (project 1806), a New Zealand Ministry of Business, Innovation, and Employment-funded research platform (

Author information




KMH, DW, GL and HST designed the study and collected the samples. HST and DW prepared samples for sequencing and nutrient measurements. HST and DW processed and analysed the data. HST, KMH and GL wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Kim Marie Handley.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Supplementary tables.

Additional file 2.

Supplementary results and discussion, and figures.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tee, H.S., Waite, D., Lear, G. et al. Microbial river-to-sea continuum: gradients in benthic and planktonic diversity, osmoregulation and nutrient cycling. Microbiome 9, 190 (2021).

Download citation