Metagenomic insights into the profile of antibiotic resistomes in a large drinking water reservoir

Health risk ranking of antibiotic resistance genes in the Yangtze River

Antibiotic resistance is an escalating global health concern, exacerbated by the pervasive presence of antibiotic resistance genes (ARGs) in natural environments. The Yangtze River, the world's third-longest river, traversing areas with intense human activities, presents a unique ecosystem for studying the impact of these genes on human health. Here, we explored ARGs in the Yangtze River, examining 204 samples from six distinct habitats of approximately 6000 km of the river, including free-living and particle-associated settings, surface and bottom sediments, and surface and bottom bank soils. Employing shotgun sequencing, we generated an average of 13.69 Gb reads per sample. Our findings revealed a significantly higher abundance and diversity of ARGs in water-borne bacteria compared to other habitats. A notable pattern of resistome coalescence was observed within similar habitat types. In addition, we developed a framework for ranking the risk of ARG and a corresponding method for calculating the risk index. Applying them, we identified water-borne bacteria as the highest contributors to health risks, and noted an increase in ARG risks in particle-associated bacteria correlating with heightened anthropogenic activities. Further analysis using a weighted ARG risk index pinpointed the Chengdu-Chongqing and Yangtze River Delta urban agglomerations as regions of elevated health risk. These insights provide a critical new perspective on ARG health risk assessment, highlighting the urgent need for strategies to mitigate the impact of ARGs on human health and to preserve the ecological and economic sustainability of the Yangtze River for future human use.

When antibiotics encounter microplastics in aquatic environments: Interaction, combined toxicity, and risk assessments

Antibiotics and microplastics (MPs), known as emerging pollutants, are bound to coexist in aquatic environments due to their widespread distribution and prolonged persistence. To date, few systematic summaries are available for the interaction between MPs and antibiotics in aquatic ecosystems, and a comprehensive reanalysis of their combined toxicity is also needed. Based on the collected published data, we have analyzed the source and distribution of MPs and antibiotics in global aquatic environments, finding their coexistence occurs in a lot of study sites. Accordingly, the presence of MPs can directly alter the environmental behavior of antibiotics. The main influencing factors of interaction between antibiotics and MPs have been summarized in terms of the characteristics of MPs and antibiotics, as well as the environmental factors. Then, we have conducted a meta-analysis to evaluate the combined toxicity of antibiotics and MPs on aquatic organisms and the related toxicity indicators, suggesting a significant adverse effect on algae, and inapparent on fish and daphnia. Finally, the environmental risk assessments for antibiotics and MPs were discussed, but unfortunately the standardized methodology for the risk assessment of MPs is still challenging, let alone assessment for their combined toxicity. This review provides insights into the interactions and environment risks of antibiotics and MPs in the aquatic environment, and suggests perspectives for future research.

Spread of antibiotic resistance genes in drinking water reservoirs: Insights from a deep metagenomic study using a curated database

The exploration of antibiotic resistance genes (ARGs) in drinking water reservoirs is an emerging field. Using a curated database, we enhanced the ARG detection and conducted a comprehensive analysis using 2.2 Tb of deep metagenomic sequencing data to determine the distribution of ARGs across 16 drinking water reservoirs and associated environments. Our findings reveal a greater diversity of ARGs in sediments than in water, underscoring the importance of extensive background surveys. Crucial ARG carriers-specifically Acinetobacter, Pseudomonas, and Mycobacterium were identified in drinking water reservoirs. Extensive analysis of the data uncovered a considerable concern for drinking water safety, particularly in regions reliant on river sources. Mobile genetic elements have been found to contribute markedly to the propagation of ARGs. The results of this research suggest that the establishment of drinking water reservoirs for supplying raw water may be an effective strategy for alleviating the spread of water-mediated ARGs.

Profiles of phage in global hospital wastewater: Association with microbial hosts, antibiotic resistance genes, metal resistance genes, and mobile genetic elements

Hospital wastewater (HWW) is known to host taxonomically diverse microbial communities, yet limited information is available on the phages infecting these microorganisms. To fill this knowledge gap, we conducted an in-depth analysis using 377 publicly available HWW metagenomic datasets from 16 countries across 4 continents in the NCBI SRA database to elucidate phage-host dynamics and phage contributions to resistance gene transmission. We first assembled a metagenomic HWW phage catalog comprising 13,812 phage operational taxonomic units (pOTUs). The majority of these pOTUs belonged to the Caudoviricetes order, representing 75.29 % of this catalog. Based on the lifestyle of phages, we found that potentially virulent phages predominated in HWW. Specifically, 583 pOTUs have been predicted to have the capability to lyse 81 potentially pathogenic bacteria, suggesting the promising role of HWW phages as a viable alternative to antibiotics. Among all pOTUs, 1.56 % of pOTUs carry 108 subtypes of antibiotic resistance genes (ARGs), 0.96 % of pOTUs carry 76 subtypes of metal resistance genes (MRGs), and 0.96 % of pOTUs carry 22 subtypes of non-phage mobile genetic elements (MGEs). Predictions indicate that certain phages carrying ARGs, MRGs, and non-phage MGEs could infect bacteria hosts, even potential pathogens. This suggests that phages in HWW may contribute to the dissemination of resistance-associated genes in the environment. This meta-analysis provides the first global catalog of HWW phages, revealing their correlations with microbial hosts and pahge-associated ARGs, MRG, and non-phage MGEs. The insights gained from this research hold promise for advancing the applications of phages in medical and industrial contexts.

The relationship between water quality and the microbial virulome and resistome in urban streams in Brazil

Urban streams that receive untreated domestic and hospital waste can transmit infectious diseases and spread drug residues, including antimicrobials, which can then increase the selection of antimicrobial-resistant bacteria. Here, water samples were collected from three different urban streams in the state of São Paulo, Brazil, to relate their range of Water Quality Indices (WQIs) to the diversity and composition of aquatic microbial taxa, virulence genes (virulome), and antimicrobial resistance determinants (resistome), all assessed using untargeted metagenome sequencing. There was a predominance of phyla Proteobacteria, Actinobacteria, and Bacteroidetes in all samples, and Pseudomonas was the most abundant detected genus. Virulence genes associated with motility, adherence, and secretion systems were highly abundant and mainly associated with Pseudomonas aeruginosa. Furthermore, some opportunistic pathogenic genera had negative correlations with WQI. Many clinically relevant antimicrobial resistance genes (ARGs) and efflux pump-encoding genes that confer resistance to critically important antimicrobials were detected. The highest relative abundances of ARGs were β-lactams and macrolide-lincosamide-streptogramin. No statistically supported relationship was detected between the abundance of virulome/resistome and collection type/WQI. On the other hand, total solids were a weak predictor of gene abundance patterns. These results provide insights into various microbial outcomes given urban stream quality and point to its ecological complexity. In addition, this study suggests potential consequences for human health as mediated by aquatic microbial communities responding to typical urban outputs.

Ranking the risk of antibiotic resistance genes by metagenomic and multifactorial analysis in hospital wastewater systems

Antimicrobial resistance poses a serious threat to human health. Hospital wastewater system (HWS) is an important source of antibiotic resistance genes (ARGs). The risk of ARGs in HWS is still an under-researched area. In this study, we collected publicly metagenomic datasets of 71 hospital wastewater samples from 18 hospitals in 13 cities. A total of 9838 contigs were identified to carry 383 unique ARGs across all samples, of which 2946 contigs were plasmid-like sequences. Concurrently, the primary hosts of ARGs within HWS were found to be Escherichia coli and Klebsiella pneumoniae. To further evaluate the risk of each ARG subtype, we proposed a risk assessment framework based on the importance of corresponding antibiotics as defined by the WHO and three other indicators - ARG abundance (A), mobility (M), and host pathogenicity (P). Ninety ARGs were identified as R1 ARGs having high-risk scores, which meant having a high abundance, high mobility, and carried by pathogens in HWS. Furthermore, 25% to 49% of genomes from critically important pathogens accessed from NCBI carried R1 ARGs. A significantly higher number of R1 ARGs was carried by pathogens in the effluents of municipal wastewater treatment plants from NCBI, highlighting the role of R1 ARGS in accelerating health and environmental risks.

Historical trajectories of antibiotics resistance genes assessed through sedimentary DNA analysis of a subtropical eutrophic lake

Investigating the occurrence of antibiotic-resistance genes (ARGs) in sedimentary archives provides opportunities for reconstructing the distribution and dissemination of historical (i.e., non-anthropogenic origin) ARGs. Although ARGs in freshwater environments have attracted great attention, historical variations in the diversity and abundance of ARGs over centuries to millennia remain largely unknown. In this study, we investigated the vertical change patterns of bacterial communities, ARGs and mobile genetic elements (MGEs) found in sediments of Lake Chenghai spanning the past 600 years. Within resistome preserved in sediments, 177 ARGs subtypes were found with aminoglycosides and multidrug resistance being the most abundant. The ARG abundance in the upper sediment layers (equivalent to the post-antibiotic era since the 1940s) was lower than those during the pre-antibiotic era, whereas the ARG diversity was higher during the post-antibiotic era, possibly because human-induced lake eutrophication over the recent decades facilitated the spread and proliferation of drug-resistant bacteria. Statistical analysis suggested that MGEs abundance and the bacterial community structure were significantly correlated with the abundance and diversity of ARGs, suggesting that the occurrence and distribution of ARGs may be transferred between different bacteria by MGEs. Our results provide new perspectives on the natural history of ARGs in freshwater environments and are essential for understanding the temporal dynamics and dissemination of ARGs.

