Environmental DNA metabarcoding reveals the impacts of anthropogenic pollution on multitrophic aquatic communities across an urban river of western China.

Anthropogenic activities are intensively affecting the structure and function of biological communities in river ecosystems. The effects of anthropogenic pollution on single-trophic community have been widely explored, but their effects on the structures and co-occurrence patterns of multitrophic communities remain largely unknown. In this study, we collected 13 water samples from the Neijiang River in Chengdu City of China, and identified totally 2352 bacterial, 207 algal, 204 macroinvertebrate, and 33 fish species based on the eDNA metabarcoding to systematically investigate the responses of multitrophic communities to environmental stressors. We observed significant variations in bacterial, algal, and macroinvertebrate community structures (except fish) with the pollution levels in the river. Network analyses indicated a more intensive interspecific co-occurrence pattern at high pollution level. Although taxonomic diversity of the multitrophic communities varied insignificantly, phylogenetic diversities of fish and algae showed significantly positive and negative associations with the pollution levels, respectively. We demonstrated the primary role of environmental filtering in driving the structures of bacteria, algae, and macroinvertebrates, while the fish was more controlled by dispersal limitation. Nitrogen was identified as the most important factor impacting the multitrophic community, where bacterial composition was mostly associated with NO3--N, algal spatial differentiation with TN, and macroinvertebrate and fish with NH4+-N. Further partial least-squares path model confirmed more important effect of environmental variables on the relative abundance of bacteria and algae, while macroinvertebrate and fish communities were directly driven by the algae-mediated pathway in the food web. Our study highlighted the necessity of integrated consideration of multitrophic biodiversity for riverine pollution management, and emphasized the importance of controlling nitrogen inputs targeting a healthy ecosystem.

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.

Roles of dissolved organic matter (DOM) in shaping the distribution pattern of heavy metal in the Yangtze River.

Dissolved organic matter (DOM) strongly influences the solid-liquid partitioning and migration characteristics of heavy metals, yet little is known about the metal distribution and risk with the participation of DOM in large riverine systems. This study investigated the spatiotemporal distribution of 14 heavy metals and DOM along the entire Yangtze River (over 6000 km), and highlighted the critical roles of DOM in regulating the environmental behaviors of heavy metals. Significant spatial variations of metal contents were observed, with the river source and lower reach remarkably different from the upper-middle reaches. Heavy metals in the Yangtze River were mainly from the natural sources with minor anthropogenic disturbance. We found DOM could promote the conversion of metals from solid to liquid phase and DOM with higher aromaticity showed higher metal affinities. Although low ecological risks were observed in the Yangtze River, potential risks of metal leaching warrant attention, especially for As, Cd and Sb in the middle-lower reaches with higher DOM content and aromaticity. This study established a source-to-sea investigative approach to evaluate the influences of DOM features on metal partitioning, which is crucial for the risk control and sustainable management of large rivers.

Deciphering the dynamics of metal and antibiotic resistome profiles under different metal(loid) contamination levels.

Metal(loid) contaminations pose considerable threats to ecological security and public health, yet little is known about the dynamics of metal resistance genes (MRGs) and antibiotic resistance genes (ARGs) under different metal(loid) contamination levels. Here, we provided a systematic investigation of MRGs and ARGs in three zones (Zones I, II, and III) with different metal(loid) contamination levels across an abandoned sewage reservoir. More diverse MRGs and ARGs were detected from the high-contaminated Zone I and the moderate-contaminated Zone II, while the abundant MGEs (mobile genetic elements) potentially enhanced the horizontal gene transfer potential and the resistome diversity in Zone I. Particularly, resistome hosts represented by Thiobacillus, Ramlibacter, and Dyella were prevalent in Zone II, promoting the vertical gene transfer of MRGs and ARGs. The highest health risk of ARGs was predicted for Zone I (about 7.58% and 0.48% of ARGs classified into Rank I and Rank II, respectively), followed by Zone II (2.11% and 0%) and Zone III (0% and 0%). However, the ARGs co-occurring with MRGs might exhibit low proportions and low health risks (all were Rank IV) in the three zones. Overall, these findings uncovered the dynamic responses of resistomes and their hosts to different metal(loid) contamination levels, contributing to formulating accurate management and bioremediation countermeasures for various metal(loid) contaminated environments.

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.

