Henk C. den Bakker

Listeria monocytogenes lineages: Genomics, evolution, ecology, and phenotypic characteristics

Listeria monocytogenes consists of at least 4 evolutionary lineages (I, II, III, and IV) with different but overlapping ecological niches. Most L. monocytogenes isolates seem to belong to lineages I and II, which harbor the serotypes more commonly associated with human clinical cases, including serotype 1/2a (lineage II) and serotypes 1/2b and 4b (lineage I). Lineage II strains are common in foods, seem to be widespread in the natural and farm environments, and are also commonly isolated from animal listeriosis cases and sporadic human clinical cases. Most human listeriosis outbreaks are associated with lineage I isolates though. In addition, a number of studies indicate that, in many countries, lineage I strains are overrepresented among human isolates, as compared to lineage II strains. Lineage III and IV strains on the other hand are rare and predominantly isolated from animal sources. The apparent differences in the distribution of strains representing the L. monocytogenes lineages has lead to a number of studies aimed at identifying phenotypic differences among the different lineages. Interestingly, lineage II isolates seem to carry more plasmids than lineage I isolates and these plasmids often confer resistance to toxic metals and possibly other compounds that may be found in the environment. Moreover, lineage II isolates seem to be more resistant to bacteriocins than lineage I isolates, which probably confers an advantage in environments where bacteriocin-producing organisms are abundant. A large number of lineage II isolates and strains have been shown to be virulence-attenuated due to premature stop codon mutations in inlA and mutations in prfA. A subset of lineage I isolates carry a listeriolysin S hemolysin, which is not present in isolates belonging to lineages II, III, or IV. While lineage II isolates also show higher recombination rates than lineage I isolates, possibly facilitating adaptation of lineage II strains to diverse environments, lineage I isolates are clonal and show a low prevalence of plasmids and IS elements, suggesting that lineage I isolates may have mechanisms that limit the acquisition of foreign DNA by horizontal gene transfer. Diversifying selection has also been shown to have played an important role during evolution of the L. monocytogenes lineages and during divergence of L. monocytogenes from the non-pathogenic species L. innocua. Overall evidence thus suggests that the 4 L. monocytogenes lineages identified so far represent distinct ecologic, genetic, and phenotypic characteristics, which appear to affect their ability to be transmitted through foods and to cause human disease. Further insights into the ecology, evolution, and characteristics of these lineages will thus not only provide an improved understanding of the evolution of this foodborne pathogen, but may also facilitate improved control of foodborne listeriosis.

Listeria marthii sp. nov., isolated from the natural environment, Finger Lakes National Forest.

Four isolates (FSL S4-120(T), FSL S4-696, FSL S4-710, and FSL S4-965) of Gram-positive, motile, facultatively anaerobic, non-spore-forming bacilli that were phenotypically similar to species of the genus Listeria were isolated from soil, standing water and flowing water samples obtained from the natural environment in the Finger Lakes National Forest, New York, USA. The four isolates were closely related to one another and were determined to be the same species by whole genome DNA-DNA hybridization studies (>82 % relatedness at 55 degrees C and >76 % relatedness at 70 degrees C with 0.0-0.5 % divergence). 16S rRNA gene sequence analysis confirmed their close phylogenetic relatedness to Listeria monocytogenes and Listeria innocua and more distant relatedness to Listeria welshimeri, L. seeligeri, L. ivanovii and L. grayi. Phylogenetic analysis of partial sequences for sigB, gap, and prs showed that these isolates form a well-supported sistergroup to L. monocytogenes. The four isolates were sufficiently different from L. monocytogenes and L. innocua by DNA-DNA hybridization to warrant their designation as a new species of the genus Listeria. The four isolates yielded positive reactions in the AccuProbe test that is purported to be specific for L. monocytogenes, did not ferment L-rhamnose, were non-haemolytic on blood agar media, and did not contain a homologue of the L. monocytogenes virulence gene island. On the basis of their phenotypic characteristics and their genotypic distinctiveness from L. monocytogenes and L. innocua, the four isolates should be classified as a new species within the genus Listeria, for which the name Listeria marthii sp. nov. is proposed. The type strain of L. marthii is FSL S4-120(T) (=ATCC BAA-1595(T) =BEIR NR 9579(T) =CCUG 56148(T)). L. marthii has not been associated with human or animal disease at this time.

