MLVA as a tool for public health surveillance of human Salmonella Typhimurium: prospective study in Belgium and evaluation of MLVA loci stability

Isolation and molecular characterization of the Salmonella Typhimurium orphan phage Arash

The current threat of multidrug resistant strains necessitates development of alternatives to antibiotics such as bacteriophages. This study describes the isolation and characterization of a novel Salmonella Typhimurium phage 'Arash' from hospital wastewater in Leuven, Belgium. Arash has a myovirus morphology with a 95 nm capsid and a 140 nm tail. The host range of Arash is restricted to its isolation host. Approximately 86% of the phage particles are adsorbed to a host cell within 10 min. Arash has latent period of 65 min and burst size of 425 PFU/cell. Arash has a dsDNA genome of 180,819 bp with GC content of 53.02% with no similarities to any characterized phages, suggesting Arash as a novel species in the novel 'Arashvirus' genus. Arash carries no apparent lysogeny-, antibiotic resistance- nor virulence-related genes. Proteome analysis revealed 116 proteins as part of the mature phage particles of which 27 could be assigned a function. Therefore, the present findings shed light on the morphological, microbiological and genomic characteristics of Arash and suggest its potential application as therapeutic and/or biocontrol agent.

Application of a new multi-locus variable number tandem repeat analysis (MLVA) scheme for the seasonal investigation of Cryptosporidium parvum cases in Wales and the northwest of England, spring 2022

The protozoan Cryptosporidium parvum is an important cause of gastroenteritis in humans and livestock, and cryptosporidiosis outbreaks are common. However, a multi-locus genotyping scheme is not widely adopted. We describe the further development and application of a seven-locus multi-locus variable number of tandem repeats analysis (MLVA) scheme. From 28th March to 31st July 2022, confirmed C. parvum stools (n = 213) from cryptosporidiosis patients (cases) in Wales (n = 95) and the north west of England (n = 118) were tested by MLVA. Typability (defined as alleles identified at all seven loci in a sample) was 81.2% and discriminatory power estimated by Hunter Gaston Discriminatory Index was 0.99. A MLVA profile was constructed from the alleles, expressed in chromosomal order. Profiles were defined as simple (single allele at each locus) or mixed (more than one allele at any locus). A total of 161 MLVA profiles were identified; 13 were mixed, an additional 38 simple profiles contained null records, and 110 were complete simple profiles. A minimum spanning tree was constructed of simple MLVA profiles and those identical at all seven loci defined genetic clusters of cases (here, null records were considered as an allele); 77 cases formed 25 clusters, ranging from two to nine (mode = two) cases. The largest cluster, following epidemiological investigation, signalled a newly-identified outbreak. Two other cases with mixed profiles that contained the outbreak alleles were included in the outbreak investigation. In another epidemiologically-identified outbreak of six initial cases, MLVA detected two additional cases. In a third, small outbreak of three cases, identical MLVA profiles strengthened the microbiological evidence. Review of the performance characteristics of the individual loci and of the seven-locus scheme suggested that two loci might be candidates for review, but a larger dataset over a wider geographical area and longer timeframe will help inform decision-making about the scheme by user laboratories and stakeholders (such as public health agencies). This MLVA scheme is straightforward in use, fast and cheap compared to sequence-based methods, identifies mixed infections, provides an important tool for C. parvum surveillance, and can enhance outbreak investigations and public health action.

Genetic diversity and delineation of Salmonella Agona outbreak strains by next generation sequencing, Bavaria, Germany, 1993 to 2018

Background

In 2017, a food-borne Salmonella Agona outbreak caused by infant milk products from a French supplier occurred in Europe. Simultaneously, S. Agona was detected in animal feed samples in Bavaria.

Aim

Using next generation sequencing (NGS) and three data analysis methods, this study’s objectives were to verify clonality of the Bavarian feed strains, rule out their connection to the outbreak, explore the genetic diversity of Bavarian S. Agona isolates from 1993 to 2018 and compare the analysis approaches employed, for practicality and ability to delineate outbreaks caused by the genetically monomorphic Agona serovar.

Methods

In this observational retrospective study, three 2017 Bavarian feed isolates were compared to a French outbreak isolate and 48 S. Agona isolates from our strain collections. The later included human, food, feed, veterinary and environmental isolates, of which 28 were epidemiologically outbreak related. All isolates were subjected to NGS and analysed by: (i) a publicly available species-specific core genome multilocus sequence typing (cgMLST) scheme, (ii) single nucleotide polymorphism phylogeny and (iii) an in-house serovar-specific cgMLST scheme. Using additional international S. Agona outbreak NGS data, the cluster resolution capacity of the two cgMLST schemes was assessed.

Results

We could prove clonality of the feed isolates and exclude their relation to the French outbreak. All approaches confirmed former Bavarian epidemiological clusters.

Conclusion

Even for S. Agona, species-level cgMLST can produce reasonable resolution, being standardisable by public health laboratories. For single samples or homogeneous sample sets, higher resolution by serovar-specific cgMLST or SNP genotyping can facilitate outbreak investigations.

High genetic similarity between non-typhoidal Salmonella isolated from paired blood and stool samples of children in the Democratic Republic of the Congo

BACKGROUND

Non-typhoidal Salmonella (NTS) serotypes Typhimurium and Enteritidis are a major cause of bloodstream infections in children in sub-Saharan Africa but their reservoir is unknown. We compared pairs of NTS blood and stool isolates (with the same NTS serotype recovered in the same patient) for genetic similarity.

METHODS

Between November 2013 and April 2017, hospital-admitted children (29 days to 14 years) with culture-confirmed NTS bloodstream infections were enrolled in a cross-sectional study at Kisantu Hospital, DR Congo. Stool cultures for Salmonella were performed on a subset of enrolled children, as well as on a control group of non-febrile hospital-admitted children. Pairs of blood and stool NTS isolates were assessed for genetic similarity by multiple-locus variable-number of tandem repeats (MLVA) and genomics analysis.

RESULTS

