Whole Genome and Core Genome Multilocus Sequence Typing and Single Nucleotide Polymorphism Analyses of Listeria monocytogenes Isolates Associated with an Outbreak Linked to Cheese, United States, 2013.

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 3 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 results 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 nonimplicated food source in 2012. This unrelated isolate differed from the outbreak isolates by only 7 to 14 SNPs, and as a result, the minimum spanning tree from the whole genome analyses and certain variant calling approach and phylogenetic algorithm for core genome-based analyses could not provide differentiation between unrelated isolates. Our data also suggest that SNP/allele counts should always be combined with WGS clustering analysis generated by phylogenetically meaningful algorithms on a sufficient number of isolates, and the SNP/allele threshold alone does not provide 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 and 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 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 analysis 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 confidence levels during outbreak investigations.

Listeria monocytogenes in Stone Fruits Linked to a Multistate Outbreak: Enumeration of Cells and Whole-Genome Sequencing.

In 2014, the identification of stone fruits contaminated with Listeria monocytogenes led to the subsequent identification of a multistate outbreak. Simultaneous detection and enumeration of L. monocytogenes were performed on 105 fruits, each weighing 127 to 145 g, collected from 7 contaminated lots. The results showed that 53.3% of the fruits yielded L. monocytogenes (lower limit of detection, 5 CFU/fruit), and the levels ranged from 5 to 2,850 CFU/fruit, with a geometric mean of 11.3 CFU/fruit (0.1 CFU/g of fruit). Two serotypes, IVb-v1 and 1/2b, were identified by a combination of PCR- and antiserum-based serotyping among isolates from fruits and their packing environment; certain fruits contained a mixture of both serotypes. Single nucleotide polymorphism (SNP)-based whole-genome sequencing (WGS) analysis clustered isolates from two case-patients with the serotype IVb-v1 isolates and distinguished outbreak-associated isolates from pulsed-field gel electrophoresis (PFGE)-matched, but epidemiologically unrelated, clinical isolates. The outbreak-associated isolates differed by up to 42 SNPs. All but one serotype 1/2b isolate formed another WGS cluster and differed by up to 17 SNPs. Fully closed genomes of isolates from the stone fruits were used as references to maximize the resolution and to increase our confidence in prophage analysis. Putative prophages were conserved among isolates of each WGS cluster. All serotype IVb-v1 isolates belonged to singleton sequence type 382 (ST382); all but one serotype 1/2b isolate belonged to clonal complex 5. IMPORTANCE: WGS proved to be an excellent tool to assist in the epidemiologic investigation of listeriosis outbreaks. The comparison at the genome level contributed to our understanding of the genetic diversity and variations among isolates involved in an outbreak or isolates associated with food and environmental samples from one facility. Fully closed genomes increased our confidence in the identification and comparison of accessory genomes. The diversity among the outbreak-associated isolates and the inclusion of PFGE-matched, but epidemiologically unrelated, isolates demonstrate the high resolution of WGS. The prevalence and enumeration data could contribute to our further understanding of the risk associated with Listeria monocytogenes contamination, especially among high-risk populations.

Tracing Origins of the Salmonella Bareilly Strain Causing a Food-borne Outbreak in the United States.

BACKGROUND: Using a novel combination of whole-genome sequencing (WGS) analysis and geographic metadata, we traced the origins of Salmonella Bareilly isolates collected in 2012 during a widespread food-borne outbreak in the United States associated with scraped tuna imported from India. METHODS: Using next-generation sequencing, we sequenced the complete genome of 100 Salmonella Bareilly isolates obtained from patients who consumed contaminated product, from natural sources, and from unrelated historically and geographically disparate foods. Pathogen genomes were linked to geography by projecting the phylogeny on a virtual globe and produced a transmission network. RESULTS: Phylogenetic analysis of WGS data revealed a common origin for outbreak strains, indicating that patients in Maryland and New York were infected from sources originating at a facility in India. CONCLUSIONS: These data represent the first report fully integrating WGS analysis with geographic mapping and a novel use of transmission networks. Results showed that WGS vastly improves our ability to delimit the scope and source of bacterial food-borne contamination events. Furthermore, these findings reinforce the extraordinary utility that WGS brings to global outbreak investigation as a greatly enhanced approach to protecting the human food supply chain as well as public health in general.

Prevalence and Level of Listeria monocytogenes in Ice Cream Linked to a Listeriosis Outbreak in the United States.

A most-probable-number (MPN) method was used to enumerate Listeria monocytogenes in 2,320 commercial ice cream scoops manufactured on a production line that was implicated in a 2015 listeriosis outbreak in the United States. The analyzed samples were collected from seven lots produced in November 2014, December 2014, January 2015, and March 2015. L. monocytogenes was detected in 99% (2,307 of 2,320) of the tested samples (lower limit of detection, 0.03 MPN/g), 92% of which were contaminated at <20 MPN/g. The levels of L. monocytogenes in these samples had a geometric mean per lot of 0.15 to 7.1 MPN/g. The prevalence and enumeration data from an unprecedented large number of naturally contaminated ice cream products linked to a listeriosis outbreak provided a unique data set for further understanding the risk associated with L. monocytogenes contamination for highly susceptible populations.