Thermokarst lakes are hotspots of antibiotic resistance genes in permafrost regions on the Qinghai-Tibet Plateau

Antibiotic resistance genes (ARGs) are natural products and emerging pollutants in remote environments, including permafrost regions that are rapidly thawing due to climate warming. We investigated the role of thermokarst lakes (including sediment and water) in reserving ARGs compared to permafrost soils across the permafrost regions on the Qinghai-Tibet Plateau. As intrinsically connected distinct environments, permafrost soil, lake sediment, and lake water harbored 1239 ARGs in total, while a considerable number of same ARGs (683 out of 1239) concurrently presented in all these environments. Soil and sediment had a higher number of ARGs than water. Multidrug resistance genes were the most diverse and abundant in all three environments, where cls, ropB, mdfA, fabI, and macB were the top five most abundant ARGs while with different orders. Soil and sediment had similar ARG profiles, and the alpha and beta diversity of ARGs in sediment were positively correlated with that in soil. The beta diversity of ARG profiles between sediment and soil was highly contributed by turnover component (89%). However, turnover and nestedness components were almost equality contributed (46%-54%) to the beta diversity of ARG profiles between soil and water as well as between sediment and water. The results suggested that thermokarst lake sediments might inherit the ARGs in permafrost soils. Water ARGs are the subset of soil ARGs and sediment ARGs to a certain degree with species turnover playing a significant role. When accounting the ARGs in sediment and water together, thermokarst lakes had a significantly higher number of ARGs than permafrost soils, suggesting that thermokarst lakes act as the hotspots of ARGs in permafrost regions. These findings are disturbing especially due to the fact that tremendous number of thermokarst lakes are forming under accelerating climate change.

Prevalence of antibiotic and metal resistance genes in phytoremediated cadmium and zinc contaminated soil assisted by chitosan and Trichoderma harzianum

Heavy metal in soil have been shown to be toxic with high concentrations and acts as selective pressure on both bacterial metal and antibiotic resistance determinants, posing a serious risk to public health. In cadmium (Cd) and zinc (Zn) contaminated soil, chitosan (Chi) and Trichoderma harzianum (Tri) were applied alone and in combination to assist phytoremediation by Amaranthus hypochondriacus L. Prevalence of antibiotic and metal resistance genes (ARGs and MRGs) in the soil was also evaluated using metagenomic approach. Results indicated that the phytoremediation of Cd and Zn contaminated soil was promoted by Chi, and Tri further reinforced this effect, along with the increased availability of Cd and Zn in soil. Meanwhile, combination of Chi and Tri enhanced the prevalence of ARGs (e.g., multidrug and β-lactam resistance genes) and maintained a high level of MRGs (e.g., chromium, copper) in soil. Soil available Zn and Cd fractions were the main factors contributing to ARGs profile by co-selection, while boosted bacterial hosts (e.g., Mitsuaria, Solirubrobacter, Ramlibacter) contributed to prevalence of most MRGs (e.g., Cd). These findings indicate the potential risk of ARGs and MRGs propagation in phytoremediation of metal contaminated soils assisted by organic and biological agents.

Impacts of net cages on pollutant accumulation and its consequence on antibiotic resistance genes (ARGs) dissemination in freshwater ecosystems: Insights for sustainable urban water management

There has been increasing interest in the role of human activities in disseminating antibiotic-resistance genes (ARGs) in aquatic ecosystems. However, the influence of pollutant accumulation on anthropogenic pollutant-ARG synergistic actions is limited. This study explored the association of net cages with the propagation of anthropogenic pollutants and their consequences for influencing the enrichment of ARGs using high-throughput metagenomic sequencing. We showed that net cages could substantially impact the ecology of freshwater systems by enhancing i) ARG diversity and the tendency for ARG-horizontal gene transfer and ii) the overlap of mobile genetic elements (MGEs) with biocide-metal resistance genes (BMRGs) and ARGs. These findings suggested that the cotransfer of these three genetic determinants would be favored in net cage plots and that nonantibiotic factors such as metal(loid)s, particularly iron (Fe), displayed robust selective pressures on ARGs exerted by the net cage. The resistome risk scores of net cage sediments and biofilms were higher than those from off-net cage plots, indicating that the net cage-origin antibiotic resistome should be of great concern. The combination of deterministic and stochastic processes acting on bacterial communities could explain the higher ARG variations in cage plots (8.2%) than in off-cage plots (3.4%). Moreover, MGEs and pollutants together explained 43.3% of the total variation in ARG communities, which was higher than that of off-cage plots (8.8%), considering pollutants, environmental variables, MGEs, and assembly processes. These findings will inform the development of policies and guidelines to more effectively limit the spread of antimicrobial resistance and achieve the goal of sustainability in freshwater systems in urban areas.

Occurrence, fate and potential health risks of antibiotic resistomes in a constructed wetlands-reservoir ecosystem for drinking water source improvement

The development of effective and feasible engineering technologies to control the transmission of antibiotic resistance genes (ARGs) and pathogenic antibiotic-resistant bacteria (PARB) form drinking water sources is urgently needed for ensuring drinking water safety. In this study, metagenomic analysis was applied to systematically explore the full profiles, removal, and potential health risks of antibiotic resistomes in a large constructed wetlands-reservoir ecosystem (CWs-R) for drinking water source improvement. A total of 343 ARG subtypes belonging to 18 ARG types were identified from water and sediment samples in the CWs-R ecosystem, with an average abundance of 0.339 copies/cell, and bacitracin and multidrug resistance genes were the predominant ARG types in the water and sediment, respectively. The CWs-R ecosystem showed an excellent removal efficiency of ARGs and mobile genetic elements (MGEs) in water, with the total removal rate reaching 64.82 % and 77.09 %, respectively, among which the emergent plant zone and ecological storage unit played major roles. The metagenomic assembly tracked many mobile ARGs and opportunistic pathogens in the CWs-R ecosystem and identified 19 contigs as ARG-carrying pathogens, including Staphylococcus aureus, Salmonella enterica, Escherichia coli, and Klebsiella pneumonia. Overall, the CWs-R ecosystem has an important role in reducing the potential public health risks posed by antibiotic resistomes in drinking water sources but still cannot fully eliminate them. Therefore, we further classified water and sediment samples in the CWs-R ecosystem and identified potential ARGs and PARB indicators based on the metagenomic analysis results by considering the potential for horizontal transfer of ARGs to opportunistic pathogens. Taken together, this work demonstrates the CWs-R ecosystem as an economical and feasible engineering technology to reduce the dissemination of antibiotic resistomes in the drinking water source, provides useful information for monitoring and controlling antibiotic resistance in similar water sources, and ensures biosafety of drinking water.

Antibiotic resistant genes profile in the surface water of subtropical drinking water river-reservoir system

To comprehensively understand antibiotic resistant genes (ARGs) profile in the subtropical drinking water river-reservoir system, this study selected Dongzhen river-reservoir system in Mulan Creek as object to investigate the spatial-temporal characteristics of ARGs diversity, bacterial host and resistance mechanism, and to analyze the key environmental factors driving ARGs profile variation. The results indicated that a total of 440 ARGs were detected in the target system, and the ARGs distribution pattern in the reservoir was attributed to autologous evolution or the comprehensive influence of feeding river system. The predominant bacterial host at different sites showed similar variations to dominated ARGs, and Proteobacteria, Actinobacteria and Bacteroidetes harbored most ARGs at phylum level, which showed the highest proportions of 74%, 37% and 35%, respectively. Antibiotic efflux was the primary resistance mechanism in all samples from wet season (45%-60%), yet the samples from dry season exhibited multiple resistance mechanisms, including inactivation (37%-52%), efflux (44%), and target alteration (43%). The total relative abundances of ARGs in the target system ranged from 0.89 × 10-2 to 1.71 × 10-2, and seasonal variation had a more significant influence on ARGs abundance than spatial variation (R = 0.68, P < 0.01). Environmental factors analysis indicated that the concentrations of nitrite nitrogen and total organic carbon were significant factors explaining ARGs number and various resistance mechanism proportions (P < 0.01), accounting for 48.7% and 61.1% of the variation, respectively; ammonia nitrogen concentration, total organic carbon concentration, temperature and pH were the significant influence factors on the relative abundance of ARGs (P < 0.05), with standardized regression weights of 0.700, 1.414, 1.447, and 1.727, respectively. In summary, in the surface water of the target system, ARGs diversity was primarily driven by ARGs horizontal transfer and antibiotics biosynthesis. Nutrients mainly promoted ARGs abundance by providing abundant energy, rather than increasing bacterial reproductive capacity.

A novel culture-enriched metagenomic sequencing strategy effectively guarantee the microbial safety of drinking water by uncovering the low abundance pathogens

