Metagenomic profiling of microbial composition and antibiotic resistance determinants in Puget Sound

Metagenomic profiling of antibiotic resistance genes in Red Sea brine pools

Antibiotic resistance (AR) is an alarming global health concern, causing an annual death rate of more than 35,000 deaths in the US. AR is a natural phenomenon, reported in several pristine environments. In this study, we report AR in pristine Red Sea deep brine pools. Antimicrobial resistance genes (ARGs) were detected for several drug classes with tetracycline and macrolide resistance being the most abundant. As expected, ARGs abundance increased in accordance with the level of human impact with pristine Red Sea samples having the lowest mean ARG level followed by estuary samples, while activated sludge samples showed a significantly higher ARG level. ARG hierarchical clustering grouped drug classes for which resistance was detected in Atlantis II Deep brine pool independent of the rest of the samples. ARG abundance was significantly lower in the Discovery Deep brine pool. A correlation between integrons and ARGs abundance in brine pristine samples could be detected, while insertion sequences and plasmids showed a correlation with ARGs abundance in human-impacted samples not seen in brine pristine samples. This suggests different roles of distinct mobile genetic elements (MGEs) in ARG distribution in pristine versus human-impacted sites. Additionally, we showed the presence of mobile antibiotic resistance genes in the Atlantis II brine pool as evidenced by the co-existence of integrases and plasmid replication proteins on the same contigs harboring predicted multidrug-resistant efflux pumps. This study addresses the role of non-pathogenic environmental bacteria as a silent reservoir for ARGs, and the possible horizontal gene transfer mechanism mediating ARG acquisition.

Profile and actual transmissibility of Carbapenem resistance genes: Intracellular and extracellular DNA in hospital wastewater

The current worldwide spread of carbapenem resistance genes (CRGs) has posed a major public health threat, which continues to grow in severity. Hospital wastewaters (HWWs) are major reservoirs for antibiotic resistance genes, while resistomes in HWWs are still poorly characterized when it comes to CRGs. We comprehensively characterized the profile and actual transmissibility of extracellular CRGs (eCRGs) and intracellular CRGs (iCRGs) in HWWs for the first time. In this study, CRGs showed similar relative abundance in treated and untreated HWWs. Meanwhile, HWWs treatments led to the enrichment of bla_IMP-8, probably attributed to the promotion of Novosphingobium and Prosthecobacter after treatment. To evaluate the transmission potential of CRGs, extracellular and intracellular carbapenem-resistant plasmids were captured from HWWs by transformation and conjugation, respectively. We found an interesting phenomenon regarding the transmission characteristics of CRGs: bla_KPC-carrying plasmids could only be captured by transformation, while bla_NDM-carrying plasmids were captured by conjugation. Further experiments showed that HWW treatments increased the conjugation ability of bla_NDM. In conclusion, our study demonstrated that HWWs are significant reservoirs of CRGs and various CRGs exhibit different modes of transmission in HWWs. CRGs cannot be removed by membrane bioreactor and chlorine disinfection. An urgent need is to develop more efficient wastewater treatments to limit CRG dissemination.

Sediment prokaryotic assembly, methane cycling, and ammonia oxidation potentials in response to increasing antibiotic pollution at shrimp aquafarm

Antibiotic pollution poses serious threats to public health and ecological processes. However, systematic research regarding the interactive effects of increasing nutrient and antibiotic pollutions on the prokaryotic community, particularly taxa that contribute to greenhouse gas emissions, is lacking. By exploring the complex interactions that occur between interkingdom bacteria and archaea, biotic and abiotic factors, the responses of sediment prokaryotic assembly were determined along a significant antibiotic pollution gradient. Bacterial and archaeal communities were primarily governed by sediment antibiotic pollution, ammonia, phosphate, and redox potential, which further affected enzyme activities. The two communities nonlinearly responded to increasing antibiotic pollution, with significant tipping points of 3.906 and 0.979 mg/kg antibiotics, respectively. The combined antibiotic concentration-discriminatory taxa of bacteria and archaea accurately (98.0% accuracy) diagnosed in situ antibiotic concentrations. Co-abundance analysis revealed that the methanogens, methanotrophs, sulfate-reducing bacteria, and novel players synergistically contributed to methane cycling. Antibiotic pollution caused the dominant role of ammonia-oxidizing archaea in ammonia oxidation at these alkaline sediments. Collectively, the significant tipping points and bio-indicators afford indexes for regime shift and quantitative diagnosis of antibiotic pollution, respectively. Antibiotic pollution could expedite methane cycling and mitigate nitrous oxide yield, which are previously unrecognized ecological effects. These findings provide new insights into the interactive biological and ecological consequences of increasing nutrient and antibiotic pollutions.

Metagenomic insights into surface water microbial communities of a South Asian mangrove ecosystem

Estuaries are one of the most productive ecosystems and their productivity is maintained by resident microbial communities. Recent alterations driven by climate change have further escalated these stressors leading to the propagation of traits such as antibiotic resistance and heavy metal resistance in microbial communities. Surface water samples from eleven stations along the Thakuran and Matla estuaries of the Sundarbans Biosphere Reserve (SBR) of Sundarbans mangrove located in South Asia were sampled in monsoon (June) 2019 to elucidate resident microbial communities based on Nanopore sequencing. Metagenomic analyses revealed the widespread dominance of Proteobacteria across all the stations along with a high abundance of Firmicutes. Other phyla, including Euryarchaeota, Thaumarchaeota, Actinobacteria, Bacteroidetes and Cyanobacteria showed site-specific trends in abundance. Further taxonomic affiliations showed Gammaproteobacteria and Alphaproteobacteria to be dominant classes with high abundances of Bacilli in SBR_Stn58 and SBR_Stn113. Among the eukaryotic communities, the most abundant classes included Prasinophyceae, Saccharyomycetes and Sardariomycetes. Functional annotation showed metabolic activities such as carbohydrate, amino acid, nitrogen and phosphorus metabolisms to be uniformly distributed across all the studied stations. Pathways such as stress response, sulphur metabolism and motility-associated genes appeared in low abundances in SBR. Functional traits such as antibiotic resistance showed overwhelming dominance of genes involved in multidrug resistance along with widespread resistance towards commonly used antibiotics including Tetracycline, glycopeptide and aminoglycoside. Metal resistance genes including arsenic, nickel and copper were found in comparable abundances across the studied stations. The prevalence of ARG and MRG might indicate presence of pollutants and hint toward deteriorating ecosystem health status of Sundarbans mangrove.