A total of 299 children with NTS grown from blood cultures (Typhimurium 68.6%, Enteritidis 30.4%, other NTS 1.0%) had a stool sample processed; in 105 (35.1%) of them NTS was detected (Typhimurium 70.5%, Enteritidis 25.7%, other NTS 3.8%). A total of 87/105 (82.9%) pairs of blood and stool NTS isolates were observed (representing 29.1% of the 299 children). Among 1598 controls, the proportion of NTS stool excretion was 2.1% (p < 0.0001). MLVA types among paired isolates were identical in 82/87 (94.3%) pairs (27.4% of the 299 children; 61/66 (92.4%) in Typhimurium and 21/21 (100%) in Enteritidis pairs). Genomics analysis confirmed high genetic similarity within 41/43 (95.3%) pairs, showing a median SNP difference of 1 (range 0-77) and 1 (range 0-4) for Typhimurium and Enteritidis pairs respectively. Typhimurium and Enteritidis isolates belonged to sequence types ST313 lineage II and ST11 respectively.

CONCLUSION

Nearly 30% of children with NTS bloodstream infection showed stool excretion of an NTS isolate with high genetic similarity, adding to the evidence of humans as a potential reservoir for NTS.

Multidrug-Resistant Salmonella enterica Serovar Rissen Clusters Detected in Azores Archipelago, Portugal

Gastrointestinal infections caused by nontyphoidal Salmonella (NTS) remain one of the main causes of foodborne illness worldwide. Within the multiple existing Salmonella enterica serovars, the serovar Rissen is rarely reported, particularly as a cause of human salmonellosis. Between 2015 and 2017, the Portuguese National Reference Laboratory of Gastrointestinal Infections observed an increase in the number of clinical cases caused by multidrug-resistant (MDR) S. enterica serovar Rissen, particularly from the Azores archipelago. In the present study, we analyzed by whole genome sequencing (WGS) all clinical, animal, food, and environmental isolates received up to 2017 in the Portuguese Reference Laboratories. As such, through a wgMLST-based gene-by-gene analysis, we aimed to identify potential epidemiological clusters linking clinical and samples from multiple sources, while gaining insight into the genetic diversity of S. enterica serovar Rissen. We also investigated the genetic basis driving the observed multidrug resistance. By integrating 60 novel genomes with all publicly available serovar Rissen genomes, we observed a low degree of genetic diversity within this serovar. Nevertheless, the majority of Portuguese isolates showed high degree of genetic relatedness and a potential link to pork production. An in-depth analysis of these isolates revealed the existence of two major clusters from the Azores archipelago composed of MDR isolates, most of which were resistant to at least five antimicrobials. Considering the well-known spread of MDR between gastrointestinal bacteria, the identification of MDR circulating clones should constitute an alert to public health authorities. Finally, this study constitutes the starting point for the implementation of the "One Health" approach for Salmonella surveillance in Portugal.

In vivo passage of Salmonella Typhimurium results in minor mutations in the bacterial genome and increases in vitro invasiveness

Eggs and raw or undercooked egg-containing food items are frequently identified as the bacterial source during epidemiolocal investigation of Salmonella outbreaks. Multi-locus variable number of tandem repeats analysis (MLVA) is a widely used Salmonella typing method enabling the study of diversity within populations of the same serotype. In vivo passage, however, has been linked with changes in MLVA type and more broadly the Salmonella genome. We sought to investigate whether in vivo passage through layer hens had an effect on MLVA type as well as the bacterial genome and whether any mutations affected bacterial virulence. Layer hens were infected with either Salmonella Typhimurium DT9 (03-24-11-11-523) as part of a single infection or were co-infected with an equal amount of Salmonella Mbandaka. Salmonella shedding in both single and co-infected birds was variable over the course of the 16-week experiment. Salmonella Typhimurium and Salmonella Mbandaka were identified in feces of co-infected birds. Salmonella colonies isolated from fecal samples were subtyped using MLVA. A single change in SSTR-6 was observed in Salmonella Typhimurium strains isolated from co-infected birds. Isolates of Salmonella Typhimurium of both the parent (03-24-11-11-523) and modified (03-24-12-11-523) MLVA type were sequenced and compared with the genome of the parent strain. Sequence analysis revealed that in vivo passaging resulted in minor mutation events. Passaged isolates exhibited significantly higher invasiveness in cultured human intestinal epithelial cells than the parent strain. The microevolution observed in this study suggests that changes in MLVA may arise more commonly and may have clinical significance.

Rise and fall of outbreak-specific clone inside endemic pulsotype of Salmonella 4,[5],12:i:-; insights from high-resolution molecular surveillance in Emilia-Romagna, Italy, 2012 to 2015

Background and aimEpidemiology of human non-typhoid salmonellosis is characterised by recurrent emergence of new clones of the pathogen over time. Some clonal lines of Salmonella have shaped epidemiology of the disease at global level, as happened for serotype Enteritidis or, more recently, for Salmonella 4,[5],12:i:-, a monophasic variant of serotype Typhimurium. The same clonal behaviour is recognisable at sub-serotype level where single outbreaks or more generalised epidemics are attributable to defined clones. The aim of this study was to understand the dynamics of a clone of Salmonella 4,[5],12:i:- over a 3-year period (2012-15) in a province of Northern Italy where the clone caused a large outbreak in 2013. Furthermore, the role of candidate outbreak sources was investigated and the accuracy of multilocus variable-number tandem repeat analysis (MLVA) was evaluated. Methods: we retrospectively investigated the outbreak through whole genome sequencing (WGS) and further monitored the outbreak clone for 2 years after its conclusion. Results: The study showed the transient nature of the clone in the population, possibly as a consequence of its occasional expansion in a food-processing facility. We demonstrated that important weaknesses characterise conventional typing methods applied to clonal pathogens such as Salmonella 4,[5],12:i:-, namely lack of accuracy for MLVA and inadequate resolution power for PFGE to be reliably used for clone tracking. Conclusions: The study provided evidence for the remarkable prevention potential of whole genome sequencing used as a routine tool in systems that integrate human, food and animal surveillance.

Assessment and Comparison of Molecular Subtyping and Characterization Methods for Salmonella

