Inversions over the terminus region in Salmonella and Escherichia coli: IS200s as the sites of homologous recombination inverting the chromosome of Salmonella enterica serovar typhi

Genomic rearrangements (duplications and inversions) in enteric bacteria such as Salmonella enterica serovar Typhimurium LT2 and Escherichia coli K12 are frequent (10(-3) to 10(-5)) in culture, but in wild-type strains these genomic rearrangements seldom survive. However, inversions commonly survive in the terminus of replication (TER) region, where bidirectional DNA replication terminates; nucleotide sequences from S. enterica serovar Typhimurium LT2, S. enterica serovar Typhi CT18, E. coli K12, and E. coli O157:H7 revealed genomic inversions spanning the TER region. Assuming that S. enterica serovar Typhimurium LT2 represents the ancestral genome structure, we found an inversion of 556 kb in serovar Typhi CT18 between two of the 25 IS200 elements and an inversion of about 700 kb in E. coli K12 and E. coli O157:H7. In addition, there is another inversion of 500 kb in E. coli O157:H7 compared with E. coli K12. PCR analysis confirmed that all S. enterica serovar Typhi strains tested, but not strains of other Salmonella serovars, have an inversion at the exact site of the IS200 insertions. We conclude that inversions of the TER region survive because they do not significantly change replication balance or because they are part of the compensating mechanisms to regain chromosome balance after it is disrupted by insertions, deletions, or other inversions.     

Characterization of a T7-like lytic bacteriophage (phiSG-JL2) of Salmonella enterica serovar gallinarum biovar gallinarum

PhiSG-JL2 is a newly discovered lytic bacteriophage infecting Salmonella enterica serovar Gallinarum biovar Gallinarum but is nonlytic to a rough vaccine strain of serovar Gallinarum biovar Gallinarum (SG-9R), S. enterica serovar Enteritidis, S. enterica serovar Typhimurium, and S. enterica serovar Gallinarum biovar Pullorum. The phiSG-JL2 genome is 38,815 bp in length (GC content, 50.9%; 230-bp-long direct terminal repeats), and 55 putative genes may be transcribed from the same strand. Functions were assigned to 30 genes based on high amino acid similarity to known proteins. Most of the expected proteins except tail fiber (31.9%) and the overall organization of the genomes were similar to those of yersiniophage phiYeO3-12. phiSG-JL2 could be classified as a new T7-like virus and represents the first serovar Gallinarum biovar Gallinarum phage genome to be sequenced. On the basis of intraspecific ratios of nonsynonymous to synonymous nucleotide changes (Pi[a]/Pi[s]), gene 2 encoding the host RNA polymerase inhibitor displayed Darwinian positive selection. Pretreatment of chickens with phiSG-JL2 before intratracheal challenge with wild-type serovar Gallinarum biovar Gallinarum protected most birds from fowl typhoid. Therefore, phiSG-JL2 may be useful for the differentiation of serovar Gallinarum biovar Gallinarum from other Salmonella serotypes, prophylactic application in fowl typhoid control, and understanding of the vertical evolution of T7-like viruses.     

The evolving genome of Salmonella enterica serovar Pullorum

Salmonella enterica serovar Pullorum is a fowl-adapted bacterial pathogen that causes dysentery (pullorum disease). Host adaptation and special pathogenesis make S. enterica serovar Pullorum an exceptionally good system for studies of bacterial evolution and speciation, especially regarding pathogen-host interactions and the acquisition of pathogenicity. We constructed a genome map of S. enterica serovar Pullorum RKS5078, using I-CeuI, XbaI, AvrII, and SpeI and Tn10 insertions. Pulsed-field gel electrophoresis was employed to separate the large DNA fragments generated by the endonucleases. The genome is 4,930 kb, which is similar to most salmonellas. However, the genome of S. enterica serovar Pullorum RKS5078 is organized very differently from the majority of salmonellas, with three major inversions and one translocation. This extraordinary genome structure was seen in most S. enterica serovar Pullorum strains examined, with different structures in a minority of S. enterica serovar Pullorum strains. We describe the coexistence of different genome structures among the same bacteria as genomic plasticity. Through comparisons with S. enterica serovar Typhimurium, we resolved seven putative insertions and eight deletions ranging in size from 12 to 157 kb. The genomic plasticity seen among S. enterica serovar Pullorum strains supported our hypothesis about its association with bacterial evolution: a large genomic insertion (157 kb in this case) disrupted the genomic balance, and rebalancing by independent recombination events in individual lineages resulted in diverse genome structures. As far as the structural plasticity exists, the S. enterica serovar Pullorum genome will continue evolving to reach a further streamlined and balanced structure.     