Distinctive signatures of pathogenic and antibiotic resistant potentials in the hadal microbiome

BACKGROUND

Hadal zone of the deep-sea trenches accommodates microbial life under extreme energy limitations and environmental conditions, such as low temperature, high pressure, and low organic matter down to 11,000 m below sea level. However, microbial pathogenicity, resistance, and adaptation therein remain unknown. Here we used culture-independent metagenomic approaches to explore the virulence and antibiotic resistance in the hadal microbiota of the Mariana Trench.

RESULTS

The results indicate that the 10,898 m Challenger Deep bottom sediment harbored prosperous microbiota with contrasting signatures of virulence factors and antibiotic resistance, compared with the neighboring but shallower 6038 m steep wall site and the more nearshore 5856 m Pacific basin site. Virulence genes including several famous large translocating virulence genes (e.g., botulinum neurotoxins, tetanus neurotoxin, and Clostridium difficile toxins) were uniquely detected in the trench bottom. However, the shallower and more nearshore site sediment had a higher abundance and richer diversity of known antibiotic resistance genes (ARGs), especially for those clinically relevant ones (e.g., fosX, sul1, and TEM-family extended-spectrum beta-lactamases), revealing resistance selection under anthropogenic stresses. Further analysis of mobilome (i.e., the collection of mobile genetic elements, MGEs) suggests horizontal gene transfer mediated by phage and integrase as the major mechanism for the evolution of Mariana Trench sediment bacteria. Notably, contig-level co-occurring and taxonomic analysis shows emerging evidence for substantial co-selection of virulence genes and ARGs in taxonomically diverse bacteria in the hadal sediment, especially for the Challenger Deep bottom where mobilized ARGs and virulence genes are favorably enriched in largely unexplored bacteria.

CONCLUSIONS

This study reports the landscape of virulence factors, antibiotic resistome, and mobilome in the sediment and seawater microbiota residing hadal environment of the deepest ocean bottom on earth. Our work unravels the contrasting and unique features of virulence genes, ARGs, and MGEs in the Mariana Trench bottom, providing new insights into the eco-environmental and biological processes underlying microbial pathogenicity, resistance, and adaptative evolution in the hadal environment.

Microbiome and antibiotic resistance profiling in submarine effluent-receiving coastal waters in Croatia

Wastewater treatment plant (WWTP) effluents are pointed as hotspots for the introduction of both commensal and pathogenic bacteria as well as their antibiotic resistance genes (ARGs) in receiving water bodies. For the first time, the effect of partially treated submarine effluents was explored at the bottom and surface of the water column to provide a comprehensive overview of the structure of the microbiome and associated AR, and to assess environmental factors leading to their alteration. Seawater samples were collected over a 5-month period from submarine outfalls in central Adriatic Sea, Croatia. 16S rRNA amplicon sequencing was used to establish taxonomic and resistome profiles of the bacterial communities. The community differences observed between the two discharge areas, especially in the abundance of Proteobacteria and Firmicutes, could be due to the origin of wastewaters treated in WWTPs and the limiting environmental conditions such as temperature and nutrients. PICRUSt2 analysis inferred the total content of ARGs in the studied microbiomes and showed the highest abundance of resistance genes encoding multidrug efflux pumps, such as MexAB-OprM, AcrEF-TolC and MdtEF-TolC, followed by the modified peptidoglycan precursors, transporter genes encoding tetracycline, macrolide and phenicol resistance, and the bla operon conferring β-lactam resistance. A number of pathogenic genera introduced by effluents, including Acinetobacter, Arcobacter, Bacteroides, Escherichia-Shigella, Klebsiella, Pseudomonas, and Salmonella, were predicted to account for the majority of efflux pump-driven multidrug resistance, while Acinetobacter, Salmonella, Bacteroides and Pseudomonas were also shown to be the predominant carriers of non-efflux ARGs conferring resistance to most of nine antibiotic classes. Taken together, we evidenced the negative impact of submarine discharges of treated effluents via alteration of physico-chemical characteristics of the water column and enrichment of bacterial community with nonindigenous taxa carrying an arsenal of ARGs, which could contribute to the further propagation of the AR in the natural environment.

Metagenomic and Recombination Analyses of Antimicrobial Resistance Genes from Recreational Waters of Black Sea Coastal Areas and Other Marine Environments Unveil Extensive Evidence for Their both Intrageneric and Intergeneric Transmission across Genetically Very Diverse Microbial Communities

Microbial communities of marine coastal recreation waters have become large reservoirs of AMR genes (ARGs), contributing to the emergence and transmission of various zoonotic, foodborne and other infections that exhibit resistance to various antibiotics. Thus, it is highly imperative to determine ARGs assemblages as well as mechanisms and trajectories of their transmission across these microbial communities for our better understanding of the evolutionary trends of AMR (AMR). In this study, using metagenomics approaches, we screened for ARGs in recreation waters of the Black Sea coastal areas of the Batumi City (Georgia). Also, a large array of the recombination detection algorithms of the SplitsTree, RDP4, and GARD was applied to elucidate genetic recombination of ARGs and trajectories of their transmission across various marine microbial communities. The metagenomics analyses of sea water samples, obtained from across the above marine sites, could identify putative ARGs encoding for multidrug resistance efflux transporters mainly from the Major Facilitator and Resistance Nodulation Division superfamilies. The data, generated by SplitsTree (fit ≥95.619; bootstrap values ≥ 95; Phi p ≤ 0.0494), RDP4 (p ≤ 0.0490), and GARD, provided strong statistical evidence not only for intrageneric recombination of these ARGs, but also for their intergeneric recombination across fairly large and diverse microbial communities of marine environment. These bacteria included both human pathogenic and nonpathogenic species, exhibiting collectively the genera of Vibrio, Aeromonas, Synechococcus, Citromicrobium, Rhodobacteraceae, Pseudoalteromonas, Altererythrobacter, Erythrobacter, Altererythrobacter, Marivivens, Xuhuaishuia, and Loktanella. The above nonpathogenic bacteria are strongly suggested to contribute to ARGs transmission in marine ecosystems.

Surveillance for Antibiotic-Resistant E. coli in the Salish Sea Ecosystem

