plsMD: A plasmid reconstruction tool from short-read assemblies

Abstract

While whole genome sequencing (WGS) has become a cornerstone of antimicrobial resistance (AMR) surveillance, the reconstruction of plasmid sequences from short-read WGS data remains a challenge due to repetitive sequences and assembly fragmentation. Current computational tools for plasmid identification and binning, such as PlasmidFinder, cBAR, PlasmidSPAdes, and Mob-recon, have limitations in reconstructing full plasmid sequences, hindering downstream analyses like phylogenetic studies and AMR gene tracking. To address this gap, we present plsMD, a tool designed for full plasmid reconstruction from short-read assemblies. plsMD integrates Unicycler assemblies with replicon and full plasmid sequence databases (PlasmidFinder, MOB-typer and PLSDB) to guide plasmid reconstruction through a series of contig manipulations. Using two datasets, one established benchmark dataset used in previous benchmarking studies and another novel dataset consisting of newly sequenced bacterial isolates, plsMD outperformed existing tools in both. In the benchmark dataset, it achieved excellent recall, precision, and F1 scores of 91.3%, 95.5%, and 92.0%, respectively. In the novel dataset, it achieved good recall, precision, and F1 scores of 77.6%, 88.9%, and 74.5%, respectively. plsMD supports two usage modalities: single-sample analysis for plasmid reconstruction and gene annotation, and batch-sample analysis for phylogenetic investigations of plasmid transmission. This computational tool represents a significant advancement in plasmid analysis, offering a robust solution for utilizing existing short-read WGS data to study plasmid-mediated AMR spread and evolution. Key pointsO_LIAccurate plasmid reconstruction from short-read assemblies, surpassing existing binning-based tools. C_LIO_LIReplicon-guided approach enables detection of divergent plasmids. C_LIO_LISupports single and batch sample analysis, enabling gene annotation as well as plasmid transmission and evolutionary studies. C_LI Bibliographical noteThe authors are members of the Genomics and Metagenomics Program at Childrens Cancer Hospital Egypt, working in cancer genomics, bioinformatics, and antimicrobial resistance research. Graphical abstractplsMD uses Unicycler assemblies to identify plasmid replicons via PlasmidFinder and MOB-typer, and aligns the assemblies to PLSDB. It refines alignments, selects reference plasmids, and reconstructs full plasmid sequences. Circular plasmids without replicons are identified separately. plsMD outputs plasmid and non-plasmid FASTA files. Two workflows are supported: single-sample (plasmid/non-plasmid separation, annotation of AMR, VF, IS, and replicons) and batch-sample (grouping plasmids by replicon, MAFFT alignment, rotation, and phylogenetic tree construction). Validation on two datasets showed notably better performance than other benchmarking tools. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=150 SRC="FIGDIR/small/643493v2_ufig1.gif" ALT="Figure 1"> View larger version (31K): org.highwire.dtl.DTLVardef@719537org.highwire.dtl.DTLVardef@fb761dorg.highwire.dtl.DTLVardef@1076b72org.highwire.dtl.DTLVardef@1f2fd9_HPS_FORMAT_FIGEXP M_FIG C_FIG