Most efforts to understand the biology of Vibrio cholerae have focused on a single group, the pandemic-generating lineage harbouring the strains responsible for all known cholera pandemics. Consequently, little is known about the diversity of this species in its native aquatic environment. To understand the differences in the V. cholerae populations inhabiting in regions with varying history of cholera cases and how that might influence the abundance of pandemic strains, a comparative analysis of population composition was performed. Little overlap was found in lineage compositions between those in Dhaka (cholera endemic) located in the Ganges delta, and of Falmouth (no known history of cholera), a small coastal town on the US East Coast. The most striking difference was the presence of a group of related lineages at high abundance in Dhaka which was completely absent from Falmouth. Phylogenomic analysis revealed that these lineages form a cluster at the base of the phylogeny of V. cholerae species, sufficiently differentiated genetically and phenotypically to form a novel species. Strains from this species have been anecdotally isolated from around the world and were isolated as early as 1916 from a British soldier in Egypt suffering from choleraic diarrhoea. In 1935 Gardner and Venkatraman unofficially referred to a member of this group as Vibrio paracholerae. In recognition of this earlier designation, we propose the name Vibrio paracholerae, sp. nov. for this bacterium. Genomic analysis suggests a link with human populations for this novel species and substantial interaction with its better-known sister species.Importance Cholera continues to remain a major public health threat around the globe. Understanding the ecology, evolution and environmental adaptation of the causative agent Vibrio cholerae and tracking the emergence of novel lineages with pathogenic potential are essential to combat the problem. In this study, we investigated the population dynamics of Vibrio cholerae in an inland locality which is known as endemic for cholera and compared with that of a cholera free coastal location. We found the consistent presence of the pandemic generating V. cholerae in cholera-endemic Dhaka and an exclusive presence of a lineage phylogenetically distinct from other V. cholerae. Our study suggests that this lineage represents a novel species having pathogenic potential and a human link to its environmental abundance. The possible association with human population, co-existence and interaction with toxigenic V. cholerae in the natural environment make this potential human pathogen an important subject for future studies.Competing Interest StatementThe authors have declared no competing interest.

Retracing microbial emergence and spread is essential to understanding the evolution and dynamics of pathogens. The bacterial foodborne pathogen Listeria monocytogenes clonal complex 1 (Lm-CC1) is the most prevalent clonal group associated with listeriosis, and is strongly associated with cattle and dairy products. Here we analysed 2,021 Lm-CC1 isolates collected from 40 countries, since the first Lm isolation to the present day, to define its evolutionary history and population dynamics. Our results suggest that Lm-CC1 spread worldwide from North America following the Industrial Revolution through two waves of expansion, coinciding with the transatlantic livestock trade in the second half of the 19th century and the rapid growth of cattle farming in the 20th century. Lm-CC1 then firmly established at a local level, with limited inter-country spread. This study provides an unprecedented insight into Lm-CC1 phylogeography and dynamics and can contribute to effective disease surveillance to reduce the burden of listeriosis.Competing Interest StatementThe authors have declared no competing interest.