The food industry is facing a major transition regarding methods for confirmation, characterization, and subtyping of Salmonella. Whole-genome sequencing (WGS) is rapidly becoming both the method of choice and the gold standard for Salmonella subtyping; however, routine use of WGS by the food industry is often not feasible due to cost constraints or the need for rapid results. To facilitate selection of subtyping methods by the food industry, we present: (i) a comparison between classical serotyping and selected widely used molecular-based subtyping methods including pulsed-field gel electrophoresis, multilocus sequence typing, and WGS (including WGS-based serovar prediction) and (ii) a scoring system to evaluate and compare Salmonella subtyping assays. This literature-based assessment supports the superior discriminatory power of WGS for source tracking and root cause elimination in food safety incident; however, circumstances in which use of other subtyping methods may be warranted were also identified. This review provides practical guidance for the food industry and presents a starting point for further comparative evaluation of Salmonella characterization and subtyping methods.

Traditional Salmonella Typhimurium typing tools (phage typing and MLVA) are sufficient to resolve well-defined outbreak events only

Between 1991 and 2014 the per capita notification rate of salmonellosis in Australia increased from 31.9 to 69.7 cases per 100,000 people. Salmonella Typhimurium accounted for nearly half the human cases until the end of 2014. In this study, we used cluster analysis tools to compare S. Typhimurium isolates from a chicken-meat study with those reported to the National Enteric Pathogen Surveillance System (NEPSS) from the coincident human and non-human populations. There was limited phage type diversity within all populations and a lack of specificity of MLVA profiling within phage types. The chicken-meat study isolates were not significantly clustered with the human cases and at least 7 non-human sources, based on typing profiles (PT/MLVA combination), could be implicated as a source of human cases during the same period. In the absence of a strong surveillance system representative of all putative sources, MLVA and phage typing alone or in combination are insufficient to identify the source of human cases.

Salmonella spp. transmission in a vertically integrated poultry operation: Clustering and diversity analysis using phenotyping (serotyping, phage typing) and genotyping (MLVA)

The transmission of Salmonella enterica within a vertically integrated poultry operation was investigated longitudinally over an 18-month period (2013–2014). Thirty six percent of all samples collected (1503 of 4219) were positive for salmonellae with seven Salmonella enterica subsp. enterica serovars, and one Salmonella enterica subsp. salamae serovar detected. Both Salmonella enterica subsp. enterica serovars Infantis and Typhimurium were detected in all locations sampled. Salmonella Typhimurium was the most frequently detected serovar (63% of serotyped samples) with 8 phage types (PT) and 41 multiple-locus variable-number tandem-repeats analysis (MLVA) profiles identified. The most frequently identified phage types were PT135a and DT135. A total of 62 PT/MLVA combinations were observed. MLVA profiles 03-14-10-09-525 and 03-15-11-11-525 were the most frequently identified and 83% of the isolates shared at least one MLVA profile with an isolate from another phage type. The use of phage typing and MLVA profiling, on their own or in combination, were insufficient to understand the complexity of the epidemiological relationships between locations within this production system. Despite the high level of apparent diversity, cluster analysis was unable to differentiate the transmission pathways of all S. Typhimurium variants detected within the integrated enterprise. Using additional epidemiological information, the parent breeder rearing site was identified as the most likely point of introduction of two S. Typhimurium isolates into the production system with subsequent dissemination to the broiler flocks via the hatchery. This complexity is unable to be resolved in the absence of intensive sampling programs at all generations of the production system.

Comparison of SNP-based subtyping workflows for bacterial isolates using WGS data, applied to Salmonella enterica serotype Typhimurium and serotype 1,4,[5],12:i:

Whole genome sequencing represents a promising new technology for subtyping of bacterial pathogens. Besides the technological advances which have pushed the approach forward, the last years have been marked by considerable evolution of the whole genome sequencing data analysis methods. Prior to application of the technology as a routine epidemiological typing tool, however, reliable and efficient data analysis strategies need to be identified among the wide variety of the emerged methodologies. In this work, we have compared three existing SNP-based subtyping workflows using a benchmark dataset of 32 Salmonella enterica subsp. enterica serovar Typhimurium and serovar 1,4,[5],12:i:- isolates including five isolates from a confirmed outbreak and three isolates obtained from the same patient at different time points. The analysis was carried out using the original (high-coverage) and a down-sampled (low-coverage) datasets and two different reference genomes. All three tested workflows, namely CSI Phylogeny-based workflow, CFSAN-based workflow and PHEnix-based workflow, were able to correctly group the confirmed outbreak isolates and isolates from the same patient with all combinations of reference genomes and datasets. However, the workflows differed strongly with respect to the SNP distances between isolates and sensitivity towards sequencing coverage, which could be linked to the specific data analysis strategies used therein. To demonstrate the effect of particular data analysis steps, several modifications of the existing workflows were also tested. This allowed us to propose data analysis schemes most suitable for routine SNP-based subtyping applied to S. Typhimurium and S. 1,4,[5],12:i:-. Results presented in this study illustrate the importance of using correct data analysis strategies and to define benchmark and fine-tune parameters applied within routine data analysis pipelines to obtain optimal results.

Use of multiple-locus variable-number of tandem repeats analysis (MLVA) to investigate genetic diversity of Salmonella enterica subsp. enterica serovar Typhimurium isolates from human, food, and veterinary sources

Salmonella enterica subspecies enterica serovar Typhimurium is the most common zoonotic pathogen in Bulgaria. To allow efficient outbreak investigations and surveillance in the food chain, accurate and discriminatory methods for typing are needed. This study evaluated the use of multiple-locus variable-number of tandem repeats analysis (MLVA) and compared results with antimicrobial resistance (AMR) determinations for 100 S. Typhimurium strains isolated in Bulgaria during 2008-2012 (50 veterinary/food and 50 human isolates). Results showed that isolates were divided into 80 and 34 groups using MLVA and AMR, respectively. Simpson's index of diversity was determined to 0.994 ± 0.003 and 0.945 ± 0.012. The most frequently encountered MLVA profiles were 3-11-9-NA-211 (n = 5); 3-12-9-NA-211 (n = 3); 3-12-11-21-311 (n = 3); 3-17-10-NA-311 (n = 3); 2-20-9-7-212 (n = 3); and 2-23-NA-NA-111 (n = 3). No clustering of isolates related to susceptibility/resistance to antimicrobials, source of isolation, or year of isolation was observed. Some MLVA types were found in both human and veterinary/food isolates, indicating a possible route of transmission. A majority (83%) of the isolates were found to be resistant against at least one antimicrobial and 44% against ≥4 antimicrobials. Further studies are needed to verify MLVA usefulness over a longer period of time and with more isolates, including outbreak strains.

