Search

Pathogenic Vibrio parahaemolyticus ( Vp) ( tdh+/ trh+) represent a small percentage of environmental Vp populations, and very little is known about this subpopulation. Repetitive extragenic palindromic PCR and multilocus sequence analysis revealed heterogeneity among 41 Vp containing thermostable direct hemolysin ( tdh) and tdh-related hemolysin ( trh) that were isolated from Mississippi coastal environments from October 2006 to April 2007. There was no source-specific sequestering in oysters, water, or sediment. [ABSTRACT FROM AUTHOR]

Abstract Background Insertion elements (IS) are known to play an important role in the evolution and genomic diversification of Escherichia coli O157:H7 lineages. In particular, IS629 has been found in multiple copies in the E. coli O157:H7 genome and is one of the most prevalent IS in this serotype. It was recently shown that the lack of O157 antigen expression in two O rough E. coli O157:H7 strains was due to IS629 insertions at 2 different locations in the gne gene that is essential for the O antigen biosynthesis. Results The comparison of 4 E. coli O157:H7 genome and plasmid sequences showed numerous IS629 insertion sites, although not uniformly distributed among strains. Comparison of IS629s found in O157:H7 and O55:H7 showed the presence of at least three different IS629 sub-types. O157:H7 strains carry IS629 elements sub-type I and III whereby the ancestral O55:H7 carries sub-type II. Analysis of strains selected from various clonal groups defined on the E. coli O157:H7 stepwise evolution model showed that IS629 was not observed in sorbitol fermenting O157 (SFO157) clones that are on a divergent pathway in the emergence of O157:H7. This suggests that the absence of IS629 in SFO157 strains probably occurred during the divergence of this lineage, albeit it remains uncertain if it contributed, in part, to their divergence from other closely related strains. Conclusions The highly variable genomic locations of IS629 in O157:H7 strains of the A6 clonal complex indicates that this insertion element probably played an important role in genome plasticity and in the divergence of O157:H7 lineages.

Abstract Background The isolation of DNA-free RNA is a crucial step in the reverse transcription PCR (RT-PCR). Every RNA extraction procedure results in RNA samples contaminated with genomic DNA, which can cause false-positive outcomes in highly sensitive applications, including a recently developed quantitative real-time PCR (RT-qPCR) assay that targets invA mRNA for the detection of live Salmonella cells. The assay of this specific mRNA can be used to indicate the presence of live, as opposed to dead, cells of Salmonella enterica in a food matrix. Findings We evaluated the ability of five RNA extraction kits to produce RNA preparations from exponentially growing Salmonella cells. The acceptability of the preparations for use in downstream applications such as RT-qPCR was judged in terms of the total amount of RNA recovered, the integrity of the RNA molecules, and minimal content of DNA. The five kits produced RNA preparations that differed markedly in yield, integrity of the Salmonella RNA and the amount of contaminant DNA. The greatest RNA recovery was achieved with the MasterPure kit; however, the preparation contained high levels of genomic DNA. The UltraClean extraction kit gave a low level of RNA recovery with a poor level of integrity. The RNeasy Mini, RiboPure and PureLink extraction kits produced high-quality, DNA-free RNA suitable for Salmonella detection by RT-qPCR. Conclusions We showed that the RNeasy Mini and PureLink RNA extraction kits were the most suitable for the detection of Salmonella invA mRNA by RT-qPCR. The use of these two kits will greatly reduce the frequency of false-positive results and might allow fast RT-qPCR determination of invA mRNA produced by viable Salmonella in food samples.

We present the closed genome sequence of the Clostridium botulinum BT-22100019 strain isolated from the stool specimen of an infant diagnosed with botulism. With 4.33-Mb genome size and 28.0% G + C content, the bont/B1 gene encoded for botulinum neurotoxin serotype B was found on a 262 kb plasmid arranged in a ha + orfx - cluster., Competing Interests: The authors declare no conflict of interest.

Antibiotic resistance in Citrobacter freundii is a public health concern. This study evaluated the closed genome of a C. freundii isolated from the stool of a hospitalized patient initially related to a Salmonella outbreak. Confirmation of the isolate was determined by whole-genome sequencing. Nanopore sequencing was performed using a MinION with a Flongle flow cell. Assembly using SPAdes and Unicycler yielded a closed genome annotated by National Center for Biotechnology Information Prokaryotic Genome Annotation Pipeline. Genomic analyses employed MLST 2.0, ResFinder4.1, PlasmidFinder2.1, and VFanalyzer. Phylogenetic comparison utilized the Center for Food Safety and Applied Nutrition (CFSAN)-single nucleotide polymorphism pipeline and Genetic Algorithm for Rapid Likelihood Inference. Antimicrobial susceptibility was tested by broth microdilution following Clinical and Laboratory Standards Institute criteria. Multi-locus sequence type in silico analysis assigned the C. freundii as sequence type 64 and the blaCMY-41 gene was detected in resistome investigation. The susceptibility to antibiotics, determined using Sensititre® plates, revealed resistance to aztreonam, colistin, cefoxitin, amoxicillin/clavulanic acid, sulfisoxazole, ampicillin, and streptomycin. The genetic relatedness of the C. freundii CFSAN077772 with publicly available C. freundii genomes revealed a close relationship to a C. freundii SRR1186659, isolated in 2009 from human stool in Tanzania. In addition, C. freundii CFSAN077772 is nested in the same cluster with C. freundii clinical strains isolated in Denmark, Mexico, Myanmar, and Canada, suggesting a successful intercontinental spread., (© The Author(s) 2024. Published by Oxford University Press on behalf of Applied Microbiology International.)