Comparative genomics of the bacterial genus Listeria: Genome evolution is characterized by limited gene acquisition and limited gene loss

BackgroundThe bacterial genus Listeria contains pathogenic and non-pathogenic species, including the pathogens L. monocytogenes and L. ivanovii, both of which carry homologous virulence gene clusters such as the prfA cluster and clusters of internalin genes. Initial evidence for multiple deletions of the prfA cluster during the evolution of Listeria indicates that this genus provides an interesting model for studying the evolution of virulence and also presents practical challenges with regard to definition of pathogenic strains.ResultsTo better understand genome evolution and evolution of virulence characteristics in Listeria, we used a next generation sequencing approach to generate draft genomes for seven strains representing Listeria species or clades for which genome sequences were not available. Comparative analyses of these draft genomes and six publicly available genomes, which together represent the main Listeria species, showed evidence for (i) a pangenome with 2,032 core and 2,918 accessory genes identified to date, (ii) a critical role of gene loss events in transition of Listeria species from facultative pathogen to saprotroph, even though a consistent pattern of gene loss seemed to be absent, and a number of isolates representing non-pathogenic species still carried some virulence associated genes, and (iii) divergence of modern pathogenic and non-pathogenic Listeria species and strains, most likely circa 47 million years ago, from a pathogenic common ancestor that contained key virulence genes.ConclusionsGenome evolution in Listeria involved limited gene loss and acquisition as supported by (i) a relatively high coverage of the predicted pan-genome by the observed pan-genome, (ii) conserved genome size (between 2.8 and 3.2 Mb), and (iii) a highly syntenic genome. Limited gene loss in Listeria did include loss of virulence associated genes, likely associated with multiple transitions to a saprotrophic lifestyle. The genus Listeria thus provides an example of a group of bacteria that appears to evolve through a loss of virulence rather than acquisition of virulence characteristics. While Listeria includes a number of species-like clades, many of these putative species include clades or strains with atypical virulence associated characteristics. This information will allow for the development of genetic and genomic criteria for pathogenic strains, including development of assays that specifically detect pathogenic Listeria strains.

Rapid Whole-Genome Sequencing for Surveillance of Salmonella enterica Serovar Enteritidis

For Salmonella enterica serovar Enteritidis, 85% of isolates can be classified into 5 pulsed-field gel electrophoresis (PFGE) types. However, PFGE has limited discriminatory power for outbreak detection. Although whole-genome sequencing has been found to improve discrimination of outbreak clusters, whether this procedure can be used in real-time in a public health laboratory is not known. Therefore, we conducted a retrospective and prospective analysis. The retrospective study investigated isolates from 1 confirmed outbreak. Additional cases could be attributed to the outbreak strain on the basis of whole-genome data. The prospective study included 58 isolates obtained in 2012, including isolates from 1 epidemiologically defined outbreak. Whole-genome sequencing identified additional isolates that could be attributed to the outbreak, but which differed from the outbreak-associated PFGE type. Additional putative outbreak clusters were detected in the retrospective and prospective analyses. This study demonstrates the practicality of implementing this approach for outbreak surveillance in a state public health laboratory.

Identification and Characterization of Psychrotolerant Sporeformers Associated with Fluid Milk Production and Processing

