The 3Cs provide a novel concept of bacterial species: messages from the genome as illustrated by Salmonella

Similar yet different: phylogenomic analysis to delineate Salmonella and Citrobacter species boundaries

Background

Salmonella enterica is a leading cause of foodborne illness worldwide resulting in considerable public health and economic costs. Testing for the presence of this pathogen in food is often hampered by the presence of background microflora that may present as Salmonella (false positives). False positive isolates belonging to the genus Citrobacter can be difficult to distinguish from Salmonella due to similarities in their genetics, cell surface antigens, and other phenotypes. In order to understand the genetic basis of these similarities, a comparative genomic approach was used to define the pan-, core, accessory, and unique coding sequences of a representative population of Salmonella and Citrobacter strains.

Results

Analysis of the genomic content of 58 S. enterica strains and 37 Citrobacter strains revealed the presence of 31,130 and 1540 coding sequences within the pan- and core genome of this population. Amino acid sequences unique to either Salmonella (n = 1112) or Citrobacter (n = 195) were identified and revealed potential niche-specific adaptations. Phylogenetic network analysis of the protein families encoded by the pan-genome indicated that genetic exchange between Salmonella and Citrobacter may have led to the acquisition of similar traits and also diversification within the genera.

Conclusions

Core genome analysis suggests that the Salmonella enterica and Citrobacter populations investigated here share a common evolutionary history. Comparative analysis of the core and pan-genomes was able to define the genetic features that distinguish Salmonella from Citrobacter and highlight niche specific adaptations.

Genetic boundaries delineate the potential human pathogen Salmonella bongori into discrete lineages: divergence and speciation

BACKGROUND

Salmonella bongori infect mainly cold-blooded hosts, but infections by S. bongori in warm-blooded hosts have been reported. We hypothesized that S. bongori might have diverged into distinct phylogenetic lineages, with some being able to infect warm-blooded hosts.

RESULTS

To inspect the divergence status of S. bongori, we first completely sequenced the parakeet isolate RKS3044 and compared it with other sequenced S. bongori strains. We found that RKS3044 contained a novel T6SS encoded in a pathogenicity island-like structure, in addition to a T6SS encoded in SPI-22, which is common to all S. bongori strains so far reported. This novel T6SS resembled the SPI-19 T6SS of the warm-blooded host infecting Salmonella Subgroup I lineages. Genomic sequence comparisons revealed different genomic sequence amelioration events among the S. bongori strains, including a unique CTAG tetranucleotide degeneration pattern in RKS3044, suggesting non-overlapping gene pools between RKS3044 and other S. bongori lineages/strains leading to their independent accumulation of genomic variations. We further proved the existence of a clear-cut genetic boundary between RKS3044 and the other S. bongori lineages/strains analyzed in this study.

CONCLUSIONS

The warm-blooded host-infecting S. bongori strain RKS3044 has diverged with distinct genomic features from other S. bongori strains, including a novel T6SS encoded in a previously not reported pathogenicity island-like structure and a unique genomic sequence degeneration pattern. These findings alert cautions about the emergence of new pathogens originating from non-pathogenic ancestors by acquiring specific pathogenic traits.

Distinct evolutionary origins of common multi-drug resistance phenotypes in Salmonella typhimurium DT104: a convergent process for adaptation under stress

Antimicrobial resistance makes pathogenic bacteria hard to control, but little is known about the general processes of resistance gain or loss. Here, we compared distinct S. typhimurium DT104 strains resistant to zero, two, five, or more of the tested antimicrobials. We found that common resistance phenotypes could be encoded by distinct genes, on SGI-1 or plasmid. We also demonstrated close clonality among all the tested non-resistant and differently resistant DT104 strains, demonstrating dynamic acquisition or loss (by total deletion or gradual decaying of multi-drug resistance gene clusters) of the genetic traits. These findings reflect convergent processes to make the bacteria resistant to multiple antimicrobials by acquiring the needed traits from stochastically available origins. When the antimicrobial stress is absent, the resistance genes may be dropped off quickly, so the bacteria can save the cost for maintaining unneeded genes. Therefore, this work reiterates the importance of strictly controlled use of antimicrobials.

Differential degeneration of the ACTAGT sequence among Salmonella: a reflection of distinct nucleotide amelioration patterns during bacterial divergence

