Single and Combined Effects of Phenanthrene and Silver Nanoparticles on Denitrification Processes in Coastal Marine Sediments

The increasing production and utilization of polycyclic aromatic hydrocarbons (PAHs) and commercial silver nanoparticles (AgNPs) have raised concerns about their potential environmental release, with coastal sediments as a substantial sink. To better understanding the effects of these contaminants on denitrification processes in coastal marine sediments, a short-term exposure simulation experiment was conducted. We investigated the effects of single and combined contamination of phenanthrene (Phe) and AgNPs on denitrification processes in a coastal marine sediment. Results showed that all contaminated treatment groups had different degrees of inhibitory effect on denitrification activity, denitrifying enzyme activity, total bacteria count and denitrifying genes. The inhibitory effect sequence of each treatment group was combined treatment > AgNPs treatment > Phe treatment. Moreover, the inhibitory effects of denitrifying genes were much larger than that of total bacteria count, indicating that the pollutants had specific toxic effects on denitrifying bacteria. The sequence of sensitivity of three reduction process to pollutants was N2O > NO2− > NO3−. All contaminated treatment groups could increase NO3−, NO2− and N2O accumulation. Furthermore, according to the linear relationship between functional gene or reductase and denitrification process, we also found that the abundance of denitrifying genes could better predict the influence of Phe and AgNPs on sediment denitrification than the denitrifying bacterial diversity. In addition, at the genus level, the community structure of nirS- and nosZ-type denitrifying bacteria changed dramatically, while changes at the phylum level were comparatively less pronounced. Single and combined contamination of Phe and AgNPs could reduce the dominance of Pseudomonas, which may lead to a potential slow-down in the degradation of Phe and inhibition of denitrification, especially the combined contamination. Overall, our study revealed that combined contamination of Phe and AgNPs could lead to an increase in NO3−, NO2− and N2O accumulation in coastal sediment, which poses a risk of eutrophication in coastal areas, exacerbates the greenhouse effect and has adverse effects on global climate change.