Diversity of genome structure in Salmonella enterica serovar Typhi populations

The genomes of most strains of Salmonella and Escherichia coli are highly conserved. In contrast, all 136 wild-type strains of Salmonella enterica serovar Typhi analyzed by partial digestion with I-CeuI (an endonuclease which cuts within the rrn operons) and pulsed-field gel electrophoresis and by PCR have rearrangements due to homologous recombination between the rrn operons leading to inversions and translocations. Recombination between rrn operons in culture is known to be equally frequent in S. enterica serovar Typhi and S. enterica serovar Typhimurium; thus, the recombinants in S. enterica serovar Typhi, but not those in S. enterica serovar Typhimurium, are able to survive in nature. However, even in S. enterica serovar Typhi the need for genome balance and the need for gene dosage impose limits on rearrangements. Of 100 strains of genome types 1 to 6, 72 were only 25.5 kb off genome balance (the relative lengths of the replichores during bidirectional replication from oriC to the termination of replication [Ter]), while 28 strains were less balanced (41 kb off balance), indicating that the survival of the best-balanced strains was greater. In addition, the need for appropriate gene dosage apparently selected against rearrangements which moved genes from their accustomed distance from oriC. Although rearrangements involving the seven rrn operons are very common in S. enterica serovar Typhi, other duplicated regions, such as the 25 IS200 elements, are very rarely involved in rearrangements. Large deletions and insertions in the genome are uncommon, except for deletions of Salmonella pathogenicity island 7 (usually 134 kb) from fragment I-CeuI-G and 40-kb insertions, possibly a prophage, in fragment I-CeuI-E. The phage types were determined, and the origins of the phage types appeared to be independent of the origins of the genome types.     

A genome map of Salmonella enterica serovar Agona: numerous insertions and deletions reflecting the evolutionary history of a human pathogen

Salmonella enterica serovar Agona is an important zoonotic pathogen, causing serious human illness worldwide, but knowledge about its genetics and evolution, especially regarding the genomic events that might have contributed to the formation of S . Agona as an important pathogen, is lacking. As a first step toward understanding this pathogen and characterizing its genomic differences with other salmonellae, we constructed a physical map of S . Agona in strain SARB1 using I-CeuI, XbaI, AvrII and Tn 10 insertions with pulsed-field gel electrophoresis techniques. On the 4815-kb genomic map, we located 82 genes, revealed one inversion of about 1000 kb and resolved seven deletions and seven insertions ranging from 10 to 67 kb relative to the genome of Salmonella typhimurium LT2. These genomic features clearly distinguish S . Agona from other previously analyzed salmonellae and provide clues to the molecular basis for its genomic divergence. Additionally, these kinds of physical maps, combined with emerging high-speed sequencing technologies, such as the Solexa or SOLiD techniques, which require a pre-existing high-resolution physical map such as the S . Agona map reported here, will play important roles in genomic comparative studies of bacteria involving large numbers of strains.     

Genomic diversification among archival strains of Salmonella enterica serovar typhimurium LT7

To document genomic changes during long periods of storage, we analyzed Salmonella enterica serovar Typhimurium LT7, a mutator strain that was previously reported to have higher rates of mutations compared to other serovar Typhimurium strains such as LT2. Upon plating directly from sealed agar stabs that had been stocked at room temperature for up to four decades, many auxotrophic mutants derived from LT7 gave rise to colonies of different sizes. Restreaking from single colonies consistently yielded colonies of diverse sizes even when we repeated single-colony isolation nine times. Colonies from the first plating had diverse genomic changes among and even within individual vials, including translocations, inversions, duplications, and point mutations, which were detected by rare-cutting endonuclease analysis with pulsed-field gel electrophoresis. Interestingly, even though the colony size kept diversifying, all descendents of the same single colonies from the first plating had the same sets of detected genomic changes. We did not detect any colony size or genome structure diversification in serovar Typhimurium LT7 stocked at -70 degrees C or in serovar Typhimurium LT2 stocked either at -70 degrees C or at room temperature. These results suggest that, although colony size diversification occurred during rapid growth, all detected genomic changes took place during the storage at room temperature and were carried over to their descendents without further changes during rapid growth in rich medium. We constructed a genomic cleavage map on the LT7 strain that had been stocked at -70 degrees C and located all of the detected genomic changes on the map. We speculated on the significance of mutators for survival and evolution under environmentally stressed conditions.     

