River restoration changes distributions of antibiotics, antibiotic resistance genes, and microbial community

Potential of nature-based solutions to reduce antibiotics, antimicrobial resistance, and pathogens in aquatic ecosystems. a critical review

This comprehensive scientific review evaluates the effectiveness of nature-based solutions (NBS) in reducing antibiotics (ABs), combating antimicrobial resistance (AMR), and controlling pathogens in various aquatic environments at different river catchment levels. It covers conventional and innovative treatment wetland configurations for wastewater treatment to reduce pollutant discharge into the aquatic ecosystems as well as exploring how river restoration and saltmarshes can enhance pollutant removal. Through the analysis of experimental studies and case examples, the review shows NBS's potential for providing sustainable and cost-effective solutions to improve the health of aquatic ecosystems. It also evaluates the use of diagnostic indicators to predict NBS effectiveness in removing specific pollutants such as ABs and AMR. The review concludes that NBS are feasible for addressing the new challenges stemming from human activities such as the presence of ABs, AMR and pathogens, contributing to a better understanding of NBS, highlighting success stories, addressing knowledge gaps, and providing recommendations for future research and implementation.

Detection of Pharmaceutical Contamination in Amphipods of Lake Baikal by the HPLC-MS Method

Pollution by active ingredients is one of the most significant and widespread forms of pollution on Earth. Medicines can have a negative impact on ecosystems, and contamination can have unpredictable consequences. An urgent and unexplored task is to study the Lake Baikal ecosystem and its organisms for the presence of trace concentrations of active pharmaceutical ingredients. Our study aimed to conduct a qualitative analysis of active pharmaceutical ingredients, and quantitative analysis of ibuprofen in endemic amphipods of Lake Baikal, using methods of high-performance liquid chromatography and mass spectrometry (HPLC-MS). Acetylsalicylic acid (aspirin), ibuprofen, acetaminophen, azithromycin, dimetridazole, metronidazole, amikacin, spiramycin, and some tetracycline antibiotics were detected in the studied littoral amphipods. We also detected different annual loads of active pharmaceutical ingredients on amphipods. Using the multiple reaction monitoring (MRM) mode mentioned in GOST International Technical Standards, we detected molecules, fragmented as amikacin, chlortetracycline, doxycycline, oxytetracycline, dimetridazole, metronidazole and spiramycin. Thus, we first revealed that invertebrates of Lake Baikal can uptake pharmaceutical contaminants in the environment.

Response characteristics and functional predictions of soil microorganisms to heavy metals, antibiotics, and their resistance genes originating from different animal farms amended with Herbaspirillum huttiense

Current understanding is limited regarding technologies that use biochar and microorganisms to simultaneously treat soils contaminated with both veterinary antibiotics (VAs) and heavy metals (HMs) from different animal farms. The contributions of the keystone taxa and their similarities from different animal farms under VA and HM stresses before and after soil remediation should be further investigated as well. An innovative treatment of Herbaspirillum huttiense (HHS1) inoculated waste fungus chaff-based (WFCB) biochar was designed for immobilization of copper (Cu) and zinc (Zn), and the removal of oxytetracycline (OTC), enrofloxacin (ENR), and a subsequent reduction in their resistance genes in soils from pig, cow, and chicken farms. Roles of indigenous microorganisms which can treat soils contaminated with VAs and HMs were summarized. Results showed that available Cu and Zn were reduced by 19.5% and 28.1%, respectively, while 49.8% of OTC and 85.1% of ENR were removed by WFCB-HHS1. The decrease in ENR improved overall microbial community diversity, and the increases in genera HHS1, Pedobacter, Flavobacterium and Aequorivita, along with the decreases of genera Bacillus, Methylobacter, and Fermentimonas were indirectly favorable to treat HMs and VAs in soils from different animal farms. Bacterial communities in different animal farm soils were predominantly influenced by stochastic processes. The regulations of functional genes associated with metabolism and environmental information processing, which contribute to HM and VA defense, were altered when using WFCB-HHS1. Furthermore, the spread of their antibiotic resistance genes was restricted.

Antibiotic resistance genes and mobile genetic elements in different rivers: The link with antibiotics, microbial communities, and human activities

