Exploring variation of antibiotic resistance genes in activated sludge over a four-year period through a metagenomic approach

Impact of dairy manure pre-application treatment on manure composition, soil dynamics of antibiotic resistance genes, and abundance of antibiotic-resistance genes on vegetables at harvest

Manuring ground used for crop production is an important agricultural practice. Should antibiotic-resistant enteric bacteria carried in the manure be transferred to crops that are consumed raw, their consumption by humans or animals will represent a route of exposure to antibiotic resistance genes. Treatment of manures prior to land application is a potential management option to reduce the abundance of antibiotic resistance genes entrained with manure application. In this study, dairy manure that was untreated, anaerobically digested, mechanically dewatered or composted was applied to field plots that were then cropped to lettuce, carrots and radishes. The impact of treatment on manure composition, persistence of antibiotic resistance gene targets in soil following application, and distribution of antibiotic resistance genes and bacteria on vegetables at harvest was determined. Composted manure had the lowest abundance of antibiotic resistance gene targets compared to the other manures. There was no significant difference in the persistence characteristics of antibiotic resistance genes following land application of the various manures. Compared to unmanured soil, antibiotic resistance genes were detected more frequently in soil receiving raw or digested manure, whereas they were not in soil receiving composted manure. The present study suggests that vegetables grown in ground receiving raw or digested manure are at risk of contamination with manure-borne antibiotic resistant bacteria, whereas vegetables grown in ground receiving composted manure are less so.

Microbial community analysis of bulk sludge/cake layers and biofouling-causing microbial consortia in a full-scale aerobic membrane bioreactor

Pyrosequencing was used to investigate biofouling-causing microbial consortia at the community level in bulk sludge and cake layers within a full-scale membrane bioreactor (MBR). The analysis revealed Chao's estimates of total operational taxonomic units (OTUs) of 1726, 1806, and 1362 for bulk sludge, cake outer layer, and cake inner layer, respectively. The bulk sludge and cake outer layer OTUs clustered together, whereas the cake inner layer OTUs formed a separate group, indicating that environmental conditions affected the microbial community composition within the MBR. Bacteroidetes, Proteobacteria, and Chloroflexi were the dominant phyla in both the bulk sludge and the cake layers. Comparison at the genus level showed twelve distinct genera in the cake layers that were absent in bulk sludge. Twenty distinct genera were recorded in the inner cake layer. Those genera are likely the microbial colonization pioneers in full-scale membrane bioreactors.

Microbial diversity of a full-scale UASB reactor applied to poultry slaughterhouse wastewater treatment: integration of 16S rRNA gene amplicon and shotgun metagenomic sequencing

The 16S rRNA gene amplicon and whole-genome shotgun metagenomic (WGSM) sequencing approaches were used to investigate wide-spectrum profiles of microbial composition and metabolic diversity from a full-scale UASB reactor applied to poultry slaughterhouse wastewater treatment. The data were generated by using MiSeq 2 × 250 bp and HiSeq 2 × 150 bp Illumina sequencing platforms for 16S amplicon and WGSM sequencing, respectively. Each approach revealed a distinct microbial community profile, with Pseudomonas and Psychrobacter as predominant genus for the WGSM dataset and Clostridium and Methanosaeta for the 16S rRNA gene amplicon dataset. The virome characterization revealed the presence of two viral families with Bacteria and Archaea as host, Myoviridae, and Siphoviridae. A wide functional diversity was found with predominance of genes involved in the metabolism of acetone, butanol, and ethanol synthesis; and one-carbon metabolism (e.g., methanogenesis). Genes related to the acetotrophic methanogenesis pathways were more abundant than methylotrophic and hydrogenotrophic, corroborating the taxonomic results that showed the prevalence of the acetotrophic genus Methanosaeta. Moreover, the dataset indicated a variety of metabolic genes involved in sulfur, nitrogen, iron, and phosphorus cycles, with many genera able to act in all cycles. BLAST analysis against Antibiotic Resistance Genes Database (ARDB) revealed that microbial community contained 43 different types of antibiotic resistance genes, some of them were associated with growth chicken promotion (e.g., bacitracin, tetracycline, and polymyxin).

Shifts in microbial community structure and function in surface waters impacted by unconventional oil and gas wastewater revealed by metagenomics

Unconventional oil and gas (UOG) production produces large quantities of wastewater with complex geochemistry and largely uncharacterized impacts on surface waters. In this study, we assessed shifts in microbial community structure and function in sediments and waters upstream and downstream from a UOG wastewater disposal facility. To do this, quantitative PCR for 16S rRNA and antibiotic resistance genes along with metagenomic sequencing were performed. Elevated conductivity and markers of UOG wastewater characterized sites sampled downstream from the disposal facility compared to background sites. Shifts in overall high level functions and microbial community structure were observed between background sites and downstream sediments. Increases in Deltaproteobacteria and Methanomicrobia and decreases in Thaumarchaeota were observed at downstream sites. Genes related to dormancy and sporulation and methanogenic respiration were 18-86 times higher at downstream, impacted sites. The potential for these sediments to serve as reservoirs of antimicrobial resistance was investigated given frequent reports of the use of biocides to control the growth of nuisance bacteria in UOG operations. A shift in resistance profiles downstream of the UOG facility was observed including increases in acrB and mexB genes encoding for multidrug efflux pumps, but not overall abundance of resistance genes. The observed shifts in microbial community structure and potential function indicate changes in respiration, nutrient cycling, and markers of stress in a stream impacted by UOG waste disposal operations.

Functional Metagenomics as a Tool for Identification of New Antibiotic Resistance Genes from Natural Environments

Antibiotic resistance has become a major concern for human and animal health, as therapeutic alternatives to treat multidrug-resistant microorganisms are rapidly dwindling. The problem is compounded by low investment in antibiotic research and lack of new effective antimicrobial drugs on the market. Exploring environmental antibiotic resistance genes (ARGs) will help us to better understand bacterial resistance mechanisms, which may be the key to identifying new drug targets. Because most environment-associated microorganisms are not yet cultivable, culture-independent techniques are essential to determine which organisms are present in a given environmental sample and allow the assessment and utilization of the genetic wealth they represent. Metagenomics represents a powerful tool to achieve these goals using sequence-based and functional-based approaches. Functional metagenomic approaches are particularly well suited to the identification new ARGs from natural environments because, unlike sequence-based approaches, they do not require previous knowledge of these genes. This review discusses functional metagenomics-based ARG research and describes new possibilities for surveying the resistome in environmental samples.

Polycyclic aromatic hydrocarbons (PAHs) enriching antibiotic resistance genes (ARGs) in the soils

The prevalence of antibiotic resistance genes (ARGs) in modern environment raises an emerging global health concern. In this study, soil samples were collected from three sites in petrochemical plant that represented different pollution levels of polycyclic aromatic hydrocarbons (PAHs). Metagenomic profiling of these soils demonstrated that ARGs in the PAHs-contaminated soils were approximately 15 times more abundant than those in the less-contaminated ones, with Proteobacterial being the preponderant phylum. Resistance profile of ARGs in the PAHs-polluted soils was characterized by the dominance of efflux pump-encoding ARGs associated with aromatic antibiotics (e.g., fluoroquinolones and acriflavine) that accounted for more than 70% of the total ARGs, which was significantly different from representative sources of ARG pollution due to wide use of antibiotics. Most of ARGs enriched in the PAHs-contaminated soils were not carried by plasmids, indicating the low possibilities of them being transferred between bacteria. Significant correlation was observed between the total abundance of ARGs and that of Proteobacteria in the soils. Proteobacteria selected by PAHs led to simultaneously enriching of ARGs carried by them in the soils. Our results suggested that PAHs could serve as one of selective stresses for greatly enriching of ARGs in the human-impacted environment.