Differences in gene content between Salmonella enterica serovar Enteritidis isolates and comparison to closely related serovars Gallinarum and Dublin

Salmonella enterica serovar Enteritidis is often transmitted into the human food supply through eggs of hens that appear healthy. This pathogen became far more prevalent in poultry following eradication of the fowl pathogen S. enterica serovar Gallinarum in the mid-20th century. To investigate whether changes in serovar Enteritidis gene content contributed to this increased prevalence, and to evaluate genetic heterogeneity within the serovar, comparative genomic hybridization was performed on eight 60-year-old and nineteen 10- to 20-year-old serovar Enteritidis strains from various hosts, using a Salmonella-specific microarray. Overall, almost all the serovar Enteritidis genomes were very similar to each other. Excluding two rare strains classified as serovar Enteritidis in the Salmonella reference collection B, only eleven regions of the serovar Enteritidis phage type 4 (PT4) chromosome (sequenced at the Sanger Center) were absent or divergent in any of the other serovar Enteritidis strains tested. The more recent isolates did not have consistent differences from 60-year-old field isolates, suggesting that no large genomic additions on a whole-gene scale were needed for serovar Enteritidis to become more prevalent in domestic fowl. Cross-hybridization of phage genes on the array with related genes in the examined genomes grouped the serovar Enteritidis isolates into two major lineages. Microarray comparisons of the sequenced serovar Enteritidis PT4 to isolates of the closely related serovars Dublin and Gallinarum (biovars Gallinarum and Pullorum) revealed several genomic areas that distinguished them from serovar Enteritidis and from each other. These differences in gene content could be useful in DNA-based typing and in understanding the different phenotypes of these related serovars.     

DNA amplification and rearrangements in archival Salmonella enterica serovar Typhimurium LT2 cultures

Variations in genome size and gene order were observed in archival Salmonella enterica serovar Typhimurium cultures stored for over 40 years. In one strain, microarray analysis revealed a large, stable amplification. PCR analysis of the same strain revealed a genomic duplication that underwent a translocation. Other strains had smaller duplications and deletions. These results demonstrate that storage in stabs over time at room temperature not only allows for further bacterial growth but also may produce an environment that selects for a variety of mutations, including genomic rearrangements.     

Variation between pathogenic serovars within Salmonella pathogenicity islands

Although four of the five Salmonella pathogenicity islands (SPIs) have been characterized in detail for Salmonella enterica serovar Typhimurium, and the fifth has been characterized for Salmonella enterica serovar Dublin, there have been limited studies to examine them in detail in a range of pathogenic serovars of S. enterica. The aim of this study was to examine these regions, shown to be crucial in virulence, in pathogenic serovars to identify any major deletions or insertions that may explain variation in virulence and provide further understanding of the elements involved in the evolution of these regions. Multiple strains of each of the 13 serovars were compared by Southern blot hybridization using a series of probes that together encompassed the full length of all five SPIs. With the exception of serovar Typhimurium, all strains of the same serovar were identical in all five SPIs. Those serovars that differed from serovar Typhimurium in SPI-1 to SPI-4 and from serovar Dublin in SPI-5 were examined in more detail in the variant regions by PCR, and restriction endonuclease digestion and/or DNA sequencing. While most variation in hybridization patterns was attributable to loss or gain of single restriction endonuclease cleavage sites, three regions, in SPI-1, SPI-3, and SPI-5, had differences due to major insertions or deletions. In SPI-1 the avrA gene was replaced by a 200-base fragment in three serovars, as reported previously. In SPI-5, two serovars had acquired an insertion with similarity to the pagJ and pagK genes between pipC and pipD. In SPI-3 the genes sugR and rhuM were deleted in most serovars and in some were replaced by sequences that were very similar to either the Escherichia coli fimbrial operon, flanked by two distinct insertion sequence elements, or to the E. coli retron phage PhiR73. The distribution of these differences suggests that there have been a number of relatively recent horizontal transfers of genes into S. enterica and that in some cases the same event has occurred in multiple lineages of S. enterica. Thus, it seems that insertion sequences and retron phages are likely to be involved in continuing evolution of the pathogenicity islands of pathogenic Salmonella serovars.     

