REVIEW article

Front. Microbiol., 02 December 2022

Sec. Evolutionary and Genomic Microbiology

Volume 13 - 2022 | https://doi.org/10.3389/fmicb.2022.1066995

Metagenomic-based surveillance systems for antibiotic resistance in non-clinical settings

  • 1. Delft Bioinformatics Lab, Delft University of Technology, Delft, Netherlands

  • 2. Department of Biotechnology, Delft University of Technology, Delft, Netherlands

  • 3. Infectious Disease and Microbiome Program, Broad Institute of MIT and Harvard, Cambridge, MA, United States

Article metrics

View details

59

Citations

12,9k

Views

7,5k

Downloads

Abstract

The success of antibiotics as a therapeutic agent has led to their ineffectiveness. The continuous use and misuse in clinical and non-clinical areas have led to the emergence and spread of antibiotic-resistant bacteria and its genetic determinants. This is a multi-dimensional problem that has now become a global health crisis. Antibiotic resistance research has primarily focused on the clinical healthcare sectors while overlooking the non-clinical sectors. The increasing antibiotic usage in the environment – including animals, plants, soil, and water – are drivers of antibiotic resistance and function as a transmission route for antibiotic resistant pathogens and is a source for resistance genes. These natural compartments are interconnected with each other and humans, allowing the spread of antibiotic resistance via horizontal gene transfer between commensal and pathogenic bacteria. Identifying and understanding genetic exchange within and between natural compartments can provide insight into the transmission, dissemination, and emergence mechanisms. The development of high-throughput DNA sequencing technologies has made antibiotic resistance research more accessible and feasible. In particular, the combination of metagenomics and powerful bioinformatic tools and platforms have facilitated the identification of microbial communities and has allowed access to genomic data by bypassing the need for isolating and culturing microorganisms. This review aimed to reflect on the different sequencing techniques, metagenomic approaches, and bioinformatics tools and pipelines with their respective advantages and limitations for antibiotic resistance research. These approaches can provide insight into resistance mechanisms, the microbial population, emerging pathogens, resistance genes, and their dissemination. This information can influence policies, develop preventative measures and alleviate the burden caused by antibiotic resistance.

Introduction

Antibiotic-resistant bacterial infections have contributed to 4.95 million deaths worldwide, with 1.27 million of those directly resulting from antimicrobial resistance (AMR; Gigante et al., 2022; Papavarnavas and Mendelson, 2022). It is the leading cause of death worldwide and affects high, low, and middle-income countries (Pokharel et al., 2019; Gupta and Gupta, 2021). Antibiotic-resistant bacterial infections are rising exponentially, making them harder to treat (Dušan et al., 2016; Robinson et al., 2016). This continuation will lead to an estimated 10 million deaths annually by 2050 if the burden of AMR is not alleviated (Van Boeckel et al., 2015).

For most of human history, bacteria have caused infectious diseases and contributed to high mortality rates (MacLean and San, 2019). To counteract this, antibiotics were routinely prescribed to treat and prevent such diseases. This saved millions of lives each year. Due to its success, antibiotics were used in animal husbandry and agricultural practices with the added benefits of growth promotion (Kang et al., 2018). Unfortunately, this success was limited to the emergence of antibiotic-resistant bacteria (ARB) and its effects on human health and longevity (Fair and Tor, 2014; MacLean and San, 2019; Kumar et al., 2020).

The continuous and unwarranted use of antibiotics in humans, animals, and agricultural areas exerts selection pressure on the bacteria found in these sectors. This gives rise to bacteria resistant to multiple antibiotics classes leading to untreatable infections (Davies and Davies, 2010; Cabello and Godfrey, 2016).

The ability of bacteria to resist antibiotics is embedded in their evolutionary history, which has led to various phenotypic and genotypic resistance mechanisms (Pornsukarom and Thakur, 2017; Power, 2019). These include the inactivation, modification, degradation, and expulsion of the antibiotic or its target site, therefore protecting the bacteria (Capita and Alonso-Calleja, 2013; Blair et al., 2015; Munita and Arias, 2016; Kumar, 2017; Kang et al., 2018). It is also possible that bacteria gain resistance through mutations, changing the nature of proteins expressed in the bacterial organism, or by horizontal gene transfer phenomena, which is the exchange of genetic material by the use of mobile genetic elements (MGEs) between bacterial strains (Pornsukarom and Thakur, 2017; MacLean and San, 2019). This allows bacteria to collect numerous resistant traits and ultimately become multi-drug resistant (MDR; Munita and Arias, 2016; von Wintersdorff et al., 2016; Kumar, 2017; Power, 2019).

The ineffectiveness of antibiotics threatens the ability to treat common bacterial infections and minor injuries. Since first-line antibiotics are deemed sub-par, more expensive treatment, alternative therapy and additional care is necessary. This can result in prolonged infections and an increase in morbidity and mortality rates. This puts an economic burden on families, society, and modern medicine (Gandra et al., 2014; Friedman et al., 2016; Aslam et al., 2018; Power, 2019).

For over a decade, antibiotic resistance has been studied using the traditional culture-dependent approach. This involves the use of artificial conditions to culture microorganisms. Unfortunately, some microorganisms cannot be cultured by this approach (Schmieder and Edwards, 2012; Nowrotek et al., 2019). Approximately 80% of bacterial species in the human gut and 99% of environmental bacterial species remain uncultured (Lagier et al., 2012; Nowrotek et al., 2019). Factors such as slow growth, microbial competition, specific growth requirements, and environmental stressors can affect the culturing process. This leaves limited information regarding gene transfer between these bacterial communities (Garmendia et al., 2012; Schmieder and Edwards, 2012; Suenaga, 2012; Nowrotek et al., 2019; Moralez et al., 2021).

Whole metagenome DNA sequencing can describe the genomes of the total microbial community found in nature. This approach enables the study of culturable and non-culturable bacteria by bypassing the need for isolation and laboratory cultivation of microorganisms. DNA directly isolated from the environmental sample can broaden the understanding of the structure, gene/species richness and distribution, and the functional and metabolic potential of a microbial community (Suenaga, 2012).

Applying this approach to antibiotic resistance in different microbial communities can identify known and novel resistance genes and mobile genetic elements, i.e., plasmids, integrons, transposons, and phage’s (Partridge et al., 2018; Stalder et al., 2019; Hall et al., 2022; Yao et al., 2022). This information is the stepping stone to making new policies for infection and prevention control measures, thereby reducing the incidence of infection and optimizing the use of antibiotics by health professionals, in the healthcare, animal, and agricultural industries. In addition, it will provide improved awareness and understanding of antibiotic resistance as a whole.

Unfortunately, AMR monitoring systems focus on clinical and the public health sectors as a representation of AMR as a whole. Organizations such as WHO and the implementation of the Global Antimicrobial Surveillance System (GLASS) have relied on the use of traditional culturable methods for clinical AMR research. This creates a biased outlook on AMR as non-clinical sectors, and non-culturable bacteria are excluded. This review will focus on AMR in non-clinical sectors and the use of bioinformatics as an aid in this research.

DNA sequencing technologies used in antibiotic resistance research

Microbiomes are complex and require sequencing technologies to be comprehensive enough to capture the present representative sequences from species. Several sequencing platforms can be used to study AMR in clinical and non-clinical sectors. They are useful in detecting antibiotic resistance genes (ARGs), virulence factors, MGEs, and diversity of the microbial community (Forbes et al., 2017). The advancements in metagenomics have been driven by next-generation sequencing technologies (NGS) i.e., second and third-generation sequencing technologies (Table 1).

Table 1

2nd Generation3rd Generation
DefinitionNGS technology where the nucleic acid is fragmented and amplified followed by sequencing of many different clusters of short reads taking place in parallel where base detection is monitored by different platforms such as IlluminaNGS technology allows for long-read single-molecular sequencing by monitoring the incorporation of fluorescently labelled nucleotides or by monitoring an electric single as nucleic acid is fed through a nano-sized pore
Advantages+ High reproducibility Requires less DNA/RNA (ng) Strain typing Resistance determinants such as single nucleotide polymorphism level analysis Detection of chromosomal mutations leading to AMR Detection of ARGs and allelic variants+ Suitable for pathogen identification High sensitivity for microorganism detection Detection of ARGs Longer reads allow for the detection and mapping of mobile genetic elements in multi-drug-resistant strains Provides genetic context
LimitationsResults in days therefore cannot be used as a rapid diagnostic tool Complex analysis due to high volumes of data and short reads Unable to link ARGs to genetic context Short reads result in fragmented assemblies around repetitive regions associated with ARGs or mobile genetic elementsHigh cost Resistance determinants such as single nucleotide polymorphisms cannot be analysed Chromosomal mutations leading to AMR cannot be detected Allelic variants cannot be detected A high concentration of input nucleic acid is required
Platforms
  • Life Technologies (Ion Torrent and Ion Proton) Illumina (HiSeq, MiSeq, GenomeAnalyser)

  • Pacific Biosciences (PacBio RS) Oxford Nanopore (MinION, PromethION, GridION)

Second-generation versus third-generation sequencing technologies for antibiotic susceptibility testing.

Second-generation sequencing has been the most widely used sequencing technology for microbial genomics. Second-generation sequencing technologies, i.e., Illumina and Ion Torrent, produce short reads from 100 up to 300 base pairs with high accuracy (~0.1% error rate; De Maio et al., 2019; Yee and Simner, 2019). These short reads provide information on species population, evolutionary relationships, allelic variations, SNPs, and the detection of specific genes, i.e., ARGs (Boolchandani et al., 2019; De Maio et al., 2019; Yee and Simner, 2019). Although short reads provide information on the genomic content of the bacterial isolates in a sample, it does not provide information on the genomic structure, i.e., which genes are present on chromosomes and MGEs such as plasmids which are crucial in understanding the dissemination of ARGs (George et al., 2017; De Maio et al., 2019; Yee and Simner, 2019; Maguire et al., 2020).

Third-generation technologies overcome this limitation by providing a better view of the genetic structure as longer reads (on average 1–100 Kbp) span most of the repetitive sequences, linking ARGs to MGEs, e.g., plasmids (Boolchandani et al., 2019; Yee and Simner, 2019). This long-read technology is advantageous for studying AMR in bacterial isolates, and metagenomic samples as the complexity of assembly is reduced (De Maio et al., 2019). Unfortunately, the high error rate makes it difficult to identify specific allelic variants or SNPs in chromosomal genes leading to AMR (Yee and Simner, 2019). On the other hand, Pacbio and ONT are easily accessible and able to produce raw data in real-time, making them a good tool for the rapid diagnostics (Che et al., 2019). Long reads are also a practical method to determine the species/gene richness, distribution, and functional potential of a microbial community (Garmendia et al., 2012; Crofts et al., 2017).

While third-generation technologies are successful in AMR research, it is limited by the high error rate and low accuracy (Yee and Simner, 2019). To overcome this, a combination of short-reads (i.e., Illumina) with long-reads (i.e., ONT) is a promising way to generate fully resolved and accurate bacterial genome assemblies which characterize antibiotic resistance genes on plasmids and genomes from environmental samples (Liu et al., 2020; Chen et al., 2021; Moralez et al., 2021). This combination allows for better assembly of complex genomes as long reads provide information on the structure of the genome, and short reads can be used to correct errors in long reads (De Maio et al., 2019).

These sequencing technologies directly characterize the microbiomes in humans, animals, and environmental samples, providing high-throughput sequencing of either whole genomes, targeted amplicons, or whole metagenomes (Di Bella et al., 2013).

Whole-genome sequencing (WGS) uses DNA sequencing technologies to sequence the entire genome of isolated pure organisms to characterize genomic variants (Table 2; Collineau et al., 2019). WGS overcomes the limitations of traditional culture-dependent approaches and phenotypic tests, i.e., disk diffusion for AMR (Balloux et al., 2018). By opting to use WGS, more information on the bacterial genetic determinants conveying resistance and their association with mobile genetic elements can be established (Balloux et al., 2018; Oniciuc et al., 2018). This contributes to outbreak detection, infection control, and epidemiological surveillance, all of which are an essential part of the AMR surveillance (Collineau et al., 2019).

Table 2

Whole-genome sequencingMeta-proteomicsWhole-metagenome sequencingQuantitative PCR
DefinitionThe sequencing of an entire genome of a single organism is obtained from a culture to identify species/ strains, genes, and mutations associated with AMR.Identification and quantification of proteins conferring AMR and Multi-drug resistance (MDR) to localized MGEs from microbial communities.A culture-independent approach is used for the identification of all ARGs/ mutations in all organisms present in complex microbial communities.Method for the detection, quantification, and typing of specific microbial species/ strains or ARGs/ markers.
TechniqueDNA is extracted from bacterial culture and sequenced to generate FastQ reads.The generation of peptides and high-resolution mass spectra analysis to identify and quantify proteins.DNA is extracted directly from a sample and sequenced to generate reads representing the entire microbial community.An enzymatic reaction in combination with fluorophores (SYBR Green and TaqMan) is used to amplify a specific target of interest and quantify gene expression levels.
Advantages+ Reproducible High resolution Standardization of analysis methods Rapid turn-around time High specificity Read lengths of 1 – 50KbAccurate prediction Large-scale protein identification Rapid turn-around time+ High resolution Accurate predictionFast High-throughput detection and quantification of target DNA sequences avoid cross-contamination Sensitivity Specificity
Limitations- Requires cultured organism Computationally demanding Cost- False positives from large datasets False discovery rates in large datasets Limited reproducibility- Difficult to differentiate between host and pathogen Expensive Computationally demanding- Incapable of distinguishing between live and dead cells
Outcome
  • Microbial typing and tracing Ability to predict AMR from only genomic data Discovery of novel resistance genes or mutations

Host-pathogen interaction Virulence factors Antibiotic resistance mechanisms Identification of functions involved in biological pathways
  • Species-or strain level identification AMR Virulence potential Detect Mobile genetic elements Discovery of viruses

Determine the presence of specific genes and alleles AMR profiling Toxin production Typing of strains and isolates

Metagenomic sequencing vs. other techniques used to study qualitatively and quantitatively antibiotic resistance.

Amplicon sequencing uses marker genes such as 16S (prokaryotes) or 18S rRNA (eukaryotes) genes. Both approaches require the DNA extracted directly from the microbial community and are subjected to either direct sequencing or amplification via polymerase chain reaction (PCR; Table 2). In both instances, the sequences are queried against a reference database (Table 3). This is a reliable process to identify the composition of a microbial community down genus-level within a sample (Ranjan et al., 2016; Forbes et al., 2017; Oniciuc et al., 2018).

Table 3

General databasesDescriptionUpdateWebsiteReference
ResFinder and ResFinderFGA curated database that identifies acquired genes and/ or finds chromosomal mutations.Apr-21https://cge.cbs.dtu.dk//services/ResFinder/(Zankari et al., 2012; Xie et al., 2018)
CARDManually curated database and an ontology resource that provides information on ARGs and resistance mechanisms.Jul-21https://card.mcmaster.ca/(Ahmed et al., 2007; Jia et al., 2016; Chen et al., 2017; Guo et al., 2017)
SARG (v2)ARG database derived from ARDB, CARD, and NCBI-NR databases for characterization and quantification of ARGs.Apr-21http://smile.hku.hk/SARGs(Yin et al., 2018)
FARMECurated database focused on environmentally derived metagenomes conferring resistance.Mar-19http://staff.washington.edu/jwallace/farme/(Wallace et al., 2017)
MEGARESDatabase containing sequence data for approximately 8,000 manually curated ARGs.Oct-19https://megares.meglab.org/(Lakin et al., 2017)
MUSTARDDatabase of 6,059 AMR determinants from 20 families in the human gut microbiota and curated sets of genes. Identified via functional metagenomics.Sep-17http://mgps.eu/Mustard/(Ruppé et al., 2019)
RESFAMSA profile HMM-based curated database to confirm antibiotic resistance function and organized by ontology.Jan-15http://www.dantaslab.org/resfams/(Gibson et al., 2015)
NDAROCurated and collated data from multiple databases. Contains 5,804 sequences.https://www.ncbi.nlm.nih.gov/bioproject/PRJNA313047
Specific DatabasesDescriptionUpdateWebsiteReference
LacEDCurated database providing information on mutation, sequences, and structures of TEM and SHV β-lactamases.TEMLacED: 2017 SHVED:2010http://www.laced.uni-stuttgart.de/(Thai et al., 2009)
CBMARDatabase identifying and characterizing novel β-lactamases based on Amber classification.Sep-14http://proteininformatics.org/mkumar/lactamasedb/(Srivastava et al., 2014)
BLDBA manually curated database of AR enzymes classified by class, family, and subfamilyJul-21http://bldb.eu/(Naas et al., 2017)

