Plasmid metagenome reveals high levels of antibiotic resistance genes and mobile genetic elements in activated sludge

The overuse or misuse of antibiotics has accelerated antibiotic resistance, creating a major challenge for the public health in the world. Sewage treatment plants (STPs) are considered as important reservoirs for antibiotic resistance genes (ARGs) and activated sludge characterized with high microbial density and diversity facilitates ARG horizontal gene transfer (HGT) via mobile genetic elements (MGEs). However, little is known regarding the pool of ARGs and MGEs in sludge microbiome. In this study, the transposon aided capture (TRACA) system was employed to isolate novel plasmids from activated sludge of one STP in Hong Kong, China. We also used Illumina Hiseq 2000 high-throughput sequencing and metagenomics analysis to investigate the plasmid metagenome. Two novel plasmids were acquired from the sludge microbiome by using TRACA system and one novel plasmid was identified through metagenomics analysis. Our results revealed high levels of various ARGs as well as MGEs for HGT, including integrons, transposons and plasmids. The application of the TRACA system to isolate novel plasmids from the environmental metagenome, coupled with subsequent high-throughput sequencing and metagenomic analysis, highlighted the prevalence of ARGs and MGEs in microbial community of STPs.     

Metagenomic analysis reveals the prevalence and persistence of antibiotic- and heavy metal-resistance genes in wastewater treatment plant

The increased antibiotic resistance among microorganisms has resulted into growing interest for investigating the wastewater treatment plants (WWTPs) as they are reported to be the major source in the dissemination of antibiotic resistance genes (ARGs) and heavy metal resistance genes (HMRGs) in the environment. In this study, we investigated the prevalence and persistence of ARGs and HMRGs as well as bacterial diversity and mobile genetic elements (MGEs) in influent and effluent at the WWTP in Gwangju, South Korea, using high-throughput sequencing based metagenomic approach. A good number of broad-spectrum of resistance genes (both ARG and HMRG) were prevalent and likely persistent, although large portion of them were successfully removed at the wastewater treatment process. The relative abundance of ARGs and MGEs was higher in effluent as compared to that of influent. Our results suggest that the resistance genes with high abundance and bacteria harbouring ARGs and MGEs are likely to persist more through the treatment process. On analyzing the microbial community, the phylum Proteobacteria, especially potentially pathogenic species belonging to the genus Acinetobacter, dominated in WWTP. Overall, our study demonstrates that many ARGs and HMRGs may persist the treatment processes in WWTPs and their association to MGEs may contribute to the dissemination of resistance genes among microorganisms in the environment.     

High-Temperature-Induced Transposition of Insertion Elements in Burkholderia multivorans ATCC 17616

An efficient and quantitative method to analyze the transposition of various insertion sequence (IS) elements in Burkholderia multivorans ATCC 17616 was devised. pGEN500, a plasmid carrying a Bacillus subtilis-derived sacB gene, was introduced into ATCC 17616 cells, and 25% of their sucrose-resistant derivatives were found to carry various IS elements on pGEN500. A PCR-based experimental protocol, in which a mixture of several specific primer pairs was used, revealed that pGEN500 captured, in addition to five previously reported IS elements (IS401, IS402, IS406, IS407, and IS408), three novel IS elements, ISBmu1, ISBmu2, and ISBmu3. The global transposition frequency of these IS elements was enhanced more than sevenfold under a high-temperature condition (42°C) but not under oxidative stress or starvation conditions. To our knowledge, this is the first report demonstrating the elevated transposition activities of several IS elements at a high temperature. The efficient experimental protocol developed in this study will be useful in quantitatively and simultaneously investigating various IS elements, as well as in capturing novel functional mobile elements from a wide variety of bacteria.     

Interplay between the cell envelope and mobile genetic elements shapes gene flow in populations of the nosocomial pathogen Klebsiella pneumoniae