Phenotypic and Genotypic Eligible Methods for Salmonella Typhimurium Source Tracking

Salmonellosis is one of the most common causes of foodborne infection and a leading cause of human gastroenteritis. Throughout the last decade, Salmonella enterica serotype Typhimurium (ST) has shown an increase report with the simultaneous emergence of multidrug-resistant isolates, as phage type DT104. Therefore, to successfully control this microorganism, it is important to attribute salmonellosis to the exact source. Studies of Salmonella source attribution have been performed to determine the main food/food-production animals involved, toward which, control efforts should be correctly directed. Hence, the election of a ST subtyping method depends on the particular problem that efforts must be directed, the resources and the data available. Generally, before choosing a molecular subtyping, phenotyping approaches such as serotyping, phage typing, and antimicrobial resistance profiling are implemented as a screening of an investigation, and the results are computed using frequency-matching models (i.e., Dutch, Hald and Asymmetric Island models). Actually, due to the advancement of molecular tools as PFGE, MLVA, MLST, CRISPR, and WGS more precise results have been obtained, but even with these technologies, there are still gaps to be elucidated. To address this issue, an important question needs to be answered: what are the currently suitable subtyping methods to source attribute ST. This review presents the most frequently applied subtyping methods used to characterize ST, analyses the major available microbial subtyping attribution models and ponders the use of conventional phenotyping methods, as well as, the most applied genotypic tools in the context of their potential applicability to investigates ST source tracking.

Salmonella Typhimurium DT193 and DT99 are present in great and blue tits in Flanders, Belgium

Endemic infections with the common avian pathogen Salmonella enterica subspecies enterica serovar Typhimurium (Salmonella Typhimurium) may incur a significant cost on the host population. In this study, we determined the potential of endemic Salmonella infections to reduce the reproductive success of blue (Cyanistes caeruleus) and great (Parus major) tits by correlating eggshell infection with reproductive parameters. The fifth egg of each clutch was collected from nest boxes in 19 deciduous forest fragments. Out of the 101 sampled eggs, 7 Salmonella Typhimurium isolates were recovered. The low bacterial prevalence was reflected by a similarly low serological prevalence in the fledglings. In this study with a relatively small sample size, presence of Salmonella did not affect reproductive parameters (egg volume, clutch size, number of nestlings and number of fledglings), nor the health status of the fledglings. However, in order to clarify the impact on health and reproduction a larger number of samples have to be analyzed. Phage typing showed that the isolates belonged to the definitive phage types (DT) 193 and 99, and multi-locus variable number tandem repeat analysis (MLVA) demonstrated a high similarity among the tit isolates, but distinction to human isolates. These findings suggest the presence of passerine-adapted Salmonella strains in free-ranging tit populations with host pathogen co-existence.

Enrichment, Amplification, and Sequence-Based Typing of Salmonella enterica and Other Foodborne Pathogens

Detection of Salmonella enterica in foods typically involves microbiological enrichment, molecular-based assay, and subsequent isolation and identification of a pure culture. This is ideally followed by strain typing, which provides information critical to the investigation of outbreaks and the attribution of their sources. Pulsed-field gel electrophoresis is the "gold standard" for S. enterica strain typing, but its limitations have encouraged the search for alternative methods, including whole genome sequencing. Both methods typically require a pure culture, which adds to the cost and turnaround time. A more rapid and cost-effective method with sufficient discriminatory power would benefit food industries, regulatory agencies, and public health laboratories. To address this need, a novel enrichment, amplification, and sequence-based typing (EAST) approach was developed involving (i) overnight enrichment and total DNA preparation, (ii) amplification of polymorphic tandem repeat-containing loci with electrophoretic detection, and (iii) DNA sequencing and bioinformatic analysis to identify related strains. EAST requires 3 days or less and provides a strain resolution that exceeds serotyping and is comparable to pulsed-field gel electrophoresis. Evaluation with spiked ground turkey demonstrated its sensitivity (with a starting inoculum of ≤1 CFU/g) and specificity (with unique or nearly unique alleles relative to databases of >1,000 strains). In tests with unspiked retail chicken parts, 3 of 11 samples yielded S. enterica -specific PCR products. Sequence analysis of three distinct typing targets (SeMT1, SeCRISPR1, and SeCRISPR2) revealed consistent similarities to specific serotype Schwarzengrund, Montevideo, and Typhimurium strains. EAST provides a time-saving and cost-effective approach for detecting and typing foodborne S. enterica , and postenrichment steps can be commercially outsourced to facilitate its implementation. Initial studies with Listeria monocytogenes and Shiga toxigenic Escherichia coli suggest that EAST can be extended to these foodborne pathogens.

Evolution of Variable Number Tandem Repeats and Its Relationship with Genomic Diversity in Salmonella Typhimurium

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing human infections in Australia and many other countries. A total of 12,112 S. Typhimurium isolates from New South Wales were analyzed by multi-locus variable number of tandem repeat (VNTR) analysis (MLVA) using five VNTRs from 2007 to 2014. We found that mid ranges of repeat units of 8–14 in VNTR locus STTR5, 6–13 in STTR6, and 9–12 in STTR10 were always predominant in the population (>50%). In vitro passaging experiments using MLVA type carrying extreme length alleles found that the majority of long length alleles mutated to short ones and short length alleles mutated to longer ones. Both data suggest directional mutability of VNTRs toward mid-range repeats. Sequencing of 28 isolates from a newly emerged MLVA type and its five single locus variants revealed that single nucleotide variation between isolates with up to two MLVA differences ranged from 0 to 12 single nucleotide polymorphisms (SNPs). However, there was no relationship between SNP and VNTR differences. A population genetic model of the joint distribution of VNTRs and SNPs variations was used to estimate the mutation rates of the two markers, yielding a ratio of 1 VNTR change to 6.9 SNP changes. When only one VNTR repeat difference was considered, the majority of pairwise SNP difference between isolates were 4 SNPs or fewer. Based on this observation and our previous findings of SNP differences of outbreak isolates, we suggest that investigation of S. Typhimurium community outbreaks should include cases of 1 repeat difference to increase sensitivity. This study offers new insights into the short-term VNTR evolution of S. Typhimurium and its application for epidemiological typing.