E. coli was isolated from the Salish Sea (Puget Sound) ecosystem, including samples of marine and fresh water, and wildlife dependent on this environment. E. coli isolates were assessed for phenotypic and genotypic resistance to antibiotics. A total of 305 E. coli isolates was characterized from samples collected from: marine water obtained in four quadrants of the Salish Sea; select locations near beaches; fresh water from streams near marine beaches; and fecal samples from harbor porpoises (Phocoena phocoena), harbor seals (Phoca vitulina), river otters (Lontra canadensis), and English sole (Parophrys vetulus). Isolates were evaluated using antimicrobial susceptibility typing, whole-genome sequencing, fumC, and multilocus sequence typing. Resistance and virulence genes were identified from sequence data. Of the 305 isolates from Salish Sea samples, 20 (6.6%) of the E. coli were intermediate, and 31 (10.2%) were resistant to ≥1 class of antibiotics, with 26.9% of nonsusceptible (resistant and intermediate resistant) E. coli isolates from marine mammals and 70% from river otters. The proportion of nonsusceptible isolates from animals was significantly higher than samples taken from marine water (p < 0.0001). A total of 196 unique STs was identified including 37 extraintestinal pathogenic E. coli (ExPEC)-associated STs [ST10, ST38, ST58, ST69, ST73, ST117, ST131, and ST405]. The study suggests that animals may be potential sentinels for antibiotic-resistant and ExPEC E. coli in the Salish Sea ecosystem.

Environmental Spread of Antibiotic Resistance

Antibiotic resistance represents a global health concern. Soil, water, livestock and plant foods are directly or indirectly exposed to antibiotics due to their agricultural use or contamination. This selective pressure has acted synergistically to bacterial competition in nature to breed antibiotic-resistant (AR) bacteria. Research over the past few decades has focused on the emergence of AR pathogens in food products that can cause disease outbreaks and the spread of antibiotic resistance genes (ARGs), but One Health approaches have lately expanded the focus to include commensal bacteria as ARG donors. Despite the attempts of national and international authorities of developed and developing countries to reduce the over-prescription of antibiotics to humans and the use of antibiotics as livestock growth promoters, the selective flow of antibiotic resistance transmission from the environment to the clinic (and vice-versa) is increasing. This review focuses on the mechanisms of ARG transmission and the hotspots of antibiotic contamination resulting in the subsequent emergence of ARGs. It follows the transmission of ARGs from farm to plant and animal food products and provides examples of the impact of ARG flow to clinical settings. Understudied and emerging antibiotic resistance selection determinants, such as heavy metal and biocide contamination, are also discussed here.

The relationship between depression and risk of metabolic syndrome: a meta-analysis of observational studies

INTRODUCTION

The link between metabolic syndrome and depression has always been controversial. Different studies that have examined the relationship between metabolic syndrome and depression have reported different results. Therefore, the goal of the present study was to examine the association between depression and MetS by meta-analysis.

METHODS

Embase, Scopus, PubMed, and ISI were searched for publications in English from January 1990 to February 2020. Search included cohort and cross-sectional studies aimed at examining the association between depression and MetS. The risk of bias was assessed by Newcastle-Ottawa Scale. Heterogeneity and publication bias were tested, subgroup analysis and meta-regression were conducted.

RESULTS

49 studies with total sample size 399,494 were analyzed. Results indicated the odds of MetS was higher in depressed compared to non-depressed individuals [OR: 1.48; 95 %CI: 1.33-1.64) vs. (OR: 1.38; 95 %CI: 1.17-1.64)]. For cross-sectional studies, depressed patients in Europe (OR = 1.35; 95 %CI: 1.47-1.99) were at higher odds of MetS compared to those in America and Asia. For cohort studies, depressed patients in America (OR = 1.46; 95 %CI: 1.16-1.84) were at higher odds of MetS than those in Europe. Cross-sectional studies indicated women with depression were at higher odds of MetS (OR = 1.95; 95 %CI: 1.38-2.74) compared to men. In both types of studies, the odds of MetS decreased with age.

CONCLUSIONS

Metabolic syndrome is more common in depressed compared to non-depressed individuals.

Antibiotic resistome from the One-Health perspective: understanding and controlling antimicrobial resistance transmission

The concept of the antibiotic resistome was introduced just over a decade ago, and since then, active resistome studies have been conducted. In the present study, we describe the previously established concept of the resistome, which encompasses all types of antibiotic resistance genes (ARGs), and the important findings from each One-Health sector considering this concept, thereby emphasizing the significance of the One-Health approach in understanding ARG transmission. Cutting-edge research methodologies are essential for deciphering the complex resistome structure in the microbiomes of humans, animals, and the environment. Based on the recent achievements of resistome studies in multiple One-Health sectors, future directions for resistome research have been suggested to improve the understanding and control of ARG transmission: (1) ranking the critical ARGs and their hosts; (2) understanding ARG transmission at the interfaces of One-Health sectors; (3) identifying selective pressures affecting the emergence, transmission, and evolution of ARGs; and (4) elucidating the mechanisms that allow an organism to overcome taxonomic barriers in ARG transmission.

Structure of Bacterial Community with Resistance to Antibiotics in Aquatic Environments. A Systematic Review

Aquatic environments have been affected by the increase in bacterial resistant to antibiotics. The aim of this review is to describe the studies carried out in relation to the bacterial population structure and antibiotic resistance genes in natural and artificial water systems. We performed a systematic review based on the PRISMA guideline (preferred reporting items for systematic reviews and meta-analyzes). Articles were collected from scientific databases between January 2010 and December 2020. Sixty-eight papers meeting the inclusion criteria, i.e., "reporting the water bacterial community composition", "resistance to antibiotics", and "antibiotic resistance genes (ARG)", were evaluated according to pre-defined validity criteria. The results indicate that the predominant phyla were Firmicutes and Bacteroidetes in natural and artificial water systems. Gram-negative bacteria of the family Enterobacteraceae with resistance to antibiotics are commonly reported in drinking water and in natural water systems. The ARGs mainly reported were those that confer resistance to β-lactam antibiotics, aminoglycosides, fluoroquinolones, macrolides and tetracycline. The high influence of anthropogenic activity in the environment is evidenced. The antibiotic resistance genes that are mainly reported in the urban areas of the world are those that confer resistance to the antibiotics that are most used in clinical practice, which constitutes a problem for human and animal health.

Source tracking of antibiotic resistance genes in the environment - Challenges, progress, and prospects