Introduction
Nitrogen (N) is a crucial and indispensable element for organisms and ecosystems [1].Excessive use of nitrogenous fertilizers can lead to not only an overabundance in nitrogen discharge but also eutrophication, resulting in the deterioration of water quality, decreased biodiversity and destruction of aquatic habitats [2].In the coastal ecosystem, denitrification is a significant process, where bacteria transform nitrate into nitrite and then to nitrogen gas and other gaseous forms of nitrogen, playing a vital role in controlling the fate of nitrogen in this environment [3].Functional genes such as amoA (encode ammonia onooxygenase), nxrA (nitrite oxidase), narG/napA (nitrate reductase), nirK/nirS (nitrite reductase), norB (nitric oxide reductase) and nosZ (nitrous oxide reductase) primarily drive this multi-step process [4].During denitrification, an intermediate gas, N 2 O, impacts the atmosphere by acting as a potent ozone-depleting substance and a greenhouse gas [5,6].It is crucial to note that the global warming potential of N 2 O is 310 times higher than carbon dioxide [3].
Polycyclic aromatic hydrocarbons (PAHs) belong to a group of persistent organic pollutants (POPs) and are known for their toxicity, which stems from their potential hazardous properties, recalcitrance and widespread occurrence.These compounds are generated through both natural processes and human activities, leading to their wide distribution in the environment.The PAHs content in the sediments of different types of water varies, and the PAHs content of sediments in ports, bays and estuaries is usually higher [7].
Silver nanoparticles (AgNPs) are becoming increasingly important nanomaterials with broad applications in both traditional industries and emerging fields, primarily due to their unique antibacterial properties.It was estimated that global consumption of AgNPs will exceed 360 tons per year [8,9].Some researchers predicted the growth rate of AgNPs concentration in various types of sediments in Europe, the United States and Switzerland by using a probabilistic model, and the ranges of the growth rates were 978-8593 ∆ng kg −1 y −1 , 153-1638 ∆ng kg −1 y −1 and 965-10184 ∆ng kg −1 y −1 , respectively [10].With the extensive utilization of AgNPs, a considerable amount of AgNPs infiltrates rivers through various pathways, eventually accumulating in estuarine areas.Then, the majority of AgNPs tend to settle into sediment, while less than 5% persist in seawater [11].With the increase in the production and emission of PAHs and AgNPs products, the presence of PAHs and AgNPs in coastal areas poses a substantial threat to the environment and aquatic organisms due to their persistent and bio-accumulative nature in the environment [10,12,13].
Recent studies have focused on the impacts of PAHs and AgNPs on denitrifying bacteria and denitrification processes in the soil and sediment.These studies suggested that exposure to PAHs and AgNPs, discretely, not only reduced the abundance and diversity of denitrifying genes (such as narG, nirK, nirS, nrfA, nosZ and so on) in soil and sediments but also lowered the potential denitrification activity (PDA), further impacting the community structure and function of denitrifying bacteria [8,9,[14][15][16].As a consequence, the incomplete denitrification can inhibit inorganic nitrogen decomposition and cause large N 2 O emission, which might lead to eutrophication, global warming and stratospheric ozone destruction [17].Jin et al. (2022) [18] revealed that a mixture of PAHs promoted PAHs functional genes while inhibiting denitrification functional genes.Additionally, some researchers indicated that AgNPs have adverse effects on the related bacteria and functional genes [19,20].For instance, the addition of 10 mg/L AgNPs could cause an inhibitory effect on the abundance of functional genes such as nirK and nirS, as well as the predominant bacteria Paracoccus, causing a halt in the release rate of N 2 O and N 2 , thereby increasing the accumulation of inorganic nitrogen [21].AgNPs, to a large extent, could also decrease the NO 3 − reduction rate in the soil.The degree of inhibition is closely related to the concentration of silver nanoparticles added [22].Moreover, the presence of AgNPs could inhibit the activity of nitric oxide reductase (NOR) and nitrous oxide reductase (NOS) and the abundance of corresponding encoding genes of denitrification process, norB and nosZ genes, in the sediments of hypertrophic and mesotrophic lakes [23].A single exposure to 10 mg/L of AgNPs has been proven to cause an obvious toxicity on the activity and abundance of denitrifying bacteria in the first 30 days in sediments of the Taihu Lake, China.The toxicity was reflected by the decrease in nicotinamide adenine dinucleotide (NADH) formation, electron transport system (ETS) activity, NO 2 − reductase activity (NIR), N 2 O reductase activity (NOS) and nirK gene copy number, ultimately leading to a significant decline in denitrification rate [24].However, the above studies were all focused on single effects of PAHs or AgNPs on the denitrification processes in the soil or sediment.Previ-ous studies have shown that, in addition to their inherent toxicity, nanomaterials in the environment can also interact with different types of pollutants, thereby affecting their biological effects in the environment [12].Dai et al. (2013) [25] showed that the combined contamination of multi-wall carbon nanotubes and nano-zinc oxide with phenanthrene (Phe), a PAH, could induce the generation of free radicals in the brain and liver tissues of crucian and increase the enrichment of Phe in the brain, eye and embryo of fish.Schwab et al. (2013) [26] found that carbon nanomaterials can adsorb the herbicide diuron, and their combined exposure can enhance the toxic effect of single contamination on Chlorella vulgaris.The combined exposure of nano-microplastics and PAHs could also increase the toxicity effects on Chlamydomonas reinhardtii via cell membrane damage and modifications in the physiological system.We have already explored the acute (24 h) toxicity effects of PAHs or AgNPs exposure on the denitrification processes in estuarine and marine sediments [27,28].However, research on their combined effects on denitrification remains limited.Given the widespread presence of PAHs and engineered AgNPs in the aquatic system, it is critical to investigate their combined effects on denitrification.
Pollutants tend to be more concentrated and complex in coastal areas.Nanomaterials and co-pollutants form new complexes through physical or chemical interactions, and the differences in morphological properties can change the toxic effects of nanomaterials and co-pollutants [29].The coastal marine sediment of Jiaozhou Bay (JZB), where several rivers meet the sea, was selected as the sampling site.Situated in northern China (35 • 43 -36 • 18 N, 120 • 04 -120 • 23 E), the JZB is home to Qingdao Port and has experienced significant alterations in the ecosystem due to human activities, including bank re-alignment, port construction, shipyard establishment and the development of fishing harbors.Moreover, the ecosystem of the JZB is influenced by inputs from agricultural runoff, as well as urban and industrial sewage, which contribute nutrients and waste from residential, industrial and agricultural areas, further impacting its environmental health.Therefore, the JZB is also a main gathering area of AgNPs.Meanwhile, as a typical representative of threeringed PAHs, the detection rate of Phe in the JZB and the East China Sea and other sea areas has increased in recent years [30,31].Herein, the JZB was selected as the coastal area, while Phe and AgNPs represented PAHs and nanomaterials, respectively.Through short-term exposure simulation experiments, this study aimed to evaluate the response of denitrifying capacity, the denitrification enzyme activity, the abundance of denitrifying genes and the community structure and diversity of denitrifying bacteria when exposed to single and combined contamination of Phe and AgNPs in a typical coastal marine sediment.In addition, the inhibition mechanism of Phe and AgNPs on denitrification processes and the denitrifying genes were also discussed.This study provides valuable insights into the effects of single and combined contamination of Phe and AgNPs on denitrification under more environmentally relevant conditions.

Sample Collection
Surface sediment (0-5 cm) samples were collected from the station S (36 • 07 37 N, 120 • 11 02 E) in July 2017 (Figure 1).To prevent sediment oxidation, the sediment was carefully sealed in polyethylene bags, placed in covered buckets and stored at 4 • C in a dark, cold room without freezing until further use.Before conducting sediment characterization and the following experiments, the sediment was further sieved thoroughly to remove debris larger than 1 mm [32].Seawater samples were collected from the water column above the sediment using the water sampler.The collected seawater was stored in bottles (polyethylene acid) cleaned with 10% HCl solution and rinsed with milli-Q water.Seawater samples were immediately filtered using 0.22 µm membrane filters, then stored at 4 • C. The contents of PAHs and Ag in the sediment was determined by gas chromatographymass spectrometry and inductively coupled plasma mass spectrometry, respectively.The sediment's PAHs content was found to be 0.041 ± 0.004 mg/Kg, while the Ag content was 0.019 ± 0.002 mg/Kg.The background values of PAHs and Ag in the sediment are negligible compared to the experimental addition values.
column above the sediment using the water sampler.The collected seawater was stored in bottles (polyethylene acid) cleaned with 10% HCl solution and rinsed with milli-Q water.Seawater samples were immediately filtered using 0.22 µm membrane filters, then stored at 4 °C.The contents of PAHs and Ag in the sediment was determined by gas chromatography-mass spectrometry and inductively coupled plasma mass spectrometry, respectively.The sediment's PAHs content was found to be 0.041 ± 0.004 mg/Kg, while the Ag content was 0.019 ± 0.002 mg/Kg.The background values of PAHs and Ag in the sediment are negligible compared to the experimental addition values.

Exposure to Phe and AgNPs
The exposure concentration of Phe and AgNPs in this study was determined based on the 24 h EC50 (concentration of Phe or AgNPs in sediment causing 50% inhibition of PDA) of the pollutant acting on sediment PDA alone.The EC50 of phenanthine and nanosilver acting on sediment PDA in Jiaozhou Bay was 61.19 mg kg −1 and 147.47 mg kg −1 , respectively, through the 24 h laboratory simulated culture experiment [27,28].
Before the experiment, crystalline Phe (L01921, Alfa Aesar, Lancaster, PA, USA) was dissolved in acetone and uniformly sprayed onto dry sediment, resulting in sub-samples with Phe levels of 0 (Phe-free) and 3059.5 mg kg −1 of dry sediment weight.The sub-samples were mixed thoroughly and left overnight in the dark to allow the solvent to evaporate.Meanwhile, seawater was amended with 50 mM KNO3 and 1% glucose, respectively.The AgNPs (CAS No. 7440-22-4), purchased from Sigma-Aldrich (St. Louis, MO, USA), had an average particle diameter of 100 nm (Figure 2).The original stock suspension of AgNPs were prepared by dispersing them in the aforementioned seawater to make a concentration of 147.47 mg L −1 .The suspension then underwent ultrasonication for 1 h (20 kHz and 500 W) at 4 °C to break down aggregates before being diluted to exposure concentrations.

Experiment Setups 2.2.1. Exposure to Phe and AgNPs
The exposure concentration of Phe and AgNPs in this study was determined based on the 24 h EC 50 (concentration of Phe or AgNPs in sediment causing 50% inhibition of PDA) of the pollutant acting on sediment PDA alone.The EC 50 of phenanthine and nanosilver acting on sediment PDA in Jiaozhou Bay was 61.19 mg kg −1 and 147.47 mg kg −1 , respectively, through the 24 h laboratory simulated culture experiment [27,28].
Before the experiment, crystalline Phe (L01921, Alfa Aesar, Lancaster, PA, USA) was dissolved in acetone and uniformly sprayed onto dry sediment, resulting in sub-samples with Phe levels of 0 (Phe-free) and 3059.5 mg kg −1 of dry sediment weight.The sub-samples were mixed thoroughly and left overnight in the dark to allow the solvent to evaporate.Meanwhile, seawater was amended with 50 mM KNO 3 and 1% glucose, respectively.The AgNPs (CAS No. 7440-22-4), purchased from Sigma-Aldrich (St. Louis, MO, USA), had an average particle diameter of 100 nm (Figure 2).The original stock suspension of AgNPs were prepared by dispersing them in the aforementioned seawater to make a concentration of 147.47 mg L −1 .The suspension then underwent ultrasonication for 1 h (20 kHz and 500 W) at 4 • C to break down aggregates before being diluted to exposure concentrations.To evaluate the single and combined effects of Phe and AgNPs on denitrification, a 6-day denitrification treatment was conducted in 150 mL serum bottles.The exposure ex-

Six Days Batch Exposure Experiment
To evaluate the single and combined effects of Phe and AgNPs on denitrification, a 6-day denitrification treatment was conducted in 150 mL serum bottles.The exposure experiments were divided into four treatment groups: (1) SC (sediment without any Phe and AgNPs); (2) SP (sediment contaminated with Phe (61.19 mg kg −1 wet sediment)); (3) SA (sediment contaminated with AgNPs (147.47 mg kg −1 wet sediment)); ( 4) SJ (sediment contaminated with Phe (61.19 mg kg −1 wet sediment) and AgNPs (147.47 mg kg −1 wet sediment)).Three replicates were prepared for each group.The treatment process of each experimental group is as follows: For the SC and SP treatment groups, 1 g of 0 and 3059.5 mg kg −1 dry weight sediments were mixed with 50 g wet sediment and 50 mL filtered seawater amended with 50 mM KNO 3 and 1% glucose in a 150 mL serum bottle.For the SA treatment group, 1 g of 3059.5 mg kg −1 dry sediment was mixed with 50 g wet sediment and 50 mL AgNPs stock suspensions in a 150 mL serum bottle.For the SJ treatment group, 1 g of 3059.5 mg kg −1 dry weight sediment was mixed with 50 g wet sediment and 50 mL AgNPs stock suspensions in a 150 mL serum bottle.The stock suspension was sonicated for 30 min at 4 • C before use.All serum bottles were hermetically sealed with a butyl stopper and aluminum crimp, then purged with N 2 for 15 min to remove O 2 .The bottles were incubated at a constant temperature of 25 • C with continuous shaking at 80 rpm in the dark for 6 days.After 6 days of exposure, all bottles were shaken to homogenize the samples, then 5 mL sediment slurries were collected in 10 mL sterilized centrifuge tube for DNA extraction.Additionally, 5 mL sediment slurries were mixed with 0.2 mL of 1 g L −1 HgCl and collected in 10 mL sterilized centrifuge tubes, then centrifuged with 0.22 µm filter to measure the content of NO 3 − and NO 2 − .Six 5 mL sediment slurries were taken from each parallel sample for subsequent denitrification activity culture experiments and determination of denitrifying reductase activity.Furthermore, 2 g wet sediment was transferred to the sterilized centrifuge tubes containing 5 mL of saline solution (0.22 µm filtered, 9 g L −1 NaCl, 200 mL of 0.22 µm filtered, 12.5% v/v Tween 80) and fixed with 100 µL of formaldehyde (0.22 µm filtered, 4% v/v) for total cell count analysis using flow cytometry.

Denitrification Activity
The denitrification rate was measured by the acetylene (C 2 H 2 ) inhibition technique described by Sørensen (1978) [33].Six 5 mL pre-incubated sediment slurries from each replicate of all the above-mentioned treatment groups were collected by centrifuging at 1500 rpm and washed thrice with clean seawater to remove any potential NO3 − residue.These sediment slurries were then shaken and incubated at room temperature overnight to ensure uniform microbial activity across all treatment groups.Then, the sediment was transferred into 15 mL serum bottle with 5 mL filtered seawater (amended with 300 mM KNO 3 and 2 mM glucose).All serum bottles were hermetically sealed with a butyl stopper and aluminum crimp, then purged with N 2 for 15 min to remove O 2 .Each treatment group was divided into three experimental groups, and three replicates were set for each experimental group.The first groups were filled with acetylene (C 2 H 2 , 20%, v/v) to prevent N 2 O from being further reduced to N 2 .Gas samples (6 mL) were collected from each serum bottle of the second groups after vigorous shaking, then injected into 10 mL vacuum bottles for later N 2 O analysis at time zero.All serum bottles of the first and third groups were shacked at 80 rpm in the dark for 4 h at a constant temperature of 25 • C.After 4 h, gas samples (6 mL) were collected from each serum bottle and injected into a 10 mL evacuated serum vials for later N 2 O analysis.Subsequently, the sediment slurries were treated with 0.2 mL of 1 g L −1 HgCl and collected in 10 mL sterilized centrifuge tubes, followed by centrifugation and filtration through a 0.22 µm filter for the measurement of NO 3 − and NO 2 − .Denitrification rate and N 2 O accumulation were calculated by quantifying N 2 O produced with and without C 2 H 2 , respectively.N 2 O reduction via denitrification was determined as the difference between the N 2 O produced with and without C 2 H 2 [34].The NO 3 − reduction rate was calculated as the difference between concentrations of N compounds measured at T 0 and T 4h .Similarly, the NO 2 − reduction rate was calculated by subtracting the NO 3 − reduction from the NO 2 − residues.

Determination of Denitrifying Reductase
The measurements of key enzyme activities can provide insights into the potential influences of contaminations [35,36].Six 5 mL pre-incubated sediment slurries from each replicate of all the above-mentioned treatments were transferred to 10 mL sterile centrifuge tube.Cells were harvested by centrifugation (12,000 rpm for 10 min), washed thrice with 0.1 M phosphate-buffered saline (PBS) (pH 7.4) to eliminate any possible NO 3 − residue and resuspended in the same buffer at 4 • C. Crude cell extracts were prepared by sonication of cells for 1 min at 4 • C and 20 kHz, followed by centrifugation at 12,000 rpm for 10 min at 4 • C, then the extracts were immediately used for determination of enzyme activities.
NO 3 − and NO 2 − reductase (NAR and NIR) activity: 2 mL of mixed culture medium was added to an anaerobic bottle.The mixed culture medium contained 10 mM PBS, 10 mM methyl viologen, 5 mm Na 2 S 2 O 4 and 1 mM NaNO 3 or NaNO 2 .Then, 0.1 mL of the above-mentioned supernatant was added to the bottle.The mixture was incubated at room temperature (25 • C) for 30 min under anaerobic conditions, and the NAR and NIR activity were calculated through the reduction of NO 3 − and NO 2 − [37,38].The measurement of nitrous oxide reductase (N 2 OR) activity was based on the methyl viologen method [39], and the improved process was as follows: 1.5 mL of 20 mM Tris HCl buffer (pH 7.0), 0.2 mL of 10 mM methyl viologen and 0.1 mL of the above-mentioned supernatant were added to the 3.5 mL anaerobic cuvette, which was filled with nitrogen.A total of 50 µL of 50 mM Na 2 S 2 O 4 was injected into the anaerobic cuvette by a micro syringe, which can provide reducing electrons to methyl viologen and stabilize the OD600 of the mixed reaction solution at 1.0-1.2.Then, a micro syringe was used to inject 0.1 mL of N 2 Osaturated solution into the anaerobic cuvette, and the cuvette was cultured in the ultraviolet spectrophotometer at room temperature (25 • C).In addition, the fluorescence value was measured at the same time.The change curve of OD600 was drawn to calculate the average slope, which indicates the N 2 OR activity (N 2 OR activity (U) = average slope/molar extinction coefficient of substrate) The molar extinction coefficient of substrate methyl viologen was 5.8 mM in this study.

Total Bacterial Count
The subsamples of slurries collected were treated as follows: The samples were stirred at 150 rpm for 15 min to ensure proper mixing and then subjected to sonication for 30 s at a low intensity (0.5 cycle, 20% amplitude).The sonication process was performed to disrupt any cell aggregates and ensure a more accurate analysis of bacterial cells in the cleaned suspensions.The cleaned bacterial cells in the suspension were subsequently analyzed by a BD Accuri C6 flow cytometer (Becton, Dickinson and Company, Franklin Lakes, NJ, USA) [27,40].

DNA Extraction and Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR)
The genomic DNA was extracted from 1 g wet weight of sediment using a Fast DNA SPIN kit for soil (SK8233, Sangon Biotech Co., Ltd., Shanghai, China) following the manufacturer's instructions and stored at −80 • C for later use in RT-qPCR.
The nirG, nirS and nosZ genes were selected for denitrifying bacteria.The primer pairs 1960m2F and 2050m2R-GC [39], Cd3aF and R3cdR-GC [41] and F2 and R2 [42] were used to amplify narG, nirS and nosZ gene fragments, respectively.The product length of the three genes were 110, 406 and 267 bp, respectively.Quantification of denitrifying genes was performed by a Roche LightCycler 480 (Roche, Basel, Switzerland).The 25 µL reaction mixtures contained 10 µL of SYBR Green Fast qPCR Master Mix (Roche, Basel, Switzerland), 0.4 µL of each 20 mM primer, 7.2 µL ddH 2 0 and 2 µL of template DNA.Thermal cycling conditions for nirS and nosZ genes were as follows: pre-incubation at 95 • C for 3 min, 45 cycles consisting of denaturation at 94 • C for 7 s, annealing at 57 • C for 10 s and extension at 72 • C for 30 s.After the cycling, a melting curve analysis was conducted from 60-95 • C with a heating rate of 0.11 • C s −1 and continuous fluorescence measurement.Finally, the samples were cooled down to 40 • C. The gene expression measurement was performed by RT-qPCR in triplicate, and the mean value was used for final analysis.The relative abundance of each gene was analyzed by the 2 − CT method [43].To ensure accuracy, all qPCR reactions, including the unknown samples, were performed in triplicate.In order to account for any contamination or background noise, no-template control (NTC) treatments were included in all experimental runs.

