Salmonella enterica vs Escherichia coli Core Synteny and Pathogenicity Islands

Salmonella enterica and Escherichia coli are sister genera of facultatively anaerobic, peritrichously flagellated, Gram-negative enteric rods within the family Enterobacteriaceae (order Enterobacterales) that diverged from a common ancestor approximately 100 to 140 million years ago—roughly contemporaneous with the evolutionary radiation of placental mammals (Eutheria). Despite >100 million years of independent evolution and ecological divergence (E. coli as a lactose-fermenting commensal of the lower mammalian colon versus S. enterica as a lactose-negative, hydrogen sulfide H₂S-producing intracellular pathogen across reptiles, birds, and mammals), whole-genome alignments between S. enterica serovar Typhimurium LT2 (4.86 Mb, 52.2% GC) and E. coli K-12 MG1655 (4.64 Mb, 50.8% GC) demonstrate remarkable chromosomal collinearity (synteny) across ~3,700 shared orthologous core housekeeping genes (~80–85% average nucleotide identity ANI).

Rather than gradual point-mutational drift of their shared metabolic backbone, the divergence of Salmonella from Escherichia was driven by episodic horizontal gene transfer (HGT) of large, AT-rich genomic islands integrated into chromosomal tRNA genes: specifically, Salmonella Pathogenicity Island 1 (SPI-1, ~40 kb), which was acquired at the base of the genus Salmonella (S. enterica and S. bongori) to enable invasion of non-phagocytic intestinal epithelial cells via the T3SS-1 injectisome; and Salmonella Pathogenicity Island 2 (SPI-2, ~40 kb), acquired subsequently by S. enterica at the tRNA-Val (valV) locus to enable survival and replication inside macrophage phagolysosomes via the T3SS-2 nanomachine. In clinical metagenomics and bioinformatics, comparing S. enterica and E. coli also provides the textbook empirical demonstration of why multi-copy 16S ribosomal RNA (rrsA–rrsH) barcoding collapses between enteric sister genera (d_K2P = 0.0078 inter-species vs 0.0091 intra-genomic E. coli paralog divergence), requiring single-copy rpoB, gyrB, or whole-contig tetranucleotide Kullback–Leibler divergence (D_KL) for diagnostic resolution.

Conservation of the oriC–ter Replichore Architecture

Across the 4.64–4.95 Mb circular chromosomes of S. enterica and E. coli, bidirectional DNA replication initiates at oriC and DnaA (~3,925 kb in K-12) and proceeds down two symmetric replichores toward the terminus (terC) bounded by Tus-binding TerA–TerJ replication forks traps. Over 85% of the ~3,700 shared orthologous genes (thrLABC → dnaK → dnaA → gyrB → rpoB–rpoC → atpIBEFHAGDC) maintain identical gene order and replichore orientation (leading-strand vs lagging-strand transcriptional bias to prevent head-on collisions between replisomes and RNA polymerase), interrupted only by a well-characterized ~500 kb symmetrical chromosomal inversion across the ter replication terminus (between rrlG and pyrD) in S. enterica serovar Typhimurium.

At the nucleotide level, orthologous core genes between S. enterica and E. coli exhibit highly constrained 2nd-codon-position nonsynonymous conservation (dN ≈ 0.02–0.05) under purifying selection, paired with substantial 3rd-codon-position synonymous divergence (dS ≈ 0.45–0.65, overall ω = dN/dS ≈ 0.04–0.08) accumulated over ~100–140 Myr according to the bacterial molecular clock.

Compositional Signature of Horizontally Acquired Pathogenicity Islands (SPI-1 to SPI-5)

Superimposed onto this conserved 51–52% GC housekeeping backbone are >15 lineage-specific Salmonella Pathogenicity Islands (SPI-1 through SPI-19) in S. enterica (mirroring the enteropathogenic LEE and uropathogenic PAI islands in pathogenic E. coli). Most SPIs exhibit four canonical signatures of horizontal acquisition:

1. Anomalous GC Content and Tetranucleotide Skew: Whereas the core S. enterica chromosome averages 52.2% GC, both SPI-1 (44.6% GC) and SPI-2 (44.7% GC) are strongly AT-rich, creating sharp local peaks (Z > +4.0 SD) in sliding-window tetranucleotide Kullback–Leibler divergence (D_KL) against the chromosomal background.
2. Site-Specific Integration at 3′ Ends of tRNA Genes: With the exception of SPI-1 (inserted between the housekeeping fhlA and mutS genes at centisome 63), SPI-2 (tRNA-Val valV at centisome 30.7), SPI-3 (tRNA-SelC selC at centisome 82), SPI-4 (tRNA-Leu leuX), and SPI-5 (tRNA-Ser serT) are integrated directly into the conserved 3′ CCA-proximal stem of chromosomal tRNA genes via bacteriophage-like tyrosine integrases.
3. Xenogeneic Silencing by H-NS and Counter-Silencing by HilD/SsrB: Because foreign AT-rich DNA (<47% GC) contains narrow minor-groove A-tracts, the nucleoid-associated histone-like protein H-NS Polymerizes across SPI-1 and SPI-2 to transcriptionally silence them under non-host conditions, preventing fitness cost until host environmental cues activate the dedicated AraC-like/two-component counter-silencers HilD (SPI-1) or SsrA–SsrB (SPI-2).