General and specific databases in AMR research.

Metagenomic sequencing involves the fragmentation, sequencing, assembly, and annotation of the total genomic DNA isolated in a given sample (Oniciuc et al., 2018; Xu et al., 2021). Similar to WGS, whole metagenome sequencing (WMS) can provide information on the entire gene content of either prokaryotic and eukaryotic organisms as well as species or strain-level identification and virulence or resistance potential (Ranjan et al., 2016; Quainoo et al., 2017; Hendriksen et al., 2019). In addition, WMS makes it possible to predict the metabolic potential of the microbial community. This comprehensive view of all gene content allows access to information on AMR that is essential for setting up control and prevention strategies in all sectors (Xu et al., 2021).

Meta-proteomics focuses on functional change in a microbiome. Assisting in the control strategies and plays a key role in identifying molecular mechanisms of bacterial pathogenesis and disease outcome determinants. It can aid in developing pathogen-specific treatment strategies that can lower the spread of AMR. Meta-proteomics is reliable for reviewing bacteria in soil, sludge, food, and the ocean (Wilmes and Bond, 2004; Bastida et al., 2014; Williams and Cavicchioli, 2014; Maier et al., 2017). Unfortunately, there are inconsistent protocols for sample preparation, inefficient bioinformatic tools, and challenges in measuring low-abundance proteins within a complex protein sample (Zhang et al., 2018). This technique, while reliable, should be applied with other sequencing technologies to have an in-depth analysis of microbial communities, antibiotic resistance genes, and host-pathogen interaction (Zhang and Figeys, 2019; Khodadadi et al., 2020).

Bioinformatic analysis of metagenomic sequences for antibiotic resistance research

The constant development of new computational pipelines is advantageous in providing an accurate depiction of the resistome from metagenomic sequencing data. These analyses rely on a variety of algorithms: quality control of DNA sequencing data, genome and metagenome assembly algorithms, read-mapping, variant detection, phylogenetics, taxonomic databases, and visualization (McArthur and Tsang, 2017). AMR databases are discussed in more detail in the next section.

The typical process to identify antibiotic resistance genes from sequenced metagenomes is based on one of two general techniques: (i) read-mapping and (ii) assembly. Each of these could be followed by binning or an annotation step (Table 4). A bioinformatics pipeline for analysing metagenomic data for antibiotic resistance research has been presented in Figure 1.

Table 4

Read-mapping approachAssembly-based approachBinning approachAnnotation
DefinitionThe detection of genes without genome assembly. The direct alignment of reads to a reference database.The identification of genes with de novo assembly profiling.The evaluation of taxonomic groups by clustering assembled sequences into individual groups that represent microbial species.The prediction of genes from either assembled or unassembled reads.
FunctionAntibiotic resistance gene discovery and taxonomic identification.Antibiotic-resistant gene discovery and taxonomic identification.Assessment of taxonomic diversity and gene association to taxonomic groups.Gene function and taxonomic identification.
TechniqueHigh-quality reads are directly aligned to a reference database using pairwise alignment tools or splitting reads into k-mers and mapping to reference databases.High-quality reads are assembled into contigs that are aligned to a reference database using BLAST (Altschul et al., 1990) or HMMscan.The process of clustering contigs into individual genome bins by machine learning methods or comparing metagenomic sequences against a reference database of genomic sequences.Sequences are classified by coding and non-coding regions to determine gene function.
Advantages+ Fast Less computationally demanding Identify genes from low abundance organisms in complex communities Resistome analysis of large datasets+ Identification of protein-coding regions Identification of known and novel resistance genes Regulatory and mobile genetic element sequences identified+ Easy Reliable+ High accuracy High sensitivity
Limitations- Lacks information required to study upstream and downstream factors of identified resistance genes Large datasets may cause false-positive predictions- Loss of information due to incomplete and poor assembly Computationally expensive Time-consuming Requires a high genome coverage- Lack of whole-genome sequences within the available databasesComputationally demanding Limited success for short reads Expensive
ToolsSRST2 (Inouye et al., 2014), ARIBA (Hunt et al., 2017), BOWTIE2 (Langmead and Salzberg, 2012), Burrows-wheeler aligner (BWA) (Li and Durbin, 2009)MEGAHIT (Li et al., 2015), MetaSpades (Nurk et al., 2017), MetaVelvet (Namiki et al., 2012), Spades (Bankevich et al., 2012), Velvet, Abyss (Simpson et al., 2009) and SOAPdenovo (Luo et al., 2012)MetaBat2 (Kang et al., 2019), Maxbin2 (Wu et al., 2016), CONCOCT (Alneberg et al., 2014), GroopM (Imelfort et al., 2014), DASTool (Sieber et al., 2018). MetaCluster (Leung et al., 2011), MEGAN (Huson et al., 2007), MGRAST (Glass et al., 2010)RAST (Aziz et al., 2008), IMG (Markowitz et al., 2009), MetaGeneAnnotator (Noguchi et al., 2008), MetaGeneMark (McHardy et al., 2007)

Various bioinformatics approaches for the prediction of AMR and MGEs.

Figure 1

Figure 1

Pathway map of antibiotic agents and antibiotic resistance dissemination between livestock, agriculture, aquaculture, and water. The movement of antibiotic resistance or antibiotic agents is indicated by solid lines. The potential routes of transmission and exposure between non-clinical sectors and antibiotic usage is indicated by the abbreviations: Ro: run-off /leaching, T: treatment, P: prevention, GP: growth promoter, H: animal house, C: cleaning, N: nutrients, I: irrigation (Created with BioRender.com). Routes between non-clinical sectors (↔) and antibiotic usage in non-clinical sectors (→) have been described by the previous studies mentioned. Water ↔ crop: (Hu et al., 2010; Gatica and Cytryn, 2013; Meena et al., 2015; Bürgmann et al., 2018; Rodríguez-Molina et al., 2019), water ↔ soil: (Hu et al., 2010; Seiler and Berendonk, 2012; Michael et al., 2013; Meena et al., 2015; Tang, 2016), water ↔ livestock: (Miyata et al., 2011; Hao et al., 2014; Van Boeckel et al., 2015; Magouras et al., 2017; Madec and Haenni, 2018; McEwen and Collignon, 2018; Opatowski et al., 2021), water ↔ aquaculture: (Seiler and Berendonk, 2012; Shah et al., 2014; Grenni et al., 2018; Pérez-Sánchez et al., 2018; Yuan et al., 2019; Preena et al., 2020a,b; Reverter et al., 2020; Sonone et al., 2020), antibiotic use → aquaculture: (Shah et al., 2014; Pérez-Sánchez et al., 2018; Stärk et al., 2019; Preena et al., 2020b), antibiotic use → livestock: (Ahmed et al., 2007; Marshall and Levy, 2011; Landers et al., 2012; Hao et al., 2014; Speksnijder et al., 2015; Van Boeckel et al., 2017; McEwen and Collignon, 2018), antibiotic use → crops: (Hirniesen et al., 2012; Pan and Chu, 2017; Sundin and Wang, 2018; Zhu et al., 2019; Gao, 2021), livestock ↔ soil and crops: (Popowska et al., 2012; Wang et al., 2015; Lin et al., 2016; Chen et al., 2017; Xie et al., 2018; Pérez-Valera et al., 2019; Zhang et al., 2019).

Read mapping allows for raw read sequences to be directly aligned to reference databases (Table 3) using pairwise tools such as MiniMap2 (Li, 2018), Bowtie2 (Langmead and Salzberg, 2012), Burrows-wheeler aligner (BWA; Langmead and Salzberg, 2012), or by splitting reads into k-mers and mapping to a reference database using k-mer alignment (KMA; Clausen et al., 2016; Figure 2). This provides information on the similarity of one sequence to another by counting the number of alignments to quantify the abundance of similar sequence patterns (Clausen et al., 2016).

In the last decade, read-mapping approaches were deemed superior to assembly approaches for AMR gene detection, i.e., reads are first assembled and mapped to a database using BLAST (Johnson et al., 2008; Hunt et al., 2017; Clausen et al., 2018). The success of an assembly-based approach is highly dependent on the quality of the assembly (Inouye et al., 2014; Clausen et al., 2016). An approach as such is problematic as genes of interest, i.e., ARG, s can be split over two or more contigs and will not be identified if the assembly is poor quality (Clausen et al., 2016, 2018).

Read-mapping based methods overcome this by using the tools above, which allow for fast mapping and alignment of raw reads against large reference genomes and entire databases (Clausen et al., 2016). This can identify genes from low abundance organisms in a complex community with speed and ease of computation. Since epidemiological databases are constantly being updated with new sequences due to natural evolution, the read-mapping approach becomes difficult as results constantly change. There is no guarantee that the read will cover a unique part of the reference sequence and will result in a tie for the best match due to the random selection (Clausen et al., 2016). Fortunately, tools like SRST2 (Inouye et al., 2014) resolve the ties by pre-and post-processing of sequences for read-mapping, which can predict the presence of ARGs in a sample (Clausen et al., 2018; Taş et al., 2021).

Identification and characterization of ARGs can also be achieved using assembly-based approaches (Table 4; Figure 2; Breitwieser et al., 2019). This method allows reads to be assembled into contigs and then queried against reference databases (Table 3). The assembly of WMS data is complicated as there is an uneven or unknown abundance of different genomes (Ghurye et al., 2016).

Unrelated genomes may also contain nearly identical DNA repeats, which could represent MGEs (Boolchandani et al., 2019). In addition, multiple individual organisms could be from the same species but may harbour small genetic differences indicating strain variants (Ghurye et al., 2016). Technical factors, e.g., library preparation, sequencing depth, and sequencing platforms, affect the accuracy in assembling WMS data into larger contigs (Escobar-Zepeda et al., 2015). Metagenomic assemblers such as MEGAHIT (Li et al., 2015), MetaSPades (Nurk et al., 2017), IDBA-UD (Peng et al., 2012), and MetaVelvet (Namiki et al., 2012) have been developed to overcome these challenges by normalizing sequencing depth, correcting read errors, and detecting and reporting genomic variants and repeats (Ghurye et al., 2016; Vosloo et al., 2021). To date, no single assembler is the best at accurately reconstructing known genomes and capturing taxonomic diversity in a metagenomic sample. Multiple assemblers are recommended to be applied to a subset of samples to determine which has the best fit (Escobar-Zepeda et al., 2015; Ghurye et al., 2016; Boolchandani et al., 2019).

Binning classification can predict the taxonomical composition using information contained in reads (Figure 2). Reads are grouped to represent an individual genome or closely related genomes (Table 4). Binning can be done either with reads or assembled sequences and employs two strategies to obtain taxonomic assignment; (i) sequence composition classification and (ii) sequence alignment against references. The first is based on looking at genomic signatures using k-mers to identify evolutionary conservation among species and uses software such as TETRA (Teeling et al., 2004), MetaClusterTA (Wang et al., 2014), and PhylophytiaS (McHardy et al., 2007). Similarly, software like MaxBin (Wu et al., 2016) and Amphora2 (Wu and Scott, 2012) use k-mer signatures but also considers gene markers, GC content, and coverage information for binning on assembled sequences or reads. Another method based on reference read alignment is the Burrows-Wheeler Transform (Burrows and Wheeler, 1994) which indexes like BWA (Li and Durbin, 2009) or Bowtie (Langmead and Salzberg, 2012) and is fast and accurate in assessing species richness and abundance in WMS by aligning them to reference genomes. Genometa (Davenport and Tümmler, 2013) is a software that looks at OTUs in the WMS data and groups according to genomic islands and operons. Taxonomic classification in long reads using the potential coding regions to search in annotated protein databases using BLAST (Johnson et al., 2008) or Megan (Huson et al., 2007). Binning achieved after the assembly of contigs can lead to the generation of partial genomes of unknown or uncultured organisms. This can be used to perform similarity-based binning of other metagenomic datasets. If binning is achieved before assembly, it can reduce the complexity of an assembly-based approach and may reduce the computational demands. However, caution should be taken to ensure the validity of the genome bins as there is a chance of false assignments (Thomas et al., 2012; Bengtsson-Palme, 2018a). In this instance, CheckM (Parks et al., 2015) can be used to evaluate the quality of bins on metagenomic samples and assess the quality of assembled metagenomic genomes by estimation of completeness and contamination based on marker-gene validation (Du and Sun, 2022).

Taxonomic assignments of reads can also be done with genome annotation. This process identifies the coding regions and their location to determine gene function. Annotation of metagenomic sequences can be done in two ways, (i) existing pipelines, e.g., RAST (Aziz et al., 2008) and IMG (Markowitz et al., 2009), can be used for assembled genomes, and (ii) annotation can be performed on an entire community and relies on unassembled reads and short contigs. WMS data is annotated by identifying genes (feature prediction) and assigning putative gene functions and taxonomic neighbours (functional annotation; Thomas et al., 2012). Prediction tools such as MetaGeneAnnotator (Noguchi et al., 2008) or MetaGeneMark use internal information to classify sequences as either coding or non-coding and are done with a low error rate of 2% (Thomas et al., 2012).

Antibiotic resistance databases

The performance of antibiotic resistance gene prediction and taxonomic identification depends on the availability of accurate databases. Two types of public databases exist, generalized AMR databases, which include a wide range of ARGs and mechanistic information, and specialized databases, which provide extensive information on specific gene families, e.g., the β-lactamase family (Hendriksen et al., 2019).

Several bioinformatic tools and databases can be used in AMR research. They focus on fast and reliable predictions of ARB and their genetic determinants in complex communities (Van Camp et al., 2020). It is important to note that the quality and type of sequence play an essential role in AMR research. The bioinformatics tool’s robustness can differ when analysing sequences of low-quality (Hendriksen et al., 2019). Any database’s value depends on maintenance, curation, and continuous updates with detailed metadata. Correctly predicting ARGs from large-scale datasets with improved accuracy depends on input parameters such as e-values, bit scores, identity, and query coverage levels (Lal Gupta et al., 2020).

General antibiotic resistance databases

General antibiotic resistance databases incorporate a wide variety of ARGs and mechanistic information.

Established tools like ResFinder (Zankari et al., 2012) accept short reads and contigs to detect the presence of acquired ARGs. This can be done for complete or partial genome sequences and uses BLAST (Johnson et al., 2008), or KMA (Clausen et al., 2018) approaches (Zankari et al., 2012; Table 3). The specificity of using BLAST to detect antibiotic resistance genes depends on the selection criteria for gene length and percentage similarity. This method requires computational expertise and could be challenging to identify multiple copies of ARGs in de novo assembled genomes (Clausen et al., 2018; Hameed and Fatima, 2021). This can be solved by using longer reads which increases the cost.

In comparison, the KMA approach was developed to map the raw reads against redundant AMR databases as the traditional BLAST method remains too slow to map raw reads directly (Clausen et al., 2018; Hendriksen et al., 2019). This approach provides accurate bacterial genome analysis and can identify ARGs present in low abundance, which might be excluded in incomplete assemblies. There is also a chance of false positives due to sequencing errors which can be overcome by setting the minimum threshold for the number of reads needed for a positive outcome (Hendriksen et al., 2019; Hameed and Fatima, 2021).

The Comprehensive Antibiotic Resistance Database (CARD; Jia et al., 2016) for the identification of resistance genes, their products, and associated phenotypes. This database allows for the identification of chromosomal mutations that confer resistance to antibiotics. Similar to ResFinder (Zankari et al., 2012), the percentage coverage and identity can be adjusted to allow for the selection of new ARGs. The Resistance Gene Identifier (RGI), a tool for de novo annotation of genes, complete genomes, or genome assembled sequences for AMR, uses CARD’s curated database to predict ARGs in either DNA or amino acid format. This is because some sequences cannot be mapped due to divergence in the nucleotide sequences; however, amino acid sequences in this instance can be translated to proteins to characterize the resistome (Xavier et al., 2016; Lakin et al., 2019). CARD uses two prediction models: a protein homolog model based on functional resistance homologs and a protein variant model that detects mutations conferring resistance (Jia et al., 2016; Imchen et al., 2020). The application of ontologies such as CARD’s Antibiotic Resistance Ontology (ARO) has organized AMR information based on molecular determinants, resistance mechanisms, individual antimicrobials, antimicrobial targets, and drug classes, therefore, making this database favourable (McArthur and Tsang, 2017).

