Low nitrous oxide concentration and spatial microbial community transition across an urban river affected by treated sewage

Urbanization significantly increases greenhouse gas emissions from a subtropical headwater stream in Southeast China

Streams are disproportionately significant contributors to increases in greenhouse gas (GHG) effluxes in river networks. In the context of global urbanization, a growing number of streams are affected by urbanization, which has been suggested to stimulate the water-air GHG emissions from fluvial systems. This study investigated the seasonal and longitudinal profiles of GHG (N2O, CH4, and CO2) concentrations of Jiuxianghe Stream, a headwater stream undergoing urbanization, and estimated its GHG diffusive fluxes and global warming potentials (GWPs) using the boundary layer method. The results showed that N2O, CH4, and CO2 concentrations in Jiuxianghe Stream were 0.45-7.19 μg L-1, 0.31-586.85 μg L-1, and 0.16-11.60 mg L-1, respectively. N2O, CH4, and CO2 concentrations in the stream showed 4.55-, 23.70-, and 7.68-fold increases from headwaters to downstream, respectively, corresponding to the forest-urban transition within the watershed. Multiple linear regression indicated that NO3--N, NH4+-N, and DOC:NO3--N accurately predicted N2O and CO2 concentrations, indicating that N nutrients were the driving factors. The Jiuxianghe Stream was a source of atmospheric GHGs with a daily GWP of 7.31 g CO2-eq m-2 d-1 on average and was significantly positively correlated with the ratio of construction land and forest in the sub-watershed. This study highlights the critical role of urbanization in amplifying GHG emissions from streams, thereby augmenting our understanding of GHG emissions from river networks. With global urbanization on the rise, streams experiencing urbanization are expected to make an unprecedentedly significant contribution to riverine GHG budgets in the future.

Hysteretic and asynchronous regime shifts of bacterial and micro-eukaryotic communities driven by nutrient loading

Nutrient pollution is pervasive in many urban rivers, while restoration measures that reduce nutrient loading but fail to improve biological communities often lack effectiveness due to the indispensable role of biota, especially multi-taxa, in enhancing ecosystem stability and function. The investigation of the response patterns of multi-taxa to the nutrient loading in urban rivers is important for the recovery of biota structure and thus ecosystem function. However, little is known about the response patterns of multi-taxa and their impact on ecosystem structure and function in urban rivers. Here, the study, from the perspective of alternative stable states theory, showed the hysteretic response of both bacterial and micro-eukaryotic communities to nutrient loading based on the field investigation and environmental DNA metabarcoding. Bistability was shown to exist in both bacterial and micro-eukaryotic communities, demonstrating that the response of microbiota to nutrient loading was a regime shifts with hysteresis. Potential analysis then indicated that the increased nutrient loading drove regime shifts in the bacterial community and the micro-eukaryotic community towards a state dominated by anaerobic bacteria and benthic Bacillariophyta, respectively. High nutrient loading was found to reduce the relative abundance of metazoan, but increase that of eukaryotic algae, which made the trophic pyramid top-lighter and bottom-heavier, probably exacerbating the degradation of ecosystem function. It should be noted that, in response to the reduced nutrient loading, the recovery threshold of micro-eukaryotic communities (nutrient loading = ∼0.5) was lower than that of bacterial communities (nutrient loading = ∼1.2), demonstrating longer hysteresis of micro-eukaryotic communities. In addition, the markedly positive correlation between the status of microbial communities and N-related enzyme activities suggested the recovery of microbial communities probably will benefit the improvement of N-cycling functionality. The obtained results provide a deep insight into the collapse and recovery trajectories of multi-trophic microbiota to the nutrient loading gradient and their impact on the N transformation potential, therefore benefiting the restoration and management of urban rivers.

Benthic clade II-type nitrous oxide reducers suppress nitrous oxide emissions in shallow lakes