Co-occurrence of antibiotic and metal resistance genes revealed in complete genome collection

The high frequency of antibiotic resistance is a global public health concern. More seriously, widespread metal pressure in the environment may facilitate the proliferation of antibiotic resistance via coselection of antibiotic resistance genes (ARGs) and metal resistance genes (MRGs). Given the lack of comprehensive understanding of the ARG and MRG coselection, in this study both abundance relationship and genetic linkage between ARGs and MRGs were rigorously investigated by performing a genomic analysis of a large complete genome collection. Many more ARGs were enriched in human-associated bacteria compared with those subjected to less anthropogenic interference. The signatures of ARG and MRG co-occurrence were much more frequent and the distance linkages between ARGs and MRGs were much more intimate in human pathogens than those less human-associated bacteria. Moreover, the co-occurrence structures in the habitat divisions were significantly different, which could be attributed to their distinct gene transfer potentials. More exogenous ARGs and MRGs on the genomes of human pathogens indicated the importance of recent resistance acquisition in resistome development of human commensal flora. Overall, the study emphasizes the potential risk associated with ARG and MRG coselection of both environmental and medical relevance.

The Composition and Spatial Patterns of Bacterial Virulence Factors and Antibiotic Resistance Genes in 19 Wastewater Treatment Plants

Bacterial pathogenicity and antibiotic resistance are of concern for environmental safety and public health. Accumulating evidence suggests that wastewater treatment plants (WWTPs) are as an important sink and source of pathogens and antibiotic resistance genes (ARGs). Virulence genes (encoding virulence factors) are good indicators for bacterial pathogenic potentials. To achieve a comprehensive understanding of bacterial pathogenic potentials and antibiotic resistance in WWTPs, bacterial virulence genes and ARGs in 19 WWTPs covering a majority of latitudinal zones of China were surveyed by using GeoChip 4.2. A total of 1610 genes covering 13 virulence factors and 1903 genes belonging to 11 ARG families were detected respectively. The bacterial virulence genes exhibited significant spatial distribution patterns of a latitudinal biodiversity gradient and a distance-decay relationship across China. Moreover, virulence genes tended to coexist with ARGs as shown by their strongly positive associations. In addition, key environmental factors shaping the overall virulence gene structure were identified. This study profiles the occurrence, composition and distribution of virulence genes and ARGs in current WWTPs in China, and uncovers spatial patterns and important environmental variables shaping their structure, which may provide the basis for further studies of bacterial virulence factors and antibiotic resistance in WWTPs.

Antimicrobial resistance surveillance in the genomic age

The loss of effective antimicrobials is reducing our ability to protect the global population from infectious disease. However, the field of antibiotic drug discovery and the public health monitoring of antimicrobial resistance (AMR) is beginning to exploit the power of genome and metagenome sequencing. The creation of novel AMR bioinformatics tools and databases and their continued development will advance our understanding of the molecular mechanisms and threat severity of antibiotic resistance, while simultaneously improving our ability to accurately predict and screen for antibiotic resistance genes within environmental, agricultural, and clinical settings. To do so, efforts must be focused toward exploiting the advancements of genome sequencing and information technology. Currently, AMR bioinformatics software and databases reflect different scopes and functions, each with its own strengths and weaknesses. A review of the available tools reveals common approaches and reference data but also reveals gaps in our curated data, models, algorithms, and data-sharing tools that must be addressed to conquer the limitations and areas of unmet need within the AMR research field before DNA sequencing can be fully exploited for AMR surveillance and improved clinical outcomes.

Metagenomic insights into ultraviolet disinfection effects on antibiotic resistome in biologically treated wastewater

High-throughput sequencing-based metagenomic approaches were used to comprehensively investigate ultraviolet effects on the microbial community structure, and diversity and abundance of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) in biologically treated wastewater. After ultraviolet radiation, some dominant genera, like Aeromonas and Halomonas, in the wastewater almost disappeared, while the relative abundance of some minor genera including Pseudomonas and Bacillus increased dozens of times. Metagenomic analysis showed that 159 ARGs within 14 types were detectable in the samples, and the radiation at 500 mJ/cm(2) obviously increased their total relative abundance from 31.68 ppm to 190.78 ppm, which was supported by quantitative real time PCR. As the dominant persistent ARGs, multidrug resistance genes carried by Pseudomonas and bacitracin resistance gene bacA carried by Bacillus mainly contributed to the ARGs abundance increase. Bacterial community shift and MGEs replication induced by the radiation might drive the resistome alteration. The findings may shed new light on the mechanism behind the ultraviolet radiation effects on antibiotic resistance in wastewater.

Can chlorination co-select antibiotic-resistance genes?

Selective pressures, such as chemical or heavy metal pollution, may co-select for bacterial antibiotic resistance in the environment. However, whether chlorination in water treatment can co-select antibiotic-resistant bacteria is controversial. In this study, high capacity quantitative polymerase chain reaction (qPCR) analysis was applied to target almost all known antibiotic-resistance genes (ARGs) (282 types) and 13 mobile genetic elements (MGEs) in bacteria detected in secondary effluents from a municipal wastewater treatment plant after chlorination. The results revealed that 125 unique ARGs were detected in non-chlorinated samples, and the number decreased (79-91 types) as the chlorine concentration was increased. Moreover, 7.49 × 10(4)-3.92 × 10(7) copies/100 ml water reduction of ARGs occurred with 4 mg Cl2/l. Considering the relative abundance of ARGs (i.e., ARG copies normalized to 16S rRNA gene copies), 119 ARGs decreased in response to chlorination, whereas only six ARGs, such as dfrA1, tetPB-03, tetPA, ampC-04, tetA-02, and erm(36), were potentially enriched by 10.90-, 10.06-, 8.63-, 6.86-, 3.77-, and 1.09-fold, respectively. Furthermore, the relative abundance of 12 detected MGEs was lower after chlorination. Therefore, chlorination was effective in reducing ARGs and MGEs rather than co-selecting them.

Metagenomic Analysis Revealing Antibiotic Resistance Genes (ARGs) and Their Genetic Compartments in the Tibetan Environment

Comprehensive profiles of antibiotic resistance genes (ARGs) and mobile genetic elements (MGEs) in a minimally impacted environment are essential to understanding the evolution and dissemination of modern antibiotic resistance. Chemical analyses of the samples collected from Tibet demonstrated that the region under investigation was almost devoid of anthropogenic antibiotics. The soils, animal wastes, and sediments were different from each other in terms of bacterial community structures, and in the typical profiles of ARGs and MGEs. Diverse ARGs that encoded resistance to common antibiotics (e.g., beta-lactams, fluoroquinolones, etc.) were found mainly via an efflux mechanism completely distinct from modern antibiotic resistome. In addition, a very small fraction of ARGs in the Tibetan environment were carried by MGEs, indicating the low potential of these ARGs to be transferred among bacteria. In comparison to the ARG profiles in relatively pristine Tibet, contemporary ARGs and MGEs in human-impacted environments have evolved substantially since the broad use of anthropogenic antibiotics.

Changes of resistome, mobilome and potential hosts of antibiotic resistance genes during the transformation of anaerobic digestion from mesophilic to thermophilic

