<figure><img src="image_1.png" /></figure>

HAL
open science

# Analysis of COMPASS, a New Comprehensive Plasmid Database Revealed Prevalence of Multireplicon and Extensive Diversity of IncF Plasmids

Pierre-Emmanuel Douarre, Ludovic Mallet, Nicolas Radomski, Arnaud Felten, Michel-Yves Mistou

► To cite this version:

Pierre-Emmanuel Douarre, Ludovic Mallet, Nicolas Radomski, Arnaud Felten, Michel-Yves Mistou. Analysis of COMPASS, a New Comprehensive Plasmid Database Revealed Prevalence of Multi-replicon and Extensive Diversity of IncF Plasmids. Frontiers in Microbiology, 2020, 11, pp.483. 10.3389/fmicb.2020.00483. anses-03821437

HAL Id: anses-03821437

https://anses.hal.science/anses-03821437

Submitted on 19 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Analysis of COMPASS, a New Comprehensive Plasmid Database Revealed Prevalence of Multireplicon and Extensive Diversity of IncF Plasmids

Pierre-Emmanuel Douarre¹*, Ludovic Mallet¹, Nicolas Radomski¹, Arnaud Felten¹ and Michel-Yves Mistou²*

¹ Agence Nationale de Sécurité Sanitaire de l'Alimentation, de l'Environnement et du Travail, Laboratory for Food Safety, Paris, France, ² INRAE, MalAGE, Université Paris-Saclay, Jouy-en-Josas, France

## OPEN ACCESS

**Edited by:**
Vasco Ariston De Carvalho
Azevedo,
Federal University of Minas Gerais,
Brazil

**Reviewed by:**
Dariusz Bartosik,
University of Warsaw, Poland
Richard John Meinersmann,
US National Poultry Research Center
(USDA-ARS), United States
Kristina Kadlec,
Independent Researcher, Wunstorf,
Germany

*Correspondence:
Pierre-Emmanuel Douarre
pierre-emmanuel.douarre@anses.fr
Michel-Yves Mistou
michel-yves.mistou@inrae.fr

**Specialty section:**
This article was submitted to
Evolutionary and Genomic
Microbiology,
a section of the journal
Frontiers in Microbiology

**Received:** 08 October 2019
**Accepted:** 05 March 2020
**Published:** 24 March 2020

**Citation:**
Douarre P-E, Mallet L, Radomski N, Felten A and Mistou M-Y (2020) Analysis of COMPASS, a New Comprehensive Plasmid Database Revealed Prevalence of Multireplicon and Extensive Diversity of IncF Plasmids. Front. Microbiol. 11:483. doi: 10.3389/fmicb.2020.00483

Plasmids are genetic elements that enable rapid adaptation and evolution by transferring genes conferring selective advantages to their hosts. Conjugative plasmids are predominantly responsible for the global dissemination of antimicrobial resistance, representing an important threat to global health. As the number of plasmid sequences grows exponentially, it becomes critical to depict the global diversity and decipher the distribution of circulating plasmids in the bacterial community. To this end, we created COMPASS, a novel and comprehensive database compiling 12,084 complete plasmids with associated metadata from 1571 distinct species isolated worldwide over more than 100 years. The curation of the database allowed us to identify identical plasmids across different bacteria revealing mainly intraspecies dissemination and rare cases of horizontal transmission. We outlined and analyzed all relevant features, plasmid properties, host range and characterized their replication and mobilization systems. After an exhaustive comparison of PlasmidFinder and MOB-typer, the MOB-typer-based analysis revealed that the current knowledge embedded in the current typing schemes fails to classify all the plasmid sequences collected in COMPASS. We were able to categorize 6828 and 5229 plasmids by replicon and MOB typing, respectively, mostly associated with *Proteobacteria* and *Firmicutes*. We then searched for the presence of multiple core genes involved in replication and propagation. Our results showed that 2403 plasmids carried multiple replicons that were distributed in 206 bacterial species. The co-integration of replicon types from different incompatibility (Inc) groups is an adaptive mechanism, which plays an important role in plasmid survival and dissemination by extending their host range. Our results highlight the crucial role of IncF alleles (present in 56% of all multireplicons) and revealed that IncH, IncR, and IncU replicons were also frequently carried in multireplicons. Here, we provided a comprehensive picture of the different IncF subtypes by identifying 20 different profiles in 849 IncF multireplicons, which were mostly associated with *Enterobacteriaceae*. These results could provide the basis for a novel IncF plasmid nomenclature based on different allelic profiles.

**Keywords:** plasmid database, host range, replicon typing, MOB typing, plasmid mobility, multireplicon, IncF plasmid
# INTRODUCTION

Plasmids are extrachromosomal fragments of DNA that can replicate with different degrees of autonomy from the host's replicative proteins and can transfer between bacterial species. They are found in the three domain of life and are widely represented in bacterial species. Plasmids can either be linear or circular, vary in size, and can represent a large proportion of the whole bacterial genome (Smillie et al., 2010; Shintani et al., 2015). These mobile genetic elements have a typical modular structure consisting of a conserved backbone region ensuring their survival or propagation and a variable accessory region that may encode host-beneficial traits (Garcillán-Barcia et al., 2011). The acquisition or deletion of adaptive modules from different phylogenetic origins begets enormous diversity among plasmids (de Toro et al., 2014; Pesesky et al., 2019). Plasmids are a keystone of horizontal gene transfer, facilitating the rapid evolution and adaptation of their hosts under changing environmental conditions. Plasmids contribute to the epidemic propagation of antibiotic resistance genes in bacterial pathogens and the dissemination of multidrug resistant plasmids is an increasing threat for modern medicine (Carattoli, 2013; von Wintersdorff et al., 2016). Some plasmids can be transferred between different bacterial species from diverse taxonomic groups and thus have a significant ecological impact. Understanding the genetic plasticity and transmission routes of plasmids is crucial in fighting against antimicrobial resistant pathogens.

The classification of plasmids using typing schemes based on replication or mobilization apparatus (replicon or MOB typing) is helpful to describe plasmid diversity but their lack of resolution limits their use for epidemiological studies (Orlek et al., 2017c).

The term “replicon” can refer to any DNA sequence that is capable of replicating as a unit such as a plasmid or can refer more precisely, like in this study to specific regions or genes encoding functions that enable the replication (Carattoli et al., 2014). Plasmids sharing one or more elements of the plasmid replication (same origin of replication or same replicon) or partitioning systems cannot stably coexist in a cell together and are defined as incompatible. Competition for replication factors also lead to competition between plasmids; meaning that plasmid having growth advantages rapidly outgrow other plasmids in the cell (Velappan et al., 2007). Traditionally, replicon typing classifies plasmids according to incompatibility (Inc) groups.

Plasmid mobility is determined by the propagation module that contains the mobilization system (relaxase and origin of transfer—*oriT*) and the mating pore formation (MPF) system. During the conjugation process, the relaxase recognizes the *oriT* (a short sequence of tens to hundreds of base pairs) and cuts the plasmid DNA at a conserved *nic* site to initiate the transfer into recipient cells while the MPF build the mating channel (Smillie et al., 2010). Thanks to their relaxase-specific properties, it is possible to assign a MOB-group to *oriT* sequences. The presence/absence of these three components allows the prediction of the plasmid mobility (Robertson and Nash, 2018).

While Whole Genome Sequencing offers a high-resolution method for studying the genetic structure and variability of plasmids, full exploitation of this massively available data is still a hurdle. Identifying and tracking circulating plasmids heavily relies on the ability to accurately assemble their whole genome sequences, which remains particularly challenging from current short reads mainly due to the presence of numerous repetitive elements, especially in large plasmids (>50 kb) (Arredondo-Alonso et al., 2017). Even though several bioinformatics tools can be applied to reconstruct plasmid sequences from short reads, a contiguous assembly is still difficult to obtain (Zhou and Xu, 2010; Lanza et al., 2014; Antipov et al., 2016; Rozov et al., 2017). This pitfall can be overcome by using long read sequencing technology such as *Oxford Nanopore* or *PacBio Technologies* (Hoffmann et al., 2017; Margos et al., 2017; Li et al., 2018).

With the advances of sequencing technologies and bioinformatics tools, the number of plasmid sequences available in public databases is growing rapidly. However, comparing newly sequenced plasmids with existing resources is cumbersome because of partial and misannotated sequences as well as a paucity in bioinformatic tools to explore the large amount of available data. The creation of a curated and comprehensive complete plasmid database is critical to integrate and elicit understanding of the available data. The first attempt was made by Orlek et al. (2017a), but the database was restricted to complete plasmids belonging to the *Enterobacteriaceae* family (Orlek et al., 2017a). More recently, a database containing all the complete plasmids from the NCBI's RefSeq database was integrated in the visualization and exploration tool pATLAS to facilitate the identification of plasmid from assemblies (Jesus et al., 2019). A similar database called PLSDB was also implemented in a user friendly and interactive web server, which offers the possibility to upload new sequences and query against the database (Galata et al., 2019). Both PLSDB and pATLAS databases allow the browsing and filtering of plasmids with their associated metadata and annotations. There are also bioinformatics tools that use plasmid databases as a reference for plasmid detection and identification. For example, the dataset that was created by Robertson and Nash (2018) and implemented in the MOB-suite software to reconstruct draft genome assemblies using known plasmid sequences.

In the present paper, we unveil the construction and the curation of a Comprehensive and Complete Plasmid Database we named “COMPASS.” This new database constituted of 12 084 plasmid sequences expands previously described resources and provides a full description of its content in terms of relevant properties (size, taxonomy, geographical origin, and year of isolation). With the objective of getting a global overview of the diversity of circulating plasmids, we characterized the database by *in silico* typing based on replication and mobilization systems and predicted the mobility from the analysis of dissemination markers. The distribution in the bacterial community was carefully examined to highlight particular features and specific associations between core genes and the plasmid host range.
# MATERIALS AND METHODS

