Front. Microbiol. Frontiers in Microbiology Front. Microbiol. 1664-302X Frontiers Media S.A. 10.3389/fmicb.2022.1062703 Microbiology Original Research Individual and combined contamination of oxytetracycline and cadmium inhibited nitrification by inhibiting ammonia oxidizers Cao Xiaoxu Zhao Wei Zhang He Lin Jitong Hu Jingying Lou Yanhong Wang Hui Yang Quangang Pan Hong * Zhuge Yuping * National Engineering Research Center for Efficient Utilization of Soil and Fertilizer Resources, College of Resources and Environment, Shandong Agricultural University, Tai’an, China

Edited by: Yongxin Lin, Fujian Normal University, China

Reviewed by: Xue Zhou, Hohai University, China; Kailou Liu, Jiangxi Institute of Red Soil, China

 *Correspondence: Hong Pan, panhong6239@163.com Yuping Zhuge, zhugeyp@sdau.edu.cn

This article was submitted to Microbiological Chemistry and Geomicrobiology, a section of the journal Frontiers in Microbiology

 01 12 2022 2022 13 1062703 06 10 2022 08 11 2022 Copyright © 2022 Cao, Zhao, Zhang, Lin, Hu, Lou, Wang, Yang, Pan and Zhuge. 2022 Cao, Zhao, Zhang, Lin, Hu, Lou, Wang, Yang, Pan and Zhuge

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

 Introduction

The large-scale development of animal husbandry and industrialization lead to more and more serious co-contamination from heavy metals and antibiotics in soils. Ecotoxic effects of residues from antibiotics and heavy metals are of increasing concern.

 Materials and Methods

In this study, oxytetracycline (OTC) and cadmium (Cd) were selected as target pollutants to evaluate the individual and combined effects on nitrification process using four different soil types sampled from North to South China through a 56-day incubation experiment.

 Results and Discussion

The results demonstrated that the contaminations of OTC and Cd, especially combined pollution had significant inhibitory effects on net nitrification rates (NNRs) as well as on AOA and AOB abundance. The toxic effects of contaminants were greatly enhanced with increasing OTC concentration. AOB was more sensitive than AOA to exogenous contaminants. And the interaction effects of OTC and Cd on ammonia oxidizers were mainly antagonistic. Furthermore, Cd contaminant (with or without OTC) had indirect effects on nitrification activity via inhibiting mineral N and AOA/AOB, while OTC alone indirectly inhibited nitrification activity by inhibiting ammonia oxidizers. The results could provide theoretical foundation for exploring the eco-environmental risks of antibiotics and heavy metals, as well as their toxic effects on nitrification processes.

 oxytetracycline cadmium nitrification ammonia-oxidizing archaea ammonia-oxidizing bacteria 2021CXGC010804 Major Science and Technology Innovation Projects Introduction

Antibiotics have become important antimicrobial drugs in recent years, and have been widely applied in medicine and agriculture for disease treatment and health protection. However, antibiotics are incompletely metabolized by humans or animals, and most antibiotics are excreted into the environment, leading to increase in antibiotic resistant bacteria (ARB) and antibiotic-resistant genes (ARGs; Knapp et al., 2010; Li et al., 2012; Wang et al., 2019b). China was a major producer and consumer of antibiotics and approximately 53,800 tons of different antibiotics were imported into environment only in 2013 (Zhang Q.-Q. et al., 2015). The worldwide consumption of antimicrobials was expected to increase from 631,511,560 to 105,596 3,605 tons from 2010 to 2030, and the consumption of antibiotics in Brazil, Russia, India, China, and South Africa would increase by 99% (Van Boeckel et al., 2015). Veterinary antibiotics were commonly used to improve animal growth efficiency. The tetracyclines such as tetracycline, oxytetracycline and chlortetracycline were widely applied to animal feed in many countries. All these contaminants tended to be excreted into the environment as part of animal excrement, such as manure and urine, and agricultural organic fertilizers (Sarmah et al., 2006; Zhang H. et al., 2015; Cao et al., 2016; Li et al., 2019; Guo et al., 2021).

Oxytetracycline (OTC) is one of the most widely used antibiotics in the world and is widely applied as feed additive in animal husbandry to stimulate animal growth (Qin et al., 2019). Besides, OTC has various functional groups, which are easily combined with soil cations or metal oxides, and thus could be adsorbed on the soil surface by cation exchange, hydrogen bond exchange, surface ligand chelation. Researchers found that the adsorption curve of OTC fit the Freundlich model, with adsorption rates close to 100%, indicating that the soil had high affinity for OTC (Tolls, 2001; Mifflin et al., 2006; Pan and Chu, 2016; He et al., 2019; Conde-Cid et al., 2020; Guo et al., 2021). In addition, some studies found that only small portions of tetracycline were ingested by animals, and nearly 60%–90% of tetracycline was released into the environment, including mud, soil and water (Liu et al., 2020). Antibiotics have direct or indirect effects on soil microbial communities, damaging the structure and function of microbial communities in various ways, resulting in the inhibition and even extinction of certain microorganisms involved in critical ecological functions (Grenni et al., 2018).

Heavy metal (loid)s contaminants in soils and plants harm food safety and threaten human health (Qin et al., 2021). Furthermore, heavy metals cannot be biodegraded in the environment, which will persist in nature and cause serious environmental challenges (Leong and Chang, 2020). Rapid industrialization and inadequate waste management have caused high incidences of toxic metal pollution in Africa and other developing countries (Okereafor et al., 2020). Similarly, the lack of effective environmental management in some nonferrous metal smelters resulted in the emissions of heavy metals in the surrounding environment and caused crop heavy metal concentrations far above national maximum allowance concentrations (Cai et al., 2019; Wang et al., 2019a).

Cadmium (Cd) was representative toxic metal that affected soil nutrient cycle and plant growth (Clemens and Ma, 2016; Rahman and Singh, 2019). Moreover, based on the results of widespread database evaluation of heavy metal concentrations in industrial and agricultural land soils, Cd, Pb and As were selected as priority control heavy metals (Yang et al., 2018). Cd was mainly existed as CdCO3 (high pH) and CdS (low pH) in soils these complexes could be transformed and be absorbed by plants, and eventually threatened human health by accumulating through the food chain (He et al., 2015; Haider et al., 2021). Huang et al. (2019) found continuous increase of soil Cd and Hg concentration in industrial areas countrywide by analyzing data published from 2005 to 2017. Hu et al. (2020) proposed that heavy metal levels exceeded national environmental background limits and increased over time in most Chinese provinces. Previous studies found that Cd was mainly derived from human activities such as mining, metallurgy, wastewater irrigation, and fertilization. Fertilizer, in particular, was found to be the major source of Cd (Shi et al., 2019). Meanwhile, the levels of heavy metals and antibiotics in organic fertilizers were significantly higher than environmental background concentrations (Zhu et al., 2013). It has been reported that metal contamination as long-term and widespread persistent selection stress could transmit antibiotic resistance genes (ARG; Baker-Austin et al., 2006; Guo et al., 2021).

Nitrogen (N) is crucial component and nutrient for the growth and development of plants. N availability depends on various chemical processes catalyzed by microorganisms, thus microorganisms play important role in N cycle (Kuypers et al., 2018). Nitrification is a microbially regulated conversion process of ammonia to nitrate via nitrite, which is essential for the global N cycle. Of which, ammonia oxidation is the first step and also rate limiting step of nitrification performed by ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA; Galloway et al., 2008; Beeckman et al., 2018; Pan et al., 2018). Recently, microorganisms capable of complete oxidation of ammonia to nitrate (comammox) have been discovered and radically challenged the conventional concept of two-step nitrification (Daims et al., 2015). The ammonia monooxygenase subunit A (amoA) gene can be used as molecular marker, which was detected in many environments and widely applied to analyze the species, abundance, community structure and evolutionary relationships of ammonia-oxidizing microorganisms (Wang J. et al., 2018; Aigle et al., 2019; Pan et al., 2021).

Previous studies have extensively reported the adverse effects of antibiotics and heavy metals on soil microorganisms and ecosystem function (Wang J. et al., 2018; Wang L. et al., 2018; Tang et al., 2020). For example, OTC inhibited protein synthesis by disrupting amino acid chain elongation at the 30S subunit of the ribosome, leading to decreased nitrification rate and bacterial growth (Kong et al., 2006). The soil microbial community function showed pronounced negative effect with increasing concentrations of OTC (Kong et al., 2006; Ma et al., 2016). Likewise, Cd exposure significantly decreased ammonification, potential nitrification activity, denitrification, as well as the amoA gene abundance and soil urease activity, which exerted a negative impact on N cycle processes (Chen et al., 2001; Bai et al., 2021; Lu et al., 2022). Zhao et al. (2020) also found that Cd inhibited soil nitrification by eliminating AOB in acidic ferrosol. Tang et al. (2020) demonstrated that addition of heavy metals alone or combined with OTC inhibited NH4+ oxidation to NO2−, then reduced the concentration of NO3−-N in sandy loam and clay loamy soils.