This study aimed to reveal how antibiotic resistance genes (ARGs) and their horizontal and vertical transfer-related items (mobilome and bacterial hosts) respond to the transformation of anaerobic digestion (AD) from mesophilic to thermophilic using one-step temperature increase. The resistomes and mobilomes of mesophilic and thermophilic sludge were investigated using metagenome sequencing, and the changes in 24 representative ARGs belonging to three categories, class 1 integron and bacterial genera during the transition period were further followed using quantitative PCR and 454-pyrosequencing. After the temperature increase, resistome abundance in the digested sludge decreased from 125.97 ppm (day 0, mesophilic) to 50.65 ppm (day 57, thermophilic) with the reduction of most ARG types except for the aminoglycoside resistance genes. Thermophilic sludge also had a smaller mobilome, including plasmids, insertion sequences and integrons, than that of mesophilic sludge, suggesting the lower horizontal transfer potential of ARGs under thermophilic conditions. On the other hand, the total abundance of 18 bacterial genera, which were suggested as the possible hosts for 13 ARGs through network analysis, decreased from 23.27% in mesophilic sludge to 11.92% in thermophilic sludge, indicating fewer hosts for the vertical expansion of ARGs after the increase in temperature. These results indicate that the better reduction of resistome abundance by thermophilic AD might be associated with the decrease of both the horizontal and vertical transferability of ARGs.

Dissemination of veterinary antibiotics and corresponding resistance genes from a concentrated swine feedlot along the waste treatment paths

Swine feedlots are an important pollution source of antibiotics and antibiotic resistance genes (ARGs) to the environment. This study investigated the dissemination of two classes of commonly-used veterinary antibiotics, namely, tetracyclines (TCs) and sulfonamides (SAs), and their corresponding ARGs along the waste treatment paths from a concentrated swine feedlot located in Beijing, China. The highest total TC and total SA concentrations detected were 166.7mgkg(-1) and 64.5μgkg(-1) in swine manure as well as 388.7 and 7.56μgL(-1) in swine wastewater, respectively. Fourteen tetracycline resistance genes (TRGs) encoding ribosomal protection proteins (RPP), efflux proteins (EFP) and enzymatic inactivation proteins, three sulfonamide resistance genes (SRGs), and two integrase genes were detected along the waste treatment paths with detection frequencies of 33.3-75.0%. The relative abundances of target ARGs ranged from 2.74×10(-6) to 1.19. The antibiotics and ARGs generally declined along both waste treatment paths, but their degree of reduction was more significant along the manure treatment path. The RPP TRGs dominated in the upstream samples and then decreased continuously along both waste treatment paths, whilst the EFP TRGs and SRGs maintained relatively stable. Strong correlations between antibiotic concentrations and ARGs were observed among both manure and wastewater samples. In addition, seasonal temperature, and integrase genes, moisture content and nutrient level of tested samples could all impact the relative abundances of ARGs along the swine waste treatment paths. This study helps understand the evolution and spread of ARGs from swine feedlots to the environment as well as assess the environmental risk arising from swine waste treatment.

High Levels of Antibiotic Resistance Genes and Their Correlations with Bacterial Community and Mobile Genetic Elements in Pharmaceutical Wastewater Treatment Bioreactors

To understand the diversity and abundance of antibiotic resistance genes (ARGs) in pharmaceutical wastewater treatment bioreactors, the ARGs in sludge from two full-scale pharmaceutical wastewater treatment plants (PWWTPs) were investigated and compared with sludge samples from three sewage treatment plants (STPs) using metagenomic approach. The results showed that the ARG abundances in PWWTP sludge ranged from 54.7 to 585.0 ppm, which were higher than those in STP sludge (27.2 to 86.4 ppm). Moreover, the diversity of ARGs in PWWTP aerobic sludge (153 subtypes) was higher than that in STP aerobic sludge (118 subtypes). In addition, it was found that the profiles of ARGs in PWWTP aerobic sludge were similar to those in STP aerobic sludge but different from those in PWWTP anaerobic sludge, suggesting that dissolve oxygen (DO) could be one of the important factors affecting the profiles of ARGs. In PWWTP aerobic sludge, aminoglycoside, sulfonamide and multidrug resistance genes were frequently detected. While, tetracycline, macrolide-lincosamide-streptogramin and polypeptide resistance genes were abundantly present in PWWTP anaerobic sludge. Furthermore, we investigated the microbial community and the correlation between microbial community and ARGs in PWWTP sludge. And, significant correlations between ARG types and seven bacterial genera were found. In addition, the mobile genetic elements (MGEs) were also examined and correlations between the ARGs and MGEs in PWWTP sludge were observed. Collectively, our results suggested that the microbial community and MGEs, which could be affected by DO, might be the main factors shaping the profiles of ARGs in PWWTP sludge.

Impacts of human activities on distribution of sulfate-reducing prokaryotes and antibiotic resistance genes in marine coastal sediments of Hong Kong

Sulfate-reducing prokaryotes (SRPs) and antibiotic resistance genes (ARGs) in sediments could be biomarkers for evaluating the environmental impacts of human activities, although factors governing their distribution are not clear yet. By using metagenomic approach, this study investigated the distributions of SRPs and ARGs in marine sediments collected from 12 different coastal locations of Hong Kong, which exhibited different pollution levels and were classified into two groups based on sediment parameters. Our results showed that relative abundances of major SRP genera to total prokaryotes were consistently lower in the more seriously polluted sediments (P-value < 0.05 in 13 of 20 genera), indicating that the relative abundance of SRPs is a negatively correlated biomarker for evaluating human impacts. Moreover, a unimodel distribution pattern for SRPs along with the pollution gradient was observed. Although total ARGs were enriched in sediments from the polluted sites, distribution of single major ARG types could be explained neither by individual sediment parameters nor by corresponding concentration of antibiotics. It supports the hypothesis that the persistence of ARGs in sediments may not need the selection of antibiotics. In summary, our study provided important hints of the niche differentiation of SRPs and behavior of ARGs in marine coastal sediment.

Monitoring and evaluation of antibiotic resistance genes in four municipal wastewater treatment plants in Harbin, Northeast China

The development and proliferation of antibiotic resistance in pathogenic and environmental microorganisms is of great concern for public health. In this study, the distribution and removal efficiency of intI1 and eight subtypes of antibiotic resistance genes (ARGs) for tetracycline, sulfonamides, beta-lactams resistance in four municipal wastewater treatment plants (WWTPs) in Harbin, which locates in Songhua River basin in cold areas of China, were monitored by real-time fluorescent quantitative PCR. The results showed that intI1 and 6 ARGs except for blaTEM and blaSHV were detected in wastewater and sludge samples and 0.3-2.7 orders of magnitude of ARGs removal efficiency in the four WWTPs were observed. The investigation on the removal of ARGs of different treatment units in one WWTP showed that the biological treatment unit played the most important role in ARGs removal (1.2-1.8 orders of magnitude), followed by UV disinfection, while primary physical treatment units can hardly remove any ARGs. Although all the WWTPs can remove ARGs effectively, ARGs concentrations are still relatively high in the effluent, their further attenuation should be investigated.

Metagenomic analysis revealed highly diverse microbial arsenic metabolism genes in paddy soils with low-arsenic contents