Transposition of the heat-stable toxin astA gene into a gifsy-2-related prophage of Salmonella enterica serovar Abortusovis

The horizontal transfer and acquisition of virulence genes via mobile genetic elements have been a major driving force in the evolution of Salmonella pathogenicity. Serovars of Salmonella enterica carry variable assortments of phage-encoded virulence genes, suggesting that temperate phages play a pivotal role in this process. Epidemic isolates of S. enterica serovar Typhimurium are consistently lysogenic for two lambdoid phages, Gifsy-1 and Gifsy-2, carrying known virulence genes. Other serovars of S. enterica, including serovars Dublin, Gallinarum, Enteritidis, and Hadar, carry distinct prophages with similarity to the Gifsy phages. In this study, we analyzed Gifsy-related loci from S. enterica serovar Abortusovis, a pathogen associated exclusively with ovine infection. A cryptic prophage, closely related to serovar Typhimurium phage Gifsy-2, was identified. This element, named Gifsy-2AO, was shown to contribute to serovar Abortusovis systemic infection in lambs. Sequence analysis of the prophage b region showed a large deletion which covers genes encoding phage tail fiber proteins and putative virulence factors, including type III secreted effector protein SseI (GtgB, SrfH). This deletion was identified in most of the serovar Abortusovis isolates tested and might be dependent on the replicative transposition of an adjacent insertion sequence, IS1414, previously identified in pathogenic Escherichia coli strains. IS1414 encodes heat-stable toxin EAST1 (astA) and showed multiple genomic copies in isolates of serovar Abortusovis. To our knowledge, this is the first evidence of intergeneric transfer of virulence genes via insertion sequence elements in Salmonella. The acquisition of IS1414 (EAST1) and its frequent transposition within the chromosome might improve the fitness of serovar Abortusovis within its narrow ecological niche.     

Detection of a Salmonella enterica serovar California strain spreading in spanish feed mills and genetic characterization with DNA microarrays

We performed an epidemiological study on Salmonella isolated from raw plant-based feed in Spanish mills. Overall, 32 different Salmonella serovars were detected. Despite its rare occurrence in humans and animals, Salmonella enterica serovar California was found to be the predominant serovar in Spanish feed mills. Different typing techniques showed that isolates of this serovar were genetically closely related, and comparative genomic hybridization using microarray technology revealed 23 S. enterica serovar Typhimurium LT2 gene clusters that are absent from serovar California.     

Genome-wide detection of spontaneous chromosomal rearrangements in bacteria

Genome rearrangements have important effects on bacterial phenotypes and influence the evolution of bacterial genomes. Conventional strategies for characterizing rearrangements in bacterial genomes rely on comparisons of sequenced genomes from related species. However, the spectra of spontaneous rearrangements in supposedly homogenous and clonal bacterial populations are still poorly characterized. Here we used 454 pyrosequencing technology and a 'split mapping' computational method to identify unique junction sequences caused by spontaneous genome rearrangements in chemostat cultures of Salmonella enterica Var. Typhimurium LT2. We confirmed 22 unique junction sequences with a junction microhomology more than 10 bp and this led to an estimation of 51 true junction sequences, of which 28, 12 and 11 were likely to be formed by deletion, duplication and inversion events, respectively. All experimentally confirmed rearrangements had short inverted (inversions) or direct (deletions and duplications) homologous repeat sequences at the endpoints. This study demonstrates the feasibility of genome wide characterization of spontaneous genome rearrangements in bacteria and the very high steady-state frequency (20-40%) of rearrangements in bacterial populations.     

Genome sequences of Salmonella enterica serovar typhimurium, Choleraesuis, Dublin, and Gallinarum strains of well- defined virulence in food-producing animals

Salmonella enterica is an animal and zoonotic pathogen of worldwide importance and may be classified into serovars differing in virulence and host range. We sequenced and annotated the genomes of serovar Typhimurium, Choleraesuis, Dublin, and Gallinarum strains of defined virulence in each of three food-producing animal hosts. This provides valuable measures of intraserovar diversity and opportunities to formally link genotypes to phenotypes in target animals.     

Direct and quantitative analysis of Salmonella enterica serovar Typhimurium using real-time PCR from artificially contaminated chicken meat