Rivers as a critical sink for antibiotic resistance genes (ARGs), and the distribution and spread of ARGs are related to environmental factors, human activities, and biotic factors (e.g. mobile genetic elements (MGEs)). However, the potential link among ARGs, microbial community, and MGEs in rivers under different antibiotic concentration and human activities remains unclear. In this study, 2 urban rivers (URs), 1 rural-urban river (RUR), and 2 rural rivers (RRs) were investigated to identify the spatial-temporal variation and driving force of ARGs. The total concentration of quinolones (QNs) was 160.1-2151 ng·g-1 in URs, 23.34-1188 ng·g-1 in RUR, and 16.39-85.98 ng·g-1 in RRs. Total population (TP), gross domestic production (GDP), sewage, industrial enterprise (IE), and IEGDP appeared significantly spatial difference in URs, RUR, and RRs. In terms of ARGs, 145-161 subtypes were detected in URs, 59-61 subtypes in RURs, and 46-79 subtypes in RRs. For MGEs, 55-60 MGEs subtypes were detected in URs, 29-30 subtypes in RUR, and 29-35 subtypes in RRs. Significantly positive correlation between MGEs and ARGs were found in these rivers. More ARGs subtypes were related to MGEs in URs than those in RUR and RRs. Overall, MGEs and QNs showed significantly direct positive impact on the abundance of ARGs in all rivers, while microbial community was significantly positive impact on the ARGs abundance in URs and RUR. The ARGs abundance in URs/RUR were directly positive influenced by microbial community/MGEs/socioeconomic elements (SEs)/QNs, while those in RRs were directly positive influenced by QNs/MGEs and indirectly positive impacted by SEs. Most QNs resistance risk showed significantly positive correlation with the abundance of ARGs types. Therefore, not only need to consider the concentration of antibiotics, but also should pay more attention to SEs and MGEs in antibiotics risk management and control.

The response of microbial compositions and functions to chronic single and multiple antibiotic exposure by batch experiment

Understanding the response of the microbial community to external disturbances such as micropollutants is vital for ecological risk evaluation. In this study, the effect of chronic antibiotic exposure on community compositions and functions was investigated by two batch experiments. The first experiment investigated the effect of chronic sulfamethoxazole (SMX) exposure, while the second investigated the combined effect of dissolved organic matter (DOM) sources and multi-antibiotic exposure. The results showed that the community responses to chronic antibiotic exposure depended on the dynamic balance among community resistance, adaptation, recovery, and selection, leading to nonlinear composition diversity variations. The disturbance strength of chronic SMX exposure increased with concentration (0.5-50 μg/L). However, complex sources and structures of coexisting organic matter might delay the disturbance by elevating metabolic activity and generating functional redundancy. Especially, when nutrient was a limiting factor, the disturbance strength by DOM source was greater than that by chronic antibiotic exposure. The resistance of abundant taxa to external distributions resulted in a low explanation of community diversity, while rare taxa played key roles in response to community variation and thereby affected community assembly. Long-term SMX exposure reduced the number of key species and favored the deterministic assembly process by 21%. However, elevated community adaptability might weaken the influence of antibiotic selection. Chronic SMX exposure elevated the relative abundance of sulfonamide resistance genes (sul1, sul2) by a factor of 1.2-4.3, while that of nitrogen-fixing genes (nifH, nifK) and the metabolic pathways related to the toluene, ethylbenzene, and dioxin degradation decreased. However, the combined influence of DOM sources and multi-antibiotic exposure barely caused the difference in the genes linking to element metabolism and drug resistance of microbial communities between blank and exposed groups. This study suggested that more concern should be given to the chronic environmental effect of organic micropollutants.

Community- and genome-based evidence for a shaping influence of redox potential on bacterial protein evolution

Despite deep interest in how environments shape microbial communities, whether redox conditions influence the sequence composition of genomes is not well known. We predicted that the carbon oxidation state (ZC) of protein sequences would be positively correlated with redox potential (Eh). To test this prediction, we used taxonomic classifications for 68 publicly available 16S rRNA gene sequence data sets to estimate the abundances of archaeal and bacterial genomes in river & seawater, lake & pond, geothermal, hyperalkaline, groundwater, sediment, and soil environments. Locally, ZC of community reference proteomes (i.e., all the protein sequences in each genome, weighted by taxonomic abundances but not by protein abundances) is positively correlated with Eh corrected to pH 7 (Eh7) for the majority of data sets for bacterial communities in each type of environment, and global-scale correlations are positive for bacterial communities in all environments. In contrast, archaeal communities show approximately equal frequencies of positive and negative correlations in individual data sets, and a positive pan-environmental correlation for archaea only emerges after limiting the analysis to samples with reported oxygen concentrations. These results provide empirical evidence that geochemistry modulates genome evolution and may have distinct effects on bacteria and archaea. IMPORTANCE The identification of environmental factors that influence the elemental composition of proteins has implications for understanding microbial evolution and biogeography. Millions of years of genome evolution may provide a route for protein sequences to attain incomplete equilibrium with their chemical environment. We developed new tests of this chemical adaptation hypothesis by analyzing trends of the carbon oxidation state of community reference proteomes for microbial communities in local- and global-scale redox gradients. The results provide evidence for widespread environmental shaping of the elemental composition of protein sequences at the community level and establish a rationale for using thermodynamic models as a window into geochemical effects on microbial community assembly and evolution.