Mobile genetic elements (MGEs) drive genetic transfers between bacteria using mechanisms that require a physical interaction with the cellular envelope. In the high-priority multidrug-resistant nosocomial pathogens (ESKAPE), the first point of contact between the cell and virions or conjugative pili is the capsule. While the capsule can be a barrier to MGEs, it also evolves rapidly by horizontal gene transfer (HGT). Here, we aim at understanding this apparent contradiction by studying the covariation between the repertoire of capsule genes and MGEs in approximately 4,000 genomes of Klebsiella pneumoniae (Kpn). We show that capsules drive phage-mediated gene flow between closely related serotypes. Such serotype-specific phage predation also explains the frequent inactivation of capsule genes, observed in more than 3% of the genomes. Inactivation is strongly epistatic, recapitulating the capsule biosynthetic pathway. We show that conjugative plasmids are acquired at higher rates in natural isolates lacking a functional capsular locus and confirmed experimentally this result in capsule mutants. This suggests that capsule inactivation by phage pressure facilitates its subsequent reacquisition by conjugation. Accordingly, capsule reacquisition leaves long recombination tracts around the capsular locus. The loss and regain process rewires gene flow toward other lineages whenever it leads to serotype swaps. Such changes happen preferentially between chemically related serotypes, hinting that the fitness of serotype-swapped strains depends on the host genetic background. These results enlighten the bases of trade-offs between the evolution of virulence and multidrug resistance and caution that some alternatives to antibiotics by selecting for capsule inactivation may facilitate the acquisition of antibiotic resistance genes (ARGs).

A study of how the complex interaction between capsules and mobile genetic elements shapes gene flow in populations of Klebsiella pneumoniae reveals that capsule inactivation by phage pressure facilitates its subsequent re-acquisition by conjugation, and this loss and re-gain process influences the gene flow towards other lineages whenever it leads to serotype changes.     

Mobilization of vancomycin resistance by transposon-mediated fusion of a VanA plasmid with an Enterococcus faecium sex pheromone-response plasmid

A striking feature of recent outbreaks of vancomycin-resistant (Vm^R) enterococci is the apparent horizontal dissemination of resistance determinants. The plasmids pHKK702 and pHKK703 from Enterococcus faecium clinical isolate R7 have been implicated in the conjugal transfer of Vm^R. pHKK702 is a 41-kb plasmid that contains an element indistinguishable from the glycopeptide-resistance transposon Tn1546. pHKK703 is an approx. 55-kb putative sex pheromone-response plasmid that is required for conjugative mobilization of pHKK702. During experiments in which strain R7 was used as a donor, a highly conjugative Vm^R transconjugant was isolated that formed constitutive cellular aggregates. Restriction analyses and DNA hybridizations revealed that the transconjugant harbored a single plasmid of approx. 92 kb and this plasmid (pHKK701) was composed of DNA from both pHKK702 and pHKK703. Results from DNA sequence analyses showed that a 39-kb composite transposon (Tn5506) from pHKK702 had inserted into pHKK703. The left end of Tn5506 contained a single insertion sequence (IS) element, IS1216V2, whereas the right end was composed of a tandem IS structure consisting of the novel 1065-bp IS1252 nested within an IS1216V1 element. Transposition of Tn5506 from pHKK702 to pHKK703 created an 8-bp target sequence duplication at the site of insertion and interrupted an ORF (ORFX) that was 91% identical to that of prgX, a gene proposed to negatively regulate sex pheromone response of the E. faecalis plasmid, pCF10. We propose that the interruption of ORFX by Tn5506 led to the constitutive cellular aggregation phenotype and thereby enhanced the efficiency with which Vm^R was transferred. Similar IS1216V-mediated transposition events may contribute to the horizontal spread of glycopeptide resistance among enterococci in nature.     

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.     

Integrons,transposons and IS elements promote diversification of multidrug resistance plasmids and adaptation of their hosts to antibiotic pollutants from pharmaceutical companies

Plasmids are important vehicles for the dissemination of antibiotic resistance genes (ARGs) among bacteria by conjugation. Here, we determined the complete nucleotide sequences of nine different plasmids previously obtained by exogenous plasmid isolation from river and creek sediments and wastewater from a pharmaceutical company. We identified six IncP/P-1ε plasmids and single members of IncL, IncN and IncFII-like plasmids. Genetic structures of the accessory regions of the IncP/P-1ε plasmids obtained implied that multiple insertions and deletions had occurred, mediated by different transposons and Class 1 integrons with various ARGs. Our study provides compelling evidence that Class 1 integrons, Tn402-like transposons, Tn3-like transposons and/or IS26 played important roles in the acquisition of ARGs across all investigated plasmids. Our plasmid sequencing data provide new insights into how these mobile genetic elements could mediate the acquisition and spread of ARGs in environmental bacteria.     

Formation of the tandem repeat (IS30)2 and its role in IS30-mediated transpositional DNA rearrangements

