Search

Objectives Implementing whole-genome sequencing (WGS) technologies in clinical microbiology laboratories can increase the amount and quality of information available for healthcare practitioners. In this study, we analysed the applicability of this method and determined the distribution of bacterial species processed in clinical settings in Denmark. Methods We performed a point-prevalence study of all bacterial isolates (n = 2,009) processed and reported in the Clinical Microbiology Laboratories in Denmark in one day in January 2018. We compared species identification as performed by classical methods (MALDI-TOF) and by bioinformatics analysis (KmerFinder and rMLST) of WGS (Illumina NextSeq) data. We compared the national point-prevalence of bacterial isolates observed in clinical settings with the research attention given to those same genera in scientific literature. Results The most prevalent bacterium was Escherichia coli isolated from urine (n = 646), followed by Staphylococcus spp. from skin or soft tissues (n = 197). The distribution of bacterial species throughout the country was not homogeneous. We observed concordance of species identification for all methods in 95.7% (n = 1,919) of isolates, furthermore obtaining concordance for 99.7% (n = 1,999) at genus level. The number of scientific publications in the country did not correlate with the number of bacterial isolates of each genera analysed in this study. Conclusions WGS technologies have the potential to be applied in clinical settings for routine diagnostics purposes. This study also showed that bioinformatics databases should be continuously improved and results from local point-prevalence surveys should not be applied at national levels without previously determining possible regional variations.

Escherichia coli sequence type 131 (ST131) is known for its contribution to multidrug resistance and the worldwide spread of this clone has become a global problem. Understanding the trends among ST131 clades will help design strategies to prevent its rapid dissemination. In this study, 72 ST131 strains were subjected to comparative genomic analysis and 64 clade C strains were compared with clade C strains reported from other regions using publicly available whole-genome sequencing data. C1 (n = 31 [48.4%]) and C2 (n = 33 [%51.5]) strains had the same prevalence in our collection, and C1-M27 (n = 22) strains were closely related, carried a unique plasmid type (F1:A2:B20), and exhibited virotype C. Removal of 11 C2 strains with varied virotype patterns and the heterogeneous IncF type identified 22 closely related virotype E/F strains with replicon type F31/F36:A4:B1, forming what we denote as the “C2-subset.” In a global context, the C2-subset constituted a distinct cluster with international virotype E strains and harbored a genomic island, GI-pheU. Association of cnf1/hlyCABD genes with 1 to 7 mobile genetic elements, mostly IS682/ISKpn37 combination within GI-pheU was identified. The C2-subset accounted for excess resistance/virulence of subclade C2 relative to C1 strains. In addition, a conserved chromosomal IS26-mediated composite transposon (IS15DIV-ISEcp1-blaCTX-M-15-WbuC cupin fold metalloprotein-Tn2-IS15DIV) was observed in the C2-subset. The local spread of the C2-subset in the hospital studied, with the carriage of higher virulence/resistance markers and a peculiar F-type plasmid, demonstrates the potential for diversification of the ST131 lineage and the emergence of subpopulations with higher survival potential to cause health care-associated outbreaks.IMPORTANCE Escherichia coli sequence type 131 (ST131) is a globally dominant multidrug-resistant clone that is commonly associated with extraintestinal infections. Specific sublineages have been shown to have emerged and spread within ST131, highlighting the complex nature of ST131 epidemiology. This study systematically compared the Iranian ST131 population to those reported from other countries and found a subpopulation harboring virotype E, a homogeneous profile of plasmid Inc-F type F31/F36:A4:B1 harboring cnf1/hemolysin genes on the genomic island GI-pheU, and up to seven mobile genetic elements (MGEs) flanking cnf1/hemolysin virulence markers. The results of this study highlight the importance of MGEs for virulence gene acquisition and the formation of new subpopulations among pandemic clones such as E. coli ST131.

Carbapenems are critically important antimicrobials for the treatment of clinical infections caused by MDR Enterobacteriaceae, with colistin reserved for those cases in which even carbapenems are ineffective. Therefore, co-localization of carbapenem and colistin resistance genes on the same plasmid is worrisome and, to date, has been described only in China in an Escherichia coli from a diseased chicken (blaNDM-4 and mcr-1 on an IncHI2/ST3 plasmid)1 and in an E. coli from a diseased pet cat (blaNDM-5 and mcr-1 on an IncX3-X4 hybrid plasmid).2 Here, we report a novel case of co-localization of blaOXA-162 carbapenemase and mcr-1 on a transferable IncHI2 plasmid in E. coli from a healthy broiler chicken in Europe.

Objectives Fully sequenced IncI1 plasmids obtained from CTX-M-1-producing Escherichia coli of human and animal origin were compared. Methods Twelve E. coli isolates sharing identical ESBL genes and plasmid multilocus STs sequenced on Illumina and MinION platforms were obtained from the Danish antimicrobial resistance surveillance programme, DANMAP. After de novo assembly, the sequences of plasmids harbouring blaCTX-M-1 were manually curated and ORFs annotated. Within-group comparisons were performed separately for the IncI1 ST3 plasmid type and the IncI1 ST7 plasmid type. The IncI1 ST3 plasmid group was obtained from 10 E. coli isolates (2 from patients with bloodstream infections, 6 from food and 2 from animals). The IncI1 ST7 plasmids originated from E. coli isolates obtained from a patient with bloodstream infection and from a pig. Sequences of IncI1 ST3 and IncI1 ST7 plasmids harbouring blaCTX-M-1 with determined origin were retrieved from GenBank and used for comparison within the respective group. Results The 10 IncI1 ST3 blaCTX-M-1 plasmids were highly similar in structure and organization with only minor plasmid rearrangements and differences in the variable region. The IncI1 ST7 blaCTX-M-1 plasmids also showed high similarity in structure and organization. The high level of similarity was also observed when including plasmids from E. coli of animal origin from Australia, Switzerland, the Netherlands and France. Conclusions This study shows broad spread of a very successful CTX-M-1-producing IncI1 type plasmid among E. coli of both human and animal origin.

Food authentication is a rapidly growing field driven by increasing public awareness of food quality and safety. Foods containing herbs are particularly prone to industrial fraud and adulteration. Several methodologies are currently used to evaluate food authenticity. DNA-based technologies have increased focus, with DNA barcoding the most widely used. DNA barcoding is based on the sequencing and comparison of orthologous DNA regions from all species in a sample, but the approach is limited by its low resolution to distinguish closely-related species. Here we developed a customised database and bioinformatics pipeline (Herbs Authenticity - GitHub) to identify herbal ingredients implemented as a metagenomics approach for plant-derived product authenticity testing. We evaluated the accuracy of the method by using publicly available plant genomes and databases to allow the construction of our customised database barcodes, which were also complemented with entries from publicly available resources (iBOL and ENA). The pipeline performance was then tested with new 47 de novo partly sequenced whole plant genomes or barcodes as query sequences. Our results show that using our mapping algorithm with the customised barcode database correctly identifies the main components of a wide range of plant-derived samples, albeit with variable additional noise across samples depending on the tested samples and barcodes. Our result also show that at the current stage the usefulness of metagenomics is limited by the availability of reference sequences and the needed sequencing depth. However, this method shows promise for evaluating the authenticity of different herbal products provided that the method is further refined to increase the qualitative and quantitative accuracy.

Background A novel Brucella strain closely related to Brucella (B.) melitensis biovar (bv) 3 was found in Croatian cattle during testing within a brucellosis eradication programme. Case presentation Standardised serological, brucellin skin test, bacteriological and molecular diagnostic screening for Brucella infection led to positive detection in one dairy cattle herd. Three isolates from that herd were identified to species level using the Bruce ladder method. Initially, two strains were typed as B. melitensis and one as B. abortus, but multiplex PCR based on IS711 and the Suis ladder showed that all of them to belong to B. melitensis, and the combination of whole-genome and multi-locus sequencing as well as Multi-Locus Variable numbers of tandem repeats Analysis (MLVA) highlighted a strong proximity within the phylogenetic branch of B. melitensis strains previously isolated from Croatia, Albania, Kosovo and Bosnia and Herzegovina. Two isolates were determined to be B. melitensis bv. 3, while the third showed a unique phylogenetic profile, growth profile on dyes and bacteriophage typing results. This isolate contained the 609-bp omp31 sequence, but not the 723-bp omp31 sequence present in the two isolates of B. melitensis bv. 3. Conclusions Identification of a novel Brucella variant in this geographic region is predictable given the historic endemicity of brucellosis. The emergence of a new variant may reflect a combination of high prevalence among domestic ruminants and humans as well as weak eradication strategies. The zoonotic potential, reservoirs and transmission pathways of this and other Brucella variants should be explored.

