- Research
- Open access
- Published:
Metagenomic approach revealed the mobility and co-occurrence of antibiotic resistomes between non-intensive aquaculture environment and human
Microbiome volume 12, Article number: 107 (2024)
Abstract
Background
Aquaculture is an important food source worldwide. The extensive use of antibiotics in intensive large-scale farms has resulted in resistance development. Non-intensive aquaculture is another aquatic feeding model that is conducive to ecological protection and closely related to the natural environment. However, the transmission of resistomes in non-intensive aquaculture has not been well characterized. Moreover, the influence of aquaculture resistomes on human health needs to be further understood. Here, metagenomic approach was employed to identify the mobility of aquaculture resistomes and estimate the potential risks to human health.
Results
The results demonstrated that antibiotic resistance genes (ARGs) were widely present in non-intensive aquaculture systems and the multidrug type was most abundant accounting for 34%. ARGs of non-intensive aquaculture environments were mainly shaped by microbial communities accounting for 51%. Seventy-seven genera and 36 mobile genetic elements (MGEs) were significantly associated with 23 ARG types (p < 0.05) according to network analysis. Six ARGs were defined as core ARGs (top 3% most abundant with occurrence frequency > 80%) which occupied 40% of ARG abundance in fish gut samples. Seventy-one ARG-carrying contigs were identified and 75% of them carried MGEs simultaneously. The qacEdelta1 and sul1 formed a stable combination and were detected simultaneously in aquaculture environments and humans. Additionally, 475 high-quality metagenomic-assembled genomes (MAGs) were recovered and 81 MAGs carried ARGs. The multidrug and bacitracin resistance genes were the most abundant ARG types carried by MAGs. Strikingly, Fusobacterium_A (opportunistic human pathogen) carrying ARGs and MGEs were identified in both the aquaculture system and human guts, which indicated the potential risks of ARG transfer.
Conclusions
The mobility and pathogenicity of aquaculture resistomes were explored by a metagenomic approach. Given the observed co-occurrence of resistomes between the aquaculture environment and human, more stringent regulation of resistomes in non-intensive aquaculture systems may be required.
Video Abstract
Introduction
Aquaculture occupies a significant portion of the global food supply, and the growing consumption of aquatic animal protein surpasses that of all other animal proteins combined. Asia represented 69% of worldwide aquaculture production, with China alone accounting for 35% [1, 2]. The demand for aquaculture animals is expected to increase by 27% by 2030 [3]. At present, in response to the increasing demand for aquatic products, the typical approach is the intensification of aquaculture production or extensive use of antibiotics [4, 5]. However, the extensive use of antibiotics will lead to the emergence and enrichment of antibiotic resistance bacteria and antibiotic resistance genes (ARGs) [6, 7]. Antibiotic resistance in aquaculture might further lead to the prevalence of aquatic animals increasing and aquaculture production reducing [1, 8]. In addition, ARGs could be transferred through the food chain into humans, posing potential health risks [1, 9].
Previous studies have found that high concentrations of antibiotics were the main driving factors of the development and transfer of resistomes in aquaculture environments [10, 11]. Even though the antibiotics use was restricted, antibiotic residues in livestock manure or irrigation water could form a certain selective pressure promoting horizontal gene transfer (HGT) through mobile genetic elements (MGEs) [12,13,14]. However, most previous studies focused on intensive aquaculture environments, where antibiotics were often present at high concentrations. Traditional intensive aquaculture often demands a large number of cultured animals in a confined space, which will promote the transfer of resistome and parasite infection, resulting in a reduction of aquaculture production [15]. Non-intensive aquaculture, which is based on privately owned fishponds with much less antibiotic use than intensive aquaculture, is the most common form of aquaculture in southern China [16]. Non-intensive aquaculture provides more space and avoids overuse of feed. The relatively extensive operation mode of non-intensive aquaculture is more likely to produce environmental interferences [17]. However, the influencing factors of the emergence and transfer of resistomes in non-intensive aquaculture environments are not yet clear. Although some ARGs have been found in non-intensive aquaculture environments [18,19,20], their potential mobility and pathogenicity have not been well evaluated.
Under the “One Health” concept, the transfer of resistomes between the environment, animals, and humans has received extensive attention [8, 21, 22]. Previous researches on the resistome transfer often focused on limited environmental components, such as animal waste and surrounding soil or water [23,24,25], fishpond sediment and water [12, 16, 19, 26], activated sludge, and influent or effluent from wastewater treatment plants [27,28,29]. There is still very little scientific understanding of resistome development in an integrated system including humans, animals, and the environment. It is thus necessary to explore the pattern of ARG transfer between complex environmental system and humans to deepen the understanding of potential health risks.
In this study, we conducted comprehensive sampling (including the chicken gut, fish gut, fishpond sediment, and water) from a non-intensive aquaculture farm and collected local human gut metagenomic data from a public database. A metagenomic approach was employed to characterize the distribution and mobility of resistomes. The results revealed the potential transfer risk of various ARGs between aquaculture environment and human. The pathogens carrying multiple ARGs shared with the aquaculture system and human gut were identified. Our study emphasizes the transfer potential as well as the health risks of aquaculture resistomes and provides a scientific basis for controlling resistomes within the One Health concept.
Materials and methods
Sampling and sequencing
Sampling activities were conducted at a non-intensive aquaculture farm in southern China (113.33 E, 22.10 N). The fish and chickens were raised on the farm without antibiotics. The chickens were fed with natural grains free of antibiotics, the composition of the fish feed could be found in Table S1, and the concentrations of antibiotics and metals in fishpond water were measured (Table S2). Chicken gut samples were collected from chickens at different growth stages with different weights. The fish gut samples of different species, including carp, tilapia, grass carp, bighead carp, and crucian carp, were collected. The fish gut specimens were squeezed out, and the contents were collected carefully to avoid contamination of the intestinal tissue [30]. Fresh chicken droppings were collected from the coop as chicken gut samples [31]. The sediment and water samples were collected at the four corners of the fishponds and then mixed separately. For each water sample, 300 mL of water was filtered using 0.22-µm microporous membranes (Millipore, MA, USA) to enrich microorganisms. All samples were stored at −80 ℃ until DNA extraction. Finally, 35 samples were collected, including 6 chicken guts, 17 fish guts, 8 sediments, and 4 water samples. The detailed sample information can be found in Table S3.
E.Z.N.A.® Water DNA Kit (Omega Bio-Tek, USA) was used for DNA extraction from water samples, and FastDNA® SPIN Kit for Soil (MP Biomedicals, USA) was used for DNA extraction of the other samples according to the manufacturers’ protocols. The extracted DNA was used to construct the library with an insert size of 350 bp and sequenced using the paired-end (PE) 150 bp strategy on the Illumina NovaSeq platform by Novogene company (Beijing, China). All the metagenomic sequencing raw data have been submitted to the Sequence Read Archive (SRA) of NCBI (PRJNA988937). The “aquaculture system samples” in downstream data analysis refers to the collection of all chicken gut, fish gut, sediment, and water samples.
To further analyze the sharing of resistomes between aquaculture system and human, 30 metagenomic data sets about the local human gut (Guangdong Province, China) were downloaded from the NCBI SRA database (Table S4). According to the reference, the volunteers did not receive any antibiotic treatment for 2 months before sampling [32]. Metagenomic data sets with size > 4.5 Gb were selected, and a total of 311.6 Gb of data was downloaded. The raw data was filtered according to the human genome reference (hg18) to eliminate human host DNA before subsequent analysis [32].
Identification and quantification of bacteria, ARGs, and MGEs
All raw reads were quality-filtered using fastp (v0.20.1) [33] to remove sequences with a quality value < 20 or ambiguous nucleotides > 10. After quality filtering, clean reads were used for taxonomic annotation and quantification by Kraken 2 (v2.0.7) [34] and Bracken (v2.0) [35]. ARGs-OAP (v3.0) [36, 37] was used for identification and quantification of ARGs. Clean reads were aligned with the structured ARG database (SARG v3.0) by BLAST + (v2.12.0) (similarity ≥ 80%, e value ≤ 1e − 7 and query coverage ≥ 75%). Finally, the abundance of ARGs was indicated as a “copy of ARG per copy of 16S rRNA gene” based on ARGs-OAP [38, 39]. The analysis of MGEs followed the similar procedure as ARGs, and the SARG database was replaced with the MGE database constructed in the previous reference [40]. This MGE database contains 278 distinct gene annotations and at least 2000 unique sequences.
Metagenomic assembly and binning
After quality filtering, clean reads were assembled using SPAdes (v3.15.3) [41] with the parameters “–meta”. The quality of contigs was evaluated by quast (v5.0.2) [42], and contigs longer than 1000 bp were used for subsequent analysis. Kraken2 (v2.0.7) and Bracken (v2.0) were employed for identifying the taxonomy of contigs. Prodigal [43] (v2.6.3) was used to predict open reading frames (ORFs) on contigs. ARGs and MGEs were searched among the predicted ORFs using BLAST + (v2.12.0) against the database of ARGs and MGEs (refer to 2.2 for database details, e value ≤ 1e − 5, identity ≥ 60%). The ARG-carrying contig (ACC) will be used for the subsequent analysis about the co-occurrence of ARGs and MGEs. The opportunistic pathogenic ACC was identified according to a published list of opportunistic pathogenic species [44].
The co-assembly results of each sample type (chicken gut, fish gut, sediment, water, and human gut) were imported to MetaWRAP (v1.3.2) [45] to recover metagenomic-assembled genomes (MAGs), and the bin_refinement module was performed to improve the quality of MAGs. The completeness and contamination of the recovered MAGs were estimated using CheckM (v1.2.0) [46] with default parameters (completeness > 50% and contamination < 10%). Only high-quality MAGs (completeness—5*contamination ≥ 50%) were retained for downstream analysis. The recovered MAGs were dereplicated using dRep (v3.4.0) [47], and the MAGs with average nucleotide identity (ANI) ≥ 95% were considered to be the same species [48, 49]. The Quant_bin module (default parameters) from MetaWRAP was used to calculate the relative abundance of non-redundant MAGs, and genome copies per million reads (GPMR) was used as the abundance unit. GTDB-Tk (v1.4.1) [50] was used to assign taxonomic classifications to the MAGs, and FastTree [51] was employed to perform phylogenetic analyses and construct phylogenetic trees, which were then visualized by iTOL (https://itol.embl.de).
Co‑occurrence analysis of ARGs and MGEs
To investigate the transfer potential of ARGs between the aquaculture system and human, gene arrangements were deciphered on ACC. The different ARG-MGE combinations that existed simultaneously between different sample groups were visualized by R package gggenes. To further analyze the co-occurrence of ARG-MGE on contigs, two indices introduced by Li et al. [52] were employed to measure their co-occurrence (incidence of encountering and average minimum distance). Briefly, the average minimum distance was calculated by adding the distance from each ARG to the nearest MGE and then dividing it by the total number of ARGs. The number of MGEs found within a certain range per ARG was counted, and the average number of MGEs found in each contig was used to determine the incidence of ARGs encountering MGEs.
The model and method of ARG source tracking
SourceTracker (v1.0) was a useful classification tool based on Bayesian and machine learning [53]. It has been verified to be effective in predicting ARG sources from complex environmental samples [54, 55]. Here, SourceTracker was used for tracking human gut ARGs from various environments, including soil, river, estuary, activated sludge, and non-intensive aquaculture system. The different environment data were collected from the NCBI SRA database (Table S5) and processed in accordance with the above methods (refer to 2.2). The abundances of ARG subtypes were used as SourceTracker input and the following parameters were set: burnin = 100, nrestarts = 10, delay = 10, α = 0.001, β = 0.01, and rarefaction_depth = 1000. The leave-one-out strategy was used to evaluate the performance of the model and the suitability of source settings as previous research described [56].
Statistical analysis and visualization
All statistical analyses were performed in R (v4.1.3) software. Data sets were normalized before analyses, and the results with p < 0.05 by t test were considered as statistically significant. The differences in ARG composition and abundance between different sample groups were revealed by principal coordinate analysis (PCoA). The visualization of network analysis was performed by Gephi (0.9.2) [57]. Variation partitioning analysis (VPA) was used to explore the effects of microbial communities and MGEs on ARGs in aquaculture systems. TBtools [58] was used to make heatmaps. These statistical analyses were performed by vegan package in R, and the ggplot2 package was used for result visualizations.
Results
ARG profile in the aquaculture environment
A total of 22 types and 576 subtypes of ARGs were detected, and the distribution of ARGs was slightly different among different sample groups (Fig. 1a). Overall, multidrug resistance genes (MDR) were relatively abundant in aquaculture system samples, accounting for about 34% of the total ARG abundance. In contrast, macrolide-lincosamide-streptogramin (MLS) and tetracycline resistance genes were more abundant in human gut samples, accounting for about 24% and 41% of the total ARG abundance. MDR and tetracycline resistance genes were relatively abundant ARG types in chicken gut counting for 30% and 35%. The abundances of MDR and bacitracin resistance genes were relatively high in sediments and water, counting for 31% and 27%. Interestingly, we observed that MDR were consistently the most abundant ARGs across the five fish species of different feeding habits (i.e., grass carp, bighead carp, tilapia, carp, and crucian carp).
The average abundance of all ARG types detected in chicken gut, fish gut, sediment, water, and human gut were 0.050 ± 0.03, 0.002 ± 0.001, 0.006 ± 0.002, 0.005 ± 0.0009, and 0.030 ± 0.04 copies of ARG per copy of 16S rRNA gene, respectively (the data were expressed as means ± standard deviation). Compared with the other sample groups, the average abundances of ARGs were higher in the chicken gut and the human gut. The abundance of ARGs in the human gut was significantly different from fish gut, sediment, and water (t test, p < 0.001) (Fig. 1b). In order to compare the composition of ARGs in different sample groups, PCoA based on Bray–Curtis distance was performed and obvious clustering between the sample groups was observed (p < 0.01) (Fig. 1c). The clustering of water and sediment indicated that their composition of ARGs was not significantly different, while the fish gut samples were clustered separately.
ARG composition and abundance profile across samples. a Heatmap of ARG relative abundance (log-transformed). Different color blocks at the top represent different sample groups. b Box plot of ARG abundance in different sample groups. Significant difference between groups is indicated by asterisk with t test (***p < 0.001). c PCoA shows the sample clustering based on ARG abundance profile (95% confidence interval)
The identification of aquaculture core ARGs
In order to better understand the connection of ARGs between aquaculture and human gut, the common and distinct ARGs between various sample groups were identified. Four hundred fifty-seven ARG subtypes were detected in the human gut, and 195 of them were uniquely carried (Fig. 2a). There were 262 ARG subtypes detected simultaneously in the human gut and aquaculture system (Fig. 2b), accounting for 57% of the ARGs carried in the human gut. It is worth noting that 85 ARG subtypes were present in all five sample groups, and 66 ARG subtypes were carried by both human gut and chicken gut (Fig. 2a), indicating certain overlap of ARGs between aquaculture system and human.
The core ARGs in aquaculture system. a UpSet diagram shows the number of ARGs shared and unique between different sample groups. Solid black points indicate ARG occurrence in the sample group, and points linked by lines indicate ARGs shared by different sample groups. The bar chart on the top shows the number of ARG unique or shared by sample groups, and the bar chart on the left shows the total number of detected ARG subtypes. b The Venn diagram shows the number of ARGs shared by the human gut and the aquaculture system. c The percent stacked histogram shows the composition of typical and core ARGs (indicated by different colors). Human gut and environmental typical ARGs defined as ARGs detected in all human gut and environment samples (fishpond sediment and water samples), respectively. Fish gut typical ARGs defined as ARGs detected in 24% fish gut samples, and chicken gut typical ARGs defined as ARGs detected in 50% chicken gut samples
The core ARGs of the aquaculture system were identified according to the ARG occurrence frequency. The negative correlation between the subtype number and occurrence frequency of detected ARGs (R2 = 0.892) demonstrated that the number of persistent ARGs in the aquaculture system was limited (Fig. S1). Notably, only 6 ARG subtypes had occurrence frequency greater than 80% (described as core ARGs), and 229 ARG subtypes had occurrence frequency less than 20% (described as transient ARGs). In addition, the average abundance of core ARGs was 5.95 × 10–3 copies of ARG per copy of 16S rRNA gene, while the average abundance of the transient ARGs was only 7.31 × 10–5 copies of ARG per copy of 16S rRNA gene which was much lower than core ARGs. Overall, the core ARGs with low diversity but high abundance were found in the aquaculture system. These core ARGs accounted for 29% ± 20% and 4% ± 2% of the ARGs composition in aquaculture systems and human gut, respectively (Fig. 2c). It is worth noting that the core ARGs occupied a substantial fraction of ARG abundance in fish gut samples (40% ± 21%).
The influencing factors of ARGs in the aquaculture system
Several results revealed the aquaculture ARGs were shaped by microbial communities and MGEs. Firstly, the abundance of ARGs was significantly correlated with the number (Fig. 3c) and abundance (Fig. 3d) of MGEs. Secondly, network analysis was used to explore the associations between ARGs, MGEs, and potential hosts (Fig. 3a). Seventy-seven genera and 36 MGEs were significantly associated with 23 ARG types (p < 0.05). The detail information of network analysis results can be found in Table S6. The results indicated that Anabaena, Aphanizomenon, Cylindrospermopsis, and Raphidiopsis were significantly correlated with 13 ARGs. Fusobacterium was observed to be associated with 8 ARGs. The insertion sequence (IS) was the dominant MGE type associated with ARGs, accounting for approximately 75% of all MGE types. Tn916 and qacEdelta were associated with 16 and 13 ARGs, respectively. Aminoglycoside, sulfonamide, and fosmidomycin resistance genes were associated with the largest number of genera and MGEs. These ARG potential hosts and the correlation between ARGs and MGEs were also observed in analysis based on contigs and MAGs (refer to results 3.4 and 3.5). Finally, the impacts of MGEs and bacteria on ARGs were further quantitatively characterized by VPA (Fig. 3b). The VPA results indicated that bacterial flora could explain 51% of ARGs profile and MGEs could explain the other 19%.
Influencing factors of ARGs in aquaculture system. a The network analysis shows the relationship between ARGs, MGEs, and bacteria (genus level) (p < 0.05). Larger nodes indicate more connections. The red lines indicate positive correlations, and the green lines indicate negative correlations. b VPA shows the interpretation rate of ARGs by MGEs and bacteria. c Linear relationship between ARG abundance and the number of MGEs (p < 0.001). d Linear relationship between ARG abundance and MGE abundance (p < 0.001) (discrete sample points f3, f4, and f15 were excluded)
In addition, Spearman’s rank coefficient revealed the correlation of core ARGs with MGEs and genera (Table S7). There was a relatively strong correlation between Thiobacillus with MLS resistance gene macB and multidrug transporter (p ≤ 0.01). Many types of IS (ISSfl3, IS679, and IS621) were significantly positively correlated with macB, multidrug transporter, and cAMP-regulatory protein (p < 0.01). The positive correlation was also identified between transposon tniA and macB. The above evidence suggested that aquaculture core ARGs might be influenced by MGEs and bacteria.
Co-occurrence of diverse ARG and MGE combinations
The transfer risk of ARGs in aquaculture systems was assessed by the incidence of co-occurrence and average minimum distance between ARGs and MGEs [52]. As the distance from ARGs increased, the incidence of MGE co-occurrence increased (Fig. 4a), and the distance between ARGs and MGEs in aquaculture system was significantly shorter than that in the human gut (p < 0.001) (Fig. 4b). In order to further reveal the co-occurrence pattern of ARGs and MGEs, we analyzed genetic locations of the two elements on metagenomic-assembled contigs retrieved from all samples. Overall, 71 ACCs were identified, among which 53 ACCs (75%) carried MGEs, 21 ACCs (30%) carried multiple ARGs, and 10 ACCs (14%) carried multiple ARGs and MGEs simultaneously. Furthermore, the co-occurrence of ARGs and MGEs shared by aquaculture systems and human gut samples were investigated, and various co-occurrence patterns of ARGs and MGEs were observed. The aminoglycoside resistance gene aadA, sulfonamide resistance gene sul1, and the quaternary ammonium compound (qac) resistance gene qacEdelta1 form a stable combination and were present in different sample groups (Fig. 4c). The qac gene was reported as a potential ancestor of class I integrons, and clinical class I integrons may contain qacEdelta gene [59, 60], hence the occurrence with the elements will increase the HGT probability of aadA and sul1. It is noteworthy that we identified two contigs W3_NODE_35346 and hg21_NODE_1841, predicted to be derived from opportunistic pathogens Pseudomonas aeruginosa and Escherichia coli, respectively.
The co-occurrence pattern of ARGs and MGEs based on contigs. a The co-occurrence incidence of ARGs and MGEs. The x-axis means the distance of ARG-MGE, and the y-axis means the co-occurrence incidence. With the increase of distance, the co-occurrence rate also increases and gradually reaches a peak. b The shortest distance of ARG-MGE in the human gut and aquaculture system (**p < 0.01 indicates significant difference). c Gene arrangement patterns of sul1-qacEdelta1 and neighborhood genes in different sample groups. Different colored contigs belong to different sample groups, and the taxonomy names colored in red indicate opportunistic pathogenic bacteria
In addition to the combination of aadA-sul1-qacEdelta1, multiple ARGs and MGEs were found to co-occur on various contigs from different samples, such as sul4-tnpA, tetA-tnpA, and tetM-Tn916 (Fig. S2). Especially, the tetracycline resistance gene tetP, one of the core ARGs, was found to be widely distributed in chicken, fish, and human gut samples (Fig. S2c). The co-occurrence of another core ARG transcriptional regulatory protein CpxR between integrase int2 and transposase tnpA was also observed in chicken and human gut samples.
The source tracking of ARGs in the human gut
To further track the source of human gut ARGs, the Bayesian-based tool SourceTracker was used, and different environment types were set as potential sources. Four hundred fifty-seven metagenomic data collected from definite environments (including non-intensive aquaculture environments from this study) were used for analysis by a leave-one-out strategy to compare predicted environmental types with definite environmental types. The results showed that all environment types were predicted with more than 80% accuracy by five independent tests (Fig. S3). The aquaculture environment type was predicted with 88% accuracy comparable to previous studies [54]. The results indicated that SourceTracker could effectively predict ARG sources from different environment types. Then, the model was used to predict the potential source of ARGs in the human gut by five independent analyses. The result showed that the aquaculture environment was the main potential source of human gut ARGs, accounting for 81% (Table S8). This result identified the potential risk of ARG transfer from the aquaculture system to the human gut.
Identification and quantification of ARG-carrying genomes (ACGs)
Recovering MAGs and detecting ACGs could provide insights into the host of ARGs and assess potential transfer risks. In total, 475 high-quality MAGs were recovered and the taxonomy of 460 non-redundant MAGs was annotated (Table S9). At the phylum level, the most abundant bacteria were Firmicutes, Bacteroidota, and Proteobacteria, and 12 bins were identified as opportunistic pathogenic bacteria (Fig. S4). The Circos plot displayed the distribution of ARGs across MAGs (Fig. 5), with MDR and bacitracin resistance genes being the two most abundant ARGs carried by MAGs. It is noteworthy that the MDR carried by MAGs contained 29 subtypes, while the bacitracin resistance genes contained only one subtype, bacA.
The heatmap further displayed the relative average abundances of ACGs in different sample groups and the types of ARGs carried (Fig. 6). Strikingly, some ACGs carried multiple ARGs which may be highly mobile. For all 81 ACGs, 23 carried multiple ARGs, 18 carried MGEs, and 6 ACGs carried multiple ARGs and MGEs simultaneously (Table S10), which were considered as genomes with high ARG enrichment and mobility, thus posing a high risk of promoting ARG dissemination. These 6 ACGs were hg_bin.40 (Fusobacterium_A sp900543175), hg_bin.256 (Eubacterium_G sp000434315), hg_bin.49 (Holdemanella porci), fg_bin.37 (Acetatifactor intestinalis), fg_bin.33 (Solirubrobacterales bacterium 67–14), and cg_bin.33 (Nannocystaceae) carrying 48, 46, 38, 7, 6, and 2 ARGs, respectively, which encode bacitracin, beta-lactam, fosfomycin, fosmidomycin, kasugamycin, MLS, MDR, polymyxin, and tetracycline resistance.
The heatmap shows the ACG abundance, mobility, and carried ARG types. a The color bar A indicates the ACG carries one or more ARGs. The color bar B indicates the ACG does not carry MGE, carries one MGE, or carries multiple MGEs. The color bar C indicates the sample group from which the ACG was recovered. The heatmap shows the average abundance of ACGs in different sample groups (data was log transformed). ACG name labels in red indicate opportunistic pathogenic bacteria. b The number of ARG types carried by ACGs. The color indicated the carried ARG number
Especially, 11 MAGs carrying core ARGs were identified (Table S11). Five MAGs carried multiple core ARGs and 7 MAGs also carried MGEs, while no MAG carrying core ARGs were identified as opportunistic pathogens. It is worth noting that hg_bin.40 carried MGEs and 4 core ARGs simultaneously, and the relative abundances of hg_bin.40 was high in both chicken gut and human gut samples (Fig. 6a), suggesting the potential risk of core-ARG transfer.
Discussion
A non-intensive aquaculture farm was selected for sampling to fill the research gap on non-intensive aquaculture resistomes. A study about the large-scale intensive aquaculture farms in Guangdong Province found that the concentrations of erythromycin in fishpond water could reach 80–1400 ng/L [16], much higher than this study (< 1 ng/L, Table S2). Despite the lack of antibiotic pressure in non-intensive aquaculture environment, multiple ARG types were also detected (Fig. 1a). It has been reported that MDR and bacitracin resistance genes were the main resistance types in freshwater aquaculture [20], which were similar to our findings (Fig. 1a). It is worth noting that although the average abundance of ARGs detected in this research was approximately 0.06 copies of ARG per copy of 16S rRNA gene (Fig. 1b), which was lower than 0.25 copies of ARG per copy of 16S rRNA gene in a previous study about integrated freshwater aquaculture farms [26], the potential transfer risk of ARGs could be still demonstrated (Fig. 4c). Despite almost zero use of antibiotics in this relatively closed non-intensive aquaculture system, fishpond irrigation water sourced from surrounding river may introduce extra antibiotic pressure. The relatively low concentrations of antibiotics and metals (Table S2) might form a combined selection pressure leading to the occurrence and transfer of ARGs [61]. These findings suggest that the mobility of ARGs needs more attention even in non-intensive aquaculture environments with low antibiotic pressure.
Previous studies have suggested that the microbial structure of fish gut was influenced by the host species and diet [62, 63]. However, the microbial structure of 5 different fish species in our research did not show significant differences (R2 = 0.33, p = 0.18) (Fig. S5a). The fish gut microbiome could profile the intestinal ARGs by affecting the hosts’ metabolism [64], while there was no significant difference between the intestinal ARGs of different fish species either (R2 = 0.23, p = 0.55) (Fig. S5b). This might be due to the fact that the water and feed received by the fish were from the same surrounding environments in our research, thus shaping the similar gut community and resistomes of different fish species. These results suggest that the fish gut microbiome and ARGs might be determined by the habitat rather than fish species [65,66,67]. Extensive large-scale research is needed to explore the contribution of habitat environments to fish gut resistomes in the future. Moreover, approximately 40% of the ARGs carried by the fish were identified as aquaculture core ARGs that indicated fish may play an important role in the transmission of aquaculture ARGs (Fig. 2c).
The core ARGs have been demonstrated to exist in a variety of environments [68,69,70]. To accurately identify core ARGs in aquaculture systems, more stringent criteria (approximately the top 3% most abundant with detection frequency greater than 80% in all aquaculture system samples) were used to define aquaculture core ARGs compared to previous studies [69, 71]. The results showed that 33% of core ARGs in aquaculture environments were MDR, similar to the composition of core ARGs in municipal wastewater treatment plants [72, 73] and farmland soils [71]. Notably, only a small fraction of the human gut ARG profile consisted of the core ARGs (4% ± 2%). Previous research indicated that the HGT of core ARGs was seldom observed due to the barriers of bacterial phylogeny and the absence of MGEs [74, 75]. However, the core ARGs were highly shared between different aquaculture system samples (29% ± 20%, Fig. 2c), and one of the aquaculture core ARG tetP was shared between the aquaculture system and human gut (Fig. S2). In this work, different samples were collected from a relatively closed aquaculture system, circulation of materials, and microbiome transfer could be frequent, which might increase the possibility of core ARGs exchange between different environmental components.
Those ARGs capable of transferring from the environment to human gut are of high health risks, which call for priority attention when conducting relevant surveillance and risk assessment [76]. In this study, 262 ARG subtypes were co-shared by aquaculture system and human gut (Fig. 2b). ARGs in aquaculture system had a positive correlation with MGEs indicating the potential mobility of aquaculture ARGs (Fig. 3c and d). The gene arrangement showed the possible transfer of different ARG-MGE combinations between aquaculture environment and human. In addition, some MGEs (intI1, tniB, and tnpA) were located in proximity of ARGs, suggesting the high mobility of these ARGs (Fig. 4c). The relatively short distance between ARGs and MGEs implied high risk of ARGs transfer, of which the mobility of those clinical important ARGs between aquaculture environment and human may raise public health concern (Fig. 4b). Furthermore, the results of source tracking showed that aquaculture might occupied 80% of the ARG contribution in the human gut (Table S8), which further indicated the potential risk of ARG transfer from aquaculture environment to the human gut. In fact, it has been reported that the resistomes could be shared between human and animal gut [9, 77,78,79]; however, there is limited evidence of resistomes transmission from environment to human due to factors like ecological barrier, phylogenetic restriction, and low mobility [74, 75]. The combination of sul1-qacEdelta1, which could transfer from aquaculture system to human gut, has not only been detected in this study (Fig. 4c), but also been reported in wastewater treatment plants [80], plastisphere [81], and drinking water [82]. These results provide the evidence for resistomes transfer potential between environments and humans, thus indicate that the high mobility of aquaculture resistomes, especially the ARG transfer from environmental bacteria to clinical pathogens raise concern of infectious untreatable [76, 83]. Furthermore, samples based on time series will help explore the specific transfer process of resistomes between the environment and human in the future.
Identification and quantification of human opportunistic pathogenic bacteria could assist in assessing relevant health risks [84]. In this study, contig hg21_NODE_1841 obtained from human gut was identified as E. coli and W3_NODE_35346 obtained from fishpond water was identified as P. aeruginosa (Fig. 4c). These contigs carried ARG-MGE and were recognized as human opportunistic pathogenic bacteria meanwhile indicating the potential pathogenicity of aquaculture resistomes. Additionally, 12 opportunistic pathogenic MAGs were obtained and 75% of them were recovered from human gut (Table S9). It is worth noting that high abundances of hg_bin.40 (Fusobacterium) were detected in both chicken gut and human gut (Fig. 6a). Fusobacterium species have been found to cause a wide range of opportunistic infections related to oral diseases and colorectal carcinoma [85]. The Fusobacterium MAG was found carried multiple ARGs and MGEs, which was similar to a previous study that Fusobacterium might carry tetracycline resistance genes [86]. The detection of same opportunistic pathogenic bacteria both in human gut and aquaculture system samples revealed the potential role for aquaculture in promoting pathogenicity of human gut microbiome. In fact, only 3% MAGs were detected as opportunistic pathogenic bacteria and 1% MAGs carried multiple ARGs and MGEs in this study. Compared to culture-based experiments, MAG-based methods are limited in only capturing those abundant bacteria in sequencing, providing indirect evidence about ARG mobility or bacterial pathogenicity [87], lacking phenotypic characteristics [77]. In the future, transcriptome data and culture-based experiments are needed to complement the metagenomics in providing a more comprehensive insight into the risks of ARGs and the host microorganisms.
Conclusion
This research explores the transmission capability of resistomes and potential risks to human health in the understudied non-intensive aquaculture environment. Multiple types of ARGs were identified with possible mobility in the aquaculture environment. The detected ARG-MGE combinations had potential transfer risk between human gut and aquaculture system. The potential transfer of opportunistic pathogens carrying multiple ARGs between aquaculture system and human gut could increase the risk of human disease. In summary, our study emphasizes the aquaculture environment even with almost zero antibiotic application as “hidden” reservoir of ARGs and indicates potential mobility of resistomes across human–environment boundaries.
Availability of data and materials
All sequencing data is available at the Sequence Read Archive (SRA) of NCBI (PRJNA988937).
References
Schar D, Zhao C, Wang Y, Larsson DGJ, Gilbert M, Van Boeckel TP. Twenty-year trends in antimicrobial resistance from aquaculture and fisheries in Asia. Nat Commun. 2021;12(1):5384.
FAO. The State of World fisheries and aquaculture 2020. Sustainability in action. Rome: FAO; 2020. p. 244. https://doi.org/10.4060/ca9229en.
Reverter M, Sarter S, Caruso D, Avarre J-C, Combe M, Pepey E, Pouyaud L, Vega-Heredía S, de Verdal H, Gozlan RE. Aquaculture at the crossroads of global warming and antimicrobial resistance. Nat Commun. 2020;11(1):1870.
Van Boeckel TP, Brower C, Gilbert M, Grenfell BT, Levin SA, Robinson TP, Teillant A, Laxminarayan R. Global trends in antimicrobial use in food animals. Proc Natl Acad Sci. 2015;112(18):5649–54.
Herrero M, Wirsenius S, Henderson B, Rigolot C, Thornton P, Havlík P, de Boer I, Gerber PJ. Livestock and the environment: what have we learned in the past decade? Annu Rev Environ Resour. 2015;40(1):177–202.
Preena PG, Swaminathan TR, Rejish Kumar VJ, Bright Singh IS. Unravelling the menace: detection of antimicrobial resistance in aquaculture. Lett Appl Microbiol. 2020;71(1):26–38.
Schar D, Klein EY, Laxminarayan R, Gilbert M, Van Boeckel TP. Global trends in antimicrobial use in aquaculture. Sci Rep. 2020;10(1):21878.
Zhao Y, Yang QE, Zhou X, Wang F-H, Muurinen J, Virta MP, Brandt KK, Zhu Y-G. Antibiotic resistome in the livestock and aquaculture industries: status and solutions. Crit Rev Environ Sci Technol. 2021;51(19):2159–96.
Cao H, Bougouffa S, Park TJ, Lau A, Tong MK, Chow KH, Ho PL. Sharing of antimicrobial resistance genes between humans and food animals. mSystems. 2022;7(6):e00775–00722.
He X, Xu Y, Chen J, Ling J, Li Y, Huang L, Zhou X, Zheng L, Xie G. Evolution of corresponding resistance genes in the water of fish tanks with multiple stresses of antibiotics and heavy metals. Water Res. 2017;124:39–48.
Shao Y, Wang Y, Yuan Y, Xie Y. A systematic review on antibiotics misuse in livestock and aquaculture and regulation implications in China. Sci Total Environ. 2021;798:149205.
Xiong W, Sun Y, Zhang T, Ding X, Li Y, Wang M, Zeng Z. Antibiotics, antibiotic resistance genes, and bacterial community composition in fresh water aquaculture environment in China. Microb Ecol. 2015;70(2):425–32.
Davies J, Davies D. Origins and evolution of antibiotic resistance. Microbiol Mol Biol Rev. 2010;74(3):417–33.
Forster SC, Liu J, Kumar N, Gulliver EL, Gould JA, Escobar-Zepeda A, Mkandawire T, Pike LJ, Shao Y, Stares MD, et al. Strain-level characterization of broad host range mobile genetic elements transferring antibiotic resistance from the human microbiome. Nat Commun. 2022;13(1):1445.
Lulijwa R, Rupia EJ, Alfaro AC. Antibiotic use in aquaculture, policies and regulation, health and environmental risks: a review of the top 15 major producers. Rev Aquac. 2020;12(2):640–63.
Chen H, Liu S, Xu X-R, Diao Z-H, Sun K-F, Hao Q-W, Liu S-S, Ying G-G. Tissue distribution, bioaccumulation characteristics and health risk of antibiotics in cultured fish from a typical aquaculture area. J Hazard Mater. 2018;343:140–8.
Naylor RL, Hardy RW, Buschmann AH, Bush SR, Cao L, Klinger DH, Little DC, Lubchenco J, Shumway SE, Troell M. A 20-year retrospective review of global aquaculture. Nature. 2021;591(7851):551–63.
Fang H, Huang K, Yu J, Ding C, Wang Z, Zhao C, Yuan H, Wang Z, Wang S, Hu J, et al. Metagenomic analysis of bacterial communities and antibiotic resistance genes in the Eriocheir sinensis freshwater aquaculture environment. Chemosphere. 2019;224:202–11.
Shen X, Jin G, Zhao Y, Shao X. Prevalence and distribution analysis of antibiotic resistance genes in a large-scale aquaculture environment. Sci Total Environ. 2020;711:134626.
Wang C, Liu X, Yang Y, Wang Z. Antibiotic and antibiotic resistance genes in freshwater aquaculture ponds in China: a meta-analysis and assessment. J Clean Prod. 2021;329:129719.
Jin L, Pruden A, Boehm AB, Alvarez PJJ, Raskin L, Kohn T, Li X. Integrating environmental dimensions of “One Health” to combat antimicrobial resistance: essential research needs. Environ Sci Technol. 2022;56(21):14871–4.
Ikhimiukor OO, Odih EE, Donado-Godoy P, Okeke IN. A bottom-up view of antimicrobial resistance transmission in developing countries. Nat Microbiol. 2022;7(6):757–65.
Yu X, Zhou ZC, Shuai XY, Lin ZJ, Liu Z, Zhou JY, Lin YH, Zeng GS, Ge ZY, Chen H. Microplastics exacerbate co-occurrence and horizontal transfer of antibiotic resistance genes. J Hazard Mater. 2023;451:131130.
Gao F-Z, He L-Y, Chen X, Chen J-L, Yi X, He L-X, Huang X-Y, Chen Z-Y, Bai H, Zhang M, et al. Swine farm groundwater is a hidden hotspot for antibiotic-resistant pathogenic Acinetobacter. ISME Commun. 2023;3(1):34.
Luo Y, Tan L, Zhang H, Bi W, Zhao L, Wang X, Lu X, Xu X, Sun R, Alvarez PJJ. Characteristics of wild bird resistomes and dissemination of antibiotic resistance genes in interconnected bird-habitat systems revealed by similarity of blaTEM polymorphic sequences. Environ Sci Technol. 2022;56(21):15084–95.
Xu C, Lv Z, Shen Y, Liu D, Fu Y, Zhou L, Liu W, Chen K, Ye H, Xia X, et al. Metagenomic insights into differences in environmental resistome profiles between integrated and monoculture aquaculture farms in China. Environ Int. 2020;144:106005.
Wang Y, Han Y, Li L, Liu J, Yan X. Distribution, sources, and potential risks of antibiotic resistance genes in wastewater treatment plant: a review. Environ Pollut. 2022;310:119870.
Yin X, Yang Y, Deng Y, Huang Y, Li L, Chan LYL, Zhang T. An assessment of resistome and mobilome in wastewater treatment plants through temporal and spatial metagenomic analysis. Water Res. 2022;209:117885.
Yu K, Li P, He Y, Zhang B, Chen Y, Yang J. Unveiling dynamics of size-dependent antibiotic resistome associated with microbial communities in full-scale wastewater treatment plants. Water Res. 2020;187:116450.
Sáenz JS, Marques TV, Barone RSC, Cyrino JEP, Kublik S, Nesme J, Schloter M, Rath S, Vestergaard G. Oral administration of antibiotics increased the potential mobility of bacterial resistance genes in the gut of the fish Piaractus mesopotamicus. Microbiome. 2019;7(1):24.
Mazhar SH, Li X, Rashid A, Su J, Xu J, Brejnrod AD, Su J-Q, Wu Y, Zhu Y-G, Zhou SG, et al. Co-selection of antibiotic resistance genes, and mobile genetic elements in the presence of heavy metals in poultry farm environments. Sci Total Environ. 2021;755:142702.
Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, Liang S, Zhang W, Guan Y, Shen D, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490(7418):55–60.
Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018;34(17):i884–90.
Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20(1):257.
Jennifer L, Florian PB, Peter T, Steven LS. Bracken: estimating species abundance in metagenomics data. PeerJ Comput Sci. 3:e104. https://doi.org/10.7717/peerj-cs.104.
Yin X, Zheng X, Li L, Zhang A-N, Jiang XT, Zhang T. ARGs-OAP v3.0: Antibiotic-resistance gene database curation and analysis pipeline optimization. Engineering. 2022;27:234–41.
Yin X, Chen X, Jiang XT, Yang Y, Li B, Shum MHH, Lam TTY, Leung GM, Rose J, Sanchez-Cid C, et al. Toward a universal unit for quantification of antibiotic resistance genes in environmental samples. Environ Sci Technol. 2023;57(26):9713–21.
Yang Y, Jiang X, Chai B, Ma L, Li B, Zhang A, Cole JR, Tiedje JM, Zhang T. ARGs-OAP: online analysis pipeline for antibiotic resistance genes detection from metagenomic data using an integrated structured ARG-database. Bioinformatics. 2016;32(15):2346–51.
Li B, Yang Y, Ma L, Ju F, Guo F, Tiedje JM, Zhang T. Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes. ISME J. 2015;9(11):2490–502.
Pärnänen K, Karkman A, Hultman J, Lyra C, Bengtsson-Palme J, Larsson DGJ, Rautava S, Isolauri E, Salminen S, Kumar H, et al. Maternal gut and breast milk microbiota affect infant gut antibiotic resistome and mobile genetic elements. Nat Commun. 2018;9(1):3891.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11(1):119.
Yi X, Liang JL, Su JQ, Jia P, Lu Jl, Zheng J, Wang Z, Feng SW, Luo ZH, Ai HX, et al. Globally distributed mining-impacted environments are underexplored hotspots of multidrug resistance genes. ISME J. 2022;16(9):2099–113.
Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome. 2018;6(1):158.
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25(7):1043–55.
Olm MR, Brown CT, Brooks B, Banfield JF. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J. 2017;11(12):2864–8.
Olm Matthew R, Crits-Christoph A, Diamond S, Lavy A. Matheus Carnevali Paula B, Banfield Jillian F: Consistent metagenome-derived metrics verify and delineate bacterial species boundaries. mSystems. 2020;5(1):e00731–00719.
Parks DH, Chuvochina M, Chaumeil P-A, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat Biotechnol. 2020;38(9):1079–86.
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2020;36(6):1925–7.
Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5(3):e9490.
Li L-G, Xia Y, Zhang T. Co-occurrence of antibiotic and metal resistance genes revealed in complete genome collection. ISME J. 2017;11(3):651–62.
Knights D, Kuczynski J, Charlson ES, Zaneveld J, Mozer MC, Collman RG, Bushman FD, Knight R, Kelley ST. Bayesian community-wide culture-independent microbial source tracking. Nat Methods. 2011;8(9):761–3.
Li L-G, Yin X, Zhang T. Tracking antibiotic resistance gene pollution from different sources using machine-learning classification. Microbiome. 2018;6(1):93.
Brown CM, Mathai PP, Loesekann T, Staley C, Sadowsky MJ. Influence of library composition on sourcetracker predictions for community-based microbial SourceTracking. Environ Sci Technol. 2019;53(1):60–8.
Wang C, Yang H, Liu H, Zhang X-X, Ma L. Anthropogenic contributions to antibiotic resistance gene pollution in household drinking water revealed by machine-learning-based source-tracking. Water Res. 2023;246:120682.
Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. In: Third international AAAI conference on weblogs and social media. 2009. p. 2009.
Chen C, Wu Y, Xia R. A painless way to customize Circos plot: from data preparation to visualization using TBtools. iMeta. 2022;1(3):e35.
Gillings Michael R. Integrons: past, present, and future. Microbiol Mol Biol Rev. 2014;78(2):257–77.
Gillings MR, Xuejun D, Hardwick SA, Holley MP, Stokes HW. Gene cassettes encoding resistance to quaternary ammonium compounds: a role in the origin of clinical class 1 integrons? ISME J. 2009;3(2):209–15.
Li Z-H, Yuan L, Wang L, Liu Q-H, Sheng G-P. Coexistence of silver ion and tetracycline at environmentally relevant concentrations greatly enhanced antibiotic resistance gene development in activated sludge bioreactor. J Hazard Mater. 2022;423:127088.
Parata L, Mazumder D, Sammut J, Egan S. Diet type influences the gut microbiome and nutrient assimilation of Genetically Improved Farmed Tilapia (Oreochromis niloticus). PLoS ONE. 2020;15(8):e0237775.
Sylvain F-É, Holland A, Bouslama S, Audet-Gilbert É, Lavoie C, Val Adalberto L, Derome N. Fish skin and gut microbiomes show contrasting signatures of host species and habitat. Appl Environ Microbiol. 2020;86(16):e00789–e720.
Yukgehnaish K, Kumar P, Sivachandran P, Marimuthu K, Arshad A, Paray BA, Arockiaraj J. Gut microbiota metagenomics in aquaculture: factors influencing gut microbiome and its physiological role in fish. Rev Aquac. 2020;12(3):1903–27.
Kim PS, Shin N-R, Lee J-B, Kim M-S, Whon TW, Hyun D-W, Yun J-H, Jung M-J, Kim JY, Bae J-W. Host habitat is the major determinant of the gut microbiome of fish. Microbiome. 2021;9(1):166.
Yuan L, Wang L, Li Z-H, Zhang M-Q, Shao W, Sheng G-P. Antibiotic resistance and microbiota in the gut of Chinese four major freshwater carp from retail markets. Environ Pollut. 2019;255:113327.
Chen C-Z, Li P, Liu L, Li Z-H. Exploring the interactions between the gut microbiome and the shifting surrounding aquatic environment in fisheries and aquaculture: a review. Environ Res. 2022;214:114202.
Yin X, Deng Y, Ma L, Wang Y, Chan LYL, Zhang T. Exploration of the antibiotic resistome in a wastewater treatment plant by a nine-year longitudinal metagenomic study. Environ Int. 2019;133:105270.
Su J-Q, An X-L, Li B, Chen Q-L, Gillings MR, Chen H, Zhang T, Zhu Y-G. Metagenomics of urban sewage identifies an extensively shared antibiotic resistome in China. Microbiome. 2017;5(1):84.
Yin X, Li L, Chen X, Liu Y-Y. Lam TT-Y, Topp E, Zhang T: Global environmental resistome: distinction and connectivity across diverse habitats benchmarked by metagenomic analyses. Water Res. 2023;235:119875.
Du S, Shen J-P, Hu H-W, Wang J-T, Han L-L, Sheng R, Wei W-X, Fang Y-T, Zhu Y-G, Zhang L-M, et al. Large-scale patterns of soil antibiotic resistome in Chinese croplands. Sci Total Environ. 2020;712:136418.
Li W, Mao F, Ng C, Jong MC, Goh SG, Charles FR, Ng OT, Marimuthu K, He Y. Gin KY-H: Population-based variations of a core resistome revealed by urban sewage metagenome surveillance. Environ Int. 2022;163:107185.
Tian L, Li Q, Cai X, Wang Y, Wang Y, Mao Y. Dynamic distribution and potential transmission of antibiotic resistance genes in activated sludge. Appl Microbiol Biotechnol. 2022;106(19):6785–97.
Munck C, Albertsen M, Telke A, Ellabaan M, Nielsen PH, Sommer MOA. Limited dissemination of the wastewater treatment plant core resistome. Nat Commun. 2015;6(1):8452.
Forsberg KJ, Patel S, Gibson MK, Lauber CL, Knight R, Fierer N, Dantas G. Bacterial phylogeny structures soil resistomes across habitats. Nature. 2014;509(7502):612–6.
Forsberg KJ, Reyes A, Wang B, Selleck EM, Sommer MOA, Dantas G. The shared antibiotic resistome of soil bacteria and human pathogens. Science. 2012;337(6098):1107–11.
Wang Y, Hu Y, Liu F, Cao J, Lv N, Zhu B, Zhang G, Gao GF. Integrated metagenomic and metatranscriptomic profiling reveals differentially expressed resistomes in human, chicken, and pig gut microbiomes. Environ Int. 2020;138:105649.
Wang C, Song Y, Tang N, Zhang G, Leclercq Sébastien O, Feng J. The shared resistome of human and pig microbiota is mobilized by distinct genetic elements. Appl Environ Microbiol. 2021;87(5):e01910–01920.
Ma L, Xia Y, Li B, Yang Y, Li L-G, Tiedje JM, Zhang T. Metagenomic assembly reveals hosts of antibiotic resistance genes and the shared resistome in pig, chicken, and human feces. Environ Sci Technol. 2016;50(1):420–7.
Bengtsson-Palme J, Hammarén R, Pal C, Östman M, Björlenius B, Flach C-F, Fick J, Kristiansson E, Tysklind M, Larsson DGJ. Elucidating selection processes for antibiotic resistance in sewage treatment plants using metagenomics. Sci Total Environ. 2016;572:697–712.
Zhao Y, Gao J, Wang Z, Cui Y, Zhang Y, Dai H, Li D. Distinct bacterial communities and resistance genes enriched by triclocarban-contaminated polyethylene microplastics in antibiotics and heavy metals polluted sewage environment. Sci Total Environ. 2022;839:156330.
Ke Y, Sun W, Jing Z, Zhao Z, Xie S. Seasonal variations of microbial community and antibiotic resistome in a suburb drinking water distribution system in a northern Chinese city. J Environ Sci. 2023;127:714–25.
Martinez JL, Coque TM, Baquero F. What is a resistance gene? Ranking risk in resistomes. Nat Rev Microbiol. 2015;13(2):116–23.
Liang J, Mao G, Yin X, Ma L, Liu L, Bai Y, Zhang T, Qu J. Identification and quantification of bacterial genomes carrying antibiotic resistance genes and virulence factor genes for aquatic microbiological risk assessment. Water Res. 2020;168:115160.
Brennan CA, Garrett WS. Fusobacterium nucleatum—symbiont, opportunist and oncobacterium. Nat Rev Microbiol. 2019;17(3):156–66.
Juricova H, Matiasovicova J, Kubasova T, Cejkova D, Rychlik I. The distribution of antibiotic resistance genes in chicken gut microbiota commensals. Sci Rep. 2021;11(1):3290.
Boolchandani M, D’Souza AW, Dantas G. Sequencing-based methods and resources to study antimicrobial resistance. Nat Rev Genet. 2019;20(6):356–70.
Acknowledgements
Thank The University of Hong Kong Shenzhen Institute of Research and Innovation for the financial support.
Funding
This work was supported by the Natural Science Foundation of Guangdong Province (2022A1515011961), National Natural Science Foundation of China (22193062, 52200070), Natural Science Foundation of Guangdong Province (2023A1515012019).
Author information
Authors and Affiliations
Contributions
T.L., L.L. and M.Y. wrote the main manuscript text. Z.T. and M.Y. designed the work. F.G. collected samples. T.L., F.G. and L.G. prepared figures and analyzed data. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The data sets of human gut used in this study are public. No ethical approval is needed.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Fig. S1.
The relationship between the mean number and abundance of detected ARGs with the occurrence frequency. Fig. S2. Gene arrangements of ARG-MGE combinations. Fig. S3. The result of source tracking by leave-one-out strategy. Fig. S4. The phylogenetic tree of MAGs. Fig. S5. The PCoA of fish gut microbiome and ARGs. ARGs: antibiotic resistance genes. MGEs: mobile genetic elements. MAGs: metagenomic-assembled genomes. PCoA: principal coordinate analysis.
Additional file 2: Table S1.
The composition of the fish feed. Table S2. The concentrations of antibiotics and metals in fishpond water. Table S3. The detailed sample information. Table S4. The information of human gut samples. Table S5. The information of different environment samples used for source tracking. Table S6. The specific network association between ARGs, MGEs and genera. Table S7. The Spearman’s rank correlation between core ARGs, MGEs and genera. Table S8. The source tracking results of human gut ARGs. Table S9. The information of 460 non-redundant MAGs. Table S10. The detailed information of 6 ACGs carried multiple ARGs and MGEs. Table S11. The detailed information of 11 MAGs carrying core ARGs.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Tian, L., Fang, G., Li, G. et al. Metagenomic approach revealed the mobility and co-occurrence of antibiotic resistomes between non-intensive aquaculture environment and human. Microbiome 12, 107 (2024). https://doi.org/10.1186/s40168-024-01824-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40168-024-01824-x