Draft Genome Sequence of Salmonella enterica subsp. enterica Serovar Give, Isolated from an Imported Chili Powder Product.

We report the genome sequence of Salmonella enterica subsp. enterica serovar Give (CFSAN012622), isolated from imported chili powder in 2014. This genome contains genes previously reported to be specific only to S. enterica serovar Enteritidis. This strain shows a unique pulsed-field gel electrophoresis (PFGE) pattern clustering with serovar Enteritidis (JEG X01.0005).

Genome-wide methylation patterns in Salmonella enterica Subsp. enterica Serovars.

The methylation of DNA bases plays an important role in numerous biological processes including development, gene expression, and DNA replication. Salmonella is an important foodborne pathogen, and methylation in Salmonella is implicated in virulence. Using single molecule real-time (SMRT) DNA-sequencing, we sequenced and assembled the complete genomes of eleven Salmonella enterica isolates from nine different serovars, and analysed the whole-genome methylation patterns of each genome. We describe 16 distinct N6-methyladenine (m6A) methylated motifs, one N4-methylcytosine (m4C) motif, and one combined m6A-m4C motif. Eight of these motifs are novel, i.e., they have not been previously described. We also identified the methyltransferases (MTases) associated with 13 of the motifs. Some motifs are conserved across all Salmonella serovars tested, while others were found only in a subset of serovars. Eight of the nine serovars contained a unique methylated motif that was not found in any other serovar (most of these motifs were part of Type I restriction modification systems), indicating the high diversity of methylation patterns present in Salmonella.

Complete Genome Sequences of Salmonella enterica Serovar Heidelberg Strains Associated with a Multistate Food-Borne Illness Investigation.

Next-generation sequencing is being evaluated for use with food-borne illness investigations, especially when the outbreak strains produce patterns that cannot be discriminated from non-outbreak strains using conventional procedures. Here we report complete genome assemblies of two Salmonella enterica serovar Heidelberg strains with a common pulsed-field gel electrophoresis pattern isolated during an outbreak investigation.

Draft Genome Sequences of Two O104:H21 Escherichia coli Isolates Causing Hemorrhagic Colitis during a 1994 Montana Outbreak Provide Insight into Their Pathogenicity.

We sequenced the genomes of two strains of O104:H21 enterohemorrhagic Escherichia coli (EHEC) isolated during an outbreak of hemorrhagic colitis in Montana in 1994. These strains carried a plasmid that contains several virulence genes not present in pO157. The genome sequences will improve phylogenetic analysis of other non-O157 E. coli strains in the future.

Association between longer duration of HIV-suppressive therapy and partial recovery of the V gamma 2 T cell receptor repertoire.

Human immunodeficiency virus (HIV) infection alters the function and repertoire of peripheral blood gamma delta T cells expressing the V gamma 2/V delta 2 T cell receptor (TCR). A cross-sectional comparison of the V gamma 2 TCR repertoire was performed for 56 HIV-infected subjects, 51 of whom were receiving highly active antiretroviral therapy (HAART), to measure changes in the V gamma 2 TCR repertoire. Longer durations of HAART were associated with partial recovery of the normal TCR repertoire, and recovery was greatest in subjects with complete virus suppression. Changes in the V gamma 2 TCR repertoire are sensitive measures for the effects of HAART and of residual HIV replication.

Effects of virus burden and chemokine expression on immunity to SHIV in nonhuman primates.

HIV-1 vaccine candidates are designed to elicit Type 1 immune responses, including cytotoxic T cells and neutralizing antibodies. The type of immune response is influenced by many factors, including the levels of antigen expression and production of cytokines or chemokines; we designed a nonhuman primate study to evaluate the influence of these factors on protective immunity. Recombinant SHIV were engineered to express macrophage inflammatory protein-1 alpha (MIP-1alpha), regulated upon activation, normal T-cell expressed and secreted (RANTES), or Lymphotactin (Ltn) in place of nef in SHIV(89.6) (SHIV(89.6-MIP-1), SHIV(89.6-RANTES), SHIV(89.6-Ltn)). The parental virus SHIV(89.6) was included because it replicates to higher titer while still not causing disease. Control groups included animals that received a recombinant SHIV with a truncated chemokine construct (SHIV(89.6-dLtn)) and unvaccinated macaques. After pathogenic challenge with SHIV(89.6pd), animals from groups that received recombinant (nef-deleted) viruses had peak viremia levels three orders of magnitude lower than unvaccinated controls and increased survival times. Animals that received the original SHIV(89.6) (nef+) were highly resistant to both intrarectal and intravenous challenge with SHIV(89.6PD), and showed no signs of disease. There were no differences in survival times comparing unvaccinated and SHIV(89.6-dLtn) (control) groups, indicating that nef deleted viruses did not provide durable protection in this model. Strongest protection was seen in animals with the highest replicating virus (SHIV(89.6)), and the lower effect on survival after SHIV(89.6) nef-deleted vaccination, likely reflects differences in replication capacity. The protective effect of nef-deleted virus was partly restored by expressing Type 1 chemokines to augment viral immunity.