Additional file 2. Textual description of phylogenetic tree.

Candida albicans and Candida glabrata are opportunistic fungal pathogens with increasing incidence worldwide and higher-than-expected prevalence in Denmark. We whole-genome sequenced yeast isolates collected from Danish Clinical Microbiology Laboratories to obtain an overview of the Candida population in the country. The majority of the 30 C. albicans isolates were found to belong to three globally prevalent clades, and, with one exception, the remaining isolates were also predicted to cluster with samples from other geographical locations. Similarly, most of the eight C. glabrata isolates were predicted to be prevalent subtypes. Antifungal susceptibility testing proved all C. albicans isolates to be susceptible to both azoles and echinocandins. Two C. glabrata isolates presented azole-resistant phenotypes, yet all were susceptible to echinocandins. There is no indication of causality between population structure and resistance phenotypes for either species.

Bacteriophages, or phages, are ubiquitous bacterial and archaeal viruses with an estimated total global population of 1031. It is well-known that wherever there are bacteria, their phage counterparts will be found, aiding in shaping the bacterial population. The present study used metagenomic data from global influent sewage in 79 cities in 60 countries to identify phages associated with bacteria and to explore their potential role in antimicrobial resistance gene (ARG) dissemination. The reads were mapped to known databases for bacteriophages and their abundances determined and correlated to geographic origin and the countries socio-economic status, as well as the abundances of bacterial species and ARG. We found that some phages were not equally distributed on a global scale, but their distribution was rather dictated by region and the socioeconomic status of the specific countries. This study provides a preliminary insight into the global and regional distribution of phages and their potential impact on the transmission of ARGs between bacteria. Moreover, the findings may indicate that phages in sewage could have adopted a lytic lifestyle, meaning that most may not be associated with bacteria and instead may be widely distributed as free-living phages, which are known to persist longer in the environment than their hosts. In addition, a significant correlation between phages and ARGs was obtained, indicating that phages may play a role in ARG dissemination. However, further analyses are needed to establish the true relationship between phages and ARGs due to a low abundance of the phages identified.

Additional file 3. Figure Phylogenetic Tree.

It has been 30 years since the initial emergence and subsequent rapid global spread of multidrug-resistant Salmonella enterica serovar Typhimurium DT104 (MDR DT104). Nonetheless, its origin and transmission route have never been revealed. We used whole-genome sequencing (WGS) and temporally structured sequence analysis within a Bayesian framework to reconstruct temporal and spatial phylogenetic trees and estimate the rates of mutation and divergence times of 315 S . Typhimurium DT104 isolates sampled from 1969 to 2012 from 21 countries on six continents. DT104 was estimated to have emerged initially as antimicrobial susceptible in ∼1948 (95% credible interval [CI], 1934 to 1962) and later became MDR DT104 in ∼1972 (95% CI, 1972 to 1988) through horizontal transfer of the 13-kb Salmonella genomic island 1 (SGI1) MDR region into susceptible strains already containing SGI1. This was followed by multiple transmission events, initially from central Europe and later between several European countries. An independent transmission to the United States and another to Japan occurred, and from there MDR DT104 was probably transmitted to Taiwan and Canada. An independent acquisition of resistance genes took place in Thailand in ∼1975 (95% CI, 1975 to 1990). In Denmark, WGS analysis provided evidence for transmission of the organism between herds of animals. Interestingly, the demographic history of Danish MDR DT104 provided evidence for the success of the program to eradicate Salmonella from pig herds in Denmark from 1996 to 2000. The results from this study refute several hypotheses on the evolution of DT104 and suggest that WGS may be useful in monitoring emerging clones and devising strategies for prevention of Salmonella infections.

Background Sequence type 131 (ST131) of Escherichia coli is a pandemic clone that drives the increasing rates of antibiotic resistance. While the pervasiveness of ST131 clade C, especially subclades C2 and C1-M27, has been demonstrated in numerous global surveys, no report about the ST131 clades and their virotypes has been published from Iran so far. Methods A collection of 73 consecutive ST131 isolates from extraintestinal specimens was investigated for determination of virotypes, antibiotic susceptibility patterns, resistance/virulence determinants, and clade subsets. Results Most of the isolates belonged to subclade C2 (33/73; 45.2%), which had the highest virulence factor (VF) scores and resistance rates, followed by C1-M27 (18; 24.6%), C1-non-M27 (14; 19.1%), and A (8; 10.9%). The distinctive profiles of subclade C2 virulence genes were revealed by principle coordinates analysis testing. The distribution of the hlyA virulence gene among subclade C2 was not uniform, so that positive strains (21; 63.6%) showed significantly higher rates of resistance (blaCTX-M-15, blaOXA-1, aac(6’)-Ib-cr, aac(6’)-Ib, aac(3)-IIa) and virulence (hra, tia/hek, K5, cnf, papGII, papC) markers and gentamicin/tobramycin resistance. Virotype C as the most common virotype (34; 46.5%) was predominant among the subclade C1 population, while virotypes E and F (21; 28.7%) were detected among subclade C2, which had the highest VF scores and aminoglycoside resistance rates. Conclusions The appearance of virotypes E and F among subclade C2 strains with higher rates of aminoglycoside resistance/virulence gene content shows the shifting dynamics of this pandemic clone in response to antibiotic selection pressure by establishing subsets with higher survival potential.

Background Typhoid fever, caused by S. enterica ser. Typhi, continues to be a substantial health burden in developing countries. Little is known of the genotypic diversity of S. enterica ser. Typhi in Zimbabwe, but this is key for understanding the emergence and spread of this pathogen and devising interventions for its control. Objectives To report the molecular epidemiology of S. enterica ser. Typhi outbreak strains circulating from 2012 to 2019 in Zimbabwe, using comparative genomics. Methods A review of typhoid cases records from 2012 to 2019 in Zimbabwe was performed. The phylogenetic relationship of outbreak isolates from 2012 to 2019 and emergence of antibiotic resistance was investigated by whole-genome sequence analysis. Results A total 22 479 suspected typhoid cases, 760 confirmed cases were reported from 2012 to 2019 and 29 isolates were sequenced. The majority of the sequenced isolates were predicted to confer resistance to aminoglycosides, β-lactams, phenicols, sulphonamides, tetracycline and fluoroquinolones (including qnrS detection). The qnrS1 gene was associated with an IncN (subtype PST3) plasmid in 79% of the isolates. Whole-genome SNP analysis, SNP-based haplotyping and resistance determinant analysis showed that 93% of the isolates belonged to a single clade represented by multidrug-resistant H58 lineage I (4.3.1.1), with a maximum pair-wise distance of 22 SNPs. Conclusions This study has provided detailed genotypic characterization of the outbreak strain, identified as S. Typhi 4.3.1.1 (H58). The strain has reduced susceptibility to ciprofloxacin due to qnrS carried by an IncN (subtype PST3) plasmid resulting from ongoing evolution to full resistance.

Salmonella Infantis is one of the five serovars most frequently causing human salmonellosis in Europe, mainly associated with poultry. A clone harbouring a conjugative plasmid of emerging S. Infantis (pESI)-like megaplasmid, carrying multidrug resistant (MDR) and extended-spectrum beta-lactamases (ESBL) genes, has spread in the Italian broiler chicken industry also causing human illness. This work is aimed at elucidating the molecular epidemiology of S. Infantis and pESI-like in Europe using whole-genome sequencing and bioinformatics analysis, and to investigate the genetic relatedness of S. Infantis clones and pESI-like from animals, meat, feed and humans provided by institutions of nine European countries. Two genotyping approaches were used: chromosome or plasmid SNP-based analysis and the minimum spanning tree (MST) algorithm based on core-genome multilocus sequence typing (cgMLST). The European S. Infantis population appeared heterogeneous, with different genetic clusters defined at core-genome level. However, pESI-like variants present in 64.1 % of the isolates were more genetically homogeneous and capable of infecting different clonal lineages in most of the countries. Two different pESI-like with ESBL genes (n=82) were observed: bla CTX-M-1-positive in European isolates and bla CTX-M-65-positive in American isolates (study outgroup). Both variants had toxin-antitoxin systems, resistance genes towards tetracyclines, trimethoprim, sulphonamides and aminoglycosides, heavy metals (merA) and disinfectants (qacEΔ). Worryingly, 66 % of the total isolates studied presented different gyrA chromosomal point mutations associated with (fluoro)quinolone resistance (MIC range 0.125–0.5 mg/L), while 18 % displayed transferable macrolide resistance mediated by mph, mef and erm(B) genes. Proper intervention strategies are needed to prevent further dissemination/transmission of MDR S. Infantis and pESI-like along the food chain in Europe.