Antibiotic resistance has become a global public health concern, rendering common infections untreatable. Given the widespread occurrence, increasing attention is being turned toward environmental pathways that potentially contribute to antibiotic resistance gene (ARG) dissemination outside the clinical realm. Studies during the past decade have clearly proved the increased ARG pollution trend along with gradient of anthropogenic interference, mainly through marker-ARG detection by PCR-based approaches. However, accurate source-tracking has been always confounded by various factors in previous studies, such as autochthonous ARG level, spatiotemporal variability and environmental resistome complexity, as well as inherent method limitation. The rapidly developed metagenomics profiles ARG occurrence within the sample-wide genomic context, opening a new avenue for source tracking of environmental ARG pollution. Coupling with machine-learning classification, it has been demonstrated the potential of metagenomic ARG profiles in unambiguously assigning source contribution. Through identifying indicator ARG and recovering ARG-host genomes, metagenomics-based analysis will further increase the resolution and accuracy of source tracking. In this review, challenges and progresses in source-tracking studies on environmental ARG pollution will be discussed, with specific focus on recent metagenomics-guide approaches. We propose an integrative metagenomics-based framework, in which coordinated efforts on experimental design and metagenomic analysis will assist in realizing the ultimate goal of robust source-tracking in environmental ARG pollution.

Whole-genome analysis of extraintestinal pathogenic Escherichia coli (ExPEC) MDR ST73 and ST127 isolated from endangered southern resident killer whales (Orcinus orca)

BACKGROUND

Limited studies have investigated the microbial diversity of wild marine mammals.

OBJECTIVES

This study characterized Escherichia coli isolates collected from fresh faecal samples of endangered southern resident killer whales (Orcinus orca) located by detection dogs.

METHODS

WGS of each strain was done to determine ST (using MLST), clonotype (C:H), antimicrobial resistance and virulence profile. Conjugation experiments were done to determine the mobility of the tet(B) tetracycline resistance gene.

RESULTS

All isolates belonged to extraintestinal pathogenic E. coli (ExPEC) clonal lineages ST73 (8/9) and ST127 (1/9), often associated with human community-acquired urinary tract disease. Clonotyping using fumC and fimH alleles showed divergence in clonal lineages, with ST73 isolates belonging to the C24:H10 clade and the ST127 isolate belonging to C14:H2. The eight ST73 isolates carried multiple acquired antibiotic resistance genes, including aadA1, sul1 and tet(B), encoding aminoglycoside, sulphonamide and tetracycline resistance, respectively. Conjugative transfer of the resistance gene tet(B) was observed for three of the eight isolates. ST127 did not carry any of these acquired resistance genes. Virulence-associated genes identified included those encoding adhesins (iha, papC, sfaS), toxins (sat, vat, pic, hlyA, cnf1), siderophores (iutA, fyuA, iroN, ireA), serum survival/protectins (iss, ompT), capsule (kpsM) and pathogenicity island marker (malX).

CONCLUSIONS

Orca whales can carry antibiotic-resistant potentially pathogenic strains of E. coli. Possible sources include contamination of the whale's environment and/or food. It is unknown whether these isolates cause disease in southern resident killer whales, which could contribute to the ongoing decline of this critically endangered population.

Freshwater environments as reservoirs of antibiotic resistant bacteria and their role in the dissemination of antibiotic resistance genes

Freshwater environments are susceptible to possible contamination by residual antibiotics that are released through different sources, such as agricultural runoffs, sewage discharges and leaching from nearby farms. Freshwater environment can thus become reservoirs where an antibiotic impact microorganisms, and is an important public health concern. Degradation and dilution processes are fundamental for predicting the actual risk of antibiotic resistance dissemination from freshwater reservoirs. This study reviews major approaches for detecting and quantifying antibiotic resistance bacteria (ARBs) and genes (ARGs) in freshwater and their prevalence in these environments. Finally, the role of dilution, degradation, transmission and the persistence and fate of ARB/ARG in these environments are also reviewed. Culture-based single strain approaches and molecular techniques that include polymerase chain reaction (PCR), quantitative polymerase chain reaction (qPCR) and metagenomics are techniques for quantifying ARB and ARGs in freshwater environments. The level of ARBs is extremely high in most of the river systems (up to 98% of the total detected bacteria), followed by lakes (up to 77% of the total detected bacteria), compared to dam, pond, and spring (<1%). Of most concern is the occurrence of extended-spectrum β-lactamase producing Enterobacteriaceae, methicillin resistant Staphylococcus aureus (MRSA) and vancomycin resistant Enterococcus (VRE), which cause highly epidemic infections. Dilution and natural degradation do not completely eradicate ARBs and ARGs in the freshwater environment. Even if the ARBs in freshwater are effectively inactivated by sunlight, their ARG-containing DNA can still be intact and capable of transferring resistance to non-resistant strains. Antibiotic resistance persists and is preserved in freshwater bodies polluted with high concentrations of antibiotics. Direct transmission of indigenous freshwater ARBs to humans as well as their transitory insertion in the microbiota can occur. These findings are disturbing especially for people that rely on freshwater resources for drinking, crop irrigation, and food in form of fish.

Deciphering extracellular antibiotic resistance genes (eARGs) in activated sludge by metagenome

