Shu-Lin Liu

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.

Salmonella paratyphi C: Genetic Divergence from Salmonella choleraesuis and Pathogenic Convergence with Salmonella typhi

Background Although over 1400 Salmonella serovars cause usually self-limited gastroenteritis in humans, a few, e.g., Salmonella typhi and S. paratyphi C, cause typhoid, a potentially fatal systemic infection. It is not known whether the typhoid agents have evolved from a common ancestor (by divergent processes) or acquired similar pathogenic traits independently (by convergent processes). Comparison of different typhoid agents with non-typhoidal Salmonella lineages will provide excellent models for studies on how similar pathogens might have evolved. Methodologies/Principal Findings We sequenced a strain of S. paratyphi C, RKS4594, and compared it with previously sequenced Salmonella strains. RKS4594 contains a chromosome of 4,833,080 bp and a plasmid of 55,414 bp. We predicted 4,640 intact coding sequences (4,578 in the chromosome and 62 in the plasmid) and 152 pseudogenes (149 in the chromosome and 3 in the plasmid). RKS4594 shares as many as 4346 of the 4,640 genes with a strain of S. choleraesuis, which is primarily a swine pathogen, but only 4008 genes with another human-adapted typhoid agent, S. typhi. Comparison of 3691 genes shared by all six sequenced Salmonella strains placed S. paratyphi C and S. choleraesuis together at one end, and S. typhi at the opposite end, of the phylogenetic tree, demonstrating separate ancestries of the human-adapted typhoid agents. S. paratyphi C seemed to have suffered enormous selection pressures during its adaptation to man as suggested by the differential nucleotide substitutions and different sets of pseudogenes, between S. paratyphi C and S. choleraesuis. Conclusions S. paratyphi C does not share a common ancestor with other human-adapted typhoid agents, supporting the convergent evolution model of the typhoid agents. S. paratyphi C has diverged from a common ancestor with S. choleraesuis by accumulating genomic novelty during adaptation to man.

Salmonella enterica Serovar Typhi Strains from Which SPI7, a 134-Kilobase Island with Genes for Vi Exopolysaccharide and Other Functions, Has Been Deleted

Salmonella enterica serovar Typhi has a 134-kb island of DNA identified as salmonella pathogenicity island 7 (SPI7), inserted between pheU and pheU (truncated), two genes for tRNA(Phe). SPI7 has genes for Vi exopolysaccharide, for type IVB pili, for putative conjugal transfer, and for sopE bacteriophage. Pulsed-field gel electrophoresis following digestion with the endonuclease I-CeuI, using DNA from a set of 120 wild-type strains of serovar Typhi assembled from several sources, identified eight strains in which the I-CeuI G fragment, which contains SPI7, had a large deletion. In addition, agglutination tests with Vi antiserum and phage typing with Vi phages show that all eight strains are Vi negative. We therefore tested these strains for deletion of SPI7 by multiplex PCR, by microarray analysis, and by sequencing of PCR amplicons. Data show that seven of the eight strains are precise deletions of SPI7: a primer pair flanking SPI7 results in a PCR amplicon containing a single pheU gene; microarrays show that all SPI7 genes are deleted. Two of the strains produce amplicons which have A derived from pheU at bp 27, while five have C derived from pheU at this position; thus, the position of the crossover which results in the deletion can be inferred. The deletion in the eighth strain, TYT1669, removes 175 kb with junction points in genes STY4465 and STY4664; the left junction of SPI7 and adjacent genes, as well as part of SPI7 including the viaB operon for Vi exopolysaccharide, was removed, while the right junction of SPI7 was retained. We propose that these deletions occurred during storage following isolation.

The XbaI‐BlnI‐CeuI genomic cleavage map of Salmonella enteritidis shows an inversion relative to Salmonella typhimurium LT2

We have established the genomic cleavage map of Salmonella enteritidis strain SSU7998 using pulsed‐field gel electrophoresis. The chromosome of 4600kb was analysed by XbaI (16 fragments), I‐CeuI (7 fragments) and BlnI (12 fragments); the genome also contains a plasmid of 60 kb. Cleavage sites of I‐CeuI, in the large subunit ribosomal RNA gene, are conserved from Salmonella typhimurium and Escherichia coli K‐12, and the XbaI and BinI sites in glt‐tRNA are also conserved, but other sites are less conserved. Transposon Tn10, located at 60 different positions in the chromosome of S. typhimurium, was transduced by bacteriophage P22 into S. enteritidis and the insertion mapped using the XbaI and BlnI sites on Tn10. Gene order in S. enteritidis is identical to S. typhimurium LT2 and similar to E. coli K‐12 except for an inversion of 815 kb, which covers the terminus region including T1 and T2. Endpoints are in the NDZs, or non‐divisible zones, in which inversion endpoints were not detected in experiments in E. coli K‐12 and S. typhimurium LT2. This inversion resembles the inversion between S. typhimurium and E. coli, but is longer at both ends.

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.