Gene Amplification, High Throughput Sequencing and Data Analysis
The V4 region of the bacterial 16S rRNA gene (about 290 nucleotides) was targeted with the primer pairs of 515f and barcoded 806r [44].Polymerase chain reaction (PCR) was performed to amplify the genes, including nirS and nosZ.The primers used for gene detection were those mentioned above.All samples were amplified in triplicate, and no-template controls were included at all steps of the process.Triplicate PCR amplicons were pooled together.The PCR reactions were carried out in 30 µL reactions with 15 µL of Phusion ® High-Fidelity PCR Master Mix (New England Biolabs Inc., Ipswich, MA, USA), 3 µL of each primer (2 µM), 10 µL of template DNA (5-10 ng) and 2 µL sterilized water.Thermal cycling consisted of initial denaturation at 98 • C for 1 min, followed by 30 cycles of denaturation at 98 • C for 10 s, annealing at 50 • C for 30 s, elongation at 72 • C for 30 s and a final step at 72 • C for 5 min.The same volume of 1× loading buffer (containing SYB green) with PCR products were mixed and electrophoresis was conducted on 2% agarose gel for detection.PCR products with a bright band were mixed in equal density ratios and purified with GeneJET Gel Extraction Kit (Thermo Scientific, Waltham, MA, USA).The purified PCR amplicons products were sequenced on the Illumina MiSeq platform (Illumina Inc., San Diego, CA, USA) at the Novogene Bioinformatics Technology Co., Ltd.(Beijing, China).
The paired-end reads from the original DNA fragments were merged by using FLASH (Version 1.2.7) [45], and UPARSE (Version 7.1), a chimera filtering approach, was employed to pick operational taxonomic units (OTUs) with 97% similarity [46] on the Bio-Linux platform.Representative sequences from each OTU were analyzed using the QIIME (Version 1.7.0)[47] software package (Quantitative Insights into Microbial Ecology), and in-house Perl scripts were used to analyze alpha (within samples) and beta (among samples) diversity analysis.We picked a representative sequence for each OTU and used the RDP classifier [48] to annotate taxonomic information for each representative sequence.The local nucleotide databases of the denitrifying genes were created by downloading DNA sequences of nirS and nosZ genes from NCBI (http://www.ncbi.nlm.nih.gov/.Accessed on 6 June 2023).
Alpha diversity was computed by rarifying the OTU table, and three metrics were calculated: Chao1, which estimates the species abundance, Observed Species, which estimates the amount of unique OTUs found in each sample and the Shannon index.Rarefaction curves were generated based on these three metrics.QIIME calculates both weighted and unweighted UniFrac, which are phylogenetic measures of beta diversity.Beta diversity was assessed using unweighted UniFrac for Principal Coordinate Analysis (PCoA) and Unweighted Pair Group Method with Arithmetic mean (UPGMA) Clustering.PCoA helps visualize complex, multidimensional data by transforming a distance matrix into a new set of orthogonal axes, demonstrating maximum variation by the first principal coordinate, second maximum by the second principal coordinate and so on.UPGMA Clustering, a hierarchical clustering method, uses average linkage to interpret the distance matrix.To further analyze the microbial diversity differences between the samples, several statistical analysis methods were applied, including T-test, MetaStat, LEfSe, Anosim and MRPP.

Effect of Single and Combined Contamination of Phe and AgNPs on Denitrification Process
The denitrification process consists of four main steps: nitrate reduction, nitrite reduction, nitric oxide reduction and nitrous oxide reduction.While the reduction rates of NO 3 − , NO 2 − , NO and N 2 O can partially reflect the denitrification ability, the denitrification rate is usually calculated by measuring the N 2 production [45].Since directly determining N 2 can be challenging, many researchers employ the acetylene inhibition technique to hinder the reduction of N 2 O to N 2 and then calculate the denitrification rate using N 2 O production as a proxy for N 2 production [14,49,50].In this paper, we used the Potential Denitrification Activity (PDA), which is measured through the acetylene inhibition technique, to represent the denitrification rate.

Nitrate and Nitrite Reduction
The control group showed the highest NO 3 − and NO 2 − reduction rates at 4.85 mg N kg −1 h −1 and 3.24 mg N kg −1 h −1 , respectively (Figure 3).The NO 3 − and NO 2 − reduction rates of the combined contaminated treatment groups were the lowest at 2.67 mg N kg −1 h −1 and 1.47 mg N kg −1 h −1 , respectively.The reduction rates of NO 3 − and NO 2 − in the Phe, AgNPs and combined contaminated treatment groups were all found to be inhibited to varying degrees compared to the control group.Specifically, the NO 3 − reduction rates in the AgNPs and combined contaminated treatment groups exhibited significant differences from the control group (ANOVA, p < 0.05), while the NO 2 − reduction rates in all three contaminated treatment groups were also significantly different from the control group (ANOVA, p < 0.05).Moreover, the inhibition rates of NO 3 − and NO 2 − reduction were both highest in the combined contaminated treatment group with 44.89% and 54.63%, respectively, followed by the AgNPs treatment group with 33.4% and 37.04%, respectively, and the Phe treatment group was the lowest with 19.60% and 24.07%, respectively.Therefore, the combined contamination of Phe and AgNPs had a stronger inhibition effect on NO 3 − and NO 2 − reduction than that of single contamination of AgNPs; meanwhile, the inhibition effect of single contamination of Phe was the weakest (ANOVA, p < 0.05).In addition, the inhibition effect of all contaminated treatment groups on NO 2 − reduction was stronger than that of NO 3 − reduction rate (ANOVA, p < 0.05).The findings indicated that the single and combined contamination of Phe and AgNPs could lead to elevated levels of NO 3 − and NO 2 − accumulation in marine environment, consequently posing an increased risk of eutrophication in the study area.
Microorganisms 2024, 12, x FOR PEER REVIEW 9 of 26 highest in the combined contaminated treatment group with 44.89% and 54.63%, respectively, followed by the AgNPs treatment group with 33.4% and 37.04%, respectively, and the Phe treatment group was the lowest with 19.60% and 24.07%, respectively.Therefore, the combined contamination of Phe and AgNPs had a stronger inhibition effect on NO3 − and NO2 − reduction than that of single contamination of AgNPs; meanwhile, the inhibition effect of single contamination of Phe was the weakest (ANOVA, p < 0.05).In addition, the inhibition effect of all contaminated treatment groups on NO2 − reduction was stronger than that of NO3 − reduction rate (ANOVA, p < 0.05).The findings indicated that the single and combined contamination of Phe and AgNPs could lead to elevated levels of NO3 − and NO2 − accumulation in marine environment, consequently posing an increased risk of eutrophication in the study area.) of different treatment groups.Note: "※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05; the same below).