Plasmids carrying two IS30 elements in the same orientation, as in the composite transposon Tn2706, are structurally unstable in Escherichia coli. A primary segregation product is formed by site-specific deletion of the sequences carried between the two IS30 elements. The resulting covalently closed replicon carries the two IS30 elements as tandem repeats separated by only 2 bp. This (IS30)₂ structure is extremely unstable, but it can nevertheless be isolated on its vector plasmid and, after purification, can be reintroduced into host cells by transformation. Among the descendants of transformants of recA ^? bacteria, replicated copies of the introduced (IS30)₂ structure are still present, together with various kinds of segregation products which provide evidence for the efficient generation of DNA rearrangements. Most abundant is the product of another site-specific recombination between two identical ends of the IS30 elements involved, which results in the presence of just one intact IS30 on the plasmid. Apart from this, and depending on the presence of appropriate targets for IS30 transposition, various transposition products of (IS30)₂ are also seen. Intramolecular reactions lead to DNA inversions and deletions with breakpoints other than IS30 ends. In intermolecular reactions inverse transposition occurs at high frequency and one also obtains simple transposition and cointegration. A mutational study revealed the requirement in cis of one intact IS30 transposase gene and of both proximal ends of the two IS30 elements concerned not only for the formation of (IS30)₂, but also for its further rearrangement reactions, including the efficient formation of site-specific deletions. A model is proposed, which postulates that (IS30)₂ intermediates play a key role in IS30 transposition pathways in which the formation of (IS30)₂ may be rate-limiting. Once this structure is formed, it gives rise to a burst of transpositional rearrangements in the subclone carrying (IS30)₂. Evolutionary implications of these findings are discussed.     

An insight into the exploration of proliferation of antibiotic resistance genes in high-fat diet induced obesity mice

Using mice as an animal model, we first demonstrated the significant proliferation of ARGs and the change of mobile genetic elements (MGEs) in high-fat diet induced obesity (DIO) mice, which the ermB and tnpA-03 genes mostly increased, illuminating that DIO could enrich the abundance of ARGs. Additionally, Lactobacillus sharply increased in the DIO mice and might contribute to the proliferation of ARGs and dramatical change of MGEs in the HFD groups. Finally, procrustes analysis showed the explanatory variables of the MGEs, the metabolites, and the microbial communities for the ARGs accounted for 94.3%, 53.4%, and 68.1%, respectively, and implying that MGEs might be the most direct factor affecting ARGs, and microbiota could be the main driver of the proliferation of ARGs in the DIO mice.     

Evidence that Tn5565, which includes the enterotoxin gene in Clostridium perfringens, can have a circular form which may be a transposition intermediate

The Clostridium perfringens enterotoxin gene is on a transposon-like element, Tn5565, integrated in the chromosome in human food poisoning strains. The flanking IS elements, IS1470 A and B, are related to IS30. The IS element found in the transposon, IS1469, is related to IS200 and has been found upstream of cpe in all Type A strains. PCR and sequencing studies from cell extracts and plasmid isolations of C. perfringens indicate that Tn5565 can form a circular form with the tandem repeat (IS1470)₂, similar to the transposition intermediates described for a number of IS elements.     

Enterobacteriaceae genome-wide analysis reveals roles for P1-like phage-plasmids in transmission of mcr-1, tetX4 and other antibiotic resistance genes

P1 -like phage-plasmids (PPs) are important gene vehicles in isolated pathogens. In this study, we conducted genome-wide and cross-species analysis of antimicrobial resistance genes (ARGs) from 35 ARG-positive P1-like PPs. LS-BSR analysis reveal that P1-like PPs had in common 7 highly variable regions and carried 48 different ARG subtypes. The most prevalent gene groups were the colistin resistance gene mcr-1 and a class 1 integron. Analysis of the flanking sequences of mcr-1 indicated an “IS30-mcr-1-ORF-IS30” as the core cluster. In particular, we found an mcr-1- and bla_CTX-M-55-coharboring large fusion P1-like PP. Also, tet(X4) was detected and flanking sequences indicated tet(X4)-bearing cluster can formed a larger size fusion plasmid mediated a wider spread via IS26 hotspots. Overall, this study demonstrated that P1-like PPs can not only mobilize a large number of ARGs in variable regions but also form larger hybrid P1-like PPs that would increase their ability to spread antimicrobial resistance.     

Efficient transposition of IS204-derived plasmids in Streptomyces coelicolor

