Search

Colistin is a last-resort drug for the treatment of infections by carbapenem-resistant Enterobacteriaceae, and the emergence of colistin resistance poses a serious clinical challenge. The aim of this study was to investigate the occurrence of colistin-resistant Escherichia coli in retail meat in Southern Italy in 2018-2020. Of 570 samples, 147 contained E. coli. Two out of 147 (1.4%) E. coli showed a non-wild-type phenotype to colistin and harboured mcr-1. mcr-1 was also detected in a wild-type isolate, resulting in a 2% mcr prevalence. mcr-1-positive isolates originated from turkey meat collected in Apulia (n = 2) and Basilicata (n = 1). A whole-genome sequencing analysis confirmed mcr-1.2 and mcr-1.1 in two and one isolate, respectively. The strains were diverse, belonging to three multi-locus sequence types (ST354, ST410, SLV of ST10) and harbouring genes mediating resistance to antimicrobials in two, six and seven classes. mcr-1 was carried by IncX4 plasmids with high nucleotide similarity to IncX4 plasmids harbouring mcr-1.2 and mcr-1.1 in Enterobacterales from different sources and geographical regions. This is the first study reporting updates on E. coli non-wild-type to colistin from retail meat in Southern Italy, highlighting the importance of phenotypic and genotypic antimicrobial resistance surveillance to contain the dissemination of mcr among E. coli.

Antimicrobial susceptibility testing (AST) should be fast and accurate, leading to proper interventions and therapeutic success. Clinical microbiology laboratories rely on phenotypic methods, but the continuous improvement and decrease in the cost of whole-genome sequencing (WGS) technologies make them an attractive alternative. Studies evaluating the performance of WGS-based prediction of antimicrobial resistance (AMR) for selected bacterial species have shown promising results. There are, however, significant gaps in the literature evaluating the applicability of WGS as a diagnostics method in real-life clinical settings against the range of bacterial pathogens experienced there. Thus, we compared standard phenotypic AST results with WGS-based predictions of AMR profiles in bacterial isolates without preselection of defined species, to evaluate the applicability of WGS as a diagnostics method in clinical settings. We collected all bacterial isolates processed by all Danish Clinical Microbiology Laboratories in 1 day. We randomly selected 500 isolates without any preselection of species. We performed AST through standard broth microdilution (BMD) for 488 isolates (n = 6,487 phenotypic AST results) and compared results with in silico antibiograms obtained through WGS (Illumina NextSeq) followed by bioinformatics analyses using ResFinder 4.0 (n = 5,229 comparisons). A higher proportion of AMR was observed for Gram-negative bacteria (10.9%) than for Gram-positive bacteria (6.1%). Comparison of BMD with WGS data yielded a concordance of 91.7%, with discordant results mainly due to phenotypically susceptible isolates harboring genetic AMR determinants. These cases correspond to 6.2% of all isolate-antimicrobial combinations analyzed and to 6.8% of all phenotypically susceptible combinations. We detected fewer cases of phenotypically resistant isolates without any known genetic resistance mechanism, particularly 2.1% of all combinations analyzed, which corresponded to 26.4% of all detected phenotypic resistances. Most discordances were observed for specific combinations of species-antimicrobial: macrolides and tetracycline in streptococci, ciprofloxacin and β-lactams in combination with β-lactamase inhibitors in Enterobacterales, and most antimicrobials in Pseudomonas aeruginosa. WGS has the potential to be used for surveillance and routine clinical microbiology. However, in clinical microbiology settings and especially for certain species and antimicrobial agent combinations, further developments in AMR gene databases are needed to ensure higher concordance between in silico predictions and expected phenotypic AMR profiles.

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.

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.

The application of high-throughput DNA sequencing technologies (WGS) data remain an increasingly discussed but vastly unexplored resource in the public health domain of quantitative microbial risk assessment (QMRA). This is due to challenges including high dimensionality of WGS data and heterogeneity of microbial growth phenotype data. This study provides an innovative approach for modeling the impact of population heterogeneity in microbial phenotypic stress response and integrates this into predictive models inputting a high-dimensional WGS data for increased precision exposure assessment using an example of Listeria monocytogenes. Finite mixture models were used to distinguish the number of sub-populations for each of the stress phenotypes, acid, cold, salt and desiccation. Machine learning predictive models were selected from six algorithms by inputting WGS data to predict the sub-population membership of new strains with unknown stress response data. An example QMRA was conducted for cultured milk products using the strains of unknown stress phenotype to illustrate the significance of the findings of this study. Increased resistance to stress conditions leads to increased growth, the likelihood of higher exposure and probability of illness. Neglecting within-species genetic and phenotypic heterogeneity in microbial stress response may over or underestimate microbial exposure and eventual risk during QMRA. ispartof: MICROORGANISMS vol:8 issue:11 ispartof: location:Switzerland status: published

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.

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

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.

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.

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.