Extended-spectrum cephalosporinase producing (ESC) E. coli are regarded as key indicator microorganisms of antimicrobial resistance (AMR), calling for a One Health integrated global surveillance strategy. Wildlife is exposed to antibiotic contaminants and/or resistant bacteria that have been released into the environment, potentially acting as reservoirs and spreaders of resistance genes as well as sentinels of anthropogenic pressure. Monitoring AMR in wildlife has become crucial in determining anthropogenic environmental impacts as well as transmission routes. In this study, we determined the occurrence and potential sources of ESC E. coli in 218 Eurasian griffon vultures (Gyps fulvus) foraging regularly on human waste disposed at a dumpsite in north-eastern Spain. Minimal inhibitory concentration for 14 different antimicrobials was performed to evaluate the phenotype of the isolates, and whole genome sequencing was carried out to investigate lineages and plasmids harbouring ESC genes. Our sequences were compared to previously published Spanish sequences of human, animal, and wildlife origin. We report a high prevalence of CTX-M-15, as well as the presence of other resistance genes such as OXA-10, CTX-M-27, and CTX-M-65 which are rarely described in European livestock, suggesting a human origin. The isolates also carried a diverse range of additional AMR genes for a broad spectrum of drug families, with the majority being multi-drug resistant. The phylogenomic analyses suggests the transmission of high-risk lineages from humans to vultures, with 49 % of our isolates matching the most common extraintestinal pathogenic E. coli (ExPEC) lineages described in humans worldwide, including ST131, ST10 and ST58. We conclude that anthropogenically altered habitats, such as landfills, are hotspots for the acquisition and spread of high-risk ESC E. coli lineages associated with hospital infections. Measures must be implemented to limit their spread into natural environments., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier B.V.)

Here, we report the complete genome of four S. enterica Infantis isolated in Costa Rica from human, poultry rinse, and raw chicken meat from 2017 to 2019. All genomes belonged to ST32 and carried a 310-kb plasmid with many antimicrobial resistance genes including the bla CTX-M65 gene., Competing Interests: The authors declare no conflict of interest.

The sharing of genome sequences in online data repositories allows for large scale analyses of specific genes or gene families. This can result in the detection of novel gene subtypes as well as the development of improved detection methods. Here, we used publicly available WGS data to detect a novel Stx subtype, Stx2n in two clinical E. coli strains isolated in the USA. During this process, additional Stx2 subtypes were detected; six Stx2j, one Stx2m strain, and one Stx2o, were all analyzed for variability from the originally described subtypes. Complete genome sequences were assembled from short- or long-read sequencing and analyzed for serotype, and ST types. The WGS data from Stx2n- and Stx2o-producing STEC strains were further analyzed for virulence genes pro-phage analysis and phage insertion sites. Nucleotide and amino acid maximum parsimony trees showed expected clustering of the previously described subtypes and a clear separation of the novel Stx2n subtype. WGS data were used to design OMNI PCR primers for the detection of all known stx 1 (283 bp amplicon), stx 2 (400 bp amplicon), intimin encoded by eae (221 bp amplicon), and stx 2f (438 bp amplicon) subtypes. These primers were tested in three different laboratories, using standard reference strains. An analysis of the complete genome sequence showed variability in serogroup, virulence genes, and ST type, and Stx2 pro-phages showed variability in size, gene composition, and phage insertion sites. The strains with Stx2j, Stx2m, Stx2n, and Stx2o showed toxicity to Vero cells. Stx2j carrying strain, 2012C-4221, was induced when grown with sub-inhibitory concentrations of ciprofloxacin, and toxicity was detected. Taken together, these data highlight the need to reinforce genomic surveillance to identify the emergence of potential new Stx2 or Stx1 variants. The importance of this surveillance has a paramount impact on public health. Per our description in this study, we suggest that 2017C-4317 be designated as the Stx2n type-strain.

Background: Listeria monocytogenes can survive in cold and wet environments, such as tree fruit packing facilities and it has been implicated in outbreaks and recalls of tree fruit products. However, little is known about microbiota that co-occurs with L. monocytogenes and its stability over seasons in tree fruit packing environments. In this 2-year longitudinal study, we aimed to characterize spatial and seasonal changes in microbiota composition and identify taxa indicative of L. monocytogenes contamination in wet processing areas of three tree fruit packing facilities (F1, F2, F3)., Methods: A total of 189 samples were collected during two apple packing seasons from floors under the washing, drying, and waxing areas. The presence of L. monocytogenes was determined using a standard culturing method, and environmental microbiota was characterized using amplicon sequencing. PERMANOVA was used to compare microbiota composition among facilities over two seasons, and abundance-occupancy analysis was used to identify shared and temporal core microbiota. Differential abundance analysis and random forest were applied to detect taxa indicative of L. monocytogenes contamination. Lastly, three L. monocytogenes-positive samples were sequenced using shotgun metagenomics with Nanopore MinION, as a proof-of-concept for direct detection of L. monocytogenes' DNA in environmental samples., Results: The occurrence of L. monocytogenes significantly increased from 28% in year 1 to 46% in year 2 in F1, and from 41% in year 1 to 92% in year 2 in F3, while all samples collected from F2 were L. monocytogenes-positive in both years. Samples collected from three facilities had a significantly different microbiota composition in both years, but the composition of each facility changed over years. A subset of bacterial taxa including Pseudomonas, Stenotrophomonas, and Microbacterium, and fungal taxa, including Yarrowia, Kurtzmaniella, Cystobasidium, Paraphoma, and Cutaneotrichosporon, were identified as potential indicators of L. monocytogenes within the monitored environments. Lastly, the DNA of L. monocytogenes was detected through direct Nanopore sequencing of metagenomic DNA extracted from environmental samples., Conclusions: This study demonstrated that a cross-sectional sampling strategy may not accurately reflect the representative microbiota of food processing facilities. Our findings also suggest that specific microorganisms are indicative of L. monocytogenes, warranting further investigation of their role in the survival and persistence of L. monocytogenes. Video Abstract., (© 2023. The Author(s).)