Potential Denitrification Activity
In comparison to the control group, the PDA in different contaminated treatment groups was significantly inhibited to varying degrees (ANOVA, p < 0.05) (Figure 4).The highest (0.79 mg N kg −1 h −1 ) PDA was recorded in the control group, followed by the AgNPs and Phe treatment groups, where the PDA was 0.45 mg N kg −1 h −1 and 0.33 mg N kg −1 h −1 , with inhibition rates of 58.11% and 42.75%, respectively.The combined contami-

Potential Denitrification Activity
In comparison to the control group, the PDA in different contaminated treatment groups was significantly inhibited to varying degrees (ANOVA, p < 0.05) (Figure 4).The highest (0.79 mg N kg −1 h −1 ) PDA was recorded in the control group, followed by the AgNPs and Phe treatment groups, where the PDA was 0.45 mg N kg −1 h −1 and 0.33 mg N kg −1 h −1 , with inhibition rates of 58.11% and 42.75%, respectively.The combined contaminated treatment group exhibited the lowest reduction rate, measuring 0.23 mg N kg −1 h −1 , with a notable inhibition rate of 71.66%.Therefore, results indicated that the inhibition effect of the three contaminated treatment groups on PDA varied in the order combined contaminated treatment > AgNPs treatment > Phe treatment, which was consistent with the inhibition effect on NO 3 − and NO 2 − reduction rates, and the inhibition effect on PDA was stronger than that of NO 3 − and NO 2 − reduction rates (ANOVA, p < 0.05).

N2O Accumulation and Reduction
The change trend of the N2O reduction rate was consistent with the PDA, and the control group recorded the highest, in which the PDA was 0.76 mg N kg −1 h −1 , followed by the Phe and AgNPs treatment groups, and the combined contaminated treatment group recorded the lowest, with a PDA of 0.10 mg N kg −1 h −1 (Figure 5).Moreover, the N2O accumulation rate in the control group was the lowest (33.87 ug N kg −1 h −1 ), the N2O accumulation rates of the three contaminated treatment groups significantly increased compared to the control group (ANOVA, p < 0.05) and the N2O accumulation rate of the AgNPs treatment group was the highest at 0.14 mg N kg −1 h −1 .The N2O accumulation ratios ranged from 4.27% to 56.37%; the N2O accumulation ratio was the lowest in the control group, followed by the Phe and AgNPs treatment groups, and was highest in the combined treatment group.However, the change trend of the N2O reduction ratio is the opposite of that of the accumulation ratio, and the reduction ratios ranged from 43.63 to 95.73%.The above results imply that single and combined contamination of Phe and AgNPs could lead to the accumulation of N2O, an important greenhouse gas known for its stable chemical properties, long retention time in the atmosphere and average atmospheric lifetime of 100-150 years.Furthermore, N2O can be transported to the atmospheric stratosphere to deplete ozone through photochemical reactions, which in turn leads to ozone layer destruction [5,51].

N 2 O Accumulation and Reduction
The change trend of the N 2 O reduction rate was consistent with the PDA, and the control group recorded the highest, in which the PDA was 0.76 mg N kg −1 h −1 , followed by the Phe and AgNPs treatment groups, and the combined contaminated treatment group recorded the lowest, with a PDA of 0.10 mg N kg −1 h −1 (Figure 5).Moreover, the N 2 O accumulation rate in the control group was the lowest (33.87 ug N kg −1 h −1 ), the N 2 O accumulation rates of the three contaminated treatment groups significantly increased compared to the control group (ANOVA, p < 0.05) and the N 2 O accumulation rate of the AgNPs treatment group was the highest at 0.14 mg N kg −1 h −1 .The N 2 O accumulation ratios ranged from 4.27% to 56.37%; the N 2 O accumulation ratio was the lowest in the control group, followed by the Phe and AgNPs treatment groups, and was highest in the combined treatment group.However, the change trend of the N 2 O reduction ratio is the opposite of that of the accumulation ratio, and the reduction ratios ranged from 43.63 to 95.73%.The above results imply that single and combined contamination of Phe and AgNPs could lead to the accumulation of N 2 O, an important greenhouse gas known for its stable chemical properties, long retention time in the atmosphere and average atmospheric lifetime of 100-150 years.Furthermore, N 2 O can be transported to the atmospheric stratosphere to deplete ozone through photochemical reactions, which in turn leads to ozone layer destruction [5,51].
posite of that of the accumulation ratio, and the reduction ratios ranged from 43.63 to 95.73%.The above results imply that single and combined contamination of Phe and AgNPs could lead to the accumulation of N2O, an important greenhouse gas known for its stable chemical properties, long retention time in the atmosphere and average atmospheric lifetime of 100-150 years.Furthermore, N2O can be transported to the atmospheric stratosphere to deplete ozone through photochemical reactions, which in turn leads to ozone layer destruction [5,51].

Effect of Single and Combined Contamination of Phe and AgNPs on Denitrifying Reductase Activities
The short-term effects of single and combined contamination of Phe and AgNPs on NAR and NIR activities of sediments in the JZB are shown in Figure 6.Both NAR and NIR activities in the control group were the highest, at 36.45 mg N kg −1 h −1 and 32.45 mg N kg −1 h −1 , respectively.This was followed by the Phe and AgNPs treatment groups, where the NAR and NIR activities of the Phe treatment group were 36.45 mg N kg −1 h −1 and 22.53 mg N kg −1 h −1 , with inhibition rates of 21.34% and 30.57%, respectively.The inhibition rates of NAR and NIR activities of the AgNPs treatment group were 44.19% and 65.03%.The NAR and NIR activities of the combined contaminated treatment groups were the lowest, with 16.84 mg N kg −1 h −1 and 6.96 mg N kg −1 h −1 , and the inhibition rates were 53.80% and 72.76%, respectively.The NAR and NIR activities of three contaminated treatment groups showed significant differences compared to those of the control group (ANOVA, p < 0.05), and the inhibition effect on NIR activity was stronger than that on NAR activity (ANOVA, p < 0.05).

