Characterization of the ELPhiS prophage from Salmonella enterica serovar Enteritidis strain LK5

Phages are a primary driving force behind the evolution of bacterial pathogens by transferring a variety of virulence genes into their hosts. Similar to other bacterial genomes, the Salmonella enterica serovar Enteritidis LK5 genome contains several regions that are homologous to phages. Although genomic analysis demonstrated the presence of prophages, it was unable to confirm which phage elements within the genome were viable. Genetic markers were used to tag one of the prophages in the genome to allow monitoring of phage induction. Commonly used laboratory strains of Salmonella were resistant to phage infection, and therefore a rapid screen was developed to identify susceptible hosts. This approach showed that a genetically tagged prophage, ELPhiS (Enteritidis lysogenic phage S), was capable of infecting Salmonella serovars that are diverse in host range and virulence and has the potential to laterally transfer genes between these serovars via lysogenic conversion. The rapid screen approach is adaptable to any system with a large collection of isolates and may be used to test the viability of prophages found by sequencing the genomes of various bacterial pathogens.     

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.     

Role of Genomic Rearrangements in Producing New Ribotypes of Salmonella typhi

Salmonella typhi is the only species of Salmonella which grows exclusively in humans, in whom it causes enteric typhoid fever. Strains of S. typhi show very little variation in electrophoretic types, restriction fragment length polymorphisms, cell envelope proteins, and intervening sequences, but the same strains are very heterogeneous for ribotypes which are detected with the restriction endonuclease PstI. In addition, the genome of S. typhi has been proven to undergo genomic rearrangement due to homologous recombination between the seven copies of rrn genes. The relationship between ribotype heterogeneity and genomic rearrangement was investigated. Strains of S. typhi which belong to 23 different genome types were analyzed by ribotyping. A limited number of ribotypes were found within the same genome type group; e. g., most strains of genome type 3 belonged to only two different ribotypes, which result from recombination between rrnH and rrnG operons. Different genome type groups normally have different ribotypes. The size and identity of the PstI fragment containing each of the seven different rrn operons from S. typhi Ty2 were determined, and from these data, one can infer how genomic rearrangement forms new ribotypes. It is postulated that genomic rearrangement, rather than mutation, is largely responsible for producing the ribotype heterogeneity in S. typhi.     

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.     

Evolution of paralogous genes: Reconstruction of genome rearrangements through comparison of multiple genomes within Staphylococcus aureus

Analysis of evolution of paralogous genes in a genome is central to our understanding of genome evolution. Comparison of closely related bacterial genomes, which has provided clues as to how genome sequences evolve under natural conditions, would help in such an analysis. With species Staphylococcus aureus, whole-genome sequences have been decoded for seven strains. We compared their DNA sequences to detect large genome polymorphisms and to deduce mechanisms of genome rearrangements that have formed each of them. We first compared strains N315 and Mu50, which make one of the most closely related strain pairs, at the single-nucleotide resolution to catalogue all the middle-sized (more than 10 bp) to large genome polymorphisms such as indels and substitutions. These polymorphisms include two paralogous gene sets, one in a tandem paralogue gene cluster for toxins in a genomic island and the other in a ribosomal RNA operon. We also focused on two other tandem paralogue gene clusters and type I restriction-modification (RM) genes on the genomic islands. Then we reconstructed rearrangement events responsible for these polymorphisms, in the paralogous genes and the others, with reference to the other five genomes. For the tandem paralogue gene clusters, we were able to infer sequences for homologous recombination generating the change in the repeat number. These sequences were conserved among the repeated paralogous units likely because of their functional importance. The sequence specificity (S) subunit of type I RM systems showed recombination, likely at the homology of a conserved region, between the two variable regions for sequence specificity. We also noticed novel alleles in the ribosomal RNA operons and suggested a role for illegitimate recombination in their formation. These results revealed importance of recombination involving long conserved sequence in the evolution of paralogous genes in the genome.     

Comparative genome analysis of the high pathogenicity Salmonella Typhimurium strain UK-1