Zoonotic Salmonella causes millions of human salmonellosis infections worldwide each year. Information about the source of the bacteria guides risk managers on control and preventive strategies. Source attribution is the effort to quantify the number of sporadic human cases of a specific illness to specific sources and animal reservoirs. Source attribution methods for Salmonella have so far been based on traditional wet-lab typing methods. With the change to whole genome sequencing there is a need to develop new methods for source attribution based on sequencing data. Four European datasets collected in Denmark (DK), Germany (DE), the United Kingdom (UK) and France (FR) are presented in this descriptor. The datasets contain sequenced samples of Salmonella Typhimurium and its monophasic variants isolated from human, food, animal and the environment. The objective of the datasets was either to attribute the human salmonellosis cases to animal reservoirs or to investigate contamination of the environment by attributing the environmental isolates to different animal reservoirs., Measurement(s)sequence_assembly • DNATechnology Type(s)sequence assembly process • DNA sequencing • Paired-End SequencingFactor Type(s)geographic location • Salmonella Typhimurium isolate sourceSample Characteristic - OrganismSalmonellaSample Characteristic - LocationKingdom of Denmark • Germany • United Kingdom • French Republic Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.11816748

Salmonella Weltevreden is increasingly reported from aquatic environments, seafood, and patients in several Southeast Asian countries. Using genome-wide analysis, we characterized S. Weltevreden isolated from cultured shrimp and tilapia from Vietnam and China to study their genetic characteristics and relatedness to clinical isolates of S. Weltevreden ST-365. The phylogenetic analysis revealed up to 312 single-nucleotide polymorphism (SNP) difference between tilapia isolates, whereas isolates from shrimp were genetically more closely related. Epidemiologically unrelated isolates from Vietnam were closely related to isolates from China, e.g., 20 SNPs differences between strains 28V and 75C. In comparison with strains from other parts of the world, our environmental isolates predominantly clustered within the continental South Asia lineage, constituted mostly of strains from human stool with as low as seven SNPs difference, e.g., 30V versus Cont_ERR495254. All sequenced isolates were MLST type ST-365 and contained the major virulence-related genes encoded by the Salmonella Pathogenicity Islands 1-5. Ten of the isolates harbored the IncFII(S) plasmid similar to the virulence genes-mediated plasmid pSPCV of S. Paratyphi C, and one isolate had the IncQ1 plasmid on the same contig with strA/B, sul2, and tetA resistance genes similar to the IncQ1 type, pNUC of S. Typhimurium. A pangenomic analysis yielded 7891 genes including a core genome of 4892 genes, with a closely related accessory genome content between clinical and environmental isolates (Benjamini p > 0.05). In a search for differences that could explain the higher prevalence of S. Weltevreden in aquatic samples, genomes were compared with those of other Salmonella enterica serovars. S. Weltevreden revealed specific regions harboring glpX (Fructose-1;6-bisphosphatase; class II), rfbC (dTDP-4-dehydrorhamnose 3;5-epimerase), and cmtB (PTS Mannitol-specific cryptic phosphotransferase enzyme IIA component) involved in carbohydrate biosynthesis pathways. Our study builds grounds for future experiments to determine genes or pathways that are essential when S. Weltevreden are in aquatic environments and microbial interactions providing survival advantages to S. Weltevreden in such environments. Published version