Background Tanzania is one of seven countries with the highest disease burden caused by cholera in Africa. We studied the evolution of Vibrio cholerae O1 isolated in Tanzania during the past three decades. Methodology/Principal findings Genome-wide analysis was performed to characterize V. cholerae O1 responsible for the Tanzanian 2015–2017 outbreak along with strains causing outbreaks in the country for the past three decades. The genomes were further analyzed in a global context of 590 strains of the seventh cholera pandemic (7PET), as well as environmental isolates from Lake Victoria. All Tanzanian cholera outbreaks were caused by the 7PET lineage. The T5 sub-lineage (ctxB3) dominated outbreaks until 1997, followed by the T10 atypical El Tor (ctxB1) up to 2015, which were replaced by the T13 atypical El Tor of the current third wave (ctxB7) causing most cholera outbreaks until 2017 with T13 being phylogenetically related to strains from East African countries, Yemen and Lake Victoria. The strains were less drug resistant with approximate 10-kb deletions found in the SXT element, which encodes resistance to sulfamethoxazole and trimethoprim. Nucleotide deletions were observed in the CTX prophage of some strains, which warrants further virulence studies. Outbreak strains share 90% of core genes with V. cholerae O1 from Lake Victoria with as low as three SNPs difference and a significantly similar accessory genome, composed of genomic islands namely the CTX prophage, Vibrio Pathogenicity Islands; toxin co-regulated pilus biosynthesis proteins and the SXT-ICE element. Conclusion/Significance Characterization of V. cholerae O1 from Tanzania reveals genetic diversity of the 7PET lineage composed of T5, T10 and T13 sub-lineages with introductions of new sequence types from neighboring countries. The presence of these sub-lineages in environmental isolates suggests that the African Great Lakes may serve as aquatic reservoirs for survival of V. cholerae O1 favoring continuous human exposure., Author summary The seventh cholera pandemic has claimed >250,000 reported cases and 13,078 deaths until 2018 in Tanzania. To understand the epidemiology and to guide control, we used genomics to study V. cholerae O1 isolated in Tanzania during the past three decades. Tanzanian cholera outbreaks were caused by the T5, T10 and T13 sub-lineages of the 7PET lineage of V. cholerae O1 with some strains showing an unusual 100-bp deletion on the CTX prophage. From 1993 to 2017, most sub-lineages found in patients were also found in the aquatic environment and the close phylogenetic relationships between strains from the two niches suggest that the African Great Lakes may act as a reservoir for cholera outbreak strains. Moreover, we reported clonal transmission at regional and global scale favored by population displacements. Regional collaborative efforts are advised for effective cholera control.

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.

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.

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.

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.

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.

Introduction: Salmonella enterica serovar Typhimurium ST313 is an invasive and phylogenetically distinct lineage present in sub-Saharan Africa. We report the presence of S. Typhimurium ST313 from patients in the Democratic Republic of Congo and Nigeria. Methodology: Eighteen S. Typhimurium ST313 isolates were characterized by antimicrobial susceptibility testing, pulsed-field gel electrophoresis (PFGE), and multilocus sequence typing (MLST). Additionally, six of the isolates were characterized by whole genome sequence typing (WGST). The presence of a putative virulence determinant was examined in 177 Salmonella isolates belonging to 57 different serovars. Results: All S. Typhimurium ST313 isolates harbored resistant genes encoded by blaTEM1b, catA1, strA/B, sul1, and dfrA1. Additionally, aac(6')1aa gene was detected. Phylogenetic analyses revealed close genetic relationships among Congolese and Nigerian isolates from both blood and stool. Comparative genomic analyses identified a putative virulence fragment (ST313-TD) unique to S. Typhimurium ST313 and S. Dublin. Conclusion: We showed in a limited number of isolates that S. Typhimurium ST313 is a prevalent sequence-type causing gastrointestinal diseases and septicemia in patients from Nigeria and DRC. We found three distinct phylogenetic clusters based on the origin of isolation suggesting some spatial evolution. Comparative genomics showed an interesting putative virulence fragment (ST313-TD) unique to S. Typhimurium ST313 and invasive S. Dublin. © 2013 Leekitcharoenphon et al., SCOPUS: ar.j, info:eu-repo/semantics/published