Microbe-mediated arsenic (As) metabolism plays a critical role in global As cycle, and As metabolism involves different types of genes encoding proteins facilitating its biotransformation and transportation processes. Here, we used metagenomic analysis based on high-throughput sequencing and constructed As metabolism protein databases to analyze As metabolism genes in five paddy soils with low-As contents. The results showed that highly diverse As metabolism genes were present in these paddy soils, with varied abundances and distribution for different types and subtypes of these genes. Arsenate reduction genes (ars) dominated in all soil samples, and significant correlation existed between the abundance of arr (arsenate respiration), aio (arsenite oxidation), and arsM (arsenite methylation) genes, indicating the co-existence and close-relation of different As resistance systems of microbes in wetland environments similar to these paddy soils after long-term evolution. Among all soil parameters, pH was an important factor controlling the distribution of As metabolism gene in five paddy soils (p = 0.018). To the best of our knowledge, this is the first study using high-throughput sequencing and metagenomics approach in characterizing As metabolism genes in the five paddy soil, showing their great potential in As biotransformation, and therefore in mitigating arsenic risk to humans.

Antibiotic Resistome: Improving Detection and Quantification Accuracy for Comparative Metagenomics

The unprecedented rise of life-threatening antibiotic resistance (AR), combined with the unparalleled advances in DNA sequencing of genomes and metagenomes, has pushed the need for in silico detection of the resistance potential of clinical and environmental metagenomic samples through the quantification of AR genes (i.e., genes conferring antibiotic resistance). Therefore, determining an optimal methodology to quantitatively and accurately assess AR genes in a given environment is pivotal. Here, we optimized and improved existing AR detection methodologies from metagenomic datasets to properly consider AR-generating mutations in antibiotic target genes. Through comparative metagenomic analysis of previously published AR gene abundance in three publicly available metagenomes, we illustrate how mutation-generated resistance genes are either falsely assigned or neglected, which alters the detection and quantitation of the antibiotic resistome. In addition, we inspected factors influencing the outcome of AR gene quantification using metagenome simulation experiments, and identified that genome size, AR gene length, total number of metagenomics reads and selected sequencing platforms had pronounced effects on the level of detected AR. In conclusion, our proposed improvements in the current methodologies for accurate AR detection and resistome assessment show reliable results when tested on real and simulated metagenomic datasets.

Antibiotic resistance genes and human bacterial pathogens: Co-occurrence, removal, and enrichment in municipal sewage sludge digesters

Understanding which/how antibiotic resistance genes (ARGs) contribute to increased acquisition of resistance by pathogens in aquatic environments are challenges of profound significance. We explored the co-occurrence and removal versus enrichment of ARGs and human bacterial pathogens (HBPs) in municipal sewage sludge digesters. We combined metagenomic detection of a wide spectrum of 323 ARGs and 83 HBPs with a correlation-based statistical approach and charted a network of their co-occurrence relationships. The results indicate that most ARGs and a minor proportion of HBPs (mainly Collinsella aerofaciens, Streptococcus salivarius and Gordonia bronchialis) could not be removed by anaerobic digestion, revealing a biological risk of post-digestion sludge in disseminating antibiotic resistance and pathogenicity. Moreover, preferential co-occurrence patterns were evident within one ARG type (e.g., multidrug, beta-lactam, and aminoglycoside) and between two different ARG types (i.e., aminoglycoside and beta-lactam), possibly implicating co-effects of antibiotic selection pressure and co-resistance on shaping antibiotic resistome in sewage sludge. Unlike beta-lactam resistance genes, ARGs of multidrug and macrolide-lincosamide-streptogramin tended to co-occur more with HBPs. Strikingly, we presented evidence that the most straightforward biological origin of an ARG-species co-occurring event is a hosting relationship. Furthermore, a significant and robust HBP-species co-occurrence correlation provides a proper scenario for nominating HBP indicators (e.g., Bifidobacterium spp. are perfect indicators of C. aerofaciens; r = 0.92-0.99 and P-values < 0.01). Combined, this study demonstrates a creative and effective network-based metagenomic approach for exploring ARG hosts and HBP indicators and assessing ARGs acquisition by HBPs in human-impacted environments where ARGs and HBPs may co-thrive.

Photochemical Transformation of Aminoglycoside Antibiotics in Simulated Natural Waters

Aminoglycoside antibiotics are widely used in human therapy and veterinary medicine. We report herein a detailed study on the natural-organic-matter- (NOM-) photosensitized degradation of aminoglycosides in aqueous media under simulated solar irradiation. It appears that the direct reaction of the excited states of NOM ((3)NOM*) with aminoglycosides is minor. The contributions of reactive oxygen species (ROSs) in the bulk solutions are also unimportant, as determined by an assessment based on steady-state concentrations and bimolecular reaction rate constants in a homogeneous reaction model. The inhibition of the photodegradation by isopropamide is rationalized through competitive sorption with aminoglycosides on the NOM surface, whereas the addition of isopropanol negligibly affects degradation because it quenches HO(•) in the bulk solution but not HO(•) localized on the NOM surface where aminoglycosides reside. Therefore, a sorption-enhanced phototransformation mechanism is proposed. The sorption of aminoglycosides on NOM follows a dual-mode model involving Langmuir and linear isotherms. The steady-state concentration of HO(•) on the surface of NOM was calculated as 10(-14) M, 2 orders of magnitude higher than that in the bulk solution. This fundamental information is important in the assessment of the fate and transport of aminoglycosides in aqueous environments.

ARGs-OAP: online analysis pipeline for antibiotic resistance genes detection from metagenomic data using an integrated structured ARG-database

Motivation: Environmental dissemination of antibiotic resistance genes (ARGs) has become an increasing concern for public health. Metagenomics approaches can effectively detect broad profiles of ARGs in environmental samples; however, the detection and subsequent classification of ARG-like sequences are time consuming and have been severe obstacles in employing metagenomic methods. We sought to accelerate quantification of ARGs in metagenomic data from environmental samples.

Results: A Structured ARG reference database (SARG) was constructed by integrating ARDB and CARD, the two most commonly used databases. SARG was curated to remove redundant sequences and optimized to facilitate query sequence identification by similarity. A database with a hierarchical structure (type-subtype-reference sequence) was then constructed to facilitate classification (assigning ARG-like sequence to type, subtype and reference sequence) of sequences identified through similarity search. Utilizing SARG and a previously proposed hybrid functional gene annotation pipeline, we developed an online pipeline called ARGs-OAP for fast annotation and classification of ARG-like sequences from metagenomic data. We also evaluated and proposed a set of criteria important for efficiently conducting metagenomic analysis of ARGs using ARGs-OAP.

Availability and Implementation: Perl script for ARGs-OAP can be downloaded from https://github.com/biofuture/Ublastx_stageone. ARGs-OAP can be accessed through http://smile.hku.hk/SARGs.

Contact:  zhangt@hku.hk or tiedjej@msu.edu

Supplementary information:  Supplementary data are available at Bioinformatics online.

Effects of antibiotic resistance genes on the performance and stability of different microbial aggregates in a granular sequencing batch reactor

Antibiotic resistance genes (ARGs) have emerged as key factors in wastewater environmental contaminants and continue to pose a challenge for wastewater treatment processes. With the aim of investigating the performance of granular sludge system when treating wastewater containing a considerable amount of ARGs, a lab-scale granular sequencing batch reactor (GSBR) where flocculent and granular sludge coexisted was designed. The results showed that after inoculation of donor strain NH4(+)-N purification efficiency diminished from 94.7% to 32.8% and recovered to 95.2% after 10 days. Meanwhile, RP4 plasmid had varying effects on different forms of microbial aggregates. As the size of aggregates increased, the abundance of RP4 in sludge decreased. The residence time of RP4 in granules with particle size exceeding 0.9 mm (14 days) was far shorter than that in flocculent sludge (26 days). Therefore, our studies conclude that with increasing number of ARGs being detected in wastewater, the use of granular sludge system in wastewater treatment processes will allow the reduction of ARGs transmissions and lessen potential ecological threats.