This study was conducted to determine the phylogenies of Salmonella strains isolated from cross-sectional studies conducted at hatcheries, broiler farms, processing plants, and retail outlets (broiler production chain) in Trinidad and Tobago over 4 yr (2016-2019). Whole-genome sequencing (WGS) was used to characterize Salmonella isolates. Core genome phylogenies of 8 serovars of public health significance were analyzed for similarities in origin and relatedness. In addition, Salmonella strains isolated from human salmonellosis cases in Trinidad were analyzed for their relatedness to the isolates detected along the broiler production chain. The common source of these isolates of diverse serovars within farms, within processing plants, between processing plants and retail outlets, and among farm-processing plant-retail outlet continuum was well-supported (bootstrap value >70%) by the core genome phylogenies for the respective serovars. Also, genome analyses revealed clustering of Salmonella serovars of regional (intra-Caribbean) and international (extra-Caribbean) origin. Similarly, strains of S. Enteritidis and S. Infantis isolated from human clinical salmonellosis in 2019 from Trinidad and Tobago clustered with our processing plant isolates recovered in 2018. This study is the first phylogenetic analysis of Salmonella isolates using WGS from the broiler industry in the Caribbean region. The use of WGS confirmed the genetic relatedness and transmission of Salmonella serovars contaminating chickens in broiler processing, and retailing in the country, with zoonotic and food safety implications for humans., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Escherichia coli ST1485 strains belong to the clinically important phylogroup F and have disseminated worldwide in humans, animals, and the environment. Here, we elucidated the pathogenome of a global collection of E. coli ST1485 isolates from diverse sources retrieved from public databases and a high-quality sequenced complete genome of colistin-resistant E. coli strain CFSAN061771 isolated from raw milk cheese which designated as a reference strain. CFSAN061771 belongs to O83:H42-ST1485 pathotype and carries a conjugative ColV plasmid, pCFSAN061771_01, combining extraintestinal virulence genes (ompt, sitA, iroN, etsC, traT, cvaC, hylF, iss, tsh, mchf, iucC, iutA) with a multidrug resistance island (bla TEM-1 , aph(6)-Id, aph(3″)-Ib, sul2, dfrA14). Comparative genomic analysis revealed a high frequency of pCFSAN061771_01-like plasmids in E. coli ST1485. A notable evolutionary genetic event in E. coli ST1485 strains is the acquisition of a pCFSAN061771_02-like plasmid, which confers resistance to several antimicrobials, tellurium, and quaternary ammonium compounds. The identical virulence and antibiotic resistance profiles identified in some human and animal strains are worrisome. This is the first study to emphasize the significance of E. coli ST1485 as a global high-risk virulent and multidrug-resistant clone with zoonotic potential., (© 2022. The Author(s).)

Actinobacillus pleuropneumoniae (APP) is the causative agent of pleuropneumonia in pigs, one of the most relevant bacterial respiratory diseases in the swine industry. To date, 19 serotypes have been described based on capsular polysaccharide typing with significant virulence dissimilarities. In this study, 16 APP isolates from Spanish origin were selected to perform antimicrobial susceptibility tests and comparative genomic analysis using whole genome sequencing (WGS). To obtain a more comprehensive worldwide molecular epidemiologic analyses, all APP whole genome assemblies available at the National Center for Biotechnology Information (NCBI) at the time of the study were also included. An in-house in silico PCR approach enabled the correct serotyping of unserotyped or incorrectly serotyped isolates and allowed for the discrimination between serotypes 9 and 11. A pangenome analysis identified the presence or absence of gene clusters to be serotype specific, as well as virulence profile analyses targeting the apx operons. Antimicrobial resistance genes were correlated to the presence of specific plasmids. Altogether, this study provides new insights into the genetic variability within APP serotypes, correlates phenotypic tests with bioinformatic analyses and manifests the benefits of populated databases for a better assessment of diversity and variability of relatively unknown pathogens. Overall, genomic comparative analysis enhances the understanding of transmission and epidemiological patterns of this species and suggests vertical transmission of the pathogen, including the resistance genes, within the Spanish integrated systems. IMPORTANCE Pleuropneumonia is one of the most relevant respiratory infections in the swine industry. Despite Actinobacillus pleuropneumoniae (APP) being one of the most important pathogens in the pig production, this is the first comparative study including all available whole genome sequencing data from NCBI. Moreover, this study also includes 16 APP isolates of Spanish origin with known epidemiological relationships through vertical integrated systems. Genomic comparisons provided a deeper understanding of molecular and epidemiological knowledge between different APP serotypes. Furthermore, determination of resistance and toxin profiles allowed correlation with the presence of mobile genetic elements and specific serotype, respectively.

Salmonella enterica is an important foodborne pathogen worldwide. We used long and short-read sequencing to close genomes of eight multidrug-resistant (MDR) S. enterica strains, belonging to serovars Infantis (2), Albany, Oranienburg, I 4,[5],12:i:-, Javiana, Schwarzengrund, and Kentucky from broiler chicken farms and processing plants in Trinidad and Tobago. They also belonged to seven different sequence types (STs- 32, 292, 1510, 19, 24, 152, and 96). Among the strains, seven had demonstrated multi-drug resistance with the presence of at least three AMR genes, whereas three isolates contained the quinolone resistance gene qnr B19 in plasmids (CFSAN103840, CFSAN103854, and CFSAN103872). The extended-spectrum β-lactamase genes bla CTX-M-65 (CFSAN103796) and bla TEM-1 (CFSAN103852) were detected in this study. The genomes closed in this study will be useful for future source tracking and outbreak investigations in Trinidad and Tobago and worldwide., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Maguire, Khan, Adesiyun, Georges and Gonzalez-Escalona.)

This cross-sectional study determined the serovars, antimicrobial resistance genes, and virulence factors of Salmonella isolated from hatcheries, broiler farms, processing plants, and retail outlets in Trinidad and Tobago. Salmonella in silico serotyping detected 23 different serovars where Kentucky 20.5% (30/146), Javiana 19.2% (28/146), Infantis 13.7% (20/146), and Albany 8.9% (13/146) were the predominant serovars. There was a 76.0% (111/146) agreement between serotyping results using traditional conventional methods and whole-genome sequencing (WGS) in in silico analysis. In silico identification of antimicrobial resistance genes conferring resistance to aminoglycosides, cephalosporins, peptides, sulfonamides, and antiseptics were detected. Multidrug resistance (MDR) was detected in 6.8% (10/146) of the isolates of which 100% originated from broiler farms. Overall, virulence factors associated with secretion systems and fimbrial adherence determinants accounted for 69.3% (3091/4463), and 29.2% (1302/4463) counts, respectively. Ten of 20 isolates of serovar Infantis (50.0%) showed MDR and contained the bla CTX-M-65 gene. This is the first molecular characterization of Salmonella isolates detected along the entire broiler production continuum in the Caribbean region using WGS. The availability of these genomes will help future source tracking during epidemiological investigations associated with Salmonella foodborne outbreaks in the region and worldwide.