Antibiotic resistance genes (ARGs) including extracellular ARGs (eARGs) and intracellular ARGs (iARGs), are recognized as emerging environmental contaminants. Despite extensive efforts to profile ARGs in their "hotspots" wastewater treatment plants (WWTPs), the contribution of eARGs to antibiotic resistance spread remains unclear. Here, we applied metagenomic sequencing to investigate the distribution, mobility and microbial hosts of eARGs in activated sludge from five WWTPs. The total relative abundance eARGs ranged from 9.5 × 10^-6 to 1.3 × 10^-4, mainly encompassing elfamycin, dual drug, and aminoglycoside resistance genes. Multiple eARGs (e.g., EF-Tu, ropB, and rpsL mutants) were shared among the five WWTPs, and some clustered in the same genetic element (e.g., EF-Tu-rpsJ). eARGs were found to frequently co-localize with the eMGEs (e.g., sul1-3'CS-TnAs3, sul2-intI1-ISVsa3, and tetX-p63039), which may facilitate the mobilization of eARGs. Most eARGs likely originated from the genera Mycobacterium (6.7%), Nitrosomonas (5.3%), Steroidobacter (5.3%), Nitrospira (5.2%) and Pseudomonas (5.1%). No significant difference in the diversity, abundance, and mobility was observed between eARGs and iARGs. The host composition of eARGs and iARGs from municipal WWTPs are consistently dominated by Nitrosomonas, Steroidobacter, Nitrospira and Pseudomonas, while some differentially enriched genera (especially Nitrospira) in the hosts of iARGs compared with those of eARGs from the swine WWTP were identified. Our findings corroborate the mobile eARGs reservoir in WWTPs, thereby laying foundation for mitigating widespread antibiotic resistance.

Understanding antimicrobial discovery and resistance from a metagenomic and metatranscriptomic perspective: advances and applications

Our inability to cultivate most microorganisms, specifically bacteria, in the laboratory has for many years restricted our view and understanding of the bacterial meta-resistome in all living and nonliving environments. As a result, reservoirs, sources and distribution of antibiotic resistance genes (ARGS) and antibiotic-producers, as well as the effects of human activity and antibiotics on the selection and dissemination of ARGs were not well comprehended. With the advances made in the fields of metagenomics and metatranscriptomics, many of the hitherto little-understood concepts are becoming clearer. Further, the discovery of antibiotics such as lugdinin and lactocillin from the human microbiota, buttressed the importance of these new fields. Metagenomics and metatranscriptomics are becoming important clinical diagnostic tools for screening and detecting pathogens and ARGs, assessing the effects of antibiotics, other xenobiotics and human activity on the environment, characterizing the microbiome and the environmental resistome with lesser turnaround time and decreasing cost, as well as discovering antibiotic-producers. However, challenges with accurate binning, skewed ARGs databases, detection of less abundant and allelic variants of ARGs and efficient mobilome characterization remain. Ongoing efforts in long-read, phased- and single-cell sequencing, strain-resolved binning, chromosomal-conformation capture, DNA-methylation binning and deep-learning bioinformatic approaches offer promising prospects in reconstructing complete strain-level genomes and mobilomes from metagenomes.

The effects of marine farm-scale sequentially integrated multi-trophic aquaculture systems on microbial community composition, prevalence of sulfonamide-resistant bacteria and sulfonamide resistance gene sul1

Aquaculture, one of the most important food production practices worldwide, faces serious challenges of mitigating the detrimental impacts of intensive farming on the environment and increased prevalence of antibiotic resistance. To develop an environment-friendly aquaculture system, a land-based and farm-scale sequentially integrated multi-trophic aquaculture (IMTA) system was established for farming Chanos chanos in southwestern Taiwan. In this system, fishes are cultured in combination with organic extractive shellfish and inorganic extractive seaweed. This study aimed to evaluate the prevalence of sulfonamide-resistant bacteria, microbial community structure, and occurrence of sulfonamide resistance genes in the IMTA and traditional aquaculture systems. Water and sediment samples were collected before raising and after harvesting C. chanos. Our results showed that the occurrence of sulfonamide-resistant phenotypes in the IMTA system was comparable with that in influent seawater, while the traditional system exhibited a high sulfonamide resistance rate. Additionally, the traditional system resulted in a deviation of the bacterial community structure from that of seawater. In the water samples from the IMTA system and influent seawater, Proteobacteria and Bacteroidetes were the two dominant phyla, representing approximately 75% and 15% of the community, respectively. In the traditional system, Actinobacteria, constituting 39% of the community, was the dominant bacterial phylum. Thirty-one sulfonamide-resistant bacterial species were isolated. In conclusion, a sequentially IMTA system showed superior ability to maintain the prevalence of antibiotic resistance and the integrity of the bacterial community structure compared to the traditional farming system, representing a potentially valuable aquaculture system for preserving the sustainability of the marine environment.

Using metagenomics to investigate human and environmental resistomes

Antibiotic resistance is a global health concern declared by the WHO as one of the largest threats to modern healthcare. In recent years, metagenomic DNA sequencing has started to be applied as a tool to study antibiotic resistance in different environments, including the human microbiota. However, a multitude of methods exist for metagenomic data analysis, and not all methods are suitable for the investigation of resistance genes, particularly if the desired outcome is an assessment of risks to human health. In this review, we outline the current state of methods for sequence handling, mapping to databases of resistance genes, statistical analysis and metagenomic assembly. In addition, we provide an overview of important considerations related to the analysis of resistance genes, and recommend some of the currently used tools and methods that are best equipped to inform research and clinical practice related to antibiotic resistance.

Variability in metagenomic samples from the Puget Sound: Relationship to temporal and anthropogenic impacts

Whole-metagenome sequencing (WMS) has emerged as a powerful tool to assess potential public health risks in marine environments by measuring changes in microbial community structure and function in uncultured bacteria. In addition to monitoring public health risks such as antibiotic resistance determinants, it is essential to measure predictors of microbial variation in order to identify natural versus anthropogenic factors as well as to evaluate reproducibility of metagenomic measurements.This study expands our previous metagenomic characterization of Puget Sound by sampling new nearshore environments including the Duwamish River, an EPA superfund site, and the Hood Canal, an area characterized by highly variable oxygen levels. We also resampled a wastewater treatment plant, nearshore and open ocean sites introducing a longitudinal component measuring seasonal and locational variations and establishing metagenomics sampling reproducibility. Microbial composition from samples collected in the open sound were highly similar within the same season and location across different years, while nearshore samples revealed multi-fold seasonal variation in microbial composition and diversity. Comparisons with recently sequenced predominant marine bacterial genomes helped provide much greater species level taxonomic detail compared to our previous study. Antibiotic resistance determinants and pollution and detoxification indicators largely grouped by location showing minor seasonal differences. Metal resistance, oxidative stress and detoxification systems showed no increase in samples proximal to an EPA superfund site indicating a lack of ecosystem adaptation to anthropogenic impacts. Taxonomic analysis of common sewage influent families showed a surprising similarity between wastewater treatment plant and open sound samples suggesting a low-level but pervasive sewage influent signature in Puget Sound surface waters. Our study shows reproducibility of metagenomic data sampling in multiple Puget Sound locations while establishing baseline measurements of antibiotic resistance determinants, pollution and detoxification systems. Combining seasonal and longitudinal data across these locations provides a foundation for evaluating variation in future studies.

Stormwater loadings of antibiotic resistance genes in an urban stream

Antibiotic resistance presents a critical public health challenge and the transmission of antibiotic resistance via environmental pathways continues to gain attention. Factors driving the spread of antibiotic resistance genes (ARGs) in surface water and sources of ARGs in urban stormwater have not been well-characterized. In this study, five ARGs (sul1, sul2, tet(O), tet(W), and erm(F)) were quantified throughout the duration of three storm runoff events in an urban inland stream. Storm loads of all five ARGs were significantly greater than during equivalent background periods. Neither fecal indicator bacteria measured (E. coli or enterococci) was significantly correlated with sul1, sul2, or erm(F), regardless of whether ARG concentration was absolute or normalized to 16S rRNA levels. Both E. coli and enterococci were correlated with the tetracycline resistance genes, tet(O) and tet(W). Next-generation shotgun metagenomic sequencing was conducted to more thoroughly characterize the resistome (i.e., full complement of ARGs) and profile the occurrence of all ARGs described in current databases in storm runoff in order to inform future watershed monitoring and management. Between 37 and 121 different ARGs were detected in each stream sample, though the ARG profiles differed among storms. This study establishes that storm-driven transport of ARGs comprises a considerable fraction of overall downstream loadings and broadly characterizes the urban stormwater resistome to identify potential marker ARGs indicative of impact.

Metagenomic analysis reveals wastewater treatment plants as hotspots of antibiotic resistance genes and mobile genetic elements

The intensive use of antibiotics results in their continuous release into the environment and the subsequent widespread occurrence of antibiotic resistant bacteria (ARB), antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs). This study used Illumina high-throughput sequencing to investigate the broad-spectrum profiles of both ARGs and MGEs in activated sludge and anaerobically digested sludge from a full-scale wastewater treatment plant. A pipeline for identifying antibiotic resistance determinants was developed that consisted of four categories: gene transfer potential, ARG potential, ARGs pathway and ARGs phylogenetic origin. The metagenomic analysis showed that the activated sludge and the digested sludge exhibited different microbial communities and changes in the types and occurrence of ARGs and MGEs. In total, 42 ARGs subtypes were identified in the activated sludge, while 51 ARG subtypes were detected in the digested sludge. Additionally, MGEs including plasmids, transposons, integrons (intI1) and insertion sequences (e.g. ISSsp4, ISMsa21 and ISMba16) were abundant in the two sludge samples. The co-occurrence pattern between ARGs and microbial taxa revealed by network analysis indicated that some environmental bacteria (e.g. Clostridium and Nitrosomonas) might be potential hosts of multiple ARGs. The findings increase our understanding of WWTPs as hotspots of ARGs and MGEs, and contribute towards preventing their release into the downstream environment.