Virulence-associated genes, antimicrobial resistance and molecular typing of Salmonella Typhimurium strains isolated from swine from 2000 to 2012 in Brazil

AIMS

The aims of this study were to assess the pathogenic potential, antimicrobial resistance and genotypic diversity of Salmonella Typhimurium strains isolated in Brazil from swine (22) and the surrounding swine environment (5) from 2000 to 2012 and compare them to the profiles of 43 human strains isolated from 1983 to 2010, which had been previously studied.

METHODS AND RESULTS

The presence of 12 SPI-1, SPI-2 and plasmid genes was assessed by PCR, the antimicrobial susceptibility to 13 antimicrobials was determined by the disc diffusion assay and genotyping was performed using pulsed-field gel electrophoresis (PFGE), multiple-locus variable-number of tandem repeats analysis (MLVA) and ERIC-PCR. More than 77·8% of the swine strains carried 10 or more of the virulence markers. Ten (37%) strains isolated from swine were multi-drug resistant (MDR). All the molecular typing techniques grouped the strains in two main clusters. Some strains isolated from swine and humans were allocated together in the PFGE-B2, MLVA-A1, MLVA-B and ERIC-A1 clusters.

CONCLUSIONS

The genotyping results suggest that some strains isolated from swine and humans may descend from a common subtype and may indicate a possible risk of MDR S. Typhimurium with high frequency of virulence genes isolated from swine to contaminate humans in Brazil.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study provided new information about the pathogenic potential, antimicrobial resistance and genotypic diversity of S. Typhimurium isolates from swine origin in Brazil, the fourth largest producer of pigs worldwide.

Application of Molecular Typing Results in Source Attribution Models: The Case of Multiple Locus Variable Number Tandem Repeat Analysis (MLVA) of Salmonella Isolates Obtained from Integrated Surveillance in Denmark

Salmonella is an important cause of bacterial foodborne infections in Denmark. To identify the main animal-food sources of human salmonellosis, risk managers have relied on a routine application of a microbial subtyping-based source attribution model since 1995. In 2013, multiple locus variable number tandem repeat analysis (MLVA) substituted phage typing as the subtyping method for surveillance of S. Enteritidis and S. Typhimurium isolated from animals, food, and humans in Denmark. The purpose of this study was to develop a modeling approach applying a combination of serovars, MLVA types, and antibiotic resistance profiles for the Salmonella source attribution, and assess the utility of the results for the food safety decisionmakers. Full and simplified MLVA schemes from surveillance data were tested, and model fit and consistency of results were assessed using statistical measures. We conclude that loci schemes STTR5/STTR10/STTR3 for S. Typhimurium and SE9/SE5/SE2/SE1/SE3 for S. Enteritidis can be used in microbial subtyping-based source attribution models. Based on the results, we discuss that an adjustment of the discriminatory level of the subtyping method applied often will be required to fit the purpose of the study and the available data. The issues discussed are also considered highly relevant when applying, e.g., extended multi-locus sequence typing or next-generation sequencing techniques.

The evaluation and application of multilocus variable number tandem repeat analysis (MLVA) for the molecular epidemiological study of Salmonella enterica subsp. enterica serovar Enteritidis infection

Background

Salmonella enterica subsp. enterica serovar Enteritidis (S. Enteritidis) is one of the most prevalent Salmonella serotypes that cause gastroenteritis worldwide and the most prevalent serotype causing Salmonella infections in China. A rapid molecular typing method with high throughput and good epidemiological discrimination is urgently needed for detecting the outbreaks and finding the source for effective control of S. Enteritidis infections.

Methods

In this study, 194 strains which included 47 from six outbreaks that were well-characterized epidemiologically were analyzed with pulse field gel electrophoresis (PFGE) and multilocus variable number tandem repeat analysis (MLVA). Seven VNTR loci published by the US Center for Disease Control and Prevention (CDC) were used to evaluate and develop MLVA scheme for S. Enteritidis molecular subtyping by comparing with PFGE, and then MLVA was applied to the suspected outbreaks detection. All S. Enteritidis isolates were analyzed with MLVA to establish a MLVA database in Shenzhen, Guangdong province, China to facilitate the detection of S. Enteritidis infection clusters.

Results

There were 33 MLVA types and 29 PFGE patterns among 147 sporadic isolates. These two measures had Simpson indices of 0.7701 and 0.8043, respectively, which did not differ significantly. Epidemiological concordance was evaluated by typing 47 isolates from six epidemiologically well-characterized outbreaks and it did not differ for PFGE and MLVA. We applied the well established MLVA method to detect two S. Enteritidis foodborne outbreaks and find their sources successfully in 2014. A MLVA database of 491 S. Enteritidis strains isolated from 2004 to 2014 was established for the surveillance of clusters in the future.

Conclusions

MLVA typing of S. Enteritidis would be an effective tool for early warning and epidemiological surveillance of S. Enteritidis infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12941-016-0119-3) contains supplementary material, which is available to authorized users.

Investigation of Clostridium difficile interspecies relatedness using multilocus sequence typing, multilocus variable-number tandem-repeat analysis and antimicrobial susceptibility testing