## Creation and Curation of the COMPASS Database

Complete bacterial plasmid accessions were retrieved from the NCBI nucleotide database on 17 June 2018 using advanced optimized filters in the entrez query to exclude incomplete and misannotated plasmids. More precisely, keywords such as “plasmid,” “bacteria,” and “complete” were used in the initial query to retrieve complete plasmid records. Specific terms in the title were filtered to remove plasmid “genes, cds, origin of replication or region,” and other words such as “integron, transposon, operon, and phage” were applied to get rid of mobile genetic elements other than plasmids. Fasta and GenBank files associated with all plasmid accessions were downloaded from the NCBI database. A python script was developed to extract metadata from the GenBank files and to curate the database as described below. Duplicated records were first identified based on specific metadata. All the plasmids with the same size, topology, taxonomy, and header description were considered identical. Duplicated records were removed keeping preferentially the most recent version and RefSeq entry. Duplicated sequences sharing 100% identity and coverage were then identified using the clustering tool CD-HIT-EST (Li and Godzik, 2006). This step allowed the identification of clusters constituted of identical plasmids. One plasmid sequence was kept (as a reference) in the database while the other records were discarded. Finally, putative chromosomal sequences mislabeled as plasmids were identified by performing *in silico* Ribosomal Multilocus Sequence Typing (rMLST) analysis (Jolley et al., 2012). rMLST is an approach which indexes variation of the 53 genes encoding the bacterial ribosome protein subunits (*rps* genes) to identify rapidly the phylogenetic position of any bacterial sequence. The 53 *rps* genes, downloaded from PubMLST<fnref n="1" /> (12 February 2019) were sought in the plasmid records using BLASTn and hits with at least 95% identity and 95% coverage were considered potential matches (Jolley et al., 2012). The scripts and database of COMPASS developed in the present study can be found in the following GitHub repository: https://github.com/itsmeludo/COMPASS.

## Plasmid Typing and Predicted Mobility

*In silico* predictions of the replication gene (replicon), mobilization protein (relaxase), mate-pair formation (MPF), and the origin of transfer (*oriT*) types were obtained using the MOB-typer module from the MOB-suite package<fnref n="2" /> (Robertson and Nash, 2018). Replicon, relaxase, MPF, and *oriT* sequences were identified by BLAST using default parameters against reference databases that Robertson and Nash built for the MOB-typer tool. These databases contained 2481 plasmid-derived replication genes that corresponded to 1683 replicon types, 916 relaxase sequences classified into six MOB families (MOB<sub>F</sub>, MOB<sub>H</sub>, MOB<sub>Q</sub>, MOB<sub>C</sub>, MOB<sub>P</sub>, and MOB<sub>V</sub>), 2066 MPF proteins clustered in four groups (MPF<sub>I</sub>, MPF<sub>F</sub>, MPF<sub>G</sub>, and MPF<sub>T</sub>), and 502 *oriT* sequences.

The plasmid mobility was predicted based on the presence-absence of the three dissemination markers. A plasmid was labeled as “conjugative” if it carried a relaxase plus an MPF and as “mobilizable” if it contained at least a relaxase or an *oriT* whereas plasmids missing both a relaxase and an *oriT* were classified as “non-mobilizable” (Robertson and Nash, 2018).

Replicon typing was also performed using the conventional *in silico* tool PlasmidFinder (version October/2019)<fnref n="3" /> with recommended parameters (Carattoli et al., 2014) and the results were compared to those obtained with MOB-typer. The replicon content of the MOB-typer database was also compared to the one of PlasmidFinder using the clustering tool CD-HIT-EST-2D (i.e., 100% identity and 100% coverage) to identify specific alleles shared between the two typing resources.

## Visualization and Statistical Analysis

The plasmid size was analyzed with R scripts (R_Core_Team, 2018) and graphically represented with the ggplot2 library (Wickham, 2016). A custom world map showing the distribution and the number of plasmids isolated per country was created using the online tool Gunnmap 2⁴. The exploration of the taxonomical composition of the COMPASS database was performed using the web-based interactive visualization tool Krona (Ondov et al., 2011). The HTML link for the interactive Krona chart can be found in the COMPASS GitHub repository mentioned above. The plasmid content of our database was compared to other plasmid resources (PLSDB, pATLAS, and MOB-suite) by comparing the nucleotide databases 2 by 2 and identifying strictly identical plasmid sequence (i.e., 100% identity and 100% coverage) using the clustering tool CD-HIT-EST-2D. The integrative venny tool (Bardou et al., 2014)<fnref n="5" /> was also used to visualize the number of plasmids shared between COMPASS and the other three databases by comparing the lists of accessions. Two co-occurrence network analyses were performed independently based on the R igraph library to highlight plasmid mobilization in different species and decipher the linkages between replication genes present in multireplicon plasmids. Sankey diagrams were built using the online tool SankeyMATIC<fnref n="6" /> to show particular associations between typing results and taxonomy.

# RESULTS

## Construction and Curation of the COMPASS Database

A total of 22,425 sequences representing complete plasmids were retrieved from the NCBI nucleotide database. Duplicate records were first identified based on metadata; records presenting identical size, topology, taxonomy, and header were considered as duplicates. These filtering steps led to the identification of 9299 duplicated records. Among these, 6329 arose from

<footnote n="1">¹https://pubmlst.org/rmlst/</footnote>

<footnote n="2">²https://github.com/phac-nml/mob-suite</footnote>

<footnote n="3">³https://cge.cbs.dtu.dk/services/PlasmidFinder/</footnote>

<footnote n="4">⁴http://lert.co.nz/map/</footnote>

<footnote n="5">⁵http://jvenn.toulouse.inra.fr/app/index.html</footnote>