Metagenomic profiles of antibiotic resistance genes in paddy soils from South China

Overuse and arbitrary discarding of antibiotics have expanded antibiotic resistance reservoirs, from gut, waste water and activated sludge, to soil, freshwater and even the ocean. Based on the structured Antibiotic Resistance Genes Database and next generation sequencing, metagenomic analysis was used for the first time to detect and quantify antibiotic resistance genes (ARGs) in paddy soils from South China. A total of 16 types of ARGs were identified, corresponding to 110 ARG subtypes. The abundances and distribution pattern of ARGs in paddy soil were distinctively different from those in activated sludge and pristine deep ocean sediment, but close to those of sediment from human-impacted estuaries. Multidrug resistance genes were the most dominant type (38-47.5%) in all samples, and the ARGs detected encompassed the three major resistance mechanisms, among which extrusion by efflux pumps was predominant. Redundancy analysis (RDA) showed that pH was significantly correlated with the distribution of ARG subtypes (P < 0.05). Our results provided a broad spectrum profile of ARGs in paddy soil, indicating that ARGs are widespread in paddy soils of South China.

High-throughput quantification of antibiotic resistance genes from an urban wastewater treatment plant

Antibiotic resistance among bacteria is a growing problem worldwide, and wastewater treatment plants have been considered as one of the major contributors to the dissemination of antibiotic resistance to the environment. There is a lack of comprehensive quantitative molecular data on extensive numbers of antibiotic resistance genes (ARGs) in different seasons with a sampling strategy that would cover both incoming and outgoing water together with the excess sludge that is removed from the process. In order to fill that gap we present a highly parallel quantitative analysis of ARGs and horizontal gene transfer potential over four seasons at an urban wastewater treatment plant using a high-throughput qPCR array. All analysed transposases and two-thirds of primer sets targeting ARGs were detected in the wastewater. The relative abundance of most of the genes was highest in influent and lower in effluent water and sludge. The resistance profiles of the samples cluster by sample location with a shift from raw influent through the final effluents and dried sludge to the sediments. Wastewater discharge enriched only a few genes, namely Tn25 type transposase gene and clinical class 1 integrons, in the sediment near the discharge pipe, but those enriched genes may indicate a potential for horizontal gene transfer.

Metagenomic Assembly Reveals Hosts of Antibiotic Resistance Genes and the Shared Resistome in Pig, Chicken, and Human Feces

The risk associated with antibiotic resistance disseminating from animal and human feces is an urgent public issue. In the present study, we sought to establish a pipeline for annotating antibiotic resistance genes (ARGs) based on metagenomic assembly to investigate ARGs and their co-occurrence with associated genetic elements. Genetic elements found on the assembled genomic fragments include mobile genetic elements (MGEs) and metal resistance genes (MRGs). We then explored the hosts of these resistance genes and the shared resistome of pig, chicken and human fecal samples. High levels of tetracycline, multidrug, erythromycin, and aminoglycoside resistance genes were discovered in these fecal samples. In particular, significantly high level of ARGs (7762 ×/Gb) was detected in adult chicken feces, indicating higher ARG contamination level than other fecal samples. Many ARGs arrangements (e.g., macA-macB and tetA-tetR) were discovered shared by chicken, pig and human feces. In addition, MGEs such as the aadA5-dfrA17-carrying class 1 integron were identified on an assembled scaffold of chicken feces, and are carried by human pathogens. Differential coverage binning analysis revealed significant ARG enrichment in adult chicken feces. A draft genome, annotated as multidrug resistant Escherichia coli, was retrieved from chicken feces metagenomes and was determined to carry diverse ARGs (multidrug, acriflavine, and macrolide). The present study demonstrates the determination of ARG hosts and the shared resistome from metagenomic data sets and successfully establishes the relationship between ARGs, hosts, and environments. This ARG annotation pipeline based on metagenomic assembly will help to bridge the knowledge gaps regarding ARG-associated genes and ARG hosts with metagenomic data sets. Moreover, this pipeline will facilitate the evaluation of environmental risks in the genetic context of ARGs.

Prevalence and proliferation of antibiotic resistance genes in two municipal wastewater treatment plants

The propagation of antibiotic resistance genes (ARGs) is an emerging health concern worldwide. Thus, it is important to understand and mitigate their occurrence in different systems. In this study, 30 ARGs that confer resistance to tetracyclines, sulfonamides, quinolones or macrolides were detected in two activated sludge wastewater treatment plants (WWTPs) in northern China. Bacteria harboring ARGs persisted through all treatment units, and survived disinfection by chlorination in greater percentages than total Bacteria (assessed by 16S rRNA genes). Although the absolute abundances of ARGs were reduced from the raw influent to the effluent by 89.0%-99.8%, considerable ARG levels [(1.0 ± 0.2) × 10(3) to (9.5 ± 1.8) × 10(5) copies/mL)] were found in WWTP effluent samples. ARGs were concentrated in the waste sludge (through settling of bacteria and sludge dewatering) at (1.5 ± 2.3) × 10(9) to (2.2 ± 2.8) × 10(11) copies/g dry weight. Twelve ARGs (tetA, tetB, tetE, tetG, tetH, tetS, tetT, tetX, sul1, sul2, qnrB, ermC) were discharged through the dewatered sludge and plant effluent at higher rates than influent values, indicating overall proliferation of resistant bacteria. Significant antibiotic concentrations (2%-50% of raw influent concentrations) remained throughout all treatment units. This apparently contributed selective pressure for ARG replication since the relative abundance of resistant bacteria (assessed by ARG/16S rRNA gene ratios) was significantly correlated to the corresponding effluent antibiotic concentrations. Similarly, the concentrations of various heavy metals (which induce a similar bacterial resistance mechanism as antibiotics - efflux pumps) were also correlated to the enrichment of some ARGs. Thus, curtailing the release of antibiotics and heavy metals to sewage systems (or enhancing their removal in pre-treatment units) may alleviate their selective pressure and mitigate ARG proliferation in WWTPs.

Experimental Design and Bioinformatics Analysis for the Application of Metagenomics in Environmental Sciences and Biotechnology

Recent advances in DNA sequencing technologies have prompted the widespread application of metagenomics for the investigation of novel bioresources (e.g., industrial enzymes and bioactive molecules) and unknown biohazards (e.g., pathogens and antibiotic resistance genes) in natural and engineered microbial systems across multiple disciplines. This review discusses the rigorous experimental design and sample preparation in the context of applying metagenomics in environmental sciences and biotechnology. Moreover, this review summarizes the principles, methodologies, and state-of-the-art bioinformatics procedures, tools and database resources for metagenomics applications and discusses two popular strategies (analysis of unassembled reads versus assembled contigs/draft genomes) for quantitative or qualitative insights of microbial community structure and functions. Overall, this review aims to facilitate more extensive application of metagenomics in the investigation of uncultured microorganisms, novel enzymes, microbe-environment interactions, and biohazards in biotechnological applications where microbial communities are engineered for bioenergy production, wastewater treatment, and bioremediation.

Development of antibiotic resistance genes in microbial communities during long-term operation of anaerobic reactors in the treatment of pharmaceutical wastewater