Multilocus sequence typing (MLST), multilocus variable-number tandem-repeat analysis (MLVA) and antimicrobial susceptibility were performed on 37 animal and human C. difficile isolates belonging to 15 different PCR-ribotypes in order to investigate the relatedness of human and animal isolates and to identify possible transmission routes. MLVA identified a total of 21 different types while MLST only distinguished 12 types. Identical C. difficile strains were detected in the same animal species for PCR-ribotypes 014, 078, UCL 16U and UCL 36, irrespective of their origin or the isolation date. Non clonal strains were found among different hosts; however, a high genetic association between pig and cattle isolates belonging to PCR-ribotype 078 was revealed. MLVA also showed genetic differences that clearly distinguished human from animal strains. For a given PCR-ribotype, human and animal strains presented a similar susceptibility to the antimicrobials tested. All strains were susceptible to vancomycin, metronidazole, chloramphenicol and rifampicin, while PCR-ribotypes 078, UCL 5a, UCL 36 and UCL 103 were associated with erythromycin resistance. The data suggest a wide dissemination of clones at hospitals and breeding-farms or a contamination at the slaughterhouse, but less probability of interspecies transmission. However, further highly discriminatory genotyping methods are necessary to elucidate interspecies and zoonotic transmission of C. difficile.

Characterization of Salmonella Typhimurium isolates from domestically acquired infections in Finland by phage typing, antimicrobial susceptibility testing, PFGE and MLVA

BACKGROUND

Salmonella enterica spp. enterica serotype Typhimurium (STM) is the most common agent of domestically acquired salmonellosis in Finland. Subtyping methods which allow the characterization of STM are essential for effective laboratory-based STM surveillance and for recognition of outbreaks. This study describes the diversity of Finnish STM isolates using phage typing, antimicrobial susceptible testing, pulsed-field gel electrophoresis (PFGE) and multilocus variable-number tandem repeat analysis (MLVA), and compares the discriminatory power and the concordance of these methods.

RESULTS

A total of 375 sporadic STM isolates were analysed. The isolates were divided into 31 definite phage (DT) types, dominated by DT1 (47 % of the isolates), U277 (9 % of the isolates) and DT104 (8 % of the isolates). Of all the isolates, 62 % were susceptible to all the 12 antimicrobials tested and 11 % were multidrug resistant. Subtyping resulted in 83 different XbaI-PFGE profiles and 111 MLVA types. The three most common XbaI-PFGE profiles (STYM1, STYM7 and STYM8) and one MLVA profile with three single locus variants accounted for 56 % and 49 % of the STM isolates, respectively. The studied isolates showed a genetic similarity of more than 70 % by XbaI-PFGE. In MLVA, 71 % of the isolates lacked STTR6 and 77 % missed STTR10p loci. Nevertheless, the calculated Simpson's diversity index for XbaI-PFGE was 0.829 (95 % CI 0.792-0.865) and for MLVA 0.867 (95 % CI 0.835-0.898). However, the discriminatory power of the 5-loci MLVA varied among the phage types. The highest concordance of the results was found between XbaI-PFGE and phage typing (adjusted Wallace coefficient was 0.833 and adjusted Rand coefficient was 0.627).

CONCLUSIONS

In general, the calculated discriminatory power was higher for genotyping methods (MLVA and XbaI-PFGE) than for phenotyping methods (phage typing). Overall, comparable diversity indices were calculated for PFGE and MLVA (both DI > 0.8). However, MLVA was phage type dependent providing better discrimination of the most common phage types. Furthermore, 5-loci MLVA was a less laborious method and easier to interpret than XbaI-PFGE. Thus, the laboratory-based surveillance of the Finnish human STM infections has been conducted with a combination of phage typing, antimicrobial susceptibility testing and 5-loci MLVA since January 2014.

Defining the Core Genome of Salmonella enterica Serovar Typhimurium for Genomic Surveillance and Epidemiological Typing

Salmonella enterica serovar Typhimurium is the most common Salmonella serovar causing foodborne infections in Australia and many other countries. Twenty-one S. Typhimurium strains from Salmonella reference collection A (SARA) were analyzed using Illumina high-throughput genome sequencing. Single nucleotide polymorphisms (SNPs) in 21 SARA strains ranged from 46 to 11,916 SNPs, with an average of 1,577 SNPs per strain. Together with 47 strains selected from publicly available S. Typhimurium genomes, the S. Typhimurium core genes (STCG) were determined. The STCG consist of 3,846 genes, a set that is much larger than that of the 2,882 Salmonella core genes (SCG) found previously. The STCG together with 1,576 core intergenic regions (IGRs) were defined as the S. Typhimurium core genome. Using 93 S. Typhimurium genomes from 13 epidemiologically confirmed community outbreaks, we demonstrated that typing based on the S. Typhimurium core genome (STCG plus core IGRs) provides superior resolution and higher discriminatory power than that based on SCG for outbreak investigation and molecular epidemiology of S. Typhimurium. STCG and STCG plus core IGR typing achieved 100% separation of all outbreaks compared to that of SCG typing, which failed to separate isolates from two outbreaks from background isolates. Defining the S. Typhimurium core genome allows standardization of genes/regions to be used for high-resolution epidemiological typing and genomic surveillance of S. Typhimurium.

A novel multiple locus variable number of tandem repeat (VNTR) analysis (MLVA) method for Propionibacterium acnes

Propionibacterium acnes plays a central role in the pathogenesis of acne and is responsible for severe opportunistic infections. Numerous typing schemes have been developed that allow the identification of phylotypes, but they are often insufficient to differentiate subtypes. To better understand the genetic diversity of this species and to perform epidemiological analyses, high throughput discriminant genotyping techniques are needed. Here we describe the development of a multiple locus variable number of tandem repeats (VNTR) analysis (MLVA) method. Thirteen VNTRs were identified in the genome of P. acnes and were used to genotype a collection of clinical isolates. In addition, publically available sequencing data for 102 genomes were analyzed in silico, providing an MLVA genotype. The clustering of MLVA data was in perfect congruence with whole genome based clustering. Analysis of the clustered regularly interspaced short palindromic repeat (CRISPR) element uncovered new spacers, a supplementary source of genotypic information. The present MLVA13 scheme and associated internet database represents a first line genotyping assay to investigate large number of isolates. Particular strains may then be submitted to full genome sequencing in order to better analyze their pathogenic potential.

Ascertaining the relationship between Salmonella Typhimurium and Salmonella 4,[5],12:i:- by MLVA and inferring the sources of human salmonellosis due to the two serovars in Italy