Effect of Single and Combined Contamination of Phe and AgNPs on Denitrifying Reductase Activities
The short-term effects of single and combined contamination of Phe and AgNPs on NAR and NIR activities of sediments in the JZB are shown in Figure 6.Both NAR and NIR activities in the control group were the highest, at 36.45 mg N kg −1 h −1 and 32.45 mg N kg −1 h −1 , respectively.This was followed by the Phe and AgNPs treatment groups, where the NAR and NIR activities of the Phe treatment group were 36.45 mg N kg −1 h −1 and 22.53 mg N kg −1 h −1 , with inhibition rates of 21.34% and 30.57%, respectively.The inhibition rates of NAR and NIR activities of the AgNPs treatment group were 44.19% and 65.03%.The NAR and NIR activities of the combined contaminated treatment groups were the lowest, with 16.84 mg N kg −1 h −1 and 6.96 mg N kg −1 h −1 , and the inhibition rates were 53.80% and 72.76%, respectively.The NAR and NIR activities of three contaminated treatment groups showed significant differences compared to those of the control group (ANOVA, p < 0.05), and the inhibition effect on NIR activity was stronger than that on NAR activity (ANOVA, p < 0.05).The N2OR activity of sediments in different contaminated treatment groups are shown in Figure 7.The N2OR activity in the control group was the highest, which is 2.34 U mg −1 .Followed by the Phe and AgNPs treatment groups, with N2OR activities were 1.54 U mg −1 and 1.28 U mg −1 , and the inhibition rates were 35.44% and 45.58%, respectively.The N2OR activity of the combined contaminated treatment group was the lowest with 0.59 U mg −1 , and the inhibition rate was 74.72%.The N2OR activities of three contaminated treatment groups exhibited significant differences compared to those of the control group The N 2 OR activity of sediments in different contaminated treatment groups are shown in Figure 7.The N 2 OR activity in the control group was the highest, which is 2.34 U mg −1 .Followed by the Phe and AgNPs treatment groups, with N 2 OR activities were 1.54 U mg −1 and 1.28 U mg −1 , and the inhibition rates were 35.44% and 45.58%, respectively.The N 2 OR activity of the combined contaminated treatment group was the lowest with 0.59 U mg −1 , and the inhibition rate was 74.72%.The N 2 OR activities of three contaminated treatment groups exhibited significant differences compared to those of the control group (ANOVA, p < 0.05), and the inhibition effect on N 2 OR activity was stronger than that on NAR and NIR activities (ANOVA, p < 0.05).In general, NAR, NIR and N 2 OR activities were mostly inhibited by the combined contaminated treatment of Phe and AgNPs, followed by the AgNPs treatment, and the inhibition effect of Phe treatment was the least.However, the inhibition effect of three contaminated treatment groups on denitrification capacity (PDA, NO 3 − and NO 2 − reduction rates) varied in the order PDA > NO 3 − reduction rate> NO 2 − reduction rate.This may be related to the cascade effects of nutrients, pollution or enzymatic reactions, while reductases involved in denitrification were usually sequentially induced under anaerobic conditions [52].Single and combined contamination of Phe and AgNPs had different degrees of inhibition effects on the PDA, NAR, NIR, N2OR, NO3 − and NO2 − reduction rates of sediment.Guo et al. (2011) [14] showed that the addition of PAHs (pyrene) could reduce PDA and denitrifying process-related reductase expression in a farmland soil.Some studies showed that PAHs can promote the increase in N2O accumulation while reducing the denitrification [53,54].PAHs have the ability to diminish enzyme activity by disrupting the structure and composition of cell membranes and altering the conformation of protein molecules [55,56].Moreover, the AgNPs could reduce the reduction rate of NO3 -under anaerobic conditions, and the reduction rate of NO3 − decreased significantly with the increase in AgNPs concentration [19,22].The AgNPs entering the bacteria can attract each other with negatively charged biomacromolecules and change the three-dimensional conformation of the enzyme or replace the metal ions in the enzyme through binding with important functional groups such as oxygen, nitrogen and sulfhydryl groups of the enzyme, thus affecting the stability of the enzyme and catalytic activity of the enzyme [57,58].AgNPs can also disrupt the integrity of the cell membrane, leading to the inactivation of membrane-bound enzymes and proteins [59].In addition, AgNPs in the environment can release some dissolved silver ions (Ag + ), which gradually approach and gather around bacteria under the action of electrostatic attraction; meanwhile, AgNPs and Ag + can cause the inactivation of NIR by causing a three-dimensional conformation change [60].NAR can be classified into two types-membrane-bound nitrate reductase (NAR) and periplasmic nitrate reductase (NAP)-and both NIR and N2OR are NAPs [61].Therefore, NIR and N2OR, as NAPs, may be more sensitive to AgNPs.The combination of AgNPs and hydrocarbons could reduce microbial enzyme activity in estuary-mouth sediments [62].Moreover, a previous study demonstrated that the combined contamination of AgNPs and phenol had a stronger inhibition effect on NO3 -and NO2 -reductase activities and NO3 -and NO2 -reducing ability than its single contamination [60].In this study, the combination of Phe and AgNPs had a stronger inhibitory effect on denitrification than its single combi- Single and combined contamination of Phe and AgNPs had different degrees of inhibition effects on the PDA, NAR, NIR, N 2 OR, NO 3 − and NO 2 − reduction rates of sediment.Guo et al. (2011) [14] showed that the addition of PAHs (pyrene) could reduce PDA and denitrifying process-related reductase expression in a farmland soil.Some studies showed that PAHs can promote the increase in N 2 O accumulation while reducing the denitrification [53,54].PAHs have the ability to diminish enzyme activity by disrupting the structure and composition of cell membranes and altering the conformation of protein molecules [55,56].Moreover, the AgNPs could reduce the reduction rate of NO 3 − under anaerobic conditions, and the reduction rate of NO 3 − decreased significantly with the increase in AgNPs concentration [19,22].The AgNPs entering the bacteria can attract each other with negatively charged biomacromolecules and change the three-dimensional conformation of the enzyme or replace the metal ions in the enzyme through binding with important functional groups such as oxygen, nitrogen and sulfhydryl groups of the enzyme, thus affecting the stability of the enzyme and catalytic activity of the enzyme [57,58].AgNPs can also disrupt the integrity of the cell membrane, leading to the inactivation of membranebound enzymes and proteins [59].In addition, AgNPs in the environment can release some dissolved silver ions (Ag + ), which gradually approach and gather around bacteria under the action of electrostatic attraction; meanwhile, AgNPs and Ag + can cause the inactivation of NIR by causing a three-dimensional conformation change [60].NAR can be classified into two types-membrane-bound nitrate reductase (NAR) and periplasmic nitrate reductase (NAP)-and both NIR and N 2 OR are NAPs [61].Therefore, NIR and N 2 OR, as NAPs, may be more sensitive to AgNPs.The combination of AgNPs and hydrocarbons could reduce microbial enzyme activity in estuary-mouth sediments [62].Moreover, a previous study demonstrated that the combined contamination of AgNPs and phenol had a stronger inhibition effect on NO 3 − and NO 2 − reductase activities and NO 3 − and NO 2 − reducing ability than its single contamination [60].In this study, the combination of Phe and AgNPs had a stronger inhibitory effect on denitrification than its single combination.A possible explanation is the toxic narcosis hypothesis proposed by Sikkema et al. (1994) [63], which states that the interaction between PAHs and lipophilic compounds on bacterial cell membranes can affect the integrity and permeability of cell membranes.Then, the activity of denitrifying bacteria is affected, leading to the inhibition of denitrification process.