Assessing the presence of waterborne pathogens and antibiotic resistance genes (ARGs) is crucial for managing the environmental quality of drinking water sources. However, detecting low abundance pathogens in such settings is challenging. In this study, a workflow was developed to enrich for broad spectrum pathogens from drinking water samples. A mock community was used to evaluate the effectiveness of various enrichment broths in detecting low-abundance pathogens. Monthly metagenomic surveillance was conducted in a drinking water source from May to September 2021, and water samples were subjected to five enrichment procedures for 6 h to recover the majority of waterborne bacterial pathogens. Oxford Nanopore Technology (ONT) was used for metagenomic sequencing of enriched samples to obtain high-quality pathogen genomes. The results showed that selective enrichment significantly increased the proportions of targeted bacterial pathogens. Compared to direct metagenomic sequencing of untreated water samples, targeted enrichment followed by ONT sequencing significantly improved the detection of waterborne pathogens and the quality of metagenome-assembled genomes (MAGs). Eighty-six high-quality MAGs, including 70 pathogen MAGs, were obtained from ONT sequencing, while only 12 MAGs representing 10 species were obtained from direct metagenomic sequencing of untreated water samples. In addition, ONT sequencing improved the recovery of mobile genetic elements and the accuracy of phylogenetic analysis. This study highlights the urgent need for efficient methodologies to detect and manage microbial risks in drinking water sources. The developed workflow provides a cost-effective approach for environmental management of drinking water sources with microbial risks. The study also uncovered pathogens that were not detected by traditional methods, thereby advancing microbial risk management of drinking water sources.

Wastewater treatment plant effluents exert different impacts on antibiotic resistome in water and sediment of the receiving river: Metagenomic analysis and risk assessment

Wastewater treatment plants (WWTPs) are considered as hotspots for the spread of antibiotic resistome into the environment. However, the differential contributions of WWTPs to the antibiotic resistome in the receiving river water and sediment are poorly understood. Here, based on metagenomic analysis, we found that the WWTP effluents significantly elevated the diversities and abundances of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) in the receiving river water from the Qinghai-Tibet Plateau, but showed less interference with the antibiotic resistome in sediment. Estimated by SourceTracker, WWTPs contributed 60.691.8% of ARGs in downstream river water, much higher than those for sediment (7.7568.0%). A holistic comparison of ARG risks based on analysis of ARG combination, mobility risk, ARG hosts and ARG-carrying pathogens further revealed the great impacts of WWTP effluents on downstream river water rather than sediment. Among various MGEs, tnpA exhibited the greatest potential for the dissemination of ARGs, and displayed highest co-occurrence frequency with multiple ARGs. P. aeruginosa, E. cloacae, and E. coli were identified as the critical-priority pathogens of ARG hosts. This study demonstrated the much greater impacts of WWTP effluents on the downstream water compared with sediment, which is significant for developing effective strategies to mitigate ARG risks.

Unraveling antibiotic resistomes associated with bacterial and viral communities in intertidal mudflat aquaculture area

The extensive use of antibiotics in intertidal mudflat aquaculture area has substantially increased the dissemination risk of antibiotic resistance genes (ARGs). As hosts of ARGs, bacteria and virus exert vital effects on ARG dissemination. However, the insights for the interrelationships among ARGs, bacteria, and virus have not been thoroughly explored in intertidal mudflat. Therefore, this study attempts to unravel the occurrence, dissemination, evolution, and driving mechanisms of ARGs associated with bacterial and viral communities using metagenomic sequencing in a typical intertidal mudflat. Abundant and diverse ARGs (22 types and 437 subtypes) were identified and those of ARGs were higher in spring than in autumn. It is worthy noted that virus occupied a more essential position than bacteria for ARGs dissemination through network analysis. Meanwhile, nitrogen exerted indirect effect on ARG profiles by shaping viral and bacterial diversity. According to the results of neutral and null models, deterministic processes dominated the ARG community assembly by controlling sediment nitrogen and antibiotics. Homogeneous and variable selection dominated phylogenetic turnover of ARG community, contributing 46.15% and 45.90% of the total processes, respectively. This study can hence theoretically support for the ARG pollution control and management in intertidal mudflat aquaculture area.

Occurrence, distribution and their correlation with different parameters of antibiotics and antibiotic resistance genes in lakes of China: A review

The exposure of antibiotics and antibiotic resistance genes (ARGs) as potential threats to the environment has raised global concern. This study provides discussion on the emergence and distribution of antibiotics and ARGs in lakes. The correlation of critical water quality parameters with antibiotics and ARGs are evaluated along with their integrative potential ecological risk. Sulfonamides (∼67.18 ng/L) and quinolones (∼77.62 ng/L) were the dominant antibiotics distributed in the aqueous phase, while the quinolones and tetracyclines were the primary contamination factors in the sediment phase. The temporal and spatial distribution revealed that the antibiotic concentrations were significantly lower in summer than other seasons and the lakes in Hebei and Jiangsu provinces exhibited the highest antibiotic pollution. The detection frequency and relative abundance of sul1 gene have been the highest among all detected ARGs. Moreover, ARGs in lakes were driven by several factors, with bacterial communities and mobile genetic elements that prevailed the positive distribution of ARGs. Antibiotics have been identified as critical factors in inducing the propagation of ARGs, which could be further enhanced by chemical contaminants (e.g., heavy metals and nutrients). Involving the risk assessment strategies, research attention should be paid on three antibiotics (ofloxacin, sulfamethoxazole and erythromycin) to strengthen the policy and management of Baiyangdian Lake and East Dongting Lake. This review analysis will provide in-depth understanding to the researchers and policy-makers in formulation of strategies for remediation of antibiotic contamination in the lakes.

Occurrence and prevalence of antibiotic resistance genes and pathogens in an industrial park wastewater treatment plant

Antibiotic resistance genes (ARGs) and pathogens are emerging environmental pollutants that pose a threat to human health and ecosystem. Industrial park wastewater treatment plants (WWTPs) treat large amounts of comprehensive wastewater derived from industrial production and park human activity, which is possible a source of ARGs and pathogens. Therefore, this study investigated the occurrence and prevalence of ARGs, ARGs hosts and pathogens and assesses the ARGs health risk in the biological treatment process in a large-sale industrial park WWTP using metagenomic analysis and omics-based framework, respectively. Results show that the major ARG subtypes are multidrug resistance genes (MDRGs), macB, tetA(58), evgS, novA, msbA and bcrA and the ARGs main hosts were genus Acidovorax, Pseudomonas, Mesorhizobium. In particular, all determined ARGs genus level hosts are pathogens. The total removal percentage of ARGs, MDRGs and pathogens were 12.77 %, 12.96 % and 25.71 % respectively, suggesting that the present treatment could not efficiently remove these pollutants. The relative abundance of ARGs, MDRGs and pathogens varied along biological treatment process that ARGs and MDRGs were enriched in activated sludge and pathogens were enriched in both secondary sedimentation tank and activated sludge. Among 980 known ARGs, 23 ARGs (e.g., ermB, gadX and tetM) were assigned into risk Rank I with characters of enrichment in the human-associated environment, gene mobility and pathogenicity. The results indicate that industrial park WWTPs might serve as an important source of ARGs, MDRGs, and pathogens. These observations invite further study of the origination, development, dissemination and risk assessment of industrial park WWTPs ARGs and pathogens.

A Review of Antibiotics, Antibiotic Resistant Bacteria, and Resistance Genes in Aquaculture: Occurrence, Contamination, and Transmission

Antibiotics are commonly used to prevent and control diseases in aquaculture. However, long-term/overuse of antibiotics not only leaves residues but results in the development of antibiotic resistant bacteria (ARB) and antibiotic resistance genes (ARGs). Antibiotics, ARB, and ARGs are widespread in aquaculture ecosystems. However, their impacts and interaction mechanisms in biotic and abiotic media remain to be clarified. In this paper, we summarized the detection methods, present status, and transfer mechanisms of antibiotics, ARB, and ARGs in water, sediment, and aquaculture organisms. Currently, the dominant methods of detecting antibiotics, ARB, and ARGs are UPLC-MS/MS, 16S rRNA sequencing, and metagenomics, respectively. Tetracyclines, macrolides, fluoroquinolones, and sulfonamides are most frequently detected in aquaculture. Generally, antibiotic concentrations and ARG abundance in sediment are much higher than those in water. Yet, no obvious patterns in the category of antibiotics or ARB are present in organisms or the environment. The key mechanisms of resistance to antibiotics in bacteria include reducing the cell membrane permeability, enhancing antibiotic efflux, and structural changes in antibiotic target proteins. Moreover, horizontal transfer is a major pathway for ARGs transfer, including conjugation, transformation, transduction, and vesiculation. Identifying, quantifying, and summarizing the interactions and transmission mechanisms of antibiotics, ARGs, and ARB would provide useful information for future disease diagnosis and scientific management in aquaculture.

Insight into the dynamics of drinking water resistome in China

Antibiotic resistance (AR) is a serious environmental hazard of the current age. Antibiotic resistance genes (ARGs) are the fundamental entities that spread AR in the environment. ARGs are likely to be transferred from the non-pathogenic to pathogenic microbes that might ultimately be responsible for the AR in humans and other organisms. Drinking water (DW) is the primary interaction route between ARGs and humans. Being the highest producer and consumer of antibiotics China poses a potential threat to developing superbugs and ARGs dissemination. Herein, we comprehensively seek to review the ARGs from dominant DW sources in China. Furthermore, the origin and influencing factors of the ARGs to the DW in China have been evaluated. Commonly used methods, both classical and modern, are being compiled. In addition, the risk posed and mitigation strategies of DW ARGs in China have been outlined. Overall, we believe this review would contribute to the assessment of ARGs in DW of China and their dissemination to humans and other animals and ultimately help the policymakers and scientists in the field to counteract this problem on an emergency basis.

Evaluation of environmental factors and microbial community structure in an important drinking-water reservoir across seasons