Salmonella enterica serovar Typhimurium, a gram-negative facultative rod-shaped bacterium causing salmonellosis and foodborne disease, is one of the most common isolated Salmonella serovars in both developed and developing nations. Several S. Typhimurium genomes have been completed and many more genome-sequencing projects are underway. Comparative genome analysis of the multiple strains leads to a better understanding of the evolution of S. Typhimurium and its pathogenesis. S. Typhimurium strain UK-1 (belongs to phage type 1) is highly virulent when orally administered to mice and chickens and efficiently colonizes lymphoid tissues of these species. These characteristics make this strain a good choice for use in vaccine development. In fact, UK-1 has been used as the parent strain for a number of nonrecombinant and recombinant vaccine strains, including several commercial vaccines for poultry. In this study, we conducted a thorough comparative genome analysis of the UK-1 strain with other S. Typhimurium strains and examined the phenotypic impact of several genomic differences. Whole genomic comparison highlights an extremely close relationship between the UK-1 strain and other S. Typhimurium strains; however, many interesting genetic and genomic variations specific to UK-1 were explored. In particular, the deletion of a UK-1-specific gene that is highly similar to the gene encoding the T3SS effector protein NleC exhibited a significant decrease in oral virulence in BALB/c mice. The complete genetic complements in UK-1, especially those elements that contribute to virulence or aid in determining the diversity within bacterial species, provide key information in evaluating the functional characterization of important genetic determinants and for development of vaccines.     

Identification and characterization of novel salmonella mobile elements involved in the dissemination of genes linked to virulence and transmission

The genetic diversity represented by >2,500 different Salmonella serovars provides a yet largely uncharacterized reservoir of mobile elements that can contribute to the frequent emergence of new pathogenic strains of this important zoonotic pathogen. Currently, our understanding of Salmonella mobile elements is skewed by the fact that most studies have focused on highly virulent or common serovars. To gain a more global picture of mobile elements in Salmonella, we used prediction algorithms to screen for mobile elements in 16 sequenced Salmonella genomes representing serovars for which no prior genome scale mobile element data were available. From these results, selected mobile elements underwent further analyses in the form of validation studies, comparative analyses, and PCR-based population screens. Through this analysis we identified a novel plasmid that has two cointegrated replicons (IncI1-IncFIB); this plasmid type was found in four genomes representing different Salmonella serovars and contained a virulence gene array that had not been previously identified. A Salmonella Montevideo isolate contained an IncHI and an IncN2 plasmid, which both encoded antimicrobial resistance genes. We also identified two novel genomic islands (SGI2 and SGI3), and 42 prophages with mosaic architecture, seven of them harboring known virulence genes. Finally, we identified a novel integrative conjugative element (ICE) encoding a type IVb pilus operon in three non-typhoidal Salmonella serovars. Our analyses not only identified a considerable number of mobile elements that have not been previously reported in Salmonella, but also found evidence that these elements facilitate transfer of genes that were previously thought to be limited in their distribution among Salmonella serovars. The abundance of mobile elements encoding pathogenic properties may facilitate the emergence of strains with novel combinations of pathogenic traits.     

The relaxing ori-ter balance of Mycoplasma genomes

Mycoplasma are wall-less bacteria with small genomes, which are thought to have resulted from massive genome reductive processes, during which the ori-ter balance may be disrupted. For technical difficulties, ori and ter have been located only in a few Mycoplasma strains. Using the Z curve method, we were able to locate turning points on the Mycoplasma genomes, with the minimum and maximum points co-locating with ori or ter in the reference genomes. Assuming Z curve correctly located ori and ter, we calculated the distances from ori to ter in both directions on the circular genome and calculated the ori-ter balance status. The Mycoplasma genomes were not balanced, possibly as a result of close association of Mycoplasma with hosts, where there would be no other microbes for Mycoplasma to compete with for nutrients, so fastest possible growth related to balanced genomes might not be needed by Mycoplasma, leading to a relaxing ori-ter balance.     

Whole-genome sequencing of staphylococcus haemolyticus uncovers the extreme plasticity of its genome and the evolution of human-colonizing staphylococcal species