Effect of Single and Combined Contamination of Phe and AgNPs on Total Bacterial Count and Relative Abundance of Denitrifying Genes
The short-term effects of single and combined contamination of Phe and AgNPs on the total bacteria count in the sediment of the JZB is shown in Figure 8.The total bacteria count in the control group was the highest at 1.19 × 10 9 cells g −1 .This was followed by the Phe and AgNPs treatment groups, with 9.58 × 10 8 cells g −1 and 8.78 × 10 8 cells g −1 and inhibition rates of 19.19% and 25.95%, respectively.The combined treatment group exhibited the lowest total bacterial count at 8.30 × 10 8 cells g −1 , with an inhibition rate of 30.00%.The total bacterial count in three contaminated treatment groups was significantly different from that in the control group (ANOVA, p < 0.05).Therefore, the order of the inhibitory effect of each contaminated treatment on total bacterial count was as follows: combined treatment > AgNPs treatment > Phe treatment.
Microorganisms 2024, 12, x FOR PEER REVIEW 13 of 26 The short-term effects of single and combined contamination of Phe and AgNPs on the total bacteria count in the sediment of the JZB is shown in Figure 8.The total bacteria count in the control group was the highest at 1.19 × 10 9 cells g −1 .This was followed by the Phe and AgNPs treatment groups, with 9.58 × 10 8 cells g −1 and 8.78 × 10 8 cells g −1 and inhibition rates of 19.19% and 25.95%, respectively.The combined treatment group exhibited the lowest total bacterial count at 8.30 × 10 8 cells g −1 , with an inhibition rate of 30.00%.The total bacterial count in three contaminated treatment groups was significantly different from that in the control group (ANOVA, p < 0.05).Therefore, the order of the inhibitory effect of each contaminated treatment on total bacterial count was as follows: combined treatment > AgNPs treatment > Phe treatment.narG, nirS and nosZ genes play crucial roles in the denitrification process as they are responsible for the reduction of NO3 − , NO2 − and N2O during denitrification, and their abundance directly affects the synthesis and function of related reductases.The relative abundances of the three genes in sediments of the three contaminated treatment groups are shown in Figure 9.The abundances of three denitrification functional genes in the three contaminated treatment groups were inhibited to varying degrees (ANOVA, p < 0.05).The nosZ gene exhibited the highest level of inhibition among the three denitrifying genes across all four treatment groups, with inhibition rates ranging from 53.23% to 78.19%; meanwhile, the narG gene was the least inhibited gene, and the inhibition rate ranged from 43.57% to 63.43%.The order of the inhibition effect of each contaminated treatment on the three denitrifying genes was as follows: nosZ > nirS > narG; the inhibition effect of the denitrifying genes was much larger than that of the total bacterial population (ANOVA, p < 0.05), which may indicate that contaminants may have some specific toxic effects on denitrifying bacteria rather the other bacteria.The nosZ gene of Pseudomonas stutzeri PCN-1 was more sensitive to heavy metal (Cd, Cu, Ni and Zn) contamination than the nirS gene under anaerobic conditions [64].The inhibition effect of AgNPs treatment on the nirS gene was less than that on NO2 − reductase (ANOVA, p > 0.05), but the difference was not significant; the inhibition effect of AgNPs treatment on the narG and nosZ genes was significantly greater than that on NAR and N2OR (ANOVA, p < 0.05).This shows that the impact of short-term (6 days) AgNPs exposure on the NO3 − and N2O reduction pro- narG, nirS and nosZ genes play crucial roles in the denitrification process as they are responsible for the reduction of NO 3 − , NO 2 − and N 2 O during denitrification, and their abundance directly affects the synthesis and function of related reductases.The relative abundances of the three genes in sediments of the three contaminated treatment groups are shown in Figure 9.The abundances of three denitrification functional genes in the three contaminated treatment groups were inhibited to varying degrees (ANOVA, p < 0.05).The nosZ gene exhibited the highest level of inhibition among the three denitrifying genes across all four treatment groups, with inhibition rates ranging from 53.23% to 78.19%; meanwhile, the narG gene was the least inhibited gene, and the inhibition rate ranged from 43.57% to 63.43%.The order of the inhibition effect of each contaminated treatment on the three denitrifying genes was as follows: nosZ > nirS > narG; the inhibition effect of the denitrifying genes was much larger than that of the total bacterial population (ANOVA, p < 0.05), which may indicate that contaminants may have some specific toxic effects on denitrifying bacteria rather the other bacteria.The nosZ gene of Pseudomonas stutzeri PCN-1 was more sensitive to heavy metal (Cd, Cu, Ni and Zn) contamination than the nirS gene under anaerobic conditions [64].The inhibition effect of AgNPs treatment on the nirS gene was less than that on NO 2 − reductase (ANOVA, p > 0.05), but the difference was not significant; the inhibition effect of AgNPs treatment on the narG and nosZ genes was significantly greater than that on NAR and N 2 OR (ANOVA, p < 0.05).This shows that the impact of short-term (6 days) AgNPs exposure on the NO 3 − and N 2 O reduction processes of coastal sediment is mainly through the inhibition of the related functional genes rather than the reductase activity, while the impact of NO 2 − reduction process may be through the combined inhibition of reductase activity and functional genes.Moreover, the inhibition effect of Phe treatment on narG, nirS and nosZ genes was significantly greater than that on NAR, NIR and N 2 OR (ANOVA, p < 0.05), which indicates that the impact of short-term (6 days) Phe exposure on the NO 3 − , NO 2 − and N 2 O reduction processes of coastal sediment is mainly through the inhibition of the related functional genes rather than the reductase activity.The inhibition effect of the combined treatment on the narG, nirS and nosZ genes was slightly greater than that on NAR, NIR and N 2 OR (ANOVA, p > 0.05), and the difference was not significant.This indicates that the impact of the short-term (6 days) combined Phe and AgNPs exposure on NO 3 − , NO 2 − and N 2 O reduction may occur through the combined inhibition of reductase activity and functional genes, but the toxic effect of denitrifying functional genes was slightly higher.In addition, correlation analysis of the narG, nirS and nosZ genes and corresponding reductases showed that there was a significant linear relationship between functional genes and their corresponding reductase activities (ANOVA, p < 0.05).
Microorganisms 2024, 12, x FOR PEER REVIEW 14 of 26 coastal sediment is mainly through the inhibition of the related functional genes rather than the reductase activity.The inhibition effect of the combined treatment on the narG, nirS and nosZ genes was slightly greater than that on NAR, NIR and N2OR (ANOVA, p > 0.05), and the difference was not significant.This indicates that the impact of the shortterm (6 days) combined Phe and AgNPs exposure on NO3 − , NO2 − and N2O reduction may occur through the combined inhibition of reductase activity and functional genes, but the toxic effect of denitrifying functional genes was slightly higher.In addition, correlation analysis of the narG, nirS and nosZ genes and corresponding reductases showed that there was a significant linear relationship between functional genes and their corresponding reductase activities (ANOVA, p < 0.05).Our previous research demonstrated that the inhibition effect of AgNPs on total bacterial count and denitrifying genes abundance was less than that of Phe under acute (24 h) conditions (ANOVA, p > 0.05) [27,28].However, the inhibition effect of three contaminated treatment groups on total bacterial count and relative abundance of denitrifying genes varied in the order combined treatment > AgNPs treatment > Phe treatment.The results obtained from this study suggest that the toxic effects of Phe and AgNPs on the total bacterial count and denitrifying genes abundance were different under acute and short-term contamination conditions.On one hand, it is based on the fact that the longer the AgNPs prevails in the environment, the higher the Ag + content released and the stronger the toxic effect produced.Bao et al. (2016) [65] found that after 4 and 7 days contaminated with 35 nm, 75 nm and PVP-coated AgNPs in two freshwater sediments, the total bacterial count in AgNPs treatment groups decreased to varying degrees.On the other hand, it is due to some bacteria having the ability to utilize PAHs as their carbon and energy source, thereby promoting growth [66].Some denitrifying bacteria can degrade PAHs and other hydrocarbons (aliphatic and aromatic hydrocarbons) [67,68].In addition, some studies have shown that certain heavy metals, nanomaterials and PAHs can significantly reduce the abundance of denitrifying bacteria [69][70][71][72].Our previous research demonstrated that the inhibition effect of AgNPs on total bacterial count and denitrifying genes abundance was less than that of Phe under acute (24 h) conditions (ANOVA, p > 0.05) [27,28].However, the inhibition effect of three contaminated treatment groups on total bacterial count and relative abundance of denitrifying genes varied in the order combined treatment > AgNPs treatment > Phe treatment.The results obtained from this study suggest that the toxic effects of Phe and AgNPs on the total bacterial count and denitrifying genes abundance were different under acute and short-term contamination conditions.On one hand, it is based on the fact that the longer the AgNPs prevails in the environment, the higher the Ag + content released and the stronger the toxic effect produced.Bao et al. (2016) [65] found that after 4 and 7 days contaminated with 35 nm, 75 nm and PVP-coated AgNPs in two freshwater sediments, the total bacterial count in AgNPs treatment groups decreased to varying degrees.On the other hand, it is due to some bacteria having the ability to utilize PAHs as their carbon and energy source, thereby promoting growth [66].Some denitrifying bacteria can degrade PAHs and other hydrocarbons (aliphatic and aromatic hydrocarbons) [67,68].In addition, some studies have shown that certain heavy metals, nanomaterials and PAHs can significantly reduce the abundance of denitrifying bacteria [69][70][71][72].Alpha diversity can be explained as the richness of a species community in the specific region, and it is mostly used to analyze the microbial community diversity within a sample.Diversity analysis of a single sample can reflect the richness and microbial community diversity within a sample and then evaluate the differences in the richness and diversity of the microbial communities among samples [73].In this study, by comparing the microbial diversity indices (Shannon, Simpson, Chao1 and Ace) of different contaminated treatment groups, the differences of the diversity indices of nirS and nosZ genes in coastal sediment under single and combined short-term contamination conditions of Phe and AgNPs were studied.The OTU number and four diversity indices of the nosZ gene were higher than those of the nirS gene (Table S1).Both Chao1 and Ace indices serve as indicators of community richness, where higher values suggest a greater diversity of species within the community; the Shannon index and Simpson index are used to assess community diversity, with a higher Shannon index and a lower Simpson index indicating a more diverse community [74].Single and combined contamination of Phe and AgNPs exhibited varying degrees of inhibition on the Shannon, Simpson, Chao1 and Ace indices (Table 1).Previous research showed that adding heavy metals or organic pollutants to soil can reduce microbial diversity [75,76], and the addition of AgNPs could also reduce soil microbial diversity [77].The Shannon, Chao1 and Ace indices of the nosZ gene in the three contaminated treatment groups were larger than those of the nirS gene, while the Simpson index was smaller than that of the nirS gene, indicating that the diversity and richness of the nosZ gene were greater than those of the nirS gene, which is consistent with the result of the PDA being greater than the NO 2 − reduction rate.This finding was inconsistent with the research in Taihu Lake sediments, which showed that the diversity and richness of the nirS gene were greater than those of the nosZ gene [21].Through the analysis of Shannon and Simpson indices of different contaminated treatment groups, the order of inhibitory effect on the nirS gene was Phe treatment > AgNPs treatment > combined treatment; the inhibition effect of the three contaminated treatment groups on nosZ gene varied in the order AgNPs treatment > Phe treatment > combined treatment.This was inconsistent with the order of inhibitory effects on the relative abundance of the nirS and nosZ genes in different contaminated treatment groups.Through the analysis of Chao1 and Ace indices, which reflect community richness, we found that the order of inhibitory effect on the nirS and nosZ genes was Phe treatment > combined treatment > AgNPs treatment, which was also inconsistent with the inhibition order of the relative abundance of the nirS and nosZ genes in different contaminated treatments.The results suggested that there was no obvious correlation between the diversity of denitrifying bacterial communities and denitrification function.However, it was found that the denitrifying genes abundance was strongly associated with denitrification function.The denitrifying genes abundance could better predict the effects of PAHs and AgNPs stress on sediment denitrification than diversity.Earlier studies have proven that the activities of nitrifiers and denitrifying bacteria in soil were not affected by the reduction of microbial diversity under stable environmental conditions [78,79].Guo et al. (2011) [14] also discovered that there was no substantial correlation between the PDA and the diversity of denitrifying communities, but there was a close relationship between the denitrification function and denitrifying bacterial abundance.

Beta Diversity
Beta diversity can be defined as the magnitude of community composition differences between individual samples [79,80].The weighted and unweighted UniFrac were proposed in 2005 and 2007, respectively; the former only considers the presence or absence of species, while the latter considers changes in species abundance.The UniFrac value is inversely proportional to the differences in species diversity between the samples [81,82].The weighted and unweighted UniFrac distances were calculated between the control groups and three contaminated treatment groups based on the nirS and nosZ genes.The analysis revealed that the largest distances were observed between the control group and the combined treatment group.Specifically, the distances for the nirS gene were 0.113 (weighted) and 0.551 (unweighted), while for the nosZ gene, the distances were 0.165 (weighted) and 0.631 (unweighted).This was followed by the control group and the AgNPs treatment group; the weighted and unweighted UniFrac distances between the control group and the Phe treatment group were the smallest (Figure 10).The difference between the community diversity of the two genes of the contaminated treatment groups and the control group varied in the order combined treatment > AgNPs treatment > Phe treatment.In addition, the weighted and unweighted UniFrac distances of the two genes between the combined treatment group and the Phe treatment group were higher than those of the AgNPs treatment group, indicating that the difference in community diversity of the two genes between the combined treatment group and the AgNPs treatment group was smaller than that of the Phe treatment group.correlation between the PDA and the diversity of denitrifying communities, but there was a close relationship between the denitrification function and denitrifying bacterial abundance.