The composition of microbial communities varies in water and sediments, and changes in environmental factors have major effects on microbiomes. Here, we characterized variations in microbial communities and physicochemical factors at two sites in a large subtropical drinking water reservoir in southern China. The microbiomes of all sites, including the diversity and abundance of microbial species, were determined via metagenomics, and the relationships between microbiomes and physicochemical factors were determined via redundancy analysis. The dominant species in sediment and water samples differed; Dinobryon sp. LO226KS and Dinobryon divergens were dominant in sediment samples, whereas Candidatus Fonsibacter ubiquis and Microcystis elabens were dominant in water. The diversity was also significantly different in microbial alpha diversity between water and sediment habitats (p < 0.01). The trophic level index (TLI) was the major factor affecting the microbial community in water samples; Mycolicibacterium litorale and Mycolicibacterium phlei were significantly positively related to TLI. Furthermore, we also studied the distribution of algal toxin-encoding genes and antibiotic-resistant genes (ARGs) in the reservoir. It found that water samples contained more phycotoxin genes, with the cylindrospermopsin gene cluster most abundant. We found three genera highly related to cylindrospermopsin and explored a new cyanobacteria Aphanocapsa montana that may produce cylindrospermopsin based on the correlation through network analysis. The multidrug resistance gene was the most abundant ARG, while the relationship between ARGs and bacteria in sediment samples was more complicated than in water. The results of this study enhance our understanding of the effects of environmental factors on microbiomes. In conclusion, research on the properties, including profiles of algal toxin-encoding genes and ARGs, and microbial communities can aid water quality monitoring and conservation.

The reasons for the spatial and media distribution variations of ARGs in a typical semi-enclosed bay

Antibiotic resistance genes (ARGs) are considered a newly emerging contaminant. This study aimed to investigate the spatial and media distribution patterns of ARGs in Jiaozhou Bay, as well as the reasons behind these patterns. The results revealed that aminoglycoside and MLSB resistant genes predominated in all samples, and the relative abundance of ARGs ranging from 10^-6 to 10^-2, 10^-6 to 10^-3 and 10^-5 to 10^-2 copies/16S rRNA in coastal water, bay water, and sediments, respectively. The significant spatial variation of ARGs was explained by the fact that the coastal water was more susceptible to human activities, whereas environmental physicochemical factors played a crucial role in the bay water. The intrinsic reason for the media distribution variation was the different assembly processes in the two media, while the external reason was that the ARGs in the water and sediments were mainly influenced by environmental physicochemical factors and heavy metals, respectively.

The distinct resistance mechanisms of cyanobacteria and green algae to sulfamethoxazole and its implications for environmental risk assessment

Cyanobacteria and green algae are the OECD recommended test organisms for environmental toxicity assessments of chemicals. Whether the differences in these two species' responses to the identical chemical affect the assessment outcomes is a question worth investigating. Firstly, we investigated the distinct resistance mechanisms of Synechococcus sp. (cyanobacteria) and R. subcapitata (green algae) to sulfamethoxazole (SMX). The antioxidant system analysis demonstrated that R. subcapitata mainly relies on enhancing the activity of first line defense antioxidants, including superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GPx), which is the most powerful and efficient response to get rid of ROS, whereas Synechococcus sp. depends upon increasing the activity of glutathione-S-transferase (GST) and GPx to resist oxidative stress. Besides, a total 7 transformation products (TPs) of SMX were identified in R. subcapitata culture medium. The analysis of conjectural transformation pathways and the predicted toxicity indicates that R. subcapitata could relieve SMX toxicity by degrading it to low eco-toxic TPs. Additionally, we summarized numerous exposure data and assessed the environmental risk of various antibiotics, revealing an inconsistent result for the same type of antibiotic by using cyanobacteria and green algae, which is most likely due to the different resistance mechanisms. In the future, modified indicators or comprehensive assessment methods should be considered to improve the rationality of environmental toxicity assessments.

Explaining the resistomes in a megacity's water supply catchment: Roles of microbial assembly-dominant taxa, niched environments and pathogenic bacteria