Antibiotics and heavy metals commonly coexisted in the soil, thus making contamination more widespread and complex (Guo et al., 2018). OTC contained several O-functional groups, which could act as potential electron donors and complexation with metal ions (Guo et al., 2021; Zhu et al., 2022). The interactions between antibiotics and heavy metals affected each other’s environmental behavior and toxicological effects (Zhao et al., 2013). In addition, OTC and Cd showed dual effects on nitrification by promoting or inhibiting the growth and activity of ammonia oxidizing microorganisms. Bai et al. (2021) found that Cd addition at 15 mg·kg−1 promoted N mineralization rates and net N mineralization accumulation, while Cd addition at 100 mg·kg−1 inhibited N transformation in constructed wetland soils. Chen et al. (2001) discovered that Cd addition at 2 and 5 mg·kg−1 could stimulate nitrification while Cd addition at 10 and 20 mg·kg−1 inhibited nitrification in alluvial soils. Moreover, Wang L. et al. (2018) found that AOA abundance was promoted at 23 mg·kg−1 of OTC addition and was inhibited at 92 and 368 mg·kg−1 of OTC addition in brown soils. Besides, Zhao et al. (2021) reported that Cd (4 and 8 mg·kg−1) directly stimulated AOB-dominated nitrification at high NH4+ levels, leading to acceleration of N transformation in clay soils. He et al. (2018) found that AOB was more sensitive than their AOA counterparts to Cu pollution in fluvo-aquic soils. Similarly, soil nitrification activity was inhibited under the threat of lead (Pb) and Cd contamination in silty clay soils. However, Wang et al. (2019) showed that single and combined addition of sulforaphane (SM2) and Cu had persistent and significant synergistic inhibitory effects on AOB and antagonistic inhibitory effects on AOA in brunisolic soils. The toxic effects of antibiotics and heavy metals on nitrification depended on a series of factors including soil organic matter (SOM), pH, clay minerals, inorganic anions, cations, as well as the contaminants concentration (Tang et al., 2020).

Therefore, four representative soils from north to south of China varied with SOM, pH, N content, etc. were selected to comprehensively explore the effects of OTC and Cd application on nitrification in different soil types. Although antibiotics and heavy metals have been frequently detected in many environments and its combined exposure significantly threaten the N cycle, the detailed ecological toxic effects of OTC and Cd on soil nitrification were currently unclear in soils. The purpose of this study was to reveal the effects of single or combined contamination of antibiotics and heavy metals on nitrification activity in different soil types. We hypothesized the following: (1) both single and combined OTC and Cd contamination could decrease net nitrification rates (NNRs) and AOA, AOB, due to immediate toxicity of OTC and Cd for broad range of microbes; (2) the inhibitory effects of OTC and Cd coexist complexes on soil nitrification processes might be greater than the effect of their single contaminants; (3) the inhibitory effects of OTC and Cd on nitrification and ammonia oxidizers might vary with soil types because of SOM and pH.

 Materials and methods Soils collection

The surface soils (0–20 cm) were sampled and were sieved (<2 mm) from Inner Mongolia (Chestnut soils), Shandong (Brown soils), Zhejiang (Paddy soils), and Hainan (Latosols) of China in June 2021 (Figure 1), then stored at 4°C prior to the incubation experiment. The vegetation type was Leymus chinensis with some Stipa grandis and Cleistogenes squarrosa in Inner Mongolia (Chestnut soils). The cropping systems in Shandong (Brown soils), Zhejiang (rice soils), and Hainan (Latosols) were winter wheat-summer corn rotation, rice-wheat rotation and rice-melon rotation, respectively. All the tested soils were uncontaminated soils, with the Cd and OTC contents below the limit of detection. Soils were sampled in an “S” pattern, with each plot being sampled at five random locations with a 5 cm diameter auger, and then mixed and combined to form a composite sample. The composite samples were transported to the laboratory at low temperature and passed through a 2 mm sieve for soil property analysis and incubation. The soil water content was measured after 24 h at 105°C. Soil pH was measured with a pH detector at a 2.5:1 (w:v) ratio of soil to distilled water. Soil organic carbon (SOC) content was analyzed by wet digestion with H2SO4-K2Cr2O7. The soil NH4+-N and NO3−-N were determined by continuous flow analyzer after extraction with 1 M KCl. The soil basic physical and chemical properties were listed in Table 1.

Sampling sites and soils in Inner Mongolia, Shandong, Zhejiang, and Hainan provinces from north to south of China. CK, control treatment without contamination; OTC0.1, OTC addition at 0.1 mg·kg−1; OTC1, OTC addition at 1 mg·kg−1; Cd5, Cd addition at 5 mg·kg−1; OTC0.1 + Cd5, OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1; OTC1 + Cd5, OTC addition at 1 mg·kg−1, and Cd addition at 5 mg·kg−1.

Physicochemical properties of sampled soils.

Sample sites Soil types pH Total nitrogen Available phosphorous Available potassium Organic matter
(g·kg−1) (mg·kg−1) (mg·kg−1) (g·kg−1)
Inner Mongolia Chestnut soils 7.17 ± 0.01a 1.59 ± 0.06b 25.19 ± 3.23a 182.33 ± 1.45a 18.46 ± 1.31b
Shandong Brown soils 6.95 ± 0.02a 1.03 ± 0.04c 18.15 ± 3.48b 104.53 ± 2.17b 16.12 ± 1.81b
Zhejiang Paddy soils 5.59 ± 0.01b 1.65 ± 0.02a 15.7 ± 2.79c 117.50 ± 1.25b 22.3 ± 2.45a
Hainan Latosols 5.70 ± 0.01b 1.23 ± 0.04c 14.50 ± 2.42c 92.67 ± 2.71b 11.77 ± 2.94c

Values represent the means ± standard error for triplicate replicates. Different letters indicate significant differences with p < 0.05 between the groups.

The net nitrification rates (NNRs) were calculated using the following equation (Khanom et al., 2021):

N = ( N O 3 − − N ) t 2 − ( N O 3 − − N ) t 1 t 2 − t 1

Where N represents the NNRs (mg·kg−1·d−1); (NO3−-N)t2 and (NO3−-N)t1 are the concentrations of NO3−-N at times 2 and 1, respectively.

The OTC with purity >95% and the CdCl2 2.5 H2O with purity >99.95% were purchased from Aladdin Reagent Company. The OTC and Cd stock solution were prepared as aqueous solutions in sterile deionized water. Other chemicals used for the incubation were of analytical grade and purchased from Sinopharm Chemical Reagent Co., Ltd.

 Incubation experiment

Fresh soil (equivalent to 50.0 g dry weight) was pre-incubated at 40% water holding capacity in a 250 ml serum bottle for 7 days at 25°C in darkness before the incubation experiment. The levels of OTC and Cd addition were referenced to the residual concentrations in soils and animal manure according to the previous research (Li et al., 2010; Afzal et al., 2019; Guo et al., 2021). Specifically, the addition of OTC at 0.1 mg·kg−1 could be conducted to explore the toxic effects in lightly contaminated soils. OTC addition at 1 mg·kg−1 could be the threshold content of OTC contamination in combined contaminated soils and could be used to analyze the dose-effect relationship (Ma et al., 2016; Tang et al., 2020). Furthermore, the toxic effects of Cd at 5 mg·kg−1 on nitrification were proved to be more typical and reliable (Kavvadias et al., 2012; Afzal et al., 2019). Microcosms were established using seven treatments (each in triplicate) including (Ι) Control treatment without contamination (CK); (II) OTC addition at 0.1 mg·kg−1 (OTC0.1); (III) OTC addition at 1 mg·kg−1 (OTC1); (IV) Cd addition at 5 mg·kg−1 (Cd5); (V) OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1 (OTC0.1 + Cd5); (VI) OTC addition at 1 mg·kg−1 and Cd addition at 5 mg·kg−1 (OTC1 + Cd5). OTC and Cd were dissolved in sterile water and added to each serum bottle at the concentrations above. The soil moisture was monitored by weight loss and adjusted to 60% of the water holding capacity. The consistent moisture content was maintained by periodically adding deionized water to compensate for any water loss during the entire incubation period. Soils were incubated at 25 ± 2°C in darkness and were collected at 0, 14, 28, and 56 days. Approximately 3 g of soil for molecular analysis was mixed and immediately stored at −80°C and the remaining soil was used to determine other properties.

 Soil DNA extraction, quantitative PCR, and high-throughput sequencing analysis