Shallow lakes, recognized as hotspots for nitrogen cycling, contribute to the emission of the potent greenhouse gas nitrous oxide (N2O), but the current emission estimates for this gas have a high degree of uncertainty. However, the role of N2O-reducing bacteria (N2ORB) as N2O sinks and their contribution to N2O reduction in aquatic ecosystems in response to N2O dynamics have not been determined. Here, we investigated the N2O dynamics and microbial processes in the nitrogen cycle, which included both N2O production and consumption, in five shallow lakes spanning approximately 500 km. The investigated sites exhibited N2O oversaturation, with excess dissolved N2O concentrations (ΔN2O) ranging from 0.55 ± 0.61 to 53.17 ± 15.75 nM. Sediment-bound N2O (sN2O) was significantly positively correlated with the nitrate concentration in the overlying water (p < 0.05), suggesting that nitrate accumulation contributes to benthic N2O generation. High N2O consumption activity (RN2O) corresponded to low ΔN2O. In addition, a significant negative correlation was found between RN2O and nir/nosZ, showing that bacteria encoding nosZ contributed to N2O consumption in the benthic sediments. Redundancy analysis indicated that benthic functional genes effectively reflected the variations in RN2O and ∆N2O. qPCR analysis revealed that the clade II nosZ gene was more sensitive to ΔN2O than the clade I nosZ gene. Furthermore, four novel genera of potential nondenitrifying N2ORB were identified based on metagenome-assembled genome analysis. These genera, which are affiliated with clade II, lack genes responsible for N2O production. Collectively, benthic N2ORB, especially for clade II-type N2ORB, harnesses N2O consumption activity leading to low N2O emissions from shallow lakes. This study advances our knowledge of the role of benthic clade II-type N2ORB in regulating N2O emissions in shallow lakes.

Biogeography and impact of nitrous oxide reducers in rivers across a broad environmental gradient on emission rates

Microbial communities that reduce nitrous oxide (N2O) are divided into two clades, nosZI and nosZII. These clades significantly differ in their ecological niches and their implications for N2O emissions in terrestrial environments. However, our understanding of N2O reducers in aquatic systems is currently limited. This study investigated the relative abundance and diversity of nosZI- and nosZII-type N2O reducers in rivers and their impact on N2O emissions. Our findings revealed that stream sediments possess a high capacity for N2O reduction, surpassing N2O production under high N2O/NO3- ratio conditions. This study, along with others in freshwater systems, demonstrated that nosZI marginally dominates more often in rivers. While microbes containing either nosZI and nosZII were crucial in reducing N2O emissions, the net contribution of nosZII-containing microbes was more significant. This can be attributed to the nir gene co-occurring more frequently with the nosZI gene than with the nosZII gene. The diversity within each clade also played a role, with nosZII species being more likely to function as N2O sinks in streams with higher N2O concentrations. Overall, our findings provide a foundation for a better understanding of the biogeography of stream N2O reducers and their effects on N2O emissions.

Labile dissolved organic matter (DOM) and nitrogen inputs modified greenhouse gas dynamics: A source-to-estuary study of the Yangtze River

Although rivers are increasingly recognized as essential sources of greenhouse gases (GHG) to the atmosphere, few systematic efforts have been made to reveal the drivers of spatiotemporal variations of dissolved GHG (dGHG) in large rivers under increasing anthropogenic stress and intensified hydrological cycling. Here, through a source-to-estuary survey of the Yangtze River in March (spring) and October (autumn) of 2018, we revealed that labile dissolved organic matter (DOM) and nitrogen inputs remarkably modified the spatiotemporal distribution of dGHG. The average partial pressure of CO2 (pCO2), CH4 and N2O concentrations of all sampling sites in the Yangtze River were 1015 ± 225 μatm, and 87.5± 36.5 nmol L-1, and 20.3 ± 6.6 nmol L-1, respectively, significantly lower than the global average. In terms of longitudinal and seasonal variations, higher GHG concentrations were observed in the middle-lower reach in spring. The dominant drivers of spatiotemporal variations in dGHG were labile, protein-like DOM components and nitrogen level. Compared with the historical data of dGHG from published literature, we found a significant increase in N2O concentrations in the Yangtze River during 2004-2018, and the increasing trend was consistent with the rising riverine nitrogen concentrations. Our study emphasized the critical roles of labile DOM and nitrogen inputs in driving the spatial hotspots, seasonal variations and annual trends of dGHG. These findings can contribute to constraining the global GHG budget estimations and controls of GHG emission in large rivers in response to global change.

Elucidating the links between N2O dynamics and changes in microbial communities following saltwater intrusions

Saltwater intrusion in estuarine ecosystems alters microbial communities as well as biogeochemical cycling processes and has become a worldwide problem. However, the impact of salinity intrusion on the dynamics of nitrous oxide (N2O) and associated microbial community are understudied. Here, we conducted field microcosms in a tidal estuary during different months (December, April and August) using dialysis bags, and microbes inside the bags encountered a change in salinity in natural setting. We then compared N2O dynamics in the microcosms with that in natural water. Regardless of incubation environment, saltwater intrusion altered the dissolved N2O depending on the initial saturation rates of N2O. While the impact of saltwater intrusion on N2O dynamics was consistent across months, the dissolved N2O was higher in summer than in winter. The N-related microbial communities following saltwater intrusion were dominated by denitrifers, with fewer nitrifiers and bacterial taxa involved in dissimilatory nitrate reduction to ammonium. While denitrification was a significant driver of N2O dynamics in the studied estuary, nitrifier-involved denitrification contributed to the additional production of N2O, evidenced by the strong associations with amoA genes and the abundance of Nitrospira. Higher N2O concentrations in the field microcosms than in natural water limited N2O consumption in the former, given the lack of an association with nosZ gene abundance. The differences in the N2O dynamics observed between the microcosms and natural water could be that the latter comprised not only indigenous microbes but also those accompanied with saltwater intrusion, and that immigrants might be functionally rich individuals and able to perform N transformation in multiple pathways. Our work provides the first quantitative assessment of in situ N2O concentrations in an estuary subjected to a saltwater intrusion. The results highlight the importance of ecosystem size and microbial connectivity in the source-sink dynamics of N2O in changing environments.