Antibiotic resistance genes (ARGs) in drinking water sources suggest the possible presence of resistant microorganisms that jeopardize human health. However, explanations for the presence of specific ARGs in situ are largely unknown, especially how their prevalence is affected by local microbial ecology, taxa assembly and community-wide gene transfer. Here, we characterized resistomes and bacterial communities in the Taipu River catchment, which feeds a key drinking water reservoir to a global megacity, Shanghai. Overall, ARG abundances decreased significantly as the river flowed downstream towards the reservoir (P < 0.01), whereas the waterborne bacteria assembled deterministically (|βNRI| > 2.0) as a function of temperature and dissolved oxygen conditions with the assembly-dominant taxa (e.g. Ilumatobacteraceae and Cyanobiaceae) defining local resistomes (P < 0.01, Cohen's D = 4.22). Bacterial hosts of intragenomic ARGs stayed at the same level across the catchment (60 ∼ 70 genome copies per million reads). Among them, the putative resistant pathogens (e.g. Burkholderiaceae) carried mixtures of ARGs that exhibited high transmission probability (transfer counts = 126, P < 0.001), especially with the microbial assembly-dominant taxa. These putative resistant pathogens had densities ranging form 3.0 to 4.0 × 10⁶ cell/L, which was more pronouncedly affected by resistome and microbial assembly structures than environmental factors (SEM, std-coeff β = 0.62 vs. 0.12). This work shows that microbial assembly and resistant pathogens play predominant roles in prevelance and dissemination of resistomes in receiving water, which deserves greater attention in devisng control strategies for reducing in-situ ARGs and resistant strains in a catchment.

Antibiotic resistance genes in landfill leachates from seven municipal solid waste landfills: Seasonal variations, hosts, and risk assessment

Landfills are reservoir of antibiotics and antibiotic resistance. Antibiotic resistance would transport to the environment through landfill leachate, posing threaten to the environment. However, long term monitoring on antibiotic resistance genes in landfill leachate transportation is limited. Furthermore, antibiotic resistance gene hosts and their risk assessment are lacking. In this study, we investigated the seasonal variation of ARGs sulI, tetO and tetW in seven Chinese municipal solid waste landfill leachates over two years (2017-2018) by quantitative polymerase chain reaction. We also evaluated the associated bacterial hosts and their risk levels based on metagenomics and omics-based framework for assessing the health risk of antimicrobial resistance genes, respectively. Because sulI, tetO and tetW are abundant and the most frequently detected ARGs in global landfill system, they are selected as target ARGs. Results showed that the relative content of target ARGs in 2017 was 100 times higher than that in 2018, suggesting ARGs attenuation. The hosts of sulI were phyla of Lentisphaerae and Proteobacteria, whereas the hosts of tetO and tetW were Bacteroidetes and Firmicutes. Remarkably, the host species include pathogenic bacterium (Salmonella enterica, Labilibaculum filiforme, Bacteroidales bacterium, Anaeromassilibacillus senegalensis, and Pseudochrobactrum sp. B5). ARGs tetO and tetW belong to the Rank II level with characters of enrichment in the human-associated environment and gene mobility, and sulI ranked as Rank VI. In addition, among 1210 known ARGs in the landfill leachate, 78 ARGs belonged to risk Rank I (enrichment in human-associated environment, gene mobility and pathogenicity), demonstrating high health risk of landfill system. These results demonstrate that antibiotic resistance in landfill and landfill leachate have high health risk and the kind of ARGs with high abundance in human-associated environment, gene mobility and pathogenicity should be paid more attention.

Deterministic assembly process dominates bacterial antibiotic resistome in wastewater effluents receiving river

Antibiotic resistance has become a concerning global health challenge, such as the dissemination of bacteria and genes between humans and the environments. Wastewater treatment plants (WWTPs) effluents, as significant reservoirs for antimicrobial resistant bacteria and antibiotic resistance genes (ARGs), pose critical risks to public health. However, whether wastewater effluent prominently contributes to the abundance of ARGs and their community assembly processes in receiving river has yet been unclear. Here we investigated the effects of the effluent discharge on the ARGs and their associate microbial community in the receiving river (Qinhuai River, Nanjing) of upstream and 2000 m downstream of one WWTPs discharge point. Results revealed that the total antibiotic concentrations of all sediment samples ranged from 37.86 to 76.11 µg/kg dw, while antibiotic concentrations and ARG abundances in the river near the wastewater discharge site were significantly higher than that of the downstream receiving river. The metagenomic assembly obtained 245 ARGs associated with 19 antibiotic types in the receiving river. Network analyses confirmed that Proteobacteria, Firmicutes, Acidobacteria, and Bacteroides were the key phylum and positively correlated with the antibiotic resistome. Additionally, the bacterial pathogens of the receiving river were identified as the most frequent strains of clinically relevant antibacterial resistance, such as Streptococcus pneumoniae, Pseudomonas aeruginosa, and Acinetobacter baumannii. Using null modeling analysis to determine the quantification of ecological processes, the results showed that heterogeneous environmental selection (81.81%) was a dominate role of the ecological mechanisms determining the ARG community reconstruction in the receiving river. Our results may contribute to control the environmental dissemination of antimicrobial resistance risks in aquatic environments.

Study on microbes and antibiotic resistance genes in karst primitive mountain marshes - A case study of Niangniang Mountain in Guizhou, China

Previous research on antibiotic resistance genes and microorganisms centered on those in urban sewage treatment plants, breeding farms, hospitals and others with serious antibiotic pollution. However, at present, there are evident proofs that antibiotic resistance genes (ARGs) indeed exist in a primitive environment hardly without any human's footprints. Accordingly, an original karst mountain swamp ecosystem in Niangniang Mountain, Guizhou, China, including herbaceous swamp, shrub swamp, sphagnum bog and forest swamp, was selected to analyze the physical and chemical parameters of sediments. Moreover, microbial compositions, functions, as well as their connections with ARGs were assayed and analyzed using metagenomic technology. The results showed that there was no significant difference in the dominant microorganisms and ARGs in the four marshes, in which the dominant bacteria phyla were Proteobacteria (37.82 %), Acidobacteriota (22.17 %) and Actinobacteriota (20.64 %); the dominant archaea Euryarchaeota. (1.00 %); and the dominant eukaryotes Ascomycota (0.07 %), with metabolism as their major functions. Based on the ARDB database, the number of ARGs annotated reached 209 including 30 subtypes, and the dominant ARGs were all Bacitracin resistance genes (bacA, 84.77 %). In terms of the diversity of microorganisms and ARGs, the herbaceous swamp ranked the top, and the shrub swamp were at the bottom. Correlation analysis between microorganisms and resistance genes showed that, apart from aac2ic, macB, smeE, tetQ, and tetL, other ARGs were positively correlated with microorganisms. Among them, baca coexisted with microorganisms. Pearson correlation analysis results showed that contrary to ARGs, microorganisms were more affected by environmental factors.

Antibiotic Toxicity Isolated and as Binary Mixture to Freshwater Algae Raphidocelis subcapitata: Growth Inhibition, Prediction Model, and Environmental Risk Assessment

Antibiotics in aqueous environments can have extremely adverse effects on non-targeted organisms. However, many research projects have only focused on the toxicological evaluation of individual antibiotics in various environments. In the present work, individual and binary mixture toxicity experiments have been conducted with the model organism Raphidocelis subcapitata (R. subcapitata), and a mixture concentration-response curve was established and contrasted with the estimated effects on the basis of both the concentration addition (CA) and the independent action (IA) models. In addition, different risk assessment methods were used and compared to evaluate the environmental risk of binary mixtures. The toxic ranking of the selected antibiotics to R. subcapitata was erythromycin (ERY) > sulfamethoxazole (SMX) > sulfamethazine (SMZ). In general, the conclusion of this study is that the adverse effects of binary mixtures are higher than the individual antibiotics. The CA model and RQSTU are more suitable for toxicity prediction and risk assessment of binary mixtures. This study reveals the potential ecological risks that antibiotics and their mixtures may pose to water ecosystems, thus providing scientific information for environmental quality regulation.

Metagenomic analysis reveals the diversity and distribution of antibiotic resistance genes in thermokarst lakes of the Yellow River Source Area

Thermokarst lakes form as the results of ice-rich permafrost thawing and act as important water resources in cold regions. However, the distributions of antibiotic resistance genes (ARGs) in thermokarst lakes are far less studied. Using metagenomic sequencing approach, we provided the first study to document ARGs in thermokarst lakes of the Yellow River Source Area on the Qinghai-Tibet Plateau (QTP). The results revealed that both sediment and water of the thermokarst lakes harbor diverse ARGs. Multidrug resistance genes were the most diverse, while rifamycin resistance genes were the most abundant with rpoB2 and rpoB genes having the highest proportion. Sediment samples contained more ARGs than water samples, but their composition differed between the two types of samples. However, the composition variations of sediment and water ARGs were closely correlated. The Sorensen dissimilarities of ARGs were controlled by strong turnover processes in sediment samples, and by turnover and nestedness in water samples. High contributions of nestedness were found between sediment and water samples. Moreover, ARGs in water had more significant relationships with environmental variables than that in sediment. Given the role of thermokarst lakes as important water resources in permafrost landscape, as well as intensifying influences of climate change and anthropogenetic activities, thermokarst lakes could bring potential ARG risks, warranting further investigation and evaluation.

Frequency of occurrence and habitat selection shape the spatial variation in the antibiotic resistome in riverine ecosystems in eastern China

BACKGROUND

Riverine ecosystems are one of the most important reservoirs of antibiotic resistance genes (ARGs) in the environment, but the occurrence and controlling factors of ARG distribution in different habitats of riverine ecosystems remain poorly understood. In this study, a metagenomic approach was used to characterize ARG types and their abundance in different habitats (rhizosphere soil, surface bulk soil, bottom bulk soil, and sediment) of riverine ecosystems in eastern China. Sampling sites were located along different rivers of eastern China, which are geographically isolated. Differences in bacterial communities, mobile genetic elements (MGEs), pattern and intensity of human activities, climate, and other environmental factors at the sampling sites and habitats were expected to affect ARG occurrence.

RESULTS

ARGs were observed with high variations in diversity (44-206 subtypes) and abundance (6.85-105.68 ×/Gb). There were significant south-north differences in ARG occurrence in the same habitat, except for surface bulk soil. And the significant difference was found in ARGs among four southern habitats. South-north differences in ARGs of the same habitat were mainly attributed to the combination of different occurrence frequencies and habitat selections of ARGs. Differences in ARG profiles among the four habitats in the south and the north were both mainly attributed to the different occurrence frequencies of ARGs. Bacterial communities and MGEs (Mobile genetic elements) could account for the observed variance in the resistome of riverine ecosystems across eastern China. The co-occurrences of specific ARGs with bacterial communities and MGEs were more frequent at the northern sampling sites than in the south, and co-occurrence patterns (i.e. ARGs and bacterial communities or ARGs and MGEs) varied between the habitats. Moreover, building land in all habitats, except bulk soils, showed significant positive correlations with ARG abundance.

CONCLUSION

This study reveals a high variance in the resistome of riverine ecosystems in eastern China and its controlling factors. We appeal to the importance of assessment of ARGs in the riverine ecosystem and the need for future prevention and intervention of ARG spread.

Metagenomic insights into the antibiotic resistome in freshwater and seawater from an Antarctic ice-free area

The comprehensive profiles of antibiotic resistance genes (ARGs) in the Antarctic water environments and their potential health risks are not well understood. The present study characterized the bacterial community compositions and ARG profiles of freshwater (11 samples) and seawater (28 samples) around the Fildes Region (an ice-free area in Antarctica) using a shotgun metagenomic sequencing approach for the first time. There were significant differences in the compositions of the bacterial community and ARG profiles between freshwater and seawater. In the 39 water samples, 114 ARG subtypes belonging to 15 ARG types were detectable. In freshwater, the dominant ARGs were related to multidrug and rifamycin resistance. In seawater, the dominant ARGs were related to peptide, multidrug, and beta-lactam resistance. Both the bacterial community compositions and ARG profiles were significantly related to certain physicochemical properties (e.g., pH, salinity, NO₃^-). Procrustes analysis revealed a significant correlation between the bacterial community compositions and ARG profiles of freshwater and seawater samples. A total of 31 metagenome-assembled genomes (MAGs) carrying 35 ARG subtypes were obtained and identified. The results will contribute to a better evaluation of the ARG contamination in relation to human health in the Antarctic aquatic environments.

Metagenomic binning and assembled genome analysis revealed the distinct composition of resistome and mobilome in the Ili River

Rivers play an important role in receiving and transporting the resistome among different environmental compartments. However, the difference in resistome and mobilome between the water and sediment and their underlying mechanisms were still poorly understood. In this study, the Ili River, an important water source in the arid area of Central Asia, was selected as the studied target. The comprehensive profile of resistome and mobilome and their host in water and sediment were studied based on metagenomic binning and assembled genome (MAG) analysis. The relative abundance of resistome and mobilome in sediment were 28.0 - 67.8 × /Gb and 46.5 - 121.1 × /Gb, respectively, which were significantly higher than those in water (23.1 - 52.8 ×/Gb and 25.3 - 67.7 ×/Gb). Multidrug and macrolides-lincosamides-streptogramin (MLS) resistance genes were the main ARG types in both water and sediment from relative abundance. Transposases dominated the relative abundance of mobilome, followed by insert elements and integrases. Strong correlations were found between the relative abundance of resistome and mobilome (r > 0.6 and p < 0.01) in both water and sediment, indicating the mobilome played an important role in the propagation of resistome in the Ili River. The main hosts for multidrug resistance genes via MAG analysis differed in water (Alphaproteobacteria and Gammaproteobacteria) and sediment (Gammaproteobacteria). Distinct compositions of resistome and mobilome existed between water and sediment in the Ili River. Specificity-occupancy analysis of the differential resistome and mobilome showed that occurrence frequencies and habitat selections of the differential ARGs shaped the resistome of water and sediment. In contrast, habitat was the main driver that shaped the mobilome in the Ili River.

The formation of specific bacterial communities contributes to the enrichment of antibiotic resistance genes in the soil plastisphere

Soil serves as a major reservoir of both antibiotic resistance genes (ARGs) and microplastics. However, the characteristics of the antibiotic resistome in the soil plastisphere remain largely unknown. In this study, we used metagenomic approaches to reveal the changing patterns of ARGs and the bacterial community and their associations in response to three types of microplastics (light density polyethylene, LDPE; polypropylene, PP; polystyrene, PS) using particles 550 µm or 75 µm in diameter. The total ARG abundances significantly increased in the plastisphere and varied across plastic types. The LDPE plastisphere had the highest ARG total abundance and lowest Shannon diversity index, indicating that this plastic had the most severe negative impact on soil bacterial diversity. The PP plastisphere contained higher relative abundances of the pathogenic bacteria Acinetobacter johnsonii and Escherichia coli, demonstrating the higher pathogenic risk of the microbial communities enriched in the plastisphere. Specifically, multidrug resistance genes (ceoB and MuxB) co-existed with more than four microbial taxa, increasing the potential risk of ARG spread in pathogenic bacteria. These findings implied that the plastisphere acts as a hotspot for acquiring and spreading antibiotic resistance and may have long-term negative effects on the soil ecosystem and human health.

Tidal flat aquaculture pollution governs sedimentary antibiotic resistance gene profiles but not bacterial community based on metagenomic data

Coastal tidal flats are intersection zones between terrestrial and marine environments and are considered repositories of pollutants from anthropogenic activities (e.g., fishery and aquaculture). Specifically, the prevalence of antibiotics and antibiotic resistance genes (ARGs) in coastal aquaculture environments pose critical threats to estuarine ecosystems. However, the contribution of aquaculture to the occurrence and abundance of ARGs and community assemblies has not been fully explored in tidal flat zones. Thus, we investigated ARGs profiles, ARG-carrying host bacteria, and their associate microbial community in the Dongtai and Sheyang tidal flat aquaculture regions of Jiangsu, China using metagenomic assembly methods. The antibiotic concentrations in the sediment samples ranged from nd to 35.50 ng/g dw, and the antibiotic pollution in the Dongtai tidal flat was more severe than in the Sheyang tidal flats. Metagenomic assembly indicated that a total of 247 ARG subtypes associated with ARG 33 types were characterized across all samples and their abundance in the Dongtai region exceeded that in the Sheyang region. Meanwhile, 21 bacteria in the tidal flat aquaculture were identified as ARG-carrying pathogens, including Escherichia coli, Vibrio fluvialis, and Staphylococcus aureus. Using neutral and null modeling analysis to determine the community ecological processes, the results revealed bacterial and ARG communities were generally dominated by stochastic and deterministic processes, respectively. The above results suggested that aquaculture pollution was contributed to shape ARG profiles in tidal flats. The observed deterministic processes affecting the ARG community in tidal flat aquaculture also provides an effective foundation to control the risks of environmental antibiotic resistance through reducing aquaculture antibiotic usage.

Linkage and driving mechanisms of antibiotic resistome in surface and ground water: Their responses to land use and seasonal variation

Antibiotic resistance in the environment, mostly mediated by antibiotic resistance genes (ARGs), has posed a threat to ecological and human health. Contamination of surface water and groundwater with ARGs has become a serious environmental concern. However, the distinctions and similarities across ARG profiles, the various ecological processes associated with ARGs, the driving mechanisms for ARG profiles in surface water and groundwater, and how they respond to land use and seasonal variation remain unknown. To tackle these issues, the contamination of ARGs in surface water and groundwater in central China was investigated using metagenomic technology. The results indicated that seasonal changes in ARG abundance and diversity were inconsistent across surface water and groundwater, and that the relationship between ARGs in surface water and groundwater was stronger during the rainy season. Land use had a greater effect on ARGs in surface water than in groundwater and was stronger during the dry season than during the rainy season. More interestingly, the ideal buffer zones with the greatest impact of land use on the ARGs of surface water and groundwater had distinct radii: 1500 m for both dry and rainy seasons in surface water, and 1000 m for dry season and 500 m for rainy season in groundwater. Furthermore, stochastic mechanisms mediated by mobile gene elements (MGEs) contribute significantly more to ARG assemblages than deterministic processes, particularly in groundwater. Furthermore, our results also showed that ARG enrichment in microbial communities was host- dependent, and the risk of ARGs in groundwater was greater both during the rainy season and dry season. In conclusion, the findings have improved our understanding of the relationship and driving mechanisms of ARGs in surface and ground water, as well as their responses to land use and seasonal variation, which may be beneficial for limiting ARG pollution in a watershed with high levels of anthropogenic activity.

Antibiotic resistomes in water supply reservoirs sediments of central China: main biotic drivers and distribution pattern

Water supply reservoirs form one of the critical drinking water resources. Their water quality directly affects human health. However, reservoir sediments have not received adequate attention in antibiotic resistance genes (ARGs) dissemination, though they reflect long-term ARGs contamination of water supply reservoirs. Moreover, the physicochemical parameters in water supply reservoir sediments are generally better than those in the other media. Thus, the main ARGs biotic drivers of the media would demonstrate their unique characteristics. In this study, sediment samples were collected from 10 water supply reservoirs in central China, and the antibiotic resistomes were determined with the metagenomic method. As revealed from the results, 174 ARGs (18 ARG types) were detected in the reservoir sediment. Besides, multidrug-, sulfonamide-, and vancomycin-ARGs were the dominant ARGs in the sediment samples. The macrolide-resistant Microcystis was prevalent (100% detection frequency with 0.35% average percentage) in reservoir sediments and posed potential risks to human health. Furthermore, the results of the Mantel test and VPA demonstrated that mobile genetic elements (MGEs) were the more essential biotic drivers in ARG contents of reservoir sediments rather than the bacteria community.

Metagenome-Assembled Genomes Reveal Mechanisms of Carbohydrate and Nitrogen Metabolism of Schistosomiasis-Transmitting Vector Biomphalaria Glabrata

Biomphalaria glabrata transmits schistosomiasis mansoni which poses considerable risks to hundreds of thousands of people worldwide, and is widely used as a model organism for studies on the snail-schistosome relationship. Gut microbiota plays important roles in multiple aspects of host including development, metabolism, immunity, and even behavior; however, detailed information on the complete diversity and functional profiles of B. glabrata gut microbiota is still limited. This study is the first to reveal the gut microbiome of B. glabrata based on metagenome-assembled genome (MAG). A total of 28 gut samples spanning diet and age were sequenced and 84 individual microbial genomes with ≥ 70% completeness and ≤ 5% contamination were constructed. Bacteroidota and Proteobacteria were the dominant bacteria in the freshwater snail, unlike terrestrial organisms harboring many species of Firmicutes and Bacteroidota. The microbial consortia in B. glabrata helped in the digestion of complex polysaccharide such as starch, hemicellulose, and chitin for energy supply, and protected the snail from food poisoning and nitrate toxicity. Both microbial community and metabolism of B. glabrata were significantly altered by diet. The polysaccharide-degrading bacterium Chryseobacterium was enriched in the gut of snails fed with high-digestibility protein and high polysaccharide diet (HPHP). Notably, B. glabrata as a mobile repository can escalate biosafety issues regarding transmission of various pathogens such as Acinetobacter nosocomialis and Vibrio parahaemolyticus as well as multiple antibiotic resistance genes in the environment and to other organisms. IMPORTANCE The spread of aquatic gastropod Biomphalaria glabrata, an intermediate host of Schistosoma mansoni, exacerbates the burden of schistosomiasis disease worldwide. This study provides insights into the importance of microbiome for basic biological activities of freshwater snails, and offers a valuable microbial genome resource to fill the gap in the analysis of the snail-microbiota-parasite relationship. The results of this study clarified the reasons for the high adaptability of B. glabrata to diverse environments, and further illustrated the role of B. glabrata in accumulation of antibiotic resistance in the environment and spread of various pathogens. These findings have important implications for further exploration of the control of snail dissemination and schistosomiasis from a microbial perspective.

Cross-regional scale pollution of freshwater biofilms unveiled by antibiotic resistance genes

A comprehensive global profile of the distribution of ARGs in freshwater biofilms is lacking. We utilized metagenomic approaches to reveal the diversity, abundance, transferability and hosts of ARGs in 96 freshwater biofilm samples from 38 sampling sites across four countries. The abundant ARGs were associated with bacitracin, multidrug, polymyxin macrolide-lincosamide-streptogramin (MLS) aminoglycoside, β-lactam, chloramphenicol, sulfonamide and tetracycline resistance, consistent with the spectrum of antibiotics commonly used in human or veterinary medicine. As expected, the resistome in freshwater biofilm habitats was significantly influenced by geographical location and human footprint. Based on the co-occurrence pattern revealed by network analysis, mdtC, kdpE, and emrB were proposed as ARG indicators in freshwater biofilms that can be used to evaluate the abundance of 46 other co-occurring ARG subtypes quantitatively. Metagenomic assembly analysis revealed that the identified ARGs were hosted by more than 46 bacterial phyla, including various pathogens, which greatly expands the knowledge of resistome diversity in freshwater biofilms. Our study points to the central roles of biofilms in harbouring ARGs. The results could enhance understanding the distribution of ARGs in freshwater habitats, thereby strengthening the global environmental risk assessment and management of ARGs.

Metagenomic insights into the profile of antibiotic resistomes in sediments of aquaculture wastewater treatment system

To meet the rapidly growing global demand for aquaculture products, large amounts of antibiotics were used in aquaculture, which might accelerate the evolution of antibiotic-resistant bacteria (ARB) and the propagation of antibiotic genes (ARGs). In our research, we revealed the ARGs profiles, their co-occurrence with mobile genetic elements (MGEs), and potential hosts in sediments of a crab pond wastewater purification system based on metagenomic analysis. The residual antibiotic seems to increase the propagation of ARGs in the crab pond, but there was no clear relationship between a given antibiotic type and the corresponding resistance genes. The effect of aquaculture on sediment was not as profound as that of other anthropogentic activities, but increased the relative abundance of sulfonamide resistance gene. A higher abundance of MGEs, especially plasmid, increased the potential ARGs dissemination risk in crab and purification ponds. Multidrug and sulfonamide resistance genes had greater potential to transfer because they were more frequently carried by MGEs. The horizontal gene transfer was likely to occur among a variety of microorganisms, and various ARGs hosts including Pseudomonas, Acinetobacter, Escherichia, and Klebsiella were identified. Bacterial community influenced the composition of ARG hosts, and Proteobacteria was the predominant hosts. Overall, our study provides novel insights into the environmental risk of ARGs in sediments of aquaculture wastewater treatment system.

Antibiotic resistome in a large urban-lake drinking water source in middle China: Dissemination mechanisms and risk assessment

The increasing pollution of urban drinking water sources by antibiotic resistance genes (ARGs) threatens human health worldwide. However, the distribution and influencing factors of ARGs, especially how to reveal the risks of ARGs in this environment remains unclear. Hence, Chaohu Lake was selected as an example to investigate the characteristics of ARGs and explore the interactions among physicochemical factors, microorganisms, and ARGs by metagenomic approach. In this work, 75 ARG subtypes with an average of 30.4 × /Gb (ranging from 15.2 ×/Gb to 57.9 ×/Gb) were identified, and multidrug and bacA were most frequent in Chaohu Lake. Non-random co-occurrence patterns and potential host bacteria of ARGs were revealed through co-occurrence networks. Microbial community and mobile genetic elements (MGEs) were the major direct factors in ARG profiles. The dissemination of ARGs was mainly driven by plasmids. Considering the interactions among MGEs, human bacterial pathogens, and ARGs, antibiotic resistome risk index (ARRI) was proposed to manifest the risks of ARGs. Overall, our work systemically investigated the composition and associated factors of ARGs and built ARRI to estimate the potential risks of ARGs in a typical urban drinking water source, providing an intuitive indicator for managing similar lakes.

Effect of Nano-Silver on Formation of Marine Snow and the Underlying Microbial Mechanism

Roller experiments were conducted to explore the effect of nano-silver on the formation of marine snow and the underlying microbial mechanism. With the increasing concentration of nano-Ag from 1 ng/L to 1 mg/L, the formation and aggregation of marine snow particles were solidly suppressed in a dose-dependent pattern. Moreover, the formed marine snows tended to be thinner fibrous particles with smaller size and increased edge smoothness and compactness in the presence of nano-Ag. The microbial analyses indicated that nano-Ag not only inhibited the development of biomass but also changed the species composition and functional profile of the microbial community. Nano-Ag obviously inhibited most of the abundant species, except for some myxobacteria, which is unfavorable for the microbial community stability. For the microbial functions, some major biological processes including the growth, metabolic, and cellular processes were also inhibited by the high dosage of nano-Ag. The strong microbial inhibition of nano-Ag would contribute to the suppression on the formation of marine snow. Specifically, the function genes of extracellular polymeric substance synthesis and secretion were significantly reduced by nano-Ag, which might be the key and straight microbial factor in suppressing the formation of marine snow.

Fate, mobility, and pathogenicity of drinking water treatment plant resistomes deciphered by metagenomic assembly and network analyses

Antibiotic resistance genes (ARGs) have been regarded as emerging environmental contaminants. The profile of resistome (collection of all ARGs) in drinking water and its fate during drinking water treatment remain unclear. This study applied metagenomic assembly combined with network analysis to decipher the profile, mobility, host, and pathogenicity of resistomes in two full-scale drinking water treatment plants (DWTPs), each applying conventional treatment and advanced treatment of ozonation followed by biological activated carbon filtration. In source waters and effluents of each treatment process collected from both DWTPs, 215 ARGs belonging to 20 types were detected with total concentration ranging from 6.30 ± 1.83 to 5.20 ± 0.26 × 10⁴ copies/mL. Both the conventional and advanced DWTPs were revealed to effectively reduce the concentration of total ARGs, with the average removal efficiency of 3.61-log₁₀ and 2.21-log₁₀, respectively. Multiple statistical analyses (including network analysis) indicated drinking water resistome correlated tightly with mobile gene elements (MGEs) and bacterial community, with the latter acting as the premier driver of resistome alteration in DWTPs. Further analysis of ARG-carrying contigs (ACCs) assembled from drinking water metagenomes (i) tracked down potential bacterial hosts of ARGs (e.g., Proteobacteria phylum as the major pool of resistome), (ii) provided co-localization information of ARGs and MGEs (e.g., MacB-E7196 plasmid1), and (iii) identified ARG-carrying human pathogens (e.g., Enterococcus faecium and Ralstonia pickettii). This work firstly determined the concentration, mobility incidence, and pathogenicity incidence of DWTP resistomes, based on which the actual health risk regarding antibiotic resistance could be quantitatively assessed in further study, providing a useful direction for decision-making concerning the risk control of ARGs in DWTPs.

Response of microbial nitrogen transformation processes to antibiotic stress in a drinking water reservoir

Effects of antibiotics on microbial nitrogen transformation processes in natural aquatic ecosystems are largely unknown. In this study, we utilized the ¹⁵N stable isotope tracers and metagenomic sequencing to identify how antibiotics drive nitrogen transformation processes in Danjiangkou Reservoir, which is the largest artificial drinking water reservoir in China. We retrieved 51 nitrogen functional genes, and found that the highest abundances of nitrate reduction and denitrification-related genes occurred in dissimilatory nitrogen transformation pathways. ¹⁵N-labelling analysis substantiated that denitrification was the main pathway for nitrogen removal, accounting for 57.1% of nitrogen loss. Nitrogen functional genes and antibiotic resistance genes co-occurred in Danjiangkou Reservoir, and they were mainly carried by the denitrifying bacteria such as Rhodoferax, Polaromonas, Limnohabitans, Pararheinheimera, Desulfobulbus, and Pseudopelobacter. Genome annotation revealed that antibiotic deactivation, Resistance-Nodulation-Division and facilitator superfamily efflux pumps were responsible for the multiple-resistance to antibiotics in these bacteria. Moreover, antibiotics showed non-significant effects on nitrogen transformation processes. It is speculated that denitrifying bacteria harboring ARGs played crucial roles in protecting nitrogen transformation from low-level antibiotics stress in the reservoir. Our results highlight that denitrifying bacteria are important hosts of ARGs, which provides a novel perspective for evaluating the effects of antibiotics on nitrogen cycle in natural aquatic ecosystems.

Vertical and horizontal distributions of clinical antibiotic resistance genes and bacterial communities in Danjiangkou Reservoir, China

The Danjiangkou Reservoir is an important water source for the middle route of the South-to-North Water Diversion Project in China. The current anthropogenic pollution of this reservoir is a great public health concern. Therefore, the horizontal and vertical distributions of seven clinical antibiotic resistance genes (ARGs) and bacterial communities in the Danjiangkou Reservoir were investigated using qPCR and next-generation sequencing, respectively. The average relative abundance of ARGs was 3.01 × 10^-4 to 1.90 × 10^-3 and 4.66 × 10^-4 to 1.85 × 10^-3 in horizontal and vertical profiles, respectively. There was a significant difference in the vertical composition of ARGs, which was caused by different media (i.e. water column vs. sediment). No significant differences in the composition of ARGs were found in the horizontal profile. Proteobacteria was the most abundant phylum, followed by Actinobacteria in the Danjiangkou Reservoir. The beta diversity pattern of the microbial communities in the vertical profile was consistent with that of the ARGs. Moreover, a significant difference in the horizontal composition of the bacterial communities among these water columns was found. All of these factors have resulted in noticeably different co-occurrence patterns of ARGs and bacterial communities between water columns and surface sediment samples. ARGs were closely associated with Proteobacteria and Chloroflexi in the sediment samples, indicating potential anthropogenic pollution in the Danjiangkou Reservoir. Although there was no significant correlation between the occurrence of ARGs and 11 opportunistic pathogens, our results point to potential risks for the development of multi-resistant pathogens due to the simultaneous presence of ARGs and pathogens in the study area. These results provide a good basis for thorough ecological evaluation and remediation of the Danjiangkou Reservoir.

The key environmental influencing factors for the change of sediment bacterial community and antibiotics resistance genes in a long-term polluted lake, China

In order to investigate the key environmental influencing factors for the change of sediment bacterial community structure (BCS) and antibiotics resistance genes (ARGs) in a long-term Quinolone antibiotics (QNs) and heavy metals (HMs) polluted lake, 16S rRNA MiSeq High-throughput Sequencing and Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) analysis methods were applied. Baiyangdian lake was subdivided into three habitats: (1) Habitat 1: greatly influenced by municipal wastewater; (2) Habitat 2: mainly impacted by aquaculture sewage and domestic wastewater; and (3) Habitat 3: with the least human activities. Through One-way ANOVA analysis, the results showed that most of QNs and HMs showed significant difference among three habitats. Both the highest richness and diversity indices of bacterial community appeared in Habitat 3. The abundance of Multidrug, Phenicol, Aminoglycoside, Teracycline, and Quinolone ARGs exhibited the highest values in Habitat 1, while the abundance of Macrolide-Lincosamide-Streptogramin (MLS), Rifamycinm, and Sulfonamide ARGs appeared the highest values in Habitat 2. The result of redundancy analysis exhibited that 68.8% and 93.8% of the change in BCS and ARGs can be explained by environmental factors, respectively. Pb (explained 37.5% (p = 0.002)) and Fleroxacin (FLE) (explained 51.3% (p = 0.026)) were the most important factors for the variation of BCS and ARGs, respectively. Therefore, the effects of antibiotics and HMs on BCS and ARGs should be simultaneously paid more attention. Furthermore, the ARGs results by PICRUSt were similar to the results by metagenomic shotgun sequencing analysis, thus the PICRUSt analysis method can be used in the future research.