Core genome multilocus sequence typing (cgMLST) has gained popularity in recent years in epidemiological research and subspecies level classification. cgMLST retains the intuitive nature of traditional MLST but offers much greater resolution by utilizing significantly larger portions of the genome. Here, we introduce a cgMLST scheme for Vibrio cholerae, a bacterium abundant in marine and freshwater environments and the etiologic agent of cholera. A set of 2,443 core genes ubiquitous in V. cholerae were used to analyze a comprehensive dataset of 1,262 clinical and environmental strains collected from 52 countries, including 65 newly sequenced genomes in this study. We established a sublineage threshold based on 133 allelic differences that creates clusters nearly identical to traditional MLST types, providing backwards compatibility to new cgMLST classifications. We also defined an outbreak threshold based on seven allelic differences that is capable of identifying strains from the same outbreak and closely related isolates which could give clues on outbreak origin. Using cgMLST, we confirmed the South Asian origin of modern epidemics and identified clustering affinity among sublineages of environmental isolates from the same geographic origin. Advantages of this method are highlighted by direct comparison with existing classification methods, such as MLST and single nucleotide polymorphism-based methods. cgMLST outperforms all existing methods in terms of resolution, standardization, and ease-of-use. We anticipate this scheme will serve as a basis for a universally applicable and standardized classification system for V. cholerae research and epidemiological surveillance in the future. This cgMLST scheme is publicly available on PubMLST (https://pubmlst.org/vcholerae/).IMPORTANCE Toxigenic Vibrio cholerae of the O1 and O139 serogroups are the causative agent of cholera, an acute diarrheal disease that plagued the world for centuries, if not millennia. Here, we introduce a core genome multilocus sequence typing (cgMLST) scheme for V. cholerae. Using cgMLST, we established an outbreak threshold that can efficiently identify outbreak related strains and potential sources of introduction. We also defined a sublineage threshold that is similar to traditional MLST sequence type which will provide context to this new typing method by relating it to previous MLST results. cgMLST outperforms all existing methods in terms of resolution, standardization, and ease-of-use, making this scheme the most suitable method for V. cholerae typing and surveillance worldwide.

Toxigenic Vibrio cholerae serogroup O1 is the etiologic agent of the disease cholera, and strains of this serogroup are responsible for pandemics. A few other serogroups have been found to carry cholera toxin genes-most notably, O139, O75, and O141-and public health surveillance in the United States is focused on these four serogroups. A toxigenic isolate was recovered from a case of vibriosis from Texas in 2008. This isolate did not agglutinate with any of the four different serogroups' antisera (O1, O139, O75, or O141) routinely used in phenotypic testing and did not display a rough phenotype. We investigated several hypotheses that might explain the recovery of this potential nonagglutinating (NAG) strain using whole-genome sequencing analysis and phylogenetic methods. The NAG strain formed a monophyletic cluster with O141 strains in a whole-genome phylogeny. Furthermore, a phylogeny of ctxAB and tcpA sequences revealed that the sequences from the NAG strain also formed a monophyletic cluster with toxigenic U.S. Gulf Coast (USGC) strains (O1, O75, and O141) that were recovered from vibriosis cases associated with exposures to Gulf Coast waters. A comparison of the NAG whole-genome sequence showed that the O-antigen-determining region of the NAG strain was closely related to those of O141 strains, and specific mutations were likely responsible for the inability to agglutinate. This work shows the utility of whole-genome sequence analysis tools for characterization of an atypical clinical isolate of V. cholerae originating from a USGC state. IMPORTANCE Clinical cases of vibriosis are on the rise due to climate events and ocean warming (1, 2), and increased surveillance of toxigenic Vibrio cholerae strains is now more crucial than ever. While traditional phenotyping using antisera against O1 and O139 is useful for monitoring currently circulating strains with pandemic or epidemic potential, reagents are limited for non-O1/non-O139 strains. With the increased use of next-generation sequencing technologies, analysis of less well-characterized strains and O-antigen regions is possible. The framework for advanced molecular analysis of O-antigen-determining regions presented herein will be useful in the absence of reagents for serotyping. Furthermore, molecular analyses based on whole-genome sequence data and using phylogenetic methods will help characterize both historical and novel strains of clinical importance. Closely monitoring emerging mutations and trends will improve our understanding of the epidemic potential of Vibrio cholerae to anticipate and rapidly respond to future public health emergencies.

Retracing microbial emergence and spread is essential to understanding the evolution and 40 dynamics of pathogens. The bacterial foodborne pathogen Listeria monocytogenes clonal 41 complex 1 (Lm-CC1) is the most prevalent clonal group associated with listeriosis, and is 42 strongly associated with cattle and dairy products. Here we analysed 2,021 Lm-CC1 43 isolates collected from 40 countries, since the first Lm isolation to the present day, to 44 define its evolutionary history and population dynamics. Our results suggest that Lm-CC1 45 spread worldwide from North America following the Industrial Revolution through two 46 waves of expansion, coinciding with the transatlantic livestock trade in the second half of 47 the 19th century and the rapid growth of cattle farming in the 20th century. Lm-CC1 then 48 firmly established at a local level, with limited inter-country spread. This study provides 49 an unprecedented insight into Lm-CC1 phylogeography and dynamics and can contribute 50 to effective disease surveillance to reduce the burden of listeriosis.

Background: While rotavirus vaccine programs effectively protect against severe rotavirus gastroenteritis, rotavirus vaccine strains have been identified in the stool of vaccinated children and their close contacts suffering from acute gastroenteritis. The prevalence of vaccine strains, the emergence of vaccine-derived strains and their role in acute gastroenteritis are not well studied. Methods: We developed a Locked Nucleic Acid Reverse Transcription real-time PCR assay (LNA-RTqPCR) to detect the monovalent rotavirus vaccine (RV1) Rotarix non-structural protein 2 in children with acute gastroenteritis and healthy controls and validated it using sequence confirmed RV1 strains. The association between RV1-derived strains and gastroenteritis was determined using logistic regression. Results: The new assay exhibited 100% (95%CI: 91.7%, 100%) diagnostic sensitivity and 99.4% (95%CI: 96.2%, 100%) diagnostic specificity, with a detection limit of 9.86 copies/reaction and qPCR efficiency of 99.7%. Using this assay, we identified the presence of RV1-derived NSP2 sequences in 7.7% of rotavirus gastroenteritis cases and 98.6% of rotavirus positive healthy children (94.4% had previously received the RV1). Among gastroenteritis cases, those whose stool contained RV1-derived strains had milder gastroenteritis symptoms compared to that of natural rotavirus infections. We observed no significant association between RV1-derived strains and gastroenteritis (OR 0.98; 95% CI 0.60, 1.72). Conclusion: Our study demonstrated that the new assay is suitable for monitoring RV1-derived rotavirus strain circulation and that the RV1-derived strains are not associated with development of gastroenteritis symptoms.

Most efforts to understand the biology of Vibrio cholerae have focused on a single group, the pandemic-generating lineage harbouring the strains responsible for all known cholera pandemics. Consequently, little is known about the diversity of this species in its native aquatic environment. To understand the differences in the V. cholerae populations inhabiting regions with a history of cholera cases and those lacking such a history, a comparative analysis of population composition was performed. Little overlap was found in lineage compositions between those in Dhaka (cholera endemic) located in the Ganges delta, and of Falmouth (no known history of cholera), a small coastal town on the United States east coast. The most striking difference was the presence of a group of related lineages at high abundance in Dhaka which was completely absent from Falmouth. Phylogenomic analysis revealed that these lineages form a cluster at the base of the phylogeny for the V. cholerae species, sufficiently differentiated genetically and phenotypically to form a novel species. A retrospective search revealed that strains from this species have been anecdotally found from around the world and were isolated as early as 1916 from a British soldier in Egypt suffering from choleraic diarrhoea. In 1935 Gardner and Venkatraman unofficially referred to a member of this group as Vibrio paracholerae. In recognition of this earlier designation, we propose the name Vibrio paracholerae sp. nov. for this bacterium. Genomic analysis suggests a link with human populations for this novel species and substantial interaction with its better-known sister species. Importance Cholera continues to remain a major public health threat around the globe. Understanding the ecology, evolution, and environmental adaptation of the causative agent Vibrio cholerae and tracking the emergence of novel lineages with pathogenic potential are essential to combat the problem. In this study, we investigated the population dynamics of Vibrio cholerae in an inland locality which is known as endemic for cholera and compared with that of a cholera free coastal location. We found the consistent presence of the pandemic generating V. cholerae in cholera-endemic Dhaka and an exclusive presence of a lineage phylogenetically distinct from other V. cholerae. Our study suggests that this lineage represents a novel species having pathogenic potential and a human link to its environmental abundance. The possible association with human population, co-existence and interaction with toxigenic V. cholerae in the natural environment make this potential human pathogen an important subject for future studies.

In July 2019, we investigated a cluster of Yersinia enterocolitica cases affecting a youth summer camp and nearby community in northeastern Pennsylvania. After initial telephone interviews with camp owners and community members, we identified pasteurized milk from a small dairy conducting on-site pasteurization, Dairy A, as a shared exposure. We conducted site visits at the camp and Dairy A where we collected milk and other samples. Samples were cultured for Y. enterocolitica. Clinical and nonclinical isolates were compared using molecular subtyping. We performed case finding, conducted telephone interviews for community cases, and conducted a cohort study among adult camp staff by administering an online questionnaire. In total, we identified 109 Y. enterocolitica cases. Consumption of Dairy A milk was known for 37 (34%); of these, Dairy A milk was consumed by 31 (84%). Dairy A had shipped 214 gallons of pasteurized milk in 5 weekly shipments to the camp by mid-July. Dairy A milk was the only shared exposure identified between the camp and community. Y. enterocolitica was isolated from Dairy A unpasteurized milk samples. Five clinical isolates from camp members, two clinical isolates from community members, and nine isolates from unpasteurized milk were indistinguishable by whole-genome sequencing. The risk for yersinosis among camp staff who drank Dairy A milk was 5.3 times the risk for those who did not (95% confidence interval: 1.6-17.3). Because Dairy A only sold pasteurized milk, pasteurized milk was considered the outbreak source. We recommend governmental agencies and small dairies conducting on-site pasteurization collaborate to develop outbreak prevention strategies.

Core genome multilocus sequence typing (cgMLST) has gained popularity in recent years in epidemiological research and subspecies level classification. cgMLST retains the intuitive nature of traditional MLST but offers much greater resolution by utilizing significantly larger portions of the genome. Here, we introduce a cgMLST scheme for Vibrio cholerae, a bacterium abundant in marine and freshwater environments and the etiologic agent of cholera. A set of 2,443 core genes ubiquitous in V. cholerae were used to analyze a comprehensive dataset of 1,262 clinical and environmental strains collected from 52 countries, including 65 newly sequenced genomes in this study. We established a sublineage threshold based on 133 allelic differences that creates clusters nearly identical to traditional MLST types, providing backwards compatibility to new cgMLST classifications. We also defined an outbreak threshold based on seven allelic differences that is capable of identifying strains from the same outbreak and closely related isolates which could give clues on outbreak origin. Using cgMLST, we confirmed the South Asian origin of modern epidemics and identified clustering affinity among sublineages of environmental isolates from the same geographic origin. Advantages of this method are highlighted by direct comparison with existing classification methods, such as MLST and single nucleotide polymorphism-based methods. cgMLST outperforms all existing methods in terms of resolution, standardization, and ease-of-use. We anticipate this scheme will serve as a basis for a universally applicable and standardized classification system for V. cholerae research and epidemiological surveillance in the future. This cgMLST scheme is publicly available on PubMLST (https://pubmlst.org/vcholerae/).IMPORTANCEToxigenic Vibrio cholerae of the O1 and O139 serogroups are the causative agents of cholera, an acute diarrheal disease that plagued the world for centuries, if not millennia. Here, we introduce a core genome multilocus sequence typing scheme for V. cholerae Using this scheme, we have standardized the definition for subspecies-level classification, facilitating global collaboration in the surveillance of V. cholerae In addition, this typing scheme allows for quick identification of outbreak-related isolates that can guide subsequent analyses, serving as an important first step in epidemiological research. This scheme is also easily scalable to analyze thousands of isolates at various levels of resolution making it an invaluable tool for large-scale ecological and evolutionary analyses.

Invasive Cronobacter infections among infants are associated with severe neurologic disabilities and death. Early Cronobacter reports typically featured hospitalized and preterm infants and recognized contaminated powdered infant formula (PIF) as a transmission vehicle. To clarify recent epidemiology, we reviewed all cases of bloodstream infection or meningitis among infants that were reported to the Centers for Disease Control and Prevention and in the literature (1961–2018; n = 183). Most infants were neonates (100/150 [67%]); 38% (42/112) died, and 79% (81/102) had reported recent PIF consumption. In the final quarter of the study period (2004–2018), case counts were significantly higher (global average 8.7 cases/year); among US cases, significantly higher proportions occurred among full-term (56% [27/48]) and nonhospitalized (78% [42/54]) infants. PIF contamination, most commonly from opened containers, was identified in 30% (21/71) of investigations. Our findings reaffirm the need to promote safer alternatives for infant feeding, particularly among neonates. | Contaminated powdered infant formula from opened containers is the most commonly identified transmission vehicle. | eng

Yersinia enterocolitica (YE) bioserotype 1B/O:8 (YE 1B/O:8) was identified in routine culture of a variety of zoo species housed at Omaha's Henry Doorly Zoo and Aquarium (OHDZA) from April to July 2011. Animal cases representing 12 species had YE detected from 34 cases during routine fecal monitoring and/or during postmortem examination: Coquerel's sifakas (Propithecus coquereli, two cases), black & white (BW) ruffed lemurs (Varecia variegata variegata, six cases), red ruffed lemurs (Varecia rubra, seven cases), white handed gibbon (Hylobates lar albimana, one case), black lemurs (Eulemur macaco, three cases), mongoose lemurs (Eulemur mongoz, two cases), African hunting dogs (Lycaon pictus, five cases), agile gibbons (Hylobates agilis, three cases), siamangs (Hylobates syndactylus, two cases), colobus monkey (Colobus angolensis palliates, one case), argus pheasant (Argusianus argus, one case), and orangutan (Pongo pygmaeus, one case). Most species were not symptomatic; however, three symptomatic cases in Coquerel's sifakas (two) and a white handed gibbon (one) showed clinical signs of diarrhea and lethargy that resulted in death for the Coquerel's sifakas. One unexpected death also occurred in a BW ruffed lemur. To the authors' knowledge, this is the first report of YE 1B/O:8 in such a large variety of zoo species. The source of the YE could not be identified, prompting the initiation of a diseases surveillance program to prevent further cases for the species that are sensitive to YE. To date, no additional cases have been identified, thus suggesting a single introduction of the YE 1B/O:8 strain into the zoo environment.

Populations of the bacterium Vibrio cholerae consist of dozens of distinct lineages, with primarily (but not exclusively) members of the pandemic generating (PG) lineage capable of causing the diarrheal disease cholera. Assessing composition and temporal dynamics of such populations requires extensive isolation efforts and thus only rarely covers large geographic areas or timeframes exhaustively. We developed a culture-independent amplicon sequencing strategy based on the protein-coding gene viuB (vibriobactin utilization) to study the structure of a V. cholerae population over the course of a summer. We show that the 26 co-occurring V. cholerae lineages continuously compete for limited space on nutrient-rich particles where only few of them can grow to large numbers. Differential abundance of lineages between locations and size-fractions associated with a particle-attached or free-swimming lifestyle could reflect adaptation to various environmental niches. In particular, a major V. cholerae lineage occasionally grow to large numbers on particles but remain undetectable using isolation-based methods, indicating selective culturability for some members of the species. We thus demonstrate that isolation-based studies may not accurately reflect the structure and complex dynamics of V. cholerae populations and provide a scalable high-throughput method for both epidemiological and ecological approaches to studying this species. This article is protected by copyright. All rights reserved.

We sequenced the genomes of eight isolates from various regions of the United States. These isolates form a monophyletic cluster clearly related to but distinct from Vibrio cholerae. Phylogenetic and genomic analyses suggest that they represent a basal lineage highly divergent from Vibrio cholerae or a novel species.

We are reporting whole-genome sequences of nine Vibrio sp. isolates closely related to the waterborne human pathogen Vibrio cholerae. These isolates were recovered from sources, including human samples, from different regions of the United States. Genome analysis suggests that this group of isolates represents a highly divergent basal V. cholerae lineage or a closely related novel species.

Latin America has experienced two of the largest cholera epidemics in modern history; one in 1991 and the other in 2010. However, confusion still surrounds the relationships between globally circulating pandemic Vibrio cholerae clones and local bacterial populations. We used whole-genome sequencing to characterize cholera across the Americas over a 40-year time span. We found that both epidemics were the result of intercontinental introductions of seventh pandemic El Tor V. cholerae and that at least seven lineages local to the Americas are associated with disease that differs epidemiologically from epidemic cholera. Our results consolidate historical accounts of pandemic cholera with data to show the importance of local lineages, presenting an integrated view of cholera that is important to the design of future disease control strategies.

The seventh cholera pandemic has heavily affected Africa, although the origin and continental spread of the disease remain undefined. We used genomic data from 1070 Vibrio cholerae O1 isolates, across 45 African countries and over a 49-year period, to show that past epidemics were attributable to a single expanded lineage. This lineage was introduced at least 11 times since 1970, into two main regions, West Africa and East/Southern Africa, causing epidemics that lasted up to 28 years. The last five introductions into Africa, all from Asia, involved multidrug-resistant sublineages that replaced antibiotic-susceptible sublineages after 2000. This phylogenetic framework describes the periodicity of lineage introduction and the stable routes of cholera spread, which should inform the rational design of control measures for cholera in Africa.

These guidelines are intended for use by healthcare professionals who care for children and adults with suspected or confirmed infectious diarrhea. They are not intended to replace physician judgement regarding specific patients or clinical or public health situations. This document does not provide detailed recommendations on infection prevention and control aspects related to infectious diarrhea.

Vibrio sp. strain 2521-89 is an environmental isolate from lake water in New Mexico, USA. Average nucleotide identity, in silico DNA-DNA hybridization, and core genome single-nucleotide polymorphism (SNP)-based phylogenetic analysis suggest that this may be a potentially novel species that is closely related to Vibrio cholerae.

In April 2016, a female infant was born via Cesarean delivery at 29 estimated gestational weeks and had a birthweight of 3 pounds (1,405 grams). Her clinical course in the neonatal intensive care unit was unremarkable until she developed signs of sepsis at age 21 days. Cultures of blood and cerebrospinal fluid yielded Cronobacter sakazakii, a gram-negative pathogenic bacillus. Despite treatment with ampicillin and cefepime, she developed seizures; brain imaging revealed liquefaction necrosis of the entire left cerebral hemisphere and right frontal lobe. The infant developed spastic cerebral palsy and global developmental delay and required a ventriculoperitoneal shunt and a gastrostomy feeding tube. | The infant had been fed pasteurized donor human milk and expressed maternal milk (EMM) during the first week after birth; thereafter, she received EMM mixed with a commercial liquid human milk fortifier. Maternal milk was expressed using a dedicated bedside hospital breast pump and the mother’s personal breast pump throughout the infant’s hospitalization. The infant did not receive powdered infant formula products.

Epidemiological findings of a listeriosis outbreak in 2013 implicated Hispanic-style cheese produced by Company A, and pulsed-field gel electrophoresis (PFGE) and whole genome sequencing (WGS) were performed on clinical isolates and representative isolates collected from Company A cheese and environmental samples during the investigation. The results strengthened the evidence for cheese as the vehicle. Surveillance sampling and WGS three months later revealed that the equipment purchased by Company B from Company A yielded an environmental isolate highly similar to all outbreak isolates. The whole genome and core genome multilocus sequence typing and single nucleotide polymorphism (SNP) analyses were compared to demonstrate the maximum discriminatory power obtained by using multiple analyses, which were needed to differentiate outbreak-associated isolates from a PFGE-indistinguishable isolate collected in a non-implicated food source in 2012. This unrelated isolate differed from the outbreak isolates by only 7 to 14 SNPs, and as a result, minimum spanning tree by the whole genome analyses and certain variant calling approach and phylogenetic algorithm for core genome-based analyses could not provide the differentiation between unrelated isolates. Our data also suggest that SNP/allele counts should always be combined with WGS clustering generated by phylogenetically meaningful algorithms on sufficient number of isolates, and SNP/allele threshold alone is not sufficient evidence to delineate an outbreak. The putative prophages were conserved across all the outbreak isolates. All outbreak isolates belonged to clonal complex 5 and serotype 1/2b, had an identical inlA sequence, which did not have premature stop codons.IMPORTANCE In this outbreak, multiple analytical approaches were used for maximum discriminatory power. A PFGE-matched, epidemiologically unrelated isolate had high genetic similarity to the outbreak-associated isolates, with as few as only 7 SNP differences. Therefore, the SNP/allele threshold should not be used as the only evidence to define the scope of an outbreak. It is critical that the SNP/allele counts be complemented by WGS clustering generated by phylogenetically meaningful algorithms to distinguish outbreak-associated isolates from epidemiologically unrelated isolates. Careful selection of a variant calling approach and phylogenetic algorithm is critical for core genome-based analyses. The whole genome-based analyses were able to construct the highly resolved phylogeny needed to support the findings of the outbreak investigation. Ultimately, epidemiologic evidence and multiple WGS analyses should be combined to increase the confidence in outbreak investigations.

Whole apples have not been previously implicated in outbreaks of foodborne bacterial illness. We investigated a nationwide listeriosis outbreak associated with caramel apples. We defined an outbreak-associated case as an infection with one or both of two outbreak strains of Listeria monocytogenes highly related by whole-genome multilocus sequence typing (wgMLST) from 1 October 2014 to 1 February 2015. Single-interviewer open-ended interviews identified the source. Outbreak-associated cases were compared with non-outbreak-associated cases and traceback and environmental investigations were performed. We identified 35 outbreak-associated cases in 12 states; 34 (97%) were hospitalized and seven (20%) died. Outbreak-associated ill persons were more likely to have eaten commercially produced, prepackaged caramel apples (odds ratio 326.7, 95% confidence interval 32.2-3314). Environmental samples from the grower's packing facility and distribution-chain whole apples yielded isolates highly related to outbreak isolates by wgMLST. This outbreak highlights the importance of minimizing produce contamination with L. monocytogenes. Investigators should perform single-interviewer open-ended interviews when a food is not readily identified.

Diagnostic testing for foodborne pathogens relies on culture-based techniques that are not rapid enough for real-time disease surveillance and do not give a quantitative picture of pathogen abundance or the response of the natural microbiome. Powerful sequence-based, culture-independent approaches such as shotgun metagenomics could sidestep these limitations, and potentially reveal a pathogen-specific signature on the microbiome that would have implications not only for diagnostics but also for better understanding disease progression and pathogen ecology. However, metagenomics have not yet been validated for foodborne pathogen detection. Toward closing these gaps, we applied shotgun metagenomics to stool samples collected from two geographically isolated (Alabama and Colorado) foodborne outbreaks, where the etiologic agents were identified as distinct strains of Salmonella enterica serovar Heidelberg by culture-dependent methods. Metagenomic investigations were consistent with the culture-based findings and revealed, in addition, the in-situ abundance and level of intra-population diversity of the pathogen, the possibility for co-infections with Staphylococcus aureus, and significant shifts in the gut microbiome during infection relative to reference healthy samples. Additionally, we designed our bioinformatics pipeline to deal with several challenges associated with analysis of clinical samples such as the high frequency of co-eluting human DNA sequences and assessment of the virulence potential of pathogens. Comparisons of these results to those of other studies revealed that in several cases of diarrheal outbreaks -but not all- the disease and healthy states of the gut microbial community might be distinguishable, opening new possibilities for diagnostics. IMPORTANCE STATEMENT: Diagnostic testing for enteric pathogens has relied for decades on culture-based techniques but a total of 38.4 million cases of foodborne illness per year cannot be attributed to specific causes. This study describes new culture-independent metagenomic approaches and the associated bioinformatics approaches to detect and type the causative agents of microbial disease with unprecedented accuracy, opening new possibilities for future development of health technologies and diagnostics. Our tools and approaches should be applicable to other microbial diseases in addition to foodborne diarrhea.

Listeria monocytogenes (Lm) is a major human foodborne pathogen. Numerous Lm outbreaks have been reported worldwide and associated with a high case fatality rate, reinforcing the need for strongly coordinated surveillance and outbreak control. We developed a universally applicable genome-wide strain genotyping approach and investigated the population diversity of Lm using 1,696 isolates from diverse sources and geographical locations. We define, with unprecedented precision, the population structure of Lm, demonstrate the occurrence of international circulation of strains and reveal the extent of heterogeneity in virulence and stress resistance genomic features among clinical and food isolates. Using historical isolates, we show that the evolutionary rate of Lm from lineage I and lineage II is low ( approximately 2.5 x 10-7 substitutions per site per year, as inferred from the core genome) and that major sublineages (corresponding to so-called 'epidemic clones') are estimated to be at least 50-150 years old. This work demonstrates the urgent need to monitor Lm strains at the global level and provides the unified approach needed for global harmonization of Lm genome-based typing and population biology.

Four Vibrio spp. isolates from the historical culture collection at the Centers for Disease Control and Prevention, obtained from human blood specimens (n = 3) and river water (n = 1), show characteristics distinct from those of isolates of the most closely related species, Vibrio navarrensis and Vibrio vulnificus, based on phenotypic and genotypic tests. They are specifically adapted to survival in both freshwater and seawater, being able to grow in rich media without added salts as well as salinities above that of seawater. Phenotypically, these isolates resemble V. navarrensis, their closest known relative with a validly published name, but the group of isolates is distinguished from V. navarrensis by the ability to utilize L-rhamnose. Average nucleotide identity and percent DNA-DNA hybridization values obtained from the pairwise comparisons of whole genome sequences of these isolates to V. navarrensis range from 95.4-95.8% and 61.9-64.3%, respectively, suggesting that the group represents a different species. Phylogenetic analysis of the core genome, including four protein-coding housekeeping genes (pyrH, recA, rpoA, and rpoB), places these four isolates into their own monophyletic clade, distinct from V. navarrensis and V. vulnificus. Based on these differences, we propose these isolates belong to a novel Vibrio species. The name Vibrio cidicii sp. nov. is proposed for these isolates; strain LMG 29267T (= CIP 111013T = 2756-81T), isolated from river water, is the type strain.

We report on a case of listeriosis in a patient who probably consumed a prepackaged romaine lettuce-containing product recalled for Listeria monocytogenes contamination. Although definitive epidemiological information demonstrating exposure to the specific recalled product was lacking, the patient reported consumption of a prepackaged romaine lettuce-containing product of either the recalled brand or a different brand. A multinational investigation found that patient and food isolates from the recalled product were indistinguishable by pulsed-field gel electrophoresis and were highly related by whole genome sequencing, differing by four alleles by whole genome multilocus sequence typing and by five high-quality single nucleotide polymorphisms, suggesting a common source. To our knowledge, this is the first time prepackaged lettuce has been identified as a likely source for listeriosis. This investigation highlights the power of whole genome sequencing, as well as the continued need for timely and thorough epidemiological exposure data to identify sources of foodborne infections.

Listeria monocytogenes(Lm) causes severe foodborne illness (listeriosis). Previous molecular subtyping methods, such as pulsed-field gel electrophoresis (PFGE), were critical in detecting outbreaks that led to food safety improvements and declining incidence, but PFGE provides limited genetic resolution. A multiagency collaboration began performing real-time, whole-genome sequencing (WGS) on all U.S.Lmisolates from patients, food, and the environment in September 2013, posting sequencing data into a public repository. Compared with the year before the project began, WGS, combined with epidemiologic and product trace-back data, detected more listeriosis clusters and solved more outbreaks (2 outbreaks in pre-WGS year, 5 in WGS year 1, and 9 in year 2). Whole-genome multilocus sequence typing and single nucleotide polymorphism analyses provided equivalent phylogenetic relationships relevant to investigations; results were most useful when interpreted in context of epidemiological data. WGS has transformed listeriosis outbreak surveillance and is being implemented for other foodborne pathogens.

The bacterial pathogen Vibrio cholerae can occupy both the human gut and aquatic reservoirs, where it may colonize chitinous surfaces that induce expression of factors for three phenotypes: chitin utilization, DNA uptake by natural transformation, and constitutive contact-dependent bacterial killing indicative of the Type VI secretion system (T6SS). In this study, we surveyed a diverse set of 53 isolates from different geographic locales collected over the past century from human clinical and environmental specimens for each phenotype outlined above. The set included pandemic isolates of serogroup O1, as well as several serogroup O139 and non-O1/non-O139 strains. We found that while chitin utilization was common, only 22.6% of all isolates tested were proficient at chitin-induced natural transformation, suggesting that transformation is expendable. Constitutive contact-dependent killing of Escherichia coli prey, which is indicative of a functional T6SS, was rare among clinical isolates (only 4 of 29) but common among environmental isolates (22 of 24). These results bolster the pathoadaptive model that tight regulation of T6SS-mediated bacterial killing is beneficial in a human host, whereas, constitutive killing by environmental isolates may give a competitive advantage in natural settings. Future sequence analysis of this set of diverse isolates may identify previously unknown regulators and structural components for both natural transformation and T6SS.

Listeriosis is a serious foodborne infection that disproportionately affects elderly adults, pregnant women, newborns, and immunocompromised individuals. Diagnosis is made by culturing Listeria monocytogenes from sterile body fluids or products of conception. This report describes investigations of two listeriosis pseudo-outbreaks caused by contaminated laboratory media made from sheep blood.

We used whole genome sequencing to determine evolutionary relationships among 20 outbreak-associated clinical isolates of Listeria monocytogenes serotypes 1/2a and 1/2b. Isolates from six of eleven outbreaks fell outside of the clonal groups or 'epidemic clones' that have been previously associated with outbreaks, suggesting that epidemic potential may be widespread in L. monocytogenes and is not limited to the recognized epidemic clones. Pairwise comparisons between epidemiologically-related isolates within clonal complexes showed that genome-level variation differed by two orders of magnitude between different comparisons, and the distribution of point mutations (core versus accessory genome) also varied. In addition, genetic divergence between one closely related pair of isolates from a single outbreak was driven primarily by changes in phage regions. The evolutionary analysis showed the changes could be attributed to horizontal gene transfer; members of the diverse bacterial community found in the production facility could have served as the source of novel genetic material at some point in the production chain. The results raise the question of how to best utilize information contained within the accessory genome in outbreak investigations. The full magnitude and complexity of genetic changes revealed by genome sequencing could not be discerned from traditional subtyping methods and the results demonstrate the challenges of interpreting genetic variation among isolates recovered from a single outbreak. Epidemiological information remains critical for proper interpretation of nucleotide and structural diversity among isolates recovered during outbreaks, and will remain so until we understand more about how various population histories influence genetic variation.

Vibrio metoecus is the closest relative of Vibrio cholerae, the causative agent of the potent diarrheal disease cholera. Although the pathogenic potential of this new species is yet to be studied in depth, it has been co-isolated with V. cholerae in coastal waters and found in clinical specimens in the USA. We used these two organisms to investigate the genetic interaction between closely related species in their natural environment. The genomes of twenty V. cholerae and four V. metoecus strains isolated from a brackish coastal pond on the US east coast, as well as four clinical V. metoecus strains were sequenced and compared to reference strains. Whole genome comparison shows 86-87% average nucleotide identity (ANI) in their core genes between the two species. On the other hand, the chromosomal integron, which occupies approximately 3% of their genomes, shows higher conservation in ANI between species than any other region of their genomes. The ANI of 93-94% observed in this region is not significantly greater within than between species, meaning that it does not follow species boundaries. V. metoecus does not encode toxigenic V. cholerae major virulence factors, the cholera toxin and toxin co-regulated pilus. However, some of the pathogenicity islands found in pandemic V. cholerae were either present in the common ancestor it shares with V. metoecus, or acquired by clinical and environmental V. metoecus in partial fragments. The virulence factors of V. cholerae are therefore both more ancient and more widespread than previously believed. There is high interspecies recombination in the core genome, which has been detected in 24% of the single-copy core genes, including genes involved in pathogenicity. V. metoecus was six times more often the recipient of DNA from V. cholerae as it was the donor, indicating a strong bias in the direction of gene transfer in the environment.