Evaluation of whole genome sequencing for outbreak detection of Salmonella enterica

Salmonella enterica is a common cause of minor and large food borne outbreaks. To achieve successful and nearly 'real-time' monitoring and identification of outbreaks, reliable sub-typing is essential. Whole genome sequencing (WGS) shows great promises for using as a routine epidemiological typing tool. Here we evaluate WGS for typing of S. Typhimurium including different approaches for analyzing and comparing the data. A collection of 34 S. Typhimurium isolates was sequenced. This consisted of 18 isolates from six outbreaks and 16 epidemiologically unrelated background strains. In addition, 8 S. Enteritidis and 5 S. Derby were also sequenced and used for comparison. A number of different bioinformatics approaches were applied on the data; including pan-genome tree, k-mer tree, nucleotide difference tree and SNP tree. The outcome of each approach was evaluated in relation to the association of the isolates to specific outbreaks. The pan-genome tree clustered 65% of the S. Typhimurium isolates according to the pre-defined epidemiology, the k-mer tree 88%, the nucleotide difference tree 100% and the SNP tree 100% of the strains within S. Typhimurium. The resulting outcome of the four phylogenetic analyses were also compared to PFGE revealing that WGS typing achieved the greater performance than the traditional method. In conclusion, for S. Typhimurium, SNP analysis and nucleotide difference approach of WGS data seem to be the superior methods for epidemiological typing compared to other phylogenetic analytic approaches that may be used on WGS. These approaches were also superior to the more classical typing method, PFGE. Our study also indicates that WGS alone is insufficient to determine whether strains are related or un-related to outbreaks. This still requires the combination of epidemiological data and whole genome sequencing results

Multiplex PCR for detection of plasmid-mediated colistin resistance determinants, mcr-1, mcr-2, mcr-3, mcr-4 and mcr-5 for surveillance purposes

Background and aim: Plasmid-mediated colistin resistance mechanisms have been identified worldwide in the past years. A multiplex polymerase chain reaction (PCR) protocol for detection of all currently known transferable colistin resistance genes (mcr-1 to mcr-5, and variants) in Enterobacteriaceae was developed for surveillance or research purposes. Methods: We designed four new primer pairs to amplify mcr-1, mcr-2, mcr-3 and mcr-4 gene products and used the originally described primers for mcr-5 to obtain a stepwise separation of ca 200 bp between ampli-cons. The primer pairs and amplification conditions allow for single or multiple detection of all currently described mcr genes and their variants present in Enterobacteriaceae. The protocol was validated testing 49 European Escherichia coli and Salmonella isolates of animal origin. Results: Multiplex PCR results in bovine and porcine isolates from Spain, Germany, France and Italy showed full concordance with whole genome sequence data. The method was able to detect mcr-1, mcr-3 and mcr-4 as singletons or in different combinations as they were present in the test isolates. One new mcr-4 variant, mcr-4.3, was also identified. Conclusions: This method allows rapid identification of mcr-positive bacteria and overcomes the challenges of phenotypic detection of colistin resistance. The multiplex PCR should be particularly interesting in settings or laboratories with limited resources for performing genetic analysis as it provides information on the mechanism of colistin resistance without requiring genome sequencing. © 2018, European Centre for Disease Prevention and Control (ECDC). All rights reserved

Virulence related sequences: insights provided by comparative genomics of Streptococcus uberis of differing virulence

Background: Streptococcus uberis, a Gram-positive, catalase-negative member of the family Streptococcaceae is an important environmental pathogen responsible for a significant proportion of subclinical and clinical bovine intramammary infections. Currently, the genome of only a single reference strain (0140J) has been described. Here we present a comparative analysis of complete draft genome sequences of an additional twelve S. uberis strains. Results: Pan and core genome analysis revealed the core genome common to all strains to be 1,550 genes in 1,509 orthologous clusters, complemented by 115-246 accessory genes present in one or more S. uberis strains but absent in the reference strain 0140J. Most of the previously predicted virulent genes were present in the core genome of all 13 strains but gene gain/loss was observed between the isolates in CDS associated with clustered regularly interspaced short palindromic repeats (CRISPRs), prophage and bacteriocin production. Experimental challenge experiments confirmed strain EF20 as non-virulent; only able to infect in a transient manner that did not result in clinical mastitis. Comparison of the genome sequence of EF20 with the validated virulent strain 0140J identified genes associated with virulence, however these did not relate clearly with clinical/non-clinical status of infection. Conclusion: The gain/loss of mobile genetic elements such as CRISPRs and prophage are a potential driving force for evolutionary change. This first “whole-genome” comparison of strains isolated from clinical vs non-clinical intramammary infections including the type virulent vs non-virulent strains did not identify simple gene gain/loss rules that readily explain, or be confidently associated with, differences in virulence. This suggests that a more complex dynamic determines infection potential and clinical outcome not simply gene content