HIV-mediated gammadelta T cell depletion is specific for Vgamma2+ cells expressing the Jgamma1.2 segment.

Circulating Vgamma2/Vdelta2(+) T cells, normally constituting 3-6% of all CD3(+) T cells in blood, are severely depleted after HIV infection. The mechanism(s) for Vgamma2/Vdelta2(+) T cell depletion are unknown, partly because these cells are CD4(-) and resistant to HIV infection. To determine whether this cell depletion was general for all Vgamma2(+) cells or specific for an individual subset, we analyzed the Vgamma2 repertoire and found consistent differences between HIV(+) and uninfected control samples. The change in Vgamma2 repertoire was the result of preferentially depleting only those Vgamma2 cells that express the Jgamma1.2 segment. The specific loss of Vgamma2-Jgamma1.2(+) cells was polyclonal, as the Vgamma subset retained normal diversity even after HIV infection, and loss occurred without significant changes in the paired chain (Vdelta2) repertoire, or in the alternate Vdelta1 chain repertoire. Specific depletion of Vgamma2-Jgamma1.2/Vdelta2 T cells is the first evidence of a common, T cell receptor-dependent cell loss in HIV disease and it provides a clear example of bystander cell depletion. Vgamma2-Jgamma1.2/Vdelta2 T cells mediate potent responses to microbial pathogens including HIV, and loss of this subset is an important aspect of AIDS pathogenesis.

Strains of Escherichia coli O157:H7 differ primarily by insertions or deletions, not single-nucleotide polymorphisms.

Escherichia coli O157:H7 (O157) strains demonstrate varied pulsed-field gel electrophoresis patterns following XbaI digestion, which enable epidemiological surveillance of this important human pathogen. The genetic events underlying PFGE differences between strains, however, are not defined. We investigated the mechanisms for strain variation in O157 by recovering and examining nucleotide sequences flanking each of the XbaI restriction enzyme sites in the genome. Our analysis demonstrated that differences between O157 strains were due to discrete insertions or deletions that contained the XbaI sites polymorphic between strains rather than single-nucleotide polymorphisms in the XbaI sites themselves. These insertions and deletions were found to be uniquely localized within the regions of the genome that are specific to O157 compared to E. coli K-12 (O islands), suggesting that strain-to-strain variation occurs in these O islands. These results may be utilized to devise novel strain-typing tools for this pathogen.

Polymorphic amplified typing sequences provide a novel approach to Escherichia coli O157:H7 strain typing.

Escherichia coli O157:H7 (O157) strains are commonly typed by pulsed-field gel electrophoresis (PFGE) following digestion of genomic DNA with the restriction enzyme XbaI. We have shown that O157 strains differ from each other by a series of discrete insertions or deletions, some of which contain recognition sites for XbaI, suggesting that these insertions and deletions are responsible for the differences in PFGE patterns. We have devised a new O157 strain typing protocol, polymorphic amplified typing sequences (PATS), based on this information. We designed PCR primer pairs to amplify genomic DNA flanking each of 40 individual XbaI sites in the genomes of two O157 reference strains. These primer pairs were tested with 44 O157 isolates, 2 each from 22 different outbreaks of infection. Thirty-two primer pairs amplified identical fragments from all 44 isolates, while eight primer pairs amplified regions that were polymorphic between isolates. The isolates could be differentiated solely on the basis of which of the eight polymorphic amplicons was detected. PATS correctly identified 21 of 22 outbreak pairs as being identical or highly related, whereas PFGE correctly identified 14 of the 22 outbreak pairs as being identical or highly related; PATS was also able to type isolates from three outbreaks that were untypeable by PFGE. However, PATS was less sensitive than PFGE in discriminating between outbreaks. These data suggest that typing by PATS may provide a simple procedure for strain typing of O157 and other bacteria and that further evaluation of the utility of this method for epidemiologic investigations is warranted.

Urease of enterohemorrhagic Escherichia coli: evidence for regulation by fur and a trans-acting factor.

Recent genomic analyses of Escherichia coli O157:H7 strain EDL933 revealed two loci encoding urease gene homologues (ureDABCEFG), which are absent in nonpathogenic E. coli strain K-12. This report demonstrates that the cloned EDL933 ure gene cluster is capable of synthesizing urease in an E. coli DH5alpha background. However, when the gene fragment is transformed back into the native EDL933 background, the enzymatic activity of the cloned determinants is undetectable. We speculate that an unidentified trans-acting factor in enterohemorrhagic E. coli (EHEC) is responsible for this regulation of ure expression. In addition, Fur-like recognition sites are present in three independent O157:H7 isolates upstream of ureD and ureA. Enzymatic assays confirmed a difference in urease expression of cloned EHEC ure clusters in E. coli MC3100Deltafur. Likewise, interruption of fur in O157:H7 isolate IN1 significantly diminished urease activity. We propose that, similar to the function of Fur in regulating the acid response of Salmonella enterica serovar Typhimurium, it modulates urease expression in EHEC, perhaps contributing to the acid tolerance of the organism.