In order to study functional gene expression in Streptomyces coelicolor, a mini-transposon encoding the apramycin resistance gene aac(3)IV within its inverted repeat (IR) boundaries was constructed based on IS204, which was previously identified in the genome of Nocardia asteroides YP21. The mini-transposon and IS204 transposase gene were then put on a kanamycin-resistant conjugative plasmid pDZY101 that can only replicate in Escherichia coli. After mating with S. coelicolor A3(2) M145, resistant colonies arose efficiently on both apramycin and kanamycin plates. Plasmid rescue indicated that entire plasmids were inserted into the M145 genome with cleavage at an inverted repeat junction formed by the right inverted repeat (IRR) and the last 18 bp of the transposase gene, while the left inverted repeat (IRL) was untouched. Southern blot analysis of the mutants using an aac(3)IV gene probe showed that transposition of plasmid pDZY101 was genetically stable, with a single-copy insertion within the S. coelicolor M145 genome. Several mutagenesis libraries of S. coelicolor M145 were constructed using plasmid pDZY101 derivatives and the transposon insertion site was determined. The correlation between novel mutant phenotypes and previously uncharacterized genes was established and these transposon locations were widely scattered around the genome.     

Evidence that the insertion events of IS2 transposition are biased towards abrupt compositional shifts in target DNA and modulated by a diverse set of culture parameters

Insertion specificity of mobile genetic elements is a rather complex aspect of DNA transposition, which, despite much progress towards its elucidation, still remains incompletely understood. We report here the results of a meta-analysis of IS2 target sites from genomic, phage, and plasmid DNA and find that newly acquired IS2 elements are consistently inserted around abrupt DNA compositional shifts, particularly in the form of switch sites of GC skew. The results presented in this study not only corroborate our previous observations that both the insertion sequence (IS) minicircle junction and target region adopt intrinsically bent conformations in IS2, but most interestingly, extend this requirement to other families of IS elements. Using this information, we were able to pinpoint regions with high propensity for transposition and to predict and detect, de novo, a novel IS2 insertion event in the 3′ region of the gfp gene of a reporter plasmid. We also found that during amplification of this plasmid, process parameters such as scale, culture growth phase, and medium composition exacerbate IS2 transposition, leading to contamination levels with potentially detrimental clinical effects. Overall, our findings provide new insights into the role of target DNA structure in the mechanism of transposition of IS elements and extend our understanding of how culture conditions are a relevant factor in the induction of genetic instability.     

Two Amino Acid Residues of Transposase Contributing to Differential Transposability of IS1 Elements in Escherichia coli

Escherichia coli W3110 contains four types of IS1 elements in the chromosome. Using an insertion element entrapping system, we collected 116 IS1 plasmid insertion mutants, which resulted from a minimum of 26 independent IS1 insertion events. All of them had insertions of IS1 of the IS1A (IS1E and IS1G) type. Inspection of the transposase sequences of the four IS1 types and the IS1 of the resistance plasmid R100 showed that two amino acid residues, His-193 and Leu-217 of transposase, might contribute to differential transposability of IS1 elements in W3110. The two amino acid residues of the transposase in IS1A (IS1E and IS1G) were altered separately by site-directed mutagenesis, and each mutant was found to mediate transposition at a frequency about 30-fold lower than that of IS1A (IS1E and IS1G). Thus, the assumption that His-193 and Leu-217 of transposase contribute to differential transposability of IS1 elements in W3110 was confirmed.     

Conjugative Interaction Induces Transposition of ISPst9 in Pseudomonas stutzeri AN10

ISPst9 is an ISL3-like insertion sequence (IS) that was recently described in the naphthalene-degrading organism Pseudomonas stutzeri strain AN10. In this paper we describe a novel strong IS regulation stimulus; transposition of ISPst9 is induced in all P. stutzeri AN10 cells after conjugative interaction with Escherichia coli. Thus, we observed that in all P. stutzeri AN10 cells that received genetic material by conjugation the ISPst9 genomic dose and/or distribution was changed. Furthermore, ISPst9 transposition was also observed when P. stutzeri AN10 cells were put in contact with the plasmidless conjugative strain E. coli S17-1λ_pir, but not when they were put in contact with E. coli DH5α (a nonconjugative strain). The mechanism of ISPst9 transposition was analyzed, and transposition was shown to proceed by excision from the donor DNA using a conservative mechanism, which generated 3- to 10-bp deletions of the flanking DNA. Our results indicate that ISPst9 transposes, forming double-stranded DNA circular intermediates consisting of the IS and a 5-bp intervening DNA sequence probably derived from the ISPst9 flanking regions. The kinetics of IS circle formation are also described.     