Shifts in the high-resolution spatial distribution of dissolved N2O and the underlying microbial communities and processes in the Pearl River Estuary

Estuaries are significant sources of the ozone-depleting greenhouse gas N2O. However, owing to large spatial heterogeneity and discrete measurements, N2O emissions from estuaries are considerably uncertain. Microbial processes are disputed in terms of the dominant N2O production under severe human disturbance. Herein, combining real-time and high-resolution measurements with bioinformatics analysis, we accurately mapped the consecutive two-dimensional N2O distribution in the Pearl River Estuary (PRE), China, and revealed its underlying microbial mechanisms. Both the horizontal and vertical distributions of N2O concentrations varied greatly at fine scales. Supersaturated N2O concentrations (9.1 to 132.2 nmol/L) in the surface water decreased along the estuarine salinity gradient, with several emission hotspots scattering upstream. The vertical N2O distribution showed marked differences from complete mixing upstream to incomplete mixing downstream, with constant or changeable concentrations with increasing depth. Furthermore, spatially varied denitrifying and nitrifying microorganisms controlled the N2O production and distribution in the PRE, with denitrification playing the dominant role. The nirK-type and nirS-type denitrifying bacteria were the primary producers of N2O in the water and sediment columns, respectively. In addition, substrate concentration (NO3- and DOC) regulated N2O production by affecting key microbial processes, while physical influences (water-mass mixing and salt wedges) reshaped N2O distribution. With these information, a conceptual model of estuarine N2O production and distribution was constructed to generalize the possible biochemical processes under environmental constraints, which could provide insights into the N2O biogeochemical cycle and emission mitigation from a mechanistic perspective.

How nutrient loading leads to alternative stable states in microbially mediated N-cycle pathways: A new insight into bioavailable nitrogen removal in urban rivers

Excessive nutrients have disrupted pathways of microbial-mediated nitrogen (N) cycle in urban rivers and caused bioavailable N to accumulate in sediments, while remedial actions sometimes fail to recover degraded river ecosystems even when environmental quality has been improved. It is not sufficient to revert the ecosystem to its original healthy state by restoring the pre-degradation environmental conditions, as explained by alternative stable states theory. Understanding the recovery of disrupted N-cycle pathways from the perspective of alternative stable states theory can benefit effective river remediation. Previous studies have found alternative microbiota states in rivers; however, the existence and implications of alternative stable states in microbial-mediated N-cycle pathway remain unclear. Here, high-throughput sequencing and N-related enzyme activities measurement were combined in the field investigation to provide empirical evidence for the bi-stability in microbially mediated N-cycle pathways. According to the behavior of bistable ecosystems, the existence of alternative stable states in microbial-mediated N-cycle pathway have been shown, and nutrient loading, mainly total nitrogen and total phosphorus, were identified as key driver of regime shifts. In addition, potential analysis revealed that reducing nutrient loading shifted the N-cycle pathway to a desirable state characterized by high ammonification and nitrification, probably avoiding the accumulation of ammonia and organic N. It should be noted that the improvement of microbiota status can facilitate the recovery of the desirable pathway state according to the relationship between microbiota states and N-cycle pathway states. Keystone species, including Rhizobiales and Sphingomonadales, were discerned by network analysis, and the increase in their relative abundance may facilitate the improvement of microbiota status. The obtained results suggested that the nutrient reduction should be combined with microbiota management to benefit the bioavailable N removal in urban rivers, therefore providing a new insight into alleviating adverse effects of the nutrient loading on urban rivers.

Factors affecting N2O fluxes from heavy metal-contaminated mangrove soils in a subtropical estuary