Soil type shapes the antibiotic resistome profiles of long-term manured soil

Animal manure fertilization facilitates the proliferation and dissemination of antibiotic resistance genes (ARGs) in soil, posing high risks to humans and ecosystem health. Although studies suggest that soil types could shape the ARG profiles in greenhouse soil, there is still a lack of comparative studies on the fate of ARGs in different types of manured soils under field trials. Thus, a metagenomic approach was used to decipher the fate of ARGs in 12-year long-term fertilized (inorganic fertilizer, compost manure and a mix of them) acidic, near-neutral and alkaline soils. A total of 408 unique ARG subtypes with multidrug, glycopeptide, beta-lactam and aminoglycoside resistance genes were identified as the most universal ARG types in all soil samples. Genes conferred to beta-lactam was the predominant ARG type in all the manure-amended soils. Genomic and statistical analyses showed that manure application caused the enrichment of 98 and 91 ARG subtypes in acidic and near-neutral soils, respectively, and 8 ARG subtypes in alkaline soil. The abundances of Proteobacteria (acidic and near-neutral soils) and Actinobacteria (alkaline soil), which are the potential hosts of ARGs, were clearly increased in manured soils. Random forest modelling and Pearson correlation analysis revealed that the soil properties (pH and bio-available Zn) and mobile genetic elements had considerable impacts on the transmission of ARGs. A structural equation model further indicated that soil types shaped the ARG profiles by significantly (P < 0.01) influencing the soil properties, bacterial abundance and bacterial diversity, where bacterial abundance was the major factor influencing the ARG profiles. This study systematically explored the mechanisms shaping the ARG profiles of long-term manured soils, and this information could support strategies to manage the dissemination of ARGs in different soil types.