Other curated databases include; MEGARes (Lakin et al., 2017), ResFams (Gibson et al., 2015), NDARO (Sayers et al., 2021), FARME (Wallace et al., 2017), and MUSTARD (Ruppé et al., 2019), all of which are curated for the detection of AMR determinants in large metagenomic datasets (Arango Argoty, 2019; Arango-Argoty et al., 2020; Imchen et al., 2020; Lal Gupta et al., 2020). Interestingly, FARME, a functional antibiotic resistance metagenomic element database, focuses on AMR gene elements from environmental samples (soil, faecal matter, wastewater treatment plants, oral and aquatic biomes) rather than individual ARGs obtained from cultured clinical isolates, providing more access and analysis of non-clinical sectors. This database is compiled from publicly available DNA sequences from 30 functional metagenomic projects and their corresponding predicted proteins conferring AMR (Wallace et al., 2017). FARME contains over seven times the number of non-redundant protein sequences as compared to CARD (Wallace et al., 2017; Lal Gupta et al., 2020). It provides information on regulatory elements, MGEs, and predicted proteins flanking ARGs, which are conserved between functional metagenomic AMR sequences from soil biomes and pathogenic clinical isolates (Forsberg et al., 2012; Wallace et al., 2017). This allows for better insight into AMR in unculturable bacteria found in non-clinical settings (Wallace et al., 2017; Bengtsson-Palme, 2018b; Imchen et al., 2020).

ResFams is a curated database of protein families linked to their profile Hidden Markov Models (HMMs) associated with AMR function are trained using unique AMR protein sequences obtained from CARD, LaCED, and Lahey databases (Lal Gupta et al., 2020). This platform primarily focuses on pathogen-associated ARGs and provides a comprehensive view of resistomes in the environment and the evolution of resistant pathogens (Lal Gupta et al., 2020). ResFams were evaluated with functional metagenomic datasets and demonstrated improved sensitivity, and identified 64% more ARGs in soil and human gut microbiomes compared to BLAST-based searches of CARD and ARDB databases. This increased sensitivity is due to the HMM-based analysis. HMM, models are specific models constructed based on observed sequence variation across genes/protein families and capture possible variation for the families (Xavier et al., 2016). Although this approach is sensitive and can detect distant matches that BLAST-based approaches cannot, ResFams is computationally expensive and requires the user to provide local computational resources to run HMM-based searches (Xavier et al., 2016; Lakin et al., 2019).

Specific antibiotic resistance databases

Specific databases are created to meet specific needs and provide extensive information about a specific gene family. So far, only specialized databases are available for β-lactamase.

The β-lactamase database (BLAD) is a catalogue of resistance patterns from all classes of β-lactamases collected from NCBI and published data with crystal structures of proteins from the Protein Data Bank (PDB; Berman et al., 2000). Metagenomic sequences can be queried against BLAD to obtain basic information about β-lactamases, search for the gene of interest, identify resistance patterns and search and analyse the 3D structure of the β-lactamase (Danishuddin et al., 2013; Naas et al., 2017).

The CBMAR β-lactamase database (Srivastava et al., 2014) includes detailed biochemical and molecular information to understand known and unknown β-lactamase genes. This provides a cache of information on nucleotides, 3D protein structures, sequence alignments, and mutation profiles.

Similarly, LacED (Thai et al., 2009) is a specialized platform that contains mutational and structural data for TEM and SHV β-lactamases. This is extremely important to understand the evolution of multidrug-resistant pathogens as β-lactamase evolution, and the emergence of new enzymes affects the now available treatment for infections. Furthermore, these resistant genes are present on plasmids which help their spread to different biomes, challenging the treatment options in both clinical and non-clinical sectors (Dandachi et al., 2019).

Antibiotic resistance in non-clinical sectors

AMR is a natural phenomenon that predates the use of antibiotics. A wide array of novel and clinically characterized ARGs have been detected in environmental samples ranging from pristine environments to agricultural soil (Yang et al., 2019). The major driving force of AMR is the mis-use of antibiotics while other factors such as poor infrastructure, i.e., hygiene and sanitation play a role in maintaining AMR (Collignon et al., 2018; Essack and Sartorius, 2018). Approximately more than 60% of total antibiotic use occurs outside the field of human medicine (Prestinaci et al., 2015; Kumar et al., 2020; Polianciuc et al., 2020). The main non-clinical sectors that are involved in the development of AMR are animal production (food-animals and aquaculture), agriculture (plants), and the environmental (water and soil) compartment. These sectors are interconnected and facilitate the spread of pathogenic ARB within and between them and ultimately to humans. Table 5 summarizes the findings done on AMR in various non-clinical reservoirs through the use of metagenomics.

Table 5

ReservoirsIntroduction of antibioticsStudyReferences
SoilOrganic Manure/agricultural run-offIdentification of unknown tetracycline and sulphonamide resistance genes in forest and grassland soil.(Willms et al., 2019)
Manure applicationNon-manure amended soil displayed a larger proportion of antibiotic resistance bacteria but carried fewer genes.(Popowska et al., 2012)
Manure application /organic compostAntibiotic resistance genes from manure amended soil accounted for 70% of total resistance genes.(Su et al., 2014)
Antibiotic-producing soil bacteriaThirteen antibiotic resistance genes and two bifunctional proteins conferring resistance to ceftazidime and β-lactamase.(Donato et al., 2010)
Antibiotic-producing soil bacteriaTotal of eleven ARGs in soil with little human influence.(Torres-Cortés et al., 2011)
Fertilizer/ manure applicationTetracycline-resistant bacteria in fertilized soil were three times higher than manure amended soil.(Kang et al., 2018)
A fertilizer with heavy metalsARGs could potentially increase by the use of fertilizer containing heavy metals.(Lin et al., 2016)
Manure applicationManure-borne bacteria contribute to the increase of ARGs whereas indigenous soil bacteria prevent the dissemination of ARGs from manure to soil.(Chen et al., 2017)
Pollution from a petrochemical plantThe soil was more abundant in ARGs than less contaminated soil.(Chen et al., 2017)
PlantsManure applicationARGs on lettuce and endive were grown in manure amended soil. ARGs were detected in the endophytes and phyllosphere of the plants.(Wang et al., 2015)
Conventionally produced (manure application) / organic fertilizersA higher abundance of diverse ARGs and MDR genes were present in the endophyte and phyllosphere of lettuce grown organically.(Zhu et al., 2017)
Manure application /fertilizerMore than 50 unique ARGs were detected in the phyllosphere of maize after manure application.(Chen et al., 2018)
Compost/ improper sanitation/ run-off from farms2015 to 2020, five outbreaks linked to green leafy vegetables, two linked to sprouts, and one linked a chopped salad kit in the USA.(CDC, 2021)
Food-animalsAnimal feed / infection preventionTetracycline resistant E. coli strain present in the gut of chickens receiving food supplemented with tetracycline.(Marshall and Levy, 2011)
Animal feed/growth promoter / indirect and direct contactARBs were found in each processing step of the beef production chain.(Yang et al., 2016)
Animal feedA total of 495 bacterial species, 50 ARGs were detected in the organic animal feed.(Eckstrom et al., 2019)
Therapeutic and non-therapeuticAntibiotics used in food animals leave antibiotic residues in eggs, milk, and meat products.(Oliveira et al., 2020)
Animal feedFaecal samples from large-scale swine farms showed a total of 146 ARGs.(Zhao et al., 2018)
Animal feedTaxonomic composition analysis showed a decrease of Proteobacteria and an increase of Actinobacteria in the gut microbiome of male broilers.(She et al., 2018)
WaterEffluent dischargeThe abundance of pathogens and ARGs increased in the effluent after wastewater treatment.(Tang et al., 2016)
ChlorinationARGs are found in high amounts in drinking water suggesting chlorination concentrates ARGs.(Shi et al., 2013)
Run-off / dischargeEffluent from municipal wastewater treatment plants and pharmaceutical manufacturing plants releases antibiotics into natural water bodies.(Guo et al., 2017)
Reuse of water / potable waterAn increase in ARGs in reclaimed potable water samples and associations between 193 ARGs and plasmid-associated genes.(Garner et al., 2018)
Sewage run-offThe distinct taxonomic composition of bacterial species found in sewage and the sea suggests sewage is diluted in environmental freshwater.(Fresia et al., 2019)
Aqua-cultureAnimal feedFifty-one different ARGs conferring resistance against 24 different antibiotic types with MDR genes located on plasmid sequences.(Tyagi et al., 2019)
Animal feed /therapeuticTetracycline resistance genes were found in 81% of samples taken from a Chilean salmon farm.(Shah et al., 2014)
Animal feed/ therapeuticAn abundance of tetracycline, sulphonamide, and beta-lactam resistance genes which strongly correlated to silver and mercury resistance genes were found indicating co-resistance and co-regulation.(Yuan et al., 2019)

Studies focusing on the introduction of antibiotics and AMR in each non-clinical sector.

Food-animals

Each year, half of the antibiotics produced are used for industrial farm animal production (Schmieder and Edwards, 2012). The global demand for animal protein has directly driven the use of antibiotics in livestock. The food-animal industry is a hotspot for antibiotic usage, estimated to increase by 67% in highly populated countries (Van Boeckel et al., 2015). Antibiotics are administered as metaphylactic and prophylactic treatment for the whole flock or herd to prevent the spread of infectious diseases. Furthermore, it is used as a growth-promoting agent across the globe (Van Boeckel et al., 2015; Founou et al., 2016; Robinson et al., 2016).

Historically, the negative effects of antibiotic usage in livestock production have been overlooked because of the drive to keep up with consumer demand (Robinson et al., 2016; Ramtahal et al., 2022). This has led to the prolonged use and dosage, creating an ideal condition for the emergence of ARB and ARGs (Marshall and Levy, 2011; Madec and Haenni, 2018).

AMR can spread through direct and indirect contact, which was confirmed by finding the same ARBs in the poultry farm and the caretaker (Wee et al., 2020). Antibiotics and their residues can be found in most connecting areas such as soil and plants, and animal products which can lead to health risks when consumed due to the direct toxicity, allergic reactions, carcinogenic effects, and disturbance of the beneficial microbiota in children, those who are pregnant, the elderly, and people living with immunocompromised people (Oliveira et al., 2020).

Areas such as the EU and United States have made efforts to ban the use of antibiotics as growth enhancers. However, they are still widely used in South Africa, Brazil, India, China, and Russia (Van Boeckel et al., 2015; Robinson et al., 2016). This strategy is based on the assumption that susceptible strains can outnumber AMR microorganisms if the advantage of becoming resistant is decreased (Andersson and Hughes, 2011). This was supported by a post-ban study conducted in Denmark and Norway. Interestingly, even though ARB decreased, the broilers that were studied were still colonized by ARB years after the ban came into effect. It is such that the ban on antibiotics is still up for debate as decreasing antibiotic use could increase the frequency of bacterial infections (Holmes et al., 2016; Andersson and Hughes, 2019).

The World Health Organization has recommended that farmers and the food industry stop using antibiotics for the growth and prevention of disease in healthy animals. Farms should consider the employment of proper preventive measures and farming practices to limit the risk of bacterial infections instead of the routine antibiotic dosing of animals to prevent disease. If this is done correctly, we could see the decreased use of antibiotics and the alleviation of (Marshall and Levy, 2011; Robinson et al., 2016; Aslam et al., 2018; Van Gompel et al., 2020).

Water

Water represents one of the most important bacterial habitats and is a major pathway and reservoir for disseminating microorganisms and AMR. This is due to the spread of AMR by humans, and animals, or contamination of the environment by pathogens (Hutton and Chase, 2017; Miłobedzka et al., 2021).

Grey water

Wastewater treatment plants (WWTPs) collect chemical pollutants, including antibiotics, which are disposed of from households, hospitals, and factories. The WWTP aims to remove contaminants before disposing water in natural environments, e.g., streams, rivers, and lakes. While the reuse and reclamation of water through WWTPs reduce water shortage, they assist in the spread and emergence of AMR (Michael et al., 2013). This impacts the economy and society of countries affected by AMR (Bürgmann et al., 2018; Prüss-Ustün et al., 2019). It could be possible that proper sanitation and safe drinking water in the high-, middle-and low-income countries affect the emergence and spread of AMR. Although the removal of contaminants is successful, some ARB are not removed by the WWTP and are released with the effluent, which can spread to other natural environments (Meena et al., 2015).

Since WWTPs receive a diverse range of ARB and their genes from different sources. This ARB accumulates in the WWTP, and these high concentrations provide a selective pressure that facilitates the emergence of new ARB strains (Bürgmann et al., 2018). This is seen in wastewater treatment that uses biological treatment processes such as activated sludge, creating a favourable environment for ARB, ARGs, and their transmission (Guo et al., 2017). Contradictorily, it was hypothesized that the WWTP is efficient in removing ARBs which was seen in a study of the Dutch (Pallares-Vega et al., 2019). It could be argued that AMR in WWTP decreased since the Netherlands decreased antibiotic use as a whole. While this can be observed, antibiotic residues are not fully removed from the WWTP and persist in lower concentration when released with the effluent which serves as a selection pressure for ARBs (Keen and Fugère, 2017). WWTPs should revise their treatment options to manage AMR concerns while also providing sanitary water to communities. Policies about sludge and biosolid disposals are needed along with better infrastructure which will help decrease the spread of AMR to other sectors (Bürgmann et al., 2018; Kraemer et al., 2019).

Blue/green water

Freshwater environments are among the natural environments susceptible to contamination with antibiotics. These are released through different sources such as agricultural run-off, sewage discharge, and leaching from farms (Hutton et al., 2007; Hutton and Chase, 2017). The combination of antibiotics with a high density of bacteria provides a favourable environment for the development of ARGs (Marti et al., 2014). ARB carrying resistance genes can persist in freshwater and ultimately return to humans and animals by horizontal gene transfer (HGT; Chen et al., 2017; Mirza et al., 2020). Freshwater serves as a source of drinking water, recreational purposes and for agricultural practices, which use 70% of all freshwater, e.g., aquifers, streams, and lakes (Bürgmann et al., 2018). Continuous use of freshwater encourages the spread of AMR and the emergence of new pathogens, which leads to an increased risk of infection and prolonged and untreatable infections in humans, animals, and the environment (Nnadozie and Odume, 2019).

To date (Table 5), there is evidence showing that the consumption of contaminated water poses health risks (Jjemba et al., 2010; Shi et al., 2013; Sanganyado and Gwenzi, 2019) and that the WWTP has the potential to contribute to the dissemination of AMR in receiving rivers (Tang et al., 2016; Guo et al., 2017). Unfortunately, information on antibiotics in streams, lakes, beaches, pond water, surface and groundwater, and drinking water is scarce as it is presumed that antibiotic concentrations are naturally low in these areas. However, low antibiotic concentrations still select for ARB (Nnadozie and Odume, 2019). Studies also focus on ARB and genes at discharge points, i.e., receiving river or effluent, of the WWTP. At this point, the concentration of AMR is high but declines further away from the WWTP, downstream, which could be a result of degradation through biological or non-biological processes, uptake by aquatic microorganisms through HGT, diluted, absorption onto particulate matter, or transportation into different water systems (Anyaduba, 2016; Mirza et al., 2020). More focus should be on AMR within each freshwater environment, its link to different sectors, AMR after the WWTP discharge points and in potable water to gain a more holistic view (Mirza et al., 2020).

Aquaculture

Aquaculture refers to the breeding, rearing, and harvesting of aquatic organisms in different water environments. This term is commonly used to describe fish farming. Aquaculture contributes to more than half of the world’s seafood consumption and production and has increased by 6% a year since 2001 (Rocha et al., 2022). The aquacultural sector is a complex interconnected system influenced by the environment, cultural, human, and economic factors (Brunton et al., 2019; Schar et al., 2021).

Like livestock production, antibiotics are used in aquaculture as a therapeutic agent to prevent and treat bacterial outbreaks. High concentrations are used due to high stocking densities, and the lack of individual treatment (Yukgehnaish et al., 2020).