ABSTRACT Psychrotolerant spore-forming bacteria represent a major challenge to the goal of extending the shelf life of pasteurized dairy products. The objective of this study was to identify prominent phylogenetic groups of dairy-associated aerobic sporeformers and to characterize representative isolates for phenotypes relevant to growth in milk. Analysis of sequence data for a 632-nucleotide fragment of rpoB showed that 1,288 dairy-associated isolates (obtained from raw and pasteurized milk and from dairy farm environments) clustered into two major divisions representing (i) the genus Paenibacillus (737 isolates, including the species Paenibacillus odorifer, Paenibacillus graminis, and Paenibacillus amylolyticus sensu lato) and (ii) Bacillus (n = 467) (e.g., Bacillus licheniformis sensu lato, Bacillus pumilus, Bacillus weihenstephanensis) and genera formerly classified as Bacillus (n = 84) (e.g., Viridibacillus spp.). When isolates representing the most common rpoB allelic types (ATs) were tested for growth in skim milk broth at 6°C, 6/9 Paenibacillus isolates, but only 2/8 isolates representing Bacillus subtypes, grew >5 log CFU/ml over 21 days. In addition, 38/40 Paenibacillus isolates but only 3/47 Bacillus isolates tested were positive for β-galactosidase activity (including some isolates representing Bacillus licheniformis sensu lato, a common dairy-associated clade). Our study confirms that Paenibacillus spp. are the predominant psychrotolerant sporeformers in fluid milk and provides 16S rRNA gene and rpoB subtype data and phenotypic characteristics facilitating the identification of aerobic spore-forming spoilage organisms of concern. These data will be critical for the development of detection methods and control strategies that will reduce the introduction of psychrotolerant sporeformers and extend the shelf life of dairy products.

Lineage specific recombination rates and microevolution in Listeria monocytogenes

BackgroundThe bacterium Listeria monocytogenes is a saprotroph as well as an opportunistic human foodborne pathogen, which has previously been shown to consist of at least two widespread lineages (termed lineages I and II) and an uncommon lineage (lineage III). While some L. monocytogenes strains show evidence for considerable diversification by homologous recombination, our understanding of the contribution of recombination to L. monocytogenes evolution is still limited. We therefore used STRUCTURE and ClonalFrame, two programs that model the effect of recombination, to make inferences about the population structure and different aspects of the recombination process in L. monocytogenes. Analyses were performed using sequences for seven loci (including the house-keeping genes gap, prs, purM and ribC, the stress response gene sigB, and the virulence genes actA and inlA) for 195 L. monocytogenes isolates.ResultsSequence analyses with ClonalFrame and the Sawyers test showed that recombination is more prevalent in lineage II than lineage I and is most frequent in two house-keeping genes (ribC and purM) and the two virulence genes (actA and inlA). The relative occurrence of recombination versus point mutation is about six times higher in lineage II than in lineage I, which causes a higher genetic variability in lineage II. Unlike lineage I, lineage II represents a genetically heterogeneous population with a relatively high proportion (30% average) of genetic material imported from external sources. Phylograms, constructed with correcting for recombination, as well as Tajimas D data suggest that both lineages I and II have suffered a population bottleneck.ConclusionOur study shows that evolutionary lineages within a single bacterial species can differ considerably in the relative contributions of recombination to genetic diversification. Accounting for recombination in phylogenetic studies is critical, and new evolutionary models that account for the possibility of changes in the rate of recombination would be required. While previous studies suggested that only L. monocytogenes lineage I has experienced a recent bottleneck, our analyses clearly show that lineage II experienced a bottleneck at about the same time, which was subsequently obscured by abundant homologous recombination after the lineage II bottleneck. While lineage I and lineage II should be considered separate species from an evolutionary viewpoint, maintaining single species name may be warranted since both lineages cause the same type of human disease.

A Whole-Genome Single Nucleotide Polymorphism-Based Approach To Trace and Identify Outbreaks Linked to a Common Salmonella enterica subsp. enterica Serovar Montevideo Pulsed-Field Gel Electrophoresis Type

ABSTRACT In this study, we report a whole-genome single nucleotide polymorphism (SNP)-based evolutionary approach to study the epidemiology of a multistate outbreak of Salmonella enterica subsp. enterica serovar Montevideo. This outbreak included 272 cases that occurred in 44 states between July 2009 and April 2010. A case-control study linked the consumption of salami made with contaminated black and red pepper to the outbreak. We sequenced, on the SOLiD System, 47 isolates with XbaI PFGE pattern JIXX01.0011, a common pulsed-field gel electrophoresis (PFGE) pattern associated with isolates from the outbreak. These isolates represented 20 isolates collected from human sources during the period of the outbreak and 27 control isolates collected from human, food, animal, and environmental sources before the outbreak. Based on 253 high-confidence SNPs, we were able to reconstruct a tip-dated molecular clock phylogeny of the isolates and to assign four human isolates to the actual outbreak. We developed an SNP typing assay to rapidly discriminate between outbreak-related cases and non-outbreak-related cases and tested this assay on an extended panel of 112 isolates. These results suggest that only a very small percentage of the human isolates with the outbreak PFGE pattern and obtained during the outbreak period could be attributed to the actual pepper-related outbreak (20%), while the majority (80%) of the putative cases represented background cases. This study demonstrates that next-generation-based SNP typing provides the resolution and accuracy needed for outbreak investigations of food-borne pathogens that cannot be distinguished by currently used subtyping methods.

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