For quantitative PCR assay of Salmonella enterica serovar Typhimurium in food samples, a real-time PCR method was developed, based on DNA genome equivalent. Specific primers and probe designed based on the STM4497 gene of S. Typhimurium LT2 showed the specificity to S. Typhimurium. Threshold cycle (Ct) values of real-time PCR were obtained from a quantitative standard curve with genomic DNA of Salmonella Typhimurium. In addition, the recovery of S. Typhimurium inoculated artificially to chicken samples with 4.5x105 to 4.5x100 CFU/ml was evaluated by using real-time PCR and plate-count methods. Result showed that the number of cells calculated from the real-time PCR method had good correlation with that of the plate-count method. This real-time PCR method could be applicable to the detection and quantification of S. Typhimurium in food samples.     

Identification of specific gene sequences conserved in contemporary epidemic strains of Salmonella enterica

Genetic elements specific to recent and contemporary epidemic strains of Salmonella enterica were identified using comparative genomic analysis. Two epidemic multidrug-resistant (MDR) strains, MDR Salmonella enterica serovar Typhimurium definitive phage type 104 (DT104) and cephalosporin-resistant MDR Salmonella enterica serovar Newport, and an epidemic pansusceptible strain, Salmonella serovar Typhimurium DT160, were subjected to Salmonella gene microarray and suppression subtractive hybridization analyses. Their genome contents were compared with those of coexisting sporadic strains matched by serotype, geographic and temporal distribution, and host species origin. These paired comparisons revealed that epidemic strains of S. enterica had specific genes and gene regions that were shared by isolates of the same subtype. Most of these gene sequences are related to mobile genetic elements, including phages, plasmids, and plasmid-like and transposable elements, and some genes may encode proteins conferring growth or survival advantages. The emergence of epidemic MDR strains may therefore be associated with the presence of fitness-associated genetic factors in addition to their antimicrobial resistance genes.     

The SopEPhi phage integrates into the ssrA gene of Salmonella enterica serovar Typhimurium A36 and is closely related to the Fels-2 prophage

Salmonella spp. are enteropathogenic gram-negative bacteria that use a large array of virulence factors to colonize the host, manipulate host cells, and resist the host's defense mechanisms. Even closely related Salmonella strains have different repertoires of virulence factors. Bacteriophages contribute substantially to this diversity. There is increasing evidence that the reassortment of virulence factor repertoires by converting phages like the GIFSY phages and SopEPhi may represent an important mechanism in the adaptation of Salmonella spp. to specific hosts and to the emergence of new epidemic strains. Here, we have analyzed in more detail SopEPhi, a P2-like phage from Salmonella enterica serovar Typhimurium DT204 that encodes the virulence factor SopE. We have cloned and characterized the attachment site (att) of SopEPhi and found that its 47-bp core sequence overlaps the 3' terminus of the ssrA gene of serovar Typhimurium. Furthermore, we have demonstrated integration of SopEPhi into the cloned attB site of serovar Typhimurium A36. Sequence analysis of the plasmid-borne prophage revealed that SopEPhi is closely related to (60 to 100% identity over 80% of the genome) but clearly distinct from the Fels-2 prophage of serovar Typhimurium LT2 and from P2-like phages in the serovar Typhi CT18 genome. Our results demonstrate that there is considerable variation among the P2-like phages present in closely related Salmonella spp.     

Adaptation of multilocus sequencing for studying variation within a major clone: evolutionary relationships of Salmonella enterica serovar Typhimurium

Serovar Typhimurium of Salmonella enterica is a model organism for studies of pathogenesis that exhibits phage-type variation and variation in host range and virulence, but in a recent study showed no sequence variation in four genes, indicating the clonal nature of this serovar. We determined the relationships of 46 Typhimurium isolates of nine phage types using mutational changes detected either by matching AFLP (amplified fragment length polymorphism) fragments to computer-modeled LT2 AFLP fragments or by sequencing intergenic regions. Fifty-one polymorphic sites were detected, which gave a single phylogenetic tree. Comparison with genome sequences of five other serovars, Typhi, Paratyphi A, Gallinarum, Enteritidis, and Pullorum, enabled determination of the root of the tree. Only two parallel events were observed, giving high confidence in the tree branching order. The mutation-based tree provided a high level of consistency and a clear lineage for the Typhimurium isolates studied. This enabled us to show that for seven of the nine phage types used, the isolates studied have a single origin, but that two phage types clearly have more than one independent origin. We found that sequencing intergenic regions provides a good strategy for detection of mutational polymorphisms and study of phylogenetic relationships of closely related isolates and would be applicable to many other species.