Insight into the spatiotemporal distribution of antibiotic resistance genes in estuarine sediments during long-term ecological restoration

In this study, we aimed to investigate the long-term spatiotemporal changes in hydrodynamics, antibiotics, nine typical subtypes of antibiotic resistance genes (ARGs), class 1 integron gene (intI1), and microbial communities in the sediments of a semi-enclosed estuary during ecological restoration with four treatment stages (influent (#1), effluent of the biological treatment area (#2), oxic area (#3), and plant treatment area (#4)). Ecological restoration of the estuary reduced common pollutants (nitrogen and phosphorus) in the water, whereas variations in ARGs showed noticeable seasonal and spatial features. The absolute abundance of ARGs at sampling site #2 considerably increased in autumn and winter, while it significantly increased at sampling site #3 in spring and summer. The strong intervention of biological treatment (from #1 to #2) and aerators (from #2 to #3) in the estuary substantially affected the distribution of ARGs and dominant antibiotic-resistant bacteria (ARB). The dominant ARB (Thiobacillus) in estuarine sediments may have low abundance but important dissemination roles. Meanwhile, redundancy and network analysis revealed that the microbial communities and intl1 were key factors related to ARG dissemination, which was affected by spatial and seasonal ecological restoration. A positive correlation between low flow velocity and certain ARGs (tetM, tetW, tetA, sul2, and ermC) was observed, implying that flow optimization should also be considered in future ecological restoration to remediate ARGs. Furthermore, the absolute abundance of ARGs can be utilized as an index to evaluate the removal capacity of ARGs by estuarine restoration.

Alizarin enhancement of the abundance of ARGs and impacts on the microbial community in water

Alizarin, a dyestuff from herbs, showed effective inhibition effects on pathogenic bacteria, and thus has been frequently used in the world as the main alternative to antibiotics in the treatment of inflammations and pathogen infections. However, it was unclear whether alizarin played key a role in antibiotic-induced antibiotic-resistant gene (ARG) alterations and impacted microbial community shifts in aquatic environments. In this study, the effects of alizarin or co-exposure of alizarin with antibiotics on the fate of ARGs, class 1 integron-integrase gene (intI1), and microbial populations in lake water were investigated, and the potential hosts for ARGs were analyzed. The results showed that the absolute abundance of 16s rRNA gene, ARGs (tetA, tetC, and qnrS), and intI1 were increased during the treatment of alizarin. The combination of alizarin and antibiotics was superior to alizarin in its ability to promote population growth of bacteria and induce ARGs. Additionally, alizarin more significantly altered the community composition of microorganisms in water, which resulted in differences in bacterial communities and functions.

3D porous carbon-embedded nZVI@Fe2O3 nanoarchitectures enable prominent performance and recyclability in antibiotic removal

Overcoming the instability and poor recyclability during the practical applications of contaminant scavengers is a challenging topic. Herein, a three-dimensional (3D) interconnected carbon aerogel (nZVI@Fe₂O₃/PC) embedding a core-shell nanostructure of nZVI@Fe₂O₃ was elaborately designed and fabricated via an in-situ self-assembly process. The porous carbon with 3D network architecture exhibits strong adsorption towards various antibiotic contaminants in water, where the stably embedded nZVI@Fe₂O₃ nanoparticles not only serve as magnetic seeds for recycling, but also avoid the shedding and oxidation of nZVI in the adsorption process. As a result, nZVI@Fe₂O₃/PC efficiently captures sulfamethoxazole (SMX), sulfamethazine (SMZ), ciprofloxacin (CIP), tetracycline (TC) and other antibiotics in water. In particular, an excellent adsorptive removal capacity of 329 mg g^-1 and a rapid capture kinetics (99% of removal efficiency in 10 min) under a wide pH adaptability (2-8) are achieved using nZVI@Fe₂O₃/PC as an SMX scavenger. nZVI@Fe₂O₃/PC displays exceptional long-term stability given that it shows excellent magnetic property after it is stored in water solution for 60 d, making it an ideal stable scavenger for contaminants in an etching-resistant and efficient manner. This work would also provide a general strategy to develop other stable iron-based functional architectures for efficient catalytic degradation, energy conversion and biomedicine.

Metagenomic analysis revealed sources, transmission, and health risk of antibiotic resistance genes in confluence of Fenhe, Weihe, and Yellow Rivers

Rivers are important vectors and reservoirs of antibiotics resistance genes (ARGs). Information regarding transmission and health risk of ARGs in river confluence is still lacking. In this study, metagenomics was used to distinguish contributions of human activities on ARGs and human pathogenic bacteria (HPB) in confluence of Fenhe, Weihe, and Yellow Rivers. Bacitracin resistance gene and bacA were the highest in all rivers, with 1.86 × 10^-2-7.26 × 10^-2 and 1.79 × 10^-2-9.12 × 10^-2 copies/16S rRNA copies, respectively. River confluence significantly increased the abundance of ARGs, especially at the confluence of three rivers with the highest 1.53 × 10^-1 copies/16S rRNA copies. Antibiotic efflux and antibiotic target alteration were the dominant resistant mechanisms in three rivers. ARGs profiles were influenced by multiple factors, with the contributions of various factors ranked as microbial communities > physicochemical factors > human activities > mobile genetic elements (MGEs). Notably, human activities and animal feces were important potential contributors of ARGs in the Weihe River and Yellow River. Transposons, as the main MGEs in three rivers, played important roles in ARGs transfer. The confluence of three rivers had the highest abundance of MGEs with the greatest transfer potentials, and therefore exhibiting the largest exposure risk of ARGs with 232.4 copies/cap·d. Furthermore, correlations of ARGs, MGEs, and HPB in different rivers were constructed via co-occurrence modes to systematically illustrate the health risks of ARGs. This study firstly unveiled the transmission and health risk of ARGs in river confluence, providing supports for ARGs control in watershed.

Structure and function response of bacterial communities towards antibiotic contamination in hyporheic zone sediments

Bacterial communities are crucial for processing and degrading contaminants in hyporheic zones (HZ). However, the effects of antibiotics on HZ bacterial communities have seldom been addressed. Here, using MiSeq 16S amplicon sequencing technology, the effects of acute exposure to Enrofloxacin, Sulfathiazole, Tetracycline hydrochloride, and Penicillin V potassium on HZ bacterial communities were investigated. Results revealed that HZ sediment communities responded differently to different classes of antibiotics, reflecting the distinct selection stress of antibiotics on HZ bacterial communities. Besides, HZ communities from the locations with more severe antibiotic contamination backgrounds (∼150 μg kg^-1) were more resistant towards antibiotic treatment. Compared with small/non-significant changes in HZ community diversity and composition treated with ng L^-1∼ug L^-1 level antibiotics compared to the control group, treatments with antibiotics over mg L^-1 level significantly reduced the diversity and changed the structures of HZ bacterial communities, and enhanced the resistance of the community to antibiotics by enriching antibiotic resistant bacteria. The exposure to mg L^-1 level antibiotics also changed community functions by restricting the growth of functional bacteria, such as ammonia oxidizing bacteria (AOB) Nitrosomonas, resulting in ammonia accumulation in sediments. The results implied that at field-relevant concentrations, there was no or minor effect of antibiotics on HZ bacterial community structure and functions, and only those areas with high antibiotic concentrations would have effects.

Influence of anthropogenic disturbances on antibiotic resistance gene distributions along the Minjiang River in Southeast China

River-reservoir systems have become ubiquitous among modern global aquatic environments due to the widespread construction of dams. However, little is known of antibiotic resistance gene (ARG) distributions in reservoir-river systems experiencing varying degrees of anthropogenic impacts. Here, the diversity, abundance, and spatial distribution of ARGs were comprehensively characterized along the main stem of the Minjiang River, a typical subtropic reservoir-river system in Southeast China using high-throughput quantitative PCR. A total of 252 ARG subtypes were detected from twelve sampling sites that were dominated by aac(3)-Via, followed by czcA, blaTEM, and sul1. Urban river waters (sites S9-S12) harbored more diverse ARGs than did the reservoir waters (sites S1-S7), indicating more serious antibiotic resistance pollution in areas with larger population densities. Dam construction could reduce the richness and absolute abundance of ARGs from upstream (site S7) to downstream (site S8). Urban river waters also harbored a higher proportion of mobile genetic elements (MGEs), suggesting that intensive human activities may promote ARG horizontal gene transfers. The mean relative abundance of Proteobacteria that could promote antibiotic resistance within microbial communities was also highest in urban river waters. Variance partitioning analysis indicated that MGEs and bacterial communities could explain 67.33%, 44.7%, and 90.29% of variation in selected ARGs for the entire watershed, aquaculture waters, and urban river waters, respectively. These results further suggest that urban rivers are ideal media for the acquisition and spread of ARGs. These findings provide new insights into the occurrence and potential mechanisms determining the distributions of ARGs in a reservoir-river system experiencing various anthropogenic disturbances at the watershed scale.

Evaluation of water quality at national scale from 2011 to 2021: Advances and challenges

More environmental policies and larger investments in protecting the aquatic environment in China have been made in the last decade than previously. It is important to assess how this will affect river water quality. Here, changes in water quality in China between 2011 and 2021 are assessed. Water bodies meeting class III or better defined in the Chinese Environmental Quality Standards for Surface Water (GB3838-2002) were labeled WQI, water bodies meeting class V or better but below class III were labeled WQII, and water bodies below class V were labeled WQIII. The percentage of WQI water bodies increased from 66.1 % in 2011 to 81.0 % in 2021, and the percentages of WQII and WQIII water bodies decreased between 2011 and 2021. The percentage of WQI water bodies increased more quickly and the percentage WQIII water bodies decreased more quickly after 2017 than between 2011 and 2016. The percentages of WQI water bodies in the Northwest River Basin (RB), Pearl RB, Southeast RB, Southwest RB, and Yangtze RB were >80 %, and were higher than the percentages of WQI water bodies in the other five RBs. The percentages of WQI and WQII water bodies increased but the percentage of WQIII water bodies decreased in the Hai RB. The percentage of WQI water bodies increased but the percentages of WQII and WQIII water bodies decreased in the Huai RB, Liao RB, Yangtze RB, and Yellow RB. The river monitoring capacity increased and pollution sources, particularly point sources, became more controlled, and this improved river water quality. River management in China has passed the first stage of controlling pollution sources after 10 years of centralized management. The next stage should be focused on strengthening control of non-point sources of pollution and rehabilitating ecological systems to improve river health.

Antibiotic resistance genes in Chishui River, a tributary of the Yangtze River, China: Occurrence, seasonal variation and its relationships with antibiotics, heavy metals and microbial communities

The large-scale use and release of antibiotics may create selective pressure on antibiotic resistance genes (ARGs), causing potential harm to human health. River ecosystems have long been considered repositories of antibiotics and ARGs. Therefore, the distribution characteristics and seasonal variation in antibiotics and ARGs in the surface water of the main stream and tributaries of the Chishui River were studied. The concentrations of antibiotics in the dry season and rainy season were 54.18-425.74 ng/L and 66.57-256.40 ng/L, respectively, gradually decreasing along the river direction. The results of antibiotics in the dry season and rainy season showed that livestock and poultry breeding were the main sources in the surface water of the Chishui River basin. Risk assessments indicated high risk levels of OFL in both seasons. In addition, analysis of ARGs and microbial community diversity showed that sul1 and sul3 were the main ARGs in the two seasons. The highest abundance of ARGs was 7.70 × 10⁷ copies/L, and intl1 was significantly positively correlated with all resistance genes (p< 0.01), indicating that it can significantly promote the transmission of ARGs. Proteobacteria were the dominant microorganisms in surface water, with a higher average abundance in the dry season (60.64 %) than in the rainy season (39.53 %). Finally, correlation analyses were performed between ARGs and antibiotics, microbial communities and heavy metals. The results showed that there was a significant positive correlation between ARGs and most microorganisms and heavy metals (p< 0.01), indicating that occurrence and transmission in the environment are influenced by various environmental factors and cross-selection. In conclusion, the persistent residue and transmission of ARGs and their transfer to pathogens are a great threat to human health and deserve further study and attention.

Antibiotic resistome and its driving factors in an urban river in northern China

Urban rivers dynamically interfered by anthropogenic activities are considered as a vital reservoir of antibiotic resistance genes (ARGs). Here, a total of 198 ARGs and 12 mobile genetic elements (MGEs) were profiled in water and sediment from the Chaobai river, Beijing. The total abundances of ARGs (1.01 × 10⁶-4.58 × 10⁸ copies/L in water and 2.92 × 10⁶-3.34 × 10⁹ copies/g in sediment), which were dominated by beta-lactamase genes, exhibited significant seasonal variations (p < 0.05). Significant linear correlations between the total abundances of ARGs and MGEs were observed in both water and sediment (p < 0.01). Variance partitioning analysis disclosed that environmental variables (i.e., water temperature (WT), dissolved oxygen (DO), nutrients, metals, etc.) and antibiotics were the main contributors to the variations of ARGs and MGEs, and explained 55-80 % and 27-67 % of the total variations in ARGs and MGEs, respectively. The partial least-squares path model revealed the ARG abundances in water and sediment were affected by environmental variables and antibiotics both directly and indirectly but by MGEs directly. Moreover, random forest algorithm explored that WT, Ni, DO, Co, and polyether and macrolide antibiotics were the main drivers (>10 %) of ARGs dissemination in water, whereas the transposase genes of Tp614, tnpA, and IS613 were the main drivers of ARGs dissemination in both water and sediment. This study provides a comprehensive understanding of the driving factors for the ARGs dissemination in an urban river, which is of great significance for risk management of antibiotic resistome.

Environmental risks caused by livestock and poultry farms to the soils: Comparison of swine, chicken, and cattle farms

The lack of treatment systems for pollutants in family-livestock and poultry sites results in large amounts of untreated manure and urine being directly discharged to environment. The risks from veterinary antibiotic (VA) and heavy metal (HM) exposure in the rural environment require further research. In this investigation, 221 samples (feed, manure, surface soil, soil profiles, water, and plant) were collected from 41 livestock and poultry farms (swine, chichen, and cattle). Copper (Cu), zinc (Zn), oxytetracycline (OTC), and enrofloxacin (ENR) were frequently detected in the samples. Metals and VAs in sandy loam soils were more inclined to migrate to deep layers than those in loam soils. Copper and Zn in the polluted soils mainly existed in available forms, which facilitated their migration to deep soil layers. In this study, OTC was also observed to migrate more easily to deeper soil layers than ENR due to its relatively high pKa value. Eighteen antibiotic resistance genes (ARGs) and 5 metal resistance genes (MRGs) along with one mobile genetic element (MGE) occurred in the soils at 80 cm depth. Luteimonas, Clostridium_sensu_stricto_1, and Rhodanobacter were dominant genera detected in the soil samples from different sites, which might increase migration of ARGs or degradation of VAs. An ecological risk assessment suggested that VAs posed threats to the growth of Triticum aestivum L, Cucumis sativus L, and Brassiaca chinensis L. Remediation techniques including biochar/modified biochar, anaerobic digestion, and manure composting should be developed urgently for joint HM and VA pollution.

Joint effects of bacterium and biochar in remediation of antibiotic-heavy metal contaminated soil and responses of resistance gene and microbial community

Soils containing both veterinary antibiotics (VAs) and heavy metals necessitate effective remediation approaches, and microbial and molecular levels of the results should be further examined. Here, a novel material combining waste fungus chaff-based biochar (WFCB) and Herbaspirillum huttiense (HHS1) was established to immobilize copper (Cu) and zinc (Zn) and degrade oxytetracycline (OTC) and enrofloxacin (ENR). Results showed that the combined material exhibited high immobilization of Cu (85.5%) and Zn (64.4%) and great removals of OTC (41.9%) and ENR (40.7%). Resistance genes including tet(PB), tetH, tetR, tetS, tetT, tetM, aacA/aphD, aacC, aadA9, and czcA were reduced. Abundances of potential hosts of antibiotic resistance genes (ARGs) including phylum Proteobacteria and genera Brevundimonas and Rhodanobacter were altered. Total phosphorus and pH were the factors driving the VA degrading microorganisms and potential hosts of ARGs. The combination of WFCB and HHS1 can serve as an important bioresource for immobilizing heavy metals and removing VAs in the contaminated soil.