A global survey of Salmonella plasmids and their associations with antimicrobial resistance

Plasmids are the primary vector for horizontal transfer of antimicrobial resistance (AMR) within bacterial populations. We applied the MOB-suite, a toolset for reconstructing and typing plasmids, to 150 767 publicly available Salmonella whole-genome sequencing samples covering 1204 distinct serovars to produce a large-scale population survey of plasmids based on the MOB-suite plasmid nomenclature. Reconstruction yielded 183 017 plasmids representing 1044 primary MOB-clusters and 830 potentially novel MOB-clusters. Replicon and relaxase typing were able to type 83.4 and 58 % of plasmids, respectively, compared to 99.9 % for MOB-clusters. Within this work, we developed an approach to characterize the horizonal transfer of MOB-clusters and AMR genes across different serotypes, as well as the diversity of MOB-cluster associations with AMR genes. Aggregating conjugative mobility predictions provided by the MOB-suite and their corresponding serovar entropy demonstrated that non-mobilizable plasmids were associated with fewer serotypes compared to mobilizable or conjugative MOB-clusters. The host-range predictions for MOB-clusters also showed differences between the mobility classes, with mobilizable MOB-clusters accounting for 88.3 % of the multi-phyla (broad-host-range) predictions compared to 3 and 8.6 % for conjugative and non-mobilizable, respectively. A total of 296 (22 %) of identified MOB-clusters were associated with at least one resistance gene, indicating that the majority of Salmonella plasmids are not involved in AMR dissemination. Shannon entropy analysis of horizontal transfer of AMR genes across serovars and MOB-clusters demonstrated that horizonal transfer of genes is higher between serovars compared to transfer between different MOB-clusters. In addition to the population structure characterization based on primary MOB-clusters, we characterized a multi-plasmid outbreak responsible for the global dissemination of bla (CMY-2) across different serotypes using higher resolution MOB-suite secondary cluster codes. The plasmid characterization approach developed here can be applied to different organisms to identify plasmids and genes which pose high risks for horizontal transfer.