Staphylococcus haemolyticus is an opportunistic bacterial pathogen that colonizes human skin and is remarkable for its highly antibiotic-resistant phenotype. We determined the complete genome sequence of S.haemolyticus to better understand its pathogenicity and evolutionary relatedness to the other staphylococcal species. A large proportion of the open reading frames in the genomes of S.haemolyticus, Staphylococcus aureus, and Staphylococcus epidermidis were conserved in their sequence and order on the chromosome. We identified a region of the bacterial chromosome just downstream of the origin of replication that showed little homology among the species but was conserved among strains within a species. This novel region, designated the "oriC environ," likely contributes to the evolution and differentiation of the staphylococcal species, since it was enriched for species-specific nonessential genes that contribute to the biological features of each staphylococcal species. A comparative analysis of the genomes of S.haemolyticus, S.aureus, and S.epidermidis elucidated differences in their biological and genetic characteristics and pathogenic potentials. We identified as many as 82 insertion sequences in the S.haemolyticus chromosome that probably mediated frequent genomic rearrangements, resulting in phenotypic diversification of the strain. Such rearrangements could have brought genomic plasticity to this species and contributed to its acquisition of antibiotic resistance.     

Typing of Salmonella enterica subsp. enterica serovars Choleraesuis,Typhimurium, Dublin and laboratory strains of Escherichia coli using subtracted restriction fingerprinting (SRF)

The aim of the present study was to evaluate the usefulness of a new typing technique called subtracted restriction fingerprinting (SRF) for bacterial strain and isolate discrimination. The technique was applied to isolates of Salmonella enterica subsp. enterica (S.) serovars Choleraesuis, Typhimurium, Dublin and to two laboratory strains of E. coli. SRF is based on the selective removal of excess fragments from a restriction digest using magnetic particles. Subsequently, the remaining subset of restriction fragments can easily be analyzed with a conventional agarose gel. Larger fragments are preferentially removed by SRF. This results in an even distribution of bands within each electrophoretic lane and significantly improves scoring. The high discriminatory index for (S.) Choleraesuis (D = 0.914) illustrated the suitability of SRF for genome typing.     

Persistent bacterial infections: the interface of the pathogen and the host immune system

Persistent bacterial infections involving Mycobacterium tuberculosis, Salmonella enterica serovar Typhi (S. typhi) and Helicobacter pylori pose significant public-health problems. Multidrug-resistant strains of M. tuberculosis and S. typhi are on the increase, and M. tuberculosis and S. typhi infections are often associated with HIV infection. This review discusses the strategies used by these bacteria during persistent infections that allow them to colonize specific sites in the host and evade immune surveillance. The nature of the host immune response to this type of infection and the balance between clearance of the pathogen and avoidance of damage to host tissues are also discussed.     

The chromosome of Salmonella paratyphi A is inverted by recombination between rrnH and rrnG.

Salmonella paratyphi A, a human-adapted bacterial pathogen, causes paratyphoid enteric fever. We established the genome map of strain ATCC 9150 by the use of four endonucleases, XbaI, I-CeuI, AvrII (= BlnI), and SpeI, which generated 27, 7, 19, and 38 fragments, respectively; the sum of the fragments in each case indicates a genome size of ca. 4,600 kb. With phage P22, we transduced Tn10 insertions in known genes from Salmonella typhimurium LT2 to S. paratyphi A ATCC 9150 and located these insertions on the S. paratyphi A chromosome through the XbaI and AvrII sites in Tn10 and through the increased size of the SpeI fragment bearing a Tn10. Compared with the maps of other Salmonella species, the S. paratyphi A genomic map showed two major differences: (i) an insertion of about 100 kb of DNA between rrnH/G and proB and (ii) an inversion of half the genome between rrnH and rrnG, postulated to be due to homologous recombination between the rrn genes. We propose that during the evolution of S. paratyphi A, the first rearrangement event was the 100-kb insertion, which disrupted the chromosomal balance between oriC and the termination of replication, forcing the rrnH/G inversion to restore the balance. The insertion and the inversion are both present in all 10 independent wild-type S. paratyphi A strains tested.     

Genomic Comparison of the Closely-Related Salmonella enterica Serovars Enteritidis,Dublin and Gallinarum