The invasome of Salmonella Dublin as revealed by whole genome sequencing

Background Salmonella enterica serovar Dublin is a zoonotic infection that can be transmitted from cattle to humans through consumption of contaminated milk and milk products. Outbreaks of human infections by S. Dublin have been reported in several countries including high-income countries. A high proportion of S. Dublin cases in humans are associated with invasive disease and systemic illness. The genetic basis of virulence in S. Dublin is not well characterized. Methods Whole genome sequencing was applied to a set of clinical invasive and non-invasive S. Dublin isolates from different countries in order to characterize the putative genetic determinants involved in the virulence and invasiveness of S. Dublin in humans. Results We identified several virulence factors that form the bacterial invasome and may contribute to increasing bacterial virulence and pathogenicity including mainly Gifsy-2 prophage, two different type 6 secretion systems (T6SSs) harbored by Salmonella pathogenicity islands; SPI-6 and SPI-19 respectively and virulence genes; ggt and PagN. Although Vi antigen and the virulence plasmid have been reported previously to contribute to the virulence of S. Dublin we did not detect them in all invasive isolates indicating that they are not the main virulence determinants in S. Dublin. Conclusion Several virulence factors within the genome of S. Dublin might contribute to the ability of S. Dublin to invade humans’ blood but there were no genomic markers that differentiate invasive from non-invasive isolates suggesting that host immune response play a crucial role in the clinical outcome of S. Dublin infection

The Role of the st313-td Gene in Virulence of Salmonella Typhimurium ST313

Multidrug-resistant Salmonella enterica serovar Typhimurium ST313 has emerged in sub-Saharan Africa causing severe infections in humans. Therefore, it has been speculated that this specific sequence type, ST313, carries factors associated with increased pathogenicity. We assessed the role in virulence of a gene with a yet unknown function, st313-td, detected in ST313 through comparative genomics. Additionally, the structure of the genomic island ST313-GI, harbouring the gene was determined. The gene st313-td was cloned into wild type S. Typhimurium 4/74 (4/74-C) as well as knocked out in S. Typhimurium ST313 02-03/002 (Δst313-td) followed by complementation (02-03/002-C). Δst313-td was less virulent in mice following i.p. challenge than the wild type and this phenotype could be partly complemented in trans, indicating that st313-td plays a role during systemic infection. The gene st313-td was shown not to affect invasion of cultured epithelial cells, while the absence of the gene significantly affects uptake and intracellular survival within macrophages. The gene st313-td was proven to be strongly associated to invasiveness, harboured by 92.5% of S. Typhimurium blood isolates (n = 82) and 100% of S. Dublin strains (n = 50) analysed. On the contrary, S. Typhimurium isolates of animal and food origin (n = 82) did not carry st313-td. Six human, non-blood isolates of S. Typhimurium from Belarus, China and Nepal harboured the gene and belonged to sequence types ST398 and ST19. Our data showed a global presence of the st313-td gene and in other sequence types than ST313. The gene st313-td was shown to be expressed during logarithmic phase of growth in 14 selected Salmonella strains carrying the gene. This study reveals that st313-td plays a role in S. Typhimurium ST313 pathogenesis and adds another chapter to understanding of the virulence of S. Typhimurium and in particular of the emerging sequence type ST313

Global monitoring of antimicrobial resistance based on metagenomics analyses of urban sewage

Antimicrobial resistance (AMR) is a serious threat to global public health, but obtaining representative data on AMR for healthy human populations is difficult. Here, we use meta-genomic analysis of untreated sewage to characterize the bacterial resistome from 79 sites in 60 countries. We find systematic differences in abundance and diversity of AMR genes between Europe/North-America/Oceania and Africa/Asia/South-America. Antimicrobial use data and bacterial taxonomy only explains a minor part of the AMR variation that we observe. We find no evidence for cross-selection between antimicrobial classes, or for effect of air travel between sites. However, AMR gene abundance strongly correlates with socio-economic, health and environmental factors, which we use to predict AMR gene abundances in all countries in the world. Our findings suggest that global AMR gene diversity and abundance vary by region, and that improving sanitation and health could potentially limit the global burden of AMR. We propose metagenomic analysis of sewage as an ethically acceptable and economically feasible approach for continuous global surveillance and prediction of AMR.Peer reviewe