The current picture of human salmonellosis shows Salmonella Typhimurium and S. 4,[5],12:i:- as the most common serovars in Italy. The aims of this study were to investigate the genetic relationship between these serovars, as well as to test the possibility of inferring sources of human salmonellosis due to S. Typhimurium and S. 4,[5],12:i:- by using multilocus variable-number tandem repeat analysis (MLVA) subtyping data. Single isolates from 268 human sporadic cases and 325 veterinary isolates (from pig, cattle, chicken, and turkey) collected over the period 2009-2011 were typed by MLVA, and the similarities of MLVA profiles were investigated using different analytical approaches. Results showed that isolates of S. 4,[5],12:i:- were more clonal compared to S. Typhimurium and that clones of both serovars from different non-human sources were very close to those which were responsible for human infections, suggesting that source attribution by MLVA typing should be possible. However, using the Asymmetric Island Model it was not possible to obtain a confident ranking of sources responsible for human infections based on MLVA profiles. The source assignments provided by the model could have been jeopardized by the high heterogeneity found within each source and the negligible divergence between sources as well as by the limited source data available, especially for some species.

Outbreak of Salmonella enterica serovar Typhimurium phage type DT41 in Danish poultry production

Salmonella enterica subspecies enterica serovar Typhimurium (S. Typhimurium) is one of the most prevalent serovars in Europe - where both poultry and poultry related products are common sources of human salmonellosis. Due to efficient control programs, the prevalence of S. Typhimurium in Danish poultry production is very low. Despite this, during the past decades there has been a reoccurring problem with infections with S. Typhimurium phage type DT41 in the Danish poultry production without identifying a clear source. In the end of 2013 and beginning of 2014 an increased isolation of S. Typhimurium DT41 was noted mainly in this production, but also in other samples. To investigate this is in more detail, 47 isolates from egg layers (n=5, 1 flock), broilers (n=33, 13 flocks), broiler breeding flocks and hatches (n=5; 2 flocks and 1 environmental hatchery sample), feed (n=1), poultry slaughter house (n=3, environmental sample and meat) were typed with multi locus variable number of tandem repeat analysis (MLVA) and pulsed-field gel electrophoresis (PFGE) to investigate the epidemiology of the outbreak. Based on PFGE results isolates were divided into four groups (Simpson's index of diversity (DI)=0.24±0.15). Due to the low DI, PFGE was not sufficient to provide information to unravel the outbreak. Based on MLVA typing the DT41 (42/47 isolates) and the RDNC isolates (5/47) were split into nine groups (DI=0.65±0.14). When a maximum divergence at one locus was permitted these could be gathered into four groups. Using this criterion, combined with epidemiological information, a spread of one type from broiler breeders to broilers and further to the poultry slaughter house was plausible. In conclusion, although it could be concluded that a spread within the broiler production pyramid had taken place the source of the sudden increase of S. Typhimurium DT41 remains unclear. To investigate this in more detail, further studies using whole genome sequencing to obtain a higher discriminatory strength and including isolates from a longer period of time and from various sources are in progress.

Extensive genetic variability linked to IS26 insertions in the fljB promoter region of atypical monophasic variants of Salmonella enterica serovar Typhimurium

Fifty-nine monophasic Salmonella enterica serovar Typhimurium isolates, collected in Belgium during the period from 2008 to 2011, have been serotyped as 4,[5]:i:- and shown to harbor an fljB coding sequence. The genetic differences between these strains and phenotypically biphasic Salmonella Typhimurium were analyzed through PCR and DNA sequencing. Genetic alterations in the fljB promoter region affecting expression of the phase 2 flagellin were observed in 53 isolates. Other genetic events in the invertible region carrying the fljB promoter were observed in 2 isolates. For the remaining 4 isolates, no molecular differences with a reference biphasic Salmonella Typhimurium strain could be observed. Next-generation sequencing of one representative isolate affected in the fljB promoter region revealed a 26-kb IS26 composite transposon insertion along with a local genomic rearrangement. Several other IS26 element-mediated alterations of this genomic region were observed. This group of monophasic Salmonella Typhimurium isolates was genetically heterogeneous, as revealed by multilocus variable-number tandem-repeat analysis (MLVA), PCR, and sequencing. Pigs and pork represented a major source of such monophasic isolates in Belgium, as reported in other countries. Three out of 5 isolates of human origin presented genetic profiles identical to those of food isolates, demonstrating the pathogenic potential of the newly characterized variants and potential dissemination along the food chain. This study highlighted the key role played by IS26 insertions in the loss of phase 2 flagellin expression and the subsequent generation of multiple monophasic variant lineages from biphasic Salmonella Typhimurium ancestors.

Multi locus variable-number tandem repeat (MLVA) typing tools improved the surveillance of Salmonella enteritidis: a 6 years retrospective study

Surveillance of Salmonella enterica subsp. enterica serovar Enteritidis is generally considered to benefit from molecular techniques like multiple-locus variable-number of tandem repeats analysis (MLVA), which allow early detection and confinement of outbreaks. Here, a surveillance study, including phage typing, antimicrobial susceptibility testing and MLVA on 1,535 S. Enteritidis isolates collected between 2007 and 2012, was used to evaluate the added value of MLVA for public health surveillance in Belgium. Phage types PT4, PT8, PT21, PT1, PT6, PT14b, PT28 and PT13 dominate the Belgian S. Enteritidis population. The isolates of S. Enteritidis were most frequently susceptible to all antibiotics tested. 172 different MLVA profiles were detected, of which 9 frequent profiles included 67.2% of the S. Enteritidis population. During a serial passage experiment on selected isolates to investigate the in vitro stability of the 5 MLVA loci, no variations over time were observed indicating that the MLVA profiles were stable. The MLVA profile of isolates originating from different outbreaks in the Democratic Republic of the Congo (DRC) between 2010 and 2011 were distinct from any of the MLVA profiles found in Belgian isolates throughout the six year observational period and demonstrates that MLVA improves public health surveillance of S. Enteritidis. However, MLVA should be complemented with other subtyping methods when investigating outbreaks is caused by the most common MLVA profile.

Trends in serotype distribution and antimicrobial susceptibility in Salmonella enterica isolates from humans in Belgium, 2009 to 2013

The Belgian National Reference Centre for Salmonella received 16,544 human isolates of Salmonella enterica between January 2009 and December 2013. Although 377 different serotypes were identified, the landscape is dominated by S. enterica serovars Typhimurium (55%) and Enteritidis (19%) in a ratio which is inverse to European Union averages. With outbreaks of Salmonella serotypes Ohio, Stanley, and Paratyphi B variant Java as prime examples, 20 serotypes displayed significant fluctuations in this 5-year period. Typhoid strains account for 1.2% of Belgian salmonellosis cases. Large-scale antibiotic susceptibility analyses (n = 4,561; panel of 12 antibiotics) showed declining resistance levels in S. Enteritis and Typhimurium isolates for 8 and 3 tested agents, respectively. Despite low overall resistance to ciprofloxacin (4.4%) and cefotaxime (1.6%), we identified clonal lineages of Salmonella serotypes Kentucky and Infantis displaying rising resistance against these clinically important drugs. Quinolone resistance is mainly mediated by serotype-specific mutations in GyrA residues Ser83 and Asp87 (92.2% not wild type), while an additional ParC_Ser80Ile mutation leads to ciprofloxacin resistance in 95.5% S. Kentucky isolates, which exceeds European averages. Plasmid-mediated quinolone resistance (PMQR) alleles qnrA1 (n = 1), qnrS (n = 9), qnrD1 (n = 4), and qnrB (n = 4) were found in only 3.0% of 533 isolates resistant to nalidixic acid. In cefotaxime-resistant isolates, we identified a broad range of Ambler class A and C β-lactamase genes (e.g., bla(SHV-12), blaTEM-52, bla(CTX-M-14), and bla(CTX-M-15)) commonly associated with members of the family Enterobacteriaceae. In conclusion, resistance to fluoroquinolones and cefotaxime remains rare in human S. enterica, but clonal resistant serotypes arise, and continued (inter)national surveillance is mandatory to understand the origin and routes of dissemination thereof.

Characterization of an unusual Salmonella phage type DT7a and report of a foodborne outbreak of salmonellosis

Salmonella enterica subsp. enterica serovar 4,[5],12,i:- is a monophasic variant of Salmonella Typhimurium and its occurrence has markedly increased in several European countries in the last ten years. In June 2011, an outbreak of Salmonella 4,[5],12,i:- was reported among attendees of a wedding reception in the North-East of Italy. The source of this outbreak was identified as a cooked pork product served during the wedding reception. All Salmonella isolates from humans and the contaminated pork products were identified as Salmonella 4,[5],12,i:- and phage typed as DT7a. Afterwards, the farm where the pigs were raised was identified and sampled, and Salmonella Typhimurium was isolated from swine fecal samples. Despite the difference in serovar, these Salmonella Typhimurium isolates were also phage typed as DT7a. In the present study, Salmonella isolates from animals, humans and pork products during the outbreak investigation were subtyped by pulsed-field gel electrophoresis (PFGE), Multiple-Locus Variable number tandem repeats Analysis (MLVA), and resistance patterns, aiming to identify the most suitable subtyping methods to characterize isolates associated with this outbreak. In addition, a collection of epidemiologically unrelated strains of Salmonella 4,[5],12,i:- and Salmonella Typhimurium sharing the same phage type (DT7a) was similarly characterized in order to investigate their genetic relationship. This study provides a first snapshot of a rare Salmonella phage type, DT7a, associated with both Salmonella 4,[5],12,i:- and Salmonella Typhimurium. Moreover, the study demonstrated that in this specific context MLVA could be a reliable tool to support outbreak investigations as well as to assess the genetic relatedness among Salmonella isolates.

Analysis of Salmonella enterica serovar Typhimurium variable-number tandem-repeat data for public health investigation based on measured mutation rates and whole-genome sequence comparisons

Variable-number tandem repeats (VNTRs) mutate rapidly and can be useful markers for genotyping. While multilocus VNTR analysis (MLVA) is increasingly used in the detection and investigation of food-borne outbreaks caused by Salmonella enterica serovar Typhimurium (S. Typhimurium) and other bacterial pathogens, MLVA data analysis usually relies on simple clustering approaches that may lead to incorrect interpretations. Here, we estimated the rates of copy number change at each of the five loci commonly used for S. Typhimurium MLVA, during in vitro and in vivo passage. We found that loci STTR5, STTR6, and STTR10 changed during passage but STTR3 and STTR9 did not. Relative rates of change were consistent across in vitro and in vivo growth and could be accurately estimated from diversity measures of natural variation observed during large outbreaks. Using a set of 203 isolates from a series of linked outbreaks and whole-genome sequencing of 12 representative isolates, we assessed the accuracy and utility of several alternative methods for analyzing and interpreting S. Typhimurium MLVA data. We show that eBURST analysis was accurate and informative. For construction of MLVA-based trees, a novel distance metric, based on the geometric model of VNTR evolution coupled with locus-specific weights, performed better than the commonly used simple or categorical distance metrics. The data suggest that, for the purpose of identifying potential transmission clusters for further investigation, isolates whose profiles differ at one of the rapidly changing STTR5, STTR6, and STTR10 loci should be collapsed into the same cluster.

Molecular typing of monophasic Salmonella 4,[5]:i:- strains isolated in Belgium (2008-2011)

To assess the distribution of Salmonella 4,[5]:i:- subtypes in the Belgian food chain and compare it to the subtypes associated with human infections, a molecular assessment was initiated. Two hundred fifty-three Salmonella isolates serotyped as 4,[5]:i:- during the period 2008-2011 in Belgium and originating from animal productions, food or human clinical samples were analysed by a specific duplex PCR. One hundred ninety-four isolates (76.7%) fit the profile of a S. Typhimurium monophasic variant as defined by the European Food Safety Authority. The other isolates possessed but did not express the phase II flagellin gene (23.3%). Multiple Locus Variable Number of Tandem Repeats Analysis (MLVA) revealed many but closely related profiles in the fljB-negative S. Typhimurium monophasic variant isolates. Some MLVA types were associated with both human and animal isolates but no unique source of human contamination could be demonstrated.