Substructure within Salmonella enterica subsp. enterica isolates from Australian wildlife

EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

The definition of bacterial species is traditionally a taxonomic issue while bacterial populations are identified by population genetics. These assignments are species specific, and depend on the practitioner. Legacy multilocus sequence typing is commonly used to identify sequence types (STs) and clusters (ST Complexes). However, these approaches are not adequate for the millions of genomic sequences from bacterial pathogens that have been generated since 2012. EnteroBase (http://enterobase.warwick.ac.uk) automatically clusters core genome MLST allelic profiles into hierarchical clusters (HierCC) after assembling annotated draft genomes from short-read sequences. HierCC clusters span core sequence diversity from the species level down to individual transmission chains. Here we evaluate HierCC's ability to correctly assign 100 000s of genomes to the species/subspecies and population levels for Salmonella, Escherichia, Clostridoides, Yersinia, Vibrio and Streptococcus. HierCC assignments were more consistent with maximum-likelihood super-trees of core SNPs or presence/absence of accessory genes than classical taxonomic assignments or 95% ANI. However, neither HierCC nor ANI were uniformly consistent with classical taxonomy of Streptococcus. HierCC was also consistent with legacy eBGs/ST Complexes in Salmonella or Escherichia and with O serogroups in Salmonella. Thus, EnteroBase HierCC supports the automated identification of and assignment to species/subspecies and populations for multiple genera. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

BASELINE HEALTH PARAMETERS FOR A NEWLY ESTABLISHED POPULATION OF LONG-NOSED POTOROO (POTOROUS TRIDACTYLUS) AT BOODEREE NATIONAL PARK, AUSTRALIA

Over two field seasons during 2014-15, 35 long-nosed potoroos (Potorous tridactylus) were captured in state forests in South Eastern New South Wales for translocation to Booderee National Park, Jervis Bay Territory, Australia. Animals were anesthetized for physical examination and collection of samples to assess general health and screen for select diseases identified during a disease risk assessment. Morphologic, hematologic, and biochemical parameters were determined, and parasites were identified where possible. Trypanosoma gilletti, Trypanosoma vegrandis, and novel genotypes most similar to a Trypanosoma wallaby-derived isolate (ABF) were identified from blood samples by PCR; the first time Trypanosoma has been described in this species. Also reported is the first confirmation of the Australian paralysis tick, Ixodes holocyclus, from the long-nosed potoroo. Surveillance showed that Cryptococcus sp. may form part of the normal nasal flora for long-nosed potoroo. Salmonella enterica serotype Dublin and Salmonella enterica subsp. enterica was identified from rectal swabs of otherwise healthy animals. The data provide baseline health and disease parameters for this newly established population and the source population and will inform future translocation and conservation management activities. These data expand current knowledge on aspects of the biology and microbiology of the long-nosed potoroo, both locally and nationally.

In vitro invasiveness and antimicrobial resistance of Salmonella enterica subspecies isolated from wild and captive reptiles

Reptiles are carriers of Salmonella and can intermittently shed bacteria in their faeces. Contact with snakes and lizards is a source of human salmonellosis. Here, two populations of reptiles, wild and captive were surveyed for Salmonella. One hundred thirty wild-caught reptiles were sampled for Salmonella including 2 turtle, 9 snake and 31 lizard species. Fifty-two of 130 (40%) animals were Salmonella positive: one of 5 (20%) turtles, 7 of 14 (50%) snakes and 44 of 111 (39.6%) lizards. One hundred twenty-two reptiles were sampled from a zoo collection including 1 turtle, 6 tortoise, 9 lizard, 14 snake and 1 crocodile species. Forty-two of 122 (34.4%) captive reptiles sampled were Salmonella positive. Salmonella was most commonly isolated from lizards and snakes. Fifteen serotypes were identified from zoo and 19 from wild-caught reptiles and most were members of subspecies enterica (I), salamae (II), arizonae (IIIa) or diarizonae (IIIb). Antimicrobial susceptibility testing was conducted on all Salmonella isolates; only two exhibited resistance, a Salmonella subsp. (II) ser. 21:z₁₀ :z₆ (Wandsbek) isolate cultured from a wild-caught reptile and a Salmonella Typhimurium DT120 isolated from a captive snake. The invasive capacity of reptile-associated Salmonella strains into cultured human intestinal epithelial (Caco2) and mouse macrophages cell lines (J774A.1) was also investigated. All isolates were invasive into both cell lines. Significant (P ≤ 0.001) variability in invasiveness into polarized Caco2 cells was observed. Salmonella Eastbourne exhibited the highest invasiveness into Caco2 cells and Salmonella Chester the lowest, with mean per cent recoveries of 19.99 ± 0.32 and 1.23 ± 0.30, respectively. Invasion into J774A.1 macrophages was also variable but was not significant. Salmonella subsp. II ser. 17:g,t:- (Bleadon) exhibited the highest invasiveness into J774A.1 with a mean per cent recovery of 10.19 ± 0.19. Thus, reptile-associated Salmonellae are likely to have different capacities to cause disease in humans.

Genomic diversity of Salmonella enterica -The UoWUCC 10K genomes project

Background: Most publicly available genomes of Salmonella enterica are from human disease in the US and the UK, or from domesticated animals in the US. Methods: Here we describe a historical collection of 10,000 strains isolated between 1891-2010 in 73 different countries. They encompass a broad range of sources, ranging from rivers through reptiles to the diversity of all S. enterica isolated on the island of Ireland between 2000 and 2005. Genomic DNA was isolated, and sequenced by Illumina short read sequencing. Results: The short reads are publicly available in the Short Reads Archive. They were also uploaded to EnteroBase, which assembled and annotated draft genomes. 9769 draft genomes which passed quality control were genotyped with multiple levels of multilocus sequence typing, and used to predict serovars. Genomes were assigned to hierarchical clusters on the basis of numbers of pair-wise allelic differences in core genes, which were mapped to genetic Lineages within phylogenetic trees. Conclusions: The University of Warwick/University College Cork (UoWUCC) project greatly extends the geographic sources, dates and core genomic diversity of publicly available S. enterica genomes. We illustrate these features by an overview of core genomic Lineages within 33,000 publicly available Salmonella genomes whose strains were isolated before 2011. We also present detailed examinations of HC400, HC900 and HC2000 hierarchical clusters within exemplar Lineages, including serovars Typhimurium, Enteritidis and Mbandaka. These analyses confirm the polyphyletic nature of multiple serovars while showing that discrete clusters with geographical specificity can be reliably recognized by hierarchical clustering approaches. The results also demonstrate that the genomes sequenced here provide an important counterbalance to the sampling bias which is so dominant in current genomic sequencing.

Genomic diversity of Salmonella enterica -The UoWUCC 10K genomes project

Background: Most publicly available genomes of Salmonella enterica are from human disease in the US and the UK, or from domesticated animals in the US. Methods: Here we describe a historical collection of 10,000 strains isolated between 1891-2010 in 73 different countries. They encompass a broad range of sources, ranging from rivers through reptiles to the diversity of all S. enterica isolated on the island of Ireland between 2000 and 2005. Genomic DNA was isolated, and sequenced by Illumina short read sequencing. Results: The short reads are publicly available in the Short Reads Archive. They were also uploaded to EnteroBase, which assembled and annotated draft genomes. 9769 draft genomes which passed quality control were genotyped with multiple levels of multilocus sequence typing, and used to predict serovars. Genomes were assigned to hierarchical clusters on the basis of numbers of pair-wise allelic differences in core genes, which were mapped to genetic Lineages within phylogenetic trees. Conclusions: The University of Warwick/University College Cork (UoWUCC) project greatly extends the geographic sources, dates and core genomic diversity of publicly available S. enterica genomes. We illustrate these features by an overview of core genomic Lineages within 33,000 publicly available Salmonella genomes whose strains were isolated before 2011. We also present detailed examinations of HC400, HC900 and HC2000 hierarchical clusters within exemplar Lineages, including serovars Typhimurium, Enteritidis and Mbandaka. These analyses confirm the polyphyletic nature of multiple serovars while showing that discrete clusters with geographical specificity can be reliably recognized by hierarchical clustering approaches. The results also demonstrate that the genomes sequenced here provide an important counterbalance to the sampling bias which is so dominant in current genomic sequencing.

Salmonella enterica Serovar Panama, an Understudied Serovar Responsible for Extraintestinal Salmonellosis Worldwide

In recent years nontyphoidal Salmonella has emerged as one of the pathogens most frequently isolated from the bloodstream in humans. Only a small group of Salmonella serovars cause this systemic infection, known as invasive nontyphoidal salmonellosis. Here, we present a focused minireview on Salmonella enterica serovar Panama, a serovar responsible for invasive salmonellosis worldwide. S Panama has been linked with infection of extraintestinal sites in humans, causing septicemia, meningitis, and osteomyelitis. The clinical picture is often complicated by antimicrobial resistance and has been associated with a large repertoire of transmission vehicles, including human feces and breast milk. Nonhuman sources of S Panama involve reptiles and environmental reservoirs, as well as food animals, such as pigs. The tendency of S Panama to cause invasive disease may be linked to certain serovar-specific genetic factors.

The diversity, evolution and ecology of Salmonella in venomous snakes

BACKGROUND

Reptile-associated Salmonella bacteria are a major, but often neglected cause of both gastrointestinal and bloodstream infection in humans globally. The diversity of Salmonella enterica has not yet been determined in venomous snakes, however other ectothermic animals have been reported to carry a broad range of Salmonella bacteria. We investigated the prevalence and diversity of Salmonella in a collection of venomous snakes and non-venomous reptiles.

METHODOLOGY/PRINCIPLE FINDINGS

We used a combination of selective enrichment techniques to establish a unique dataset of reptilian isolates to study Salmonella enterica species-level evolution and ecology and used whole-genome sequencing to investigate the relatedness of phylogenetic groups. We observed that 91% of venomous snakes carried Salmonella, and found that a diverse range of serovars (n = 58) were carried by reptiles. The Salmonella serovars belonged to four of the six Salmonella enterica subspecies: diarizonae, enterica, houtanae and salamae. Subspecies enterica isolates were distributed among two distinct phylogenetic clusters, previously described as clade A (52%) and clade B (48%). We identified metabolic differences between S. diarizonae, S. enterica clade A and clade B involving growth on lactose, tartaric acid, dulcitol, myo-inositol and allantoin.

SIGNIFICANCE

We present the first whole genome-based comparative study of the Salmonella bacteria that colonise venomous and non-venomous reptiles and shed new light on Salmonella evolution. Venomous snakes examined in this study carried a broad range of Salmonella, including serovars which have been associated with disease in humans such as S. Enteritidis. The findings raise the possibility that venomous snakes could be a reservoir for Salmonella serovars associated with human salmonellosis.

Salmonella enterica isolates from Western Australian rangeland goats remain susceptible to critically important antimicrobials

This study investigated faecal carriage and antimicrobial resistance (AMR) of Salmonella enterica recovered from rangeland goats. Faecal samples (n = 400) were collected at slaughter from four consignments of goats (n = 100 samples per consignment), each from one of four localities in Western Australia. Carriage of Salmonella spp. was detected in 106 samples (26.5%; 95% CI 22.4-31.0%). The rate of faecal carriage for each consignment ranged between 23-30%. PCR assays targeting the STM2755 and STM4497 genes revealed 84.9% (90/106) of the isolates were of serovar Typhimurium. Salmonella Chester (11/106, 10.4%) and S. Saintpaul (5/106, 4.7%) were characterised at invA and ompF genes. Antimicrobial susceptibility testing demonstrated that 84.0% of isolates were susceptible to all tested (n = 13) antimicrobials. Resistance was identified to azithromycin (14.2%), tetracycline (10.4%), ampicillin (5.7%), amoxicillin-clavulanate and cefoxitin (3.8%), trimethoprim/sulfamethoxazole (1.9%), gentamicin and streptomycin (0.9%). No isolate was resistant to four or more antimicrobials, or to critically important antimicrobials such as fluoroquinolones and extended spectrum cephalosporins. This is the first study reporting AMR in Salmonella isolates from Australian rangeland goats. The rate of detection of AMR was very low, some resistance to low-importance drugs was present in the Salmonella population, despite the absence of active selection pressure.

Spatial Variation and Survival of Salmonella enterica Subspecies in a Population of Australian Sleepy Lizards (Tiliqua rugosa)

The life cycles of many enteric bacterial species require a transition between two very distinct environments. Their primary habitat is the gastrointestinal tract of the host, while their secondary habitat, during transmission from one host to another, consists of environments external to the host, such as soil, water, and sediments. Consequently, both host and environmental factors shape the genetic structure of enteric bacterial populations. This study examined the distribution of four Salmonella enterica subspecies in a population of sleepy lizards, Tiliqua rugosa, in a semiarid region of South Australia. The lizards living within the 1,920-m by 720-m study site were radio tracked, and their enteric bacteria were sampled at regular intervals throughout their active seasons in the years 2001, 2002, and 2006. Four of the six subspecies of S. enterica were present in this population and were nonrandomly distributed among the lizards. In particular, S. enterica subsp. diarizonae was restricted to lizards living in the most shaded parts of the study site with an overstorey of Casuarina trees. Experiments undertaken to investigate the survival of S. enterica cells under seminatural conditions revealed that cell survival decreased with increased exposure to elevated temperatures and UV light. Among the three S. enterica subspecies tested, S. enterica subsp. diarizonae consistently had an average expected life span that was shorter than that observed for the other two subspecies. There was no indication in the data that there was any competitive dominance hierarchy among the S. enterica subspecies within individual hosts. Thus, the nonrandom distribution of S. enterica subspecies in this population of lizards appears to be driven by their different survival characteristics in the external environment.

A Comparison of Subtyping Methods for Differentiating Salmonella enterica Serovar Enteritidis Isolates Obtained from Food and Human Sources

Purpose:

To evaluate the abilities of these subtyping methods, we distinguished Salmonella Enteritidis (S. Enteritidis) isolated from food products and human clinical samples between 2009 and 2010 in Seoul using five subtyping methods.

Methods:

We determined the subtypes of 20 S. Enteritidis isolates from food and human sources using phage typing, antimicrobial susceptibility, pulsed-field gel electrophoresis (PFGE), repetitive sequence-based PCR (rep-PCR), and multi-locus sequence typing (MLST).

Results:

A total of 20 tested isolates were differentiated into six antimicrobial susceptibility patterns, three different phage types, four different PFGE profiles, seven rep-PCR patterns, and one MLST type. Food isolates were considerably more susceptible to antibiotics than human isolates. We were best able to discriminate among S. Enteritidis isolates using rep-PCR, and obtained the highest Simpson’s diversity index of 0.82, whereas other methods produced indices that were less than 0.71. PFGE pattern appeared to be more related to antimicrobial resistance and phage types of S. Enteritidis isolates than rep-PCR. MLST revealed identical alleles in all isolates at all seven loci examined, indicating no resolution.

Conclusion:

The results of this study suggest that rep-PCR provided the best discriminatory power for phenotypically similar S. Enteritidis isolates of food and human origins, whereas the discriminatory ability of MLST may be problematic because of the high sequence conservation of the targeted genes.

Multilocus sequence typing as a replacement for serotyping in Salmonella enterica

Salmonella enterica subspecies enterica is traditionally subdivided into serovars by serological and nutritional characteristics. We used Multilocus Sequence Typing (MLST) to assign 4,257 isolates from 554 serovars to 1092 sequence types (STs). The majority of the isolates and many STs were grouped into 138 genetically closely related clusters called eBurstGroups (eBGs). Many eBGs correspond to a serovar, for example most Typhimurium are in eBG1 and most Enteritidis are in eBG4, but many eBGs contained more than one serovar. Furthermore, most serovars were polyphyletic and are distributed across multiple unrelated eBGs. Thus, serovar designations confounded genetically unrelated isolates and failed to recognize natural evolutionary groupings. An inability of serotyping to correctly group isolates was most apparent for Paratyphi B and its variant Java. Most Paratyphi B were included within a sub-cluster of STs belonging to eBG5, which also encompasses a separate sub-cluster of Java STs. However, diphasic Java variants were also found in two other eBGs and monophasic Java variants were in four other eBGs or STs, one of which is in subspecies salamae and a second of which includes isolates assigned to Enteritidis, Dublin and monophasic Paratyphi B. Similarly, Choleraesuis was found in eBG6 and is closely related to Paratyphi C, which is in eBG20. However, Choleraesuis var. Decatur consists of isolates from seven other, unrelated eBGs or STs. The serological assignment of these Decatur isolates to Choleraesuis likely reflects lateral gene transfer of flagellar genes between unrelated bacteria plus purifying selection. By confounding multiple evolutionary groups, serotyping can be misleading about the disease potential of S. enterica. Unlike serotyping, MLST recognizes evolutionary groupings and we recommend that Salmonella classification by serotyping should be replaced by MLST or its equivalents.

Microbiologists have used serological and nutritional characteristics to subdivide pathogenic bacteria for nearly 100 years. These subdivisions in Salmonella enterica are called serovars, some of which are thought to be associated with particular diseases and epidemiology. We used MultiLocus Sequence-based Typing (MLST) to identify clusters of S. enterica isolates that are related by evolutionary descent. Some clusters correspond to serovars on a one to one basis. But many clusters include multiple serovars, which is of public health significance, and most serovars span multiple, unrelated clusters. Despite its broad usage, serological typing of S. enterica has resulted in confusing systematics, with a few exceptions. We recommend that serotyping for strain discrimination of S. enterica be replaced by a DNA-based method, such as MLST. Serotyping and other non-sequence based typing methods are routinely used for detecting outbreaks and to support public health responses. Moving away from these methods will require a major shift in thinking by public health microbiology laboratories as well as national and international agencies. However, a transition to the routine use of MLST, supplemented where appropriate by even more discriminatory sequence-based typing methods based on entire genomes, will provide a clearer picture of long-term transmission routes of Salmonella, facilitate data transfer and support global control measures.

Genome sequencing reveals diversification of virulence factor content and possible host adaptation in distinct subpopulations of Salmonella enterica

Background

Divergence of bacterial populations into distinct subpopulations is often the result of ecological isolation. While some studies have suggested the existence of Salmonella enterica subsp. enterica subclades, evidence for these subdivisions has been ambiguous. Here we used a comparative genomics approach to define the population structure of Salmonella enterica subsp. enterica, and identify clade-specific genes that may be the result of ecological specialization.

Results

Multi-locus sequence analysis (MLSA) and single nucleotide polymorphisms (SNPs) data for 16 newly sequenced and 30 publicly available genomes showed an unambiguous subdivision of S. enterica subsp. enterica into at least two subpopulations, which we refer to as clade A and clade B. Clade B strains contain several clade-specific genes or operons, including a β-glucuronidase operon, a S-fimbrial operon, and cell surface related genes, which strongly suggests niche specialization of this subpopulation. An additional set of 123 isolates was assigned to clades A and B by using qPCR assays targeting subpopulation-specific SNPs and genes of interest. Among 98 serovars examined, approximately 20% belonged to clade B. All clade B isolates contained two pathogenicity related genomic islands, SPI-18 and a cytolethal distending toxin islet; a combination of these two islands was previously thought to be exclusive to serovars Typhi and Paratyphi A. Presence of β-glucuronidase in clade B isolates specifically suggests an adaptation of this clade to the vertebrate gastrointestinal environment.

Conclusions

S. enterica subsp. enterica consists of at least two subpopulations that differ specifically in genes involved in host and tissue tropism, utilization of host specific carbon and nitrogen sources and are therefore likely to differ in ecology and transmission characteristics.