Campylobacter is the leading cause of human bacterial gastroenteritis worldwide and has a major impact on global public health. Whole Genome Sequencing (WGS) is a powerful tool applied in the study of foodborne pathogens. The objective of the present study was to apply WGS to determine the genetic diversity, virulence factors and determinants of antimicrobial resistance of the populations of C. jejuni and C. coli in Peru. A total of 129 Campylobacter strains (108 C. jejuni and 21 C. coli ) were sequenced using Illumina Miseq platform. In silico MLST analysis identified a high genetic diversity among those strains with 30 sequence types (STs), several of them within 11 clonal complexes (CC) for C. jejuni , while the strains of C. coli belonged to a single CC with 8 different STs. Phylogeny analysis showed that Peruvian C. jejuni strains were divided into 2 clades with 5 populations, while C. coli formed a single clade with 4 populations. Furthermore, in silico analyses showed the presence of several genes associated with adherence, colonization and invasion among both species: cadF (83.7%), jlpA (81.4%), racR (100%), dnaJ (83.7%), pebA (83.7%), pldA (82.1%), porA (84.5%), ceuE (82.9%), ciaB (78.3%), iamB (86.8%), and flaC (100%). The majority (82.9%) of the Campylobacter strains carried the cdtABC operon which code for cytolethal distending toxin (CDT). Half of them (50.4%) carried genes associated with the presence of T6SS, while the frequency of genes associated with T4SS were relatively low (11.6%). Genetic markers associated with resistance to quinolones, tetracycline ( tetO, tetW/N/W ), beta-lactamases ( bla oxa-61 ), macrolides (A2075G in 23S rRNA) were found in 94.5, 21.7, 66.7, 6.2, 69.8, and 18.6% of strains, respectively. The cmeABC multidrug efflux operon was present in 78.3% of strains. This study highlights the importance of using WGS in the surveillance of emerging pathogens associated with foodborne diseases, providing genomic information on genetic diversity, virulence mechanisms and determinants of antimicrobial resistance. The description of several Campylobacter genotypes having many virulence factors and resistance to quinolones and tetracyclines circulating in Peru provides important information which helps in the monitoring, control and prevention strategies of this emerging pathogen in our country., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Quino, Caro-Castro, Hurtado, Flores-León, Gonzalez-Escalona and Gavilan.)

Leveraging systems biology approaches, we illustrate how metabolically distinct species of Clostridia protect against or worsen Clostridioides difficile infection in mice by modulating the pathogen's colonization, growth, and virulence to impact host survival. Gnotobiotic mice colonized with the amino acid fermenter Paraclostridium bifermentans survive infection with reduced disease severity, while mice colonized with the butyrate-producer, Clostridium sardiniense, succumb more rapidly. Systematic in vivo analyses revealed how each commensal alters the gut-nutrient environment to modulate the pathogen's metabolism, gene regulatory networks, and toxin production. Oral administration of P. bifermentans rescues conventional, clindamycin-treated mice from lethal C. difficile infection in a manner similar to that of monocolonized animals, thereby supporting the therapeutic potential of this commensal species. Our findings lay the foundation for mechanistically informed therapies to counter C. difficile disease using systems biology approaches to define host-commensal-pathogen interactions in vivo., Competing Interests: Declaration of interests L.B. and G.K.G. are co-inventors on patents for C. difficile microbiota therapeutics. L.B., G.K.G., and A.L.S. are SAB members and hold stock in ParetoBio. G.K.G.is an SAB member and holds stock in Kaleido, Inc. A.L.S. is a co-owner of ExArca Pharmaceuticals. The remaining authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Antimicrobial resistance (AMR) is a significant public health threat. With the rise of affordable whole genome sequencing, in silico approaches to assessing AMR gene content can be used to detect known resistance mechanisms and potentially identify novel mechanisms. To enable accurate assessment of AMR gene content, as part of a multi-agency collaboration, NCBI developed a comprehensive AMR gene database, the Bacterial Antimicrobial Resistance Reference Gene Database and the AMR gene detection tool AMRFinder. Here, we describe the expansion of the Reference Gene Database, now called the Reference Gene Catalog, to include putative acid, biocide, metal, stress resistance genes, in addition to virulence genes and species-specific point mutations. Genes and point mutations are classified by broad functions, as well as more detailed functions. As we have expanded both the functional repertoire of identified genes and functionality, NCBI released a new version of AMRFinder, known as AMRFinderPlus. This new tool allows users the option to utilize only the core set of AMR elements, or include stress response and virulence genes, too. AMRFinderPlus can detect acquired genes and point mutations in both protein and nucleotide sequence. In addition, the evidence used to identify the gene has been expanded to include whether nucleotide or protein sequence was used, its location in the contig, and presence of an internal stop codon. These database improvements and functional expansions will enable increased precision in identifying AMR genes, linking AMR genotypes and phenotypes, and determining possible relationships between AMR, virulence, and stress response.

Salmonella enterica is a highly important foodborne pathogen worldwide. We report the complete genome sequence of a sequence type 14 Salmonella enterica serotype Senftenberg strain carrying the mcr-9 gene in a plasmid isolated from broken chicken eggshells in Trinidad and Tobago, obtained by using a combination of long- and short-read sequencing.

Background: Whole genome sequencing of cultured pathogens is the state of the art public health response for the bioinformatic source tracking of illness outbreaks. Quasimetagenomics can substantially reduce the amount of culturing needed before a high quality genome can be recovered. Highly accurate short read data is analyzed for single nucleotide polymorphisms and multi-locus sequence types to differentiate strains but cannot span many genomic repeats, resulting in highly fragmented assemblies. Long reads can span repeats, resulting in much more contiguous assemblies, but have lower accuracy than short reads., Results: We evaluated the accuracy of Listeria monocytogenes assemblies from enrichments (quasimetagenomes) of naturally-contaminated ice cream using long read (Oxford Nanopore) and short read (Illumina) sequencing data. Accuracy of ten assembly approaches, over a range of sequencing depths, was evaluated by comparing sequence similarity of genes in assemblies to a complete reference genome. Long read assemblies reconstructed a circularized genome as well as a 71 kbp plasmid after 24 h of enrichment; however, high error rates prevented high fidelity gene assembly, even at 150X depth of coverage. Short read assemblies accurately reconstructed the core genes after 28 h of enrichment but produced highly fragmented genomes. Hybrid approaches demonstrated promising results but had biases based upon the initial assembly strategy. Short read assemblies scaffolded with long reads accurately assembled the core genes after just 24 h of enrichment, but were highly fragmented. Long read assemblies polished with short reads reconstructed a circularized genome and plasmid and assembled all the genes after 24 h enrichment but with less fidelity for the core genes than the short read assemblies., Conclusion: The integration of long and short read sequencing of quasimetagenomes expedited the reconstruction of a high quality pathogen genome compared to either platform alone. A new and more complete level of information about genome structure, gene order and mobile elements can be added to the public health response by incorporating long read analyses with the standard short read WGS outbreak response.