Insertion sequences enrichment in extreme Red sea brine pool vent

Mobile genetic elements are major agents of genome diversification and evolution. Limited studies addressed their characteristics, including abundance, and role in extreme habitats. One of the rare natural habitats exposed to multiple-extreme conditions, including high temperature, salinity and concentration of heavy metals, are the Red Sea brine pools. We assessed the abundance and distribution of different mobile genetic elements in four Red Sea brine pools including the world's largest known multiple-extreme deep-sea environment, the Red Sea Atlantis II Deep. We report a gradient in the abundance of mobile genetic elements, dramatically increasing in the harshest environment of the pool. Additionally, we identified a strong association between the abundance of insertion sequences and extreme conditions, being highest in the harshest and deepest layer of the Red Sea Atlantis II Deep. Our comparative analyses of mobile genetic elements in secluded, extreme and relatively non-extreme environments, suggest that insertion sequences predominantly contribute to polyextremophiles genome plasticity.

Metagenomic analysis of bacterial community composition and antibiotic resistance genes in a wastewater treatment plant and its receiving surface water

The presence of pathogenic bacteria and the dissemination of antibiotic resistance genes (ARGs) may pose big risks to the rivers that receive the effluent from municipal wastewater treatment plants (WWTPs). In this study, we investigated the changes of bacterial community and ARGs along treatment processes of one WWTP, and examined the effects of the effluent discharge on the bacterial community and ARGs in the receiving river. Pyrosequencing was applied to reveal bacterial community composition including potential bacterial pathogen, and Illumina high-throughput sequencing was used for profiling ARGs. The results showed that the WWTP had good removal efficiency on potential pathogenic bacteria (especially Arcobacter butzleri) and ARGs. Moreover, the bacterial communities of downstream and upstream of the river showed no significant difference. However, the increase in the abundance of potential pathogens and ARGs at effluent outfall was observed, indicating that WWTP effluent might contribute to the dissemination of potential pathogenic bacteria and ARGs in the receiving river.

Laying a Community-Based Foundation for Data-Driven Semantic Standards in Environmental Health Sciences

BACKGROUND

Despite increasing availability of environmental health science (EHS) data, development, and implementation of relevant semantic standards, such as ontologies or hierarchical vocabularies, has lagged. Consequently, integration and analysis of information needed to better model environmental influences on human health remains a significant challenge.

OBJECTIVES

We aimed to identify a committed community and mechanisms needed to develop EHS semantic standards that will advance understanding about the impacts of environmental exposures on human disease.

METHODS

The National Institute of Environmental Health Sciences sponsored the "Workshop for the Development of a Framework for Environmental Health Science Language" hosted at North Carolina State University on 15-16 September 2014. Through the assembly of data generators, users, publishers, and funders, we aimed to develop a foundation for enabling the development of community-based and data-driven standards that will ultimately improve standardization, sharing, and interoperability of EHS information.

DISCUSSION

Creating and maintaining an EHS common language is a continuous and iterative process, requiring community building around research interests and needs, enabling integration and reuse of existing data, and providing a low barrier of access for researchers needing to use or extend such a resource.

CONCLUSIONS

Recommendations included developing a community-supported web-based toolkit that would enable a) collaborative development of EHS research questions and use cases, b) construction of user-friendly tools for searching and extending existing semantic resources, c) education and guidance about standards and their implementation, and d) creation of a plan for governance and sustainability.

CITATION

Mattingly CJ, Boyles R, Lawler CP, Haugen AC, Dearry A, Haendel M. 2016. Laying a community-based foundation for data-driven semantic standards in environmental health sciences. Environ Health Perspect 124:1136-1140; http://dx.doi.org/10.1289/ehp.1510438.

Metagenomic analysis of sediments under seaports influence in the Equatorial Atlantic Ocean

Maritime ports are anthropogenic interventions capable of causing serious alterations in coastal ecosystems. In this study, we examined the benthic microbial diversity and community structure under the influence of two maritime ports, Mucuripe (MUC) and Pecém (PEC), at Equatorial Atlantic Ocean in Northeast Brazil. Those seaports differ in architecture, time of functioning, cargo handling and contamination. The microbiomes from MUC and PEC were also compared in silico to 11 other globally distributed marine microbiomes. The comparative analysis of operational taxonomic units (OTUs) retrieved by PCR-DGGE showed that MUC presents greater richness and β diversity of Bacteria and Archaea than PEC. In line with these results, metagenomic analysis showed that MUC and PEC benthic microbial communities share the main common bacterial phyla found in coastal environments, although can be distinguish by greater abundance of Cyanobacteria in MUC and Deltaproteobacteria in PEC. Both ports differed in Archaea composition, being PEC port sediments dominated by Thaumarchaeota. The microbiomes showed little divergence in their potential metabolic pathways, although shifts on the microbial taxonomic signatures involved in nitrogen and sulphur metabolic pathways were observed. The comparative analysis of different benthic marine metagenomes from Brazil, Australia and Mexico grouped them by the geographic location rather than by the type of ecosystem, although at phylum level seaport sediments share a core microbiome constituted by Proteobacteria, Cyanobacteria, Actinobacteria, Tenericuteres, Firmicutes, Bacteriodetes and Euryarchaeota. Our results suggest that multiple physical and chemical factors acting on sediments as a result of at least 60years of port operation play a role in shaping the benthic microbial communities at taxonomic level, but not at functional level.