When bacteria diverge, they need to adapt to the new environments, such as new hosts or different tissues of the same host, by accumulating beneficial genomic variations, but a general scenario is unknown due to the lack of appropriate methods. Here we profiled the ACTAGT sequence and its degenerated forms (i.e., hexa-nucleotide sequences with one of the six nucleotides different from ACTAGT) in Salmonella to estimate the nucleotide amelioration processes of bacterial genomes. ACTAGT was mostly located in coding sequences but was also found in several intergenic regions, with its degenerated forms widely scattered throughout the bacterial genomes. We speculated that the distribution of ACTAGT and its degenerated forms might be lineage-specific as a consequence of different selection pressures imposed on ACTAGT at different genomic locations (in genes or intergenic regions) among different Salmonella lineages. To validate this speculation, we modelled the secondary structures of the ACTAGT-containing sequences conserved across Salmonella and many other enteric bacteria. Compared to ACTAGT at conserved regions, the degenerated forms were distributed throughout the bacterial genomes, with the degeneration patterns being highly similar among bacteria of the same phylogenetic lineage but radically different across different lineages. This finding demonstrates biased amelioration under distinct selection pressures among the bacteria and provides insights into genomic evolution during bacterial divergence.

Conserved intergenic sequences revealed by CTAG-profiling in Salmonella: thermodynamic modeling for function prediction

Highly conserved short sequences help identify functional genomic regions and facilitate genomic annotation. We used Salmonella as the model to search the genome for evolutionarily conserved regions and focused on the tetranucleotide sequence CTAG for its potentially important functions. In Salmonella, CTAG is highly conserved across the lineages and large numbers of CTAG-containing short sequences fall in intergenic regions, strongly indicating their biological importance. Computer modeling demonstrated stable stem-loop structures in some of the CTAG-containing intergenic regions, and substitution of a nucleotide of the CTAG sequence would radically rearrange the free energy and disrupt the structure. The postulated degeneration of CTAG takes distinct patterns among Salmonella lineages and provides novel information about genomic divergence and evolution of these bacterial pathogens. Comparison of the vertically and horizontally transmitted genomic segments showed different CTAG distribution landscapes, with the genome amelioration process to remove CTAG taking place inward from both terminals of the horizontally acquired segment.

Characterization of Klebsiella sp. strain S1: a bacterial producer of secoisolariciresinol through biotransformation

Secoisolariciresinol (SECO) is a lignan of potential therapeutic value for diseases such as cancer, but its use has been limited by the lack of ideal production methods, even though its precursors are abundant in plants, such as flaxseeds. Here, we report the characterization of a bacterial strain, S1, isolated from the human intestinal flora, which could produce secoisolariciresinol by biotransformation of precursors in defatted flaxseeds. This bacterium was a Gram-negative and facultatively anaerobic straight rod without capsules. Biochemical assays showed that it was negative for production of oxidase, lysine decarboxylase, ornithine decarboxylase, arginine dihydrolase, and β-glucolase. The G + C content of genomic DNA was 57.37 mol%. Phylogenetic analysis by 16S rRNA and rpoB gene sequences demonstrated S1's close relatedness to Klebsiella. No homologues were found for wzb or wzc (capsular genes), which may explain why Klebsiella sp. strain S1 does not have the capsule and was isolated from a healthy human individual. Based on the percentages of homologous genes with identical nucleotide sequences between the bacteria in comparison, we found that clear-cut genetic boundaries had been formed between S1 and any other Klebsiella strains compared, dividing them into distinct phylogenetic lineages. This work demonstrates that the intestinal Klebsiella, well known as important opportunistic pathogens prevalent in potentially fatal nosocomial infections, may contain lineages that are particularly beneficial to the human health.

Complete genome sequence of Salmonella enterica subspecies arizonae str. RKS2983

Salmonella arizonae (also called Salmonella subgroup IIIa) is a Gram-negative, non-spore-forming, motile, rod-shaped, facultatively anaerobic bacterium. S. arizonae strain RKS2983 was isolated from a human in California, USA. S. arizonae lies somewhere between Salmonella subgroups I (human pathogens) and V (also called S. bongori; usually non-pathogenic to humans) and so is an ideal model organism for studies of bacterial evolution from non-human pathogen to human pathogens. We hence sequenced the genome of RKS2983 for clues of genomic events that might have led to the divergence and speciation of Salmonella into distinct lineages with diverse host ranges and pathogenic features. The 4,574,836 bp complete genome contains 4,203 protein-coding genes, 82 tRNA genes and 7 rRNA operons. This genome contains several characteristics not reported to date in Salmonella subgroup I or V and may provide information about the genetic divergence of Salmonella pathogens.

CTAG-containing cleavage site profiling to delineate Salmonella into natural clusters

Background