Wavelet to predict bacterial ori and ter: a tendency towards a physical balance.

BackgroundChromosomal DNA replication in bacteria starts at the origin (ori) and the two replicores propagate in opposite directions up to the terminus (ter) region. We hypothesize that the two replicores need to reach ter at the same time to maintain a physical balance; DNA insertion would disrupt such a balance, requiring chromosomal rearrangements to restore the balance. To test this hypothesis, we needed to demonstrate that ori and ter are in a physical balance in bacterial chromosomes. Using wavelet analysis, we documented GC skew, AT skew, purine excess and keto excess on the published bacterial genomic sequences to locate the turning (minimum and maximum) points on the curves. Previously, the minimum point had been supposed to correlate with ori and the maximum to correlate with ter.ResultsWe observed a strong tendency of the bacterial chromosomes towards a physical balance, with the minima and maxima corresponding to the known or putative ori and ter and being about half chromosome separated in most of the bacteria studied. A nonparametric method based on wavelet transformation was employed to perform significance tests for the predicted loci.ConclusionsThe wavelet approach can reliably predict the ori and ter regions and the bacterial chromosomes have a strong tendency towards a physical balance between ori and ter.

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.

Gene Decay in Shigella as an Incipient Stage of Host-Adaptation

Background Many facultative bacterial pathogens have undergone extensive gene decay processes, possibly due to lack of selection pressure during evolutionary conversion from free-living to intracellular lifestyle. Shigella, the causative agents of human shigellosis, have arisen from different E. coli-like ancestors independently by convergent paths. As these bacteria all have lost large numbers of genes by mutation or deletion, they can be used as ideal models for systematically studying the process of gene function loss in different bacteria living under similar selection pressures. Methodologies/Principal Findings We compared the sequenced Shigella genomes and re-defined decayed genes (pseudogenes plus deleted genes) in these bacteria. Altogether, 85 genes are commonly decayed in the five analyzed Shigella strains and 1456 genes are decayed in at least one Shigella strain. Genes coding for carbon utilization, cell motility, transporter or membrane proteins are prone to be inactivated. Decayed genes tend to concentrate in certain operons rather than distribute averagely across the whole genome. Genes in the decayed operon accumulated more non-synonymous mutations than the rest genes and meanwhile have lower expression levels. Conclusions Different Shigella lineages underwent convergent gene decay processes, and inactivation of one gene would lead to a lesser selection pressure for the other genes in the same operon. The pool of superfluous genes for Shigella may contain at least two thousand genes and the gene decay processes may still continue in Shigella until a minimum genome harboring only essential genes is reached.

The Genome of Salmonella enterica Serovar Gallinarum: Distinct Insertions/Deletions and Rare Rearrangements

Salmonella enterica serovar Gallinarum is a fowl-adapted pathogen, causing typhoid fever in chickens. It has the same antigenic formula (1,9,12:--:--) as S. enterica serovar Pullorum, which is also adapted to fowl but causes pullorum disease (diarrhea). The close relatedness but distinct pathogeneses make this pair of fowl pathogens good models for studies of bacterial genomic evolution and the way these organisms acquired pathogenicity. To locate and characterize the genomic differences between serovar Gallinarum and other salmonellae, we constructed a physical map of serovar Gallinarum strain SARB21 by using I-CeuI, XbaI, and AvrII with pulsed-field gel electrophoresis techniques. In the 4,740-kb genome, we located two insertions and six deletions relative to the genome of S. enterica serovar Typhimurium LT2, which we used as a reference Salmonella genome. Four of the genomic regions with reduced lengths corresponded to the four prophages in the genome of serovar Typhimurium LT2, and the others contained several smaller deletions relative to serovar Typhimurium LT2, including regions containing srfJ, std, and stj and gene clusters encoding a type I restriction system in serovar Typhimurium LT2. The map also revealed some rare rearrangements, including two inversions and several translocations. Further characterization of these insertions, deletions, and rearrangements will provide new insights into the molecular basis for the specific host-pathogen interactions and mechanisms of genomic evolution to create a new pathogen.

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.