In the fight to limit the global spread of antibiotic resistance, the assembly of environmental metagenomes has the potential to provide rich contextual information (e.g., taxonomic hosts, carriage on mobile genetic elements) about antibiotic resistance genes (ARG) in the environment. However, computational challenges associated with assembly can impact the accuracy of downstream analyses. This work critically evaluates the impact of assembly leveraging short reads, nanopore MinION long-reads, and a combination of the two (hybrid) on ARG contextualization for ten environmental metagenomes using seven prominent assemblers (IDBA-UD, MEGAHIT, Canu, Flye, Opera-MS, metaSpades and HybridSpades). While short-read and hybrid assemblies produced similar patterns of ARG contextualization, raw or assembled long nanopore reads produced distinct patterns. Based on an in-silico spike-in experiment using real and simulated reads, we show that low to intermediate coverage species are more likely to be incorporated into chimeric contigs across all assemblers and sequencing technologies, while more abundant species produce assemblies with a greater frequency of inversions and insertion/deletions (indels). In sum, our analyses support hybrid assembly as a valuable technique for boosting the reliability and accuracy of assembly-based analyses of ARGs and neighboring genes at environmentally-relevant coverages, provided that sufficient short-read sequencing depth is achieved.
- Award ID(s):
- 2013998
- NSF-PAR ID:
- 10389123
- Date Published:
- Journal Name:
- Frontiers in Bioengineering and Biotechnology
- Volume:
- 10
- ISSN:
- 2296-4185
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Background There is concern that the microbially rich activated sludge environment of wastewater treatment plants (WWTPs) may contribute to the dissemination of antibiotic resistance genes (ARGs). We applied long-read (nanopore) sequencing to profile ARGs and their neighboring genes to illuminate their fate in the activated sludge treatment by comparing their abundance, genetic locations, mobility potential, and bacterial hosts within activated sludge relative to those in influent sewage across five WWTPs from three continents.
Results The abundances (gene copies per Gb of reads, aka gc/Gb) of all ARGs and those carried by putative pathogens decreased 75–90% from influent sewage (192-605 gc/Gb) to activated sludge (31-62 gc/Gb) at all five WWTPs. Long reads enabled quantification of the percent abundance of ARGs with mobility potential (i.e., located on plasmids or co-located with other mobile genetic elements (MGEs)). The abundance of plasmid-associated ARGs decreased at four of five WWTPs (from 40–73 to 31–68%), and ARGs co-located with transposable, integrative, and conjugative element hallmark genes showed similar trends. Most ARG-associated elements decreased 0.35–13.52% while integrative and transposable elements displayed slight increases at two WWTPs (1.4–2.4%). While resistome and taxonomic compositions both shifted significantly, host phyla for chromosomal ARG classes remained relatively consistent, indicating vertical gene transfer via active biomass growth in activated sludge as the key pathway of chromosomal ARG dissemination.
Conclusions Overall, our results suggest that the activated sludge process acted as a barrier against the proliferation of most ARGs, while those that persisted or increased warrant further attention.
-
more » « less
Characterization of antibiotic resistance genes (ARGs) from high-throughput sequencing data of metagenomics and cultured bacterial samples is a challenging task, with the need to account for both computational (e.g., string algorithms) and biological (e.g., gene transfers, rearrangements) aspects. Curated ARG databases exist together with assorted ARG classification approaches (e.g., database alignment, machine learning). Besides ARGs that naturally occur in bacterial strains or are acquired through mobile elements, there are chromosomal genes that can render a bacterium resistant to antibiotics through point mutations, i.e., ARG variants (ARGVs). While ARG repositories also collect ARGVs, there are only a few tools that are able to identify ARGVs from metagenomics and high throughput sequencing data, with a number of limitations (e.g., pre-assembly,
a posteriori verification of mutations, or specification of species). In this work we present thek -mer, i.e., strings of fixed lengthk , ARGV analyzer – KARGVA – an open-source, multi-platform tool that provides: (i) anad hoc , large ARGV database derived from multiple sources; (ii) input capability for various types of high-throughput sequencing data; (iii) a three-way, hash-based,k -mer search setup to process data efficiently, linkingk -mers to ARGVs,k -mers to point mutations, and ARGVs tok -mers, respectively; (iv) a statistical filter on sequence classification to reduce type I and II errors. On semi-synthetic data, KARGVA provides very high accuracy even in presence of high sequencing errors or mutations (99.2 and 86.6% accuracy within 1 and 5% base change rates, respectively), and genome rearrangements (98.2% accuracy), with robust performance onad hoc false positive sets. On data from the worldwide MetaSUB consortium, comprising 3,700+ metagenomics experiments, KARGVA identifies more ARGVs than Resistance Gene Identifier (4.8x) and PointFinder (6.8x), yet all predictions are below the expected false positive estimates. The prevalence of ARGVs is correlated to ARGs but ecological characteristics do not explain well ARGV variance. KARGVA is publicly available athttps://github.com/DataIntellSystLab/KARGVA under MIT license. -
more » « less
Abstract The amount of data produced by genome sequencing experiments has been growing rapidly over the past several years, making compression important for efficient storage, transfer and analysis of the data. In recent years, nanopore sequencing technologies have seen increasing adoption since they are portable, real-time and provide long reads. However, there has been limited progress on compression of nanopore sequencing reads obtained in FASTQ files since most existing tools are either general-purpose or specialized for short read data. We present NanoSpring, a reference-free compressor for nanopore sequencing reads, relying on an approximate assembly approach. We evaluate NanoSpring on a variety of datasets including bacterial, metagenomic, plant, animal, and human whole genome data. For recently basecalled high quality nanopore datasets, NanoSpring, which focuses only on the base sequences in the FASTQ file, uses just 0.35–0.65 bits per base which is 3–6
$$\times$$ $$\times$$ https://github.com/qm2/NanoSpring . -
null (Ed.)Antimicrobial resistance is a well-documented public health concern. The role that drinking water distribution pipes have as sources of antibiotic resistance genes (ARGs) is not well known. Metals are a known stressor for antibiotic resistance development, implying that aging metal-pipe infrastructure could be a source of ARGs. The objective of this study was to determine if ARGs, metal resistance genes (MRGs), and intI 1 were pervasive across various pipe biofilm sample types (biomass surfaces, pipe surfaces, corrosion tubercles, and under corrosion tubercles) and if the resistance genes associated with particular microbial taxa. Eight sample types in triplicate ( n = 24) were taken from inside a >100 year-old, six ft. section of a full-scale chloraminated cast iron drinking water main. Droplet digital PCR (ddPCR) was employed as a novel approach to quantify ARGs in pipes from full-scale drinking water distribution systems (DWDS) because it yielded higher detection frequencies than quantitative PCR (qPCR). Illumina sequencing was employed to characterize the microbial community based on 16S rRNA genes. ARGs and MRGs were detected in all 24 pipe samples. Every sample contained targeted genes. Interestingly, the mean absolute abundances of ARGs and MRGs only varied by approximately one log value across sample types, but the mean relative abundances (copy numbers normalized to 16S rRNA genes) varied by over two log values. The ARG and MRGs concentrations were not significantly different between sample types, despite significant changes in dominant microbial taxa. The most abundant genera observed in the biofilm communities were Mycobacterium (0.2–70%), and β-lactam resistance genes bla TEM , bla SHV , and the integrase gene of class 1 integrons ( intI 1) were positively correlated with Mycobacterium . The detection of ARGs, MRGs, and class 1 integrons across all sample types within the pipe indicates that pipes themselves can serve as sources for ARGs in DWDS. Consequently, future work should investigate the role of pipe materials as well as corrosion inhibitors to determine how engineering decisions can mitigate ARGs in drinking water that stem from pipe materials.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/fbioe.2022.1016408