DNA was extracted from each soil sample (0.5 g) after incubation for 0 and 56 days with the Fast DNA kit (MP Biomedicals, United States) based on the manufacturer’s instruction. The extracted soil DNA samples were stored at −80°C for further analysis. The abundance of the functional genes was determined using an ABI Q5 real-time PCR system (Applied Biosystems, California, United States). The 20 μl reaction system included 10 μl of SYBR qPCR Master Mix (Vazyme), 0.4 μl primers, 1 μl of DNA template, and 8.6 μl of double-distilled water. The PCR primers were detailed in Table 2. The reaction conditions of AOA and AOB were as follows: pre-denaturation at 95°C for 30 s, denaturation at 95°C for 10 s, then 30 s at annealing temperatures (55°C for AOA-amoA gene or 57°C for AOB-amoA gene), extension at 72°C for 30 s, and the cycles were all 40. The amplification efficiency of functional gene copies was 88%–110% and R2 values were between 0.990 and 0.999.

Primers and conditions used in this study.

Primer name Primer sequence (5′-3′) Target gene References
Arch-amoAF STA ATG GTC TGG CTT AGA CG archaeal amoA gene Francis et al. (2005)
Arch-amoAR GCG GCC ATC CAT CTG TAT GT
amoA-2R CCC CTC KGS AAA GCC TTC TTC bacterial amoA gene Rotthauwe et al. (1997)