The Salmonella enterica serovars Enteritidis, Dublin, and Gallinarum are closely related but differ in virulence and host range. To identify the genetic elements responsible for these differences and to better understand how these serovars are evolving, we sequenced the genomes of Enteritidis strain LK5 and Dublin strain SARB12 and compared these genomes to the publicly available Enteritidis P125109, Dublin CT 02021853 and Dublin SD3246 genome sequences. We also compared the publicly available Gallinarum genome sequences from biotype Gallinarum 287/91 and Pullorum RKS5078. Using bioinformatic approaches, we identified single nucleotide polymorphisms, insertions, deletions, and differences in prophage and pseudogene content between strains belonging to the same serovar. Through our analysis we also identified several prophage cargo genes and pseudogenes that affect virulence and may contribute to a host-specific, systemic lifestyle. These results strongly argue that the Enteritidis, Dublin and Gallinarum serovars of Salmonella enterica evolve by acquiring new genes through horizontal gene transfer, followed by the formation of pseudogenes. The loss of genes necessary for a gastrointestinal lifestyle ultimately leads to a systemic lifestyle and niche exclusion in the host-specific serovars.     

Web-based visualization tools for bacterial genome alignments

With the increase in the flow of sequence data, both in contigs and whole genomes, visual aids for comparison and analysis studies are becoming imperative. We describe three web-based tools for visualizing alignments of bacterial genomes. The first, called Enteric, produces a graphical, hypertext view of pairwise alignments between a reference genome and sequences from each of several related organisms, covering 20 kb around a user-specified position. Insertions, deletions and rearrangements relative to the reference genome are color-coded, which reveals many intriguing differences among genomes. The second, Menteric, computes and displays nucleotide-level multiple alignments of the same sequences, together with annotations of ORFs and regulatory sites, in a 1 kb region surrounding a given address. The third, a Java-based viewer called Maj, combines some features of the previous tools, and adds a zoom-in mechanism. We compare the Escherichia coli K-12 genome with the partially sequenced genomes of Klebsiella pneumoniae, Yersinia pestis, Vibrio cholerae, and the Salmonella enterica serovars Typhimurium, Typhi and Paratyphi A. Examination of the pairwise and multiple alignments in a region allows one to draw inferences about regulatory patterns and functional assignments. For example, these tools revealed that rffH, a gene involved in enterobacterial common antigen (ECA) biosynthesis, is partly deleted in one of the genomes. We used PCR to show that this deletion occurs sporadically in some strains of some serovars of S.enterica subspecies I but not in any strains tested from six other subspecies. The resulting cell surface diversity may be associated with selection by the host immune response.     

Application of Random Amplified Polymorphic DNA analysis to differentiate strains of Salmonella typhi and other Salmonella species

A Random Amplified Polymorphic DNA (RAPD) fingerprinting method was developed to differentiate isolates of Salmonella serotype typhi ( S. typhi ) and other Salmonella isolates. A panel of five primers was used to examine 63 isolates of Salm. typhi , including 56 strains isolated in Taiwan and seven strains obtained abroad. Twenty-one RAPD types were revealed using the RAPD fingerprinting method. An RAPD with primer 6032 yielded a polymorphism in a 350 bp fragment that differentiated the attenuated vaccine strain Salm. typhi Ty21a from the rest of the Salm. typhi strains. Strains of Salm. typhi were divided into five types with primer D14307. Primer D14307 also proved capable of discrimination among 65 other Salmonella isolates representing 42 different serotypes. The bacterial DNA used in this RAPD protocol was obtained using a commercially available DNA extraction kit (GeneReleaser). The DNA of various strains of Salmonella from this simple extraction procedure could be discriminated within a few hours using the RAPD technique.     

Genotypic homogeneity of multidrug resistant s. Typhimurium infecting distinct adult and childhood susceptibility groups in Blantyre, Malawi