SPI-1 (Inv/Spa Injectisome): Luminal Sensing and Actin Membrane Ruffling

Salmonella Pathogenicity Island 1 (SPI-1, ~40 kb, encoding 39 genes) is present in both species of the genus (S. enterica and the cold-blooded reptile commensal S. bongori) and encodes the Type III Secretion System 1 (T3SS-1) molecular syringe (Inv/Spa/Prg structural needle complex) alongside its secreted translocons (SipB, SipC, SipD), chaperones (SicA, InvB), and effectors (SopB, SopE, SopE2, SipA, SptP).

When Salmonella reaches the distal ileum of the small intestine, environmental cues characteristic of the intestinal lumen—hyperosmolarity (~300 mOsm), low oxygen tension (microaerobiosis), neutral-to-alkaline pH (pH ~7.4), and short-chain fatty acids (butyrate/acetate gradient)—activate a feed-forward transcriptional cascade of three AraC-family regulators (HilD → HilC/RtsA → HilA), which relieves H-NS silencing and induces T3SS-1 assembly. Upon contacting microfold (M) cells or enterocytes, T3SS-1 injects:

• SopE and SopE2: Bacterial guanine nucleotide exchange factors (GEFs) that directly activate host Rho-family GTPases (Rac1 and Cdc42), coupled with the inositol phosphatase SopB (SigD) and actin-bundling protein SipA, triggering Arp2/3-dependent actin cytoskeletal rearrangements and apical membrane ruffling (macropinocytosis) that engulf the bacterium into a membrane-bound Salmonella-Containing Vacuole (SCV).
• SptP: A bifunctional bacterial GAP (GTPase-activating protein) and tyrosine phosphatase with a longer intracellular half-life than SopE that subsequently deactivates Rac1/Cdc42, restoring normal host enterocyte apical microvilli architecture after bacterial entry.

SPI-2 (Ssa/Sse Injectisome): Phagolysosomal Acidification Sensing and Systemic Dissemination

Whereas S. bongori lacks SPI-2 and remains restricted to the intestinal mucosa of reptiles, the evolutionary acquisition of the 40 kb SPI-2 island at the valV (tRNA-Val) locus enabled Salmonella enterica to survive inside vertebrate macrophages and cause systemic enteric fever (Typhoid fever in S. enterica serovar Typhi) and bacteremia.

Once S. enterica is phagocytosed by a lamina propria macrophage or resides inside an epithelial SCV, vacuolar ATPases (V-ATPase) acidify the lumen (pH 4.5–5.5) while depleting luminal Mg2+ (< 10–50 μM), Ca2+, and inorganic phosphate (PO4^3−). These intravacuolar cues rapidly shut down SPI-1 and activate the two-component sensor kinase systems PhoP–PhoQ (sensing low Mg2+ and antimicrobial peptides) and EnvZ–OmpR (sensing acidic pH), which activate the SPI-2-encoded sensor kinase SsrA (SpiR) and response regulator SsrB to displace H-NS and transcribe the T3SS-2 (Ssa apparatus, SseB/C/D translocon, and >25 Sse/Sif/PipB effectors). Projected across the SCV membrane into the host cytosol, T3SS-2 effectors (SifA, SseF, SseG, SpiC) recruit kinesin/dynein motors (SKIP) to extend tubular Salmonella-Induced Filaments (Sifs) along microtubules, tether the SCV to the trans-Golgi network for nutrient interception, and exclude NADPH phagocyte oxidase (NOX2) and inducible nitric oxide synthase (iNOS) vesicles from the SCV.

Intragenomic 16S Paralog Variance (rrnA–rrnH) Exceeds Inter-Species Divergence

A foundational finding of comparative microbial systematics is that short-read (V3–V4) and even full-length (1,542 nt) 16S ribosomal RNA sequencing fails to reliably separate Salmonella enterica from Escherichia coli (as well as Shigella and Citrobacter). Both E. coli and S. enterica carry seven rRNA operons (rrnA, rrnB, rrnC, rrnD, rrnE, rrnG, and rrnH) distributed around the chromosome (five on one replichore and two on the other, oriented codirectionally with DNA replication).

Because the seven 16S rRNA gene copies (rrsA–rrsH) within a single E. coli chromosome undergo incomplete concerted evolution (gene conversion) and harbor microheterogeneities concentrated in hypervariable stem-loops V1, V2, and V6 (helices h6 and h21, H_loop = 1.68 bits), the maximum intragenomic Kimura two-parameter (K2P) divergence among rrsA–rrsH within a single E. coli genome reaches d_K2P = 0.0091 (~14 nucleotide differences across 1,542 bp). Meanwhile, the minimum inter-genus 16S distance between E. coli and S. enterica is only d_K2P = 0.0078 (~99.2% sequence identity)—yielding a negative barcoding gap (gap = 0.0078 − 0.0091 = −0.0013) and splitting a single clonal E. coli genome into 4 to 6 distinct 100%-identity Amplicon Sequence Variants (ASVs, Hd = 0.895).