Whole genome sequencing (WGS) provides essential public health information and is used worldwide for pathogen surveillance, epidemiology, and source tracking. Foodborne pathogens are often sequenced using rapid library preparation chemistries based on transposon technology; however, this method may miss random segments of genomes that can be important for accurate downstream analyses. As new technologies become available, it may become possible to achieve better overall coverage. Here we compare the sequence quality obtained using libraries prepared from the Nextera XT and Nextera DNA Prep (Illumina, San Diego, CA) chemistries for 31 Shiga toxin-producing Escherichia coli (STEC) O121:H19 strains, which had been isolated from flour during a 2016 outbreak. The Nextera DNA Prep gave superior performance metrics including sequence quality, assembly quality, uniformity of genome coverage, and virulence gene identification, among other metrics. Comprehensive detection of virulence genes is essential for making educated assessments of STECs virulence potential. The phylogenetic SNP analysis did not show any differences in the variants detected by either library preparation method which allows isolates prepared from either library method to be analysed together. Our comprehensive comparison of these chemistries should assist researchers wishing to improve their sequencing workflow for STECs and other genomic risk assessments., Competing Interests: The authors have declared that no competing interests exist.

Campylobacter jejuni and Campylobacter coli are the leading cause of human gastroenteritis in the industrialized world and an emerging threat in developing countries. The incidence of campylobacteriosis in South America is greatly underestimated, mostly due to the lack of adequate diagnostic methods. Accordingly, there is limited genomic and epidemiological data from this region. In the present study, we performed a genome-wide analysis of the genetic diversity, virulence, and antimicrobial resistance of the largest collection of clinical C. jejuni and C. coli strains from Chile available to date (n = 81), collected in 2017-2019 in Santiago, Chile. This culture collection accounts for more than one third of the available genome sequences from South American clinical strains. cgMLST analysis identified high genetic diversity as well as 13 novel STs and alleles in both C. jejuni and C. coli. Pangenome and virulome analyses showed a differential distribution of virulence factors, including both plasmid and chromosomally encoded T6SSs and T4SSs. Resistome analysis predicted widespread resistance to fluoroquinolones, but low rates of erythromycin resistance. This study provides valuable genomic and epidemiological data and highlights the need for further genomic epidemiology studies in Chile and other South American countries to better understand molecular epidemiology and antimicrobial resistance of this emerging intestinal pathogen., Competing Interests: The authors have declared that no competing interests exist.