BackgroundDivergence 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.ResultsMulti-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.ConclusionsS. 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.

Listeria floridensis sp. nov., Listeria aquatica sp. nov., Listeria cornellensis sp. nov., Listeria riparia sp. nov. and Listeria grandensis sp. nov., from agricultural and natural environments.

Sampling of agricultural and natural environments in two US states (Colorado and Florida) yielded 18 Listeria-like isolates that could not be assigned to previously described species using traditional methods. Using whole-genome sequencing and traditional phenotypic methods, we identified five novel species, each with a genome-wide average BLAST nucleotide identity (ANIb) of less than 85% to currently described species. Phylogenetic analysis based on 16S rRNA gene sequences and amino acid sequences of 31 conserved loci showed the existence of four well-supported clades within the genus Listeria; (i) a clade representing Listeria monocytogenes, L. marthii, L. innocua, L. welshimeri, L. seeligeri and L. ivanovii, which we refer to as Listeria sensu stricto, (ii) a clade consisting of Listeria fleischmannii and two newly described species, Listeria aquatica sp. nov. (type strain FSL S10-1188(T) = DSM 26686(T) = LMG 28120(T) = BEI NR-42633(T)) and Listeria floridensis sp. nov. (type strain FSL S10-1187(T) = DSM 26687(T) = LMG 28121(T) = BEI NR-42632(T)), (iii) a clade consisting of Listeria rocourtiae, L. weihenstephanensis and three novel species, Listeria cornellensis sp. nov. (type strain TTU A1-0210(T) = FSL F6-0969(T) = DSM 26689(T) = LMG 28123(T) = BEI NR-42630(T)), Listeria grandensis sp. nov. (type strain TTU A1-0212(T) = FSL F6-0971(T) = DSM 26688(T) = LMG 28122(T) = BEI NR-42631(T)) and Listeria riparia sp. nov. (type strain FSL S10-1204(T) = DSM 26685(T) = LMG 28119(T) = BEI NR- 42634(T)) and (iv) a clade containing Listeria grayi. Genomic and phenotypic data suggest that the novel species are non-pathogenic.

A Population Genetics-Based and Phylogenetic Approach to Understanding the Evolution of Virulence in the Genus Listeria

ABSTRACT The genus Listeria includes (i) the opportunistic pathogens L. monocytogenes and L. ivanovii, (ii) the saprotrophs L. innocua, L. marthii, and L. welshimeri, and (iii) L. seeligeri, an apparent saprotroph that nevertheless typically contains the prfA virulence gene cluster. A novel 10-loci multilocus sequence typing scheme was developed and used to characterize 67 isolates representing six Listeria spp. (excluding L. grayi) in order to (i) provide an improved understanding of the phylogeny and evolution of the genus Listeria and (ii) use Listeria as a model to study the evolution of pathogenicity in opportunistic environmental pathogens. Phylogenetic analyses identified six well-supported Listeria species that group into two main subdivisions, with each subdivision containing strains with and without the prfA virulence gene cluster. Stochastic character mapping and phylogenetic analysis of hly, a gene in the prfA cluster, suggest that the common ancestor of the genus Listeria contained the prfA virulence gene cluster and that this cluster was lost at least five times during the evolution of Listeria, yielding multiple distinct saprotrophic clades. L. welshimeri, which appears to represent the most ancient clade that arose from an ancestor with a prfA cluster deletion, shows a considerably lower average sequence divergence than other Listeria species, suggesting a population bottleneck and a putatively different ecology than other saprotrophic Listeria species. Overall, our data suggest that, for some pathogens, loss of virulence genes may represent a selective advantage, possibly by facilitating adaptation to a specific ecological niche.