Defence systems and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) is a fundamental process in prokaryotic evolution, contributing significantly to diversification and adaptation. HGT is typically facilitated by mobile genetic elements (MGEs), such as conjugative plasmids and phages, which often impose fitness costs on their hosts. However, a considerable number of bacterial genes are involved in defence mechanisms that limit the propagation of MGEs, suggesting they may actively restrict HGT. In our study, we investigated whether defence systems limit HGT by examining the relationship between the HGT rate and the presence of 73 defence systems across 12 bacterial species. We discovered that only six defence systems, three of which were different CRISPR-Cas subtypes, were associated with a reduced gene gain rate at the species evolution scale. Hosts of these defence systems tend to have a smaller pangenome size and fewer phage-related genes compared to genomes without these systems. This suggests that these defence mechanisms inhibit HGT by limiting prophage integration. We hypothesize that the restriction of HGT by defence systems is species-specific and depends on various ecological and genetic factors, including the burden of MGEs and the fitness effect of HGT in bacterial populations.     

Air pollution could drive global dissemination of antibiotic resistance genes

Antibiotic-resistant pathogens pose a significant threat to human health. Several dispersal mechanisms have been described, but transport of both microbes and antibiotic resistance genes (ARGs) via atmospheric particles has received little attention as a pathway for global dissemination. These atmospheric particles can return to the Earth’s surface via rain or snowfall, and thus promote long-distance spread of ARGs. However, the diversity and abundance of ARGs in fresh snow has not been studied and their potential correlation with particulate air pollution is not well explored. Here, we characterized ARGs in 44 samples of fresh snow from major cities in China, three in North America, and one in Europe, spanning a gradient from pristine to heavily anthropogenically influenced ecosystems. High-throughput qPCR analysis of ARGs and mobile genetic elements (MGEs) provided strong indications that dissemination of ARGs in fresh snow could be exacerbated by air pollution, severely increasing the health risks of both air pollution and ARGs. We showed that snowfall did effectively spread ARGs from point sources over the Earth surface. Together our findings urge for better pollution control to reduce the risk of global dissemination of antibiotic resistance genes.Subject terms: Environmental sciences, Air microbiology     

Complete Nucleotide Sequence of pSK41: Evolution of Staphylococcal Conjugative Multiresistance Plasmids

The 46.4-kb nucleotide sequence of pSK41, a prototypical multiresistance plasmid from Staphylococcus aureus, has been determined, representing the first completely sequenced conjugative plasmid from a gram-positive organism. Analysis of the sequence has enabled the identification of the probable replication, maintenance, and transfer functions of the plasmid and has provided insights into the evolution of a clinically significant group of plasmids. The basis of deletions commonly associated with pSK41 family plasmids has been investigated, as has the observed insertion site specificity of Tn552-like β-lactamase transposons within them. Several of the resistance determinants carried by pSK41-like plasmids were found to be located on up to four smaller cointegrated plasmids. pSK41 and related plasmids appear to represent a consolidation of antimicrobial resistance functions, collected by a preexisting conjugative plasmid via transposon insertion and IS257-mediated cointegrative capture of other plasmids.     

Mobility and Generation of Mosaic Non-Autonomous Transposons by Tn3-Derived Inverted-Repeat Miniature Elements (TIMEs)

Functional transposable elements (TEs) of several Pseudomonas spp. strains isolated from black shale ore of Lubin mine and from post-flotation tailings of Zelazny Most in Poland, were identified using a positive selection trap plasmid strategy. This approach led to the capture and characterization of (i) 13 insertion sequences from 5 IS families (IS3, IS5, ISL3, IS30 and IS1380), (ii) isoforms of two Tn3-family transposons – Tn5563a and Tn4662a (the latter contains a toxin-antitoxin system), as well as (iii) non-autonomous TEs of diverse structure, ranging in size from 262 to 3892 bp. The non-autonomous elements transposed into AT-rich DNA regions and generated 5- or 6-bp sequence duplications at the target site of transposition. Although these TEs lack a transposase gene, they contain homologous 38-bp-long terminal inverted repeat sequences (IRs), highly conserved in Tn5563a and many other Tn3-family transposons. The simplest elements of this type, designated TIMEs (Tn3 family-derived Inverted-repeat Miniature Elements) (262 bp), were identified within two natural plasmids (pZM1P1 and pLM8P2) of Pseudomonas spp. It was demonstrated that TIMEs are able to mobilize segments of plasmid DNA for transposition, which results in the generation of more complex non-autonomous elements, resembling IS-driven composite transposons in structure. Such transposon-like elements may contain different functional genetic modules in their core regions, including plasmid replication systems. Another non-autonomous element “captured” with a trap plasmid was a TIME derivative containing a predicted resolvase gene and a res site typical for many Tn3-family transposons. The identification of a portable site-specific recombination system is another intriguing example confirming the important role of non-autonomous TEs of the TIME family in shuffling genetic information in bacterial genomes. Transposition of such mosaic elements may have a significant impact on diversity and evolution, not only of transposons and plasmids, but also of other types of mobile genetic elements.