The bacterial genus Salmonella contains thousands of serotypes that infect humans or other hosts, causing mild gastroenteritis to potentially fatal systemic infections in humans. Pathogenically distinct Salmonella serotypes have been classified as individual species or as serological variants of merely one or two species, causing considerable confusion in both research and clinical settings. This situation reflects a long unanswered question regarding whether the Salmonella serotypes exist as discrete genetic clusters (natural species) of organisms or as phenotypic (e.g. pathogenic) variants of a single (or two) natural species with a continuous spectrum of genetic divergence among them. Our recent work, based on genomic sequence divergence analysis, has demonstrated that genetic boundaries exist among Salmonella serotypes, circumscribing them into clear-cut genetic clusters of bacteria.

Methodologies/Principal Findings

To further test the genetic boundary concept for delineating Salmonella into clearly defined natural lineages (e.g., species), we sampled a small subset of conserved genomic DNA sequences, i.e., the endonuclease cleavage sites that contain the highly conserved CTAG sequence such as TCTAGA for XbaI. We found that the CTAG-containing cleavage sequence profiles could be used to resolve the genetic boundaries as reliably and efficiently as whole genome sequence comparisons but with enormously reduced requirements for time and resources.

Conclusions

Profiling of CTAG sequence subsets reflects genetic boundaries among Salmonella lineages and can delineate these bacteria into discrete natural clusters.

Genomic comparison of Salmonella typhimurium DT104 with non-DT104 strains

DT104 emerged as a new branch of Salmonella typhimurium with resistance to multiple antimicrobials. To reveal some general genomic features of DT104 for clues of evolutionary events possibly associated with the emergence of this relatively new type of this pathogen, we mapped 11 independent DT104 strains and compared them with non-DT104 S. typhimurium strains. We found that all 11 DT104 strains contained three insertions absent in non-DT104 strains, i.e., the previously reported ST104, ST104B and ST64B. However, SGI-1, a genomic island known to be responsible for DT104 multidrug resistance, was not present in all DT104 strains examined in this study: one DT104 strain did not contain SGI-1 but carried a 144 kb plasmid, suggesting possible evolutionary relationships between the two DNA elements in the development of antimicrobial resistance.

Defining natural species of bacteria: clear-cut genomic boundaries revealed by a turning point in nucleotide sequence divergence

BACKGROUND

Bacteria are currently classified into arbitrary species, but whether they actually exist as discrete natural species was unclear. To reveal genomic features that may unambiguously group bacteria into discrete genetic clusters, we carried out systematic genomic comparisons among representative bacteria.

RESULTS

We found that bacteria of Salmonella formed tight phylogenetic clusters separated by various genetic distances: whereas over 90% of the approximately four thousand shared genes had completely identical sequences among strains of the same lineage, the percentages dropped sharply to below 50% across the lineages, demonstrating the existence of clear-cut genetic boundaries by a steep turning point in nucleotide sequence divergence. Recombination assays supported the genetic boundary hypothesis, suggesting that genetic barriers had been formed between bacteria of even very closely related lineages. We found similar situations in bacteria of Yersinia and Staphylococcus.

CONCLUSIONS

Bacteria are genetically isolated into discrete clusters equivalent to natural species.

Identification of genes to differentiate closely related Salmonella lineages

Background

Salmonella are important human and animal pathogens. Though highly related, the Salmonella lineages may be strictly adapted to different hosts or cause different diseases, from mild local illness like gastroenteritis to fatal systemic infections like typhoid. Therefore, rapid and accurate identification of Salmonella is essential for timely and correct diagnosis of Salmonella infections. The current identification methods such as 16S rRNA sequencing and multilocus sequence typing are expensive and time consuming. Additionally, these methods often do not have sufficient distinguishing resolution among the Salmonella lineages.

Methodologies/Principal Findings

We compared 27 completely sequenced Salmonella genomes to identify possible genomic features that could be used for differentiation of individual lineages. We concatenated 2372 core genes in each of the 27 genomes and constructed a neighbor-joining tree. On the tree, strains of each serotype were clustered tightly together and different serotypes were unambiguously separated with clear genetic distances, demonstrating systematic genomic divergence among the Salmonella lineages. We made detailed comparisons among the 27 genomes and identified distinct sets of genomic differences, including nucleotide variations and genomic islands (GIs), among the Salmonella lineages. Two core genes STM4261 and entF together could unambiguously distinguish all Salmonella lineages compared in this study. Additionally, strains of a lineage have a common set of GIs and closely related lineages have similar sets of GIs.

Conclusions

Salmonella lineages have accumulated distinct sets of mutations and laterally acquired DNA (e.g., GIs) in evolution. Two genes entF and STM4261 have diverged sufficiently among the Salmonella lineages to be used for their differentiation. Further investigation of the distinct sets of mutations and GIs will lead to novel insights into genomic evolution of Salmonella and greatly facilitate the elucidation of pathogeneses of Salmonella infections.