A 1-year field monitoring program was carried out to observe seasonal variation in N₂O fluxes at two typical mangrove wetlands in a subtropical estuary. The soils in the island-type mangrove wetland had a higher level of heavy metal(loid) contamination and a lower level of salinity compared to the small bay-type mangrove wetland. While there was a high level of similarity in the seasonal variation pattern of N₂O fluxes between the two investigated sites with both being significantly higher in summer than in other seasons, the average of N₂O fluxes in the island-type mangrove wetland was 7.19 μg·m^-2·h^-1， which tended to be lower compared to the small bay-type mangrove wetland (15.63 μg·m^-2·h^-1). Overall, N₂O flux was closely related to soil-borne heavy metal(loid)s, showing a trend to decrease with increasing concentration of these heavy metal(loid)s. The N₂O fluxes increased with decreasing abundance of either denitrifiers or nitrifiers. But the opposite was observed for the anammox bacteria present in the soils. The anammox bacteria were more sensitive to heavy metal(loid) stress but more tolerated high salinity encountered in the investigated soils compared to the denitrifiers or nitrifiers. It appears that anammox reactions mediated by anammox bacteria played a key role in affecting the spatial variation in N₂O fluxes from the mangrove soils in the study area. And an increased level of ammonium in soils tended to promote the activity of anammox bacteria and consequently enhanced N₂O emission from the mangrove soils.

Comammox biogeography subject to anthropogenic interferences along a high-altitude river

The recent discovery of comammox Nitrospira performing complete ammonia oxidation to nitrate has overturned the long-held dogma of two-step nitrification on Earth, yet little is known about the effect of urbanization interference on their distribution. Using gene-centric metagenomics, we provided the first blueprints about comammox community, biogeography, and environmental drivers along a high-elevation (> 2000 m) river flowing through the largest city on the vulnerable Qinghai-Tibetan Plateau. Our study confirmed a wide presence and diversity of yet-uncultured comammox clade B across wet and dry seasons, with average 3.0 and 2.0 times as abundant as clade-A amoA genes in water and sediments, respectively. Species identified from freshwater and drinking water treatment plants dominated the comammox guilds (58∼100%), suggesting this plateau river shared a similar comammox assemblage with the above habitat types. Compared with the urban area harboring more abundant canonical Nitrospira identified in wastewater (average 24%), the upstream suburban reach had a smaller human population but larger proportions of comammox in ammonia-oxidizing prokaryotes (24∼72% of abundances) and Nitrospira sublineages I/II. Higher contents of nitrate and nitrite in water, and antibiotics in water and sediments, may restrain comammox niches in nitrifiers over the urban area. Further random forest analysis revealed that lincosamides and quinolones were the most important antibiotic predictors for the niche differentiations between comammox and canonical nitrifiers in water, while macrolides for those in sediments. Finally, by incubation experiments, we demonstrated higher activity contributions of benthic comammox in the suburban area (36.2∼92.8% of potential ammonia-oxidation rates) than in the urban reach, and that the contribution variation had significant negative relations with macrolides and their major components. Overall, this study highlighted that anthropogenic activities hampered the advantage of riverine complete nitrifiers over the canonical two-step ones.

Enhanced microbial nitrification-denitrification processes in a subtropical metropolitan river network

Urban rivers are hotspots of regional nitrogen (N) pollution and N transformations. Previous studies have reported that the microbial community of urban rivers was different from that of natural rivers. However, how microbial community affects N transformations in the urban rivers is still unclear. In this study, we employed N nutrients-related isotope technology (includes natural-abundance isotopes survey and isotope-labeling method) and bioinformatics methods (includes 16S rRNA high-throughput sequencing and quantitative PCR analysis) to investigate the major N transformations, microbial communities as well as functional gene abundances in a metropolitan river network. Our results suggested that the bacterial community structure in the highly urbanized rivers was characterized by higher richness, less complexity and increased abundances of nitrification and denitrifying bacterium compared to those in the suburban rivers. These differences were mainly caused by high sewage discharge and N loadings. In addition, the abundances of nitrifier gene (amoA) and denitrifier genes (nirK and nirS) were significantly higher in the highly urbanized rivers (2.36 × 10³, 7.43 × 10⁷ and 2.28 × 10⁷ copies·mL^-1) than that in the suburban rivers (0.43 × 10³, 2.18 × 10⁷ and 0.99 × 10⁷ copies·mL^-1). These changes in microbes have accelerated nitrification-denitrification processes in the highly urbanized rivers as compared to those in the suburban rivers, which was evidenced by environmental isotopes and the rates of nitrification (10.52 vs. 0.03 nmol·L^-1·h^-1) and denitrification (83.31 vs. 22.49 nmol·g^-1·h^-1). Overall, this study concluded that the excess exogenous N has significantly shaped the specific aquatic bacterial communities, which had a potential for enhancing nitrification-denitrification processes in the highly urbanized river network. This study provides a further understanding of microbial N cycling in urban river ecosystems and expands the combined application of isotopic technology and bioinformatics methods in studying biogeochemical cycling.