Timothy M. Ghaly
Ian T. Paulsen
Ammara Sajjad1
Sasha G. Tetu
Michael R. Gillings- 1Department of Biological Sciences, Macquarie University, Sydney, NSW, Australia
- 2Department of Molecular Sciences, Macquarie University, Sydney, NSW, Australia
- 3ARC Centre of Excellence in Synthetic Biology, Macquarie University, Sydney, NSW, Australia
Acinetobacter species are emerging as major nosocomial pathogens, aided by their ability to acquire resistance to all classes of antibiotics. A key factor leading to their multi-drug resistance phenotypes is the acquisition of a wide variety of mobile genetic elements, particularly large conjugative plasmids. Here, we characterize a family of 21 multi-drug resistance mega-plasmids in 11 different Acinetobacter species isolated from various locations across the globe. The plasmid family exhibits a highly dynamic and diverse accessory genome, including 221 antibiotic resistance genes (ARGs) that confer resistance to 13 classes of antibiotics. We show that plasmids isolated within the same geographic region are often evolutionarily divergent members of this family based on their core-genome, yet they exhibit a more similar accessory genome. Individual plasmids, therefore, can disseminate to different locations around the globe, where they then appear to acquire diverse sets of accessory genes from their local surroundings. Further, we show that plasmids from several geographic regions were enriched with location-specific functional traits. Together, our findings show that these mega-plasmids can transmit across species boundaries, have the capacity for global dissemination, can accumulate a diverse suite of location-specific accessory genes, and can confer multi-drug resistance phenotypes of significant concern for human health. We therefore highlight this previously undescribed plasmid family as a serious threat to healthcare systems worldwide. These findings also add to the growing concern that mega-plasmids are key disseminators of antibiotic resistance and require global surveillance.
Introduction
Acinetobacter is a diverse genus of bacteria that commonly occur in soil and water, and in association with animals and plants (Fondi et al., 2010). Several species are emerging as major nosocomial pathogens (Almasaudi, 2018). In particular, Acinetobacter baumannii, Acinetobacter nosocomialis, Acinetobacter pitti, Acinetobacter ursingii, Acinetobacter haemolyticus, and Acinetobacter calcoaceticus are now a serious threat to human health (Antunes et al., 2014; Chusri et al., 2014; Weber et al., 2016; Castro-Jaimes et al., 2020; Rivera-Izquierdo et al., 2020). Clinical isolates are often characterized by their multi-drug resistance, with some isolates being resistant to all classes of antibiotics (Göttig et al., 2014). One of the key features leading to their multi-drug resistance is the acquisition of a wide variety of mobile DNA elements (Fournier et al., 2006; Hamidian et al., 2014, 2016; Lean and Yeo, 2017). In particular, Acinetobacter strains are known to harbor a strikingly diverse pool of plasmids (Fournier et al., 2006; Fondi et al., 2010; Merino et al., 2010; Lean and Yeo, 2017; Salto et al., 2018).
The evolutionary and ecological dynamics of Acinetobacter plasmids are particularly interesting as their host species can colonize a wide variety of environments and display diverse metabolic activities, even among closely related strains (Hawkey and Munday, 2004; Bach and Gutnick, 2006; Albarracín et al., 2012; Garcia-Garcera et al., 2017). Such variety in niche adaptation is largely driven by the open pangenome of their plasmids (Fondi et al., 2010). Plasmids, in general, can rapidly diversify through a series of recombinatorial and transpositional events that generate novel mosaic elements (Stokes and Gillings, 2011). These mosaic elements are characterized by a nested structure of diverse mobile elements, including transposons and “clinical” forms of integrons (Ghaly et al., 2017; Gillings, 2017), each of which are prolific at acquiring diverse antibiotic resistance genes (ARGs; Babakhani and Oloomi, 2018; Ghaly et al., 2020). Humans, largely through the use of various selective agents, such as heavy metals, disinfectants, and antibiotics, are increasing the frequency of such occurrences (Gillings and Stokes, 2012; Gillings et al., 2018). The same selective agents also upregulate the horizontal transmission mechanisms of these mosaic elements, leading to the wide-spread dissemination of their genetic cargo (Beaber et al., 2004; Schreiber et al., 2013; Ghaly and Gillings, 2018). In clinical settings, this has led to the emergence of multi-drug resistant pathogens that are well adapted to a nosocomial lifestyle.
An extensive range of resistance genes in Acinetobacter are plasmid-borne (Cerezales et al., 2020). In particular, large plasmids of more than 100 kb, referred to as mega-plasmids, are known to carry arrays of multiple resistance genes (Hamidian et al., 2016; Nigro and Hall, 2017). Mega-plasmids, in general, are characterized by their mosaic structure and carry genetic regions that originate from various ancestral sources (Pesesky et al., 2019). Recently, a family of mega-plasmids has been identified as a major disseminator of multi-drug resistance among Pseudomonas species (Cazares et al., 2020). These mega-plasmids carried large arrays of complex and dynamic resistance regions, suggesting that they play a key role in the spread of resistance.
Global surveillance of such mega-plasmids using a whole-genome sequencing approach has been highlighted as a top priority for public health management (Cazares et al., 2020). However, fully sequencing these mega-plasmids with short-read data remains difficult as they contain long repeat regions (Arredondo-Alonso et al., 2017). Long-read sequencing, however, allows these mega-plasmids to be fully resolved (Botelho et al., 2019). The accuracy of assemblies can be improved by combining short-read and long-read sequencing approaches.
Here, we use a combination of PacBio long-read and Illumina short-read sequencing to fully assemble four multi-resistance mega-plasmids resident in both clinical and environmental Acinetobacter strains. We collected an additional 17 publicly available complete plasmid sequences within this family, isolated from various locations across the world. This previously undescribed mega-plasmid family, found in 11 different Acinetobacter species so far, shares a tight core-genome, yet varies extensively in its accessory genome. By characterizing its pangenome, we show that this plasmid family can rapidly disperse around the globe, where it subsequently accumulates a diverse suite of niche-adaptive accessory genes.
Materials and Methods
Sampling, DNA Extraction, and Sequencing
The novel plasmids described in this study were recovered from Acinetobacter species. Three were isolated from Acinetobacter strains inhabiting the digestive tracts of wild-caught prawns, harvested from the Australian east coast (pR4WN-1BD1, pR4WN-12CE1, and pR4WN-E10B), and one from an ICU patient at Westmead Hospital, Sydney, Australia (pWM98B). The method by which each Acinetobacter host was isolated from the wild-caught prawns (Gillings et al., 2009; Sajjad et al., 2011) and the Westmead Hospital patient (Valenzuela et al., 2007) have been previously described. These bacterial isolates were of interest because they showed sequence similarity in the genomic landscape flanking a previously described class 1 integron (Gillings et al., 2009). To investigate this further, DNA was extracted from pure cultures using a standard phenol:chloroform protocol (Sambrook and Russell, 2006). Whole genome de novo sequencing was performed using both long-read and 100 bp paired-end sequencing on the PacBio® RS and the Illumina® HiSeq 2500 systems, respectively. Sequencing was carried out at the Macrogen Sequencing Facility (Seoul, South Korea).
Plasmid Assembly and Polishing
PacBio long-reads were assembled with Flye v2.7.1 (Kolmogorov et al., 2019) with three iterations of Flye’s in-built polisher (parameters: –plasmids -i 3). The Flye assemblies were then further polished with the HiSeq paired-end reads using Pilon v1.23 (Walker et al., 2014). Polishing with Pilon first involved mapping all paired-end reads to the PacBio assemblies using BWA v0.7.12-r1039 (Li and Durbin, 2009). The BWA mapping was then used by Pilon to polish the assemblies. Several rounds of polishing (3–4 iterations) were implemented until no further changes were made to the previous round. Polished contigs were then identified to be circular using Circlator v1.5.5 (Hunt et al., 2015), which outputs a linearized version of each circular contig (parameters: all -- assembler canu). Circularized contigs were then re-polished with the paired-end reads using Pilon. All assemblies required a further two rounds of Pilon polishing until no changes were made to the previous round.
Taxonomic Labeling and Collecting Related Plasmids
Contigs representing the complete chromosome from each assembly were assigned a taxonomy by best BLASTn match using BLAST v2.7.1 (Madden, 2013) against all complete bacterial genomes in RefSeq (downloaded: June 25, 2020). Taxonomic identification was then confirmed with Kraken v1.0 (Wood and Salzberg, 2014; parameters: -- threads 24).
From each of the four samples, a complete sequence was also assembled for the plasmid carrying the integron of interest. Sequence coverage statistics for each of these were generating by mapping sequence reads to each plasmid assembly using Minimap2 v2.17-r941 (Li, 2018) with the parameters (-ax map-pb) and (-ax sr) for PacBio and Illumina mapping, respectively. Sequence coverage details were then extracted from each read mapping using BEDTools v2.27.1 (Quinlan and Hall, 2010; parameters: genomeCoveragBed).
To identify all related plasmids that had previously been sequenced and submitted to NCBI, each plasmid was used to query the complete nr/nt database as of June 26, 20201 using BLASTn. All plasmids that were considered related had a BLAST query cover greater than 64% and a percent identity greater than 98.5%. This resulted in an additional 17 plasmids for inclusion in the detailed analyses below (GenBank accessions: CP010351, CP029396, CP032285, CP033130, CP033531, CP033569, CP038010, CP042365, CP042557, CP043053, CP043309, CP048828, KX426227, MH220285 – MH220287, and MK134375). Plasmid alignment visualizations were generated using Chromatiblock v0.4.2 (Sullivan and van Bakel, 2019) and AliTV v1.0.6 (Ankenbrand et al., 2017).
To determine to which Rep family the plasmids belonged, we identified replication initiation proteins by annotation with InterProScan v5.44-79.0 (Jones et al., 2014; parameters: -goterms) and by BLASTP alignments with representative Rep proteins from every known Acinetobacter Rep family (available as Supplementary Figure S2 from Salgado-Camargo et al., 2020).
Plasmid Pan‐ and Core-Genome Analyses
Gene and protein sequences for all 21 plasmids were predicted using Prokka v1.13 (parameters: -- kingdom Bacteria). Plasmid pan‐ and core-genome analyses were carried out using Perl scripts available in the GET_HOMOLOGUES v20092018 software package (Contreras-Moreira and Vinuesa, 2013; Vinuesa and Contreras-Moreira, 2015). For the plasmid pan-genome, this involved generating two pangenome matrices based on the COGtriangles (Kristensen et al., 2010) and OMCL (Li et al., 2003) clustering algorithms using the get_homologues.pl script. The parameters (-X -c -n 8 -t 0 -A -G) and (-X -c -n 8 -t 0 -A -G -M) were implemented, respectively. Based on the intersection of the two sets of clustered sequences, a consensus pangenome matrix was obtained using the compare_clusters.pl script (parameters: -m -T). The identified plasmid core-genome was based on the intersection of the BDBH, COGtriangles, and OMCL clustering algorithms using the get_homologues.pl script with parameters (-t 21 -e), (-G t0), and (-M -t 0), respectively.
The plasmid pangenome was divided into four occupancy classes. The core-genome was made up of the “soft core” (Kaas et al., 2012), which consisted of the genes found in >95% of plasmids, and the “core,” which contained genes found in all plasmids. The accessory genome was made up of the “cloud” (Wolf et al., 2012), which contained the rare genes found in only one to two plasmids, and the “shell” (Wolf et al., 2012), which comprised the remaining, moderately conserved genes.
The primary Perl script, get_homologues.pl with the above parameters, also estimated the minimum core-genome size for this family of plasmids. The estimate was based on 10 random samples of the 21 plasmids with a fitted curve following the function proposed by Tettelin et al. (2005).
To generate a maximum-likelihood tree that clustered plasmids based on the presence and absence of all accessory genes in the pangenome, the pangenome matrix, reduced to a binary format, was used. The tree was generated using IQ-TREE v1.6.12 (Nguyen et al., 2015; Hoang et al., 2018) with 1,000 bootstrap replicates (parameters: -alrt 1000 -bb 1000). The resulting tree was annotated with FigTree v.1.4.3 (Rambaut, 2017) with midpoint root.
Phylogenetic Analysis
The evolutionary relationship between plasmids was inferred from the plasmid core-genome. Appropriate phylogenetic loci were identified using the GET_PHYLOMARKERS software package (Vinuesa et al., 2018) with default parameters. The GET_PHYLOMARKERS uses the single-copy plasmid core-genome clusters identified by GET_HOMOLOGUES and removes recombinant sequences as well as loci that produce outlier phylogenies. This resulted in eight high-quality phylogenetic markers. These included genes that encode a TraX conjugative transfer protein, an integrase, three membrane-associated proteins, and three hypothetical proteins of unknown functions. The nucleotide sequences were concatenated and aligned using MAFFT v7.271 (Katoh and Standley, 2013; paramters: -localpair -maxiterate 1000). The best-fit substitution model to suit the alignment was determined using ModelFinder (Kalyaanamoorthy et al., 2017), and a maximum-likelihood tree was generated using IQ-TREE v1.6.12 (Nguyen et al., 2015; Hoang et al., 2018) with 1,000 bootstrap replicates (parameters: -m MFP -alrt 1000 -bb 1000). The resulting tree was annotated with FigTree v.1.4.3 (Rambaut, 2017) with midpoint root.
To test whether plasmids have transmitted horizontally across species boundaries, a co-phylogeny was generated to compare host cell and plasmid phylogenies. We used a whole-genome-based phylogeny of the host chromosomes using the GET_PHYLOMARKERS software package (Vinuesa et al., 2018) as described above. Five of the plasmids retrieved from GenBank did not have a corresponding host chromosome sequence publicly available and were thus excluded from the co-phylogenetic analysis. A co-phylogeny comparing the plasmids and their respective host strains was produced using the R package phytools v0.7-47 (Revell, 2012), which rotates the nodes of both trees to optimize congruence between the tips.
Tree topologies were quantitatively compared using a normalized PH85 tree topological distance (nPH85; Geoghegan et al., 2017). This was used to compare the two plasmid trees based on their core and accessory genomes, respectively, as well as to compare the plasmid and host co-phylogenies. The nPH85 distance quantifies the number of bipartitions that are not shared between the two trees. This is done by randomizing the tip labels for one of the trees, and calculating the PH85 for each replicate to generate a null distribution of the tree topology distance. The nPH85 metric ranges from 0, for identical trees, to 1 for trees that have no clades in common. We calculated nPH85 using the R package NELSI v0.2 (Ho et al., 2015) with 1,000 randomizations, which is shown to be robust even for very large trees (Geoghegan et al., 2017).
Gene Ontology Enrichment Analysis
A gene ontology (GO) enrichment analysis was done to investigate what functional traits were prevalent among the plasmids from each geographic region. To do this, all proteins encoded by the plasmid pangenome were annotated based on the terms from the GO database (Ashburner et al., 2000; The Gene Ontology Consortium, 2019) using InterProScan v5.44-79.0 (Jones et al., 2014; parameters: -goterms). Proteins from plasmids in each geographic region were tested for GO terms that were significantly enriched relative to the complete plasmid pangenome using WEGO v2.0 (Ye et al., 2018).
Screening for Antibiotic, Biocide and Metal Resistance Genes, and IS Elements
Antibiotic resistance genes were detected using ABRicate v0.8 (Seemann, 2018) using default parameters. To determine which ARGs were part of integron gene cassettes, IntegronFinder v1.5.1 (Cury et al., 2016) was used (parameters: -local_max -func_annot). We used CARD (Alcock et al., 2020) to confirm the class of antibiotic to which each gene conferred resistance. To screen for metal and biocide resistance genes, the complete manually-curated BactMet2 database (Pal et al., 2014) was used. All proteins were aligned against the BacMet2 database using DIAMOND v0.8.33.95 (Buchfink et al., 2015) with a cut-off alignment criteria of e-value 1 × 10−5, 80% similarity, and 90% subject cover length (parameters: --evalue 0.00005 --id 80 --subject-cover 90 --max-target-seqs 1).
Sequence Data Availability
The complete sequences of the four novel plasmids described here are available in GenBank as accessions MT742180 – MT742183. The complete chromosome sequences of their respective hosts are also available in GenBank as accessions CP059078 – CP059081.
Results and Discussion
Plasmid Family and Host Species
Using a combination of PacBio and Illumina sequencing, we generated complete assemblies of four novel mega-plasmids, all of which were recovered from Acinetobacter species. Each plasmid had a PacBio sequence coverage > 300x and an Illumina sequence coverage > 500x, except for pWM98B (46x Illumina coverage). Three of these (pR4WN-1BD1, pR4WN-12 CE1, and pR4WN-E10B) were isolated from the digestive tracts of wild-caught prawns and one (pWM98B) from an ICU patient. Another 17 related Acinetobacter plasmids were identified in the GenBank database, isolated from various locations across the world and were included in all subsequent analyses. The 21 mega-plasmids were ranged from 250 to 400 kb in size and were harbored by 11 different Acinetobacter species (see Supplementary Table S1 for more detailed plasmid information and taxonomy of host bacteria).
We have identified four of the 21 plasmids (pXBB1-9, pAHTJR1, p34AB, and pOXA23_010062) to be part of the Rep-3 superfamily group, based on their replication initiation (Rep) protein. Intriguingly, no Rep protein could be identified in any of the other plasmids, suggesting that it is not needed for successful replication. Similarly, other Acinetobacter plasmids, notably pRAY-type plasmids, have no known Rep protein (Hamidian et al., 2012; Lean and Yeo, 2017). It has been proposed that these rep-lacking Acinetobacter plasmids may use their host machinery for DNA melting and primase activity during plasmid replication, as is the case in ColE1-type plasmids (Lean and Yeo, 2017). It is possible that these mega-plasmids have acquired this replication strategy and subsequently lost their rep gene.
Plasmid Pan‐ and Core-Genome Analyses
The plasmid pangenome, comprising 1,303 homologous gene clusters, was divided into four occupancy classes (see the section “Methods” for definitions), where the “core” and “soft core” classes make up the core-genome and the “shell” and “cloud” classes make up the accessory genome (Figure 1A). The vast majority of genes carried by this plasmid family are accessory genes (Figures 1A,C), with a minimum set of 40 genes comprising the plasmid core-genome (Figures 1B,C). Individual plasmids carried 95–340 accessory genes each (Supplementary Table S1). It is thus evident that these mega-plasmids have the capacity to accumulate large and diverse sets of accessory genes that drives their diversification.
Figure 1. Pan‐ and core-genomic analyses of the mega-plasmid family. (A) The plasmid pangenome can be divided into occupancy classes based on the number of plasmids contributing to each homologous gene cluster. The majority of genes belong to the accessory genome, constituting the cloud and shell classes. (B) The estimate of the coregenome size (genes) for this plasmid family was generated from 10 random samples of the 21 plasmids. The fitted curve follows the function proposed by Tettelin et al. (2005). (C) A pangenome matrix showing the presence (red) and absence (cream) of genes among the plasmids. The Acinetobacter species in which each plasmid was resident in are annotated on the right. Plasmids share a tight core-genome, yet, vary extensively in their accessory genomes.
The plasmids consist of a number of core conserved genomic blocks separated by highly variable regions (Figure 2A). In particular, there are two “hotspots” of considerable genomic variability, labeled as “Hotspot 1” and “Hotspot 2” (Figures 2A,B). Hotspots 1 and 2 ranged from 7,000–140,000 bp and 15,000–90,000 bp in length, respectively. Alignments of the two hotspots showed a significant number of insertions, deletions, and large-scale rearrangements (Figure 2B; Supplementary Figures S1, S2). We speculate that these highly dynamic regions represent locations were genomic complexity can arise without interfering with core genes or regulatory networks. Interestingly, ARGs were overrepresented in Hotspot 2. This region, which on average constitutes less than 20% of each plasmid, contains greater than 65% of all ARGs (Supplementary Figures S1, S2).
Figure 2. Alignment of the mega-plasmid family. (A) A global alignment view of the complete plasmid sequences in this family. Each row represents a single plasmid. Core genomic blocks are visualized as vertically aligned solid rectangles that are colored according to their position in the genome. Non-core blocks are visualized as patterned rectangles. Sequences unique to a single plasmid are depicted as solid black lines. Note that, there are two “hotspots” of considerable genomic variability, labeled as “Hotspot 1” and “Hotspot2.” (B) A local alignment view of each hotspot region from five example plasmids. Homologous regions are shown by links that are colored according to their percent nucleotide identity. Locations of antibiotic resistance genes (ARGs) and transposases are annotated on each plasmid. Interestingly, Hotspot 2, on average, constitutes less than 20% of each plasmid, yet contains greater than 65% of all ARGs among these plasmids. Note that, these hotspots have undergone a significant amount of insertions, deletions, and rearrangements. To view the alignments of all 21 plasmids for hotspots 1 and 2, see Supplementary Figures S1, S2, respectively.
Plasmid Phylogeny
The evolutionary relationship between plasmids was inferred from the plasmid core-genome using a maximum-likelihood method (Figure 3). In addition, plasmids were clustered based on the presence and absence of all 1,262 accessory genes, also using a maximum-likelihood method (Figure 3). Here, we show that plasmid phylogeny based on core genes does not cluster well according to geographic region. However, based on the presence and absence of all accessory genes, the plasmids do cluster according to their isolation location. Quantitative comparison between the two tree topologies was measured using an nPH85 metric (Geoghegan et al., 2017). The nPH85 metric ranges from 0, for identical trees, to 1, for trees that have no clade bipartitions in common. The accessory and core-genome tree comparison returned an nPH85 value of 0.84, indicating a high-level of incongruence.
Figure 3. Plasmid core-genome phylogram and accessory genome dendrogram. (left) A maximum-likelihood tree showing the inferred plasmid phylogeny based on the plasmid core-genome and (right) a maximum-likelihood tree based on the presence and absence of all 1,262 accessory genes in the plasmid pangenome. Branch labels signify SH-aLRT support (%)/bootstrap support (%). Plasmids are named in the format of “plasmid-ID”_”Isolation-source” and are colored based on the geographic region from which they were sampled. Together, the two trees show that plasmids from the same geographic region are not necessarily evolutionarily closely related, yet they exhibit a more similar accessory genome.
Together, these data show that plasmids from the same geographic region are not necessarily evolutionarily closely related, yet they exhibit a more similar accessory genome. This strongly indicates that (a) the plasmids can rapidly disseminate across the globe, given that closely related sister plasmids can be isolated from different parts of the world, and (b) these plasmids preferentially acquire their repertoire of accessory genes from their local environment.
In contrast, the pBT2436-like mega-plasmid family that is driving the dissemination of multi-drug resistance in Pseudomonas (Cazares et al., 2020), exhibit similar topologies in their core and accessory genome trees, both of which cluster according to geographic location. This suggests that the Acinetobacter mega-plasmids described here have a comparatively greater capacity for global transmission and local acquisition of accessory genes.
Further, we show that the topologies of plasmid and host phylogenetic trees are highly incongruent (nPH85 = 1.0; Figure 4). The host and plasmid phylogenies thus share no clade bipartitions in common. Such co-phylogenetic analyses can be used to distinguish between a history of host-plasmid co-divergence (co-phylogenetic congruence) and cross-species transmission (co-phylogenetic incongruence). Our findings provide evidence that cross-species horizontal transmission among Acinetobacter species is highly prevalent for this plasmid family.
Figure 4. A co-phylogeny linking the host strain (left) and plasmid (right) phylogenies. Connections link host strains with their resident plasmid. Nodes have been rotated to maximize congruency between the two trees. Host phylogeny was inferred from the host core genomes using a maximum-likelihood method. The inferred plasmid phylogeny was based on the plasmid core-genome also using a maximum-likelihood approach. The high level of incongruence between the phylogenies (nPH85 = 1.0) indicates that cross-species transmission is highly prevalent among this family of plasmids.
Location-Specific Enrichment of Functional Traits
Since the accessory genomes of the plasmids cluster according to geographic region (Figure 3), we investigated whether there were any functional traits that were prevalent among the plasmids from each region. For this, we annotated proteins based on terms from the GO database. We found that plasmids from four of the seven regions had GO terms that were significantly enriched relative to the plasmid pangenome (Chi-squared tests, p < 0.05; Figure 5).
Figure 5. Geographic regions that had plasmids enriched with functional traits relative to the plasmid pangenome. Functions are based on the gene ontology (GO) database and include only those that were significantly enriched compared to the plasmid pangenome (Chi-squared tests, p < 0.05).
From these results, we can infer the broad-scale selection pressures within each region. For example, plasmids from Western China were largely enriched with genes involved in detoxification, metabolism, and cyclic compound binding (Figure 5). We speculate that these are traits that might be selected in environments of greater pollution. This is particularly significant as the Western China plasmids are polyphyletic (Figure 3), but they share common accessory traits, presumably a result of co-evolution by acquiring accessory genes from the local environment.
Genetic Cargo of Clinical Relevance
Collectively, the 21 plasmids carried 221 ARGs. All plasmids carried between 3 and 24 ARGs each, except for the plasmid isolated from Siberia (pALWED1.1), which only had one (tetracycline resistance). Together, the 221 genes encoded 35 different proteins known to confer resistance to 13 classes of antibiotics (Figure 6, full list of all ARGs in Supplementary Table S2). Of these, 31% (n = 67) were part of gene cassettes associated with class 1 integrons. In total, there were 20 class 1 integrons found in 14 of the plasmids. Consequently, class 1 integrons play a significant role in the accumulation of antibiotic resistance determinants in this plasmid family.
Figure 6. Diversity and abundance of ARGs. The number of resistance genes (color scale bar) is shown for each plasmid (rows). Plasmid names are colored according to the geographic region from which they were sampled, as described in Figure 3. The Acinetobacter species in which each plasmid was resident in are annotated on the right. Resistance genes are grouped according to the class of antibiotic that they are known to confer resistance to.
In particular, several of the class 1 integrons were associated with the same flanking miniature inverted-repeat transposable elements (MITEs; Figure 7). MITEs are non-autonomous mobile elements that transpose via molecular machinery provided in trans (Delihas, 2008). Fifteen of the 21 plasmids carried an identical pair of MITEs, each sharing the same insertion site (Figure 7A). Of these, 13 of the MITE pairs were flanking at least one class 1 integron carrying multiple ARG cassettes. Collectively, the MITEs were associated with 15 different ARGs, two virulence factors, and three biocide/metal resistance genes. There were eight unique MITE variants that we grouped into five different MITE “types” based on their cargo of resistance genes (Figure 7B). Interestingly, the same MITE type often occurred in multiple plasmid lineages and could also be present on plasmids from different geographic regions (Figure 7B). This suggests that the MITE-resistance region is highly dynamic and transmissible. Indeed, the same MITE-integron complex has also been detected within a clinical A. baumannii strain isolated from a human urine sample in Portugal (Domingues et al., 2011). The Portuguese MITE structure, however, was inserted into a different genomic context, further suggesting its high potential for transmission among Acinetobacter strains.
Figure 7. Genetic map of the miniature inverted-repeated transposable element (MITE)-integron complex present among 15 of the 21 plasmids. The insertion site of the MITEs is shown in (A). The structure of each MITE variant is shown in (B) and is grouped according to MITE “type.” The plasmids that carry each MITE variant are listed above/below each genetic map. Plasmid names are colored according to the geographic region from which they were sampled, as described in Figure 3. Directional blocks represent genes, while orange arrows represent MITEs as annotated in the top right corner of the figure. All maps are drawn to scale. To the right of each MITE type is the inferred plasmid core-genome phylogeny, highlighting each plasmid that carries the corresponding MITE type. Note that, the same MITE type often occurred in multiple plasmid lineages and could also be present on plasmids from different geographic regions. This suggests that this MITE-resistance region is highly dynamic and transmissible.
In addition, the plasmids carried 108 genes conferring resistance to bacterial biocides and heavy metals. Together, these genes confer resistance to six heavy metals (Hg, Cu, Au, Cr, Ni, and Co) and nine classes of biocides (acridines, azines, xanthenes, organosulfates, organomercury compounds, quaternary ammonium compounds, biguanides, diamindines, and paraquats). All of the detected metal/biocide resistance genes are part of efflux systems. For a full list of all biocide and metal resistance genes, see Supplementary Table S3.
Conclusion
The Acinetobacter mega-plasmids characterized here appear to be highly proficient in acquiring niche-adaptive accessory genes. Thus, the significant collection of ARGs observed within many of these is likely a selective response to exposure to one or more antibiotics. This is of serious concern, as global antibiotic usage continues to increase, particularly for antibiotics of “last-resort” (Klein et al., 2018). As an often-overlooked consequence of this, antibiotics are polluting the environment, radiating from areas of human populations and agricultural areas (Kümmerer, 2001; Campagnolo et al., 2002). Indeed, environmental concentrations of antibiotics are now often observed within the range of biological and evolutionary significance (Chow et al., 2021). Such intense and ubiquitous antibiotic selection surrounding dense human populations is likely driving the accumulation of diverse ARGs by these plasmids across the globe. Further, it is likely that each plasmid can accrue different suites of resistance genes, reflecting the specific conditions and selection pressures in their particular local environment. Thus, our findings highlight this plasmid family as a serious threat to successful bacterial infection control and human health, particularly for nosocomial infections. These findings also add to the growing concern that mega-plasmids are key disseminators of antibiotic resistance and require wide-spread surveillance (Cazares et al., 2020).
We have characterized a family of multi-drug resistance Acinetobacter mega-plasmids. We present evidence that these plasmids are highly transmissible, have the capacity for dispersing at a global-scale, can accumulate a vast cargo of niche-adaptive accessory genes, and can confer multi-drug resistance phenotypes of significant concern for human health. These findings add to a growing idea in the literature that points to the importance of viewing the antibiotic resistance crisis from a mobile element-centred outlook as opposed to a host cell-centric point of view (Lopatkin et al., 2017; Stevenson et al., 2017; Ghaly and Gillings, 2018; Pinilla-Redondo et al., 2018). Thus, expanding our therapeutic focus to consider ways of interfering with the transmission of mobile DNAs may be provide an additional, fruitful means of combating the spread of antibiotic resistance (Ghaly and Gillings, 2018).
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
TG contributed to sample processing and collection, contributed to the conception and design of analyses, performed the analyses, wrote the original draft of manuscript, and contributed to the final editing of manuscript. IP contributed to the conception and design of analyses and contributed to the final editing of manuscript. AS contributed to sample processing and collection and contributed to the final editing of manuscript. ST contributed to the conception and design of analyses and contributed to the final editing of manuscript. MG contributed to sample processing and collection, contributed to the conception and design of analyses, performed funding acquisition, and contributed to the final editing of manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by Australian Research Council Discovery Grants DP0877196 and DP130103839.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Jon Iredell for the Acinetobacter nosocomialis WM98B isolate. TG would like to thank Mary and Saoirse Ghaly for loving support.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.605952/full#supplementary-material
Footnotes
References
Albarracín, V. H., Pathak, G. P., Douki, T., Cadet, J., Borsarelli, C. D., Gärtner, W., et al. (2012). Extremophilic Acinetobacter strains from high-altitude lakes in Argentinean Puna: remarkable UV-B resistance and efficient DNA damage repair. Orig. Life Evol. Biosph. 42, 201–221. doi: 10.1007/s11084-012-9276-3
Alcock, B. P., Raphenya, A. R., Lau, T. T., Tsang, K. K., Bouchard, M., Edalatmand, A., et al. (2020). CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 48, D517–D525. doi: 10.1093/nar/gkz935
Almasaudi, S. B. (2018). Acinetobacter spp. as nosocomial pathogens: epidemiology and resistance features. Saudi J. Biol. Sci. 25, 586–596. doi: 10.1016/j.sjbs.2016.02.009
Ankenbrand, M. J., Hohlfeld, S., Hackl, T., and Förster, F. (2017). AliTV—interactive visualization of whole genome comparisons. PeerJ Comput. Sci. 3:e116. doi: 10.7717/peerj-cs.116
Antunes, L., Visca, P., and Towner, K. J. (2014). Acinetobacter baumannii: evolution of a global pathogen. Pathog. Dis. 71, 292–301. doi: 10.1111/2049-632X.12125
Arredondo-Alonso, S., Willems, R. J., van Schaik, W., and Schürch, A. C. (2017). On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data. Microb Genom. 3:e000128. doi: 10.1099/mgen.0.000128
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Babakhani, S., and Oloomi, M. (2018). Transposons: the agents of antibiotic resistance in bacteria. J. Basic Microbiol. 58, 905–917. doi: 10.1002/jobm.201800204
Bach, H., and Gutnick, D. L. (2006). Novel polysaccharide—protein-based amphipathic formulations. Appl. Microbiol. Biotechnol. 71, 34–38. doi: 10.1007/s00253-005-0149-9
Beaber, J. W., Hochhut, B., and Waldor, M. K. (2004). SOS response promotes horizontal dissemination of antibiotic resistance genes. Nature 427, 72–74. doi: 10.1038/nature02241
Botelho, J., Lood, C., Partridge, S. R., van Noort, V., Lavigne, R., Grosso, F., et al. (2019). Combining sequencing approaches to fully resolve a carbapenemase-encoding megaplasmid in a Pseudomonas shirazica clinical strain. Emerg. Microbes Infect. 8, 1186–1194. doi: 10.1080/22221751.2019.1648182
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Campagnolo, E. R., Johnson, K. R., Karpati, A., Rubin, C. S., Kolpin, D. W., Meyer, M. T., et al. (2002). Antimicrobial residues in animal waste and water resources proximal to large-scale swine and poultry feeding operations. Sci. Total Environ. 299, 89–95. doi: 10.1016/s0048-9697(02)00233-4
Castro-Jaimes, S., Bello-López, E., Velázquez-Acosta, C., Volkow-Fernández, P., Lozano-Zarain, P., Castillo Ramírez, S., et al. (2020). Chromosome architecture and gene content of the emergent pathogen Acinetobacter haemolyticus. Front. Microbiol. 11:926. doi: 10.3389/fmicb.2020.00926
Cazares, A., Moore, M. P., Hall, J. P., Wright, L. L., Grimes, M., Emond-Rhéault, J. -G., et al. (2020). A megaplasmid family driving dissemination of multidrug resistance in Pseudomonas. Nat. Commun. 11:1370. doi: 10.1038/s41467-020-15081-7
Cerezales, M., Xanthopoulou, K., Wille, J., Krut, O., Seifert, H., Gallego, L., et al. (2020). Mobile genetic elements harboring antibiotic resistance determinants in Acinetobacter baumannii isolates from Bolivia. Front. Microbiol. 11:919. doi: 10.3389/fmicb.2020.00919
Chow, L. K., Ghaly, T. M., and Gillings, M. R. (2021). A survey of sub-inhibitory concentrations of antibiotics in the environment. J. Environ. Sci. 99, 21–27. doi: 10.1016/j.jes.2020.05.030
Chusri, S., Chongsuvivatwong, V., Rivera, J. I., Silpapojakul, K., Singkhamanan, K., McNeil, E., et al. (2014). Clinical outcomes of hospital-acquired infection with Acinetobacter nosocomialis and Acinetobacter pittii. Antimicrob. Agents Chemother. 58, 4172–4179. doi: 10.1128/AAC.02992-14
Contreras-Moreira, B., and Vinuesa, P. (2013). GET_HOMOLOGUES, a versatile software package for scalable and robust microbial pangenome analysis. Appl. Environ. Microbiol. 79, 7696–7701. doi: 10.1128/AEM.02411-13
Cury, J., Jové, T., Touchon, M., Néron, B., and Rocha, E. P. (2016). Identification and analysis of integrons and cassette arrays in bacterial genomes. Nucleic Acids Res. 44, 4539–4550. doi: 10.1093/nar/gkw319
Delihas, N. (2008). Small mobile sequences in bacteria display diverse structure/function motifs. Mol. Microbiol. 67, 475–481. doi: 10.1111/j.1365-2958.2007.06068.x
Domingues, S., Nielsen, K. M., and da Silva, G. J. (2011). The blaIMP-5-carrying integron in a clinical Acinetobacter baumannii strain is flanked by miniature inverted-repeat transposable elements (MITEs). J. Antimicrob. Chemother. 66, 2667–2668. doi: 10.1093/jac/dkr327
Fondi, M., Bacci, G., Brilli, M., Papaleo, M. C., Mengoni, A., Vaneechoutte, M., et al. (2010). Exploring the evolutionary dynamics of plasmids: the Acinetobacter pan-plasmidome. BMC Evol. Biol. 10:59. doi: 10.1186/1471-2148-10-59
Fournier, P. -E., Vallenet, D., Barbe, V., Audic, S., Ogata, H., Poirel, L., et al. (2006). Comparative genomics of multidrug resistance in Acinetobacter baumannii. PLoS Genet. 2:e7. doi: 10.1371/journal.pgen.0020007
Garcia-Garcera, M., Touchon, M., Brisse, S., and Rocha, E. P. (2017). Metagenomic assessment of the interplay between the environment and the genetic diversification of Acinetobacter. Environ. Microbiol. 19, 5010–5024. doi: 10.1111/1462-2920.13949
Geoghegan, J. L., Duchêne, S., and Holmes, E. C. (2017). Comparative analysis estimates the relative frequencies of co-divergence and cross-species transmission within viral families. PLoS Pathog. 13:e1006215. doi: 10.1371/journal.ppat.1006215
Ghaly, T. M., Chow, L., Asher, A. J., Waldron, L. S., and Gillings, M. R. (2017). Evolution of class 1 integrons: mobilization and dispersal via food-borne bacteria. PLoS One 12:e0179169. doi: 10.1371/journal.pone.0179169
Ghaly, T. M., Geoghegan, J. L., Tetu, S. G., and Gillings, M. R. (2020). The peril and promise of integrons: beyond antibiotic resistance. Trends Microbiol. 28, 455–464. doi: 10.1016/j.tim.2019.12.002
Ghaly, T. M., and Gillings, M. R. (2018). Mobile DNAs as ecologically and evolutionarily independent units of life. Trends Microbiol. 26, 904–912. doi: 10.1016/j.tim.2018.05.008
Gillings, M. R. (2017). Class 1 integrons as invasive species. Curr. Opin. Microbiol. 38, 10–15. doi: 10.1016/j.mib.2017.03.002
Gillings, M. R., Labbate, M., Sajjad, A., Giguère, N. J., Holley, M. P., and Stokes, H. (2009). Mobilization of a Tn402-like class 1 integron with a novel cassette array via flanking miniature inverted-repeat transposable element-like structures. Appl. Environ. Microbiol. 75, 6002–6004. doi: 10.1128/AEM.01033-09
Gillings, M. R., and Stokes, H. (2012). Are humans increasing bacterial evolvability? Trends Ecol. Evolution 27, 346–352. doi: 10.1016/j.tree.2012.02.006
Gillings, M., Westoby, M., and Ghaly, T. (2018). Pollutants that replicate: xenogenetic DNAs. Trends Microbiol. 26, 975–977. doi: 10.1016/j.tim.2018.08.003
Göttig, S., Gruber, T. M., Higgins, P. G., Wachsmuth, M., Seifert, H., and Kempf, V. A. (2014). Detection of pan drug-resistant Acinetobacter baumannii in Germany. J. Antimicrob. Chemother. 69, 2578–2579. doi: 10.1093/jac/dku170
Hamidian, M., Ambrose, S. J., and Hall, R. M. (2016). A large conjugative Acinetobacter baumannii plasmid carrying the sul2 sulphonamide and strAB streptomycin resistance genes. Plasmid 87, 43–50. doi: 10.1016/j.plasmid.2016.09.001
Hamidian, M., Holt, K. E., Pickard, D., Dougan, G., and Hall, R. M. (2014). A GC1 Acinetobacter baumannii isolate carrying AbaR3 and the aminoglycoside resistance transposon TnaphA6 in a conjugative plasmid. J. Antimicrob. Chemother. 69, 955–958. doi: 10.1093/jac/dkt454
Hamidian, M., Nigro, S. J., and Hall, R. M. (2012). Variants of the gentamicin and tobramycin resistance plasmid pRAY are widely distributed in Acinetobacter. J. Antimicrob. Chemother. 67, 2833–2836. doi: 10.1093/jac/dks318
Hawkey, P. M., and Munday, C. J. (2004). Multiple resistance in gram-negative bacteria. Rev. Med. Microbiol. 15, 51–61. doi: 10.1097/01.revmedmi.0000131427.13353.7f
Ho, S. Y., Duchêne, S., and Duchêne, D. (2015). Simulating and detecting autocorrelation of molecular evolutionary rates among lineages. Mol. Ecol. Resour. 15, 688–696. doi: 10.1111/1755-0998.12320
Hoang, D. T., Chernomor, O., von Haeseler, A., Minh, B. Q., and Vinh, L. S. (2018). UFBoot2: improving the ultrafast bootstrap approximation. Mol. Biol. Evol. 35, 518–522. doi: 10.1093/molbev/msx281
Hunt, M., De Silva, N., Otto, T. D., Parkhill, J., Keane, J. A., and Harris, S. R. (2015). Circlator: automated circularization of genome assemblies using long sequencing reads. Genome Biol. 16:294. doi: 10.1186/s13059-015-0849-0
Jones, P., Binns, D., Chang, H. -Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kaas, R. S., Friis, C., Ussery, D. W., and Aarestrup, F. M. (2012). Estimating variation within the genes and inferring the phylogeny of 186 sequenced diverse Escherichia coli genomes. BMC Genomics 13:577. doi: 10.1186/1471-2164-13-577
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K., von Haeseler, A., and Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285
Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010
Klein, E. Y., van Boeckel, T. P., Martinez, E. M., Pant, S., Gandra, S., Levin, S. A., et al. (2018). Global increase and geographic convergence in antibiotic consumption between 2000 and 2015. Proc. Natl. Acad. Sci. 115, E3463–E3470. doi: 10.1073/pnas.1717295115
Kolmogorov, M., Yuan, J., Lin, Y., and Pevzner, P. A. (2019). Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 37, 540–546. doi: 10.1038/s41587-019-0072-8
Kristensen, D. M., Kannan, L., Coleman, M. K., Wolf, Y. I., Sorokin, A., Koonin, E. V., et al. (2010). A low-polynomial algorithm for assembling clusters of orthologous groups from intergenomic symmetric best matches. Bioinformatics 26, 1481–1487. doi: 10.1093/bioinformatics/btq229
Kümmerer, K. (2001). Drugs in the environment: emission of drugs, diagnostic aids and disinfectants into wastewater by hospitals in relation to other sources—a review. Chemosphere 45, 957–969. doi: 10.1016/S0045-6535(01)00144-8
Lean, S. S., and Yeo, C. C. (2017). Small, enigmatic plasmids of the nosocomial pathogen, Acinetobacter baumannii: good, bad, who knows? Front. Microbiol. 8:1547. doi: 10.3389/fmicb.2017.01547
Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100. doi: 10.1093/bioinformatics/bty191
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics 25, 1754–1760. doi: 10.1093/bioinformatics/btp324
Li, L., Stoeckert, C. J., and Roos, D. S. (2003). OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189. doi: 10.1101/gr.1224503
Lopatkin, A. J., Meredith, H. R., Srimani, J. K., Pfeiffer, C., Durrett, R., and You, L. (2017). Persistence and reversal of plasmid-mediated antibiotic resistance. Nat. Commun. 8:1689. doi: 10.1038/s41467-017-01532-1
Madden, T. (2013). “The BLAST sequence analysis tool” in The NCBI handbook [internet]. 2nd Edn. eds. J. Beck, D. Benson, J. Coleman, M. Hoeppner, M. Johnson, and D. Maglott (Bethesda, MD, USA: National Center for Biotechnology Information (US)).
Merino, M., Acosta, J., Poza, M., Sanz, F., Beceiro, A., Chaves, F., et al. (2010). OXA-24 carbapenemase gene flanked by XerC/XerD-like recombination sites in different plasmids from different Acinetobacter species isolated during a nosocomial outbreak. Antimicrob. Agents Chemother. 54, 2724–2727. doi: 10.1128/AAC.01674-09
Nguyen, L. -T., Schmidt, H. A., von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Nigro, S. J., and Hall, R. M. (2017). A large plasmid, pD46-4, carrying a complex resistance region in an extensively antibiotic-resistant ST25 Acinetobacter baumannii. J. Antimicrob. Chemother. 72, 3496–3498. doi: 10.1093/jac/dkx287
Pal, C., Bengtsson-Palme, J., Rensing, C., Kristiansson, E., and Larsson, D. J. (2014). BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 42, D737–D743. doi: 10.1093/nar/gkt1252
Pesesky, M. W., Tilley, R., and Beck, D. A. (2019). Mosaic plasmids are abundant and unevenly distributed across prokaryotic taxa. Plasmid 102, 10–18. doi: 10.1016/j.plasmid.2019.02.003
Pinilla-Redondo, R., Cyriaque, V., Jacquiod, S., Sørensen, S. J., and Riber, L. (2018). Monitoring plasmid-mediated horizontal gene transfer in microbiomes: recent advances and future perspectives. Plasmid 99, 56–67. doi: 10.1016/j.plasmid.2018.08.002
Quinlan, A. R., and Hall, I. M. (2010). BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. doi: 10.1093/bioinformatics/btq033
Rambaut, A. (2017). FigTree-version 1.4.3, a graphical viewer of phylogenetic trees. Available at: http://tree.bio.ed.ac.uk/software/figtree (Accessed October 5, 2016).
Revell, L. J. (2012). Phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217–223. doi: 10.1111/j.2041-210x.2011.00169.x
Rivera-Izquierdo, M., Valero-Ubierna, C., and Masso-Guijarro, P. (2020). Nosocomial bacteraemia caused by Acinetobacter ursingii in neonatology critical care unit in Spain: a rising challenge. Pediatr. Infect. Dis. J. 39, e28–e29. doi: 10.1097/INF.0000000000002562
Sajjad, A., Holley, M. P., Labbate, M., Stokes, H., and Gillings, M. R. (2011). Preclinical class 1 integron with a complete Tn402-like transposition module. Appl. Environ. Microbiol. 77, 335–337. doi: 10.1128/AEM.02142-10
Salgado-Camargo, A. D., Castro-Jaimes, S., Gutierrez-Rios, R. -M., Lozano, L. F., Altamirano-Pacheco, L., Silva-Sanchez, J., et al. (2020). Structure and evolution of Acinetobacter baumannii plasmids. Front. Microbiol. 11:1283. doi: 10.3389/fmicb.2020.01283
Salto, I. P., Tejerizo, G. T., Wibberg, D., Pühler, A., Schlüter, A., and Pistorio, M. (2018). Comparative genomic analysis of Acinetobacter spp. plasmids originating from clinical settings and environmental habitats. Sci. Rep. 8:7783. doi: 10.1038/s41598-018-26180-3
Sambrook, J., and Russell, D. W. (2006). Purification of nucleic acids by extraction with phenol:chloroform. Cold Spring Harb Protoc. 2006:pdb.prot4455. doi: 10.1101/pdb.prot4455
Schreiber, F., Szekat, C., Josten, M., Sahl, H. -G., and Bierbaum, G. (2013). Antibiotic-induced autoactivation of IS256 in Staphylococcus aureus. Antimicrob. Agents Chemother. 57, 6381–6384. doi: 10.1128/AAC.01585-13
Seemann, T. (2018). ABRicate: mass screening of contigs for antimicrobial and virulence genes. Department of Microbiology and Immunology, The University of Melbourne, Melbourne, Australia. Available at: https://github.com/tseemann/abricate (Accessed October 5, 2016).
Stevenson, C., Hall, J. P. J., Harrison, E., Wood, A., and Brockhurst, M. A. (2017). Gene mobility promotes the spread of resistance in bacterial populations. ISME J. 11, 1930–1932. doi: 10.1038/ismej.2017.42
Stokes, H. W., and Gillings, M. R. (2011). Gene flow, mobile genetic elements and the recruitment of antibiotic resistance genes into gram-negative pathogens. FEMS Microbiol. Rev. 35, 790–819. doi: 10.1111/j.1574-6976.2011.00273.x
Sullivan, M. J., and van Bakel, H. (2019). Chromatiblock: scalable whole-genome visualization of structural differences in prokaryotes. bioRxiv [Preprint]. doi: 10.1101/800920
Tettelin, H., Masignani, V., Cieslewicz, M. J., Donati, C., Medini, D., Ward, N. L., et al. (2005). Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial “pan-genome.” Proc. Natl. Acad. Sci. 102, 13950–13955. doi: 10.1073/pnas.0506758102
The Gene Ontology Consortium (2019). The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 47, D330–D338. doi: 10.1093/nar/gky1055
Valenzuela, J. K., Thomas, L., Partridge, S. R., Van Der Reijden, T., Dijkshoorn, L., and Iredell, J. (2007). Horizontal gene transfer in a polyclonal outbreak of carbapenem-resistant Acinetobacter baumannii. J. Clin. Microbiol. 45, 453–460. doi: 10.1128/JCM.01971-06
Vinuesa, P., and Contreras-Moreira, B. (2015). “Robust identification of orthologues and paralogues for microbial pan-genomics using GET_HOMOLOGUES: a case study of pIncA/C plasmids” in Bacterial pangenomics: Methods and protocols (methods in molecular biology). eds. A. Mengoni, M. Galardini, and M. Fondi (New York, USA: Springer Nature), 203–232.
Vinuesa, P., Ochoa-Sánchez, L. E., and Contreras-Moreira, B. (2018). GET_PHYLOMARKERS, a software package to select optimal orthologous clusters for phylogenomics and inferring pan-genome phylogenies, used for a critical geno-taxonomic revision of the genus Stenotrophomonas. Front. Microbiol. 9:771. doi: 10.3389/fmicb.2018.00771
Walker, B. J., Abeel, T., Shea, T., Priest, M., Abouelliel, A., Sakthikumar, S., et al. (2014). Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One 9:e112963. doi: 10.1371/journal.pone.0112963
Weber, B. S., Harding, C. M., and Feldman, M. F. (2016). Pathogenic Acinetobacter: from the cell surface to infinity and beyond. J. Bacteriol. 198, 880–887. doi: 10.1128/JB.00906-15
Wolf, Y. I., Makarova, K. S., Yutin, N., and Koonin, E. V. (2012). Updated clusters of orthologous genes for Archaea: a complex ancestor of the Archaea and the byways of horizontal gene transfer. Biol. Direct 7:46. doi: 10.1186/1745-6150-7-46
Wood, D. E., and Salzberg, S. L. (2014). Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 15:R46. doi: 10.1186/gb-2014-15-3-r46
Keywords: multi-drug resistance, pathogens, mobile genetic elements, nosocomial infections, plasmid pangenome
Citation: Ghaly TM, Paulsen IT, Sajjad A, Tetu SG and Gillings MR (2020) A Novel Family of Acinetobacter Mega-Plasmids Are Disseminating Multi-Drug Resistance Across the Globe While Acquiring Location-Specific Accessory Genes. Front. Microbiol. 11:605952. doi: 10.3389/fmicb.2020.605952
Edited by:
Santiago Castillo Ramírez, National Autonomous University of Mexico, MexicoReviewed by:
Xiaoting Hua, Zhejiang University, ChinaOlga Maria Perez-Carrascal, Université de Montréal, Canada
Copyright © 2020 Ghaly, Paulsen, Sajjad, Tetu and Gillings. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Timothy M. Ghaly, dGltb3RoeS5naGFseUBtcS5lZHUuYXU=