Antibiotics are mixed with feed before being administered to animals or can be directly applied to the aquatic environment. This leads to the dispersal and leaching of antibiotics into the environment (Brunton et al., 2019). The gut microbiota of fish is also affected by the continuous use of antibiotics which alter the benefit of host-microbiota relationships allowing the gut microbiota to become resistant to antibiotics that are administered (Tyagi et al., 2019). Aquaculture is also affected by toxic materials such as silver and mercury, which co-resist and co-regulate with antibiotics. These toxic materials accumulate in fish bodies and the aquatic food web. This toxicity can be passed to humans and animals through consumption. It can affect the nervous system leading to death, as discussed by (Sonone et al., 2020), who reviewed the various heavy metal sources, i.e., agricultural activities, electronic waste, mining, industrial effluents, power plants, and biomedical waste and its role in degrading the aquaculture population, causing physical deformities in organisms and polluting the aquatic environment.

Several studies (Table 5) have shown that the excess use of antibiotics has led to the emergence of AMR. It cannot be ignored that the aquaculture industry is integrated with sewage, industrial wastewater, and land agriculture as manure and other agricultural residues are used in fish food. Areas in Western Europe and North America have banned antibiotics as a growth promoter, but it is still used as a therapeutic agent in fish food. Since 70 to 80% of antibiotics given to fish are excreted into the water, the entire body of water is exposed. This leads to leaching into different water environments and sediments, which provide long-term selection pressure in the aquatic environment, favoring horizontal gene transfer and the spread of AMR, as discussed by (Taylor et al., 2011; Watts et al., 2017; Mirza et al., 2020). Such policies on aquacultural practices should be reviewed to decrease the incidence of ARB and infections (Marti et al., 2014; Pérez-Sánchez et al., 2018; Preena et al., 2020a; Zhao et al., 2020).

Soil

Microbes in soil are of great importance as they break down organic matter, recycle nutrients, bioremediation, and produce antibiotics (Chandra and Kumar, 2017). Variation in biotic, e.g., plants, animals and other bacteria, and abiotic, e.g., water, soil and atmosphere, conditions of soil cause the residing microbes to adapt and develop strategies for survival and successful reproduction. Antibiotic production is the most powerful adaptation strategy from soil microbes to inhibit the growth of competing microbes (Chandra and Kumar, 2017). This strategy has led to the natural development of AMR. Antibiotics and disinfectants used in medicine, agriculture, and aquaculture have been the driver of AMR through manure/fertiliser application on soil, irrigation, and run-off. Consequences of these are the emergence of new ARGs and the spread of pathogenic ARB to other environments, humans, and animals (Martínez, 2008; Wang et al., 2015; Chandra and Kumar, 2017).

Although AMR research has shown that WWTP effluents can still contain a high level of ARGs and bacteria after disinfection, there is not much research on the sludge and biosolids applied to soil and how this changes the bacterial community and resistome (Anyaduba, 2016; Nnadozie and Odume, 2019). The application of manure/fertiliser/sludge to soil should also be looked into to determine if there is any difference in soil and crop quality as they are used in agricultural and remediation practices. Studies have shown the prevalence of AMR in soil and how bioremediation and pollutants have facilitated its spread (Table 5).

Plants

The agricultural industry has been valued at an estimated US$ 3.2 trillion worldwide and accounts for the largest share of the GDP and employment in developing countries and underdeveloped nations. Unfortunately, the agriculture industry has been suffering due to population growth, pest resistance, and the burden of natural resources (Raman, 2017).

Manure is commonly used as fertilizer on vegetable farmlands to enhance crop yield. This enhances AMR as 30–90% of antibiotics consumed by animals are excreted in faeces which is used as manure. Using sewage sludge and water for irrigation contributes to rising AMR in agriculture (Minden et al., 2017; Grenni et al., 2018; Spielmeyer, 2018). Crops can take up antibiotics through water transport and passive absorption (Hu et al., 2010). Antibiotic-resistant bacteria found in soil can reach the interior of fruits and vegetables either through germination seeds in contaminated soil or the direct transmission of bacteria from the soil to the plant via irrigation. Both these routes can cause the bacteria to colonize in the roots and the edible parts of the plant (Hirniesen et al., 2012; Alegbeleye et al., 2018; Sundin and Wang, 2018).

These edible plants are consumed by animals or sold as ready-to-eat products such as bagged salads or pre-cut vegetables, which may enhance the proliferation and survival of ARB (Jung et al., 2014; Rajwar et al., 2016; Wadamori et al., 2017; Mir et al., 2018). This is due to many factors such as hygiene standards, irrigation water, temperature, storage, pH, soil, manure, and antimicrobial agents. Since these products rarely go through a heating step before consumption, it increases the risk of consumers being exposed to high numbers of different pathogenic ARB (Mir et al., 2018). This poses a threat to society as consumers who opt for ready-to-eat products are faced with the risk of acquiring foodborne illnesses which could be untreatable with antibiotics. Table 5 summarizes studies showing that manure, fertilizers, compost, and irrigation are all routes for the transmission and dissemination of AMR to plants.

Antibiotic resistance and its transmission

The use and misuse of antibiotics in humans, animals, and the environment have been linked to the emergence of ARB in each of these sectors. Antibiotic residues, ARB, and ARGs, are spread to different environments (Miłobedzka et al., 2021). Most antibiotics are excreted in animal urine and faeces (Jayalakshmi et al., 2017). This is introduced to the environment directly by manure amendment, which provides nutrients for crops, or through water which is the primary link between humans, animals, and nature (Miłobedzka et al., 2021).

Since antibiotics are used in the food animal industry, ARB can spread via run-off, which can go into lakes, streams, and rivers. Water from these areas is used for irrigation for agricultural practices. This leads to the spread of ARB from food animals to water, soil, and crops (Figure 1). This interconnectedness facilitates the exchange of AMR and allows for the emergence of new genetic determinants. This emergence and spread can occur by HGT (Schmieder and Edwards, 2012; Pan and Chu, 2017; Sarkar et al., 2018; Mahmood et al., 2019).

Figure 2

Figure 2

Workflow for determining antibiotic resistance genes, mobile genetic elements, horizontal gene transfer events and taxonomic composition using metagenomic data. The workflow is indicated by solid lines with example tools and databases annotated in the boxes. The complete lists are described in Tables 3, 4.

Horizontal gene transfer

ARB exists across animal, human and environmental sectors. The evolution of AMR and the dynamics of ARGs spread across these sectors are critical for predicting emerging pathogens and controlling AMR dissemination. HGT is primarily responsible for rapid AMR spread and dissemination of ARGs between bacteria and across species (Sun et al., 2020; Vinayamohan et al., 2022). AMR can be established either vertically by point mutations or horizontally by acquiring MGEs such as plasmids and transposons (Vinayamohan et al., 2022).

HGT is a process that introduces variation in the bacterial genome under natural selection. It is an essential process that allows bacteria to adapt to environmental changes such as exposure to antibiotics (Hall et al., 2020).

HGT predates the production and use of antibiotics. Unfortunately, the increased usage of antibiotics has put selective pressure on bacterial strains. This has sped up the process of bacteria becoming antibiotic-resistant by various genotypic and phenotypic resistance mechanisms. This diversification can be seen in the resistance to antibiotics of different classes (Power, 2019).

Genes that confer resistance can either be intrinsic or extrinsic (Ali et al., 2018). Intrinsic resistance is the bacterium’s innate ability to resist the antimicrobial agent’s activity. This occurs through structural or functional characteristics. Extrinsic or acquired resistance occurs when a previously susceptible bacterium obtains the ability to resist the activity of the antimicrobial agent. This can result from the mutation of genes involved in normal physiological processes, the acquisition of foreign resistance genes, or a combination of both (Blair et al., 2015; von Wintersdorff et al., 2016; Kumar, 2017).

Pathogenic bacteria can acquire resistance via HGT. This has caused AMR to spread from commensal and environmental bacteria to pathogenic ones. ARGs are often carried in MGEs, e.g., plasmids, transposons, or integrons which act as vectors for transferring genetic information between bacterial cells. The three classical pathways of HGT are conjugation, transformation, and transduction. These HGT mechanisms lead to increased bacteria fitness, which is essential for survival in the presence of antibiotics (Capita and Alonso-Calleja, 2013; Munita and Arias, 2016; von Wintersdorff et al., 2016; Kumar, 2017; Vinayamohan et al., 2022).

Recent studies have identified other mechanisms by which DNA can transfer between hosts. Gene transfer agents are DNA-containing particles, similar to phages, but cannot carry genes for particle production (Redfield and Soucy, 2018). Genes can also be transferred between bacteria that form intercellular connections by nanotubes or membrane fusion (Hall et al., 2017, 2020). Some bacteria can also release DNA-containing membrane-bound vesicles that carry genetic information to new hosts. Interestingly, some of these mechanisms are not under bacterial control but are controlled by semi-autonomous segments of DNA (Hall et al., 2017).

Bioinformatic approaches such as MetaCHIP (Song et al., 2019) and WAAFLE (WAAFLE, 2022) can be used to detect HGT events in metagenomic sequences. These processes involve an all-against-all BLASTN or BLASTX of genes within assembled contigs, and potential HGT events are determined based on genes with the best hits in other taxonomic groups (Douglas and Langille, 2019). These approaches are limited by sequencing technologies and assemblers, which often fail to assemble long regions with high sequence similarity (Song et al., 2019). It can be assumed that inferring HGT events will become easier as the assembly of metagenomic sequences improves (Douglas and Langille, 2019).

Mobile genetic elements

Specialized MGEs mediate HGT between bacteria. These play an important role in bacterial ecology and evolution. Several factors impact HGT in bacteria, such as gene expression, protein connectivity, and biochemical properties; however, the range of genes carried on MGEs remains unknown (Rodríguez-Beltrán et al., 2020).

MGEs are segments of DNA that encode enzymes and other proteins that mediate the movement of DNA within or between bacterial cells (Hall et al., 2017; Emamalipour et al., 2020). MGEs carry genes other than those necessary for transfer and replication (Hall et al., 2017). These accessory genes do not play a role in their vertical or horizontal transmission but affect the success of MGEs. The range of accessory genes encoded by MGEs and their ability to be phenotypically expressed in different genetic backgrounds are critical in its evolution (Rodríguez-Beltrán et al., 2020).

The most important elements that play a role in HGT are the conjugative and mobilized elements. Mobilized elements contain all genetic information required to transfer from one bacterium to another, while conjugative elements use conjugation functions of elements such as plasmids and transposons to transfer from one host to another. These MGEs can carry multi-drug resistant (MDR) plasmids and give rise to MDR bacteria via HGT (Vinayamohan et al., 2022). Bacteriophages also play a role in the spread of DNA by the transduction process. Some elements can translocate to new sites on a genome but cannot transfer to a new bacterial host. These include the transposons and mobile integrons (Van Hoek et al., 2011; Hall et al., 2017).

Detecting mobile genetic elements

The ability of MGEs containing ARGs to spread is controlled by the genetic elements and host factors. The genetic elements and host factors control the ability of MGEs containing ARG, but most important is the selective pressure in the environment. When antibiotics are present in the environment, there is strong selective pressure on the spread of resistance. Elements that promote resistance will be selected for, and those stopping the spread of mobile elements will be selected against (von Wintersdorff et al., 2016; Hall et al., 2017; Partridge et al., 2018). To control the spread of AMR, it is important to understand MGEs and the ecology of the environments in which they spread easily (Van Hoek et al., 2011; Slizovskiy et al., 2020).

Several databases such as ACLAME (Leplae et al., 2004, 2010), ISFinder (Siguier et al., 2006), INTEGRALL (Moura et al., 2009), and ICEberg (Liu et al., 2019) and web-based tools such as ISSaga2 (Varani et al., 2011), MobilomeFinder (Ou et al., 2007), oriTfinder (Li et al., 2018), TnpPred (Riadi et al., 2012) and MobileElementFinder (Johansson et al., 2021) were developed to identify MGEs in metagenomic datasets (Peng et al., 2021). These web-based tools and databases range from detecting plasmids, integrons, transposons, prophages, and insertion sequences.

Overall, there is no single solution for detecting MGEs as these databases, and web-based tools are incomplete and biased towards pathogens studied extensively. This makes it difficult and time-consuming as most web-based tools and databases need to be evaluated first and used in combination to provide an all-inclusive view on MGEs associated with AMR (Rankin et al., 2011; Brito, 2021). This information can be applied to clinical and non-clinical sectors to potentially influence policies and strategies which can help control the spread of AMR (Sun et al., 2020).

Outlook

Antibiotic resistance affects animals, plants, the environment and human health (Malagón-Rojas and Lagos, 2020). Monitoring and understanding the prevalence, mechanisms and spread of AMR is a priority in non-clinical sectors (Boolchandani et al., 2019). Since measures against AMR are focused on the clinical health sector, the environment is under-represented which allows for significant gaps in data, research and control strategies therefore failing to link the effects of AMR from one sector to the other (Malagón-Rojas and Lagos, 2020).

Metagenomic sequencing allows for the detection of known and novel ARGs, HGT events, and associated MGEs in complex communities providing in-depth information on their prevalence, distribution, and transmission in non-clinical sectors such as the animal industry, agriculture, and the environment (Forbes et al., 2017).

While metagenomics is promising, its implementation in non-clinical sectors is still in its early stages. Metagenomic sequencing creates a vast amount of data and requires specialized bioinformatic expertise, which makes this a costly approach (Forbes et al., 2017; Oniciuc et al., 2018). To date, only a few countries and laboratories have the resources and expertise to use metagenomic approaches as a surveillance system for AMR. In addition, bioinformatics methodologies need to be standardized and constantly updated and curated to allow for an accurate comparison between various samples, as the selection of bioinformatic tools and AMR databases does have a significant impact on the results (Hendriksen et al., 2019; Lal Gupta et al., 2020).

Nevertheless, using bioinformatic tools and databases can assist in identifying risky practices, effects of antibiotic usage, and hotspots of AMR. This can facilitate the design of new policies to control the spread of AMR between clinical and non-clinical sectors (Shen and Zhao, 2018; Nijsingh et al., 2019). Such knowledge and practices are urgent as these sectors, including the public health sector, are faced with a future of untreatable infections, increasingly costly medicinal treatment, a higher cost of living, and an increased mortality rate (Bengtsson and Greko, 2014; Littmann et al., 2015; Nijsingh et al., 2019).

Conclusion

The application of metagenomic and bioinformatic approaches to AMR research can provide fast and reliable predictions of AMR and antibiotic use in various non-clinical sectors. These hold a great promise for understanding AMR molecularly, predicting outbreaks and transmission, and emerging pathogens. This information can lead to better policymaking in each sector and decrease the incidence of infections as conditions for the animal industry, agriculture, and the environmental sectors are improved.

Funding

This work is based on the research supported by wholly/in part by the National Research Foundation of South Africa (Grant Numbers: 120192).

Publisher’s note

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

Statements

Author contributions

SP and TA conceived, designed and wrote the manuscript. AU and DC-F edited and proof-read the manuscript. All authors contributed to the article and approved the submitted version.

Conflict of interest

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