Diagnostic Rescue: Single-Copy rpoB/gyrB, invA PCR, and Whole-Contig D_KL

To resolve S. enterica from E. coli with 100% specificity in clinical microbiology and food-safety genomics, three orthogonal molecular assays replace unphased 16S amplicons:

1. Single-Copy Housekeeping Genes (rpoB and gyrB): Because rpoB (RNA polymerase β-subunit, 4,029 bp) and gyrB (DNA gyrase B-subunit, 2,415 bp) exist as strictly single-copy chromosomal genes (intragenomic variance = 0.0000) and accumulate synonymous 3rd-codon substitutions unrestricted by rRNA stem-loop secondary-structure base pairing, inter-species divergence jumps 9-fold to 14-fold (d_K2P = 0.068 for rpoB and d_K2P = 0.114 for gyrB), producing a wide positive barcoding gap with 100% bootstrap support.
2. Lineage-Specific SPI-1 invA** Virulence Target:** Diagnostic PCR and loop-mediated isothermal amplification (LAMP) assays target the invA gene (encoding the inner-membrane export apparatus protein of T3SS-1 inside SPI-1), which is conserved across >99.5% of all >2,600 Salmonella serovars and absent from E. coli.
3. Long-Read Operon Phasing (rrnA–rrnH) and Tetranucleotide D_KL: In shotgun metagenomics, phasing full 5 kb 16S–ITS–23S operons with flanking chromosomal synteny or scanning contigs for the 7.4× tetranucleotide D_KL elevation of SPI-1/SPI-2 separates Salmonella from commensal Enterobacteriaceae.

Why is Shigella classified taxonomically as four species (S. dysenteriae, S. flexneri, S. boydii, S. sonnei) when phylogenetically it lies entirely inside Escherichia coli?

Whole-genome phylogenetics demonstrates that the four named Shigella species are not a separate genus or even a single monophyletic clade; rather, they are polyphyletic lineages of Escherichia coli (ANI > 97–98%, d_K2P < 0.015 at rpoB/gyrB) that independently acquired the large ~210 kb pINV virulence plasmid (encoding the Mxi–Spa Type III Secretion System and IpaBCD invasion antigens) while losing flagellar motility (flhDC) and lysine decarboxylase (cadA). The genus name Shigella is retained in clinical medicine purely to prevent diagnostic confusion between commensal E. coli and bacillary dysentery (shigellosis).

Why did Salmonella bongori acquire SPI-1 but never acquire SPI-2?

Salmonella bongori diverged from S. enterica after their common ancestor acquired SPI-1 (T3SS-1), which suffices for invading intestinal epithelial cells and colonizing the gastrointestinal tract of poikilothermic (cold-blooded) reptiles (lizards, snakes, tortoises). The subsequent horizontal acquisition of SPI-2 (T3SS-2) at the valV (tRNA-Val) locus occurred exclusively in the stem lineage of Salmonella enterica, conferring the ability to survive inside the hotter (37–42 °C), highly oxidative phagolysosomes of endothermic mammalian and avian macrophages.

How does Salmonella enterica use tetrathionate respiration (ttrABC) inside the inflamed gut to outcompete commensal E. coli?

When SPI-1 and SPI-2 effectors trigger acute host intestinal inflammation (neutrophil influx), host reactive oxygen species (ROS) oxidize luminal thiosulfate (S2O3^2−, produced by colonic epithelium detoxifying microbiota-derived H₂S) into tetrathionate (S4O6^2−). Whereas commensal E. coli and strict anaerobes (Bacteroides, Clostridia) cannot use tetrathionate, S. enterica carries the SPI-2-regulated ttrSR–ttrBCA tetrathionate reductase operon, enabling it to use S4O6^2− as a terminal anaerobic electron acceptor to respire ethanolamine and 1,2-propanediol and bloom to >50% of the luminal population.

1. McClelland, M., Sanderson, K. E., Spieth, J., et al. (2001). "Complete genome sequence of Salmonella enterica serovar Typhimurium LT2." Nature, 413(6858), 852–856. https://doi.org/10.1038/35101614
2. Ochman, H., & Wilson, A. C. (1987). "Evolution in bacteria: evidence for a universal substitution rate in cellular genomes." Journal of Molecular Evolution, 26(1–2), 74–86. https://doi.org/10.1007/BF02111283
3. Hensel, M., Shea, J. E., Gleeson, C., Jones, M. D., Dalton, E., & Holden, D. W. (1995). "Simultaneous identification of bacterial virulence genes by negative selection." Science, 269(5222), 400–403. https://doi.org/10.1126/science.7618105
4. Winter, S. E., Thiennimitr, P., Winter, M. G., et al. (2010). "Gut inflammation provides a respiratory electron acceptor for Salmonella." Nature, 467(7314), 426–429. https://doi.org/10.1038/nature09415