Prevention of the emergence and spread of foodborne diseases is an important prerequisite for the improvement of public health. Source attribution models link sporadic human cases of a specific illness to food sources and animal reservoirs. With the next generation sequencing technology, it is possible to develop novel source attribution models. We investigated the potential of machine learning to predict the animal reservoir from which a bacterial strain isolated from a human salmonellosis case originated based on whole‐genome sequencing. Machine learning methods recognize patterns in large and complex data sets and use this knowledge to build models. The model learns patterns associated with genetic variations in bacteria isolated from the different animal reservoirs. We selected different machine learning algorithms to predict sources of human salmonellosis cases and trained the model with Danish Salmonella Typhimurium isolates sampled from broilers (n = 34), cattle (n = 2), ducks (n = 11), layers (n = 4), and pigs (n = 159). Using cgMLST as input features, the model yielded an average accuracy of 0.783 (95% CI: 0.77–0.80) in the source prediction for the random forest and 0.933 (95% CI: 0.92–0.94) for the logit boost algorithm. Logit boost algorithm was most accurate (valid accuracy: 92%, CI: 0.8706–0.9579) and predicted the origin of 81% of the domestic sporadic human salmonellosis cases. The most important source was Danish produced pigs (53%) followed by imported pigs (16%), imported broilers (6%), imported ducks (2%), Danish produced layers (2%), Danish produced cattle and imported cattle (

Antimicrobial resistance (AMR) in bacteria is a complex subject, why one need to look at this phenomenon from a wider and holistic perspective. The extensive use of the same antimicrobial classes in human and veterinary medicine as well as horticulture is one of the main drivers for the AMR selection. Here, we applied shotgun metagenomics to investigate the AMR epidemiology in several animal species including farm animals, which are often exposed to antimicrobial treatment opposed to an unique set of wild animals that seems not to be subjected to antimicrobial pressure. The comparison of the domestic and wild animals allowed to investigate the possible anthropogenic impact on AMR spread. Inclusion of animals with different feeding behaviors (carnivores, omnivores) enabled to further assess which AMR genes that thrives within the food chain. We tested fecal samples not only of intensively produced chickens, turkeys, and pigs, but also of wild animals such as wild boars, red foxes, and rodents. A multi-directional approach mapping obtained sequences to several databases provided insight into the occurrence of the different AMR genes. The method applied enabled also analysis of other factors that may influence AMR of intestinal microbiome such as diet. Our findings confirmed higher levels of AMR in farm animals than in wildlife. The results also revealed the potential of wildlife in the AMR dissemination. Particularly in red foxes, we found evidence of several AMR genes conferring resistance to critically important antimicrobials like quinolones and cephalosporins. In contrast, the lowest abundance of AMR was observed in rodents originating from natural environment with presumed limited exposure to antimicrobials. Shotgun metagenomics enabled us to demonstrate that discrepancies between AMR profiles found in the intestinal microbiome of various animals probably resulted from the different antimicrobial exposure, habitats, and behavior of the tested animal species.

The sources and modes of transmission of non-typhoidal Salmonella particularly zoonotic transmission are poorly understood in Africa. This study compared phenotypic and genotypic characteristics of Salmonellae isolated from cattle and humans. Faecal samples of diarrhoeic patients (n = 234), and a healthy population (n = 160), beef cattle at slaughter (n = 250), farms (n = 72) and market (n = 100) were cultured for salmonellae and serotyping and antimicrobial susceptibility were determined. Whole-genome sequence typing (WGST) of selected isolates and bioinformatic analysis were used to identify the multilocus sequence type (MLST), plasmid replicons, antimicrobial resistance genes and genetic relatedness by single nucleotide polymorphism (SNP) analysis. The Salmonella isolates, diarrhoeic patients (n = 17), healthy population (n = 13), cattle (abattoir, n = 67; farms, n = 10; market n = 5), revealed 49 serovars; some serovars were common to humans and cattle. Rare serovars were prevalent: Colindale (cattle and humans); Rubislaw and Bredeney (humans); and Dublin, Give, Eastbourne, Hadar, Marseille, Sundsvall, Bergen, Ekotedo, Carno and Ealing (cattle). The sequence types (ST) include ST 584, ST 198, ST 562 and ST 512 for S. Colindale, S. Kentucky S. Rubislaw and S. Urbana, respectively. Clonal cluster shared by cattle and human WGST isolates was not found. Antimicrobial resistance rates were generally low and towards only chloramphenicol, ampicillin, gentamicin, ciprofloxacin, tetracycline and streptomycin, range 2.7% (chloramphenicol) to 8.9% (streptomycin). Multiply resistant isolates included serovars Kentucky, 4,5,12:i:- and Typhimurium. The study presents a baseline description of the prevalence, serotypes, antimicrobial resistance phenotypes and genetic relatedness of Salmonella isolated from healthy and diarrhoeic humans, and cattle at harvest, on farm and at market. Cattle are a reservoir of diverse salmonellae with shared serovars with humans, but WGST does not support zoonotic transmission. Further study with larger samples is recommended to determine whether epidemiological link exists between cattle and humans.

Bacteriophages represent an alternative solution to control bacterial infections. When interacting, bacteria and phage can evolve, and this relationship is described as antagonistic coevolution, a pattern that does not fit all models. In this work, the model consisted of a microcosm of Salmonella enterica serovar Enteritidis and &phi, San23 phage. Samples were taken for 12 days every 48 h. Bacteria and phage samples were collected, and isolated bacteria from each time point were challenged against phages from previous, contemporary, and subsequent time points. The phage plaque tests, with the genomics analyses, showed a mutational asymmetry dynamic in favor of the bacteria instead of antagonistic coevolution. This is important for future phage-therapy applications, so we decided to explore the population dynamics of Salmonella under different conditions: pressure of one phage, a combination of phages, and phages plus an antibiotic. The data from cultures with single and multiple phages, and antibiotics, were used to create a mathematical model exploring population and resistance dynamics of Salmonella under these treatments, suggesting a nonlethal, growth-inhibiting antibiotic may decrease resistance to phage-therapy cocktails. These data provide a deep insight into bacterial dynamics under different conditions and serve as additional criteria to select phages and antibiotics for phage-therapy.

Salmonella enterica subspecies enterica serotype Choleraesuis is a swine adapted serovar. S. Choleraesuis variant Kunzendorf is responsible for the majority of outbreaks among pigs. S. Choleraesuis is rare in Europe, although there have been serious outbreaks in pigs including two outbreaks in Denmark in 1999-2000 and 2012-2013. Here, we elucidate the epidemiology, possible transmission routes and sources, and clonality of European S. Choleraesuis isolates including the Danish outbreak isolates. A total of 102 S. Choleraesuis isolates from different European countries and the United States, covering available isolates from the last two decades were selected for whole genome sequencing. We applied a temporally structured sequence analysis within a Bayesian framework to reconstruct a temporal and spatial phylogenetic tree. MLST type, resistance genes, plasmid replicons, and accessory genes were identified using bioinformatics tools. Fifty-eight isolates including 11 out of 12 strains from wild boars were pan-susceptible. The remaining isolates carried multiple resistance genes. Eleven different plasmid replicons in eight plasmids were determined among the isolates. Accessory genes were associated to the identified resistance genes and plasmids. The European S. Choleraesuis was estimated to have emerged in ∼ 1837 (95% credible interval, 1733-1983) with the mutation rate of 1.02 SNPs/genome/year. The isolates were clustered according to countries and neighbor countries. There were transmission events between strains from the United States and European countries. Wild boar and pig isolates were genetically linked suggesting cross-border transmission and transmission due to a wildlife reservoir. The phylogenetic tree shows that multiple introductions were responsible for the outbreak of 2012-2013 in Denmark, and suggests that poorly disinfected vehicles crossing the border into Denmark were potentially the source of the outbreak. Low levels of single nucleotide polymorphisms (SNPs) differences (0-4 SNPs) can be observed between clonal strains isolated from different organs of the same animal. Proper disinfection of livestock vehicles and improved quality control of livestock feed could help to prevent future spread of S. Choleraesuis or other more serious infectious diseases such as African swine fever (ASF) in the European pig production system.

The occurrence of toxigenic Vibrio cholerae O1 during a non- outbreak period in Lake Victoria was studied and genetic characteristics for environmental persistence and relatedness to pandemic strains were assessed. We analyzed 360 samples of carps, phytoplankton and water collected in 2017 during dry and rainy seasons in the Tanzanian basin of Lake Victoria. Samples were tested using PCR (ompW and ctxA) with DNA extracted from bacterial isolates and samples enriched in alkaline peptone water. Isolates were screened with polyvalent antiserum O1 followed by antimicrobial susceptibility testing. Whole genome sequencing and bioinformatics tools were employed to investigate the genomic characteristics of the isolates. More V. cholerae positive samples were recovered by PCR when DNA was obtained from enriched samples than from isolates (69.0% vs. 21.3%, p < 0.05), irrespectively of season. We identified ten V. cholerae O1 among 22 ctxA-positive isolates. Further studies are needed to serotype the remaining ctxA-positive non-O1 strains. Sequenced strains belonged to El Tor atypical biotype of V. cholerae O1 of MLST ST69 harboring the seventh pandemic gene. Major virulence genes, ctxA, ctxB, zot, ace, tcpA, hlyA, rtxA, ompU, toxR, T6SS, alsD, makA and pathogenicity islands VPI-1, VPI-2, VSP-1, and VSP-2 were found in all strains. The strains contained Vibrio polysaccharide biosynthesis enzymes, the mshA gene and two-component response regulator proteins involved in stress response and autoinducers for quorum sensing and biofilm formation. They carried the SXT integrative conjugative element with phenotypic and genotypic resistance to aminoglycoside, sulfamethoxazole, trimethoprim, phenicol, and quinolones. Strains contained a multidrug efflux pump component and were resistant to toxic compounds with copper homeostasis and cobalt-zinc-cadmium resistance proteins. The environmental strains belonged to the third wave of the seventh pandemic and most are genetically closely related to recent outbreak strains from Tanzania, Kenya, and Uganda with as low as three SNPs difference. Some strains have persisted longer in the environment and were more related to older outbreak strains in the region. V. cholerae O1 of outbreak potential seem to persist in Lake Victoria through interactions with fish and phytoplankton supported by the optimum water parameters and intrinsic genetic features enhancing survival in the aquatic environment.

Whole-genome sequencing (WGS) has become a significant tool in investigating foodborne disease outbreaks and some countries have incorporated WGS into national food control systems. However, WGS poses technical challenges that deter developing countries from incorporating it into their food safety management system. A rapid scoping review was conducted, followed by a focus group session, to understand the current situation regarding the use of WGS for foodborne disease surveillance and food monitoring at the global level and identify key limiting factors for developing countries in adopting WGS for their food control systems. The results showed that some developed nations routinely use WGS in their food surveillance systems resulting in more precise understanding of the causes of outbreaks. In developing nations, knowledge of WGS exists in the academic/research sectors; however, there is limited understanding at the government level regarding the usefulness of WGS for food safety regulatory activities. Thus, incorporation of WGS is extremely limited in most developing nations. While some countries lack the capacity to collect and analyze the data generated from WGS, the most significant technical gap in most developing countries is in data interpretation using bioinformatics. The gaps in knowledge and capacities between developed and developing nations regarding use of WGS likely introduce an inequality in international food trade, and thus, relevant international organizations, as well as the countries that are already proficient in the use of WGS, have significant roles in assisting developing nations to be able to fully benefit from the technology and its applications in food safety management.

Next‐generation sequencing (NGS) data present an untapped potential to improve microbial risk assessment (MRA) through increased specificity and redefinition of the hazard. Most of the MRA models do not account for differences in survivability and virulence among strains. The potential of machine learning algorithms for predicting the risk/health burden at the population level while inputting large and complex NGS data was explored with Listeria monocytogenes as a case study. Listeria data consisted of a percentage similarity matrix from genome assemblies of 38 and 207 strains of clinical and food origin, respectively. Basic Local Alignment (BLAST) was used to align the assemblies against a database of 136 virulence and stress resistance genes. The outcome variable was frequency of illness, which is the percentage of reported cases associated with each strain. These frequency data were discretized into seven ordinal outcome categories and used for supervised machine learning and model selection from five ensemble algorithms. There was no significant difference in accuracy between the models, and support vector machine with linear kernel was chosen for further inference (accuracy of 89% [95% CI: 68%, 97%]). The virulence genes FAM002725, FAM002728, FAM002729, InlF, InlJ, Inlk, IisY, IisD, IisX, IisH, IisB, lmo2026, and FAM003296 were important predictors of higher frequency of illness. InlF was uniquely truncated in the sequence type 121 strains. Most important risk predictor genes occurred at highest prevalence among strains from ready‐to‐eat, dairy, and composite foods. We foresee that the findings and approaches described offer the potential for rethinking the current approaches in MRA.

BACKGROUND: Reliable phenotypic antimicrobial susceptibility testing can be a challenge in clinical settings in low- and middle-income countries. WGS is a promising approach to enhance current capabilities.AIM: To study diversity and resistance determinants and to predict and compare resistance patterns from WGS data of Acinetobacter baumannii with phenotypic results from classical microbiological testing at a tertiary care hospital in Tanzania.METHODS AND RESULTS: MLST using Pasteur/Oxford schemes yielded eight different STs from each scheme. Of the eight, two STs were identified to be global clones 1 (n = 4) and 2 (n = 1) as per the Pasteur scheme. Resistance testing using classical microbiology determined between 50% and 92.9% resistance across all drugs. Percentage agreement between phenotypic and genotypic prediction of resistance ranged between 57.1% and 100%, with coefficient of agreement (κ) between 0.05 and 1. Seven isolates harboured mutations at significant loci (S81L in gyrA and S84L in parC). A number of novel plasmids were detected, including pKCRI-309C-1 (219000 bp) carrying 10 resistance genes, pKCRI-43-1 (34 935 bp) carrying two resistance genes and pKCRI-49-1 (11 681 bp) and pKCRI-28-1 (29 606 bp), each carrying three resistance genes. New ampC alleles detected included ampC-69, ampC-70 and ampC-71. Global clone 1 and 2 isolates were found to harbour ISAba1 directly upstream of the ampC gene. Finally, SNP-based phylogenetic analysis of the A. baumannii isolates revealed closely related isolates in three clusters.CONCLUSIONS: The validity of the use of WGS in the prediction of phenotypic resistance can be appreciated, but at this stage is not sufficient for it to replace conventional antimicrobial susceptibility testing in our setting.

Bacteriophages represent an alternative solution to control bacterial infections. When interacting, bacteria and phage can evolve, and this relationship is described as antagonistic coevolution, a pattern that does not fit all models. In this work, the model consisted of a microcosm of Salmonella enterica serovar Enteritidis and φSan23 phage. Samples were taken for 12 days every 48 h. Bacteria and phage samples were collected; and isolated bacteria from each time point were challenged against phages from previous, contemporary, and subsequent time points. The phage plaque tests, with the genomics analyses, showed a mutational asymmetry dynamic in favor of the bacteria instead of antagonistic coevolution. This is important for future phage-therapy applications, so we decided to explore the population dynamics of Salmonella under different conditions: pressure of one phage, a combination of phages, and phages plus an antibiotic. The data from cultures with single and multiple phages, and antibiotics, were used to create a mathematical model exploring population and resistance dynamics of Salmonella under these treatments, suggesting a nonlethal, growth-inhibiting antibiotic may decrease resistance to phage-therapy cocktails. These data provide a deep insight into bacterial dynamics under different conditions and serve as additional criteria to select phages and antibiotics for phage-therapy.

Wastewaters serve as important hot spots for antimicrobial resistance and monitoring can be used to analyse the abundance and diversity of antimicrobial resistance genes at the level of large bacterial and human populations. In this study, whole genome sequencing of beta-lactamase-producing Escherichia coli and metagenomic analysis of whole-community DNA were used to characterize the occurrence of antimicrobial resistance in hospital, municipal and river waters in the city of Brno (Czech Republic). Cefotaxime-resistant E. coli were mainly extended-spectrum beta-lactamase (ESBL) producers (95.6%, n = 158), of which the majority carried blaCTX-M (98.7%; n = 151) and were detected in all water samples except the outflow from hospital wastewater treatment plant. A wide phylogenetic diversity was observed among the sequenced E. coli (n = 78) based on the detection of 40 sequence types and single nucleotide polymorphisms (average number 34,666 ± 15,710) between strains. The metagenomic analysis revealed a high occurrence of bacterial genera with potentially pathogenic members, including Pseudomonas, Escherichia, Klebsiella, Aeromonas, Enterobacter and Arcobacter (relative abundance >50%) in untreated hospital and municipal wastewaters and predominance of environmental bacteria in treated and river waters. Genes encoding resistance to aminoglycosides, beta-lactams, quinolones and macrolides were frequently detected, however blaCTX-M was not found in this dataset which may be affected by insufficient sequencing depth of the samples. The study pointed out municipal treated wastewater as a possible source of multi-drug resistant E. coli and antimicrobial resistance genes for surface waters. Moreover, the combination of two different approaches provided a more holistic view on antimicrobial resistance in water environments. The culture-based approach facilitated insight into the dynamics of ESBL-producing E. coli and the metagenomics shows abundance and diversity of bacteria and antimicrobial resistance genes vary across water sites.

Listeria monocytogenes is a ubiquitous bacterium that may cause the foodborne illness listeriosis. Only a small amount of data about the population genetic structure of strains isolated from food is available. This study aimed to provide an accurate view of the L. monocytogenes food strain population in France. From 1999 to 2014, 1,894 L. monocytogenes strains were isolated from food at the French National Reference Laboratory for L. monocytogenes and classified according to the five risk food matrices defined by the European Food Safety Authority (EFSA). A total of 396 strains were selected on the basis of different pulsed-field gel electrophoresis (PFGE) clusters, serotypes, and strain origins and typed by multilocus sequence typing (MLST), and the MLST results were supplemented with MLST data available from Institut Pasteur, representing human and additional food strains from France. The distribution of sequence types (STs) was compared between food and clinical strains on a panel of 675 strains. High congruence between PFGE and MLST was found. Out of 73 PFGE clusters, the two most prevalent corresponded to ST9 and ST121. Using original statistical analysis, we demonstrated that (i) there was not a clear association between ST9 and ST121 and the food matrices, (ii) serotype IIc, ST8, and ST4 were associated with meat products, and (iii) ST13 was associated with dairy products. Of the two major STs, ST121 was the ST that included the fewest clinical strains, which might indicate lower virulence. This observation may be directly relevant for refining risk analysis models for the better management of food safety. IMPORTANCE This study showed a very useful backward compatibility between PFGE and MLST for surveillance. The results enabled better understanding of the population structure of L. monocytogenes strains isolated from food and management of the health risks associated with L. monocytogenes food strains. Moreover, this work provided an accurate view of L. monocytogenes strain populations associated with specific food matrices. We clearly showed that some STs were associated with food matrices, such as meat, meat products, and dairy products. We opened the way to source attribution modeling in order to quantify the relative importance of the main food matrices.

A Corrigendum onColistin Resistance Mediated by mcr-1 in ESBL-Producing, Multidrug Resistant Salmonella Infantis in Broiler Chicken Industry, Italy (2016–2017)by Carfora, V., Alba, P., Leekitcharoenphon, P., Ballarò, D., Cordaro, G., Di Matteo, P., et al. (2018) Front. Microbiol. 9:1880. doi: 10.3389/fmicb.2018.01880In the original article, there was an error in the Materials and Methods, subsection Isolates. The four S. Infantis isolates originated from broilers (n = 2) and broiler meat samples (n = 2).A correction has been made to Materials and Methods, subsection Isolates: Four multidrug resistant (MDR) S. Infantis, displaying a colistin MIC value ≥ 4 mg/L, were detected among 324 S. Infantis isolates collected in the frame of antimicrobial resistance (AMR) monitoring activities conducted from 2001 to 2017 by the National Reference Laboratory for Antimicrobial Resistance (NRL-AR) and screened for antimicrobial susceptibility. The four S. Infantis isolates originated from broilers (n = 3) and from broiler meat sample (n = 1) (Supplementary Table 1).The authors apologize for this error and state that this does not change the scientific conclusionsof the article in any way.The original article has been updated.

Objectives This study compares the genome of an ST131 CMY-2-producing Escherichia coli isolate from a Danish patient with other ST131 CMY-2-producing E. coli isolates of both human and animal origin. Methods In 2016, an ST131 CMY-2-producing E. coli isolate (ESBL20160056) was obtained from a patient with a bloodstream infection. The genome of the ESBL20160056 isolate was compared with genomes from six ST131 CMY-2-producing E. coli isolates obtained from broiler meat imported to Denmark, 15 ST131 CMY-2-producing E. coli isolates obtained from Enterobase (http://enterobase.warwick.ac.uk) and two ST131 CMY-2-producing E. coli from European collaborators. The plasmid from ESBL20160056 was sequenced using a MinION Mk1B (Oxford Nanopore Technologies). Results The E. coli isolate from the Danish patient clustered together with 13 other fimH22 ST131 CMY-2-producing E. coli isolates in a distinct clade. The clade consisted of genomes from six E. coli isolates from humans collected in Denmark, Spain, Cambodia and the USA, six E. coli isolates obtained from broiler meat samples imported to Denmark from France, the Netherlands and Germany, and two E. coli isolates obtained from broilers in Belgium and Luxembourg. The 101.5 kb plasmid with blaCMY-2 from ESBL20160056 had an IncI1 replicon and belonged to ST12 using the plasmid MLST scheme. In total, 10 of the 14 ST131 E. coli isolates belonging to the fimH22 clade carried an IncI1 ST12 plasmid with blaCMY-2. Conclusions From our data, it seems plausible that the ST131 fimH22 CMY-2-producing E. coli isolate obtained from the Danish patient could have a zoonotic broiler origin.

Colistin resistance by mobilisable mcr genes has been described in bacteria of food-animal origin worldwide, which has raised public health concerns about its potential foodborne transmission to human pathogenic bacteria. Here we provide baseline information on the molecular epidemiology of colistin-resistant, mcr-positive Escherichia coli and Salmonella isolates in food-producing animals in Italy in 2014-2015. A total 678, 861 and 236 indicator E. coli, Extended Spectrum Beta-Lactamase (ESBL)/AmpC-producing E. coli, and Salmonella isolates, respectively, were tested for colistin susceptibility. These isolates were collected according to the EU harmonized antimicrobial resistance monitoring program and are representative of at least 90 and 80% of the Italian poultry (broiler chickens and turkeys) and livestock (pigs and bovines < 12 months) production, respectively. Whole genome sequencing by Illumina technology and bioinformatics (Center for Genomic Epidemiology pipeline) were used to type 42 mcr-positive isolates by PCR. Colistin resistance was mainly observed in the ESBL/AmpC E. coli population, and was present in 25.9, 5.3, 6.5, and 3.9% of such isolates in turkeys, broilers, pigs, and bovines, respectively. Most colistin-resistant isolates (141/161, 87.5%) harbored genes of the mcr-1 group. mcr-1 was also detected in a small proportion of Salmonella isolates (3/146, 2.0%) in turkeys. Additional mcr types were mcr-3 in four ESBL-producing E. coli from bovines, and two mcr-4 in ESBL (n = 1) and indicator E. coli (n = 1) from pigs and bovines. We describe notable diversity of mcr variants with predominance of mcr-1.1 and mcr-1.2 on conjugative IncX4 plasmids in E. coli and in Salmonella serovars Typhimurium, Newport, Blockley from turkey. A new variant, mcr-1.13 was detected in the chromosome in E. coli in turkey and pig isolates. Additionally, we describe mcr-3.2 and mcr-4.3 in E. coli from bovines, and mcr-4.2 in E. coli from pigs. These findings elucidate the epidemiology of colistin resistance in food-producing animals in Italy along with its genetic background, and highlight the likelihood of mcr horizontal transfer between commensal bacteria and major food-borne pathogens (Salmonella) within the same type of productions. Thorough action and strategies are needed in order to mitigate the risk of mcr transfer to humans, in a “One Health” perspective.

Whole genome sequencing was utilized to investigate the genomic profile of Vibrio cholerae O1 strains, isolated from symptomatic patients in a low-income urban area of Dhaka, Bangladesh. Comparative genomics using bioinformatics tools were applied to identify major virulence factors, biotype and antimicrobial resistance genes in three V. cholerae O1 strains (VC-1, 2 and 3) isolated from two case patients. A phylogenetic SNP (single nucleotide polymorphism)-based analysis was conducted to infer the relatedness to V. cholerae O1 strains isolated elsewhere. The V. cholerae strains were the El Tor variant carrying ctxB1 (standard classical genotype). SNP-based global phylogeny revealed that the three isolates were strictly clonal and the closest neighbouring genomes were epidemic clones of V. cholerae O1 isolated in 2010 from cholera patients in Pakistan. All strains harboured the integrase gene of the SXT element (intSXT ), antimicrobial resistance genes for aminoglycosides, phenicol, sulphonamide and trimethoprim except VC-1 that lacked sulphonamide resistance genes. The multilocus sequence typing (MLST) revealed that the strains belonged to sequence type, ST69. The study provides knowledge on current genetic traits of clinical V. cholerae O1 circulating in urban household clusters of Bangladesh which may help in predicting emergence of new pandemic strains in Bangladesh. SIGNIFICANCE AND IMPACT OF THE STUDY Vibrio cholerae has frequently experienced genetic changes with rapid evolution of pandemic clones in the Ganges Delta region. Whole genome sequencing can reveal genetic information of current pathogenic V. cholerae in Bangladesh which includes cefotaxime genotypes, virulence factors, altered antimicrobial resistance pattern as well as mobile genetic element compared to global pandemic strains. This study data could be used in planning future surveillance strategies in Ganges Delta region by informing new epidemiology of current outbreak strains.

A total of 502 Campylobacter jejuni isolates from poultry in 12 different European countries (10 of them the largest poultry production countries in Europe) were whole genome sequenced to examine the genomic diversity of fluoroquinolone resistant (FQ‐R) and susceptible (FQ‐S) C. jejuni across the poultry producing European countries and to determine whether the emergence of fluoroquinolone resistance among C. jejuni is related to the transmission through countries or to the selection through fluoroquinolone use in the individual countries. A high genomic diversity was observed. The isolates clustered in four main clusters. All trees revealed that the isolates were clustered according to the presence/absence of the gyrA mutations causing fluoroquinolone resistance and ST‐types. The cgMLST trees of only FQ‐R and FQ‐S isolates showed that isolates from the same country of origin were distributed into multiple clusters similarly to the trees combining FQ‐R and FQ‐S isolates. The different phylogenetic methods, ranging from single nucleotide polymorphisms analysis to gene‐by‐gene approaches such as rMLST, cgMLST, wgMLST and core genome tree, provided concordant results, but it is not known which is the most accurate method for identifying the country of origin of the isolates. Allele frequency analysis of isolates under this study and a selection of previously published C. jejuni genomes in ENA showed association of geographical origin of poultry C. jejuni populations between Romania‐Poland, Italy‐Germany‐England, Portugal‐The Netherlands and USA‐Luxemburg. Allele frequency and phylogenetic analysis indicated that the isolates from Finland were genetically different from C. jejuni populations from other European countries included in this study. Trade pattern and antimicrobial use in livestock were not significantly associated with allele frequency or populations of C. jejuni, but data available to investigate these associations were limited.

Salmonella enterica subsp. enterica serotype Typhimurium sequence type 313 (ST313) is most commonly associated with invasive nontyphoidal Salmonella disease in Africa among patients with HIV infection and malignancy., Salmonella enterica subsp. enterica serotype Typhimurium sequence type 313 (ST313) is most commonly associated with invasive nontyphoidal Salmonella disease in Africa among patients with HIV infection and malignancy. Here, we report a draft genome sequence of S. Typhimurium ST313, isolated from an elderly immunosuppressed patient from India with non-Hodgkins lymphoma.

The ENGAGE project (http://www.engage-europe.eu/) was a collaboration between eight institutions across Europe. The aim was to boost the scientific cooperation to use whole genome sequencing (WGS) analysis in food safety and public health protection. ENGAGE focused on Escherichia coli (commensal E. coli) and different Salmonella spp. serotypes. A total of 3,360 genomes, 778 and 2,582 of E. coli and Salmonella, respectively, were produced. These genomes were stored and shared among partners in a temporary repository to be submitted to the European Nucleotide Archive by the end of the project. Generated genomes were used for benchmarking exercises to assess the possibility of replacing conventional typing with WGS for outbreak investigation. For the analysed strains, the benchmarking exercises showed that SPAdes assembly performed better than Velvet and that, by using different bioinformatics tools, WGS Salmonella serotyping and antimicrobial resistance genes detection, were largely in concordance with phenotypic data. Discrepancies were related to sequence quality and phenotype misclassification rather than to limitations of the bioinformatics tools. All partners were able to infer the expected phylogeny for the Salmonella and Campylobacter isolates in benchmarking exercises. Two WGS proficiency tests (assessing different genomic quality markers) were conducted among partners with satisfactory results. Guidelines including available bioinformatics tools and standard operating procedures (wet and dry lab) were prepared and posted online. Workshops, training courses and twinning programmes were conducted. The training focused on online, Galaxy‐based, and command line bioinformatics tools. To reach out beyond ENGAGE, an e‐learning course (17 videos) was developed and made available online. Several proof of concept projects were run and some outcomes published, e.g. the discovery of colistin resistance gene, mcr‐5. Overall, the project showed that laboratories without previous WGS experience need a period of time to implement and perform WGS for foodborne pathogens routine analysis. All developed material will remain available on the ENGAGE website.

Objective. To determine molecular epidemiology of methicillin-resistantS. aureusin Tanzania using whole genome sequencing.Methods. DNA from 33Staphylococcusspecies was recovered from subcultured archivedStaphylococcusisolates. Whole genome sequencing was performed on Illumina Miseq using paired-end2×250 bp protocol. Raw sequence data were analyzed using online tools.Results. Full susceptibility to vancomycin and chloramphenicol was observed. Thirteen isolates (43.3%) resisted cefoxitin and other antimicrobials tested. Multilocus sequence typing revealed 13 different sequence types among the 30S. aureusisolates, with ST-8 (n= seven, 23%) being the most common. Gene detection inS.aureusstains were as follows:mecA, 10 (33.3%);pvl, 5 (16.7%);tst, 2 (6.7%). The SNP difference among the six Tanzanian ST-8 MRSA isolates ranged from 24 to 196 SNPs and from 16 to 446 SNPs when using the USA300_FPR3757 or the USA500_2395 as a reference, respectively. The mutation rate was1.38×10-11SNPs/site/year or1.4×10-6SNPs/site/year as estimated by USA300_FPR3757 or the USA500_2395, respectively.Conclusion.S. aureusisolates causing infections in hospitalized patients in Moshi are highly diverse and epidemiologically unrelated. Temporal phylogenetic analysis provided better resolution on transmission and introduction of MRSA and it may be important to include this in future routines.

Furunculosis, a serious infection caused by the bacterium Aeromonas salmonicida subsp. salmonicida is common in sea-reared rainbow trout production in Denmark. Developing an effective control strategy requires knowledge of the epidemiology, as well as the genomic and virulent variability of the Danish A. salmonicida subsp. salmonicida isolates. To obtain this, the genomes of 101 A. salmonicida subsp. salmonicida, including 99 Danish isolates, one Scottish strain and the type strain NCIMB 1102, were sequenced using the Illumina HiSeq platform. Isolates were de novo assembled, examined for presence of plasmids, virulence and iron acquisition proteins, genomic islands, and antibiotic resistance genes. Single Nucleotide Polymorphisms were aligned and subjected to Bayesian temporal phylogenetic and maximum likelihood tree reconstruction using the published genome of A. salmonicida subsp. salmonicida A449 as reference. Bayesian temporal phylogenetic reconstruction suggests that four major introductions of A. salmonicida subsp. salmonicida into Denmark have occurred. The introductions correlate with the freshwater and subsequent seawater expansion of rainbow trout production. Initial transmission of the bacterium could have been from seawater to freshwater or vice versa, and most minor clades include a mixture of strains from different fresh- and seawater farms. Genomic variation of A. salmonicida subsp. salmonicida mostly appeared to be associated with their plasmids and plasmid encoded virulence factors. Nine A. salmonicida subsp. salmonicida isolates harbored worldwide known antibiotic resistance genes against several antibiotics and there is an indication that 33% of the isolates contained the genomic island AsaGEI1b. These findings not only support the usefulness of whole genome sequencing for genetic studies of homogeneous bacteria in general, but provide novel information about the Danish A. salmonicida subsp. salmonicida population, with implications for vaccine development in efforts to better protect Danish rainbow trout in the future.

Vibrio harveyi CAIM 1792 is a marine bacterial strain that causes mortality in farmed shrimp in north-west Mexico, and the identification of virulence genes in this strain is important for understanding its pathogenicity. The aim of this work was to compare the V. harveyi CAIM 1792 genome with related genome sequences to determine their phylogenic relationship and explore unique regions in silico that differentiate this strain from other V. harveyi strains. Twenty-one newly sequenced genomes were compared in silico against the CAIM 1792 genome at nucleotidic and predicted proteome levels. The proteome of CAIM 1792 had higher similarity to those of other V. harveyi strains (78%) than to those of the other closely related species Vibrio owensii (67%), Vibrio rotiferianus (63%) and Vibrio campbellii (59%). Pan-genome ORFans trees showed the best fit with the accepted phylogeny based on DNA-DNA hybridization and multi-locus sequence analysis of 11 concatenated housekeeping genes. SNP analysis clustered 34/38 genomes within their accepted species. The pangenomic and SNP trees showed that V. harveyi is the most conserved of the four species studied and V. campbellii may be divided into at least three subspecies, supported by intergenomic distance analysis. blastp atlases were created to identify unique regions among the genomes most related to V. harveyi CAIM 1792; these regions included genes encoding glycosyltransferases, specific type restriction modification systems and a transcriptional regulator, LysR, reported to be involved in virulence, metabolism, quorum sensing and motility.

Salmonella enterica serovar Choleraesuis is a porcine adapted serovar which may cause serious outbreaks in pigs. Here we describe outbreaks of salmonellosis due to S. Choleraesuis in four Danish pig farms in 2012-2013 by clinic, serology, and microbiology and compare the isolates to those of a previous outbreak in 1999-2000. The infection was in some herds associated with high mortality and a moderate to high sero-prevalence was found. In 2012-2013 the disease contributed to increased mortality but occurred concomitant with other disease problems in the herds, which likely delayed the diagnosis by up to several months. Nine isolates from the four farms in 2012-2013 and 14 isolates obtained from the outbreak in Denmark in 1999-2000 were subjected to typing using pulsed-field gel electrophoresis (PFGE). Seven isolates were selected for whole genome sequencing (WGS). The PFGE results of 23 isolates displayed five different profiles. The isolates from 2012 to 2013 revealed two distinct profiles, both different from the isolates recovered in 1999-2000. Two of the 2012-2013 farms shared PFGE profiles and had also transported pigs between them. The profile found in the two other 2012-2013 farms was indistinguishable but no epidemiological connection between these farms was found. Analysis of the number of single nucleotide polymorphisms (SNPs) from the WGS data indicated that the isolates from the farms in 2012-2013 were more closely related to each other than to isolates from the outbreak in 1999. It was therefore concluded that the infection was a new introduction and not a persistent infection since the outbreak in 1999. It may further be suggested that there were two or three independent rather than a single introduction. The re-introduction of S. Choleraesuis in Denmark emphasizes the importance of strict hygiene measures in the herds. Further investigations using WGS are now in progress on a larger collection of isolates to study clonality at European level and trace the origin of the infections.

Furunculosis, a serious infection caused by the bacterium Aeromonas salmonicida subsp. salmonicida is common in sea-reared rainbow trout production in Denmark. Developing an effective control strategy requires knowledge of the epidemiology, as well as the genomic and virulent variability of the Danish A. salmonicida subsp. salmonicida isolates. To obtain this, the genomes of 101 A. salmonicida subsp. salmonicida, including 99 Danish isolates, one Scottish strain and the type strain NCIMB 1102, were sequenced using the Illumina HiSeq platform. Isolates were de novo assembled, examined for presence of plasmids, virulence and iron acquisition proteins, genomic islands, and antibiotic resistance genes. Single Nucleotide Polymorphisms were aligned and subjected to Bayesian temporal phylogenetic and maximum likelihood tree reconstruction using the published genome of A. salmonicida subsp. salmonicida A449 as reference. Bayesian temporal phylogenetic reconstruction suggests that four major introductions of A. salmonicida subsp. salmonicida into Denmark have occurred. The introductions correlate with the freshwater and subsequent seawater expansion of rainbow trout production. Initial transmission of the bacterium could have been from seawater to freshwater or vice versa, and most minor clades include a mixture of strains from different fresh- and seawater farms. Genomic variation of A. salmonicida subsp. salmonicida mostly appeared to be associated with their plasmids and plasmid encoded virulence factors. Nine A. salmonicida subsp. salmonicida isolates harbored worldwide known antibiotic resistance genes against several antibiotics and there is an indication that 33% of the isolates contained the genomic island AsaGEI1b. These findings not only support the usefulness of whole genome sequencing for genetic studies of homogeneous bacteria in general, but provide novel information about the Danish A. salmonicida subsp. salmonicida population, with implications for vaccine development in efforts to better protect Danish rainbow trout in the future.

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. Electronic supplementary material The online version of this article (doi:10.1186/s12879-017-2628-x) contains supplementary material, which is available to authorized users.

Objective To determine molecular epidemiology of methicillin-resistant S. aureus in Tanzania using whole genome sequencing. Methods DNA from 33 Staphylococcus species was recovered from subcultured archived Staphylococcus isolates. Whole genome sequencing was performed on Illumina Miseq using paired-end 2 × 250 bp protocol. Raw sequence data were analyzed using online tools. Results Full susceptibility to vancomycin and chloramphenicol was observed. Thirteen isolates (43.3%) resisted cefoxitin and other antimicrobials tested. Multilocus sequence typing revealed 13 different sequence types among the 30 S. aureus isolates, with ST-8 (n = seven, 23%) being the most common. Gene detection in S. aureus stains were as follows: mecA, 10 (33.3%); pvl, 5 (16.7%); tst, 2 (6.7%). The SNP difference among the six Tanzanian ST-8 MRSA isolates ranged from 24 to 196 SNPs and from 16 to 446 SNPs when using the USA300_FPR3757 or the USA500_2395 as a reference, respectively. The mutation rate was 1.38 × 10−11 SNPs/site/year or 1.4 × 10−6 SNPs/site/year as estimated by USA300_FPR3757 or the USA500_2395, respectively. Conclusion S. aureus isolates causing infections in hospitalized patients in Moshi are highly diverse and epidemiologically unrelated. Temporal phylogenetic analysis provided better resolution on transmission and introduction of MRSA and it may be important to include this in future routines.

Acinetobacter johnsonii C6 originates from creosote-polluted groundwater and performs ecological and evolutionary interactions with Pseudomonas putida in biofilms. The draft genome of A. johnsonii C6 is 3.7 Mbp and was shaped by mobile genetic elements. It reveals genes facilitating the biodegradation of aromatic hydrocarbons and resistance to antimicrobials and metals.

Background/ objectives: Whole genome sequencing (WGS) has proven to be a powerful subtyping tool for foodborne pathogenic bacteria like L. monocytogenes. The interests of genome-scale analysis for national surveillance, outbreak detection or source tracking has been largely documented. The genomic data however can be exploited with many different bioinformatics methods like single nucleotide polymorphism (SNP), core-genome multi locus sequence typing (cgMLST), whole-genome multi locus sequence typing (wgMLST) or multi locus predicted protein sequence typing (MLPPST) on either core-genome (cgMLPPST) or pan-genome (wgMLPPST). Currently, there are little comparisons studies of these different analytical approaches. Our objective was to assess and compare different genomic methods that can be implemented in order to cluster isolates of L. monocytogenes. Methods: The clustering methods were evaluated on a collection of 207 L. monocytogenes genomes of food origin representative of the genetic diversity of the Anses collection. The trees were then compared using robust statistical analyses. Results: The backward comparability between conventional typing methods and genomic methods revealed a near-perfect concordance. The importance of selecting a proper reference when calling SNPs was highlighted, although distances between strains remained identical. The analysis also revealed that the topology of the phylogenetic trees between wgMLST and cgMLST were remarkably similar. The comparison between SNP and cgMLST or SNP and wgMLST approaches showed that the topologies of phylogenic trees were statistically similar with an almost equivalent clustering. Conclusion: Our study revealed high concordance between wgMLST, cgMLST, and SNP approaches which are all suitable for typing of L. monocytogenes. The comparable clustering is an important observation considering that the two approaches have been variously implemented among reference laboratories.

Cholera is still an important public health problem in several countries, including Thailand. In this study, a collection of clinical and environmental V. cholerae serogroup O1, O139, and non-O1/non-O139 strains originating from Thailand (1983 to 2013) was characterized to determine phenotypic and genotypic traits and to investigate the genetic relatedness. Using a combination of conventional methods and whole genome sequencing (WGS), 78 V. cholerae strains were identified. WGS was used to determine the serogroup, biotype, virulence, mobile genetic elements, and antimicrobial resistance genes using online bioinformatics tools. In addition, phenotypic antimicrobial resistance was determined by the minimal inhibitory concentration (MIC) test. The 78 V. cholerae strains belonged to the following serogroups O1: (n = 44), O139 (n = 16) and non-O1/non-O139 (n = 18). Interestingly, we found that the typical El Tor O1 strains were the major cause of clinical cholera during 1983-2000 with two Classical O1 strains detected in 2000. In 2004-2010, the El Tor variant strains revealed genotypes of the Classical biotype possessing either only ctxB or both ctxB and rstR while they harbored tcpA of the El Tor biotype. Thirty O1 and eleven O139 clinical strains carried CTXφ (Cholera toxin) and tcpA as well four different pathogenic islands (PAIs). Beside non-O1/non-O139, the O1 environmental strains also presented chxA and Type Three Secretion System (TTSS). The in silico MultiLocus Sequence Typing (MLST) discriminated the O1 and O139 clinical strains from other serogroups and environmental strains. ST69 was dominant in the clinical strains belonging to the 7th pandemic clone. Non-O1/non-O139 and environmental strains showed various novel STs indicating genetic variation. Multidrug-resistant (MDR) strains were observed and conferred resistance to ampicillin, azithromycin, nalidixic acid, sulfamethoxazole, tetracycline, and trimethoprim and harboured variants of the SXT elements. For the first time since 1986, the presence of V. cholerae O1 Classical was reported causing cholera outbreaks in Thailand. In addition, we found that V. cholerae O1 El Tor variant and O139 were pre-dominating the pathogenic strains in Thailand. Using WGS and bioinformatic tools to analyze both historical and contemporary V. cholerae circulating in Thailand provided a more detailed understanding of the V. cholerae epidemiology, which ultimately could be applied for control measures and management of cholera in Thailand.

As whole genome sequence data of microorganisms are becoming easily accessible and cheap to produce, a transformation of the traditional methods used for typing, phenotyping and phylogenetic analysis of microorganisms is on the way. Following the anticipation that most clinical microbiological and food safety laboratories will soon have a sequencer in use on a daily basis, there is a growing need for easy-to-use bioinformatics methods that can quickly convert the sequence data into useful information on, e.g., the type of bacteria, whether it is resistant towards any types of antibiotics, and whether it is part of an outbreak. The Center for Genomic Epidemiology, which is located at the Technical University of Denmark, has since its beginning in 2010 developed such bioinformatics methods and made them freely available as web-services. These web-services and their use is the focus of this chapter.

Acinetobacter johnsonii C6 originates from creosote-polluted groundwater and performs ecological and evolutionary interactions with Pseudomonas putida in biofilms. The draft genome of A. johnsonii C6 is 3.7 Mbp and was shaped by mobile genetic elements. It reveals genes facilitating the biodegradation of aromatic hydrocarbons and resistance to antimicrobials and metals.

Comparison of the publicly available genomes of the virulentSalmonellaentericaserovar Typhimurium (S.Typhimurium) strains SL1344, 14028s and D23580 to that of the virulence-attenuated isolate LT2 revealed the absence of a full sequence of bacteriophage ST64B in the latter. Four selected ST64B regions of unknown function (sb7–sb11,sb46,sb49–sb50andsb54) were mapped by PCR in two strain collections: (i) 310 isolates ofS.Typhimurium from human blood or stool samples, and from food, animal and environmental reservoirs; and (ii) 90 isolates belonging to other serovars. The regionsb49–sb50was found to be unique toS.Typhimurium and was strongly associated with strains isolated from blood samples (100 and 28.4 % of the blood and non-blood isolates, respectively). The region was cloned into LT2 and knocked out in SL1344, and these strains were compared to wild-type isogenic strains inin vitroassays used to predict virulence association. No difference in invasion of the Int407 human cell line was observed between the wild-type and mutated strains, but the isolate carrying the whole ST64B prophage was found to have a slightly better survival in blood. The study showed a high prevalence and a strong association between the prophage ST64B and isolates ofS.Typhimurium collected from blood, and may indicate that such strains constitute a selected subpopulation within this serovar. Further studies are indicated to determine whether the slight increase in blood survival observed in the strain carrying ST64B genes is of paramount importance for systemic infections.