Beta Diversity
Beta diversity can be defined as the magnitude of community composition differences between individual samples [79,80].The weighted and unweighted UniFrac were proposed in 2005 and 2007, respectively; the former only considers the presence or absence of species, while the latter considers changes in species abundance.The UniFrac value is inversely proportional to the differences in species diversity between the samples [81,82].The weighted and unweighted UniFrac distances were calculated between the control groups and three contaminated treatment groups based on the nirS and nosZ genes.The analysis revealed that the largest distances were observed between the control group and the combined treatment group.Specifically, the distances for the nirS gene were 0.113 (weighted) and 0.551 (unweighted), while for the nosZ gene, the distances were 0.165 (weighted) and 0.631 (unweighted).This was followed by the control group and the AgNPs treatment group; the weighted and unweighted UniFrac distances between the control group and the Phe treatment group were the smallest (Figure 10).The difference between the community diversity of the two genes of the contaminated treatment groups and the control group varied in the order combined treatment > AgNPs treatment > Phe treatment.In addition, the weighted and unweighted UniFrac distances of the two genes between the combined treatment group and the Phe treatment group were higher than those of the AgNPs treatment group, indicating that the difference in community diversity of the two genes between the combined treatment group and the AgNPs treatment group was smaller than that of the Phe treatment group.

UPGMA Cluster Analysis
UPGMA cluster analysis was conducted on the groups based on the weighted UniFrac distance matrix (Figure 11).For the nirS gene, the similarity between the AgNPs and combined treatment groups was greater than that between the Phe and combined treatment groups.For the nirS and nosZ genes, the similarity between the AgNPs and combined treatment group was large, and these two groups had little similarity with the control and Phe treatment groups, which was consistent with the results of the beta diversity.The differences demonstrated by UPGMA cluster analysis among different treatment groups might be due to the different toxic effects of Phe, AgNPs and combined contamination on denitrifying bacteria with the nirS and nosZ genes.

UPGMA Cluster Analysis
UPGMA cluster analysis was conducted on the groups based on the weighted UniFrac distance matrix (Figure 11).For the nirS gene, the similarity between the AgNPs and combined treatment groups was greater than that between the Phe and combined treatment groups.For the nirS and nosZ genes, the similarity between the AgNPs and combined treatment group was large, and these two groups had little similarity with the control and Phe treatment groups, which was consistent with the results of the beta diversity.The differences demonstrated by UPGMA cluster analysis among different treatment groups might be due to the different toxic effects of Phe, AgNPs and combined contamination on denitrifying bacteria with the nirS and nosZ genes.

UPGMA Cluster Analysis
UPGMA cluster analysis was conducted on the groups based on the weighted UniFrac distance matrix (Figure 11).For the nirS gene, the similarity between the AgNPs and combined treatment groups was greater than that between the Phe and combined treatment groups.For the nirS and nosZ genes, the similarity between the AgNPs and combined treatment group was large, and these two groups had little similarity with the control and Phe treatment groups, which was consistent with the results of the beta diversity.The differences demonstrated by UPGMA cluster analysis among different treatment groups might be due to the different toxic effects of Phe, AgNPs and combined contamination on denitrifying bacteria with the nirS and nosZ genes.The denitrifying bacterial community composition of the nirS and nosZ genes at the phylum level is shown in Figure 12.There were 13 phyla, 18 classes, 35 orders, 45 families and 59 genera of nirS-type denitrifying bacteria.The Proteobacteria was the dominant bacteria in all treatment groups.The abundance ratios of Proteobacteria in the control group, Phe treatment group, AgNPs treatment group and combined treatment group were 86.51%, 91.78%, 86.87% and 88.00%, respectively.There was little difference between the untreated groups.In addition, there were 15 phyla, 25 classes, 58 orders, 80 families and 128 genera of nosZ-type denitrifying bacteria.Proteobacteria was also the dominant bacteria in all treatment groups.The abundance ratio of Proteobacteria in the control group was the highest (82.89%), followed by the Phe and AgNPs groups (76.66% and 74.39%), and the combined treatment group was the lowest (67.48%).The dominance of Proteobacteria in nosZtype denitrifying bacteria was lower than that in nirS-type denitrifying bacteria, indicating that Proteobacteria of the nosZ gene were more sensitive to the added contaminations than those of the nirS gene.A previous study found that Pseudomonadales, Burkholderiales, Rhizobiales and Rhodospirillale belonging to Proteobacteria can degrade PAHs [14].For nirS-type denitrifying bacteria, single and combined contamination of Phe and AgNPs promoted the dominance of Actinobacteria to varying degrees but inhibited the dominance of Firmicutes and Chloroflexi to a certain extent.Furthermore, the contamination with single AgNPs and combined Phe and AgNPs could promote the dominance of Gemmatimonadetes.For nosZ-type denitrifying bacteria, single and combined contamination with Phe and AgNPs had a certain inhibitory effect on the dominance of Planctomycetes and Bacteroidetes but had a promoting effect on the dominance of Acidobacteria.The single contamination with Phe and AgNPs had a weak inhibitory effect on the dominance of Actinomycetota, while the combined contamination with Phe and AgNPs had a certain promoting effect on it.Researchers have reported that Proteobacteria, Actinobacteria and Acidobacteria are commonly dominant phyla in soils contaminated by heavy metals, and these bacteria played an important role in the bioremediation of contaminated soils [83][84][85].Denitrifying bacteria are mostly facultative anaerobic heterotrophic bacteria, including Proteobacteria, Aquificae, Deinococcus-Thermus, Actinobacteria, Bacteroidetes and Firmicutes, among which the Pseudomonas and Bacillus are ubiquitous in marine ecosystems.In addition, the denitrifying bacteria existed as Actinobacteria, Chloroflexi, Bacteroidetes and Nitrospirae [86].Fan et al. [87] showed that the main groups of denitrifying bacteria in Liaohe River estuary sediments were Proteobacteria, Firmicutes and Actinobacteria, and Proteobacteria were the dominant group; these results are similar to the conclusions of this study.The denitrifying bacterial community composition of the nirS and nosZ genes at the phylum level is shown in Figure 12.There were 13 phyla, 18 classes, 35 orders, 45 families and 59 genera of nirS-type denitrifying bacteria.The Proteobacteria was the dominant bacteria in all treatment groups.The abundance ratios of Proteobacteria in the control group, Phe treatment group, AgNPs treatment group and combined treatment group were 86.51%, 91.78%, 86.87% and 88.00%, respectively.There was little difference between the untreated groups.In addition, there were 15 phyla, 25 classes, 58 orders, 80 families and 128 genera of nosZ-type denitrifying bacteria.Proteobacteria was also the dominant bacteria in all treatment groups.The abundance ratio of Proteobacteria in the control group was the highest (82.89%), followed by the Phe and AgNPs groups (76.66% and 74.39%), and the combined treatment group was the lowest (67.48%).The dominance of Proteobacteria in nosZ-type denitrifying bacteria was lower than that in nirS-type denitrifying bacteria, indicating that Proteobacteria of the nosZ gene were more sensitive to the added contaminations than those of the nirS gene.A previous study found that Pseudomonadales, Burkholderiales, Rhizobiales and Rhodospirillale belonging to Proteobacteria can degrade PAHs [14].For nirS-type denitrifying bacteria, single and combined contamination of Phe and AgNPs promoted the dominance of Actinobacteria to varying degrees but inhibited the dominance of Firmicutes and Chloroflexi to a certain extent.Furthermore, the contamination with single AgNPs and combined Phe and AgNPs could promote the dominance of Gemmatimonadetes.For nosZ-type denitrifying bacteria, single and combined contamination with Phe and AgNPs had a certain inhibitory effect on the dominance of Planctomycetes and Bacteroidetes but had a promoting effect on the dominance of Acidobacteria.The single contamination with Phe and AgNPs had a weak inhibitory effect on the dominance of Actinomycetota, while the combined contamination with Phe and AgNPs had a certain promoting effect on it.Researchers have reported that Proteobacteria, Actinobacteria and Acidobacteria are commonly dominant phyla in soils contaminated by heavy metals, and these bacteria played an important role in the bioremediation of contaminated soils [83][84][85].Denitrifying bacteria are mostly facultative anaerobic heterotrophic bacteria, including Proteobacteria, Aquificae, Deinococcus-Thermus, Actinobacteria, Bacteroidetes and Firmicutes, among which the Pseudomonas and Bacillus are ubiquitous in marine ecosystems.In addition, the denitrifying bacteria existed as Actinobacteria, Chloroflexi, Bacteroidetes and Nitrospirae [86].Fan et al. [87] showed that the main groups of denitrifying bacteria in Liaohe River estuary sediments were Proteobacteria, Firmicutes and Actinobacteria, and Proteobacteria were the dominant group; these results are similar to the conclusions of this study.The heatmaps of the top 35 genera of nosSand nosZ-type denitrifying bacteria are shown in Figure 13.The heatmap visually represents the data, where each color block signifies the bacterial community structure at the genus level in the different treatment groups.The intensity of color reflects the degree of similarity or difference between these groups.

