Phylogenomic analysis identifies gene gains that define Salmonella enterica subspecies I

Salmonella Genomics in Public Health and Food Safety

The species Salmonella enterica comprises over 2,600 serovars, many of which are known to be intracellular pathogens of mammals, birds, and reptiles. It is now apparent that Salmonella is a highly adapted environmental microbe and can readily persist in a number of environmental niches, including water, soil, and various plant (including produce) species. Much of what is known about the evolution and diversity of nontyphoidal Salmonella serovars (NTS) in the environment is the result of the rise of the genomics era in enteric microbiology. There are over 340,000 Salmonella genomes available in public databases. This extraordinary breadth of genomic diversity now available for the species, coupled with widespread availability and affordability of whole-genome sequencing (WGS) instrumentation, has transformed the way in which we detect, differentiate, and characterize Salmonella enterica strains in a timely way. Not only have WGS data afforded a detailed and global examination of the molecular epidemiological movement of Salmonella from diverse environmental reservoirs into human and animal hosts, but they have also allowed considerable consolidation of the diagnostic effort required to test for various phenotypes important to the characterization of Salmonella. For example, drug resistance, serovar, virulence determinants, and other genome-based attributes can all be discerned using a genome sequence. Finally, genomic analysis, in conjunction with functional and phenotypic approaches, is beginning to provide new insights into the precise adaptive changes that permit persistence of NTS in so many diverse and challenging environmental niches.

Diverse distribution of Toxin-Antitoxin II systems in Salmonella enterica serovars

Type II Toxin-Antitoxin systems (TAs), known for their presence in virulent and antibiotic resistant bacterial strains, were recently identified in Salmonella enterica isolates. However, the relationships between the presence of TAs (ccdAB and vapBC) and the epidemiological and genetic features of different non-typhoidal Salmonella serovars are largely unknown, reducing our understanding of the ecological success of different serovars. Salmonella enterica isolates from different sources, belonging to different serovars and epidemiologically unrelated according to ERIC profiles, were investigated for the presence of type II TAs, plasmid content, and antibiotic resistance. The results showed the ubiquitous presence of the vapBC gene in all the investigated Salmonella isolates, but a diverse distribution of ccdAB, which was detected in the most widespread Salmonella serovars, only. Analysis of the plasmid toxin ccdB translated sequence of four selected Salmonella isolates showed the presence of the amino acid substitution R99W, known to impede in vitro the lethal effect of CcdB toxin in the absence of its cognate antitoxin CcdA. These findings suggest a direct role of the TAs in promoting adaptability and persistence of the most prevalent Salmonella serovars, thus implying a wider eco-physiological role for these type II TAs.

Genome-Scale Screening and Validation of Targets for Identification of Salmonella enterica and Serovar Prediction

Salmonella enterica is the most common foodborne pathogen worldwide, with 2,500 recognized serovars. Detection of S. enterica and its classification into serovars are essential for food safety surveillance and clinical diagnosis. The PCR method is useful for these applications because of its rapidity and high accuracy. We obtained 412 candidate detection targets for S. enterica using a comparative genomics mining approach. Gene ontology (GO) functional enrichment analysis of these candidate targets revealed that the GO term with the largest number of unigenes with known function (38 of 177, 21.5%) was significantly involved in pathogenesis (P < 10(-24)). All the candidate targets were then evaluated by PCR assays. Fifteen targets showed high specificity for the detection of S. enterica by verification with 151 S. enterica strains and 34 non-Salmonella strains. The phylogenetic trees of verified targets were highly comparable with those of housekeeping genes, especially for differentiating S. enterica strains into serovars. The serovar prediction ability was validated by sequencing one target (S9) for 39 S. enterica strains belonging to six serovars. Identical mutation sites existed in the same serovar, and different mutation sites were found in diverse serovars. Our findings revealed that 15 verified targets can be potentially used for molecular detection, and some of them can be used for serotyping of S. enterica strains.

Differential single nucleotide polymorphism-based analysis of an outbreak caused by Salmonella enterica serovar Manhattan reveals epidemiological details missed by standard pulsed-field gel electrophoresis