Molecular Methods for Assessment of Antibiotic Resistance in Agricultural Ecosystems: Prospects and Challenges

Agricultural ecosystems are of special interest for monitoring the potential for antibiotic resistance to spread through the environment and contribute to human exposure. Molecular methods, which target DNA, RNA, and other molecular components of bacterial cells, present certain advantages for characterizing and quantifying markers of antibiotic resistance and their horizontal gene transfer. These include rapid, unambiguous detection of targets; consistent results; and avoidance of culture bias. However, molecular methods are also subject to limitations that are not always clearly addressed or taken into consideration in the interpretation of scientific data. In particular, DNA-based methods do not directly assess viability or presence within an intact bacterial host, but such information may be inferred based on appropriate experimental design or in concert with complementary methods. The purpose of this review is to provide an overview of existing molecular methods for tracking antibiotic resistance in agricultural ecosystems, to define their strengths and weaknesses, and to recommend a path forward for future applications of molecular methods and standardized reporting in the literature. This will guide research along the farm-to-fork continuum and support comparability of the growing number of studies in the literature in a manner that informs management decisions and policy development.

The ocean as a global reservoir of antibiotic resistance genes

Recent studies of natural environments have revealed vast genetic reservoirs of antibiotic resistance (AR) genes. Soil bacteria and human pathogens share AR genes, and AR genes have been discovered in a variety of habitats. However, there is little knowledge about the presence and diversity of AR genes in marine environments and which organisms host AR genes. To address this, we identified the diversity of genes conferring resistance to ampicillin, tetracycline, nitrofurantoin, and sulfadimethoxine in diverse marine environments using functional metagenomics (the cloning and screening of random DNA fragments). Marine environments were host to a diversity of AR-conferring genes. Antibiotic-resistant clones were found at all sites, with 28% of the genes identified as known AR genes (encoding beta-lactamases, bicyclomycin resistance pumps, etc.). However, the majority of AR genes were not previously classified as such but had products similar to proteins such as transport pumps, oxidoreductases, and hydrolases. Furthermore, 44% of the genes conferring antibiotic resistance were found in abundant marine taxa (e.g., Pelagibacter, Prochlorococcus, and Vibrio). Therefore, we uncovered a previously unknown diversity of genes that conferred an AR phenotype among marine environments, which makes the ocean a global reservoir of both clinically relevant and potentially novel AR genes.

Does human activity impact the natural antibiotic resistance background? Abundance of antibiotic resistance genes in 21 Swiss lakes

Antibiotic resistance genes (ARGs) are emerging environmental contaminants, known to be continuously discharged into the aquatic environment via human and animal waste. Freshwater aquatic environments represent potential reservoirs for ARG and potentially allow sewage-derived ARG to persist and spread in the environment. This may create increased opportunities for an eventual contact with, and gene transfer to, human and animal pathogens via the food chain or drinking water. However, assessment of this risk requires a better understanding of the level and variability of the natural resistance background and the extent of the human impact. We have analyzed water samples from 21 Swiss lakes, taken at sampling points that were not under the direct influence of local contamination sources and analyzed the relative abundance of ARG using quantitative real-time PCR. Copy numbers of genes mediating resistance to three different broad-spectrum antibiotic classes (sulfonamides: sul1, sul2, tetracyclines: tet(B), tet(M), tet(W) and fluoroquinolones: qnrA) were normalized to copy numbers of bacterial 16S rRNA genes. We used multiple linear regression to assess if ARG abundance is related to human activities in the catchment, microbial community composition and the eutrophication status of the lakes. Sul genes were detected in all sampled lakes, whereas only four lakes contained quantifiable numbers of tet genes, and qnrA remained below detection in all lakes. Our data indicate higher abundance of sul1 in lakes with increasing number and capacity of wastewater treatment plants (WWTPs) in the catchment. sul2 abundance was rather related to long water residence times and eutrophication status. Our study demonstrates the potential of freshwater lakes to preserve antibiotic resistance genes, and provides a reference for ARG abundance from lake systems with low human impact as a baseline for assessing ARG contamination in lake water.

Methods for Detecting the Environmental Coccoid Form of Helicobacter pylori

Helicobacter pylori is recognized as the most common pathogen to cause gastritis, peptic and duodenal ulcers, and gastric cancer. The organisms are found in two forms: (1) spiral-shaped bacillus and (2) coccoid. H. pylori coccoid form, generally found in the environment, is the transformed form of the normal spiral-shaped bacillus after exposed to water or adverse environmental conditions such as exposure to sub-inhibitory concentrations of antimicrobial agents. The putative infectious capability and the viability of H. pylori under environmental conditions are controversial. This disagreement is partially due to the fact of lack in detecting the coccoid form of H. pylori in the environment. Accurate and effective detection methods of H. pylori will lead to rapid treatment and disinfection, and less human health damages and reduction in health care costs. In this review, we provide a brief introduction to H. pylori environmental coccoid forms, their transmission, and detection methods. We further discuss the use of these detection methods including their accuracy and efficiency.

Influence of a non-hospital medical care facility on antimicrobial resistance in wastewater

The global widespread use of antimicrobials and accompanying increase in resistant bacterial strains is of major public health concern. Wastewater systems and wastewater treatment plants are considered a niche for antibiotic resistance genes (ARGs), with diverse microbial communities facilitating ARG transfer via mobile genetic element (MGE). In contrast to hospital sewage, wastewater from other health care facilities is still poorly investigated. At the instance of a nursing home located in south-west Germany, in the present study, shotgun metagenomics was used to investigate the impact on wastewater of samples collected up- and down-stream in different seasons. Microbial composition, ARGs and MGEs were analyzed using different annotation approaches with various databases, including Antibiotic Resistance Ontologies (ARO), integrons and plasmids. Our analysis identified seasonal differences in microbial communities and abundance of ARG and MGE between samples from different seasons. However, no obvious differences were detected between up- and downstream samples. The results suggest that, in contrast to hospitals, sewage from the nursing home does not have a major impact on ARG or MGE in wastewater, presumably due to much less intense antimicrobial usage. Possible limitations of metagenomic studies using high-throughput sequencing for detection of genes that seemingly confer antibiotic resistance are discussed.

