Large metabolic rewiring from small genomic changes between strains of Shigella flexneri

The instability of Shigella genomes has been described, but how this instability causes phenotypic differences within the Shigella flexneri species is largely unknown and likely variable. We describe herein the genome of S. flexneri strain PE577, originally a clinical isolate, which exhibits several phenotypic differences compared to the model strain 2457T. Like many previously described strains of S. flexneri, PE577 lacks discernible, functional CRISPR and restriction-modification systems. Its phenotypic differences when compared to 2457T include lower transformation efficiency, higher oxygen sensitivity, altered carbon metabolism, and greater susceptibility to a wide variety of lytic bacteriophage isolates. Since relatively few Shigella phages have been isolated on 2457T or the previously characterized strain M90T, developing a more universal model strain for isolating and studying Shigella phages is critical to understanding both phages and phage-host interactions. In addition to phage biology, the genome sequence of PE577 was used to generate and test hypotheses of how pseudogenes in this strain-whether interrupted by degraded prophages, transposases, frameshifts, or point mutations-have led to metabolic rewiring compared to the model strain 2457T. Results indicate that PE577 can utilise the less-efficient pyruvate oxidase/acetyl-CoA synthetase (PoxB/Acs) pathway to produce acetyl-CoA, while strain 2457T cannot due to a nonsense mutation in acs, rendering it a pseudogene in this strain. Both strains also utilize pyruvate-formate lyase to oxidize formate but cannot survive with this pathway alone, possibly because a component of the formate-hydrogen lyase (fdhF) is a pseudogene in both strains.Importance Shigella causes millions of dysentery cases worldwide, primarily affecting children under five years old. Despite active research in developing vaccines and new antibiotics, relatively little is known about the variation of physiology or metabolism across multiple isolates. In this work, we investigate two strains of S. flexneri that share 98.9% genetic identity but exhibit drastic differences in metabolism, ultimately affecting the growth of the two strains. Results suggest additional strains within the S. flexneri species utilize different metabolic pathways to process pyruvate. Metabolic differences between these closely-related isolates suggest an even wider variety of differences in growth across S. flexneri and Shigella in general. Exploring this variation further may assist the development or application of vaccines and therapeutics to combat Shigella infections.