<footnote n="6">⁶http://sankeymatic.com/</footnote>
different databases (RefSeq vs GenBank, e.g., NZ_CP018683.2 vs CP018683), 22 were different versions of the same sequence (updated vs older version (e.g., CP022061.2 vs CP022061.1), and 2948 had different accession numbers (CP004864.1 vs CP007611.1). The list of redundant records removed from the initial list is available in **Supplementary Table S1**. The remaining 13,126 plasmids records were subjected to the stringent clustering tool CD-HIT-EST to remove strictly identical sequence (i.e., 100% identity and 100% coverage).

In total, 11,237 singletons (i.e., unique sequences) and 849 clusters containing from two (89%) to 25 identical plasmids were detected during this clustering step (Supplementary Table S2).

The analysis of these clusters through a manual curation revealed that for 547 clusters, the plasmids within each cluster were associated with the same bacterial strains. These duplicated records were not detected by the first curation step because the headers of their sequences were different (e.g., “complete_sequence” vs “complete_genome”; “mega_plasmid” vs “megaplasmid”; “virulence_plasmid_pCP301” vs “plasmid_pCP301”...). On the other hand, 47 clusters were missing information regarding the bacterial strains and could not be categorized. Interestingly, for 234 clusters, plasmids within each cluster were associated with the same species but different strains, emphasizing the strong species-specific plasmid dissemination. In contrast, rare cases of horizontal transmission between different species were also observed in 21 clusters. These clusters were represented in **Supplementary Figure S1**, where nodes are bacterial species and edges indicate the sharing of identical plasmid. The result highlights that most of the subgraphs (n = 19) group together bacterial species belonging to the same family. For example, the 33 kb plasmid of the cluster #7470 that was isolated in the United States from three different Enterobacteria (*Enterobacter hormaechei*, *Klebsiella oxytoca*, and *Enterobacter cloacae*) over 3 years (2010–2012). However, rare transfers between more distantly related organisms can also occur; like the 46.5 kb plasmid found in two bacteria from different orders (*Acinetobacter lwoffii* and *Klebsiella aerogenes*, cluster #6503). One reference plasmid was kept within each 849 clusters for downstream analyses and 1040 redundant sequences were discarded (**Supplementary Tables S1, S2**). Finally, chromosomal sequences were identified based on the number of RPS genes present. Out of the 12,086 plasmids, 108 contained less than 10 RPS genes while two records NZ_CP014062.1 and NZ_CP022019.1 had 55 and 54 exact matches with RPS genes of the rMLST database. The predicted taxa of these two records were *Salmonella enterica* and *Pseudomonas putida*, respectively. These two chromosomal sequences were excluded from the database that finally contains 12,084 plasmids. These two records have since been updated on the current GenBank database as “chromosome, complete genome,” which confirmed the present rMLST results.

## Comparison of Plasmid Databases

Our COMPASS database shared 11,374 plasmids (94%) with PLSDB (Galata et al., 2019), 9873 plasmids (82%) with pATLAS (Jesus et al., 2019), and 8271 plasmids (68%) with the MOB-suite database, respectively (Robertson and Nash,

2018) (**Supplementary Figure S2 and Supplementary Table S3**). Altogether, 6440 plasmids were shared between COMPASS, and the other three databases. The content of the three published databases was also compared with COMPASS to confirm the overlap between plasmid resources. MOB-suite, PLSDB, and pATLAS shared 11,610 (96%), 11,380 (83%), and 10,204 (77%) plasmids with COMPASS, respectively (**Supplementary Table S3**). The differences observed among the number of plasmids shared between two specific databases (for example, COMPASS versus MOB-suite and MOB-suite versus COMPASS) revealed the presence of duplicates. The presence of identical sequences within a database was indeed confirmed (by running CD-HIT-EST 100%) for MOB-suite (n = 3500), pATLAS (n = 381), and PLSDB (n = 7) but not for COMPASS. The list of duplicated records is available in **Supplementary Table S3**. The differences in composition could be explained by the construction and curation steps that were specific to each database. The keywords used for the original query to download plasmid records and the inclusion/exclusion criteria in the filtering steps were diverse. Duplicated records were also removed differently; redundant sequences were identified based on pair-wise distances created by the program mash for the PLSDB database while Robertson and Nash removed Refseq accessions, which contained the accession of another records (**Supplementary Methods** of the MOB-suite paper).

## Description of the COMPASS Database

In order to generate a comprehensive database, all the metadata available from the GenBank files were extracted for each plasmid and carefully examined. The information regarding the size, the topology, and the complete lineage was available for all the plasmids whereas the other metadata were partially available. Relevant metadata describing the host range and transmission route are presented below. All the other metadata can be found in **Supplementary Table S4**.

Plasmids ranged in size from 744 to 2,555,069 bp, with a median at 53,206 bp. The plasmid size distribution varied among the different phyla (**Figure 1A**). Large plasmids occurred more frequently in *Deinococcus-Thermus*, *Proteobacteria*, and *Actinobacteria* than the rest of the phyla (**Figure 1A**). The highest frequency of large plasmids was observed for *Deinococcus-Thermus*, while the largest plasmids were found in *Proteobacteria*. Indeed, 132 plasmids were larger than 1 Mb and 72% of these “megaplasmids” were isolated from *alphaproteobacteria*. The plasmids isolated from *Chlamydia* and *Tenericutes* were overall smaller than the rest of the plasmids with a median at 7504 and 3990 bp, respectively. With regard to the first and third quartile close to the median, the size of the chlamydial plasmids was very homogenous. Indeed, 70 chlamydial plasmids (92% of the phylum) exhibited a size ranging between 7415 and 7553 bp. This group included the plasmids pCTA, pCHL1, pSW2 from *Chlamydia trachomatis*, pMoPn from *Chlamydia muridarum*, and pCpGP1 from *Chlamydia caviae* (Zhong, 2017).

The metadata “Country” locating the isolation was available for 7096 plasmids (59%), tallying 126 countries with a frequency ranging from one plasmid collected in 21 different
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Description of the COMPASS database. (A) Violin plot displaying plasmid size distribution (log₁₀) among the main phyla (n = 12,084). The phylum entitled "Other" is composed of 57 plasmids from 12 minority phyla (n < 25) and 129 plasmids missing taxonomy data. (B) Krona plot showing the compositions of taxa and taxonomic ranks (n = 12,084).</figcaption></figure>

countries to 1703 plasmids collected in the United States (Supplementary Figure S3).

The metadata “Collection date” was obtained for 5819 plasmids (48%). Plasmids sequences were obtained from isolates dating from 1884 to 2018 with an exponential increase after 2000, accounting for 85% of all plasmid present in the COMPASS database (Supplementary Figure S3).

The vast majority (94%) of the plasmids were circular and among the 665 plasmids that were linear, 48% were isolated from bacteria belonging to the *Borreliaceae* family.

Our plasmid database covered 1571 distinct bacterial species, 443 genera, 189 families, 93 orders, 37 classes, and 21 phyla. *Proteobacteria* (63%) and *Firmicutes* (21%) were the most represented phyla and 32% of the bacterial species (n = 3963) belonged to the *Enterobacteriaceae* family (Figure 1B). Overall, *Escherichia coli* (n = 1672), *Klebsiella pneumoniae* (n = 1053), and *S. enterica* (n = 477) plasmids were the most prevalent species. The taxonomy distribution of COMPASS was compared at the phylum level to that of the PLSDB database and we did not find a significant difference between both distributions (Kolmogorov-Smirnov test, *p*-value = 0.997).

## Replicon Typing

Using MOB-typer, the presence of at least one replicon was detected for 9231 plasmids (76.4%) (Supplementary Table S5). Due to biases related to the constitutions of GenBank and replicon databases, the number of plasmids encoding replicons varied regarding the taxonomy (Supplementary Figure S4). For instance, 81.8% of the plasmids associated with *Firmicutes* and *Proteobacteria* contained at least one replicon whereas only 29.6 and 18.2% of plasmids from *Cyanobacteria* and *Spirochaetes* matched the replicon database, respectively. These

observations stressed the knowledge gap of plasmid biology associated with less represented taxa in sequence databases. The number of replicons detected within a plasmid varied from a single gene identified in 6828 plasmids to seven replicons identified in a 446 kb plasmid isolated from *Raoultella ornithinolytica* (Supplementary Figure S5). Interestingly, the amount of replicons within a plasmid does not correlate with the plasmid size (Pearson test: rho = -2%, *p*-value = 6.8 × 10⁻³); meaning that big plasmids do not necessarily possess several replication genes.

Among the 6828 plasmids specifically typed by the presence of a single replicon, 2469 plasmids harbored a known replicon from the ColRNAI family (n = 639) or were associated with a known Inc group (n = 1830) while 4359 plasmids carried a replicon that did not belong to a recognized family or Inc group. These plasmids were associated with 1436 replicon clusters from the MOB-typer database. While 19 Inc groups were identified, plasmids belonging to Inc10 (n = 2), IncT (n = 1), and IncW (n = 4) were rare (Table 1). The most frequently assigned Inc groups by MOB-typer were IncI (362 records) followed by IncX (309 records), IncF (also described as F-type replicons) (269 records), IncA/C (176 records), and IncP (158 records). The plasmids harboring a replicon from the Inc10, Inc18, IncF, IncH, IncI and IncL/M, IncR, IncT, IncU, IncW, IncX, and IncY were associated with a limited host range (less than five bacterial families) while ColRNAI, Inc11, Inc13, IncA/C, IncN, IncP, and IncQ plasmids were present in a broad host range (BHR) (14, 8, 11, 5, 20, and 7 bacterial families, respectively) (Figure 2). The IncL/M family has now been split into IncL, IncM1, and IncM2 subtypes (Carattoli et al., 2015). In addition, the IncA/C group have recently been shown to be compatible and are now distinguished to IncA and IncC (Ambrose et al., 2018).
<table><caption>TABLE 1 | Summary of plasmid replicon type features.</caption><thead><tr><th>Replicon type</th><th>No. of plasmids</th><th>Min (bp)</th><th>Max (bp)</th><th>Median (bp)</th><th>No. of families</th><th>Host range</th><th>Major family host (%)</th><th>Predicted mobility (%)</th><th>No. of plasmids</th><th>No. of relaxases</th><th>Major MOB type (%)</th></tr></thead><tbody><tr><td>ColRNAI</td><td>639</td><td>1308</td><td>51,662</td><td>5202</td><td>14</td><td>BHR</td><td>Enterobacteriaceae (81.69)</td><td>Mobilizable (80.13)</td><td>302</td><td>5</td><td>MOBp (49.01)</td></tr><tr><td>Inc10</td><td>2</td><td>1742</td><td>1813</td><td>1778</td><td>1</td><td>Narrow</td><td>Bacillaceae (100)</td><td>Non-mobilizable (100)</td><td>0</td><td>NA</td><td>NA</td></tr><tr><td>Inc11</td><td>44</td><td>1278</td><td>15,360</td><td>3908</td><td>8</td><td>BHR</td><td>Streptococcaceae (43.18)</td><td>Non-mobilizable (54.55)</td><td>20</td><td>2</td><td>MOBv (55)</td></tr><tr><td>Inc13</td><td>58</td><td>1643</td><td>11,801</td><td>3300</td><td>11</td><td>BHR</td><td>Bacillaceae (31.03)</td><td>Non-mobilizable (56.90)</td><td>22</td><td>1</td><td>MOBv (100)</td></tr><tr><td>Inc18</td><td>12</td><td>5842</td><td>70,706</td><td>26,725</td><td>3</td><td>Narrow</td><td>Enterococcaceae (50)</td><td>Non-mobilizable (58.33)</td><td>5</td><td>1</td><td>MOBv (100)</td></tr><tr><td>IncAVC</td><td>176</td><td>3645</td><td>233,057</td><td>158,754</td><td>8</td><td>BHR</td><td>Enterobacteriaceae (77.84)</td><td>Conjugative (92.61)</td><td>160</td><td>2</td><td>MOBH (98.75)</td></tr><tr><td>IncF*</td><td>1118</td><td>19,788</td><td>416,444</td><td>109,349</td><td>4</td><td>Narrow</td><td>Enterobacteriaceae (94)</td><td>Conjugative (68.25)</td><td>861</td><td>4</td><td>MOBF (82.11)</td></tr><tr><td>IncH</td><td>19</td><td>78,444</td><td>223,698</td><td>97,393</td><td>1</td><td>Narrow</td><td>Enterobacteriaceae (100)</td><td>Non-mobilizable (63.16)</td><td>7</td><td>2</td><td>MOBp (71.43)</td></tr><tr><td>Incl</td><td>362</td><td>9466</td><td>207,960</td><td>83,635</td><td>2</td><td>Narrow</td><td>Enterobacteriaceae (99.45)</td><td>Conjugative (93.09)</td><td>328</td><td>2</td><td>MOBp (99.39)</td></tr><tr><td>IncL/M</td><td>95</td><td>46,885</td><td>95,855</td><td>70,092</td><td>4</td><td>Narrow</td><td>Enterobacteriaceae (89.47)</td><td>Conjugative (91.58)</td><td>87</td><td>1</td><td>MOBp (100)</td></tr><tr><td>IncN</td><td>151</td><td>20,225</td><td>17,4695</td><td>54,242</td><td>5</td><td>BHR</td><td>Enterobacteriaceae (89.40)</td><td>Conjugative (97.35)</td><td>147</td><td>2</td><td>MOBF (97.96)</td></tr><tr><td>IncP</td><td>158</td><td>3269</td><td>1,499,175</td><td>60,505</td><td>20</td><td>BHR</td><td>Enterobacteriaceae (17.09)</td><td>Conjugative (76.58)</td><td>122</td><td>3</td><td>MOBp (95.90)</td></tr><tr><td>IncQ</td><td>53</td><td>6388</td><td>208,409</td><td>8300</td><td>7</td><td>BHR</td><td>Enterobacteriaceae (43.40)</td><td>Mobilizable (90.57)</td><td>48</td><td>4</td><td>MOBQ (45.83)</td></tr><tr><td>IncR</td><td>65</td><td>18,990</td><td>125,961</td><td>61,010</td><td>1</td><td>Narrow</td><td>Enterobacteriaceae (98.46)</td><td>Non-mobilizable (64.62)</td><td>223</td><td>2</td><td>MOBp (95.65)</td></tr><tr><td>IncT</td><td>1</td><td>83,698</td><td>83,698</td><td>83,698</td><td>1</td><td>Narrow</td><td>Enterobacteriaceae (100)</td><td>Mobilizable (100)</td><td>1</td><td>1</td><td>MOBp (100)</td></tr><tr><td>IncU</td><td>22</td><td>7995</td><td>84,749</td><td>38,938</td><td>3</td><td>Narrow</td><td>Aeromonadaceae (45.45)</td><td>Conjugative (50)</td><td>13</td><td>1</td><td>MOBp (100)</td></tr><tr><td>IncW</td><td>4</td><td>4233</td><td>39,924</td><td>38,971</td><td>3</td><td>Narrow</td><td>Enterobacteriaceae (50)</td><td>Conjugative (75)</td><td>4</td><td>2</td><td>MOBF (75)</td></tr><tr><td>IncX</td><td>309</td><td>6464</td><td>76,500</td><td>43,380</td><td>3</td><td>Narrow</td><td>Enterobacteriaceae (98.71)</td><td>Conjugative (95.79)</td><td>296</td><td>2</td><td>MOBp (99.66)</td></tr><tr><td>IncY</td><td>30</td><td>56,460</td><td>126,046</td><td>96,897</td><td>1</td><td>Narrow</td><td>Enterobacteriaceae (100)</td><td>Non-mobilizable (90)</td><td>3</td><td>2</td><td>MOBp (99.66)</td></tr></tbody></table>

<footnote>*IncF plasmids include all the plasmids carrying F-type replicons; single replicon (n = 269) and multireplicon (n = 849).</footnote>
However, the MOB-typer database does not differentiate between the subtypes within these two groups and therefore our analysis did not take these new distinctions into account.

As expected, Inc grouping schemes corresponded to families for which they were designed to perform plasmid typing. More than 75% of ColRNAI, IncA/C, IncF, IncH, IncI, IncL/M, IncN, IncR, IncT, IncX, and IncY plasmids were associated with *Enterobacteriaceae* family whereas Inc10-11-13-18 replicons most frequently occurred in *Firmicutes* (*Streptococcaceae*, *Bacillaceae*, *Enterococcaceae*) (Table 1 and Figure 2). Replicons clusters were detected in 146 different bacterial families and were the most prevalent in Gram-positive bacteria (393 *Bacillaceae*, 331 *Lactobacillaceae*, and 315 *Staphylococcaceae*). We observed that the plasmid size varied greatly among the different types (Figure 3). More than 80% of the plasmids from the Inc10 (n = 2), Inc11 (n = 44), Inc13 (n = 58), and the ColRNAI (n = 639) families were smaller than 10 kb (median at 1778, 3908, 3300, and 5202 bp, respectively) while 80% of the IncA/C (n = 176) and IncH (n = 19) plasmids were larger than 80 kb (median at 158,754 and 97,393 bp, respectively). Figure 3 also shows that the size of IncP (n = 158) and IncF (n = 269) plasmids were very diverse (3269 to 1,499,175 and 19,788 to 416,444 bp, respectively).

Multiple replication genes were detected in 2403 plasmids. All the known replicon types from the ColRNAI family and Inc groups were detected except for Inc10, Inc13, and IncW. The most frequent replicons were IncFII (24.5%), IncFIB (n = 12.4%), IncFIIA (10.3%), and IncFIA (3.8%) and 56% of the multireplicons (n = 1349) carried one of the four IncF alleles. For these 1349 plasmids, a co-occurrence network to visualize the frequently associated replicons highlights the central role of IncF in the multireplicons (Supplementary Figure S6). Overall, multireplicons plasmids were found in 206 bacterial species belonging to 40 bacterial families but were mostly represented in *Enterobacteriaceae* (65.9%) and in Gram-positive bacteria such as *Staphylococcaceae* (7.9%), *Bacillaceae* (3.8%), and *Lactobacillaceae* (3.58%). The combination of replicon types resulted in 539 different profiles of which 138 profiles (1189 plasmids) contained only known replicons. Among these 1189 plasmids, 849 carried only a combination of IncF replicons (20 profiles), 268 plasmids had an IncF that coexisted with another replicon from other Inc groups (84 profiles), and 72 plasmids did not harbor any IncF replicons (37 profiles). The profiles IncU-IncP, IncU-IncX, and IncA/C-IncP were predominant in the plasmids that did not possess IncF while ColRNAI, IncI, IncR, and IncQ were the replicons most frequently associated with IncF. A truncated IncQ replicon inside a resistance region carried by large IncF plasmids has been previously observed by Partridge et al. (2018). In contrast to these 340 plasmids that could not be assigned a typical type, the 849 multireplicons that only contained IncF alleles were categorized as IncF plasmids.

## IncF Plasmids

To gain further insights into the diversity of IncF plasmids, we characterized all the combination of IncF alleles in single and multireplicons and analyzed their host range. Overall, IncF plasmids (n = 1118) were detected in 50 bacterial species

mainly from *Enterobacteriaceae* (94%) and multireplicons were 3.5 times more prevalent than plasmids carrying a single IncF allele (Figure 4).

Among single replicons, IncFIIA (n = 40) and IncFIA (n = 17) replicon had a limited host range (present only in two and six species, respectively) compared to IncFIB (n = 159) and IncFII (n = 53) that were found in 20 and 23 species. IncFIIA and IncFIA occurred significantly more often in *E. coli* whereas IncFIB and IncFII were dominant in *K. pneumoniae* and *S. enterica* respectively (Figure 4). Similar results were previously found by using the original PCR-based replicon typing (PBRT) and the updated replicon sequence typing (RST) schemes (Carattoli et al., 2005; Villa et al., 2010). In contrast to the single replicons, IncFII was detected in almost all the multireplicons (n = 833) while IncFIB, IncFIIA, and IncFIA were present in 537, 478, and 131 plasmids, respectively.

The most prevalent profile was IncFII-IncFIIA (n = 233), followed by IncFII(2)-IncFIB (n = 181) and IncFII-IncFIIA-IncFIB (n = 151). The IncF multireplicons associated with the three most prevalent species *E. coli* (47%), *K. pneumoniae* (23%), and *S. enterica* (11%) were very diverse; belonging to 13, 10, and 8 different profiles, respectively. Conversely, some plasmids from other species were associated with a particular type. All the plasmids from *Shigella* (n = 62) were linked to the same profile IncFII-IncFIIA and the nine replicons from *Yersinia pestis* belonged to the IncFII(2)-IncFIB profile (Figure 4).

## Comparison Between MOB-Typer and PlasmidFinder

Using the conventional tool PlasmidFinder, replicons were only detected in 4931 plasmids (40%), of which 3349 have been assigned to specific types (Supplementary Table S6). The MOB-typer results of the 6828 plasmids associated with single replicons as well as the 849 IncF multireplicons were compared to the replicon types obtained with PlasmidFinder (Supplementary Table S6). We found a high concordance (> 90% match) between the known replicon types identified by MOB-typer and those detected by PlasmidFinder for IncA/C, IncI, IncL/M, IncN, IncQ, IncR, IncT, IncX, and IncY plasmids. However, discrepancies were observed for the other Inc groups. The 116 plasmids belonging to Inc10, 11, Inc13, and Inc18 types were either assigned to different types (n = 53) or were not detected (n = 63) by PlasmidFinder. Similarly, 84% (n = 16) of the IncH plasmids were assigned to a different replicon type (p0111). Finally, only 43% (n = 68) of the IncP plasmids were confirmed while 8% (n = 12) were classified as different types and 49% (n = 78) were not detected by PlasmidFinder. Regarding the 1118 plasmids identified by MOB-typer as IncF, PlasmidFinder also detected IncF replicons in 1058 plasmids (94.6%), but both methods identified different combination of IncF alleles. We observed a high correlation for the presence of IncFII (94.4%), IncFIA (99.2%), and IncFIB (98.3%) alleles. Nevertheless, PlasmidFinder identified the IncFIC allele (AP001918) in 37 plasmids that were not detected by MOB-typer and conversely the allele IncFIIA (AP014877) was identified by MOB-typer in 518 plasmids but not detected by PlasmidFinder. Among the 4359 plasmids that
<figure><img src="image_5.png" /></figure>

carry unknown replicons (replicon clusters), 3566 were not typed by PlasmidFinder while the other 793 were associated to 128 replicons types.

The comparison of the two replicon databases revealed that MOB-typer shared only 147 sequences with the PlasmidFinder and carry 2336 unique sequences (Supplementary Table S7). The analysis of the 147 shared replicons showed that 10 MOB-typer sequences associated with Inc11, Inc13, and Inc18 were identical to PlasmidFinder sequences that were classified as different types; thus clarifying the divergence mentioned above for Inc11-13-18 plasmids. Similarly for the IncH plasmids, the allele p0111 (AP010962) in PlasmidFinder was identical to IncH (CP021336) in MOB-typer. The comparison also showed that the IncFIIA and IncFIC alleles were unique in MOB-typer and PlasmidFinder, respectively, explaining the different allelic profiles found for IncF plasmids.

## MOB Typing

The second classification method is based on the variation of the relaxase protein. Using the relaxase database of MOB-typer, 821 mobilization protein alleles were detected in 5791 plasmids (47.9%) (Supplementary Table S5). A strong phylum-bias was observed for the detection of relaxases sequences. Indeed, 56.5

and 41.3% of plasmids from *Proteobacteria* and *Firmicutes*, respectively, harbored relaxase proteins but none was identified in *Cyanobacteria*, *Spirochaetes*, or *Tenericutes* (Supplementary Figure S4). While more than 25% of plasmids in the database contained multiple replicons (see above), a single relaxase protein was present in 90% (n = 5229) of the plasmids.

Among the 5229 plasmids that could be assigned a MOB type, the most prevalent was MOB~P~ (n = 2303) followed by MOB~F~ (n = 1422), MOB~Q~ (n = 717), MOB~H~ (n = 378), MOB~V~ (n = 269), and MOB~C~ (n = 106). In contrast to these known MOB, 34 relaxases did not belong to the 6 main families but matched sequences classified as “MOB_unknown” from the MOB-typer database. We found that the plasmids harboring a relaxase MOB~C~ or MOB~H~ had a limited host range including 9 and 14 bacterial families, respectively, and were mostly associated with *Enterobacteriaceae* (>80%). MOB~F~ plasmids, found in 43 families, were also highly represented in *Enterobacteriaceae* (67%) whereas the relaxase MOB~V~ were virtually absent in this family. In contrast to the other MOB groups, relaxase belonging to the MOB~P~ and MOB~Q~ were highly promiscuous (BHR plasmids in 80 and 47 bacterial families, respectively) (Figure 2).

Plasmid size trended within MOB types, we observed that 95% of MOB~H~-associated plasmids (n = 378) were larger than 110 kb
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Plasmid size (kb) distribution among the known replicon (n = 2469) and MOB types (n = 5195) from the COMPASS database.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Sankey diagram showing the associations between replicon types and replicon profiles with bacterial species in 1118 IncF plasmids. The number in brackets refers to the copy number of a particular IncF allele.</figcaption></figure>
while 91% of MOB<sub>V</sub> plasmids (n = 269) were smaller than 40 kb (Figure 3). The size of the other plasmids (MOB<sub>C</sub>, MOB<sub>P</sub>, and MOB<sub>Q</sub>) was more diversified.

Overall, the analysis of the co-occurrence of the replicons and relaxases in 1703 plasmids revealed that more than 95% of IncI, IncL/M, IncP, IncR, IncT, IncU, IncX, and IncY plasmids were associated with MOB<sub>P</sub> while IncA/C and IncN plasmids were linked to MOB<sub>H</sub> and MOB<sub>F</sub>, respectively. In contrast, Inc11, Inc13, and Inc18 plasmids found in *Firmicutes* encoded preferentially MOB<sub>V</sub>. The other replicon types such as ColRNAI, IncF, and IncQ were not associated with a specific type and carried 4 and 5 different relaxase families (Figure 5A and Table 1).

## Predicted Mobility

Plasmid mobility was predicted based on the presence/absence of three mobilization markers (relaxase, *oriT*, and MPF) (Supplementary Table S5) (Robertson and Nash, 2018).

Using the known *oriT* loci and their matching MOB group from the MOB-typer database, 2048 *oriT* sequences corresponding to the five MOB types present in the database (MOB<sub>C</sub>-associated *oriT* were missing) were detected in 1942 plasmids (16%). The vast majority of plasmids harbored a single *oriT* excepted for 104 plasmids carrying 2 and a single plasmid bearing 3. The most abundant *oriT* loci matched MOB<sub>P</sub> (n = 974) and MOB<sub>F</sub> (n = 482) groups while 96% of the associated plasmids belonged to *Proteobacteria* (n = 1624) and *Firmicutes* (n = 251). The predicted *oriT*-group matched their relaxase MOB-type for 72% of the plasmids, while discrepancies resulted from the absence of MOB (n = 411), a wrong match (n = 85), or the detection of an unknown type (n = 43).

Finally, 1831 MPF sequences were identified against the MOB-typer database in 4270 plasmids. In comparison to other typing components, the number of MPF per plasmid fluctuated considerably from 1 to 21 genes identified in 608 and 2 plasmids, respectively; 3662 plasmids carried more than 1 MPF and 19% of all MPF-carrying plasmids encoded more than 10 MPFs (Supplementary Figure S7). Regarding the typeability, all the sequences detected within a plasmid were from the same cluster therefore all the 4270 plasmids could be assigned to an MPF type. MPF<sub>T</sub> and MPF<sub>F</sub> were the most prevalent types (2053 and 1633 plasmids, respectively) while MPF<sub>I</sub> was detected in 495 plasmids and MPF<sub>G</sub> was rare (n = 13). Most of the MPF proteins were detected in *Proteobacteria* (86.2%) and *Firmicutes* (9.4%) while no MPF were identified in the phyla *Chlamydiae*, *Cyanobacteria*, *Spirochaetes*, and *Tenericutes*. Specific associations between MPF and replicon types were also observed. MPF<sub>F</sub> were strongly linked to IncA/C and IncN whereas MPF<sub>T</sub> and MPF<sub>I</sub> types were associated with IncX and IncL/M, respectively (Supplementary Figure S8).

Overall, looking at the three dissemination markers, we found that 5882 plasmids (49%) are missing both a relaxase and an *oriT* and were therefore classified as non-mobilizable; 2532 mobilizable plasmids (21%) carrying either a relaxase (n = 2121) or an *oriT* (n = 218) and 3670 conjugative plasmids (30%) harboring both MOB and MPF.

For the latter plasmids, a strong specificity between four MOB and two MPF types can be observed. Indeed, between 95 and 100% of MOB<sub>C</sub>, MOB<sub>Q</sub>, and MOB<sub>V</sub> were associated with MPF<sub>T</sub> while MOB<sub>H</sub> was linked to MPF<sub>F</sub>. All the MOB<sub>H</sub> plasmids (n = 378) and 82% of MOB<sub>F</sub> plasmids (n = 1175) were predicted to be conjugative whereas 97% of MOB<sub>V</sub> plasmids were mobilizable.

Regarding the replicon types, IncA/C, IncI, IncL/M, IncN, and IncX plasmids were mostly conjugative (>90% of plasmids predicted to be transferable by conjugation) while none of the ColRNAI, Inc10-11-13, and IncR plasmids were self-transferable. Most of the plasmids belonging to ColRNAI (80%) and IncQ (90%) types were mobilizable. In contrast to these transferable plasmids, 90% of the IncY plasmids were non-mobilizable (Figure 5B and Table 1).

All the plasmids belonging to *Chlamydiae*, *Spirochaetes*, and *Tenericutes*, were predicted to be non-mobilizable whereas the proportion of mobile plasmids (conjugative and mobilizable) varied from 59% in *Proteobacteria* to 15% in *Deinococcus-Thermus* (Supplementary Figure S9). Accordingly, conjugative plasmids were larger (median at 94,307 bp) than non-mobilizable plasmids (median at 36,766 bp), reflecting the genes needed to encode the conjugation functions and the ones conferring selective advantages (Supplementary Figure S10).

## DISCUSSION

### The Comprehensive and Complete Plasmid Database—COMPASS

Plasmids are ubiquitous in bacteria and have been shown to be involved in transferring genes important for survival and fitness to their hosts, enabling rapid adaptation to various selective pressures. They have been widely associated with the global dissemination of antimicrobial resistance and represent an important threat to global health (Chen et al., 2014; Hardiman et al., 2016). Despite the important contribution of plasmids to bacterial evolution, little is known about the global diversity of circulating plasmids and their host range. In order to fill this gap, we curated and compiled all the complete plasmid sequences and associated metadata that were available in the NCBI nucleotide archive to create a comprehensive database.

The initial retrieval of 22,425 plasmid records was followed by several curation steps based on their metadata and nucleotide sequences. The identification of redundant sequences using the CD-HIT-EST tool allowed us to identify clusters of identical plasmids. The analysis of these clusters showed that plasmid dissemination occurs more frequently within the same species (91.8%) rather than transferring horizontally (8.2%) between different species. In addition, these analyses showed that horizontal transmission is favored between species belonging to the same family.

Finally, our new collection composed of 12,084 plasmids covers 1571 distinct bacterial species that were isolated worldwide (126 countries) over more than 100 years (1884–2018). The diversity of the plasmid host range in our database was compared to that of 197,345 bacterial genomes (belonging to 63,240
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | (A) Sankey diagram showing the associations between replicon and MOB types in 1703 plasmids. (B) Associations between replicon types and predicted mobility in 2469 plasmids from the COMPASS database.</figcaption></figure>

different Taxonomy IDs) available in the NCBI genome report the 14 April 2019⁷. The host range of the plasmids in our database was found representative of the bacterial genome distribution present in NCBI (i.e., absence of significant differences based on the Kolmogorov-Smirnov test, *p*-value = 0.99). This result partially explained that for many complete plasmid sequences, a complete chromosome sequence from the same isolate may well have been deposited. Among the different bacterial families present in the COMPASS database, we observed an over-representation of plasmids isolated from *Enterobacteriaceae* (*n* = 3963), a family of particular interest for studying antimicrobial resistance as its members have become increasingly resistant to the antibiotics commonly used to treat bacterial infection (Carattoli, 2009; Iredell et al., 2016).

The comparison of COMPASS with previously existing databases showed a similar composition with PLSDB (94%) and pATLAS (82%) in comparison to MOB-suite that only shared 68% with our database. In contrast to the other plasmid resources, our database do not contain any duplicated records, demonstrating the strength of our curation steps.

## Plasmid Sizes and Bacterial Hosts Are Related

Plasmids covered a wide size range and both intra- and inter-phyla variations were observed. Indeed, the smallest (644 bp) and largest (2.5 Mb) plasmids were identified in *Proteobacteria* (*Candidatus Tremblaya* and *Paraburkholderia caribensis*, respectively). This phylum, especially the *Alphaproteobacteria*, harbored 94% of the megaplasmid (>1 Mb). In addition

of carrying a large chromosome of several megabases, the presence of megaplasmids is another feature that characterizes *Alphaproteobacteria* (Ettema and Andersson, 2009). In contrast, the smallest plasmids were detected in bacteria belonging to the *Tenericutes* and *Chlamydia* phyla, which also happened to bear the smallest chromosomes (Vetrovsky and Baldrian, 2013; Li and Du, 2014). In our database, almost all chlamydial plasmids (92%) corresponded to an non-conjugative 7.5 kb cryptic plasmid discovered 39 years ago by Lovett et al. (1980) which is involved in the pathogenicity and widely distributed among many different chlamydial strains (Zhong, 2017). These plasmids, although isolated from chlamydial species that infect different animal host species, from different countries and timeline are highly conserved. The lack of diversity of the chlamydial plasmids could be explained by the fact that the members of the *Chlamydia* phylum are a group of obligate intracellular bacteria therefore limiting the propensity for horizontal gene transfer (Blanc et al., 2007).

## Limitations of Current Plasmid Typing Methods

The large number of plasmids in our collection required a thorough classification to decipher plasmid diversity. Here, the *in silico* typing of all the plasmids in our database was performed using the recent reference-based tool MOB-typer (Robertson and Nash, 2018) and the conventional tool PlasmidFinder. MOB-typer provides replicon typing similar to PlasmidFinder but with the inclusion of transferability predictions based on the presence of relaxase, oriT, and MPF. Even though several studies investigated the diversity of plasmid mobilization proteins (Garcillán-Barcia et al., 2009, 2011; Fernandez-Lopez et al., 2017),

<footnote n="7">⁷ftp://ftp.ncbi.nlm.nih.gov/genomes/GENOME_REPORTS/</footnote>
MOB-typer was (at the time of the manuscript writing) the only tool providing *in silico* predictions of the plasmid mobilization system but a novel tool, MOBscan developed by Garcillan-Barcia et al. (2019) recently became available online.

The proportion of plasmids that can be assigned a replicon type is influenced by the dataset and the bioinformatics methods used. Orlek et al. (2017b) found that 85% of the curated *Enterobacteriaceae* plasmids (n = 2097) could be replicon-typed (Orlek et al., 2017b). In contrast, only 39.5% of the extensive PLSDB dataset (n = 13,789) were typed (Galata et al., 2019). These results both obtained by PlasmidFinder highlight the fact that the performance of the replicon typing is constrained by the bacterial host and most suitable for *Enterobacteriaceae* and Gram-positive species. The percentage of plasmids typed with PlasmidFinder in our study is very similar (40.8%) to the one obtained by Galata et al. (2019). In contrast, MOB-typer was able to type 76% of the plasmids present in COMPASS. Overall, few discrepancies were observed between plasmid assignments of PlasmidFinder and MOB-typer at the nucleotide level (e.g., IncH replicons). However, several plasmid replicons were identified by MOB-typer while they were not detected by PlasmidFinder (e.g., Inc10-11-Inc13-Inc18, IncFIIA, and IncP). The fact that MOB-typer database contains 2021 replicon sequences more than PlasmidFinder and that only 147 sequences are shared between the two databases could explain the observed differences. For example, MOB-typer database contains 34 different replicons associated with IncP while PlasmidFinder possesses only four. In addition, it must be emphasized that 56% of the IncP plasmids undetected by PlasmidFinder are larger than 100,000 bp. IncP plasmids have been reported to be small plasmids (Popowska and Krawczyk-Balska, 2013). However, our analysis showed that the size of IncP plasmid ranges from 3269 and 1,499,175 bp suggesting that some alleles of IncP replicon could be present in larger plasmids than previously reported.

These results highlight the need to fill gaps in plasmid classification, but also the importance of the completeness of plasmid database to have a global view of plasmid properties. With the continuous accumulation of sequencing data, it would be useful to set up a real-time, curated feeding system to maintain comprehensive typing databases. A clear limitation of the current plasmid classification scheme revealed by the present analysis of COMPASS is the frequent occurrence of multiple replicons, complicating the classification.

## Characterization of the COMPASS Database Using MOB-Typer

MOB-typer detected 9231 plasmids bearing replicons in COMPASS. Among these plasmids, 6828 plasmids carried a single replicon of which 2469 belonged to the ColRNAI family and 19 different Inc groups. Unlike replication genes, the presence of multiple relaxases within a plasmid was rare (4%) and the number of unknown relaxases was very low (n = 22). Consequently, MOB-typer was able to specifically classified 5207 plasmids into the six MOB families (MOBC-F-H-P-Q-V). The most frequent replicon families assigned by MOB-typer were

ColRNAI (n = 639), IncI (n = 362), and IncX (n = 309) while the most predominant MOB types were MOBP (n = 2303), MOBF (n = 1422), and MOBQ (n = 717). Overall, we found that half of the plasmids (51%) of the COMPASS database are predicted to be mobilizable or self-transmissible. The percentage of mobile plasmid in our study is higher than the one (39%) found by Smillie et al. (2010) suggesting that new relaxase and MPF sequences were identified in the last decade (Ramachandran et al., 2017).

Our results confirmed that the plasmid typeability differed among the different phyla and we found that plasmids from *Proteobacteria* and *Firmicutes* were the most frequently typed. Our findings regarding the taxonomic diversity of bacterial host among the different Inc groups correlated with previous studies where host ranges are generally narrow for IncF, IncH, IncI, IncT, and IncX or broad for IncA/C, IncP, and IncQ plasmids (Suzuki et al., 2010; Rozwandowicz et al., 2018). ColRNAI have been reported as narrow host range plasmids, i.e., cannot replicate in two different taxonomic classes (Smorawinska et al., 2012). However, the analysis of the COMPASS database revealed the presence of ColRNAI replicons in bacterial species belonging to five different taxonomic classes, although it was mostly found associated to *Gammaproteobacteria* (96%). This result suggests that ColRNAI have a broader replication range than previously reported and should be considered as putative BHR replicon.

Several studies previously showed that IncP plasmids can transfer and replicate in almost all Gram-negative bacteria. These conjugative plasmids are widely present in the environment and are also present in pathogenic and opportunistic bacteria. The presence of antibiotic resistance genes and their broad distribution raises more and more concern (Popowska and Krawczyk-Balska, 2013).

The analysis of COMPASS allowed to conclude at a wide scale on the preferential associations between replicons and relaxases (Garcillán-Barcia et al., 2011; Shintani et al., 2015; Rozwandowicz et al., 2018). IncI, IncL/M, IncP, and IncX replicons were mostly associated with MOBP while IncA/C and IncN were linked to MOBH and MOBF. In comparison to these plasmids that were predicted to be mainly conjugative, the replicons belonging to Inc10-11-13-18 were mostly identified on non-mobilizable plasmids harboring MOBV. Finally, the BHR ColRNAI and IncQ plasmids were not restricted to one specific relaxase and were mostly mobilizable.

Interestingly, we observed a strong size stratification associated with certain replicon-MOB types: Inc10, Inc11, and Inc13 plasmids carrying MOBV were small (median from 1778 to 3300 bp) while MOBH-associated IncA/C plasmids were large (median at 158,754 bp).

## Multireplicon Plasmids Are Widespread Among Circulating Plasmids

Multireplicon plasmids present a particular interest as they had been associated with both virulence and multidrug resistance (Carattoli et al., 2006; Hopkins et al., 2006; Johnson et al., 2007). Multiple replicons were identified in 19% of the characterized plasmids and were represented in
all major phyla. Their prevalence could indicate a propensity
to circulate between different hosts and thus be vectors
more prone to the dissemination of undesirable traits
such as antimicrobial resistance and virulence. Indeed,
it has already been described that a narrow host range
plasmid can broaden its host spectrum by cointegrating a
BHR replicon (Osborn et al., 2000; Rozwandowicz et al.,
2018). In such multireplicon plasmids, it was suggested
that one replicon is expressed due to the selective pressure
of plasmid replication while the other is free to diverge.
Thus, plasmids can evolve alternatively through replicon
acquisition and sequence divergence (Carattoli, 2009;
Villa et al., 2010).

Overall, a high degree of plasmid diversity was observed, with 539 different profiles occurring among the 2403 multireplicon plasmids examined. The analysis of these profiles revealed the crucial role of IncF replicon types. These replicons were present in more than half of the multirePLICONS where they coexist either with other IncF alleles (n = 849) or with other Inc groups (n = 500). The most frequent replicon associated with IncF was the BHR replicon ColRNAI, thus extending the narrow host spectrum of IncF plasmids.

MultirePLICONS harboring IncH, IncR, and IncU were also
common. Similarly to IncF, these three replicon types were
occurring more often in multirePLICONS (n = 159, 106, and 97)
than single replicons (n = 19, 65, and 22) and the most frequent
profile (IncU-IncP) associated a limited host range replicon
with a BHR replicon.

Extensive Diversity of IncF Plasmids and
Species-Specific Profiles

IncF plasmids represent one of the most widespread
plasmid types in clinical *Enterobacteriaceae* and have
been associated with the dissemination of relevant
antimicrobial resistance (Johnson et al., 2007; Villa et al.,
2010; Rozwandowicz et al., 2018).

Overall, 1118 plasmids (9%) belonging mostly to *Enterobacteriaceae* were classified as IncF plasmids and categorized into single (n = 269) and multirePLICons (n = 849). The occurrence of the four IncF alleles varied between the two categories and the biggest difference was observed for the IncFII replicon, which was present in 19 and 98% of the single and multirePLICons, respectively. It has been shown that IncFII are free to diverge when associated with FIA and FIB, evolving toward the formation of new compatible variants (Villa et al., 2010; Rozwandowicz et al., 2018). Thus, the acquisition and the sequence divergence of IncFII enable the replication in a broader host spectrum and contribute to the global dissemination of IncF. Preferential associations between bacterial hosts and IncF profiles were observed. For instance, *K. pneumoniae* and *S. enterica* both carried favorably IncFIB and IncFII in single replicons and these two alleles also coexisted in multirePLICons suggesting that these plasmids could have arisen from recombination of two single IncF plasmids. Regarding the mobilization system,

our results correlate with previous studies (Rozwandowicz
et al., 2018). While 68% were predicted to be conjugative, 82
and 85% of the IncF plasmids were associated with MOB<sub>F</sub>
and MPF<sub>F</sub>.

CONCLUSION

It has been shown that highly similar bacterial strains can
contain a compelling variety of different plasmids (Lanza
et al., 2014). The analysis of the diversity and host range
of circulating plasmids is therefore essential to characterize
bacterial isolates. Overall, *in silico* typing is likely to remain an
important tool for plasmid analysis, but their performance is
constrained by the diverse and dynamic nature of plasmids,
the plasmid dataset of interest, and the applied bioinformatics
methods. We showed the advantages of using MOB-typer
versus PlasmidFinder that currently does not characterize
the mobility of plasmids and only includes replicons from
*Enterobacteriaceae* and a few Gram-positive bacteria. However,
more efforts need to be done to develop a general classification
scheme for plasmids from all microbial lineages. COMPASS
plasmids could be easily integrated as a reference database
into pipelines for reconstructing plasmids and ascertaining
host range properties to new identified plasmids. The
complete classification will help deciphering the molecular
epidemiology of antimicrobial resistance by identifying
specific types and plasmid features linked to a particular
resistance gene.

DATA AVAILABILITY STATEMENT

The plasmid sequences of the COMPASS database compiled in
this project, the metadata associated for each plasmid, and the
HTML link for the interactive Krona chart are available online on
the following GitHub repository: https://github.com/itsmeludo/
COMPASS.

AUTHOR CONTRIBUTIONS

M-YM conceived and piloted the project. P-ED carried
out the analyses and wrote the manuscript. LM
developed the script for the creation and curation of the
database. LM, NR, and AF provided advice and help
for bioinformatics and genome analyses. M-YM, NR,
LM, and AF participated in the discussion and reviewed
the manuscript.

FUNDING

This study was supported by the French Agency for Food,
Environmental and Occupational Health and Safety (ANSES).
This article is also part of the European Joint Program One
Health EJP. This project has received funding from the European Union's Horizon 2020 Research and Innovation Program (H2020 EJP-AMR3-RaDAR) under Grant Agreement No. 773830.

## ACKNOWLEDGMENTS

We would like to thank the GAMeR (Génome Analyse Modélisation et Risque) team for bioinformatics development, advice on genomic analyses, computational resources, and reviewing the article. We thank Federica Palma for fruitful discussion.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.00483/full#supplementary-material

**FIGURE S1** | Network of identical plasmids isolated across different bacterial species from the COMPASS database (n = 54). Each line represents an identical plasmid shared between two bacterial species and the thickness of the line represent the number of plasmid cluster. Each color represents a bacterial family.

**FIGURE S2** | Comparison of the plasmid contents of four plasmid databases (COMPASS, PLSDB, pATLAS, and MOB-suite).

**FIGURE S3** | World map and histogram showing the distribution and the occurrence of plasmids isolated per countries and per year.

**FIGURE S4** | Detection of plasmid replicon and relaxase (in%) among the different phyla from the COMPASS database (n = 12,084). The phylum entitled "Other" is composed of 57 plasmids from 12 minority phyla (n < 25) and 129 plasmids missing taxonomy data.

**FIGURE S5** | Correlation between the number of replicons and the plasmid size from the COMPASS database (n = 9231).

**FIGURE S6** | Network of multiple replicon types (i.e., >2) isolated from different bacterial species from the COMPASS database (n = 1349). Each line represents shared replicon types harbored by plasmids and the thickness of the line represents the number of replicons types harbored by these latter.

**FIGURE S7** | Number of mating pore formation (MPF) proteins per plasmids.

**FIGURE S8** | Sankey diagram showing the associations between replicon and MPF types in 1 378 plasmids from the COMPASS database.

**FIGURE S9** | Percentage of transferable and non-mobilizable plasmids among the different phyla from the COMPASS database (n = 12,084).

**FIGURE S10** | Violin plot displaying plasmid size distribution (log₁₀) among non-mobilizable and transferable plasmids from the COMPASS database (n = 12,084).

**TABLE S1** | Curation of COMPASS database.

**TABLE S2** | CD-HIT clusters analysis.

**TABLE S3** | Comparison of plasmid databases.

**TABLE S4** | Description of COMPASS metadata.

**TABLE S5** | MOB-typer results.

**TABLE S6** | Comparison of PlasmidFinder and MOB-typer results.

**TABLE S7** | Comparison of PlasmidFinder and MOB-typer replicon databases.

## REFERENCES

Ambrose, S. J., Harmer, C. J., and Hall, R. M. (2018). Compatibility and entry exclusion of IncA and IncC plasmids revisited: IncA and IncC plasmids are compatible. *Plasmid* 9, 7–12. doi: 10.1016/j.plasmid.2018.02.002

Antipov, D., Hartwick, N., Shen, M., Raiko, M., Lapidus, A., and Pevzner, P. A. (2016). plasmidSPAdes: assembling plasmids from whole genome sequencing data. *Bioinformatics* 32, 3380–3387.

Arredondo-Alonso, S., Willems, R. J., Van Schaik, W., and Schurch, A. C. (2017). On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data. *Microb. Genom.* 3:e000128.

Bardou, P., Mariette, J., Escudié, F., Djemiel, C., and Klopp, C. (2014). jvenn an interactive Venn diagram viewer. *BMC Bioinformatics* 15:293. doi: 10.1186/1471-2105-15-293

Blanc, G., Ogata, H., Robert, C., Audic, S., Claverie, J. M., and Raoult, D. (2007). Lateral gene transfer between obligate intracellular bacteria: evidence from the Rickettsia massiliae genome. *Genome Res.* 17, 1657–1664. doi: 10.1101/gr.6742107

Carattoli, A. (2009). Resistance plasmid families in *Enterobacteriaceae*. *Antimicrob. Agents Chemother.* 53, 2227–2238. doi: 10.1128/aac.01707-08

Carattoli, A. (2013). Plasmids and the spread of resistance. *Int. J. Med. Microbiol.* 303, 298–304. doi: 10.1016/j.ijmm.2013.02.001

Carattoli, A., Bertini, A., Villa, L., Falbo, V., Hopkins, K. L., and Threlfall, E. J. (2005). Identification of plasmids by PCR-based replicon typing. *J. Microbiol. Methods* 63, 219–228. doi: 10.1016/j.mimet.2005.03.018

Carattoli, A., Miriagou, V., Bertini, A., Loli, A., Colinon, C., Villa, L., et al. (2006). Replicon typing of plasmids encoding resistance to newer beta-lactams.pdf. *Emerg. Infect. Dis.* 12, 145–148.

Carattoli, A., Seiffert, S. N., Schwendener, S., Perreten, V., and Endimiani, A. (2015). Differentiation of IncL and IncM plasmids associated with the spread of clinically relevant antimicrobial resistance. *PLoS ONE* 10:e0123063. doi: 10.1371/journal.pone.0123063

Carattoli, A., Zankari, E., Garcia-Fernandez, A., Voldby Larsen, M., Lund, O., Villa, L., et al. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. *Antimicrob. Agents Chemother.* 58, 3895–3903. doi: 10.1128/aac.02412-14

Chen, L., Chavda, K. D., Melano, R. G., Hong, T., Rojtman, A. D., Jacobs, M. R., et al. (2014). Molecular survey of the dissemination of two blaKPC-harboring IncFIA plasmids in New Jersey and New York hospitals. *Antimicrob. Agents Chemother.* 58, 2289–2294. doi: 10.1128/aac.02749-13

de Toro, M., Garcilla-Barcia, M. P., and De La Cruz, F. (2014). Plasmid diversity and adaptation analyzed by massive sequencing of *Escherichia coli* plasmids. *Microbiol. Spectr.* 2. doi: 10.1128/microbiolspec.PLAS-0031-2014

Ettema, T. J., and Andersson, S. G. (2009). The alpha-proteobacteria: the Darwin finches of the bacterial world. *Biol. Lett.* 5, 429–432. doi: 10.1098/rsbl.2008.0793

Fernandez-Lopez, R., Redondo, S., Garcillan-Barcia, M. P., and De La Cruz, F. (2017). Towards a taxonomy of conjugative plasmids. *Curr. Opin. Microbiol.* 38, 106–113. doi: 10.1016/j.mib.2017.05.005

Galata, V., Fehlmann, T., Backes, C., and Keller, A. (2019). PLSDB: a resource of complete bacterial plasmids. *Nucleic Acids Res.* 47, D195–D202.

Garcillán-Barcia, M. P., Alvarado, A., and De La Cruz, F. (2011). Identification of bacterial plasmids based on mobility and plasmid population biology. *FEMS Microbiol. Rev.* 35, 936–956. doi: 10.1111/j.1574-6976.2011.00291.x

Garcillán-Barcia, M. P., Francia, M. V., and De La Cruz, F. (2009). The diversity of conjugative relaxases and its application in plasmid classification. *FEMS Microbiol. Rev.* 33, 657–687. doi: 10.1111/j.1574-6976.2009.00168.x

Garcillán-Barcia, M. P., Redondo-Salvo, S., Vielva, L., and De La Cruz, F. (2019). "Horizontal gene transfer," in *Methods in Molecular Biology*, eds M. Boekels Gogarten, J. Peter Gogarten, and L. Olendzenski (Totowa, NJ: Humana Press), 2075.

Hardiman, C. A., Weingarten, R. A., Conlan, S., Khil, P., Dekker, J. P., Mathers, A. J., et al. (2016). Horizontal transfer of carbapenemase-encoding plasmids and
comparison with hospital epidemiology data. *Antimicrob. Agents Chemother.* 60, 4910-4919. doi: 10.1128/aac.00014-16

Hoffmann, M., Pettengill, J. B., Gonzalez-Escalona, N., Miller, J., Ayers, S. L., Zhao, S., et al. (2017). Comparative sequence analysis of multidrug-resistant inca/c plasmids from *Salmonella enterica*. *Front. Microbiol.* 8:1459. doi: 10.3389/fmicb. 2017.01459

Hopkins, K. L., Liebana, E., Villa, L., Batchelor, M., Threlfall, E. J., and Carattoli, A. (2006). Replicon typing of plasmids carrying CTX-M or CMY beta-lactamases circulating among *Salmonella* and *Escherichia coli* isolates. *Antimicrob. Agents Chemother.* 50, 3203-3206. doi: 10.1128/aac.00149-06

Iredell, J., Brown, J., and Tagg, K. (2016). Antibiotic resistance in *Enterobacteriaceae*: mechanisms and clinical implications. *BMJ* 352:h6420. doi: 10.1136/bmj.h6420

Jesus, T. F., Ribeiro-Goncalves, B., Silva, D. N., Bortolaia, V., Ramirez, M., and Carrico, J. A. (2019). Plasmid ATLAS: plasmid visual analytics and identification in high-throughput sequencing data. *Nucleic Acids Res.* 47, D188–D194.

Johnson, T. J., Wannemuehler, Y. M., Johnson, S. J., Logue, C. M., White, D. G., Doetkott, C., et al. (2007). Plasmid replicon typing of commensal and pathogenic *Escherichia coli* isolates. *Appl. Environ. Microbiol.* 73, 1976–1983. doi: 10.1128/aem.02171-06

Jolley, K. A., Bliss, C. M., Bennett, J. S., Bratcher, H. B., Brehony, C., Colles, F. M., et al. (2012). Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain. *Microbiology* 158, 1005–1015. doi: 10.1099/ mic.0.055459-0

Lanza, V. F., de Toro, M., Garcilla-N-Barcia, M. P., Mora, A., Blanco, J., Coque, T. M., et al. (2014). Plasmid Flux in *Escherichia coli* ST131 sublineages analyzed by plasmid constellation network (PLACNET), a new method for plasmid reconstruction from WGS. *PLoS Genet.* 10:e1004766. doi: 10.1371/journal.pgen. 1004766

Li, R., Xie, M., Dong, N., Lin, D., Yang, X., Wong, M. H. Y., et al. (2018). Efficient generation of complete sequences of MDR-encoding plasmids by rapid assembly of MINION barcoding sequencing data. *Gigascience* 7, 1–9.

Li, W., and Godzik, A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. *Bioinformatics* 22, 1658–1659. doi: 10.1093/bioinformatics/btl158

Li, X. Q., and Du, D. (2014). Variation, evolution, and correlation analysis of C+G content and genome or chromosome size in different kingdoms and phyla. *PLoS ONE* 9:e88339. doi: 10.1371/journal.pone.0088339

Lovett, M., Kuo, K.-K., Holmes, K., and Falkow, S. (1980). "Plasmids of the genus Chlamydia," in *Proceedings of the 11th International Congress of Chemotherapy and the 19th Interscience Conference on Antimicrobial Agents and Chemotherapy*, eds J. D. Nelson and C. Grassi (Washington, D.C: American Society for Microbiology).

Margos, G., Hepner, S., Mang, C., Marosevic, D., Reynolds, S. E., Krebs, S., et al. (2017). Lost in plasmids: next generation sequencing and the complex genome of the tick-borne pathogen *Borrelia burgdorferi*. *BMC Genomics* 18:422. doi: 10.1186/s12864-017-3804-5

Ondov, B. D., Bergman, N. H., and Phillippy, A. M. (2011). Interactive metagenomic visualization in a Web browser. *BMC Bioinformatics* 12:385. doi: 10.1186/1471-2105-12-385

Orlek, A., Phan, H., Sheppard, A. E., Doumith, M., Ellington, M., Peto, T., et al. (2017a). A curated dataset of complete *Enterobacteriaceae* plasmids compiled from the NCBI nucleotide database. *Data Brief* 12, 423–426. doi: 10.1016/j.dib. 2017.04.024

Orlek, A., Phan, H., Sheppard, A. E., Doumith, M., Ellington, M., Peto, T., et al. (2017b). Ordering the mob: Insights into replicon and MOB typing schemes from analysis of a curated dataset of publicly available plasmids. *Plasmid* 91, 42–52. doi: 10.1016/j.plasmid.2017.03.002

Orlek, A., Stoesser, N., Anjum, M. F., Doumith, M., Ellington, M. J., Peto, T., et al. (2017c). Plasmid classification in an era of whole-genome sequencing: application in studies of antibiotic resistance epidemiology. *Front. Microbiol.* 8:182. doi: 10.3389/fmicb.2017.00182

Osborn, A. M., Tatley, F. M. D. S., Steyn, L. M., Pickup, R. W., and Saunders, J. R. (2000). Mosaic plasmids and mosaic replicons evolutionary lessons from the analysis of genetic diversity in IncFII-related replicons. *Microbiology* 146, 2267–2275. doi: 10.1099/00221287-146-9-2267

Partridge, S. R., Kwong, S. M., Firth, N., and Jensen, S. O. (2018). Mobile genetic elements associated with antimicrobial resistance. *Clin. Microbiol. Rev.* 31: e88-17.

Pesesky, M. W., Tilley, R., and Beck, D. A. C. (2019). Mosaic plasmids are abundant and unevenly distributed across prokaryotic taxa. *Plasmid* 102, 10–18. doi: 10.1016/j.plasmid.2019.02.003

Popowska, M., and Krawczyk-Balska, A. (2013). Broad-host-range IncP-1 plasmids and their resistance potential. *Front. Microbiol.* 4:44. doi: 10.3389/fmicb.2013. 00044

Ramachandran, G., Miguel-Arribas, A., Abia, D., Singh, P. K., Crespo, I., Gago-Cordoba, C., et al. (2017). Discovery of a new family of relaxases in *Firmicutes bacteria*. *PLoS Genet.* 13:e1006586. doi: 10.1371/journal.pgen.100 6586

Robertson, J., and Nash, J. H. E. (2018). MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies. *Microb. Genom.* 4. doi: 10.1099/mgen.0.000206

Rozov, R., Brown Kav, A., Bogumil, D., Shterzer, N., Halperin, E., Mizrahi, I., et al. (2017). Recycler: an algorithm for detecting plasmids from de novo assembly graphs. *Bioinformatics* 33, 475–482.

Rozwandowicz, M., Brouwer, M. S. M., Fischer, J., Wagenaar, J. A., Gonzalez-Zorn, B., Guerra, B., et al. (2018). Plasmids carrying antimicrobial resistance genes in *Enterobacteriaceae*. *J. Antimicrob. Chemother.* 73, 1121–1137. doi: 10.1093/jac/ dkx488

R_Core_Team (2018). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing.

Shintani, M., Sanchez, Z. K., and Kimbara, K. (2015). Genomics of microbial plasmids: classification and identification based on replication and transfer systems and host taxonomy. *Front. Microbiol.* 6:242. doi: 10.3389/fmicb.2015. 00242

Smillie, C., Garcillan-Barcia, M. P., Francia, M. V., Rocha, E. P., and De La Cruz, F. (2010). Mobility of plasmids. *Microbiol. Mol. Biol. Rev.* 74, 434–452.

Smorawinska, M., Szulewska, M., Zaleski, P., Wawrzyniak, P., Maj, A., Plucienniczak, A., et al. (2012). Mobilizable narrow host range plasmids as natural suicide vectors enabling horizontal gene transfer among distantly related bacterial species. *FEMS Microbiol. Lett.* 326, 76–82. doi: 10.1111/j.1574- 6968.2011.02432.x

Suzuki, H., Yano, H., Brown, C. J., and Top, E. M. (2010). Predicting plasmid promiscuity based on genomic signature. *J. Bacteriol.* 192, 6045–6055. doi: 10.1128/jb.00277-10

Velappan, N., Sblattero, D., Chasteen, L., Pavlik, P., and Bradbury, A. R. (2007). Plasmid incompatibility: more compatible than previously thought? *Protein Eng. Des. Select.* 20, 309–313. doi: 10.1093/protein/gzm005

Vetrovsky, T., and Baldrian, P. (2013). The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses. *PLoS One* 8:e57923. doi: 10.1371/journal.pone.0057923

Villa, L., Garcia-Fernandez, A., Fortini, D., and Carattoli, A. (2010). Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants. *J. Antimicrob. Chemother.* 65, 2518–2529. doi: 10.1093/jac/dkq347

von Wintersdorff, C. J., Penders, J., Van Niekerk, J. M., Mills, N. D., Majumder, S., Van Alphen, L. B., et al. (2016). Dissemination of antimicrobial resistance in microbial ecosystems through horizontal gene transfer. *Front. Microbiol.* 7:173. doi: 10.3389/fmicb.2016.00173

Wickham, H. (2016). *ggplot2: Elegant Graphics for Data Analysis*. New York, NY: Springer-Verlag.

Zhong, G. (2017). Chlamydial plasmid-dependent pathogenicity. *Trends Microbiol.* 25, 141–152. doi: 10.1016/j.tim.2016.09.006

Zhou, F., and Xu, Y. (2010). cBar: a computer program to distinguish plasmid-derived from chromosome-derived sequence fragments in metagenomics data. *Bioinformatics* 26, 2051–2052. doi: 10.1093/bioinformatics/ btq299

**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.

Copyright © 2020 Douarre, Mallet, Radomski, Felten and Mistou. 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.