Whole genome sequencing (WGS) has made impressive progress in the field of molecular biology. Its most common application for public health is in the area of surveillance of food-borne diseases. WGS has the potential for providing a large amount of information, such as the identification of the strain type, the characterization of antibiotic resistance and virulence, and phylogeny. In our study, thirty-nine non-typhoidal Salmonella strains were isolated from diverse sources in Tunisia. Non-typhoidal Salmonella are among the most common pathogens contaminating food animals. The presence of virulence and antimicrobial resistance determinants in those strains were investigated using whole genome sequencing (WGS) and appropriate data analysis. The genomes were screened for several Salmonella virulence genes using the Virulence Factor Database VFDB. Twelve different virulence profiles, which correspond to the 12 identified serovars, were recognized. Several antimicrobial resistance genes were also detected: aac (6')-Iaa, sul1, tetA, bla-TEM and qnrS genes. Phylogenetic relationships among the strains were further assessed by a cgMLST analysis. The resulting phylogenetic tree consisted of several clusters consistently with the in silico multilocus sequence typing (MLST) and serotyping. Our findings demonstrated that WGS and subsequent data analysis provided an accurate tool for genetic characterization of bacterial strains compared to usual molecular typing techniques. To the best of our knowledge, this is the first report of an application of WGS for genetic characterization of food-borne Tunisian strains., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Foodborne pathogens have been implicated in illnesses worldwide. Here, we report the complete closed genome sequences of 28 bacterial strains belonging to 18 different species. These genomes belong to known foodborne pathogens. The genomes were closed by a combination of long-read and short-read sequencing.

Campylobacter species are the leading cause of gastroenteritis worldwide and an emerging threat in developing countries. Here, we report the draft whole-genome sequences of 51 Campylobacter jejuni and 12 Campylobacter coli strains isolated from patients with gastroenteritis in Santiago, Chile.

Clostridium botulinum is a Gram-positive, spore-forming anaerobic bacterium that produces botulinum neurotoxin (BoNT). Closing their genomes provides information about their neurotoxin clusters' arrangement(s) and their location (e.g., chromosome or plasmid) which cannot be assessed using draft genomes. Therefore, we tested the use of long-read sequencing (nanopore sequencing) in combination with short-read sequencing to close two toxin-producing strains. These genomes could be used by the Public Health Emergency Preparedness and Response staff during botulism outbreaks. The genomes of two toxin-producing C. botulinum strains, one from an environmental sample (83F_CFSAN034202) and the other from a clinical sample (CDC51232_CFSAN034200) were sequenced using MinION and MiSeq devices. The genomes, including the chromosomes and the plasmids, were closed by a combination of long-read and short-read sequencing. They belonged to different C. botulinum sequence types (STs), with 83F belonging to ST4 and CDC51232 to ST7. A whole genome single nucleotide polymorphism (SNP) analysis clustered these two strains with strains in lineage 2 (e.g., 6CDC297) and 4 (e.g., NCTC2916) from Group I, respectively. These two strains were also bivalent strains with the BoNTB and BoNTA4 clusters located in the larger plasmid for CDC51232, and the BoNTB and BoNTA1 clusters located both in the chromosome for 83F. Overall, this study showed the advantage of combining these two sequencing methods to obtain high quality closed C. botulinum genomes that could be used for SNP phylogenies (source tracking) as well as for fast identification of BoNT clusters and their gene arrangements., (Copyright © 2020 Gonzalez-Escalona and Sharma.)

We report transcontinental expansion of Vibrio parahaemolyticus sequence type 36 into Lima, Peru. From national collections, we identified 7 isolates from 2 different Pacific Northwest complex lineages that surfaced during 2011-2016. Sequence type 36 is likely established in environmental reservoirs. Systematic surveillance enabled detection of these epidemic isolates.

The complete genome sequences of Brevundimonas naejangsanensis strain FS1091 and Bacillus amyloliquefaciens strain FS1092, which were isolated from a commercial fresh-cut-produce-processing facility, were determined. Both FS1091 and FS1092 have one circular chromosome of approximately 3.15 and 4.24 Mb, respectively.

Colistin has become the last-line antimicrobial for the treatment of multidrug resistant (MDR) Enterobacterales in human medicine. To date, several colistin resistance genes have been described. Of them mcr -1 is disseminated worldwide in Escherichia coli of human and animal origin. The aim of this study was to characterize mcr -mediated resistance plasmids from E. coli of animal origin in Spain. From our strain collection, 70 E. coli of pig origin collected between 2005 and 2014 (10 per year, except for years 2009-2010-2013) were randomly selected and screened for the presence of mcr -genes. Additionally, 20 E. coli isolated in 2011 from white storks ( Ciconia ciconia ) from the same urban household waste landfill associated colony were also included. Whole genome sequencing of mcr -positive isolates was carried out on a MiSeq (Illumina). Hybrid whole genome sequencing strategy combining nanopore and Illumina technologies were performed in a selection of isolates to close the genomes and plasmids and identify the presence of antimicrobial resistance genes. Minimum inhibitory concentration (MIC) was used to assess the susceptibility to colistin. Mating experiments were carried out to evaluate transferability of the mcr -genes. A total of 19 mcr -1 and one mcr -4 positive isolates were detected, 15 from pigs distributed during the study period, and five from storks collected in 2011. No other mcr -variants were found. The MICs for colistin ranged between 4 and >4 mg/L. High diversity of STs were detected among the mcr-1 positive E. coli isolates, with only ST-10 shared between pigs and white storks. Except for one isolate, all were genotypic and phenotypically MDR, and five of them also harbored cephalosporin resistance genes ( bla CTX-M- 14 , bla SHV- 12 , and three bla CMY- 2 ). mcr -1 genes were mobilizable by conjugation, associated with IncX4, IncHI2, and IncI2 plasmids. In our study, mcr -1 genes have been circulating in pig farms since 2005 harbored by a variety of E. coli clones. Its persistence may be driven by co-selection since plasmids containing mcr -1 also exhibit resistance to multiple drugs used in veterinary medicine. Furthermore, this is the first report of the presence of mcr -1 gene in isolates from white storks in Spain. This finding highlights the potential importance of wildlife that forage at urban household waste landfills in the transmission and spread of colistin resistance genes., (Copyright © 2020 Migura-Garcia, González-López, Martinez-Urtaza, Aguirre Sánchez, Moreno-Mingorance, Perez de Rozas, Höfle, Ramiro and Gonzalez-Escalona.)

Here, we report the genomes of two Salmonella enterica subsp. enterica serovar Montevideo strains (CFSAN005645 and FCC0123) and two Salmonella enterica subsp. enterica serovar Senftenberg strains (FSW0104 and CFSAN087304) isolated from pistachios. The genomes were closed using a hybrid assembly method using short- and long-read sequencing technology., (Copyright © 2019 Haendiges et al.)

There is emerging evidence that food of animal origin may be responsible for the spread of multidrug resistant extraintestinal pathogenic Escherichia coli in the community. Here, we describe the emergence of colistin resistance gene, mcr-1, in a strain belonging to the dominant uropathogenic E. coli ST69 lineage. E. coli strain CFSAN061770 was isolated during monitoring of the popular Egyptian raw milk cheese, karish cheese, for the presence of colistin resistance. The complete genome of E. coli strain CFSAN061770 comprises a chromosome of 5,292,297 bp with a G + C content of 50.6%. Further, three plasmids named pEGY1-MCR-1, pEGY2 and pEGY3 of 228,947 bp, 103,234 bp and 87,012 bp were detected, respectively. Plasmid pEGY1-MCR-1 belongs to the IncHI2 incompatibility group and carries the colistin resistance mcr-1 gene flanked by two ISApl1 elements and forms a composite transposon. It mediates resistance to aminoglycosides (aadA1 and aadA2), phenicol (cmlA1 and floR), sulfonamides (sul3), and tetracycline (tet(A)), and these loci were found clustered in a multidrug resistant region. Plasmid pEGY3 carries a complex multiple resistance locus (CMR) (aph(3')-Ia, strA, strB, sul2, and bla TEM-1 ) encoding resistance to different classes of antibiotics. Interestingly, the closest plasmids to plasmid pEGY1-MCR-1 detected from the NCBI Blast search belonged to the incompatibility group IncHI2 and were from the Kingdom of Saudi Arabia and Qatar which suggests a dissemination of pEGY1-MCR-1-like plasmids in the Middle East. Most striking, and of great public health concern is that strain CFSAN061770 carries five virulence genes (iss, fimH, iutA, kpsMIII and kpsTIII) which were identified in clinical extraintestinal pathogenic E. coli. Besides that, it carries the astA gene, which codes for the enteroaggregative E. coli heat-stable toxin 1 (EAST1)., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Here, we report the genome sequences of three Salmonella enterica strains belonging to serovars Weltevreden (CFSAN047349), Saintpaul (CFSAN047351), and Thompson (CFSAN047352), isolated from river water in Sinaloa, Mexico. The genomes were closed by a combination of long-read and short-read sequencing. The strain sequence types (STs) are ST365, ST50, and ST26, respectively.

Foodborne pathogen surveillance in the United States is transitioning from strain identification using restriction digest technology (pulsed-field gel electrophoresis [PFGE]) to shotgun sequencing of the entire genome (whole-genome sequencing [WGS]). WGS requires a new suite of analysis tools, some of which have long histories in academia but are new to the field of public health and regulatory decision making. Although the general workflow is fairly standard for collecting and analyzing WGS data for disease surveillance, there are a number of differences in how the data are collected and analyzed across public health agencies, both nationally and internationally. This impedes collaborative public health efforts, so national and international efforts are underway to enable direct comparison of these different analysis methods. Ultimately, the harmonization efforts will allow the (mutually trusted and understood) production and analysis of WGS data by labs and agencies worldwide, thus improving outbreak response capabilities globally. This review provides a historical perspective on the use of WGS for pathogen tracking and summarizes the efforts underway to ensure the major steps in phylogenomic pipelines used for pathogen disease surveillance can be readily validated. The tools for doing this will ensure that the results produced are sound, reproducible, and comparable across different analytic approaches.

Nineteen Salmonella strains were isolated from 5,907 randomly selected mussel samples during a monitoring programme for the presence of Salmonella in shellfish in Galicia, north-west Spain (2012-16). Serovars, sequence type and antimicrobial resistance genes were determined through genome sequencing. Presence of the mcr -1 gene in one strain belonging to serovar Rissen and ST-469 was identified. The mcr -1 gene had not been isolated previously in environmental Salmonella isolated from mussels in Spain.

Here we report the genome sequences of two toxin-producing Clostridium botulinum strains, one environmental sample (83F) and one clinical sample (CDC51232). The genomes were closed by a combination of long-read and short-read sequencing. The strains belong to C. botulinum sequence type 4 (ST4) and ST7, respectively.

Galicia in northwestern Spain has been considered a hotspot for Vibrio parahaemolyticus infections. Infections abruptly emerged in 1998 and, over the next 15 years, were associated with large outbreaks caused by strains belonging to a single clone. We report a recent transition in the epidemiologic pattern in which cases throughout the region have been linked to different and unrelated strains. Global genome-wide phylogenetic analysis revealed that most of the pathogenic strains isolated from infections were associated with globally diverse isolates, indicating frequent episodic introductions from disparate and remote sources. Moreover, we identified that the 2 major switches in the epidemic dynamics of V. parahaemolyticus in the regions, the emergence of cases and an epidemiologic shift in 2015-2016, were associated with the rise of sea surface temperature in coastal areas of Galicia. This association may represent a fundamental contributing factor in the emergence of illness linked to these introduced pathogenic strains.

Whole genome sequencing (WGS) has been broadly used to provide detailed characterization of foodborne pathogens. These genomes for diverse species including Salmonella, Escherichia coli, Listeria, Campylobacter and Vibrio have provided great insight into the genetic make-up of these pathogens. Numerous government agencies, industry and academia have developed new applications in food safety using WGS approaches such as outbreak detection and characterization, source tracking, determining the root cause of a contamination event, profiling of virulence and pathogenicity attributes, antimicrobial resistance monitoring, quality assurance for microbiology testing, as well as many others. The future looks bright for additional applications that come with the new technologies and tools in genomics and metagenomics., (Published by Elsevier Ltd.)

Vibrio parahaemolyticus is the leading cause of seafood-related infections with illnesses undergoing a geographic expansion. In this process of expansion, the most fundamental change has been the transition from infections caused by local strains to the surge of pandemic clonal types. Pandemic clone sequence type 3 (ST3) was the only example of transcontinental spreading until 2012, when ST36 was detected outside the region where it is endemic in the U.S. Pacific Northwest causing infections along the U.S. northeast coast and Spain. Here, we used genome-wide analyses to reconstruct the evolutionary history of the V. parahaemolyticus ST36 clone over the course of its geographic expansion during the previous 25 years. The origin of this lineage was estimated to be in ~1985. By 1995, a new variant emerged in the region and quickly replaced the old clone, which has not been detected since 2000. The new Pacific Northwest (PNW) lineage was responsible for the first cases associated with this clone outside the Pacific Northwest region. After several introductions into the northeast coast, the new PNW clone differentiated into a highly dynamic group that continues to cause illness on the northeast coast of the United States. Surprisingly, the strains detected in Europe in 2012 diverged from this ancestral group around 2000 and have conserved genetic features present only in the old PNW lineage. Recombination was identified as the major driver of diversification, with some preliminary observations suggesting a trend toward a more specialized lifestyle, which may represent a critical element in the expansion of epidemics under scenarios of coastal warming. IMPORTANCE Vibrio parahaemolyticus and Vibrio cholerae represent the only two instances of pandemic expansions of human pathogens originating in the marine environment. However, while the current pandemic of V. cholerae emerged more than 50 years ago, the global expansion of V. parahaemolyticus is a recent phenomenon. These modern expansions provide an exceptional opportunity to study the evolutionary process of these pathogens at first hand and gain an understanding of the mechanisms shaping the epidemic dynamics of these diseases, in particular, the emergence, dispersal, and successful introduction in new regions facilitating global spreading of infections. In this study, we used genomic analysis to examine the evolutionary divergence that has occurred over the course of the most recent transcontinental expansion of a pathogenic Vibrio , the spreading of the V. parahaemolyticus sequence type 36 clone from the region where it is endemic on the Pacific coast of North America to the east coast of the United States and finally to the west coast of Europe.

Shellfish-transmitted Vibrio parahaemolyticus infections have recently increased from locations with historically low disease incidence, such as the Northeast United States. This change coincided with a bacterial population shift toward human-pathogenic variants occurring in part through the introduction of several Pacific native lineages (ST36, ST43, and ST636) to nearshore areas off the Atlantic coast of the Northeast United States. Concomitantly, ST631 emerged as a major endemic pathogen. Phylogenetic trees of clinical and environmental isolates indicated that two clades diverged from a common ST631 ancestor, and in each of these clades, a human-pathogenic variant evolved independently through acquisition of distinct Vibrio pathogenicity islands (VPaI). These VPaI differ from each other and bear little resemblance to hemolysin-containing VPaI from isolates of the pandemic clonal complex. Clade I ST631 isolates either harbored no hemolysins or contained a chromosome I-inserted island we call VPaIβ that encodes a type 3 secretion system (T3SS2β) typical of Trh hemolysin producers. The more clinically prevalent and clonal ST631 clade II had an island we call VPaIγ that encodes both tdh and trh and that was inserted in chromosome II. VPaIγ was derived from VPaIβ but with some additional acquired elements in common with VPaI carried by pandemic isolates, exemplifying the mosaic nature of pathogenicity islands. Genomics comparisons and amplicon assays identified VPaIγ-type islands containing tdh inserted adjacent to the ure cluster in the three introduced Pacific and most other emergent lineages that collectively cause 67% of infections in the Northeast United States as of 2016. IMPORTANCE The availability of three different hemolysin genotypes in the ST631 lineage provided a unique opportunity to employ genome comparisons to further our understanding of the processes underlying pathogen evolution. The fact that two different pathogenic clades arose in parallel from the same potentially benign lineage by independent VPaI acquisition is surprising considering the historically low prevalence of community members harboring VPaI in waters along the Northeast U.S. coast that could serve as the source of this material. This illustrates a possible predisposition of some lineages to not only acquire foreign DNA but also become human pathogens. Whereas the underlying cause for the expansion of V. parahaemolyticus lineages harboring VPaIγ along the U.S. Atlantic coast and spread of this element to multiple lineages that underlies disease emergence is not known, this work underscores the need to define the environment factors that favor bacteria harboring VPaI in locations of emergent disease., (Copyright © 2017 American Society for Microbiology.)

Background: Enterohemorrhagic Escherichia coli (EHEC) O145 are among the top non-O157 serogroups associated with severe human disease worldwide. Two serotypes, O145:H25 and O145:H28 have been isolated from human patients but little information is available regarding the virulence repertoire, origin and evolutionary relatedness of O145:H25. Hence, we sequenced the complete genome of two O145:H25 strains associated with hemolytic uremic syndrome (HUS) and compared the genomes with those of previously sequenced O145:H28 and other EHEC strains., Results: The genomes of the two O145:H25 strains were 5.3 Mbp in size; slightly smaller than those of O145:H28 and other EHEC strains. Both strains contained three nearly identical plasmids and several prophages and integrative elements, many of which differed significantly in size, gene content and organization as compared to those present in O145:H28 and other EHECs. Furthermore, notable variations were observed in several fimbrial gene cluster and intimin types possessed by O145:H25 and O145:H28 indicating potential adaptation to distinct areas of host colonization. Comparative genomics further revealed that O145:H25 are genetically more similar to other non-O157 EHEC strains than to O145:H28., Conclusion: Phylogenetic analysis accompanied by comparative genomics revealed that O145:H25 and O145:H28 evolved from two separate clonal lineages and that horizontal gene transfer and gene loss played a major role in the divergence of these EHEC serotypes. The data provide further evidence that ruminants might be a possible reservoir for O145:H25 but that they might be impaired in their ability to establish a persistent colonization as compared to other EHEC strains.

Determinants of multidrug resistance (MDR) are often encoded on mobile elements, such as plasmids, transposons, and integrons, which have the potential to transfer among foodborne pathogens, as well as to other virulent pathogens, increasing the threats these traits pose to human and veterinary health. Our understanding of MDR among Salmonella has been limited by the lack of closed plasmid genomes for comparisons across resistance phenotypes, due to difficulties in effectively separating the DNA of these high-molecular weight, low-copy-number plasmids from chromosomal DNA. To resolve this problem, we demonstrate an efficient protocol for isolating, sequencing and closing IncA/C plasmids from Salmonella sp. using single molecule real-time sequencing on a Pacific Biosciences (Pacbio) RS II Sequencer. We obtained six Salmonella enterica isolates from poultry, representing six different serovars, each exhibiting the MDR-Ampc resistance profile. Salmonella plasmids were obtained using a modified mini preparation and transformed with Escherichia coli DH10Br. A Qiagen Large-Construct kit™ was used to recover highly concentrated and purified plasmid DNA that was sequenced using PacBio technology. These six closed IncA/C plasmids ranged in size from 104 to 191 kb and shared a stable, conserved backbone containing 98 core genes, with only six differences among those core genes. The plasmids encoded a number of antimicrobial resistance genes, including those for quaternary ammonium compounds and mercury. We then compared our six IncA/C plasmid sequences: first with 14 IncA/C plasmids derived from S. enterica available at the National Center for Biotechnology Information (NCBI), and then with an additional 38 IncA/C plasmids derived from different taxa. These comparisons allowed us to build an evolutionary picture of how antimicrobial resistance may be mediated by this common plasmid backbone. Our project provides detailed genetic information about resistance genes in plasmids, advances in plasmid sequencing, and phylogenetic analyses, and important insights about how MDR evolution occurs across diverse serotypes from different animal sources, particularly in agricultural settings where antimicrobial drug use practices vary.

Vibrio parahaemolyticus is an important human foodborne pathogen whose transmission is associated with the consumption of contaminated seafood, with a growing number of infections reported over recent years worldwide. A multilocus sequence typing (MLST) database for V. parahaemolyticus was created in 2008, and a large number of clones have been identified, causing severe outbreaks worldwide (sequence type 3 [ST3]), recurrent outbreaks in certain regions (e.g., ST36), or spreading to other regions where they are nonendemic (e.g., ST88 or ST189). The current MLST scheme uses sequences of 7 genes to generate an ST, which results in a powerful tool for inferring the population structure of this pathogen, although with limited resolution, especially compared to pulsed-field gel electrophoresis (PFGE). The application of whole-genome sequencing (WGS) has become routine for trace back investigations, with core genome MLST (cgMLST) analysis as one of the most straightforward ways to explore complex genomic data in an epidemiological context. Therefore, there is a need to generate a new, portable, standardized, and more advanced system that provides higher resolution and discriminatory power among V. parahaemolyticus strains using WGS data. We sequenced 92 V. parahaemolyticus genomes and used the genome of strain RIMD 2210633 as a reference (with a total of 4,832 genes) to determine which genes were suitable for establishing a V. parahaemolyticus cgMLST scheme. This analysis resulted in the identification of 2,254 suitable core genes for use in the cgMLST scheme. To evaluate the performance of this scheme, we performed a cgMLST analysis of 92 newly sequenced genomes, plus an additional 142 strains with genomes available at NCBI. cgMLST analysis was able to distinguish related and unrelated strains, including those with the same ST, clearly showing its enhanced resolution over conventional MLST analysis. It also distinguished outbreak-related from non-outbreak-related strains within the same ST. The sequences obtained from this work were deposited and are available in the public database (http://pubmlst.org/vparahaemolyticus). The application of this cgMLST scheme to the characterization of V. parahaemolyticus strains provided by different laboratories from around the world will reveal the global picture of the epidemiology, spread, and evolution of this pathogen and will become a powerful tool for outbreak investigations, allowing for the unambiguous comparison of strains with global coverage., (Copyright © 2017 Gonzalez-Escalona et al.)