Deciphering of antibiotic resistance genes (ARGs) and potential abiotic indicators for the emergence of ARGs in an interconnected lake-river-reservoir system

This study aimed to decipher the patterns of antibiotic resistance genes (ARGs) and linkages of key abiotic indicators with ARGs in an interconnected lake-river-reservoir system. The results showed that seasonal variations in the relative abundance of ARGs and mobile gene elements (MGEs) were significant (KW, p < 0.05). ARGs representative of fecal pollution and natural environment were primarily distributed in the river and reservoir, respectively. The lake, river, and reservoir shared 54.5% of ARGs subtypes, most of which are multidrug resistance genes encoding for efflux pumps. Network results showed that ARGs conferring resistance to aminoglycoside frequently co-occurred with class 1 integrons and Limnohabitans. The resistance risks were low and associated with non-corresponding ARGs, and the highest resistance risk was caused by enrofloxacin in the Dianshan Lake. Fluorescence indices derived from two methods exhibited consistent positive correlations with abundance of individual genes (i.e. aada1 and aadA2-03) as well as total aminoglycoside resistance genes (Pearson, p < 0.05). Moreover, ARGs indicators of human and animal fecal pollution showed linkages with humic-like and fulvic-like indices (Pearson, p < 0.05). The results provide novel insights into the roles of abiotic factors on indicating dynamics of ARGs in aquatic environment impacted by anthropogenic activities.