Shotgun metagenomics reveals a wide array of antibiotic resistance genes and mobile elements in a polluted lake in India

There is increasing evidence for an environmental origin of many antibiotic resistance genes. Consequently, it is important to identify environments of particular risk for selecting and maintaining such resistance factors. In this study, we described the diversity of antibiotic resistance genes in an Indian lake subjected to industrial pollution with fluoroquinolone antibiotics. We also assessed the genetic context of the identified resistance genes, to try to predict their genetic transferability. The lake harbored a wide range of resistance genes (81 identified gene types) against essentially every major class of antibiotics, as well as genes responsible for mobilization of genetic material. Resistance genes were estimated to be 7000 times more abundant than in a Swedish lake included for comparison, where only eight resistance genes were found. The sul2 and qnrD genes were the most common resistance genes in the Indian lake. Twenty-six known and 21 putative novel plasmids were recovered in the Indian lake metagenome, which, together with the genes found, indicate a large potential for horizontal gene transfer through conjugation. Interestingly, the microbial community of the lake still included a wide range of taxa, suggesting that, across most phyla, bacteria has adapted relatively well to this highly polluted environment. Based on the wide range and high abundance of known resistance factors we have detected, it is plausible that yet unrecognized resistance genes are also present in the lake. Thus, we conclude that environments polluted with waste from antibiotic manufacturing could be important reservoirs for mobile antibiotic resistance genes.

Metagenomic frameworks for monitoring antibiotic resistance in aquatic environments

BACKGROUND

High-throughput genomic technologies offer new approaches for environmental health monitoring, including metagenomic surveillance of antibiotic resistance determinants (ARDs). Although natural environments serve as reservoirs for antibiotic resistance genes that can be transferred to pathogenic and human commensal bacteria, monitoring of these determinants has been infrequent and incomplete. Furthermore, surveillance efforts have not been integrated into public health decision making.

OBJECTIVES

We used a metagenomic epidemiology-based approach to develop an ARD index that quantifies antibiotic resistance potential, and we analyzed this index for common modal patterns across environmental samples. We also explored how metagenomic data such as this index could be conceptually framed within an early risk management context.

METHODS

We analyzed 25 published data sets from shotgun pyrosequencing projects. The samples consisted of microbial community DNA collected from marine and freshwater environments across a gradient of human impact. We used principal component analysis to identify index patterns across samples.

RESULTS

We observed significant differences in the overall index and index subcategory levels when comparing ecosystems more proximal versus distal to human impact. The selection of different sequence similarity thresholds strongly influenced the index measurements. Unique index subcategory modes distinguished the different metagenomes.

CONCLUSIONS

Broad-scale screening of ARD potential using this index revealed utility for framing environmental health monitoring and surveillance. This approach holds promise as a screening tool for establishing baseline ARD levels that can be used to inform and prioritize decision making regarding management of ARD sources and human exposure routes.

CITATION

Port JA, Cullen AC, Wallace JC, Smith MN, Faustman EM. 2014. Metagenomic frameworks for monitoring antibiotic resistance in aquatic environments. Environ Health Perspect 122:222–228; http://dx.doi.org/10.1289/ehp.1307009

Metagenomic profiles of antibiotic resistance genes (ARGs) between human impacted estuary and deep ocean sediments

Knowledge of the origins and dissemination of antibiotic resistance genes (ARGs) is essential for understanding modern resistomes in the environment. The mechanisms of the dissemination of ARGs can be revealed through comparative studies on the metagenomic profiling of ARGs between relatively pristine and human-impacted environments. The deep ocean bed of the South China Sea (SCS) is considered to be largely devoid of anthropogenic impacts, while the Pearl River Estuary (PRE) in south China has been highly impacted by intensive human activities. Commonly used antibiotics (sulfamethazine, norfloxacin, ofloxacin, tetracycline, and erythromycin) have been detected through chemical analysis in the PRE sediments, but not in the SCS sediments. In the relatively pristine SCS sediments, the most prevalent and abundant ARGs are those related to resistance to macrolides and polypeptides, with efflux pumps as the predominant mechanism. In the contaminated PRE sediments, the typical ARG profiles suggest a prevailing resistance to antibiotics commonly used in human health and animal farming (including sulfonamides, fluoroquinolones, and aminoglycosides), and higher diversity in both genotype and resistance mechanism than those in the SCS. In particular, antibiotic inactivation significantly contributed to the resistance to aminoglycosides, β-lactams, and macrolides observed in the PRE sediments. There was a significant correlation in the levels of abundance of ARGs and those of mobile genetic elements (including integrons and plasmids), which serve as carriers in the dissemination of ARGs in the aquatic environment. The metagenomic results from the current study support the view that ARGs naturally originate in pristine environments, while human activities accelerate the dissemination of ARGs so that microbes would be able to tolerate selective environmental stress in response to anthropogenic impacts.

Metagenomic profiling of antibiotic resistance genes and mobile genetic elements in a tannery wastewater treatment plant

Antibiotics are often used to prevent sickness and improve production in animal agriculture, and the residues in animal bodies may enter tannery wastewater during leather production. This study aimed to use Illumina high-throughput sequencing to investigate the occurrence, diversity and abundance of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) in aerobic and anaerobic sludge of a full-scale tannery wastewater treatment plant (WWTP). Metagenomic analysis showed that Proteobacteria, Firmicutes, Bacteroidetes and Actinobacteria dominated in the WWTP, but the relative abundance of archaea in anaerobic sludge was higher than in aerobic sludge. Sequencing reads from aerobic and anaerobic sludge revealed differences in the abundance of functional genes between both microbial communities. Genes coding for antibiotic resistance were identified in both communities. BLAST analysis against Antibiotic Resistance Genes Database (ARDB) further revealed that aerobic and anaerobic sludge contained various ARGs with high abundance, among which sulfonamide resistance gene sul1 had the highest abundance, occupying over 20% of the total ARGs reads. Tetracycline resistance genes (tet) were highly rich in the anaerobic sludge, among which tet33 had the highest abundance, but was absent in aerobic sludge. Over 70 types of insertion sequences were detected in each sludge sample, and class 1 integrase genes were prevalent in the WWTP. The results highlighted prevalence of ARGs and MGEs in tannery WWTPs, which may deserve more public health concerns.