We retrospectively analyzed a rare Salmonella enterica serovar Manhattan outbreak that occurred in Italy in 2009 to evaluate the potential of new genomic tools based on differential single nucleotide polymorphism (SNP) analysis in comparison with the gold standard genotyping method, pulsed-field gel electrophoresis. A total of 39 isolates were analyzed from patients (n=15) and food, feed, animal, and environmental sources (n=24), resulting in five different pulsed-field gel electrophoresis (PFGE) profiles. Isolates epidemiologically related to the outbreak clustered within the same pulsotype, SXB_BS.0003, without any further differentiation. Thirty-three isolates were considered for genomic analysis based on different sets of SNPs, core, synonymous, nonsynonymous, as well as SNPs in different codon positions, by Bayesian and maximum likelihood algorithms. Trees generated from core and nonsynonymous SNPs, as well as SNPs at the second and first plus second codon positions detailed four distinct groups of isolates within the outbreak pulsotype, discriminating outbreak-related isolates of human and food origins. Conversely, the trees derived from synonymous and third-codon-position SNPs clustered food and human isolates together, indicating that all outbreak-related isolates constituted a single clone, which was in line with the epidemiological evidence. Further experiments are in place to extend this approach within our regional enteropathogen surveillance system.

Genetic diversity of Salmonella pathogenicity islands SPI-5 and SPI-6 in Salmonella Newport

Salmonella enterica subspecies enterica serotype Newport is one of the common serotypes causing foodborne salmonellosis outbreaks in the United States. Salmonella Newport consists of three lineages exhibiting extensive genetic diversity. Due to the importance of Salmonella pathogenicity islands 5 and 6 (SPI-5 and SPI-6) in virulence of pathogenic Salmonella, the genetic diversity of these two SPIs may relate to different potentials of Salmonella Newport pathogenicity. Most Salmonella Newport strains from North America belong to Salmonella Newport lineages II and III. A total 28 Salmonella Newport strains of lineages II and III from diverse sources and geographic locations were analyzed, and 11 additional Salmonella genomes were used as outgroup in phylogenetic analyses. SPI-5 was identified in all Salmonella Newport strains and 146 single nucleotide polymorphisms (SNPs) were detected. Thirty-nine lineage-defining SNPs were identified, including 18 nonsynonymous SNPs. Two 40-kb genomic islands (SPI5-GI1 and SPI5-GI2) encoding bacteriophage genes were found between tRNA-ser and pipA. SPI5-GI1 was only present in Salmonella Newport multidrug-resistant strains of lineage II. SPI-6 was found in all strains but three Asian strains in Salmonella Newport lineage II, whereas the three Asian strains carried genomic island SPI6-GI1 at the same locus as SPI-6 in other Salmonella. SPI-6 exhibited 937 SNPs, and phylogenetic analysis demonstrated that clustering of Salmonella Newport isolates was a reflection of their geographic origins. The sequence diversity within SPI-5 and SPI-6 suggests possible recombination events and different virulence potentials of Salmonella Newport. The SNPs could be used as biomarkers during epidemiological investigations.

On the evolutionary history, population genetics and diversity among isolates of Salmonella Enteritidis PFGE pattern JEGX01.0004

Facile laboratory tools are needed to augment identification in contamination events to trace the contamination back to the source (traceback) of Salmonella enterica subsp. enterica serovar Enteritidis (S. Enteritidis). Understanding the evolution and diversity within and among outbreak strains is the first step towards this goal. To this end, we collected 106 new S. Enteriditis isolates within S. Enteriditis Pulsed-Field Gel Electrophoresis (PFGE) pattern JEGX01.0004 and close relatives, and determined their genome sequences. Sources for these isolates spanned food, clinical and environmental farm sources collected during the 2010 S. Enteritidis shell egg outbreak in the United States along with closely related serovars, S. Dublin, S. Gallinarum biovar Pullorum and S. Gallinarum. Despite the highly homogeneous structure of this population, S. Enteritidis isolates examined in this study revealed thousands of SNP differences and numerous variable genes (n = 366). Twenty-one of these genes from the lineages leading to outbreak-associated samples had nonsynonymous (causing amino acid changes) changes and five genes are putatively involved in known Salmonella virulence pathways. While chromosome synteny and genome organization appeared to be stable among these isolates, genome size differences were observed due to variation in the presence or absence of several phages and plasmids, including phage RE-2010, phage P125109, plasmid pSEEE3072_19 (similar to pSENV), plasmid pOU1114 and two newly observed mobile plasmid elements pSEEE1729_15 and pSEEE0956_35. These differences produced modifications to the assembled bases for these draft genomes in the size range of approximately 4.6 to 4.8 mbp, with S. Dublin being larger (∼4.9 mbp) and S. Gallinarum smaller (4.55 mbp) when compared to S. Enteritidis. Finally, we identified variable S. Enteritidis genes associated with virulence pathways that may be useful markers for the development of rapid surveillance and typing methods, potentially aiding in traceback efforts during future outbreaks involving S. Enteritidis PFGE pattern JEGX01.0004.