Biological treatment processes offer the ideal conditions in which a high diversity of microorganisms can grow and develop. The wastewater produced during these processes is contaminated with antibiotics and, as such, they provide the ideal setting for the acquisition and proliferation of antibiotic resistance genes (ARGs). This research investigated the occurrence and variation in the ARGs found during the one-year operation of the anaerobic sequencing batch reactors (SBRs) used to treat pharmaceutical wastewater that contained combinations of sulfamethoxazole-tetracycline-erythromycin (STE) and sulfamethoxazole-tetracycline (ST). The existence of eighteen ARGs encoding resistance to sulfamethoxazole (sul1, sul2, sul3), erythromycin (ermA, ermF, ermB, msrA, ereA), tetracycline (tetA, tetB, tetC, tetD, tetE, tetM, tetS, tetQ, tetW, tetX) and class Ι integron gene (intΙ 1) in the STE and ST reactors was investigated by quantitative real-time PCR. Due to the limited availability of primers to detect ARGs, Illumina sequencing was also performed on the sludge and effluent of the STE and ST reactors. Although there was good reactor performance in the SBRs, which corresponds to min 80% COD removal efficiency, tetA, tetB, sul1, sul2 and ermB genes were among those ARGs detected in the effluent from STE and ST reactors. A comparison of the ARGs acquired from the STE and ST reactors revealed that the effluent from the STE reactor had a higher number of ARGs than that from the ST reactor; this could be due to the synergistic effects of erythromycin. According to the expression of genes results, microorganisms achieve tetracycline and erythromycin resistance through a combination of three mechanisms: efflux pumping protein, modification of the antibiotic target and modifying enzymes. There was also a significant association between the presence of the class 1 integron and sulfamethoxazole resistance genes.

Limited dissemination of the wastewater treatment plant core resistome

Horizontal gene transfer is a major contributor to the evolution of bacterial genomes and can facilitate the dissemination of antibiotic resistance genes between environmental reservoirs and potential pathogens. Wastewater treatment plants (WWTPs) are believed to play a central role in the dissemination of antibiotic resistance genes. However, the contribution of the dominant members of the WWTP resistome to resistance in human pathogens remains poorly understood. Here we use a combination of metagenomic functional selections and comprehensive metagenomic sequencing to uncover the dominant genes of the WWTP resistome. We find that this core resistome is unique to the WWTP environment, with <10% of the resistance genes found outside the WWTP environment. Our data highlight that, despite an abundance of functional resistance genes within WWTPs, only few genes are found in other environments, suggesting that the overall dissemination of the WWTP resistome is comparable to that of the soil resistome.

Exploring antibiotic resistance genes and metal resistance genes in plasmid metagenomes from wastewater treatment plants

Plasmids operate as independent genetic elements in microorganism communities. Through horizontal gene transfer (HGT), they can provide their host microorganisms with important functions such as antibiotic resistance and heavy metal resistance. In this study, six metagenomic libraries were constructed with plasmid DNA extracted from influent, activated sludge (AS) and digested sludge (DS) of two wastewater treatment plants (WWTPs). Compared with the metagenomes of the total DNA extracted from the same sectors of the wastewater treatment plant, the plasmid metagenomes had significantly higher annotation rates, indicating that the functional genes on plasmids are commonly shared by those studied microorganisms. Meanwhile, the plasmid metagenomes also encoded many more genes related to defense mechanisms, including ARGs. Searching against an antibiotic resistance genes (ARGs) database and a metal resistance genes (MRGs) database revealed a broad-spectrum of antibiotic (323 out of a total 618 subtypes) and MRGs (23 out of a total 23 types) on these plasmid metagenomes. The influent plasmid metagenomes contained many more resistance genes (both ARGs and MRGs) than the AS and the DS metagenomes. Sixteen novel plasmids with a complete circular structure that carried these resistance genes were assembled from the plasmid metagenomes. The results of this study demonstrated that the plasmids in WWTPs could be important reservoirs for resistance genes, and may play a significant role in the horizontal transfer of these genes.

Profile and Fate of Bacterial Pathogens in Sewage Treatment Plants Revealed by High-Throughput Metagenomic Approach

The broad-spectrum profile of bacterial pathogens and their fate in sewage treatment plants (STPs) were investigated using high-throughput sequencing based metagenomic approach. This novel approach could provide a united platform to standardize bacterial pathogen detection and realize direct comparison among different samples. Totally, 113 bacterial pathogen species were detected in eight samples including influent, effluent, activated sludge (AS), biofilm, and anaerobic digestion sludge with the abundances ranging from 0.000095% to 4.89%. Among these 113 bacterial pathogens, 79 species were reported in STPs for the first time. Specially, compared to AS in bulk mixed liquor, more pathogen species and higher total abundance were detected in upper foaming layer of AS. This suggests that the foaming layer of AS might impose more threat to onsite workers and citizens in the surrounding areas of STPs because pathogens in foaming layer are easily transferred into air and cause possible infections. The high removal efficiency (98.0%) of total bacterial pathogens suggests that AS treatment process is effective to remove most bacterial pathogens. Remarkable similarities of bacterial pathogen compositions between influent and human gut indicated that bacterial pathogen profiles in influents could well reflect the average bacterial pathogen communities of urban resident guts within the STP catchment area.

The horizontal transfer of antibiotic resistance genes is enhanced by ionic liquid with different structure of varying alkyl chain length

Antibiotic resistance genes (ARGs) have become a global health concern. In our previous study, an ionic liquid (IL) 1-butyl-3-methylimidazolium hexafluorophosphate ([BMIm][PF6]) had been proven to facilitate the dissemination of ARGs in the environment. However, enhanced alkyl group chain length or the substitution of alkyl groups with the cation ring corresponded with increased antimicrobial effects. In this study, we investigated how different structures of ILs with 4, 6, and 8 C atoms in the longer alkyl chain on the imidazolium cations facilitated the dissemination of ARGs. The promotion of plasmid RP4 transfer frequency decreased with [CnMIM][BF4] increasing the alkyl chain length from 4 carbon atoms to 8 carbon atoms on the imidazolium cations, which is observed with [BMIM][BF4] (n = 4, 5.9 fold) > HMIM][BF4] (n = 6, 2.2 fold) > [OMIM][BF4] (n = 8, 1.7 fold). This illustrates that [CnMIM][BF4] with increasing the alkyl chain length exert decreasing ability in facilitating plasmid RP4 horizontal transfer, which is possibly related to IL-structure dependent toxicity. The IL-structure dependent plasmid RP4 transfer frequency was attributable to bacterial cell membrane permeability weaken with increasing alkyl chain length of [CnMIM][PF4], which was evidenced by flow cytometry. In freshwater microcosm, [CnMIm][BF4] promoted the relative abundance of the sulI and intI genes for 4.6 folds, aphA and traF for 5.2 folds higher than the untreated groups, promoting the propagation of ARGs in the aquatic environment. This is the first report that ILs with different structure of varying alkyl chain length facilitate horizontal transfer of plasmid RP4 which is widely distributed in the environment, and thus add the adverse effects of the environmental risk of ILs.

Metagenomic Analysis of Antibiotic Resistance Genes in Dairy Cow Feces following Therapeutic Administration of Third Generation Cephalosporin