Twenty-first century molecular methods for analyzing antimicrobial resistance in surface waters to support One Health assessments

Antimicrobial resistance (AMR) in the environment is a growing global health concern, especially the dissemination of AMR into surface waters due to human and agricultural inputs. Within recent years, research has focused on trying to understand the impact of AMR in surface waters on human, agricultural and ecological health (One Health). While surface water quality assessments and surveillance of AMR have historically utilized culture-based methods, culturing bacteria has limitations due to difficulty in isolating environmental bacteria and the need for a priori information about the bacteria for selective isolation. The use of molecular techniques to analyze AMR at the genetic level has helped to overcome the difficulties with culture-based techniques since they do not require advance knowledge of the bacterial population and can analyze uncultivable environmental bacteria. The aim of this review is to provide an overview of common contemporary molecular methods available for analyzing AMR in surface waters, which include high throughput real-time polymerase chain reaction (HT-qPCR), metagenomics, and whole genome sequencing. This review will also feature how these methods may provide information on human and animal health risks. HT-qPCR works at the nanoliter scale, requires only a small amount of DNA, and can analyze numerous gene targets simultaneously, but may lack in analytical sensitivity and the ability to optimize individual assays compared to conventional qPCR. Metagenomics offers more detailed genomic information and taxonomic resolution than PCR by sequencing all the microbial genomes within a sample. Its open format allows for the discovery of new antibiotic resistance genes; however, the quantity of DNA necessary for this technique can be a limiting factor for surface water samples that typically have low numbers of bacteria per sample volume. Whole genome sequencing provides the complete genomic profile of a single environmental isolate and can identify all genetic elements that may confer AMR. However, a main disadvantage of this technique is that it only provides information about one bacterial isolate and is challenging to utilize for community analysis. While these contemporary techniques can quickly provide a vast array of information about AMR in surface waters, one technique does not fully characterize AMR nor its potential risks to human, animal, or ecological health. Rather, a combination of techniques (including both molecular- and culture-based) are necessary to fully understand AMR in surface waters from a One Health perspective.

Reduction of iron (hydr)oxide-bound arsenate: Evidence from high depth resolution sampling of a reducing aquifer in Yinchuan Plain, China

Sediment in fluvial-deltaic plains with high-As groundwater is heterogenous but its characterization of As and Fe oxidation states lacks resolution, and is rarely attempted for aqueous and solid phases simultaneously. Here, we pair high-resolution (> 1 sample/meter) Fe extended fine-structure spectroscopy (EXAFS, n = 40) and As X-ray absorption near-edge spectroscopy (XANES, n = 49) with groundwater composition and metagenomics measurements for two sediment cores and their associated wells (n = 8) from the Yinchuan Plain in northwest China. At shallower depths, nitrate and Mn/Fe reducing sediment zones are fine textured, contain 9.6 ± 5.6 mg kg^-1 of As(V) and 2.3 ± 2.7 mg kg^-1 of As(III) with 9.1 ± 8.1 g kg^-1 of Fe(III) (hydr)oxides, with bacterial genera capable of As and Fe reduction identified. In four deeper 10-m sections, sulfate-reducing sediments are coarser and contain 2.6 ± 1.3 mg kg^-1 of As(V) and 1.1 ± 1.0 mg kg^-1 of As(III) with 3.2 ± 2.6 g kg^-1 of Fe(III) (hydr)oxides, even though groundwater As concentrations can exceed 200 μg/L, mostly as As(III). Super-enrichment of sediment As (42-133 mg kg^-1, n = 7) at shallower depth is due to redox trapping during past groundwater discharge. Active As and Fe reduction is supported by the contrast between the As(III)-dominated groundwater and the As(V)-dominated sediment, and by the decreasing sediment As(V) and Fe(III) (hydr)oxides concentrations with depth.

A metagenomic-based method to study hospital air dust resistome

As a symbol of the defense mechanisms that bacteria have evolved over time, the genes that make bacteria resist antibiotics are overwhelmingly present in the environment. Currently, bacterial antibiotic resistance genes (ARGs) in the air are a serious concern. Previous studies have identified bacterial communities and summarized putative routes of transmissions for some dominant hospital-associated pathogens from hospital indoor samples. However, little is known about the possible indoor air ARG transportation. In this study, we mainly surveyed air-conditioner air dust samples under different airflow conditions and analyzed these samples using a metagenomic-based method. The results show air dust samples exhibited a complex resistome, and the average concentration is 0.00042 copies/16S rRNA gene, which is comparable to some other environments. The hospital air-conditioners can form resistome over time and accumulate pathogens. In addition, our results indicate that the Outpatient hall is one of the main ARG transmission sources, which can distribute ARGs to other departments (explains >80% resistome). We believe that the management should focus on ARG carrier genera such as Staphylococcus, Micrococcus, Streptococcus, and Enterococcus in this hospital and our novel evidence-based network strategy proves that plasmid-mediated ARG transfer can occur frequently. Overall, these results provide insights into the characteristics of air dust resistome and possible route for how ARGs are spread in air.

Occurrence and distribution of antimicrobial resistance genes in the soil of an industrial park in China: A metagenomics survey

As zoned areas of industries, industrial parks have great impacts on the environment. Several studies have demonstrated that chemical compounds and heavy metals released from industrial parks can contaminate soil, water, and air. However, as an emerging pollutant, antimicrobial resistance genes (ARGs) in industrial parks have not yet been investigated. Here, we collected soil samples from 35 sites in an industrial park in China and applied a metagenomics strategy to profile the ARGs and virulence factors (VFs). We further compared the relative abundance of ARGs between the sites (TZ_31-35) located in a beta-lactam antimicrobial-producing factory and other sites (TZ_1-30) in this industrial park. Metagenomic sequencing and assembly generated 14, 383, 065 contigs and 17, 631, 051 open reading frames (ORFs). Taxonomy annotation revealed Proteobacteria and Actinobacteria as the most abundant phylum and class, respectively. The 32 pathogenic bacterial genera listed in the virulence factor database (VFDB) were all identified from the soil metagenomes in this industrial park. In total, 685,354 ARGs (3.89% of the ORFs) and 272,694 virulence factors (VFs) (1.55% of the ORFs) were annotated. These ARGs exhibited resistance to several critically important antimicrobials, such as rifampins, fluroquinolones, and beta-lactams. In addition, no significant difference in the relative abundance of ARGs was observed between sites TZ_31-35 and TZ_1-30, indicating that ARGs have already disseminated widely in this industrial park. The present study gave us a better understanding of the whole picture of the resistome and virulome in the soil of the industrial park and suggested that we should treat the industrial park as a whole in the surveillance and maintenance of ARGs.