The AOA-amoA and AOB-amoA genes were sequenced on the Illumina Miseq PE300 (Illumina Inc., San Diego, CA, United States) at Shanghai Majorbio Bio-Pharm Technology Co., Ltd. platform using the same primers used for qPCR. For AOA, only one sequence direction was used for the downstream analysis, since amplified AOA amoA gene sequences with the primer set ArchamoAF/ArchamoAR had a length of 651 bp and thus did not have overlap (Chen et al., 2019; Lin et al., 2020). Raw reads were quality controlled using Trimmomatic software and spliced using FLASH. The sequences were OTU clustered based on 97% similarity using UPARSE (version 7.1 http://drive5.com/uparse/). In addition, chimeric sequences were checked and rejected using UCHIME. Raw sequence data were deposited into the NCBI Sequence Read Archive (SRA) database under accession number PRJNA874156 (AOB) and PRJNA874158 (AOA).

 Statistical analysis

Average, standard error, significance test and two-factor Analysis of Variance (ANOVA) were analyzed by SPSS 21.0. Significance was analyzed by the Duncan method with p < 0.05. Spearman’s rank correlations between NNRs and AOA, AOB amoA gene abundances were calculated in SPSS 21.0. The graphic drawing was performed on Origin 2021. A structural equation model (SEM) was constructed using IBM SPSS Amos 24.0 (AMOS IBM, United States) to evaluate the effects of contaminants on nitrification. Principal coordinates analysis (PCoA) based on Bray–Curtis distance matrices was carried out with the “vegan” package of R language (R version 3.3.1; R Core Team, 2013). The types of interactions between OTC and Cd combined contamination were analyzed by using the inhibition ration (IR).

IR = CK − Treated population CK

CK represents the amoA gene copies of the soil without OTC and/or Cd.

Treated population represents the amoA gene copies of the soil with OTC and/or Cd.

 Results Individual and combined effects of OTC and Cd on soil mineral N

To understand the effect of exogenous pollutants on nitrification process, soil mineral N under all treatments was monitored throughout the 56-day incubation period (Supplementary Figures S1, S2). The NH4+-N contents gradually and significantly decreased with incubation time regardless of OTC and Cd addition or not (Supplementary Figure S1). The individual and combined addition of OTC and Cd significantly increased NH4+-N concentrations at day 14 in all soils (Supplementary Figure S1). In addition, the NH4+-N contents showed no significant changes at days 28 and 56 compared with the control treatment in chestnut soils (Supplementary Figures S1A,D). The individual and combined addition of OTC and Cd exerted no significant effect on NH4+-N concentrations at day 28, while remarkably decreased NH4+-N concentrations at day 56 in brown soils (Supplementary Figure S1B). With respect to paddy soils, NH4+-N concentrations showed no significant change by individual and combined addition of OTC and Cd at days 28 and 56. OTC1 + Cd5 addition broadly increased NH4+-N content at day 28, OTC1 drastically increased NH4+-N content at day 56 (Supplementary Figure S1C). Similarly, NH4+-N concentrations had no significant change by individual and combined addition of OTC and Cd at days 28 and 56 in latosols. The combined addition of OTC and Cd significantly decreased NH4+-N concentrations at day 28 (Supplementary Figure S1D). The NO3−-N concentrations increased continuously and significantly during the 56-day incubation period (Supplementary Figure S2). Application of OTC and Cd significantly decreased the NO3−-N concentrations. The combined application of OTC and Cd showed a greater inhibitory effect on NO3−-N concentrations than OTC or Cd alone (Supplementary Figure S2). The inhibitory effects were significant at days 14 and 28 in northern alkaline soils (chestnut and brown soils) and were only remarkable at day 14 in southern acidic soils (paddy and latosols soils; Supplementary Figures S2A–D).

NNRs were calculated to further reveal the effects of individual and combined pollution of OTC and Cd on nitrification activity (Figure 2). The pollution of OTC and Cd clearly decreased the NNRs compared with CK treatment in all soils (except brown soils; Figures 2A–D). The inhibitory effects were distinct at days 14 and 28 in northern alkaline soils (chestnut and brown soils) and were only remarkable at day 14 in southern acidic soils (paddy and latosols soils) compared with the control treatment. Additionally, the inhibitory effects were increased with increasing OTC concentration. Combined pollution of both OTC and Cd further stimulated the inhibitory effects.

Dynamics of NNRs in chestnut soils (A), brown soils (B), paddy soils (C), and latosols (D) under CK, OTC0.1, OTC1, Cd5, OTC0.1 + Cd5, and OTC1 + Cd5 treatments over the 56-day incubation. The vertical bars indicated the standard errors of the mean of triplicate samples. Different lowercase letters above the error bars indicated significant differences among treatments (p < 0.05). Different uppercase letters above the error bars indicated significant differences among sampling time points under the same treatment (p < 0.05). CK, control treatment without contamination; OTC0.1, OTC addition at 0.1 mg·kg−1; OTC1, OTC addition at 1 mg·kg−1; Cd5, Cd addition at 5 mg·kg−1; OTC0.1 + Cd5, OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1; OTC1 + Cd5, OTC addition at 1 mg·kg−1 and Cd addition at 5 mg·kg−1.

 Individual and combined effects of OTC and Cd contamination on the <italic>amoA</italic> gene abundance

Generally, the copy numbers of AOA and AOB amoA genes increased over the whole incubation time in OTC and Cd-treated soils. The abundances of AOA amoA ranged from 2.96 × 107 to 5.31 × 107, 3.95 × 105 to 5.86 × 105, 4.70× 106 to 1.18 × 107, and 2.67 × 106 to 5.78 × 106 copies g−1 d.w.s. (dry weight soil) in the chestnut, brown, paddy, and latosols soils, respectively, from day 0 to day 56 (Figure 3). The AOA amoA gene copy numbers were significantly higher in chestnut soils than in other soils. Overall, the dynamics of AOA abundance resembled that of NNRs. Individual and combined pollution of OTC and Cd significantly decreased AOA amoA gene copy numbers. And the inhibitory effects were increased with the increasing of OTC concentration, especially in combined polluted soils.

Abundance of AOA amoA gene in chestnut soils (A), brown soils (B), paddy soils (C), and latosols (D) under CK, OTC0.1, OTC1, Cd5, OTC0.1 + Cd5, and OTC1 + Cd5 treatments at day 0 and day 56. The vertical bars indicated the standard errors of the mean of triplicate samples. Different lowercase letters above the error bars indicated significant differences among treatments (p < 0.05). Different uppercase letters above the error bars indicated significant differences among sampling time points under the same treatment (p < 0.05). CK, control treatment without contamination; OTC0.1, OTC addition at 0.1 mg·kg−1; OTC1, OTC addition at 1 mg·kg−1; Cd5, Cd addition at 5 mg·kg−1; OTC0.1 + Cd5, OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1; OTC1 + Cd5, OTC addition at 1 mg·kg−1, and Cd addition at 5 mg·kg−1.

The copy numbers of the AOB amoA gene in the chestnut, brown, paddy and latosols soils ranged from 1.65 × 104 to 1.29 × 105, 2.37 × 104 to 5.97 × 104, 1.69 × 104 to 4.95 × 104, 1.24 × 104 to 1.42 × 105 copies g−1 d.w.s. from day 0 to day 56, respectively (Figure 4). The dynamics of AOB abundance followed the same trend with AOA amoA gene copy numbers. The inhibitory effects of pollution on AOB abundance followed the trend OTC1 + Cd5 > OTC0.1 + Cd5 > Cd5 > OTC1 > OTC0.1.

Abundance of AOB amoA gene in chestnut soils (A), brown soils (B), paddy soils (C), and latosols (D) under CK, OTC0.1, OTC1, Cd5, OTC0.1 + Cd5, and OTC1 + Cd5 treatments at day 0 and day 56. The vertical bars indicated the standard errors of the mean of triplicate samples. Different lowercase letters above the error bars indicated significant differences among treatments (p < 0.05). Different uppercase letters above the error bars indicated significant differences among sampling time points under the same treatment (p < 0.05). CK, control treatment without contamination; OTC0.1, OTC addition at 0.1 mg·kg−1; OTC1, OTC addition at 1 mg·kg−1; Cd5, Cd addition at 5 mg·kg−1; OTC0.1 + Cd5, OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1; OTC1 + Cd5, OTC addition at 1 mg·kg−1, and Cd addition at 5 mg·kg−1.

 Individual and combined effects of OTC and Cd contamination on AOA and AOB community structures

The principal coordinate analysis (PCoA) based on Bray-Curtis distance matrices was used to analyze the differences in the community structure of AOA and AOB among different treatments. The results showed that the AOA and AOB community patterns were significantly influenced by individual and combined pollution of OTC and Cd (Figures 5, 6). As for the AOA community, the first two principle components explained 57.65%, 60.44%, 68.54%, and 62.79% of the total variation of the AOA community composition in chestnut, brown, paddy soils, and latosols, respectively (Figure 5). The PCoA plots demonstrated distinct AOA community separations among treatments including OTC1, OTC0.1 + Cd5, OTC1 + Cd5 and other treatments (CK, Cd5, OTC0.1) along the first principle coordinates in chestnut soils and paddy soils (Figures 5A,C). Distinct separations were observed between high OTC treatments (OTC1, OTC1 + Cd5) and other treatments (CK, Cd5, OTC0.1, OTC0.1 + Cd5) along the first principle coordinates in brown soils (Figure 5B). Distinct separations were observed between low OTC treatments (OTC0.1, OTC0.1 + Cd5) and other treatments (CK, Cd5, OTC1, OTC1 + Cd5) along the first principle coordinates in latosols soils (Figure 5D).

Principal coordinate analysis (PCoA) of the Bray–Curtis distance matrix representing differences in the community structure of AOA under CK, OTC0.1, OTC1, Cd5, OTC0.1 + Cd5, and OTC1 + Cd5 treatments in chestnut soils (A), brown soils (B), paddy soils (C), and latosols (D). CK, control treatment without contamination; OTC0.1, OTC addition at 0.1 mg·kg−1; OTC1, OTC addition at 1 mg·kg−1; Cd5, Cd addition at 5 mg·kg−1; OTC0.1 + Cd5, OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1; OTC1 + Cd5, OTC addition at 1 mg·kg−1, and Cd addition at 5 mg· kg−1.

Principal coordinate analysis (PCoA) of the Bray–Curtis distance matrix representing differences in the community structure of AOB under CK, OTC0.1, OTC1, Cd5, OTC0.1 + Cd5, and OTC1 + Cd5 treatments in chestnut soils (A), brown soils (B), paddy soils (C), and latosols (D). CK, control treatment without contamination; OTC0.1, OTC addition at 0.1 mg·kg−1; OTC1, OTC addition at 1 mg·kg−1; Cd5, Cd addition at 5 mg·kg−1; OTC0.1 + Cd5, OTC addition at 0.1 mg·kg−1 and Cd addition at 5 mg·kg−1; OTC1 + Cd5, OTC addition at 1 mg·kg−1, and Cd addition at 5 mg· kg−1.

With regards to the AOB community, the first two principle components explained 60.07, 50.35, 89.71, and 79.57% of the total variation of AOB community patterns in chestnut, brown, paddy, and latosols soils, respectively (Figure 6). The PCoA plots demonstrated distinct AOB community separations between Cd treatments (with or without OTC) and other treatments (CK, OTC0.1, OTC1) along the first principle coordinates in chestnut soils (Figure 6A). Distinct AOB community separations between treatments including OTC1, OTC0.1 + Cd5, OTC1 + Cd5 and other treatments (CK, Cd5, OTC0.1) along the first principle coordinates in brown soils (Figure 6B). No significant separations were observed in AOB communities in paddy soils among various treatments (Figure 6C). In latosols, individual OTC or Cd pollution significantly separated AOB community from combined pollution treatments (OTC0.1 + Cd5, OTC1 + Cd5) along the first principle coordinates (Figure 6D).

 Individual and combined effects of OTC and Cd contamination on nitrification

The effects of soil type and pollution on NNRs, AOA and AOB abundance were analyzed by two-factor ANOVA analysis (Table 3). Either Soil type or contamination treatments exerted significant effects on AOA (p < 0.01, R2 = 0.990) and AOB (p < 0.01, R2 = 0.990). In addition, significant effects of soil types alone (p < 0.01, R2 = 0.911), and pollution alone on NNRs (p < 0.05, R2 = 0.911) were also observed. The interactive effects of soil types and pollution were significant on soil nitrification and associated functional microbes. Moreover, OTC and Cd exerted antagonistic inhibitory effects on AOB in all tested soils. As for AOA, the interaction of OTC and Cd was synergistic in Latosols soils and antagonistic in chestnut and paddy soils (Table 4).

Two-factor ANOVA based on soil types and pollution treatments.

NNRs at 56 days AOA amoA at 56 days AOB amoA at 56 days
Factor df F p df F p df F p
Soil types 3 150.868 <0.01** 3 1475.664 <0.01** 3 1396.397 <0.01**
Contamination treatments 5 2.59 <0.05* 5 30.385 <0.01** 5 124.521 <0.01**
Soil types and contamination treatments 15 1.795 <0.05* 15 11.829 <0.01** 15 7.697 <0.01**

represents significance at p < 0.05.

extremely indicates significance at p < 0.01.

df, degrees of freedom.

Interaction types of OTC and Cd on AOA and AOB at day 56.

Microbe Treatments Inhibition rate (%) Type of interaction
Chestnut soils Brown soils Paddy soils Latosols Chestnut soils Brown soils Paddy soils Latosols
AOA OTC0.1 4.04 4.61 12.28 5.01
OTC1 23.25 19.31 38.89 9.44
Cd5 30.7 19.26 54.06 16.88
OTC0.1 + Cd5 32.34 28.75 51.3 39.51 A S A S
OTC1 + Cd5 44.35 32.57 60.05 53.88 A A A S
AOB OTC0.1 16.29 7.66 42.76 9.28
OTC1 31.51 17.67 59.28 18.95
Cd5 29.48 26.65 62.15 23.57
OTC0.1 + Cd5 35.85 30.16 67.95 20.66 A A A A
OTC1 + Cd5 48.44 33.9 70.63 38.97 A A A A

If the combined inhibition rate of OTC and Cd is greater than the sum of the individual inhibition rates of each contaminant, the type of interaction is considered as synergistic effect (S), or is considered as antagonistic effect (A).

Correlation analysis revealed that NNRs showed strong and positive relationships with AOA (r = 0.546, p < 0.01; r = 0.460, p < 0.05; r = 0.382, p < 0.01) and AOB (r = 0.619, p < 0.01; r = 0.576, p < 0.01; r = 0.492, p < 0.01) under OTC pollution, Cd pollution, as well as under combined OTC and Cd contamination, respectively. SEM analysis was constructed to further investigate the direct or indirect effects of pollutants on NNRs (Figure 7). The results indicated the directly and negatively effects of combined and single OTC and Cd treatments on AOB and AOA amoA genes. The inhibitory effects of pollutants on AOB were stronger than those on AOA. Furthermore, both AOA and AOB showed a significantly positive influence on NNRs. Individual pollution of OTC inhibited NNRs by influencing ammonia oxidizers and Cd contaminant (presence or absence of OTC) had indirect effects on nitrification activity via inhibiting mineral N and ammonia oxidizers.

Structural equation models (SEM) showing the direct and indirect effects of contamination treatments on AOA, AOB, NH4+-N, and NNRs. * indicates significant differences at p < 0.05, ** indicates significant differences at p < 0.01, and *** indicates significant differences at p < 0.001. Numbers on arrows indicate significant standardized path coefficients. Positive effects are represented by red lines, while negative effects are represented by black lines. Continuous and dashed lines represent significant and nonsignificant relationships, respectively.

 Discussion OTC inhibited soil nitrification process by inhibiting ammonia oxidizers

The addition of OTC (1 mg·kg−1) alone significantly decreased NO3−-N levels compared to the control treatment in brown and paddy soils during the whole incubation period, while this effect was only observed at days 14 and 28 in chestnut and latosols soils (Supplementary Figure S2). Likewise, previous researches also revealed that 1 mg·kg−1 OTC inhibited nitrification and decreased NO3−-N concentrations in long-term non-farmed soils (Liu et al., 2020) and silt loam soils (Cao et al., 2016). The addition of 3,300 mg·kg−1 OTC significantly inhibited nitrification after incubation for 94 days in pristine sandy loam soils (Piotrowska-Seget et al., 2008). Trifonova et al. (2021) found that OTC (200 mg·kg−1) inhibited urease activity but had no significant effect on nitrification activity in sod-podzolic soils during 90 days incubation. These discrepancies could be attributed to the fact that the active N nutrient substrate could be depleted in soils, and the antibiotics degraded as the half-life values of 1 and 100 mg·kg−1 OTC were 31.5 and 33 days, respectively (Liu et al., 2020; Tang et al., 2020).

OTC1 significantly inhibited AOB and AOA in all tested soils, while OTC0.1 inhibited AOA and AOB only in paddy soils (Figures 3, 4). In addition, AOB was inhibited more strongly than AOA under OTC1 in brown and latosols soils, and under OTC0.1 in chestnut and latosols soils. Combined with the results that AOA showed remarkably higher abundance than AOB, we speculated that AOA was more tolerant than AOB to OTC pollution. Similarly, previous results reported that OTC addition as well as the individual and combined contaminations of SM2 and Cu inhibited AOB more serious than AOA, and thus AOB was more sensitive to contaminants than AOA (Wang et al., 2019; Tang et al., 2020). The possible reason might be due to the fact that archaea shared peptidoglycan-free cell walls and ether lipids membranes, which were more resistant to antibiotics (Conkle and White, 2012; Tang et al., 2020). Besides, the inhibitory effects of OTC on ammonia oxidizers could be explained by the fact that the nitrobacteria were typically Gram-negative bacteria, whose growth was inhibited by OTC (Chopra and Roberts, 2001). In contrast, some studies found that OTC addition could stimulate ammonia oxidizers rather than inhibit their growth. For example, Cao et al. (2016) revealed that 1 mg·kg−1 OTC significantly increased AOB abundance and inhibited AOA amoA gene in silt loam soils. However, Wang L. et al. (2018) proposed that AOA was facilitated by low concentrations of OTC (23 mg·kg−1) and inhibited by high concentrations (92 and 368 mg·kg−1). This could result from the antibiotic content in the soil exceeded the concentration threshold even cause complete inhibition of microbial populations. In addition, AOA and AOB share different cell structure, hence their tolerance and sensitivity to same environment are different. Besides, our results revealed that AOA (r = 0.546, p < 0.01) and AOB (r = 0.619, p < 0.01) amoA genes were strongly and positively correlated with NNRs under OTC treatments. Combined with the SEM analysis results that OTC pollution indirectly inhibited nitrification activity by inhibiting AOA and AOB abundance (Figure 7), we concluded that both AOA and AOB dominated nitrification in OTC-contaminated soils. In addition, the results revealed that the inhibition effects of OTC on nitrification were increased with the increasing concentrations of OTC, which was further supported by previous researchers who determined dose-dependence of the inhibitory effects on AOA and AOB (Wang J. et al., 2018; Wang L. et al., 2018).

 Cd inhibited soil nitrification process by influencing mineral N and ammonia oxidizers

In this study, Cd addition alone significantly inhibited the NNRs in chestnut, brown and paddy soils at 14 and 28 days (Figure 2). Moreover, Cd alone decreased the AOB amoA gene copy numbers in four soils and decreased AOA abundance only in chestnut and paddy soils (Figures 3, 4). Both correlation analysis and SEM analysis indicated that Cd inhibited nitrification activity through influencing mineral N and AOA/AOB (Figure 7). During the later stages of incubation, nitrification gradually recovered, which could be explained by the screening effect of Cd stress on Cd-resistant species, making the newly selected microorganisms more resistant to the toxic effects of Cd than the native species, resulting in the gradual recovery of soil microbial functions and enzyme activities. For example, Ruyters et al. (2010) found the gradual development of a Zn-tolerant nitrifying community, whose ability to adapt to Zn stress in the soil increased after a 12-month Zn-polluted incubation. The results that Cd contamination decreased the N cycle microbial abundance and slowed the N transformation (Oka and Uchida, 2018; Afzal et al., 2019) might arise from the fact that Cd2+ could replace the active site in ammonia monooxygenase (AMO), thus preventing AMO from oxidizing NH3 to NH2OH (Radniecki and Ely, 2008; Kapoor et al., 2015). It was also reported that Cd stress (8 mg·kg−1) produced sustained stimulatory effect on both AOB and AOA, with AOA being most strongly stimulated in the later incubation stage (Zhao et al., 2021). Some studies demonstrated that low levels of Cd (2 and 5 mg·kg−1) addition would stimulate nitrification, while inhibition effects were only found at high levels of Cd addition (10 and 20 mg·kg−1) in alluvial soils (Chen et al., 2001). Furthermore, Salam et al. (2020) reported that Cd contamination significantly affected the structure and function of soil microbial communities, resulting in massive loss of non-adapted native microorganisms, as well as changed the soil physicochemical properties. Thus, the response of nitrification process to Cd addition varied with Cd concentrations and soil types (Liu et al., 2019).

 Combined pollution of OTC and Cd inhibited soil nitrification process by influencing mineral N and ammonia oxidizers

Combined contamination of OTC and Cd obviously inhibited AOA and AOB amoA gene copies as well as NNRs in all soils. And the response time and extent of nitrification activity to the combined contamination of OTC and Cd varied with soil types. For example, combined addition of OTC (0.1 and 1 mg·kg−1) and Cd decreased the NNRs in brown soils at day 56, while the NNRs content wasn’t significantly changed in chestnut and paddy soils (Figure 2). This could be attributed to the fact that interaction effects of antibiotics and heavy metals were affected by soil properties, and the complexation of antibiotics and heavy metals varied with soils, thus altering the sorption and desorption of antibiotics and heavy metals (Xu et al., 2015; Zhang et al., 2018). It had been reported that coexistence of TC and Cd (II) reduced each other’s mobility in alkaline environment such as calcareous and saline soils, while Cd (II) had no effect on TC mobility in acidic soils (Du et al., 2021). In addition, the absorption of antibiotics by SOM can be significantly prolonged and enhanced (Kong et al., 2012; Álvarez-Esmorís et al., 2020). And SOM was the dominant adsorbent for Cd and decreased Cd availability or mobilization (Filipović et al., 2018; Chen et al., 2020). In addition, OTC0.1 significantly inhibited AOA in paddy soils while 0.1 mg·kg−1 of OTC had no significant toxic effects on AOA in other tested soils. Correlation analysis and SEM analysis further confirmed that combined pollution of OTC and Cd inhibited soil nitrification process by influencing mineral N and ammonia oxidizers (Figure 7). This was explained by antibiotics and heavy metals inhibited the replication and transcription of the function gene amoA of AOB and AOA, thus inhibiting the oxidation of NH4+-N and reducing soil NO3−-N concentrations (Kapoor et al., 2015; Tong et al., 2015; Tang et al., 2020). Hence, we speculated that the variation of pH and organic matter content in different soils could be the main reason for the different toxic effect. And an antagonism between OTC and Cd on AOB was observed in all tested soils. However, OTC and Cd had synergistic inhibitions on AOA in brown and latosols soils (Table 4). Wang J. et al. (2018) found enrofloxacin (ENR) and Cd interaction effects mainly antagonistic on AOB, which might be related to various factors such as concentrations, incubation time of mixture contaminants as well as soil properties. In addition, the combined pollution of OTC and Cd had greater negative effects on nitrification than those with the single contaminant. This resulted from the relative inhibition of dissipation and increased sorption of OTC by Cd, and OTC-Cd complexes might be formed by OTC and Cd, and the complexes showed greater toxic effects than the containments used alone. Jia et al. (2008) found that the TC-Cu2+ complexes formed in the presence of TC had greater affinity for the soil surface than Cu2+ alone, which enhanced the adsorption of Cu2+ in red soils and paddy soils. Thus, the antibiotics bound metals via multiple coordination sited to produce the complex (Zhang et al., 2012; Tong et al., 2015). This result was further supported by Tang et al. (2020), who revealed that combined pollution of antibiotics and heavy metals had greater toxicity than single pollution in soil nitrification process. Similarly, Kong et al. (2006) found that OTC and Cu also had significant negative effects on soil microbial community function, with the negative effects were more significant in the co-existence of both contaminants.

Additionally, our results also suggested that AOA amoA gene abundance was higher than AOB amoA gene abundance and the inhibitory effects of OTC and Cd on AOB were greater than on AOA in all tested soils. This was consistent with previous studies showing that the AOA amoA gene copies were higher than AOB amoA gene copies in many soils, and AOB was more sensitive than AOA to the antibiotics and heavy metals combined treatments (Leininger et al., 2006; Wang L. et al., 2018; Wang et al., 2019; Tang et al., 2020). Moreover, the combined toxicity of two pollutants was significantly greater than the toxicity of either pollutant alone and was most pronounced in the early stage during incubation, but the inhibitory effects decreased with incubation time. The results were partially consistent with previous studies (Wang J. et al., 2018; Wang L. et al., 2018), which might be due to the degradation of OTC and the reduction of Cd potency, and the concentration of antibiotics and bioavailable heavy metals at day 56 might be below the limit for inhibiting the activity of ammonia-oxidizing microorganisms. The inhibition effect of OTC and Cd on AOB in latosols soils enhanced with increasing OTC levels, as well as on NNRs in brown soils, indicated that the dose-dependent relationship also existed in the combined OTC and Cd pollution (Ma et al., 2016; Wang L. et al., 2018). In addition, due to the large differences in the properties of different antibiotics and heavy metals, only typical types of antibiotics and heavy metals were selected for this study. Therefore, the results of this study could provide a reference for the mechanism of the action of single or combined contamination of heavy metals and antibiotics on soil N cycle. Moreover, an experiment investigated the structure and form of complex contaminants after OTC and Cd exposure to soils would be considered in the near future to comprehensively understand the morphology and structure of contaminants in different soil types after OTC and Cd pollution.

 Conclusion

The effects of individual and combined OTC and Cd on soil nitrification were examined in the study. The single and combined application of OTC and Cd showed different levels of significant inhibitory effect on soil nitrification. The combined toxicity of OTC and Cd was generally greater than that from either contaminant acting alone. Moreover, the toxic effects of single and combined contaminants on nitrification were significantly enhanced with increasing concentrations of OTC, and varied with soil types because of soil pH and SOM. The interaction of OTC and Cd on AOB and AOA was mainly antagonistic. Both AOA and AOB dominated in nitrification in the individual and combined contaminated soils, and AOB was more sensitive than AOA to both OTC and Cd toxicity. OTC inhibited nitrification activity by inhibiting ammonia oxidizers, and Cd (with or without OTC) inhibited nitrification through influencing mineral N and ammonia oxidizers.

 Data availability statement

The authors confirm that the data supporting the findings of this study are available within the article and its Supplementary materials.

 Author contributions

HP and YZ provided the idea of this study. XC, HP, and YZ designed the study. XC completed most of the experimental procedures, data analysis and writing the manuscript. JH performed the correlation analysis and SEM analysis in the revised manuscript. WZ, HZ, JL, and JH analyzed the soil physicochemical properties. HP and YZ supervised the work. HP, YZ, YL, HW, and QY revised the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This research was funded by the Major Science and Technology Innovation Projects in Shandong Province (2021CXGC010804).

 Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1062703/full#supplementary-material

 References Afzal M. Yu M. Tang C. Zhang L. Muhammad N. Zhao H. . (2019). The negative impact of cadmium on nitrogen transformation processes in a paddy soil is greater under non-flooding than flooding conditions. Environ. Int. 129, 451–460. doi: 10.1016/j.envint.2019.05.058, PMID: 31154147 Aigle A. Prosser J. I. Gubry-Rangin C. (2019). The application of high-throughput sequencing technology to analysis of amoA phylogeny and environmental niche specialisation of terrestrial bacterial ammonia-oxidisers. Environ. Microbiome 14:3. doi: 10.1186/s40793-019-0342-6, PMID: 33902715 Álvarez-Esmorís C. Conde-Cid M. Fernández-Calviño D. Fernández-Sanjurjo M. J. Núñez-Delgado A. Álvarez-Rodríguez E. . (2020). Adsorption-desorption of doxycycline in agricultural soils: batch and stirred-flow-chamber experiments. Environ. Res. 186:109565. doi: 10.1016/j.envres.2020.109565, PMID: 32371275 Bai J. Yu P. Wen X. Wang W. Jia J. Wang X. (2021). Effects of cadmium addition on net nitrogen mineralization processes in the urban constructed wetland soils of a Chinese delta. Environ. Geochem. Health 43, 1155–1164. doi: 10.1007/s10653-020-00597-0, PMID: 32419088 Baker-Austin C. Wright M. S. Stepanauskas R. McArthur J. V. (2006). Co-selection of antibiotic and metal resistance. Trends Microbiol. 14, 176–182. doi: 10.1016/j.tim.2006.02.006 Beeckman F. Motte H. Beeckman T. (2018). Nitrification in agricultural soils: impact, actors and mitigation. Curr. Opin. Biotechnol. 50, 166–173. doi: 10.1016/j.copbio.2018.01.014, PMID: 29414056 Cai L.-M. Wang Q.-S. Luo J. Chen L.-G. Zhu R.-L. Wang S. . (2019). Heavy metal contamination and health risk assessment for children near a large cu-smelter in Central China. Sci. Total Environ. 650, 725–733. doi: 10.1016/j.scitotenv.2018.09.081, PMID: 30212703 Cao J. Wang C. Dou Z. Ji D. (2016). Independent and combined effects of oxytetracycline and antibiotic-resistant Escherichia coli O157: H7 on soil microbial activity and partial nitrification processes. Soil Biol. Biochem. 98, 138–147. doi: 10.1016/j.soilbio.2016.03.014 Chen C. Liang H. Gao D. (2019). Community diversity and distribution of ammonia-oxidizing archaea in marsh wetlands in the black soil zone in north-East China. Front. Environ. Sci. Eng. 13:58. doi: 10.1007/s11783-019-1146-z Chen W. Peng L. Hu K. Zhang Z. Peng C. Teng C. . (2020). Spectroscopic response of soil organic matter in mining area to Pb/cd heavy metal interaction: a mirror of coherent structural variation. J. Hazard. Mater. 393:122425. doi: 10.1016/j.jhazmat.2020.122425, PMID: 32135370 Chen Y. Wang K. Lin Q. Yang Y. (2001). Effects of heavy metals on ammonification, nitrification and denitrification in maize rhizosphere. Pedosphere 11, 115–122. Chopra I. Roberts M. (2001). Tetracycline antibiotics: mode of action, applications, molecular biology, and epidemiology of bacterial resistance. Microbiol. Mol. Biol. Rev. 65, 232–260. doi: 10.1128/MMBR.65.2.232-260.2001, PMID: 11381101 Clemens S. Ma J. F. (2016). Toxic heavy metal and metalloid accumulation in crop plants and foods. Annu. Rev. Plant Biol. 67, 489–512. doi: 10.1146/annurev-arplant-043015-112301, PMID: 27128467 Conde-Cid M. Fernández-Calviño D. Núñez-Delgado A. Fernández-Sanjurjo M. J. Arias-Estévez M. Álvarez-Rodríguez E. (2020). Estimation of adsorption/desorption Freundlich's affinity coefficients for oxytetracycline and chlortetracycline from soil properties: experimental data and pedotransfer functions. Ecotoxicol. Environ. Saf. 196:110584. doi: 10.1016/j.ecoenv.2020.110584, PMID: 32278142 Conkle J. L. White J. R. (2012). An initial screening of antibiotic effects on microbial respiration in wetland soils. J. Environ. Sci. Health A 47, 1381–1390. doi: 10.1080/10934529.2012.672315, PMID: 22571526 Daims H. Lebedeva E. V. Pjevac P. Han P. Herbold C. Albertsen M. . (2015). Complete nitrification by Nitrospira bacteria. Nature 528, 504–509. doi: 10.1038/nature16461, PMID: 26610024 Du H. Du J. Liu F. Zhang Y. Guo H. Wan D. (2021). Binding of tetracycline on soil phyllosilicates with cd (II) as affected by pH and mineral type. J. Soil. Sediment. 21, 775–783. doi: 10.1007/s11368-020-02867-x Filipović L. Romić M. Romić D. Filipović V. Ondrašek G. (2018). Organic matter and salinity modify cadmium soil (phyto)availability. Ecotoxicol. Environ. Saf. 147, 824–831. doi: 10.1016/j.ecoenv.2017.09.041, PMID: 28968923 Francis C. A. Roberts K. J. Beman J. M. Santoro A. E. Oakley B. B. (2005). Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. Proc. Natl. Acad. Sci. 102, 14683–14688. doi: 10.1073/pnas.0506625102, PMID: 16186488 Galloway J. N. Townsend A. R. Erisman J. W. Bekunda M. Cai Z. Freney J. R. . (2008). Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science 320, 889–892. doi: 10.1126/science.1136674, PMID: 18487183 Grenni P. Ancona V. Barra C. A. (2018). Ecological effects of antibiotics on natural ecosystems: a review. Microchem. J. 136, 25–39. doi: 10.1016/j.microc.2017.02.006 Guo T. Lou C. Zhai W. Tang X. Hashmi M. Z. Murtaza R. . (2018). Increased occurrence of heavy metals, antibiotics and resistance genes in surface soil after long-term application of manure. Sci. Total Environ. 635, 995–1003. doi: 10.1016/j.scitotenv.2018.04.194, PMID: 29710621 Guo H. Xue S. Nasir M. Gu J. Lv J. (2021). Impacts of cadmium addition on the alteration of microbial community and transport of antibiotic resistance genes in oxytetracycline contaminated soil. J. Environ. Sci. 99, 51–58. doi: 10.1016/j.jes.2020.04.015, PMID: 33183716 Haider F. U. Liqun C. Coulter J. A. Cheema S. A. Wu J. Zhang R. . (2021). Cadmium toxicity in plants: impacts and remediation strategies. Ecotoxicol. Environ. Saf. 211:111887. doi: 10.1016/j.ecoenv.2020.111887, PMID: 33450535 He S. He Z. Yang X. Stoffella P. J. Baligar V. C. (2015). “Chapter four - soil biogeochemistry, plant physiology, and phytoremediation of cadmium-contaminated soils” in Advances in agronomy. 134. ed. Sparks D. L. (Academic Press), 135–225. He Y. Liu C. Tang X.-Y. Xian Q.-S. Zhang J.-Q. Guan Z. (2019). Biochar impacts on sorption-desorption of oxytetracycline and florfenicol in an alkaline farmland soil as affected by field ageing. Sci. Total Environ. 671, 928–936. doi: 10.1016/j.scitotenv.2019.03.414 He H. Liu H. Shen T. Wei S. Dai J. Wang R. (2018). Influence of Cu application on ammonia oxidizers in fluvo-aquic soil. 321, 141–150. doi: 10.1016/j.geoderma.2018.01.037 Hu B. Shao S. Ni H. Fu Z. Hu L. Zhou Y. . (2020). Current status, spatial features, health risks, and potential driving factors of soil heavy metal pollution in China at province level. Environ. Pollut. 266:114961. doi: 10.1016/j.envpol.2020.114961, PMID: 32622003 Huang Y. Wang L. Wang W. Li T. He Z. Yang X. (2019). Current status of agricultural soil pollution by heavy metals in China: a meta-analysis. Sci. Total Environ. 651, 3034–3042. doi: 10.1016/j.scitotenv.2018.10.185, PMID: 30463153 Jia D.-A. Zhou D.-M. Wang Y.-J. Zhu H.-W. Chen J.-L. (2008). Adsorption and cosorption of cu(II) and tetracycline on two soils with different characteristics. Geoderma 146, 224–230. doi: 10.1016/j.geoderma.2008.05.023 Kapoor V. Li X. Elk M. Chandran K. Impellitteri C. A. Santo Domingo J. W. (2015). Impact of heavy metals on transcriptional and physiological activity of nitrifying bacteria. Environ. Sci. Technol. 49, 13454–13462. doi: 10.1021/acs.est.5b02748, PMID: 26501957 Kavvadias V. Paschalidis C. Vavoulidou E. Petropoulos D. Koriki A. (2012). Effects of soil amended with cadmium and Lead on growth, yield, and metal accumulation and distribution in parsley. Commun. Soil Sci. Plant Anal. 43, 161–175. doi: 10.1080/00103624.2012.634708 Khanom A. Azad M. A. K. Ali M. M. Ali M. Y. Biswas S. K. Rahman M. M. (2021). Plants and microbes' responses to the net nitrification rates of chemical fertilizers in vegetable soils. Appl. Soil Ecol. 158:103783. doi: 10.1016/j.apsoil.2020.103783 Knapp C. W. Dolfing J. Ehlert P. A. I. Graham D. W. (2010). Evidence of increasing antibiotic resistance gene abundances in archived soils since 1940. Environ. Sci. Technol. 44, 580–587. doi: 10.1021/es901221x, PMID: 20025282 Kong W. Li C. Dolhi J. M. Li S. He J. Qiao M. (2012). Characteristics of oxytetracycline sorption and potential bioavailability in soils with various physical–chemical properties. Chemosphere 87, 542–548. doi: 10.1016/j.chemosphere.2011.12.062, PMID: 22245075 Kong W.-D. Zhu Y.-G. Fu B.-J. Marschner P. He J.-Z. (2006). The veterinary antibiotic oxytetracycline and cu influence functional diversity of the soil microbial community. Environ. Pollut. 143, 129–137. doi: 10.1016/j.envpol.2005.11.003, PMID: 16413090 Kuypers M. M. M. Marchant H. K. Kartal B. (2018). The microbial nitrogen-cycling network. Nat. Rev. Microbiol. 16, 263–276. doi: 10.1038/nrmicro.2018.9 Leininger S. Urich T. Schloter M. Schwark L. Qi J. Nicol G. W. . (2006). Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature 442, 806–809. doi: 10.1038/nature04983, PMID: 16915287 Leong Y. K. Chang J.-S. (2020). Bioremediation of heavy metals using microalgae: recent advances and mechanisms. Bioresour. Technol. 303:122886. doi: 10.1016/j.biortech.2020.122886, PMID: 32046940 Li Z.-J. Qi W.-N. Feng Y. Liu Y.-W. Ebrahim S. Long J. (2019). Degradation mechanisms of oxytetracycline in the environment. J. Integr. Agric. 18, 1953–1960. doi: 10.1016/S2095-3119(18)62121-5 Li W. Shi Y. Gao L. Liu J. Cai Y. (2012). Occurrence of antibiotics in water, sediments, aquatic plants, and animals from Baiyangdian Lake in North China. Chemosphere 89, 1307–1315. doi: 10.1016/j.chemosphere.2012.05.079, PMID: 22698376 Li Y.-X. Xiong X. Lin C.-Y. Zhang F.-S. Wei L. Wei H. (2010). Cadmium in animal production and its potential hazard on Beijing and Fuxin farmlands. J. Hazard. Mater. 177, 475–480. doi: 10.1016/j.jhazmat.2009.12.057, PMID: 20060219 Lin Y.-X. Ye G. Ding W.-X. Hu H.-W. Zheng Y. Fan J.-B. . (2020). Niche differentiation of comammox Nitrospira and canonical ammonia oxidizers in soil aggregate fractions following 27-year fertilizations. Agric. Ecosyst. Environ. 304:107147. doi: 10.1016/j.agee.2020.107147 Liu M. Cao J. Wang C. (2020). Bioremediation by earthworms on soil microbial diversity and partial nitrification processes in oxytetracycline-contaminated soil. Ecotoxicol. Environ. Saf. 189:109996. doi: 10.1016/j.ecoenv.2019.109996, PMID: 31785943 Liu Y. Xue C. Yu S. Li F. (2019). Variations of abundance and community structure of ammonia oxidizers and nitrification activity in two Paddy soils polluted by heavy metals. Geomicrobiol. J. 36, 1–10. doi: 10.1080/01490451.2018.1471108 Lu L. Chen C. Ke T. Wang M. Sima M. Huang S. (2022). Long-term metal pollution shifts microbial functional profiles of nitrification and denitrification in agricultural soils. Sci. Total Environ. 830:154732. doi: 10.1016/j.scitotenv.2022.154732, PMID: 35346706 Ma T. Pan X. Chen L. K. Liu W. Christie P. Luo Y. . (2016). Effects of different concentrations and application frequencies of oxytetracycline on soil enzyme activities and microbial community diversity. Eur. J. Soil Biol. 76, 53–60. doi: 10.1016/j.ejsobi.2016.07.004 Mifflin A. L. Konek C. T. Geiger F. M. (2006). Tracking oxytetracyline mobility across environmental interfaces by second harmonic generation. J. Phys. Chem. B 110, 22577–22585. doi: 10.1021/jp063089p, PMID: 17092004 Oka M. Uchida Y. (2018). Heavy metals in slag affect inorganic N dynamics and soil bacterial community structure and function. Environ. Pollut. 243, 713–722. doi: 10.1016/j.envpol.2018.09.024, PMID: 30228069 Okereafor U. Makhatha M. Mekuto L. Uche-Okereafor N. Sebola T. Mavumengwana V. (2020). Toxic metal implications on agricultural soils, plants, animals, aquatic life and human health. Int. J. Environ. Res. Public Health 17:2204. doi: 10.3390/ijerph17072204, PMID: 32218329 Pan M. Chu L. M. (2016). Adsorption and degradation of five selected antibiotics in agricultural soil. Sci. Total Environ. 545-546, 48–56. doi: 10.1016/j.scitotenv.2015.12.040, PMID: 26745292 Pan H. Feng H. Liu Y. Lai C.-Y. Zhuge Y. Zhang Q. . (2021). Grazing weakens competitive interactions between active methanotrophs and nitrifiers modulating greenhouse-gas emissions in grassland soils. ISME Commun. 1:74. doi: 10.1038/s43705-021-00068-2 Pan H. Xie K. Zhang Q. Jia Z. Xu J. Di H. . (2018). Archaea and bacteria respectively dominate nitrification in lightly and heavily grazed soil in a grassland system. Biol. Fertil. Soils 54, 41–54. doi: 10.1007/s00374-017-1236-7 Piotrowska-Seget Z. Engel R. Nowak E. Kozdrój J. (2008). Successive soil treatment with captan or oxytetracycline affects non-target microorganisms. World J. Microbiol. Biotechnol. 24, 2843–2848. doi: 10.1007/s11274-008-9815-2 Qin G. Niu Z. Yu J. Li Z. Ma J. Xiang P. (2021). Soil heavy metal pollution and food safety in China: effects, sources and removing technology. Chemosphere 267:129205. doi: 10.1016/j.chemosphere.2020.129205, PMID: 33338709 Qin J. Xiong H. Ma H. Li Z. (2019). Effects of different fertilizers on residues of oxytetracycline and microbial activity in soil. Environ. Sci. Pollut. Res. 26, 161–170. doi: 10.1007/s11356-018-3603-9, PMID: 30387057 Radniecki T. S. Ely R. L. (2008). Zinc chloride inhibition of Nitrosococcus mobilis. Biotechnol. Bioeng. 99, 1085–1095. doi: 10.1002/bit.21672, PMID: 17929325 Rahman Z. Singh V. P. (2019). The relative impact of toxic heavy metals (THMs) (arsenic (as), cadmium (cd), chromium (Cr)(VI), mercury (hg), and lead (Pb)) on the total environment: an overview. Environ. Monit. Assess. 191:419. doi: 10.1007/s10661-019-7528-7, PMID: 31177337 Rotthauwe J.-H. Witzel K.-P. Liesack W. (1997). The ammonia monooxygenase structural gene amoA as a functional marker: molecular fine-scale analysis of natural ammonia-oxidizing populations. Appl. Environ. Microbiol. 63, 4704–4712. doi: 10.1128/aem.63.12.4704-4712.1997, PMID: 9406389 Ruyters S. Mertens J. Springael D. Smolders E. (2010). Stimulated activity of the soil nitrifying community accelerates community adaptation to Zn stress. Soil Biol. Biochem. 42, 766–772. doi: 10.1016/j.soilbio.2010.01.012 Salam L. B. Obayori O. S. Ilori M. O. Amund O. O. (2020). Effects of cadmium perturbation on the microbial community structure and heavy metal resistome of a tropical agricultural soil. Bioresourc. Bioprocess. 7, 1–19. doi: 10.1186/s40643-020-00314-w Sarmah A. K. Meyer M. T. Boxall A. B. A. (2006). A global perspective on the use, sales, exposure pathways, occurrence, fate and effects of veterinary antibiotics (VAs) in the environment. Chemosphere 65, 725–759. doi: 10.1016/j.chemosphere.2006.03.026, PMID: 16677683 Shi T. Zhang Y. Gong Y. Ma J. Wei H. Wu X. . (2019). Status of cadmium accumulation in agricultural soils across China (1975–2016): from temporal and spatial variations to risk assessment. Chemosphere 230, 136–143. doi: 10.1016/j.chemosphere.2019.04.208, PMID: 31103859 Tang Q. Xia L. Ti C. Zhou W. Fountain L. Shan J. . (2020). Oxytetracycline, copper, and zinc effects on nitrification processes and microbial activity in two soil types. Food Energy Secur. 9:e248. doi: 10.1002/fes3.248 R Core Team. (2013). R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, Vienna, Austria). Tolls J. (2001). Sorption of veterinary Pharmaceuticals in Soils: a review. Environ. Sci. Technol. 35, 3397–3406. doi: 10.1021/es0003021, PMID: 11563639 Tong F. Zhao Y. Gu X. Gu C. Lee C. C. (2015). Joint toxicity of tetracycline with copper (II) and cadmium (II) to Vibrio fischeri: effect of complexation reaction. Ecotoxicology 24, 346–355. doi: 10.1007/s10646-014-1383-7, PMID: 25398505 Trifonova T. Kosmacheva A. Sprygin A. Chesnokova S. Byadovskaya O. (2021). Enzymatic activity and microbial diversity of sod-Podzolic soil microbiota using 16S rRNA amplicon sequencing following antibiotic exposure. Antibiotics 10:970. doi: 10.3390/antibiotics10080970, PMID: 34439020 Van Boeckel T. P. Brower C. Gilbert M. Grenfell B. T. Levin S. A. Robinson T. P. . (2015). Global trends in antimicrobial use in food animals. Proc. Natl. Acad. Sci. 112, 5649–5654. doi: 10.1073/pnas.1503141112, PMID: 25792457 Wang P. Chen H. Kopittke P. M. Zhao F.-J. (2019a). Cadmium contamination in agricultural soils of China and the impact on food safety. Environ. Pollut. 249, 1038–1048. doi: 10.1016/j.envpol.2019.03.063, PMID: 31146310 Wang J. Wang L. Zhu L. Wang J. (2018). Individual and combined effects of enrofloxacin and cadmium on soil microbial biomass and the ammonia-oxidizing functional gene. Sci. Total Environ. 624, 900–907. doi: 10.1016/j.scitotenv.2017.12.096, PMID: 29274613 Wang L. Wang J. Zhu L. Wang J. (2018). Toxic effects of oxytetracycline and copper, separately or combined, on soil microbial biomasses. Environ. Geochem. Health 40, 763–776. doi: 10.1007/s10653-017-0022-7 Wang P. Wu D. You X. Li W. Xie B. (2019b). Distribution of antibiotics, metals and antibiotic resistance genes during landfilling process in major municipal solid waste landfills. Environ. Pollut. 255:113222. doi: 10.1016/j.envpol.2019.113222, PMID: 31563781 Wang L. Xia X. Zhang W. Wang J. Zhu L. Wang J. . (2019). Separate and joint eco-toxicological effects of sulfadimidine and copper on soil microbial biomasses and ammoxidation microorganisms abundances. Chemosphere 228, 556–564. doi: 10.1016/j.chemosphere.2019.04.165, PMID: 31055070 Xu Y. Yu W. Ma Q. Zhou H. (2015). Interactive effects of sulfadiazine and cu(II) on their sorption and desorption on two soils with different characteristics. Chemosphere 138, 701–707. doi: 10.1016/j.chemosphere.2015.07.049, PMID: 26247413 Yang Q. Li Z. Lu X. Duan Q. Huang L. Bi J. (2018). A review of soil heavy metal pollution from industrial and agricultural regions in China: pollution and risk assessment. Sci. Total Environ. 642, 690–700. doi: 10.1016/j.scitotenv.2018.06.068, PMID: 29909337 Zhang Y. Cai X. Lang X. Qiao X. Li X. Chen J. (2012). Insights into aquatic toxicities of the antibiotics oxytetracycline and ciprofloxacin in the presence of metal: complexation versus mixture. Environ. Pollut. 166, 48–56. doi: 10.1016/j.envpol.2012.03.009, PMID: 22475550 Zhang H. Luo Y. Wu L. Huang Y. Christie P. (2015). Residues and potential ecological risks of veterinary antibiotics in manures and composts associated with protected vegetable farming. Environ. Sci. Pollut. Res. 22, 5908–5918. doi: 10.1007/s11356-014-3731-9, PMID: 25354434 Zhang J. Sui Q. Tong J. Zhong H. Wang Y. Chen M. . (2018). Soil types influence the fate of antibiotic-resistant bacteria and antibiotic resistance genes following the land application of sludge composts. Environ. Int. 118, 34–43. doi: 10.1016/j.envint.2018.05.029, PMID: 29793114 Zhang Q.-Q. Ying G.-G. Pan C.-G. Liu Y.-S. Zhao J.-L. (2015). Comprehensive evaluation of antibiotics emission and fate in the river basins of China: source analysis, multimedia modeling, and linkage to bacterial resistance. Environ. Sci. Technol. 49, 6772–6782. doi: 10.1021/acs.est.5b00729, PMID: 25961663 Zhao H. Lin J. Wang X. Shi J. Dahlgren R. A. Xu J. (2021). Dynamics of soil microbial N-cycling strategies in response to cadmium stress. Environ. Sci. Technol. 55, 14305–14315. doi: 10.1021/acs.est.1c04409, PMID: 34617741 Zhao Y. Tan Y. Guo Y. Gu X. Wang X. Zhang Y. (2013). Interactions of tetracycline with cd (II), cu (II) and Pb (II) and their cosorption behavior in soils. Environ. Pollut. 180, 206–213. doi: 10.1016/j.envpol.2013.05.043, PMID: 23786789 Zhao H. Yu L. Yu M. Afzal M. Dai Z. Brookes P. . (2020). Nitrogen combined with biochar changed the feedback mechanism between soil nitrification and cd availability in an acidic soil. J. Hazard. Mater. 390:121631. doi: 10.1016/j.jhazmat.2019.121631, PMID: 31776087 Zhu Y.-G. Johnson T. A. Su J.-Q. Qiao M. Guo G.-X. Stedtfeld R. D. . (2013). Diverse and abundant antibiotic resistance genes in Chinese swine farms. Proc. Natl. Acad. Sci. 110, 3435–3440. doi: 10.1073/pnas.1222743110, PMID: 23401528 Zhu W. Zhu D. He J. Lian X. Chang Z. Guo R. . (2022). Phytoremediation of soil co-contaminated with heavy metals (HMs) and tetracyclines: effect of the co-contamination and HM bioavailability analysis. J. Soil. Sediment. 22, 2036–2047. doi: 10.1007/s11368-022-03206-y