References

  • 1

    AhmedA. M.MotoiY.SatoM.MaruyamaA.WatanabeH.FukumotoY.et al. (2007). Zoo animals as reservoirs of gram-negative bacteria harboring integrons and antimicrobial resistance genes. Appl. Environ. Microbiol.73, 66866690. doi: 10.1128/AEM.01054-07

  • 2

    AlegbeleyeO. O.SingletonI.Sant AnaA. S. (2018). Sources and contamination routes of microbial pathogens to fresh produce during field cultivation: a review. Food Microbiol.73, 177208. doi: 10.1016/j.fm.2018.01.003

  • 3

    AliA. S.Al-QahtaniA.NabiS.AltafK.SonekhiG. B.AhmedJ.et al. (2018). A short review on antibiotics and ever-changing microbial resistance mechanisms. British Journal of Pharmacy.3:1. doi: 10.5920/bjpharm.2018.XX

  • 4

    AlnebergJ.BjarnasonB. S.De BruijnI.SchirmerM.QuickJ.IjazU. Z.et al. (2014). Binning metagenomic contigs by coverage and composition. Nat. Methods, 11441146. doi: 10.1038/nmeth.3103

  • 5

    AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403410. doi: 10.1016/S0022-2836(05)80360-2

  • 6

    AnderssonD. I.HughesD. (2011). Persistence of antibiotic resistance in bacterial populations. FEMS Microbiol. Rev.35, 901911. doi: 10.1111/j.1574-6976.2011.00289.x

  • 7

    AnderssonD. I.HughesD. (2019). Selection and transmission of antibiotic-resistant bacteria. Microbial Transmission.35, 117137.

  • 8

    AnyadubaT. (2016). Fate and transport of antibiotic resistance genes in aquatic ecosystem. 116.

  • 9

    Arango ArgotyG. A. (2019). Computational tools for annotating antibiotic resistance in metagenomicUnited States: Data: Virginia Tech.

  • 10

    Arango-ArgotyG.GuronG.GarnerE.RiquelmeM. V.HeathL.PrudenA.et al. (2020). ARGminer: a web platform for the crowdsourcing-based curation of antibiotic resistance genes. Bioinformatics36, 29662973. doi: 10.1093/bioinformatics/btaa095

  • 11

    AslamB.WangW.ArshadM. I.KhurshidM.MuzammilS.RasoolM. H.et al. (2018). Antibiotic resistance: a rundown of a global crisis. Infection and drug resistance.11:1645. doi: 10.2147/IDR.S173867

  • 12

    AzizR. K.BartelsD.BestA. A.DeJonghM.DiszT.EdwardsR. A.et al. (2008). The RAST server: rapid annotations using subsystems technology. BMC Genomics9, 115. doi: 10.1186/1471-2164-9-75

  • 13

    BallouxF.BrynildsrudO. B.Van DorpL.ShawL. P.ChenH.HarrisK. A.et al. (2018). From theory to practice: translating whole-genome sequencing (WGS) into the clinic. Trends Microbiol.26, 10351048. doi: 10.1016/j.tim.2018.08.004

  • 14

    BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.19, 455477. doi: 10.1089/cmb.2012.0021

  • 15

    BastidaF.HernándezT.GarcíaC. (2014). Metaproteomics of soils from semiarid environment: functional and phylogenetic information obtained with different protein extraction methods. J. Proteome101, 3142. doi: 10.1016/j.jprot.2014.02.006

  • 16

    BengtssonB.GrekoC. (2014). Antibiotic resistance—consequences for animal health, welfare, and food production. Ups. J. Med. Sci.119, 96102. doi: 10.3109/03009734.2014.901445

  • 17

    Bengtsson-PalmeJ. (2018a). Strategies for taxonomic and functional annotation of metagenomes. Netherlands: Metagenomics: Elsevier. p. 5579.

  • 18

    Bengtsson-PalmeJ. (2018b). The diversity of uncharacterized antibiotic resistance genes can be predicted from known gene variants—but not always. Microbiome.6:125. doi: 10.1016/B978-0-08-102268-9.00003-3

  • 19

    BermanH. M.WestbrookJ.FengZ.GillilandG.BhatT. N.WeissigH.et al. (2000). The protein data bank. Nucleic Acids Res.28, 235242. doi: 10.1093/nar/28.1.235

  • 20

    BlairJ. M.WebberM. A.BaylayA. J.OgboluD. O.PiddockL. J. (2015). Molecular mechanisms of antibiotic resistance. Nat. Rev. Microbiol.13, 4251. doi: 10.1038/nrmicro3380

  • 21

    BoolchandaniM.D’SouzaA. W.DantasG. (2019). Sequencing-based methods and resources to study antimicrobial resistance. Nat. Rev. Genet.20, 356370. doi: 10.1038/s41576-019-0108-4

  • 22

    BreitwieserF. P.LuJ.SalzbergS. L. (2019). A review of methods and databases for metagenomic classification and assembly. Brief. Bioinform.20, 11251136. doi: 10.1093/bib/bbx120

  • 23

    BritoI. L. (2021). Examining horizontal gene transfer in microbial communities. Nat. Rev. Microbiol.19, 442453. doi: 10.1038/s41579-021-00534-7

  • 24

    BruntonL. A.DesboisA. P.GarzaM.WielandB.MohanC. V.HäslerB.et al. (2019). Identifying hotspots for antibiotic resistance emergence and selection, and elucidating pathways to human exposure: application of a systems-thinking approach to aquaculture systems. Sci. Total Environ.687, 13441356. doi: 10.1016/j.scitotenv.2019.06.134

  • 25

    BürgmannH.FrigonD.GazeH.WM.ManaiaC.PrudenA.et al. (2018). Water and sanitation: an essential battlefront in the war on antimicrobial resistance. FEMS Microbiol. Ecol.94:fiy101. doi: 10.1093/femsec/fiy101

  • 26

    BurrowsM.WheelerD.Technical report 124. Palo Alto, CA: Digital Equipment Corporation. (1994).

  • 27

    CabelloF. C.GodfreyH. P. (2016). Even therapeutic antimicrobial use in animal husbandry may generate environmental hazards to human health. Environ. Microbiol.18, 311313. doi: 10.1111/1462-2920.13247

  • 28

    CapitaR.Alonso-CallejaC. (2013). Antibiotic-resistant bacteria: a challenge for the food industry. Crit. Rev. Food Sci. Nutr.53, 1148. doi: 10.1080/10408398.2010.519837

  • 29

    CDC. Gov. list of selected multistate foodborne outbreak investigations. (2021). Available from:http://www.cdc.gov/foodsafety/outbreaks/multistate-outbreaks/outbreaks-list.html.

  • 30

    ChandraN.KumarS. (2017). “Antibiotics producing soil microorganisms” in Antibiotics and antibiotics resistance genes in soils. eds. MuhammadH. Z.VladimirS.VarmaA. (Berlin: Springer), 118.

  • 31

    CheY.XiaY.LiuL.LiA.YangY.ZhangT. (2019). Mobile antibiotic resistome in wastewater treatment plants revealed by Nanopore metagenomic sequencing. Microbiome.7:44. doi: 10.1186/s40168-019-0663-0

  • 32

    ChenQ.-L.AnX.-L.LiH.ZhuY.-G.SuJ.-Q.CuiL. (2017). Do manure-borne or indigenous soil microorganisms influence the spread of antibiotic resistance genes in manured soil?Soil Biol. Biochem.114, 229237. doi: 10.1016/j.soilbio.2017.07.022

  • 33

    ChenQ.-L.AnX.-L.ZhengB.-X.MaY.-B.SuJ.-Q. (2018). Long-term organic fertilization increased antibiotic resistome in phyllosphere of maize. Sci. Total Environ.645, 12301237. doi: 10.1016/j.scitotenv.2018.07.260

  • 34

    ChenY.-H.ChiangP.-W.RogozinD. Y.DegermendzhyA. G.ChiuH.-H.TangS.-L. (2021). Salvaging high-quality genomes of microbial species from a meromictic lake using a hybrid sequencing approach. Communications Biology.4:996. doi: 10.1038/s42003-021-02510-6

  • 35

    ChenB.HeR.YuanK.ChenE.LinL.ChenX.et al. (2017). Polycyclic aromatic hydrocarbons (PAHs) enriching antibiotic resistance genes (ARGs) in the soils. Environ. Pollut.220, 10051013. doi: 10.3389/fmicb.2017.01133

  • 36

    ChenZ.YuD.HeS.YeH.ZhangL.WenY.et al. (2017). Prevalence of antibiotic-resistant Escherichia coli in drinking water sources in Hangzhou City. Front. Microbiol.8:1133.

  • 37

    ClausenP. T. L. C.AarestrupF. M.LundO. (2018). Rapid and precise alignment of raw reads against redundant databases with KMA. BMC Bioinformatics.19:307. doi: 10.1186/s12859-018-2336-6

  • 38

    ClausenP. T. L. C.ZankariE.AarestrupF. M.LundO. (2016). Benchmarking of methods for identification of antimicrobial resistance genes in bacterial whole genome data. J. Antimicrob. Chemother.71, 24842488. doi: 10.1093/jac/dkw184

  • 39

    CollignonP.BeggsJ. J.WalshT. R.GandraS.LaxminarayanR. (2018). Anthropological and socioeconomic factors contributing to global antimicrobial resistance: a univariate and multivariable analysis. The Lancet Planetary Health.2, e398e405. doi: 10.1016/S2542-5196(18)30186-4

  • 40

    CollineauL.BoerlinP.CarsonC. A.ChapmanB.FazilA.HetmanB.et al. (2019). Integrating whole-genome sequencing data into quantitative risk assessment of foodborne antimicrobial resistance: a review of opportunities and challenges. Front. Microbiol.10:1107. doi: 10.3389/fmicb.2019.01107

  • 41

    CroftsT. S.GasparriniA. J.DantasG. (2017). Next-generation approaches to understand and combat the antibiotic resistome. Nat. Rev. Microbiol.15:422. doi: 10.1038/nrmicro.2017.28

  • 42

    DandachiI.ChaddadA.HannaJ.MattaJ.DaoudZ. (2019). Understanding the epidemiology of multi-drug resistant gram-negative bacilli in the Middle East using a one health approach. Front. Microbiol.10, 139. doi: 10.3389/fmicb.2019.01941

  • 43

    DanishuddinM.Hassan BaigM.KaushalL.KhanA. U. (2013). BLAD: a comprehensive database of widely circulated beta-lactamases. Bioinformatics29, 25152516. doi: 10.1093/bioinformatics/btt417

  • 44

    DavenportC. F.TümmlerB. (2013). Advances in computational analysis of metagenome sequences. Environ. Microbiol.15, 15. doi: 10.1111/j.1462-2920.2012.02843.x

  • 45

    DaviesJ.DaviesD. (2010). Origins and evolution of antibiotic resistance. Microbiol. Mol. Biol. Rev.74, 417433. doi: 10.1128/MMBR.00016-10

  • 46

    De MaioN.ShawL. P.HubbardA.GeorgeS.SandersonN. D.SwannJ.et al. (2019). Comparison of long-read sequencing technologies in the hybrid assembly of complex bacterial genomes. Microbial genomics.5:294. doi: 10.1099/mgen.0.000294

  • 47

    Di BellaJ. M.BaoY.GloorG. B.BurtonJ. P.ReidG. (2013). High throughput sequencing methods and analysis for microbiome research. J. Microbiol. Methods95, 401414. doi: 10.1016/j.mimet.2013.08.011

  • 48

    DonatoJ. J.MoeL. A.ConverseB. J.SmartK. D.BerkleinF. C.McManusP. S.et al. (2010). Metagenomic analysis of apple orchard soil reveals antibiotic resistance genes encoding predicted bifunctional proteins. Appl. Environ. Microbiol.76, 43964401. doi: 10.1128/AEM.01763-09

  • 49

    DouglasG. M.LangilleM. G. (2019). Current and promising approaches to identify horizontal gene transfer events in metagenomes. Genome Biol. Evol.11, 27502766. doi: 10.1093/gbe/evz184

  • 50

    DuY.SunF. (2022). HiCBin: binning metagenomic contigs and recovering metagenome-assembled genomes using hi-C contact maps. Genome Biol.23:63. doi: 10.1186/s13059-022-02626-w

  • 51

    DušanJ.JasperL.AnnaZ.OttoC. (2016). Antimicrobial resistance—a threat to the world’s sustainable development. Ups. J. Med. Sci.121, 159164. doi: 10.1080/03009734.2016.1195900

  • 52

    EckstromK.BarlowJ. W.ChangY.-F. (2019). Resistome metagenomics from plate to farm: the resistome and microbial composition during food waste feeding and composting on a Vermont poultry farm. PLoS One14:e0219807-e. doi: 10.1371/journal.pone.0219807

  • 53

    EmamalipourM.SeidiK.Zununi VahedS.Jahanban-EsfahlanA.JaymandM.MajdiH.et al. (2020). Horizontal gene transfer: from evolutionary flexibility to disease progression. Frontiers in Cell and Developmental Biology.8:229. doi: 10.3389/fcell.2020.00229

  • 54

    Escobar-ZepedaA.Vera-Ponce de LeónA.Sanchez-FloresA. (2015). The road to metagenomics: from microbiology to DNA sequencing technologies and bioinformatics. Front. Genet.6:348. doi: 10.3389/fgene.2015.00348

  • 55

    EssackS. Y.SartoriusB. (2018). Global antibiotic resistance: of contagion, confounders, and the COM-B model. The Lancet Planetary Health.2, e376e377. doi: 10.1016/S2542-5196(18)30187-6

  • 56

    FairR. J.TorY. (2014). Antibiotics and bacterial resistance in the 21st century. Perspectives in medicinal chemistry.6:PMCS14459. doi: 10.4137/PMC.S14459

  • 57

    ForbesJ. D.KnoxN. C.RonholmJ.PagottoF.ReimerA. (2017). Metagenomics: the next culture-independent game changer. Front. Microbiol.8:1069. doi: 10.3389/fmicb.2017.01069

  • 58

    ForsbergK. J.ReyesA.WangB.SelleckE. M.SommerM. O. A.DantasG. (2012). The shared antibiotic Resistome of soil bacteria and human pathogens. Science337, 11071111. doi: 10.1126/science.1220761

  • 59

    FounouL. L.FounouR. C.EssackS. Y. (2016). Antibiotic resistance in the food chain: a developing country-perspective. Front. Microbiol.7:1881. doi: 10.3389/fmicb.2016.01881

  • 60

    FresiaP.AnteloV.SalazarC.GiménezM.D’AlessandroB.AfshinnekooE.et al. (2019). Urban metagenomics uncover antibiotic resistance reservoirs in coastal beach and sewage waters. Microbiome.7, 19. doi: 10.1186/s40168-019-0648-z

  • 61

    FriedmanN. D.TemkinE.CarmeliY. (2016). The negative impact of antibiotic resistance. Clin. Microbiol. Infect.22, 416422. doi: 10.1016/j.cmi.2015.12.002

  • 62

    GandraS.BarterD.LaxminarayanR. (2014). Economic burden of antibiotic resistance: how much do we really know?Clin. Microbiol. Infect.20, 973980. doi: 10.1111/1469-0691.12798

  • 63

    GaoC. (2021). Genome engineering for crop improvement and future agriculture. Cells184, 16211635. doi: 10.1016/j.cell.2021.01.005

  • 64

    GarmendiaL.HernandezA.SanchezM.MartinezJ. (2012). Metagenomics and antibiotics. Clin. Microbiol. Infect.18, 2731. doi: 10.1111/j.1469-0691.2012.03868.x

  • 65

    GarnerE.ChenC.XiaK.BowersJ.EngelthalerD. M.McLainJ.et al. (2018). Metagenomic characterization of antibiotic resistance genes in full-scale reclaimed water distribution systems and corresponding potable systems. Environ. Sci. Technol.52, 61136125. doi: 10.1021/acs.est.7b05419

  • 66

    GaticaJ.CytrynE. (2013). Impact of treated wastewater irrigation on antibiotic resistance in the soil microbiome. Environ. Sci.20, 35293538. doi: 10.1007/s11356-013-1505-4

  • 67

    GeorgeS.PankhurstL.HubbardA.VotintsevaA.StoesserN.SheppardA. E.et al. (2017). Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches. Microbial genomics.3:118. doi: 10.1099/mgen.0.000118

  • 68

    GhuryeJ. S.Cepeda-EspinozaV.PopM. (2016). Focus: microbiome: metagenomic assembly: overview, challenges and applications. Yale J. Biol. Med.89:353. PMID: . PMID: .

  • 69

    GibsonM. K.ForsbergK. J.DantasG. (2015). Improved annotation of antibiotic resistance determinants reveals microbial resistomes cluster by ecology. ISME J.9, 207216. doi: 10.1038/ismej.2014.106

  • 70

    GiganteV.SatiH.BeyerP. (2022). Recent advances and challenges in antibacterial drug development. ADMET and DMPK.10, 147151. doi: 10.5599/admet.1271

  • 71

    GlassE. M.WilkeningJ.WilkeA.AntonopoulosD.MeyerF. (2010). Using the metagenomics RAST server (MG-RAST) for analyzing shotgun metagenomes. Cold Spring Harbor Protocols. 2010:pdb prot5368.

  • 72

    GrenniP.AnconaV.BarraC. A. (2018). Ecological effects of antibiotics on natural ecosystems: a review. Microchem. J.136, 2539. doi: 10.1016/j.microc.2017.02.006

  • 73

    GuoJ.LiJ.ChenH.BondP.YuanZ. (2017). Metagenomic analysis reveals wastewater treatment plants as hotspots of antibiotic resistance genes and mobile genetic elements. Water Res.123, 468478. doi: 10.1016/j.watres.2017.07.002

  • 74

    GuptaA.GuptaB. (2021). Silent slow pandemic of antimicrobial resistance, vol. 8, 354355. doi: 10.18231/j.ijca.2021.067

  • 75

    HallJ. P.BrockhurstM. A.HarrisonE. (2017). Sampling the mobile gene pool: innovation via horizontal gene transfer in bacteria. Philosophical Transactions of the Royal Society B: Biological Sciences.372:20160424. doi: 10.1098/rstb.2016.0424

  • 76

    HallJ. P.HarrisonE.BaltrusD. A. (2022). Introduction: the secret lives of microbial mobile genetic elements. The Royal Society. 377:20200460. doi: 10.1098/rstb.2020.0460

  • 77

    HallR. J.WhelanF. J.McInerneyJ. O.OuY.Domingo-SananesM. R. (2020). Horizontal gene transfer as a source of conflict and cooperation in prokaryotes. Front. Microbiol.11:1569. doi: 10.3389/fmicb.2020.01569

  • 78

    HameedS.FatimaZ. (2021). Integrated Omics approaches to infectious diseases. Berlin: Springer.

  • 79

    HaoH.ChengG.IqbalZ.AiX.HussainH. I.HuangL.et al. (2014). Benefits and risks of antimicrobial use in food-producing animals. Front. Microbiol.5:288. doi: 10.3389/fmicb.2014.00288

  • 80

    HendriksenR. S.BortolaiaV.TateH.TysonG. H.AarestrupF. M.McDermottP. F. (2019). Using genomics to track global antimicrobial resistance. Frontiers. Public Health7, 117. doi: 10.3389/fpubh.2019.00242

  • 81

    HirniesenK.SharmaM.KnielK. (2012). Human enteric pathogen internalization by root uptake into food crops. Foodborne Pathog. Dis.9, 396405. doi: 10.1089/fpd.2011.1044

  • 82

    HolmesA. H.MooreL. S. P.SundsfjordA.SteinbakkM.RegmiS.KarkeyA.et al. (2016). Understanding the mechanisms and drivers of antimicrobial resistance. Lancet387, 176187. doi: 10.1016/S0140-6736(15)00473-0

  • 83

    HuX.ZhouQ.LuoY. (2010). Occurrence and source analysis of typical veterinary antibiotics in manure, soil, vegetables and groundwater from organic vegetable bases, northern China. Environ. Pollut.158, 29922998. doi: 10.1016/j.envpol.2010.05.023

  • 84

    HuntM.MatherA. E.Sánchez-BusóL.PageA. J.ParkhillJ.KeaneJ. A.et al. (2017). ARIBA: rapid antimicrobial resistance genotyping directly from sequencing reads. Microbial genomics.3:131. doi: 10.1099/mgen.0.000131

  • 85

    HusonD. H.AuchA. F.QiJ.SchusterS. C. (2007). MEGAN analysis of metagenomic data. Genome Res.17, 377386. doi: 10.1101/gr.5969107

  • 86

    HuttonG.ChaseC. (2017). Water supply, sanitation, and hygiene. 171198.

  • 87

    HuttonG.HallerL.BartramJ. (2007). Global cost-benefit analysis of water supply and sanitation interventions. J. Water Health5, 481502. doi: 10.2166/wh.2007.009

  • 88

    ImchenM.MoopantakathJ.KumavathR.BarhD.TiwariS.GhoshP.et al. (2020). Current trends in experimental and computational approaches to combat antimicrobial resistance. Front. Genet.11:563975. doi: 10.3389/fgene.2020.563975

  • 89

    ImelfortM.ParksD.WoodcroftB. J.DennisP.HugenholtzP.TysonG. W. (2014). GroopM: an automated tool for the recovery of population genomes from related metagenomes. PeerJ.2:e603. doi: 10.7717/peerj.603

  • 90

    InouyeM.DashnowH.RavenL.-A.SchultzM. B.PopeB. J.TomitaT.et al. (2014). SRST2: rapid genomic surveillance for public health and hospital microbiology labs. Genome Med.6:90. doi: 10.1186/s13073-014-0090-6

  • 91

    JayalakshmiK.ParamasivamM.SasikalaM.TamilamT.SumithraA. (2017). Review on antibiotic residues in animal products and its impact on environments and human health. J Entomol Zool Stud.5, 14461451.

  • 92

    JiaB.RaphenyaA. R.AlcockB.WaglechnerN.GuoP.TsangK. K.et al. (2016). CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res.45:gkw1004. doi: 10.1093/nar/gkw1004

  • 93

    JjembaP. K.WeinrichL. A.ChengW.GiraldoE.LeChevallierM. W. (2010). Regrowth of potential opportunistic pathogens and algae in reclaimed-water distribution systems. Appl. Environ. Microbiol.76, 41694178. doi: 10.1128/AEM.03147-09

  • 94

    JohanssonM. H.BortolaiaV.TansirichaiyaS.AarestrupF. M.RobertsA. P.PetersenT. N. (2021). Detection of mobile genetic elements associated with antibiotic resistance in salmonella enterica using a newly developed web tool: MobileElementFinder. J. Antimicrob. Chemother.76, 101109. doi: 10.1093/jac/dkaa390

  • 95

    JohnsonM.ZaretskayaI.RaytselisY.MerezhukY.McGinnisS.MaddenT. L. (2008). NCBI BLAST: a better web interface. Nucleic Acids Res.36, W5W9. doi: 10.1093/nar/gkn201

  • 96

    JungY.JangH.MatthewsK. R. (2014). Effect of the food production chain from farm practices to vegetable processing on outbreak incidence. Microb. Biotechnol.7, 517527. doi: 10.1111/1751-7915.12178

  • 97

    KangD. D.LiF.KirtonE.ThomasA.EganR.AnH.et al. (2019). MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ.7:e7359. doi: 10.7717/peerj.7359

  • 98

    KangY.LiQ.YinZ.ShenM.ZhaoH.BaiY.et al. (2018). High diversity and abundance of cultivable tetracycline-resistant bacteria in soil following pig manure application. Sci. Rep.8:1489. doi: 10.1038/s41598-018-20050-8

  • 99

    KeenP. L.FugèreR. (2017). Antimicrobial resistance in wastewater treatment processes. New York: John Wiley & Sons.

  • 100

    KhodadadiE.ZeinalzadehE.TaghizadehS.MehramouzB.KamounahF. S.KhodadadiE.et al. (2020). Proteomic applications in antimicrobial resistance and clinical microbiology studies. Infection and Drug Resistance.13:1785. doi: 10.2147/IDR.S238446

  • 101

    KraemerS. A.RamachandranA.PerronG. G. (2019). Antibiotic pollution in the environment: from microbial ecology to public policy. Microorganisms.7:180. doi: 10.3390/microorganisms7060180

  • 102

    KumarC. S. (2017). Molecular mechanism of drug resistance: common themes. Drug Resistance in Bacteria, Fungi, Malaria, and Cancer., 2546. doi: 10.1007/978-3-319-48683-3_2

  • 103

    KumarS. B.ArnipalliS. R.ZiouzenkovaO. (2020). Antibiotics in food chain: the consequences for antibiotic resistance. Antibiotics.9:688. doi: 10.3390/antibiotics9100688

  • 104

    LagierJ. C.ArmougomF.MillionM.HugonP.PagnierI.RobertC.et al. (2012). Microbial culturomics: paradigm shift in the human gut microbiome study. Clin. Microbiol. Infect.18, 11851193. doi: 10.1111/1469-0691.12023

  • 105

    LakinS. M.DeanC.NoyesN. R.DettenwangerA.RossA. S.DosterE.et al. (2017). MEGARes: an antimicrobial resistance database for high throughput sequencing. Nucleic Acids Res.45, D574D580. doi: 10.1093/nar/gkw1009

  • 106

    LakinS. M.KuhnleA.AlipanahiB.NoyesN. R.DeanC.MuggliM.et al. (2019). Hierarchical hidden Markov models enable accurate and diverse detection of antimicrobial resistance sequences. Communications Biology.2:294. doi: 10.1038/s42003-019-0545-9

  • 107

    Lal GuptaC.Kumar TiwariR.CytrynE. (2020). Platforms for elucidating antibiotic resistance in single genomes and complex metagenomes. Environ. Int.138:105667. doi: 10.1016/j.envint.2020.105667

  • 108

    LandersT. F.CohenB.WittumT. E.LarsonE. L. (2012). A review of antibiotic use in food animals: perspective, policy, and potential. Public Health Rep.127, 422. doi: 10.1177/003335491212700103

  • 109

    LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods9, 357359. doi: 10.1038/nmeth.1923

  • 110

    LeplaeR.HebrantA.WodakS. J.ToussaintA. (2004). ACLAME: a CLAssification of Mobile genetic elements. Nucleic Acids Res.32, D45D49. doi: 10.1093/nar/gkh084

  • 111

    LeplaeR.Lima-MendezG.ToussaintA. (2010). ACLAME: a CLAssification of Mobile genetic elements, update 2010. Nucleic Acids Res.38, D57D61. doi: 10.1093/nar/gkp938

  • 112

    LeungH. C.YiuS.-M.YangB.PengY.WangY.LiuZ.et al. (2011). A robust and accurate binning algorithm for metagenomic sequences with arbitrary species abundance ratio. Bioinformatics27, 14891495. doi: 10.1093/bioinformatics/btr186

  • 113

    LiH. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics34, 30943100. doi: 10.1093/bioinformatics/bty191

  • 114

    LiH.DurbinR. (2009). Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics25, 17541760. doi: 10.1093/bioinformatics/btp324

  • 115

    LiD.LiuC.-M.LuoR.SadakaneK.LamT.-W. (2015). MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics31, 16741676. doi: 10.1093/bioinformatics/btv033

  • 116

    LiX.XieY.LiuM.TaiC.SunJ.DengZ.et al. (2018). oriTfinder: a web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements. Nucleic Acids Res.46, W229W234. doi: 10.1093/nar/gky352

  • 117

    LinH.SunW.ZhangZ.ChapmanS. J.FreitagT. E.FuJ.et al. (2016). Effects of manure and mineral fertilization strategies on soil antibiotic resistance gene levels and microbial community in a paddy–upland rotation system. Environ. Pollut.211, 332337. doi: 10.1016/j.envpol.2016.01.007

  • 118

    LittmannJ.BuyxA.CarsO. (2015). Antibiotic resistance: an ethical challenge. Int. J. Antimicrob. Agents46, 359361. doi: 10.1016/j.ijantimicag.2015.06.010

  • 119

    LiuM.LiX.XieY.BiD.SunJ.LiJ.et al. (2019). ICEberg 2.0: an updated database of bacterial integrative and conjugative elements. Nucleic Acids Res.47, D660D665. doi: 10.1093/nar/gky1123

  • 120

    LiuL.WangY.CheY.ChenY.XiaY.LuoR.et al. (2020). High-quality bacterial genomes of a partial-nitritation/anammox system by an iterative hybrid assembly method. Microbiome.8:155. doi: 10.1186/s40168-020-00937-3

  • 121

    LuoR.LiuB.XieY.LiZ.HuangW.YuanJ.et al. (2012). SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience.1:2047. doi: 10.1186/2047-217X-1-18

  • 122

    MacLeanR. C.SanM. A. (2019). The evolution of antibiotic resistance. Science365, 10821083. doi: 10.1126/science.aax3879

  • 123

    MadecJ.-Y.HaenniM. (2018). Antimicrobial resistance plasmid reservoir in food and food-producing animals. Plasmid99, 7281. doi: 10.1016/j.plasmid.2018.09.001

  • 124

    MagourasI.CarmoL. P.StärkK. D.Schüpbach-RegulaG. (2017). Antimicrobial usage and-resistance in livestock: where should we focus?Frontiers in veterinary science.4:148. doi: 10.3389/fvets.2017.00148

  • 125

    MaguireF.JiaB.GrayK. L.LauW. Y. V.BeikoR. G.BrinkmanF. S. (2020). Metagenome-assembled genome binning methods with short reads disproportionately fail for plasmids and genomic islands. Microbial genomics.6:436. doi: 10.1099/mgen.0.000436

  • 126

    MahmoodA. R.Al-HaideriH. H.HassanF. M. (2019). Detection of antibiotics in drinking water treatment plants in Baghdad City, Iraq. Advances. Public Health2019, 111. doi: 10.1155/2019/7851354

  • 127

    MaierT. V.LucioM.LeeL. H.VerBerkmoesN. C.BrislawnC. J.BernhardtJ.et al. (2017). Impact of dietary resistant starch on the human gut microbiome, metaproteome, and metabolome. MBio8, e01343e01317. doi: 10.1128/mBio.01343-17

  • 128

    Malagón-RojasJ.LagosL. (2020). From environment to clinic: the role of pesticides in antimicrobial resistance. Revista Panamericana de Salud Publica= Pan. Am. J. Public Health44:e44-e. doi: 10.26633/RPSP.2020.44

  • 129

    MarkowitzV. M.MavromatisK.IvanovaN. N.ChenI.-M. A.ChuK.KyrpidesN. C. (2009). IMG ER: a system for microbial genome annotation expert review and curation. Bioinformatics25, 22712278. doi: 10.1093/bioinformatics/btp393

  • 130

    MarshallB. M.LevyS. B. (2011). Food animals and antimicrobials: impacts on human health. Clin. Microbiol. Rev.24, 718733. doi: 10.1128/CMR.00002-11

  • 131

    MartiE.VariatzaE.BalcazarJ. L. (2014). The role of aquatic ecosystems as reservoirs of antibiotic resistance. Trends Microbiol.22, 3641. doi: 10.1016/j.tim.2013.11.001

  • 132

    MartínezJ. L. (2008). Antibiotics and antibiotic resistance genes in natural environments. Science321, 365367. doi: 10.1126/science.1159483

  • 133

    McArthurA. G.TsangK. K. (2017). Antimicrobial resistance surveillance in the genomic age. Ann. N. Y. Acad. Sci.1388, 7891. doi: 10.1111/nyas.13289

  • 134

    McEwenS. A.CollignonP. J. (2018). Antimicrobial resistance: a one health perspective. Antimicrobial resistance in bacteria from livestock and companion animals.6, 521547. doi: 10.1128/microbiolspec.ARBA-0009-2

  • 135

    McHardyA. C.MartínH. G.TsirigosA.HugenholtzP.RigoutsosI. (2007). Accurate phylogenetic classification of variable-length DNA fragments. Nat. Methods4, 6372. doi: 10.1038/nmeth976

  • 136

    MeenaV.DotaniyaM.SahaJ.PatraA. (2015). Antibiotics and antibiotic resistant bacteria in wastewater: impact on environment, soil microbial activity and human health. Afr. J. Microbiol. Res.9, 965978. doi: 10.5897/AJMR2015.7195

  • 137

    MichaelI.RizzoL.McArdellC.ManaiaC.MerlinC.SchwartzT.et al. (2013). Urban wastewater treatment plants as hotspots for the release of antibiotics in the environment: a review. Water Res.47, 957995. doi: 10.1016/j.watres.2012.11.027

  • 138

    MiłobedzkaA.FerreiraC.Vaz-MoreiraI.Calderón-FrancoD.GoreckiA.PurkrtovaS.et al. (2021). Monitoring antibiotic resistance genes in wastewater environments: the challenges of filling a gap in the one-health cycle. J. Hazard. Mater.424:127407. doi: 10.1016/j.jhazmat.2021.127407

  • 139

    MindenV.DeloyA.VolkertA. M.LeonhardtS. D.PufalG. (2017). Antibiotics impact plant traits, even at small concentrations. AoB Plants.9:010. doi: 10.1093/aobpla/plx010

  • 140

    MirS. A.ShahM. A.MirM. M.DarB. N.GreinerR.RoohinejadS. (2018). Microbiological contamination of ready-to-eat vegetable salads in developing countries and potential solutions in the supply chain to control microbial pathogens. Food Control85, 235244. doi: 10.1016/j.foodcont.2017.10.006

  • 141

    MirzaS. A.LiaqatI.AwanM. U. F.AfzaalM. (2020). Long-range transport of antibiotics and AMR/ARGs. Antibiotics and Antimicrobial Resistance Genes in the Environment: Elsevier1, 117125. doi: 10.1016/B978-0-12-818882-8.00007-3

  • 142

    MiyataM.IharaI.YoshidG.ToyodK.UmetsuK. (2011). Electrochemical oxidation of tetracycline antibiotics using a Ti/IrO2 anode for wastewater treatment of animal husbandry. Water Sci. Technol.63, 456461. doi: 10.2166/wst.2011.243

  • 143

    MoralezJ.SzenkielK.HamiltonK.PrudenA.LopatkinA. J. (2021). Quantitative analysis of horizontal gene transfer in complex systems. Curr. Opin. Microbiol.62, 103109. doi: 10.1016/j.mib.2021.05.001

  • 144

    MouraA.SoaresM.PereiraC.LeitãoN.HenriquesI.CorreiaA. (2009). INTEGRALL: a database and search engine for integrons, integrases and gene cassettes. Bioinformatics25, 10961098. doi: 10.1093/bioinformatics/btp105

  • 145

    MunitaJ. M.AriasC. A. (2016). Mechanisms of antibiotic resistance. Virulence mechanisms of bacterial pathogens.4, 481511. doi: 10.1128/microbiolspec.VMBF-0016-2015

  • 146

    NaasT.OueslatiS.BonninR. A.DabosM. L.ZavalaA.DortetL.et al. (2017). Beta-lactamase database (BLDB)–structure and function. J. Enzyme Inhib. Med. Chem.32, 917919. doi: 10.1080/14756366.2017.1344235

  • 147

    NamikiT.HachiyaT.TanakaH.SakakibaraY. (2012). MetaVelvet: an extension of velvet assembler to de novo metagenome assembly from short sequence reads. Nucleic Acids Res.40:e155-e. doi: 10.1093/nar/gks678

  • 148

    NijsinghN.MuntheC.LarssonD. J. (2019). Managing pollution from antibiotics manufacturing: charting actors, incentives and disincentives. Environ. Health18, 117. doi: 10.1186/s12940-019-0531-1

  • 149

    NnadozieC. F.OdumeO. N. (2019). Freshwater environments as reservoirs of antibiotic resistant bacteria and their role in the dissemination of antibiotic resistance genes. Environ. Pollut.254:113067. doi: 10.1016/j.envpol.2019.113067

  • 150

    NoguchiH.TaniguchiT.ItohT. (2008). MetaGeneAnnotator: detecting species-specific patterns of ribosomal binding site for precise gene prediction in anonymous prokaryotic and phage genomes. DNA Res.15, 387396. doi: 10.1093/dnares/dsn027

  • 151

    NowrotekM.JałowieckiŁ.HarniszM.PłazaG. A. (2019). Culturomics and metagenomics: in understanding of environmental resistome. Front. Environ. Sci. Eng.13:40. doi: 10.1007/s11783-019-1121-8

  • 152

    NurkS.MeleshkoD.KorobeynikovA.PevznerP. A. (2017). metaSPAdes: a new versatile metagenomic assembler. Genome Res.27, 824834. doi: 10.1101/gr.213959.116

  • 153

    OliveiraN.GonçalvesB.LeeS.OliveiraC.CorassinC. (2020). Use of antibiotics in animal production and its impact on human health. J Food Chem Nanotechnol.6, 4047. doi: 10.17756/jfcn.2020-082

  • 154

    OniciucE. A.LikotrafitiE.Alvarez-MolinaA.PrietoM.SantosJ. A.Alvarez-OrdóñezA. (2018). The present and future of whole genome sequencing (WGS) and whole metagenome sequencing (WMS) for surveillance of antimicrobial resistant microorganisms and antimicrobial resistance genes across the food chain. Genes.9:268. doi: 10.3390/genes9050268

  • 155

    OpatowskiL.OpatowskiM.VongS.TemimeL. (2021). A one-health quantitative model to assess the risk of antibiotic resistance Acquisition in Asian Populations: impact of exposure through food, water. Livestock and Humans. Risk Analysis.41, 14271446. doi: 10.1111/risa.13618

  • 156

    OuH.-Y.HeX.HarrisonE. M.KulasekaraB. R.ThaniA. B.KadiogluA.et al. (2007). MobilomeFINDER: web-based tools for in silico and experimental discovery of bacterial genomic islands. Nucleic Acids Res.35, W97W104. doi: 10.1093/nar/gkm380

  • 157

    Pallares-VegaR.BlaakH.van der PlaatsR.de Roda HusmanA. M.Hernandez LealL.van LoosdrechtM. C. M.et al. (2019). Determinants of presence and removal of antibiotic resistance genes during WWTP treatment: a cross-sectional study. Water Res.161, 319328. doi: 10.1016/j.watres.2019.05.100

  • 158

    PanM.ChuL. (2017). Fate of antibiotics in soil and their uptake by edible crops. Sci. Total Environ.599, 599600.

  • 159

    PapavarnavasN. S.MendelsonM. (2022). Optimising blood cultures: the interplay between diagnostic and antimicrobial stewardship. S. Afr. Med. J.112, 395396. doi: 10.7196/SAMJ.2022.v112i6.16557

  • 160

    ParksD. H.ImelfortM.SkennertonC. T.HugenholtzP.TysonG. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res.25:1043. doi: 10.1101/gr.186072.114

  • 161

    PartridgeS. R.KwongS. M.FirthN.JensenS. O. (2018). Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev.31, 161. doi: 10.1128/CMR.00088-17

  • 162

    PengY.LeungH. C.YiuS.ChinF. Y. (2012). IDBA-UD: A de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics28, 14201428. doi: 10.1093/bioinformatics/bts174

  • 163

    PengZ.MaoY.ZhangN.ZhangL.WangZ.HanM. (2021). Utilizing metagenomic data and Bioinformatic tools for elucidating antibiotic resistance genes in environment. Front. Environ. Sci.9:757365. doi: 10.3389/fenvs.2021.757365

  • 164

    Pérez-SánchezT.Mora-SánchezB.BalcázarJ. L. (2018). Biological approaches for disease control in aquaculture: advantages, limitations and challenges. Trends Microbiol.26, 896903. doi: 10.1016/j.tim.2018.05.002

  • 165

    Pérez-ValeraE.KyselkováM.AhmedE.SladecekF. X. J.GobernaM.ElhottováD. (2019). Native soil microorganisms hinder the soil enrichment with antibiotic resistance genes following manure applications. Sci. Rep.9:6760. doi: 10.1016/j.watres.2019.05.100

  • 166

    PokharelS.RautS.AdhikariB. (2019). Tackling antimicrobial resistance in low-income and middle-income countries. BMJ Specialist Journals. e002104. doi: 10.1136/bmjgh-2019-002104

  • 167

    PolianciucS. I.GurzăuA. E.KissB.ŞtefanM. G.LoghinF. (2020). Antibiotics in the environment: causes and consequences. Medicine and pharmacy reports.93:231. doi: 10.15386/mpr-1742

  • 168

    PopowskaM.RzeczyckaM.MiernikA.Krawczyk-BalskaA.WalshF.DuffyB. (2012). Influence of soil use on prevalence of tetracycline, streptomycin, and erythromycin resistance and associated resistance genes. Antimicrob. Agents Chemother.56, 14341443. doi: 10.1128/AAC.05766-11

  • 169

    PornsukaromS.ThakurS. (2017). Horizontal dissemination of antimicrobial resistance determinants in multiple salmonella serotypes following isolation from the commercial swine operation environment after manure application. Appl. Environ. Microbiol.83, e01503e01517. doi: 10.1128/AEM.01503-17

  • 170

    PowerM. (2019). 17 Antimicrobial resistance. In Current Therapy in Medicine of Australian Mammals. (2nd ed.). eds. VogelnestL.PortasT. (Australia: CSIRO Publishing), 285292.

  • 171

    PreenaP. G.SwaminathanT. R.KumarV. J. R.SinghI. S. B. (2020b). Antimicrobial resistance in aquaculture: a crisis for concern. Biologia75, 14971517. doi: 10.2478/s11756-020-00456-4

  • 172

    PreenaP.SwaminathanT. R.Rejish KumarV.BrightS. I. (2020a). Unravelling the menace: detection of antimicrobial resistance in aquaculture. Lett. Appl. Microbiol.71, 2638. doi: 10.1111/lam.13292

  • 173

    PrestinaciF.PezzottiP.PantostiA. (2015). Antimicrobial resistance: a global multifaceted phenomenon. Pathogens and global health.109, 309318. doi: 10.1179/2047773215Y.0000000030

  • 174

    Prüss-UstünA.WolfJ.BartramJ.ClasenT.CummingO.FreemanM. C.et al. (2019). Burden of disease from inadequate water, sanitation and hygiene for selected adverse health outcomes: an updated analysis with a focus on low-and middle-income countries. Int. J. Hyg. Environ. Health222, 765777. doi: 10.1016/j.ijheh.2019.05.004

  • 175

    QuainooS.CoolenJ. P.van HijumS. A.HuynenM. A.MelchersW. J.van SchaikW.et al. (2017). Whole-genome sequencing of bacterial pathogens: the future of nosocomial outbreak analysis. Clin. Microbiol. Rev.30, 10151063. doi: 10.1128/CMR.00016-17

  • 176

    RajwarA.SrivastavaP.SahgalM. (2016). Microbiology of fresh produce: route of contamination, detection methods, and remedy. Crit. Rev. Food Sci. Nutr.56, 23832390. doi: 10.1080/10408398.2013.841119

  • 177

    RamanR. (2017). The impact of genetically modified (GM) crops in modern agriculture: a review. GM Crops & Food.8, 195208. doi: 10.1080/21645698.2017.1413522

  • 178

    RamtahalM. A.AmoakoD. G.AkebeA. L.SomboroA. M.BesterL. A.EssackS. Y. (2022). A public health insight into salmonella in poultry in Africa: a review of the past decade: 2010–2020. Microb. Drug Resist.28, 710733. doi: 10.1089/mdr.2021.0384

  • 179

    RanjanR.RaniA.MetwallyA.McGeeH. S.PerkinsD. L. (2016). Analysis of the microbiome: advantages of whole genome shotgun versus 16S amplicon sequencing. Biochem. Biophys. Res. Commun.469, 967977. doi: 10.1016/j.bbrc.2015.12.083

  • 180

    RankinD.RochaE.BrownS. (2011). What traits are carried on mobile genetic elements, and why?Heredity106, 110. doi: 10.1038/hdy.2010.24

  • 181

    RedfieldR. J.SoucyS. M. (2018). Evolution of bacterial gene transfer agents. Front. Microbiol.9:2527. doi: 10.3389/fmicb.2018.02527

  • 182

    ReverterM.SarterS.CarusoD.AvarreJ.-C.CombeM.PepeyE.et al. (2020). Aquaculture at the crossroads of global warming and antimicrobial resistance. Nat. Commun.11, 18. doi: 10.1038/s41467-020-15735-6

  • 183

    RiadiG.Medina-MoenneC.HolmesD. S. (2012). TnpPred: a web service for the robust prediction of prokaryotic transposases. Comparative and functional genomics.2012:761. doi: 10.1155/2012/678761

  • 184

    RobinsonT. P.BuD.Carrique-MasJ.FèvreE. M.GilbertM.GraceD.et al. (2016). Antibiotic resistance is the quintessential one health issue. Trans. R. Soc. Trop. Med. Hyg.110, 377380. doi: 10.1093/trstmh/trw048

  • 185

    RobinsonT. P.WertheimH. F. L.KakkarM.KariukiS.BuD.PriceL. B. (2016). Animal production and antimicrobial resistance in the clinic. Lancet387, e1e3. doi: 10.1016/S0140-6736(15)00730-8

  • 186

    RochaC. P.CabralH. N.MarquesJ. C.GonçalvesA. M. (2022). A global overview of aquaculture food production with a focus on the Activity’s development in transitional systems—the case study of a south European country (Portugal). Journal of Marine Science and Engineering.10:417. doi: 10.3390/jmse10030417

  • 187

    Rodríguez-BeltránJ.SørumV.Toll-RieraM.de la VegaC.Peña-MillerR.SanM. Á. (2020). Genetic dominance governs the evolution and spread of mobile genetic elements in bacteria. Proc. Natl. Acad. Sci.117, 1575515762. doi: 10.1073/pnas.2001240117

  • 188

    Rodríguez-MolinaD.MangP.SchmittH.ChifiriucM. C.RadonK.WengenrothL. (2019). Do wastewater treatment plants increase antibiotic resistant bacteria or genes in the environment? Protocol for a systematic review. Syst. Rev.8:304. doi: 10.1186/s13643-019-1236-9

  • 189

    RuppéE.GhozlaneA.TapJ.PonsN.AlvarezA.-S.MaziersN.et al. (2019). Prediction of the intestinal resistome by a three-dimensional structure-based method. Nat. Microbiol.4, 112123. doi: 10.1038/s41564-018-0292-6

  • 190

    SanganyadoE.GwenziW. (2019). Antibiotic resistance in drinking water systems: occurrence, removal, and human health risks. Sci. Total Environ.669, 785797. doi: 10.1016/j.scitotenv.2019.03.162

  • 191

    SarkarD. J.MukherjeeI.ShakilN. A.RanaV. S.KaushikP.DebnathS. (2018). Antibiotics in agriculture: Use and impact. 419.

  • 192

    SayersE. W.BeckJ.BoltonE. E.BourexisD.BristerJ. R.CaneseK.et al. (2021). Database resources of the national center for biotechnology information. Nucleic Acids Res.49:D10. doi: 10.1093/nar/gkaa892

  • 193

    ScharD.ZhaoC.WangY.LarssonD.GilbertM.Van BoeckelT. P. (2021). Twenty-year trends in antimicrobial resistance from aquaculture and fisheries in Asia. Nat. Commun.12, 110. doi: 10.1038/s41467-021-25655-8

  • 194

    SchmiederR.EdwardsR. (2012). Insights into antibiotic resistance through metagenomic approaches. Future Microbiol.7, 7389. doi: 10.2217/fmb.11.135

  • 195

    SeilerC.BerendonkT. U. (2012). Heavy metal driven co-selection of antibiotic resistance in soil and water bodies impacted by agriculture and aquaculture. Front. Microbiol.3:399. doi: 10.3389/fmicb.2012.00399

  • 196

    ShahS. Q.CabelloF. C.L'Abée-LundT. M.TomovaA.GodfreyH. P.BuschmannA. H.et al. (2014). Antimicrobial resistance and antimicrobial resistance genes in marine bacteria from salmon aquaculture and non-aquaculture sites. Environ. Microbiol.16, 13101320. doi: 10.1111/1462-2920.12421

  • 197

    SheY.CaiH.LiuG. (2018). Effects of antibiotic on microflora in ileum and cecum for broilers by 16S rRNA sequence analysis. 16801691.

  • 198

    ShenQ.ZhaoF. (2018). Antibiotics and antibiotic resistance from animal manures to soil: a review. Eur. J. Soil Sci.69, 181195.

  • 199

    ShiP.JiaS.ZhangX.-X.ZhangT.ChengS.LiA. (2013). Metagenomic insights into chlorination effects on microbial antibiotic resistance in drinking water. Water Res.47, 111120. doi: 10.1016/j.watres.2012.09.046

  • 200

    SieberC. M.ProbstA. J.SharrarA.ThomasB. C.HessM.TringeS. G.et al. (2018). Recovery of genomes from metagenomes via a dereplication, aggregation and scoring strategy. Nat. Microbiol.3, 836843. doi: 10.1038/s41564-018-0171-1

  • 201

    SiguierP.PérochonJ.LestradeL.MahillonJ.ChandlerM. (2006). ISfinder: the reference Centre for bacterial insertion sequences. Nucleic Acids Res.34, D32D36. doi: 10.1093/nar/gkj014

  • 202

    SimpsonJ. T.WongK.JackmanS. D.ScheinJ. E.JonesS. J.BirolI. (2009). ABySS: a parallel assembler for short read sequence data. Genome Res.19, 11171123. doi: 10.1101/gr.089532.108

  • 203

    SlizovskiyI. B.MukherjeeK.DeanC. J.BoucherC.NoyesN. R. (2020). Mobilization of antibiotic resistance: are current approaches for colocalizing resistomes and mobilomes useful?Front. Microbiol.11:1376.

  • 204

    SongW.WemheuerB.ZhangS.SteensenK.ThomasT. (2019). MetaCHIP: community-level horizontal gene transfer identification through the combination of best-match and phylogenetic approaches. Microbiome.7, 114. doi: 10.1186/s40168-019-0649-y

  • 205

    SononeS. S.JadhavS.SankhlaM. S.KumarR. (2020). Water contamination by heavy metals and their toxic effect on aquaculture and human health through food chain. Letters in applied NanoBioScience.10, 21482166. doi: 10.33263/LIANBS102.21482166

  • 206

    SpeksnijderD.JaarsmaA.Van Der GugtenA.VerheijT. J.WagenaarJ. (2015). Determinants associated with veterinary antimicrobial prescribing in farm animals in the Netherlands: a qualitative study. Zoonoses Public Health62, 3951. doi: 10.1111/zph.12168

  • 207

    SpielmeyerA. (2018). Occurrence and fate of antibiotics in manure during manure treatments: a short review. Sustain. Chem. Pharm.9, 7686. doi: 10.1016/j.scp.2018.06.004

  • 208

    SrivastavaA.SinghalN.GoelM.VirdiJ. S.KumarM. (2014). CBMAR: a comprehensive β-lactamase molecular annotation resource. Database2014:111. doi: 10.1093/database/bau111

  • 209

    StalderT.PressM. O.SullivanS.LiachkoI.TopE. M. (2019). Linking the resistome and plasmidome to the microbiome. ISME J.13, 24372446. doi: 10.1038/s41396-019-0446-4

  • 210

    StärkK. D.PękalaA.MuellnerP. (2019). Use of molecular and genomic data for disease surveillance in aquaculture: towards improved evidence for decision making. Prev. Vet. Med.167, 190195. doi: 10.1016/j.prevetmed.2018.04.011

  • 211

    SuJ. Q.WeiB.XuC. Y.QiaoM.ZhuY. G. (2014). Functional metagenomic characterization of antibiotic resistance genes in agricultural soils from China. Environ. Int.65, 915. doi: 10.1016/j.envint.2013.12.010

  • 212

    SuenagaH. (2012). Targeted metagenomics: a high-resolution metagenomics approach for specific gene clusters in complex microbial communities. Environ. Microbiol.14, 1322. doi: 10.1111/j.1462-2920.2011.02438.x

  • 213

    SunD.JeannotK.XiaoY.KnappC. W. (2020). Horizontal gene transfer mediated bacterial antibiotic resistanceHorizontal Gene Transfer Mediated Multidrug Resistance: A Global Crisis, Published in Frontiers in microbiology 6.

  • 214

    SundinG. W.WangN. (2018). Antibiotic resistance in plant-pathogenic bacteria. Annu. Rev. Phytopathol.56, 161180. doi: 10.1146/annurev-phyto-080417-045946

  • 215

    TangP.-C. (2016). A metagenomic insight into the role of wastewater treatment plants as potential hotspot for antibiotic resistant bacteria and antibiotic resistance genes. 1378.

  • 216

    TangJ.BuY.ZhangX.-X.HuangK.HeX.YeL.et al. (2016). Metagenomic analysis of bacterial community composition and antibiotic resistance genes in a wastewater treatment plant and its receiving surface water. Ecotoxicol. Environ. Saf.132, 260269. doi: 10.1016/j.ecoenv.2016.06.016

  • 217

    TaşN.de JongA. E. E.LiY.TrublG.XueY.DoveN. C. (2021). Metagenomic tools in microbial ecology research. Curr. Opin. Biotechnol.67, 184191. doi: 10.1016/j.copbio.2021.01.019

  • 218

    TaylorN. G.Verner-JeffreysD. W.Baker-AustinC. (2011). Aquatic systems: maintaining, mixing and mobilising antimicrobial resistance?Trends Ecol. Evol.26, 278284. doi: 10.1016/j.tree.2011.03.004

  • 219

    TeelingH.WaldmannJ.LombardotT.et al. (2004). TETRA: a web-service and a stand-alone program for the analysis and comparison of tetranucleotide usage patterns in DNA sequences. BMC Bioinformatics.5:163. doi: 10.1186/1471-2105-5-163

  • 220

    ThaiQ. K.BösF.PleissJ. (2009). The lactamase engineering database: a critical survey of TEM sequences in public databases. BMC Genomics10, 19. doi: 10.1186/1471-2164-10-390

  • 221

    ThomasT.GilbertJ.MeyerF. (2012). Metagenomics-a guide from sampling to data analysis. Microbial informatics and experimentation.2, 112. doi: 10.1186/2042-5783-2-3

  • 222

    Torres-CortésG.MillánV.Ramírez-SaadH. C.Nisa-MartínezR.ToroN.Martínez-AbarcaF. (2011). Characterization of novel antibiotic resistance genes identified by functional metagenomics on soil samples. Environ. Microbiol.13, 11011114. doi: 10.1111/j.1462-2920.2010.02422.x

  • 223

    TyagiA.SinghB.ThammegowdaN. K. B.SinghN. K. (2019). Shotgun metagenomics offers novel insights into taxonomic compositions, metabolic pathways and antibiotic resistance genes in fish gut microbiome. Arch. Microbiol.201, 295303. doi: 10.1007/s00203-018-1615-y

  • 224

    Van BoeckelT. P.BrowerC.GilbertM.GrenfellB. T.LevinS. A.RobinsonT. P.et al. (2015). Global trends in antimicrobial use in food animals. Proc. Natl. Acad. Sci.112, 56495654. doi: 10.1073/pnas.1503141112

  • 225

    Van BoeckelT. P.GlennonE. E.ChenD.GilbertM.RobinsonT. P.GrenfellB. T.et al. (2017). Reducing antimicrobial use in food animals. Science357, 13501352. doi: 10.1126/science.aao1495

  • 226

    Van CampP.-J.HaslamD. B.PorolloA. (2020). Bioinformatics approaches to the understanding of molecular mechanisms in antimicrobial resistance. Int. J. Mol. Sci.21:1363. doi: 10.3390/ijms21041363

  • 227

    Van GompelL.LuikenR. E. C.HansenR. B.MunkP.BouwknegtM.HeresL.et al. (2020). Description and determinants of the faecal resistome and microbiome of farmers and slaughterhouse workers: a metagenome-wide cross-sectional study. Environ. Int.143:105939. doi: 10.1016/j.envint.2020.105939

  • 228

    Van HoekA. H.MeviusD.GuerraB.MullanyP.RobertsA. P.AartsH. J. (2011). Acquired antibiotic resistance genes: an overview. Front. Microbiol.2:203. doi: 10.3389/fmicb.2011.00203

  • 229

    VaraniA. M.SiguierP.GourbeyreE.CharneauV.ChandlerM. (2011). ISsaga is an ensemble of web-based methods for high throughput identification and semi-automatic annotation of insertion sequences in prokaryotic genomes. Genome Biol.12, 19. doi: 10.1186/gb-2011-12-3-r30

  • 230

    VinayamohanP. G.PellisseryA. J.VenkitanarayananK. (2022). Role of horizontal gene transfer in the dissemination of antimicrobial resistance in food animal production. Curr. Opin. Food Sci.47:100882. doi: 10.1016/j.cofs.2022.100882

  • 231

    von WintersdorffC. J. H.PendersJ.van NiekerkJ. M.MillsN. D.MajumderS.van AlphenL. B.et al. (2016). Dissemination of antimicrobial resistance in microbial ecosystems through horizontal gene transfer. Front. Microbiol.7. doi: 10.3389/fmicb.2016.00173

  • 232

    VoslooS.HuoL.AndersonC. L.DaiZ.SevillanoM.PintoA.et al. (2021). Evaluating de novo assembly and binning strategies for time series drinking water Metagenomes. Microbiology Spectrum.9, e01434e01421. doi: 10.1128/Spectrum.01434-21

  • 233

    WadamoriY.GooneratneR.HussainM. A. (2017). Outbreaks and factors influencing microbiological contamination of fresh produce. J. Sci. Food Agric.97, 13961403. doi: 10.1002/jsfa.8125

  • 234

    WAAFLE. (2022). [Internet]. [cited 22/07/2022].

  • 235

    WallaceJ. C.PortJ. A.SmithM. N.FaustmanE. M. (2017). FARME DB: a functional antibiotic resistance element database. Database2017:165. doi: 10.1093/database/baw165

  • 236

    WangY.LeungH. C. M.YiuS. M.ChinF. Y. L. (2014). MetaCluster-TA: taxonomic annotation for metagenomic data based on assembly-assisted binning. BMC genomics15, 19. doi: 10.1186/1471-2164-15-S1-S12

  • 237

    WangF.-H.QiaoM.ChenZ.SuJ.-Q.ZhuY.-G. (2015). Antibiotic resistance genes in manure-amended soil and vegetables at harvest. J. Hazard. Mater.299, 215221. doi: 10.1016/j.jhazmat.2015.05.028

  • 238

    WattsJ. E.SchreierH. J.LanskaL.HaleM. S. (2017). The rising tide of antimicrobial resistance in aquaculture: sources, sinks and solutions. Mar. Drugs15:158. doi: 10.3390/md15060158

  • 239

    WeeB. A.MuloiD. M.van BunnikB. A. D. (2020). Quantifying the transmission of antimicrobial resistance at the human and livestock interface with genomics. Clin. Microbiol. Infect.26, 16121616. doi: 10.1016/j.cmi.2020.09.019

  • 240

    WilliamsT. J.CavicchioliR. (2014). Marine metaproteomics: deciphering the microbial metabolic food web. Trends Microbiol.22, 248260. doi: 10.1016/j.tim.2014.03.004

  • 241

    WillmsI. M.KamranA.AßmannN. F.KroneD.BolzS. H.FiedlerF.et al. (2019). Discovery of novel antibiotic resistance determinants in Forest and grassland soil Metagenomes. Front. Microbiol.10, 111. doi: 10.3389/fmicb.2019.00460

  • 242

    WilmesP.BondP. L. (2004). The application of two-dimensional polyacrylamide gel electrophoresis and downstream analyses to a mixed community of prokaryotic microorganisms. Environ. Microbiol.6, 911920. doi: 10.1111/j.1462-2920.2004.00687.x

  • 243

    WuM.ScottA. J. (2012). Phylogenomic analysis of bacterial and archaeal sequences with AMPHORA2. Bioinformatics28, 10331034. doi: 10.1093/bioinformatics/bts079

  • 244

    WuY.-W.SimmonsB. A.SingerS. W. (2016). MaxBin 2.0: an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics32, 605607. doi: 10.1093/bioinformatics/btv638

  • 245

    XavierB. B.DasA. J.CochraneG.De GanckS.Kumar-SinghS.AarestrupF. M.et al. (2016). Consolidating and exploring antibiotic resistance gene data resources. J. Clin. Microbiol.54, 851859. doi: 10.1128/JCM.02717-15

  • 246

    XieW. Y.ShenQ.ZhaoF. (2018). Antibiotics and antibiotic resistance from animal manures to soil: a review. Eur. J. Soil Sci.69, 181195. doi: 10.1111/ejss.12494

  • 247

    XuW.ChenT.PeiY.GuoH.LiZ.YangY.et al. (2021). Characterization of shallow whole-Metagenome shotgun sequencing as a high-accuracy and low-cost method by complicated mock microbiomes. Front. Microbiol.12:380. doi: 10.3389/fmicb.2021.820380

  • 248

    YangY.AshworthA. J.WillettC.CookK.UpadhyayA.OwensP. R.et al. (2019). Review of antibiotic resistance, ecology, dissemination, and mitigation in US broiler poultry systems. Front. Microbiol.10:2639. doi: 10.3389/fmicb.2019.02639

  • 249

    YangX.NoyesN. R.DosterE.MartinJ. N.LinkeL. M.MagnusonR. J.et al. (2016). Use of metagenomic shotgun sequencing technology to detect foodborne pathogens within the microbiome of the beef production chain. Appl. Environ. Microbiol.82, 24332443. doi: 10.1128/AEM.00078-16

  • 250

    YaoY.MaddamsettiR.WeissA.HaY.WangT.WangS.et al. (2022). Intra-and interpopulation transposition of mobile genetic elements driven by antibiotic selection. Nature Ecology & Evolution.6, 555564. doi: 10.1038/s41559-022-01705-2

  • 251

    YeeR.SimnerP. J. (2019). Next-generation sequencing approaches to predicting antimicrobial susceptibility testing results. Advances in Molecular Pathology.2, 99110. doi: 10.1016/j.yamp.2019.07.008

  • 252

    YinX.JiangX.-T.ChaiB.LiL.YangY.ColeJ. R.et al. (2018). ARGs-OAP v2. 0 with an expanded SARG database and hidden Markov models for enhancement characterization and quantification of antibiotic resistance genes in environmental metagenomes. Bioinformatics34, 22632270. doi: 10.1093/bioinformatics/bty053

  • 253

    YuanL.WangL.LiZ.-H.ZhangM.-Q.ShaoW.ShengG.-P. (2019). Antibiotic resistance and microbiota in the gut of Chinese four major freshwater carp from retail markets. Environ. Pollut.255:113327. doi: 10.1016/j.envpol.2019.113327

  • 254

    YukgehnaishK.KumarP.SivachandranP.MarimuthuK.ArshadA.ParayB. A.et al. (2020). Gut microbiota metagenomics in aquaculture: factors influencing gut microbiome and its physiological role in fish. Rev. Aquac.3, 19031927. doi: 10.1111/raq.12416

  • 255

    ZankariE.HasmanH.CosentinoS.VestergaardM.RasmussenS.LundO.et al. (2012). Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother.67, 26402644. doi: 10.1093/jac/dks261

  • 256

    ZhangX.FigeysD. (2019). Perspective and guidelines for Metaproteomics in microbiome studies. J. Proteome Res.18:054. doi: 10.1021/acs.jproteome.9b00054

  • 257

    ZhangY.-J.HuH.-W.ChenQ.-L.SinghB. K.YanH.ChenD.et al. (2019). Transfer of antibiotic resistance from manure-amended soils to vegetable microbiomes. Environ. Int.130:104912. doi: 10.1016/j.envint.2019.104912

  • 258

    ZhangX.LiL.MayneJ.NingZ.StintziA.FigeysD. (2018). Assessing the impact of protein extraction methods for human gut metaproteomics. J. Proteome180, 120127. doi: 10.1016/j.jprot.2017.07.001

  • 259

    ZhaoY.SuJ.-Q.AnX.-L.HuangF.-Y.RensingC.BrandtK. K.et al. (2018). Feed additives shift gut microbiota and enrich antibiotic resistance in swine gut. Sci. Total Environ.621, 12241232. doi: 10.1016/j.scitotenv.2017.10.106

  • 260

    ZhaoY.YangQ. E.ZhouX.WangF.-H.MuurinenJ.VirtaM. P.et al. (2020). Antibiotic resistome in the livestock and aquaculture industries: status and solutions. Crit. Rev. Environ. Sci. Technol.1, 21592196. doi: 10.1080/10643389.2020.1777815

  • 261

    ZhuB.ChenQ.ChenS.ZhuY.-G. (2017). Does organically produced lettuce harbor higher abundance of antibiotic resistance genes than conventionally produced?Environ. Int.98, 152159. doi: 10.1016/j.envint.2016.11.001

  • 262

    ZhuD.GilesM.YangX.-R.DaniellT.NeilsonR.ZhuY.-G. (2019). Phyllosphere of staple crops under pig manure fertilization, a reservoir of antibiotic resistance genes. Environ. Pollut.252, 227235. doi: 10.1016/j.envpol.2019.05.098

Summary

Keywords

antibiotic resistance, metagenomics, bioinformatics, microbiomes, gene transfer

Citation

Pillay S, Calderón-Franco D, Urhan A and Abeel T (2022) Metagenomic-based surveillance systems for antibiotic resistance in non-clinical settings. Front. Microbiol. 13:1066995. doi: 10.3389/fmicb.2022.1066995

Received

11 October 2022

Accepted

09 November 2022

Published

02 December 2022

Volume

13 - 2022

Edited by

Muhammad Yasir, King Abdulaziz University, Saudi Arabia

Reviewed by

Lateef Salam, Elizade University, Nigeria; Varsha Rani Gajamer, Darjeeling Government College, India

Updates

Copyright

*Correspondence: Thomas Abeel,

This article was submitted to Evolutionary and Genomic Microbiology, a section of the journal Frontiers in Microbiology

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics