Research article

Complete genome sequence of Shigella flexneri 5b and comparison with Shigella flexneri 2a

Huan Nie^1,2* , Fan Yang^2* , Xiaobing Zhang^2* , Jian Yang^2* , Lihong Chen² , Jing Wang² , Zhaohui Xiong² , Junping Peng² , Lilian Sun² , Jie Dong² , Ying Xue² , Xingye Xu² , Shuxia Chen² , Zhijian Yao³ , Yan Shen³ and Qi Jin^1,2,4

¹  College of Biological Sciences China Agricultural University, Beijing 100094, China

²  State Key Laboratory for Molecular Virology and Genetic Engineering, Beijing 100052, China

³  National Center of Human Genome Research, Beijing 100176, China

⁴  Institute of Pathogen Biology, Chinese Academy of Medical Sciences, Beijing 100730, China

author email corresponding author email* Contributed equally

BMC Genomics 2006, 7:173doi:10.1186/1471-2164-7-173

Published:	6 July 2006

Abstract

Background

Shigella bacteria cause dysentery, which remains a significant threat to public health. Shigella flexneri is the most common species in both developing and developed countries. Five Shigella genomes have been sequenced, revealing dynamic and diverse features. To investigate the intra-species diversity of S. flexneri genomes further, we have sequenced the complete genome of S. flexneri 5b strain 8401 (abbreviated Sf8401) and compared it with S. flexneri 2a (Sf301).

Results

The Sf8401 chromosome is 4.5-Mb in size, a little smaller than that of Sf301, mainly because the former lacks the SHI-1 pathogenicity island (PAI). Compared with Sf301, there are 6 inversions and one translocation in Sf8401, which are probably mediated by insertion sequences (IS). There are clear differences in the known PAIs between these two genomes. The bacteriophage SfV segment remaining in SHI-O of Sf8401 is clearly larger than the remnants of bacteriophage SfII in Sf301. SHI-1 is absent from Sf8401 but a specific related protein is found next to the pheV locus. SHI-2 is involved in one intra-replichore inversion near the origin of replication, which may change the expression of iut/iuc genes. Moreover, genes related to the glycine-betaine biosynthesis pathway are present only in Sf8401 among the known Shigella genomes.

Conclusion

Our data show that the two S. flexneri genomes are very similar, which suggests a high level of structural and functional conservation between the two serotypes. The differences reflect different selection pressures during evolution. The ancestor of S. flexneri probably acquired SHI-1 and SHI-2 before SHI-O was integrated and the serotypes diverged. SHI-1 was subsequently deleted from the S. flexneri 5b genome by recombination, but stabilized in the S. flexneri 2a genome. These events may have contributed to the differences in pathogenicity and epidemicity between the two serotypes of S. flexneri.