Nontyphoidal Salmonella (NTS) serovars are a common cause of bacteraemia in young children and HIV-infected adults in Malawi and elsewhere in sub-Saharan Africa. These patient populations provide diverse host-immune environments that have the potential to drive bacterial adaptation and evolution. We therefore investigated the diversity of 27 multidrug resistant (MDR) Salmonella Typhimurium strains isolated over 6 years (2002-2008) from HIV-infected adults and children and HIV-uninfected children. Sequence reads from whole-genome sequencing of these isolates using the Illumina GA platform were mapped to the genome of the laboratory strain S. Typhimurium SL1344 excluding homoplastic regions that contained prophage and insertion elements. A phylogenetic tree generated from single nucleotide polymorphisms showed that all 27 strains clustered with the prototypical MDR strain D23580. There was no clustering of strains based on host HIV status or age, suggesting that these susceptible populations acquire S. Typhimurium from common sources or that isolates are transmitted freely between these populations. However, 7/14 of the most recent isolates (2006/2008) formed a distinct clade that branched off 22 SNPs away from the cluster containing earlier isolates. These data suggest that the MDR bacterial population is not static, but is undergoing microevolution which might result in further epidemiology change.     

Novel virulence gene and clustered regularly interspaced short palindromic repeat (CRISPR) multilocus sequence typing scheme for subtyping of the major serovars of Salmonella enterica subsp. enterica

Salmonella enterica subsp. enterica is the leading cause of bacterial food-borne disease in the United States. Molecular subtyping methods are powerful tools for tracking the farm-to-fork spread of food-borne pathogens during outbreaks. In order to develop a novel multilocus sequence typing (MLST) scheme for subtyping the major serovars of S. enterica subsp. enterica, the virulence genes sseL and fimH and clustered regularly interspaced short palindromic repeat (CRISPR) loci were sequenced from 171 clinical isolates from nine Salmonella serovars, Salmonella serovars Typhimurium, Enteritidis, Newport, Heidelberg, Javiana, I 4,[5],12:i:-, Montevideo, Muenchen, and Saintpaul. The MLST scheme using only virulence genes was congruent with serotyping and identified epidemic clones but could not differentiate outbreaks. The addition of CRISPR sequences dramatically improved discriminatory power by differentiating individual outbreak strains/clones. Of particular note, the present MLST scheme provided better discrimination of Salmonella serovar Enteritidis strains than pulsed-field gel electrophoresis (PFGE). This method showed high epidemiologic concordance for all serovars screened except for Salmonella serovar Muenchen. In conclusion, the novel MLST scheme described in the present study accurately differentiated outbreak strains/clones of the major serovars of Salmonella, and therefore, it shows promise for subtyping this important food-borne pathogen during investigations of outbreaks.     

Dynamics of genome rearrangement in bacterial populations

Genome structure variation has profound impacts on phenotype in organisms ranging from microbes to humans, yet little is known about how natural selection acts on genome arrangement. Pathogenic bacteria such as Yersinia pestis, which causes bubonic and pneumonic plague, often exhibit a high degree of genomic rearrangement. The recent availability of several Yersinia genomes offers an unprecedented opportunity to study the evolution of genome structure and arrangement. We introduce a set of statistical methods to study patterns of rearrangement in circular chromosomes and apply them to the Yersinia. We constructed a multiple alignment of eight Yersinia genomes using Mauve software to identify 78 conserved segments that are internally free from genome rearrangement. Based on the alignment, we applied Bayesian statistical methods to infer the phylogenetic inversion history of Yersinia. The sampling of genome arrangement reconstructions contains seven parsimonious tree topologies, each having different histories of 79 inversions. Topologies with a greater number of inversions also exist, but were sampled less frequently. The inversion phylogenies agree with results suggested by SNP patterns. We then analyzed reconstructed inversion histories to identify patterns of rearrangement. We confirm an over-representation of "symmetric inversions"-inversions with endpoints that are equally distant from the origin of chromosomal replication. Ancestral genome arrangements demonstrate moderate preference for replichore balance in Yersinia. We found that all inversions are shorter than expected under a neutral model, whereas inversions acting within a single replichore are much shorter than expected. We also found evidence for a canonical configuration of the origin and terminus of replication. Finally, breakpoint reuse analysis reveals that inversions with endpoints proximal to the origin of DNA replication are nearly three times more frequent. Our findings represent the first characterization of genome arrangement evolution in a bacterial population evolving outside laboratory conditions. Insight into the process of genomic rearrangement may further the understanding of pathogen population dynamics and selection on the architecture of circular bacterial chromosomes.