Community Structure of Denitrifying Bacteria after Single and Combined Contamination of Phenanthrene and AgNPs
AgNPs exposure can modify the microbial community composition of microorganisms in terrestrial, aquatic, marine and wastewater environments [88][89][90].PAHs with different rings could also alter the community structure of soil microorganisms [91].For nosS-type denitrifying bacteria, Pseudomonas was the most dominant in all treatment groups, and the abundance of Pseudomonas in the three contaminated treatment groups was obviously lower than that in the control group, but its abundance ratio increased.The abundance ratio of Pseudomonas in the control group was 27.53%, and the abundance ratios of Pseudomonas in the Phe, AgNPs and combined treatment groups were 40.45%, 37.45% and 43.28%, respectively.A previous study showed that the clonal proportion of Pseudomonas slightly increased before and after PAHs addition in a farmland soil [14].In addition, the abundance and abundance ratio of Bradyrhizobium in all contaminated treatment groups were obviously lower than those in the control group.The abundance and abundance ratio of Ruthenia in both the Phe and combined contaminated treatment groups decreased, while the abundance in the AgNPs treatment group decreased and the abundance ratio increased.Moreover, the abundance and abundance ratio of Chromohalobacter in the Phe and combined contaminated treatment groups increased noticeably.Some scholars studied nosS-type denitrifying bacteria in farmland soil and found that the addition of PAHs could significantly increase the clonal proportion of Pseudomonas (p < 0.05) and reduce the clonal proportion of Bradyrhizobium and Rhodobacteraceae [14], which was consistent with this study whereby the abundance ratio of Rugiella and Thioclava belonging to Rhodobacteraceae in the phenanthrene treatment group decreased noticeably.This was also demonstrated in a PAHs-contaminated soil from an oil field, which showed that Pseudomonas and Bradyrhizobium of nosS-type denitrifying bacteria had a high proportion of cloned sequences [15].In addition, the addition of AgNPs could also increase the abundance ratio of Pseudomonas among nosStype denitrifying bacteria in a lake sediment [21].For nosZ-type denitrifying bacteria, Pseudomonas was also the most dominant in all treatment groups, and the abundance of Pseudomonas in all three contaminated treatment groups was obviously lower than that in the control group, but its abundance ratio increased.The abundance ratio of Pseudomonas in the control group was 21.05%, and the abundance ratios of Pseudomonas in the Phe, AgNPs and combined treatment groups were 29.18%, 28.32% and 35.57%, respectively.Compared to the control group, the abundance and abundance ratio of Burkholderia and Variibacter in Phe treatment group were obviously lower, while Sphingomonas and Azospirillum were higher.This is due to the ability of Sphingomonas to degrade PAHs [92].A previous study also demonstrated that the AgNPs can increase the abundance ratio of Pseudomonas among nosZ-type denitrifying bacteria in a lake sediment [21].In short, both Phe and AgNPs contamination can increase the abundance ratio of Pseudomonas among nosSand nosZ-type denitrifying bacteria in coastal marine sediment, and the promotion effect of combined contamination with Phe and AgNPs on Pseudomonas, rather than the other bacteria, was stronger than their single contamination.AgNPs exposure can modify the microbial community composition of microorganisms in terrestrial, aquatic, marine and wastewater environments [88][89][90].PAHs with dif-

Environmental Implications
Our study indicated that single and combined contamination with Phe and AgNPs could inhibit denitrification activity, denitrifying enzyme activity, relative abundance of denitrifying genes and denitrifying bacterial diversity; meanwhile, it could also change denitrifying bacterial diversity structure and further proceed to change the nitrogentransforming rates, boosting the nitrogen loss in the water-sediment systems of a coastal area.Furthermore, the inhibition effect of combined contamination with Phe and AgNPs on denitrification and denitrifying bacteria was stronger than single contamination.
Although our research has provided valuable insights into assessing the single and combined effects of Phe and AgNPs on denitrification in coastal environment, further investigations are essential to form a comprehensive understanding.The toxicity and bioavailability of PAHs and AgNPs in the natural environment have been shown to closely correlate with their physicochemical properties, such as the ring content of PAHs and the size and coating of AgNPs [65,74,93].Additionally, the antibacterial properties of PAHs and AgNPs may be significantly influenced by environmental factors, including pH, dissolved oxygen levels and organic matter content [94,95].Moreover, different bacterial species with similar functions but inhabiting distinct environments may exhibit varying sensitivities to PAHs and AgNPs [96,97].Thus, the consequences of Phe and AgNPs on denitrification may vary among different coastal areas, and further explorations are needed.

Conclusions
Based on the experimental results, all contaminated treatment groups had different degrees of inhibitory effect on denitrification activity, denitrifying enzyme activity, total bacteria count and relative abundance of denitrifying genes.The inhibitory effect sequence of each treatment group was combined contamination > AgNPs contamination > Phe contamination.Moreover, the inhibitory effects of denitrifying genes were much larger than that of total bacteria count, indicating that the pollutants had specific toxic effects on denitrifying bacteria.The sequence of sensitivity of three reduction process to contaminations was Furthermore, it was found that denitrifying genes abundance could better predict the influence of Phe and AgNPs on sediment denitrification than the diversity of denitrifying bacteria.In addition, the community structure of nosSand nosZ-type denitrifying bacteria changed dramatically at the genus level but not greatly at the phylum level.Single and combined contamination with Phe and AgNPs can reduce the dominance of Pseudomonas, which may slow down the degradation of PAHs and inhibit the denitrification, especially the combined pollution.Thus, our results strongly indicate that both single and combined contamination of Phe and AgNPs have a notable suppressive effect on denitrification and denitrifying bacterial communities in a typical coastal marine sediment, with a combined contamination of Phe and AgNPs posing an increased threat to denitrifying bacteria in the short term.These results further imply that coastal areas like the JZB could be facing a toxicity risk from Phe and AgNPs all around the world.Compared with a single exposure, using repeated exposure of pollutants to investigate dynamic processes of microbial activity and function is more indicative for assessing the ecological risks of pollutants in the real environment.To sum up, although the impacts of Phe and AgNPs are notable, their concentrations in coastal areas remain lower than the exposure concentration in this study.The findings underscore the importance of maintaining the current environmental status.In addition,

Figure 1 .
Figure 1.The diagram of the sampling station in this study.

Figure 1 .
Figure 1.The diagram of the sampling station in this study.

Figure 2 .
Figure 2. Scanning electron microscopy picture of silver nanoparticles.

Figure 3 .
Figure 3. Reduction and inhibition rate of nitrate (NO 3 − ) and nitrite (NO 2 − ) of different treatment groups.Note: "※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05; the same below).

Figure 4 .
Figure 4. Potential denitrification activity and inhibition rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 4 .
Figure 4. Potential denitrification activity and inhibition rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 5 .
Figure 5. N 2 O accumulation and reduction rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Microorganisms 2024 , 26 Figure 5 .
Figure 5. N2O accumulation and reduction rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 6 .
Figure 6.Activity and inhibition rate of nitrate (NAR) and nitrite (NIR) reductase of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 6 .
Figure 6.Activity and inhibition rate of nitrate (NAR) and nitrite (NIR) reductase of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 7 .
Figure 7. Activity of N2O reductase (N2OR) and inhibition rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 7 .
Figure 7. Activity of N 2 O reductase (N 2 OR) and inhibition rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 8 .
Figure 8.Total bacterial count and inhibition rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 8 .
Figure 8.Total bacterial count and inhibition rate of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

Figure 9 .
Figure 9. Relative abundances of denitrifying genes of different treatment groups.

Figure 9 .
Figure 9. Relative abundances of denitrifying genes of different treatment groups."※" Represent the significant difference between group SC and the other three groups (ANOVA, p < 0.05).

3. 4 .
Community Diversity and Similarity of Denitrifying Bacteria after Single and Combined Exposure to Phe and AgNPs 3.4.1.Alpha Diversity

Figure 10 .
Figure 10.Thermal map of beta diversity index of nosS (a)-and nosZ (b)-type denitrifying bacteria between diffrernt groups.Note: in the same square, the upper and lower values represent weighted UniFrac and unweighted UniFrac distance, respectively.

3. 5 .
Community Structure of Denitrifying Bacteria after Single and Combined Contamination of Phenanthrene and AgNPs

Figure 10 .
Figure 10.Thermal map of beta diversity index of nosS (a)-and nosZ (b)-type denitrifying bacteria between diffrernt groups.Note: in the same square, the upper and lower values represent weighted UniFrac and unweighted UniFrac distance, respectively.

Figure 10 .
Figure 10.Thermal map of beta diversity index of nosS (a)-and nosZ (b)-type denitrifying bacteria between diffrernt groups.Note: in the same square, the upper and lower values represent weighted UniFrac and unweighted UniFrac distance, respectively.

3. 5 .
Community Structure of Denitrifying Bacteria after Single and Combined Contamination of Phenanthrene and AgNPs 3.5.1.Changes in Bacterial Community Structure at Phylum Level in Different Treatment Groups

Figure 13 .
Figure 13.Heatmap graph of nosS (a)-and nosZ (b)-type denitrifying bacteria at genus level.Note: bluer color blocks indicate lower abundance, while redder color blocks signify higher abundance.

3 − , NO 2 − 3 − 2 −
N 2 O > NO 2 − > NO 3 − .All contaminated treatment groups could increase NO and N 2 O accumulation.Furthermore, there was a significant linear relationship between functional gene abundance and reductase activity.The inhibitory effect of AgNPs on NO and N 2 O reduction was mainly through the inhibition of gene function, while the NO reduction was through the inhibition of both reductase activity and gene function.The inhibitory effect of Phe on NO 3 − , NO 2 − and N 2 O reduction was mainly through the inhibition of gene function.The inhibitory effect of the combined contamination on NO 3 − , NO 2 − and N 2 O reduction was through the inhibition of both reductase activity and gene function.

Table 1 .
Different microbial diversity indices of nosSand nosZ-type denitrifying bacteria in different treatment groups.