Different spatiotemporal dynamics, ecological drivers and assembly processes of bacterial, archaeal and fungal communities in brackish-saline groundwater.

The presence of brackish-saline groundwater (BSG) poses great harms for human health, agricultural and industrial activity. Understanding how the major environmental features in BSG determine microbiota coalescence is crucial for groundwater monitoring optimization. Based on metabarcoding analysis of 242 PCR-amplified samples, we provided the first blueprints about distinct spatiotemporal distributions, ecological drivers and assembly processes of bacterial, archaeal and fungal communities in BSG obtained from new-constructed wells at Xiong'an New Area, China. Our study demonstrated that bacterial and archaeal communities exhibited significant spatial turnovers, while fungal community displayed the most obvious seasonal variation. Environmental filtering drove bacterial compositions more than those of archaea and fungi. Total dissolved solids (TDS), one of the most critical hydrochemical factors for salinization, had a stronger effect on bacterial spatiotemporal turnover than on those of the other two taxonomic groups, while chemical oxygen demand (CODMn) was more significantly associated with prokaryotic community variations. Bacterial and archaeal taxa dominated the metacommunity network and connected closely, and TDS was mostly related to archaeal subnetwork topological features, suggesting a significant influence of TDS on species association patterns within archaea. Specific functional guilds like bacterial nitrite oxidation, anammox, and archaeal methanogenesis were enriched in lower-TDS habitats, while higher TDS favored bacterial communities involved in dark oxidation of sulfur compounds, fumarate respiration, and cellulolysis. Finally, we confirmed that bacterial and archaeal assembly processes were governed by determinism in each season, and that of fungi was more regulated by stochasticity. Higher TDS was speculated to lead bacterial assembly more deterministic and that of fungi more random. Together, these findings provided an integrate theoretical framework about the unique responses of the three life domains to brackish-saline stress, and had important implications for microbial ecological prediction in groundwater.

In-situ expressions of comammox Nitrospira along the Yangtze River.

The recent discovery of comammox Nitrospira as complete nitrifiers has significantly enriched our understanding on the nitrogen cycle, yet little is known about their metabolic transcripts in natural aquatic ecosystems. Using the genome-centric metatranscriptomics, we provided the first in-situ expression patterns of comammox Nitrospira along the Yangtze River. Our study confirmed widespread expressions of comammox Nitrospira, with the highest transcription accounting for 33.3% and 63.8% of amoA and nxrAB genes expressed in ammonia-oxidizing prokaryotes (AOPs) and Nitrospira sublineages I/II, respectively. Moreover, comammox two clades differed in nitrification, with clade A acting as the dominator to ammonia oxidation in comammox, and clade B contributing more transcripts to nitrite oxidation than to ammonia oxidation. Compared to canonical Nitrospira, comammox community had lower expressions of ammonia/nitrite transporters and nitrogen assimilatory genes, but far higher expressions in urea transport and hydrolysis, facilitating to derivation of ammonia and energy mainly through intracellular ureolytic metabolism. This suggests no need for "reciprocal-feeding" between canonical Nitrospira and AOPs in a natural river. Aerobic mixotrophy of comammox bacteria was suggested by expressions of genes coding for respiratory complexes I-V, oxidative/reductive TCA cycle, oxygen stress defenses, and transport/catabolism of simple carbohydrates and low-biosynthetic-cost amino acids. Intriguingly, significant positive correlations among expressions of ammonia monooxygenases, hydroxylamine dehydrogenase and copper-dependent nitrite reductase indicated that comammox Nitrospira had the potential of converting nitrite to nitric oxide accompanied by ammonia oxidation under low-C/N and aerobic conditions, while gene expressions in this pathway were significantly and positively associated with pH. Overall, this study illustrated novel transcriptional characteristics of comammox Nitrospira, and highlighted the necessity of reassessing their contributions to biogeochemical carbon and nitrogen cycling with perspective of in-situ meta-omics as well as culture experiments.