Although dairy manure is widely applied to land, it is relatively understudied compared to other livestock as a potential source of antibiotic resistance genes (ARGs) to the environment and ultimately to human pathogens. Ceftiofur, the most widely used antibiotic used in U.S. dairy cows, is a 3rd generation cephalosporin, a critically important class of antibiotics to human health. The objective of this study was to evaluate the effect of typical ceftiofur antibiotic treatment on the prevalence of ARGs in the fecal microbiome of dairy cows using a metagenomics approach. β-lactam ARGs were found to be elevated in feces from Holstein cows administered ceftiofur (n = 3) relative to control cows (n = 3). However, total numbers of ARGs across all classes were not measurably affected by ceftiofur treatment, likely because of dominance of unaffected tetracycline ARGs in the metagenomics libraries. Functional analysis via MG-RAST further revealed that ceftiofur treatment resulted in increases in gene sequences associated with "phages, prophages, transposable elements, and plasmids", suggesting that this treatment also enriched the ability to horizontally transfer ARGs. Additional functional shifts were noted with ceftiofur treatment (e.g., increase in genes associated with stress, chemotaxis, and resistance to toxic compounds; decrease in genes associated with metabolism of aromatic compounds and cell division and cell cycle), along with measureable taxonomic shifts (increase in Bacterioidia and decrease in Actinobacteria). This study demonstrates that ceftiofur has a broad, measureable and immediate effect on the cow fecal metagenome. Given the importance of 3rd generation cephalospirins to human medicine, their continued use in dairy cattle should be carefully considered and waste treatment strategies to slow ARG dissemination from dairy cattle manure should be explored.

Antibiotic resistome and its association with bacterial communities during sewage sludge composting

Composting is widely used for recycling of urban sewage sludge to improve soil properties, which represents a potential pathway of spreading antibiotic resistant bacteria and genes to soils. However, the dynamics of antibiotic resistance genes (ARGs) and the underlying mechanisms during sewage sludge composting were not fully explored. Here, we used high-throughput quantitative PCR and 16S rRNA gene based illumina sequencing to investigate the dynamics of ARGs and bacterial communities during a lab-scale in-vessel composting of sewage sludge. A total of 156 unique ARGs and mobile genetic elements (MGEs) were detected encoding resistance to almost all major classes of antibiotics. ARGs were detected with significantly increased abundance and diversity, and distinct patterns, and were enriched during composting. Marked shifts in bacterial community structures and compositions were observed during composting, with Actinobacteria being the dominant phylum at the late phase of composting. The large proportion of Actinobacteria may partially explain the increase of ARGs during composting. ARGs patterns were significantly correlated with bacterial community structures, suggesting that the dynamic of ARGs was strongly affected by bacterial phylogenetic compositions during composting. These results imply that direct application of sewage sludge compost on field may lead to the spread of abundant ARGs in soils.

Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes

A metagenomic approach and network analysis was used to investigate the wide-spectrum profiles of antibiotic resistance genes (ARGs) and their co-occurrence patterns in 50 samples from 10 typical environments. In total, 260 ARG subtypes belonging to 18 ARG types were detected with an abundance range of 5.4 × 10(-6)-2.2 × 10(-1) copy of ARG per copy of 16S-rRNA gene. The trend of the total ARG abundances in environments matched well with the levels of anthropogenic impacts on these environments. From the less impacted environments to the seriously impacted environments, the total ARG abundances increased up to three orders of magnitude, that is, from 3.2 × 10(-3) to 3.1 × 10(0) copy of ARG per copy of 16S-rRNA gene. The abundant ARGs were associated with aminoglycoside, bacitracin, β-lactam, chloramphenicol, macrolide-lincosamide-streptogramin, quinolone, sulphonamide and tetracycline, in agreement with the antibiotics extensively used in human medicine or veterinary medicine/promoters. The widespread occurrences and abundance variation trend of vancomycin resistance genes in different environments might imply the spread of vancomycin resistance genes because of the selective pressure resulting from vancomycin use. The simultaneous enrichment of 12 ARG types in adult chicken faeces suggests the coselection of multiple ARGs in this production system. Non-metric multidimensional scaling analysis revealed that samples belonging to the same environment generally possessed similar ARG compositions. Based on the co-occurrence pattern revealed by network analysis, tetM and aminoglycoside resistance protein, the hubs of the ARG network, are proposed to be indicators to quantitatively estimate the abundance of 23 other co-occurring ARG subtypes by power functions.

Variation of antibiotic resistance genes in municipal wastewater treatment plant with A(2)O-MBR system

The variation of five antibiotic resistance genes (ARGs)-tetG, tetW, tetX, sul1, and intI1-in a full-scale municipal wastewater treatment plant with A(2)O-MBR system was studied. The concentrations of five resistance genes both in influent and in membrane bioreactor (MBR) effluent decreased as sul1 > intI1 > tetX > tetG > tetW, and an abundance of sul1 was statistically higher than three other tetracycline resistance genes (tetG, tetW, and tetX) (p < 0.05). The concentrations of five ARGs in the influent were all higher in spring (median 10(5.81)-10(7.32) copies mL(-1)) than they were in other seasons, and tetW, tetX, and sul1 reached its lowest concentration in autumn (10(4.61)-10(6.81) copies mL(-1)). The concentration of ARGs in wastewater decreased in the anaerobic effluent and anoxic effluent, but increased in the aerobic effluent, and then sharply declined in the MBR effluent. The reduction of tetW, intI1, and sul1 was all significantly positively correlated with the reduction of 16S ribosomal DNA (rDNA) in the wastewater treatment process (p < 0.01). The concentration of ARGs (copies mg(-1)) in sludge samples increased along the treatment process, but the abundance of five ARGs (ratio of ARGs to 16S rDNA) remained the same from anaerobic to anoxic to aerobic basins, while an increment ratio in MBR was observed for all ARGs.

Metagenomics shows that low-energy anaerobic-aerobic treatment reactors reduce antibiotic resistance gene levels from domestic wastewater

Effective domestic wastewater treatment is among our primary defenses against the dissemination of infectious waterborne disease. However, reducing the amount of energy used in treatment processes has become essential for the future. One low-energy treatment option is anaerobic-aerobic sequence (AAS) bioreactors, which use an anaerobic pretreatment step (e.g., anaerobic hybrid reactors) to reduce carbon levels, followed by some form of aerobic treatment. Although AAS is common in warm climates, it is not known how its compares to other treatment options relative to disease transmission, including its influence on antibiotic resistance (AR) in treated effluents. Here, we used metagenomic approaches to contrast the fate of antibiotic-resistant genes (ARG) in anaerobic, aerobic, and AAS bioreactors treating domestic wastewater. Five reactor configurations were monitored for 6 months, and treatment performance, energy use, and ARG abundance and diversity were compared in influents and effluents. AAS and aerobic reactors were superior to anaerobic units in reducing ARG-like sequence abundances, with effluent ARG levels of 29, 34, and 74 ppm (198 ppm influent), respectively. AAS and aerobic systems especially reduced aminoglycoside, tetracycline, and β-lactam ARG levels relative to anaerobic units, although 63 persistent ARG subtypes were detected in effluents from all systems (of 234 assessed). Sulfonamide and chloramphenicol ARG levels were largely unaffected by treatment, whereas a broad shift from target-specific ARGs to ARGs associated with multi-drug resistance was seen across influents and effluents. AAS reactors show promise for future applications because they can reduce more ARGs for less energy (32% less energy here), but all three treatment options have limitations and need further study.

Investigation of fate and behavior of tetracycline in nitrifying sludge system

This study aims to investigate the fate and behavior of tetracycline (TC) in nitrifying sludge system, as well as the effects of TC dosage on sludge performance.

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.

Seasonal variability in airborne biotic contaminants in swine confinement buildings

