Comparison of the virulence plasmid genomes of two strains of Shigella which lost the ability to bind Congo red

We determined and analyzed the Shigella flexneri serotype 5 (pSF5) and S. dysenteriae serotype 1 (pSD1) virulence plasmid genomes. The total length of pSF5 is 136513 bp, including 165 open reading frames (ORFs). Of these ORFs, 133 were identified and 32 of those had no significant homology to proteins with known functions. The length of pSD1 is 182545 bp, including 224 ORFs, of which we identified 181. The remaining 43 ORFs were not significantly homologous to proteins with known functions. The insertion sequence (IS) elements are 53787 bp in pSF5, and 49616 bp in pSD1, which represents 39.4% and 27.1% of the genome, respectively. There are 22 IS element types in pSF5 and pSD1, among which we report ISEc8 and ISSbo6 for the first time in the Shigella virulence plasmid. Compared to pCP301, there are a large number of deleted genes and gene inversions in both pSF5 and pSD1. The ipa-mxi-spa locus in pSF5 is completely absent, and the genes related to the O-antigen biosynthesis are partially missing. In contrast, the above genes in pSD1 are integral, with the exception of virF. The whole genome analysis of the two plasmids shows that the loss of genes related to gene invasion or regulation also obliterates the ability of pPF5 and pSD1 to bind Congo red (Crb). Whether these genes determine the Crb function requires continued investigation. These authors contributed equally to this work.