Little is known about the seasonal dynamics of biotic contaminants in swine confinement buildings (SCBs). The biotic contaminants of seven SCBs were monitored during one visit in the winter and one during the summer. Paired-end Illumina sequencing of the 16S rRNA gene, V3 region, was used to examine seasonal shifts in bacterial community composition and diversity. The abundances of 16S rRNA genes and six tetracycline resistance genes (tetB, tetH, tetZ, tetO, tetQ, and tetW) were also quantified using real-time PCR. Bacterial abundances, community composition and diversity all showed strong seasonal patterns defined by winter peaks in abundance and diversity. Microclimatic variables of SCBs, particularly air speed, PM2.5 and total suspended particles (TSP) were found significantly correlated to abundances, community composition, and diversity of bacterial bioaerosols. Seasonal fluctuations were also observed for four tetracycline resistance genes, tetH, tetO, tetQ, and tetW. The frequency of occurrences of these resistance genes were significantly higher in samples collected during winter and was also significantly correlated with air speed, PM2.5 and TSP. Overall, our results indicate that biotic contaminants in SCBs exhibit seasonal trends, and these could be associated with the microclimatic variables of SCBs. The correlations established in the current study could be helpful in establishing better management strategies to minimize the potential health impacts on both livestock and humans working in this environment.

Evaluation of a hybrid approach using UBLAST and BLASTX for metagenomic sequences annotation of specific functional genes

The fast development of next generation sequencing (NGS) has dramatically increased the application of metagenomics in various aspects. Functional annotation is a major step in the metagenomics studies. Fast annotation of functional genes has been a challenge because of the deluge of NGS data and expanding databases. A hybrid annotation pipeline proposed previously for taxonomic assignments was evaluated in this study for metagenomic sequences annotation of specific functional genes, such as antibiotic resistance genes, arsenic resistance genes and key genes in nitrogen metabolism. The hybrid approach using UBLAST and BLASTX is 44–177 times faster than direct BLASTX in the annotation using the small protein database for the specific functional genes, with the cost of missing a small portion (<1.8%) of target sequences compared with direct BLASTX hits. Different from direct BLASTX, the time required for specific functional genes annotation using the hybrid annotation pipeline depends on the abundance for the target genes. Thus this hybrid annotation pipeline is more suitable in specific functional genes annotation than in comprehensive functional genes annotation.

Metagenomic approach reveals variation of microbes with arsenic and antimony metabolism genes from highly contaminated soil

Microbes have great potential for arsenic (As) and antimony (Sb) bioremediation in heavily contaminated soil because they have the ability to biotransform As and Sb to species that have less toxicity or are more easily removed. In this study, we integrated a metagenomic method with physicochemical characterization to elucidate the composition of microbial community and functional genes (related to As and Sb) in a high As (range from 34.11 to 821.23 mg kg⁻¹) and Sb (range from 226.67 to 3923.07 mg kg⁻¹) contaminated mine field. Metagenomic analysis revealed that microbes from 18 phyla were present in the 5 samples of soil contaminated with high As and Sb. Moreover, redundancy analysis (RDA) of the relationship between the 18 phyla and the concentration of As and Sb demonstrated that 5 phyla of microbes, i.e. Actinobacteria, Firmicutes, Nitrospirae, Tenericutes and Gemmatimonadetes were positively correlated with As and Sb concentration. The distribution, diversity and abundance of functional genes (including arsC, arrA, aioA, arsB and ACR3) were much higher for the samples containing higher As and Sb concentrations. Based on correlation analysis, the results showed a positive relationship between arsC-like (R² = 0.871) and aioA-like (R² = 0.675) gene abundance and As concentration, and indicated that intracellular As(V) reduction and As(III) oxidation could be the dominant As detoxification mechanism enabling the microbes to survive in the environment. This study provides a direct and reliable reference on the diversity of microbial community and functional genes in an extremely high concentration As- and Sb-contaminated environment.

Fate of antibiotic resistance genes in sewage treatment plant revealed by metagenomic approach

Antibiotic resistance has become a serious threat to human health. Sewage treatment plant (STP) is one of the major sources of antibiotic resistance genes (ARGs) in natural environment. High-throughput sequencing-based metagenomic approach was applied to investigate the broad-spectrum profiles and fate of ARGs in a full scale STP. Totally, 271 ARGs subtypes belonging to 18 ARGs types were identified by the broad scanning of metagenomic analysis. Influent had the highest ARGs abundance, followed by effluent, anaerobic digestion sludge and activated sludge. 78 ARGs subtypes persisted through the biological wastewater and sludge treatment process. The high removal efficiency of 99.82% for total ARGs in wastewater suggested that sewage treatment process is effective in reducing ARGs. But the removal efficiency of ARGs in sludge treatment was not as good as that in sewage treatment. Furthermore, the composition of microbial communities was examined and the correlation between microbial community and ARGs was investigated using redundancy analysis. Significant correlation between 6 genera and the distribution of ARGs were found and 5 of the 6 genera included potential pathogens. This is the first study on the fate of ARGs in STP using metagenomic analysis with high-throughput sequencing and hopefully would enhance our knowledge on fate of ARGs in STP.

High throughput profiling of antibiotic resistance genes in urban park soils with reclaimed water irrigation

Reclaimed water irrigation (RWI) in urban environments is becoming popular, due to rapid urbanization and water shortage. The continuous release of residual antibiotics and antibiotic resistance genes (ARGs) from reclaimed water could result in the dissemination of ARGs in the downstream environment. This study provides a comprehensive profile of ARGs in park soils exposed to RWI through a high-throughput quantitative PCR approach. 147 ARGs encoding for resistance to a broad-spectrum of antibiotics were detected among all park soil samples. Aminoglycoside and beta-lactam were the two most dominant types of ARGs, and antibiotic deactivation and efflux pump were the two most dominant mechanisms in these RWI samples. The total enrichment of ARGs varied from 99.3-fold to 8655.3-fold compared to respective controls. Six to 60 ARGs were statistically enriched among these RWI samples. Four transposase genes were detected in RWI samples. TnpA-04 was the most enriched transposase gene with an enrichment was up to 2501.3-fold in Urumqi RWI samples compared with control soil samples. Furthermore, significantly positive correlation was found between ARGs and transposase abundances, indicating that transposase might be involved in the propagation of ARGs. This study demonstrated that RWI resulted in the enrichment of ARGs in urban park soils.

Potentially novel copper resistance genes in copper-enriched activated sludge revealed by metagenomic analysis

In this study, we utilized the Illumina high-throughput metagenomic approach to investigate diversity and abundance of both microbial community and copper resistance genes (CuRGs) in activated sludge (AS) which was enriched under copper selective stress up to 800 mg/L. The raw datasets (~3.5 Gb for each sample, i.e., the copper-enriched AS and the control AS) were merged and normalized for the BLAST analyses against the SILVA SSU rRNA gene database and self-constructed copper resistance protein database (CuRD). Also, the raw metagenomic sequences were assembled into contigs and analyzed based on Open Reading Frames (ORFs) to identify potentially novel copper resistance genes. Among the different resistance systems for copper detoxification under the high copper stress condition, the Cus system was the most enriched system. The results also indicated that genes encoding multi-copper oxidase played a more important role than those encoding efflux proteins. More significantly, several potentially novel copper resistance ORFs were identified by Pfam search and phylogenic analysis. This study demonstrated a new understanding of microbial-mediated copper resistance under high copper stress using high-throughput shotgun sequencing technique.