Insook Son

High resolution clustering of Salmonella enterica serovar Montevideo strains using a next-generation sequencing approach

BackgroundNext-Generation Sequencing (NGS) is increasingly being used as a molecular epidemiologic tool for discerning ancestry and traceback of the most complicated, difficult to resolve bacterial pathogens. Making a linkage between possible food sources and clinical isolates requires distinguishing the suspected pathogen from an environmental background and placing the variation observed into the wider context of variation occurring within a serovar and among other closely related foodborne pathogens. Equally important is the need to validate these high resolution molecular tools for use in molecular epidemiologic traceback. Such efforts include the examination of strain cluster stability as well as the cumulative genetic effects of sub-culturing on these clusters. Numerous isolates of S. Montevideo were shot-gun sequenced including diverse lineage representatives as well as numerous replicate clones to determine how much variability is due to bias, sequencing error, and or the culturing of isolates. All new draft genomes were compared to 34 S. Montevideo isolates previously published during an NGS-based molecular epidemiological case study.ResultsIntraserovar lineages of S. Montevideo differ by thousands of SNPs, that are only slightly less than the number of SNPs observed between S. Montevideo and other distinct serovars. Much less variability was discovered within an individual S. Montevideo clade implicated in a recent foodborne outbreak as well as among individual NGS replicates. These findings were similar to previous reports documenting homopolymeric and deletion error rates with the Roche 454 GS Titanium technology. In no case, however, did variability associated with sequencing methods or sample preparations create inconsistencies with our current phylogenetic results or the subsequent molecular epidemiological evidence gleaned from these data.ConclusionsImplementation of a validated pipeline for NGS data acquisition and analysis provides highly reproducible results that are stable and predictable for molecular epidemiological applications. When draft genomes are collected at 15×-20× coverage and passed through a quality filter as part of a data analysis pipeline, including sub-passaged replicates defined by a few SNPs, they can be accurately placed in a phylogenetic context. This reproducibility applies to all levels within and between serovars of Salmonella suggesting that investigators using these methods can have confidence in their conclusions.

Genetic diversity of Arcobacter and Campylobacter on broiler carcasses during processing.

Broiler carcasses (n=325) were sampled at three sites along the processing line (prescalding, prechilling, and postchilling) in a commercial poultry processing plant during five plant visits from August to October 2004. Pulsed-field gel electrophoresis (PFGE) was used to determine the genomic fingerprints of Camospylobacter coli (n=27), Campylobacter jejuni (n=188), Arcobacter butzleri (n=138), Arcobacter cryaerophilus 1A (n=4), and A. cryaerophilus 1B (n=31) with the restriction enzymes SmaI and KpnI for Campylobacter and Arcobacter, respectively. Campylobacter species were subtyped by the Centers for Disease Control and Prevention PulseNet 24-h standardized protocol for C. jejuni. A modification of this protocol with a different restriction endonuclease (KpnI) and different electrophoresis running conditions produced the best separation of restriction fragment patterns for Arcobacter species. Both unique and common PFGE types of Arcobacter and Campylobacter strains were identified. A total of 32.8% (57 of 174) of the Arcobacter isolates had unique PFGE profiles, whereas only 2.3% (5 of 215) of the Campylobacter isolates belonged to this category. The remaining Arcobacter strains were distributed among 25 common PFGE types; only eight common Campylobacter PFGE types were observed. Cluster analysis showed no associations among the common PFGE types for either genus. Each of the eight common Campylobacter types consisted entirely of isolates from one sampling day, whereas more than half of the common Arcobacter types contained isolates from different sampling days. Our results demonstrate far greater genetic diversity for Arcobacter than for Campylobacter and suggest that the Campylobacter types are specific to individual flocks of birds processed on each sampling day.

Prevalence of Hemolysin Genes and Comparison of ehxA Subtype Patterns in Shiga toxin-producing Escherichia coli (STEC) and non-STEC Strains from Clinical, Food, and Animal Sources

ABSTRACT Shiga toxin-producing Escherichia coli (STEC) belonging to certain serogroups (e.g., O157 and O26) can cause serious conditions like hemolytic-uremic syndrome (HUS), but other strains might be equally pathogenic. While virulence factors, like stx and eae, have been well studied, little is known about the prevalence of the E. coli hemolysin genes (hlyA, ehxA, e-hlyA, and sheA) in association with these factors. Hemolysins are potential virulence factors, and ehxA and hlyA have been associated with human illness, but the significance of sheA is unknown. Hence, 435 E. coli strains belonging to 62 different O serogroups were characterized to investigate gene presence and phenotypic expression of hemolysis. We further investigated ehxA subtype patterns in E. coli isolates from clinical, animal, and food sources. While sheA and ehxA were widely distributed, e-hlyA and hlyA were rarely found. Most strains (86.7%) were hemolytic, and significantly more hemolytic (95%) than nonhemolytic strains (49%) carried stx and/or eae (P < 0.0001). ehxA subtyping, as performed by using PCR in combination with restriction fragment length polymorphism analysis, resulted in six closely related subtypes (>94.2%), with subtypes A/D being eae-negative STECs and subtypes B, C, E, and F eae positive. Unexpectedly, ehxA subtype patterns differed significantly between isolates collected from different sources (P < 0.0001), suggesting that simple linear models of exposure and transmission need modification; animal isolates carried mostly subtypes A/C (39.3%/42.9%), food isolates carried mainly subtype A (81.9%), and clinical isolates carried mainly subtype C (66.4%). Certain O serogroups correlated with particular ehxA subtypes: subtype A with O104, O113, and O8; B exclusively with O157; C with O26, O111, and O121.

Comparison of eight different agars for the recovery of clinically relevant non-O157 Shiga toxin-producing Escherichia coli from baby spinach, cilantro, alfalfa sprouts and raw milk.

The FDA Bacteriological Analytical Manual (BAM) Chapter 4a recommends several agars for isolating non-O157 Shiga toxin-producing Escherichia coli (STEC); not all have been thoroughly tested for recovering STECs from food. Using E. coli strains representing ten clinically relevant O serogroups (O26, O45, O91, O103, O104, O111, O113, O121, O128, O145) in artificially-contaminated fresh produce--bagged baby spinach, alfalfa sprouts, cilantro, and raw milk--we evaluated the performance of 8 different agars. Performance was highly dependent upon strain used and the presence of inhibitors, but not necessarily dependent on food matrix. Tellurite resistant-negative strains, O91:-, O103:H6, O104:H21, O113:H21, and O128, grew poorly on CHROMagar STEC, Rainbow agar O157, and a modified Rainbow O157 (mRB) agar. Although adding washed sheeps blood to CHROMagar STEC and mRB agars improved overall performance; however, this also reversed the inhibition of non-target bacteria provided by original formulations. Variable colony coloration made selecting colonies from Rainbow agar O157 and mRB agars difficult. Study results support a strategy using inclusive agars (e.g. L-EMB, SHIBAM) in combination with selective agars (R & F E. coli O157:H7, CHROMagar STEC) to allow for recovery of the most STECs while increasing the probability of recovering STEC in high bacterial count matrices.

Isolation of Shiga toxin-producing Escherichia coli from fresh produce using STEC heart infusion washed blood agar with mitomycin-C.

The ability to detect and isolate Shiga toxin-producing Escherichia coli (STEC) remains a major challenge for food microbiologists. Although methods based on nucleic acids and antibodies have improved detection of STECs in foods, isolation of these bacteria remains arduous. STEC isolation is necessary for matching food, environmental, and clinical isolates during outbreak investigations and for distinguishing between pathogenic and nonpathogenic organisms. STEC heart infusion washed blood agar with mitomycin-C (SHIBAM) is a modification of washed sheep blood agar prepared by adding mitomycin-C and optimizing both the washed blood and base agar to better isolate STECs. Most STEC isolates produce a zone of hemolysis on SHIBAM plates and are easily distinguishable from background microbiota. Here, we present data supporting the use of SHIBAM to isolate STECs from fresh produce. SHIBAM was tested for accuracy in identifying STECs (365 of 410 STEC strains were hemolytic, and 63 of 73 E. coli strains that did not produce Shiga toxin were not hemolytic) and for recovery from artificially inoculated fresh produce (11 of 24 romaine lettuce samples and 6 of 24 tomato samples). STEC recovery with SHIBAM agar was greatly improved when compared with recovery on Levines eosin-methylene blue agar as a reference method.

Modeling the survival kinetics of Salmonella in tree nuts for use in risk assessment

Salmonella has been shown to survive in tree nuts over long periods of time. This survival capacity and its variability are key elements for risk assessment of Salmonella in tree nuts. The aim of this study was to develop a mathematical model to predict survival of Salmonella in tree nuts at ambient storage temperatures that considers variability and uncertainty separately and can easily be incorporated into a risk assessment model. Data on Salmonella survival on raw almonds, pecans, pistachios and walnuts were collected from the peer reviewed literature. The Weibull model was chosen as the baseline model and various fixed effect and mixed effect models were fit to the data. The best model identified through statistical analysis testing was then used to develop a hierarchical Bayesian model. Salmonella in tree nuts showed slow declines at temperatures ranging from 21°C to 24°C. A high degree of variability in survival was observed across tree nut studies reported in the literature. Statistical analysis results indicated that the best applicable model was a mixed effect model that included a fixed and random variation of δ per tree nut (which is the time it takes for the first log10 reduction) and a fixed variation of ρ per tree nut (parameter which defines the shape of the curve). Higher estimated survival rates (δ) were obtained for Salmonella on pistachios, followed in decreasing order by pecans, almonds and walnuts. The posterior distributions obtained from Bayesian inference were used to estimate the variability in the log10 decrease levels in survival for each tree nut, and the uncertainty of these estimates. These modeled uncertainty and variability distributions of the estimates can be used to obtain a complete exposure assessment of Salmonella in tree nuts when including time-temperature parameters for storage and consumption data. The statistical approach presented in this study may be applied to any studies that aim to develop predictive models to be implemented in a probabilistic exposure assessment or a quantitative microbial risk assessment.

Multilaboratory validation of a luminex microbead-based suspension array for the identification of the 11 most clinically relevant Shiga toxin-producing Escherichia coli O serogroups.

Rapid and high-throughput identification and serotyping of Shiga toxin-producing Escherichia coli (STEC) O serogroups is important for detecting, investigating, and controlling STEC infection outbreaks and removing hazardous products from commerce. A Luminex microbead-based suspension array has been developed to identify the 11 most clinically relevant STEC serogroups: O26, O45, O91, O103, O104, O111, O113, O121, O128, O145, and O157. Here we present results of a blinded multilaboratory collaborative study involving 10 participants from nine laboratories using 55 unknown strains. From the total 495 analyses, two false-positive and three false-negative results were obtained, indicating the assay to be a rapid, high-throughput, and robust method for identifying clinically relevant STEC serogroups.

Comparison of different sample preparation procedures for the detection and isolation of Escherichia coli O157:H7 and Non-O157 STECs from leafy greens and cilantro

The FDA Bacteriological Analytical Manual (BAM) method for the detection/isolation of Shiga toxin-producing Escherichia coli (STEC) involves enrichment of produce rinses, blended homogenates or stomached homogenates. However, the effectiveness of rinsing produce to remove attached bacteria is largely unknown. Moreover, PCR inhibitors can be released under physical treatment. The study objective was to determine the relative effectiveness of recovery methods for STEC contaminated produce. Spinach, lettuce, and cilantro were contaminated with E. coli O157:H7 or a non-O157 STEC, subjected to both the BAM method and a soak method, and tested by real-time PCR and cultural methods. For O157:H7 and non-O157:H7 STECs, the soak method was significantly more productive than leafy green rinses. Of 320 test portions, PCR of recovered colonies confirmed 148 were positive by rinsing and 271 were positive by soaking (an 83% increase in sensitivity). For recovery of O157:H7 from cilantro, of 60 test portions, positives were 38 by soaking, 41 by stomaching, and 28 by blending. Soaking and stomaching were significantly more productive than blending, although soaking was only arithmetically superior to stomaching. Based upon these results, it is recommended that a soak method replace the current BAM procedures.

Genotyping Campylobacter jejuni by comparative genome indexing: an evaluation with pulsed-field gel electrophoresis and flaA SVR sequencing.

Comparative genome indexing (CGI) using whole-genome DNA microarrays was evaluated as a means of genotyping Campylobacter jejuni relative to two standard methods, pulsed-field gel electrophoresis (PFGE) and flaA short variable region sequencing (flaA SVR typing). Thirty-six geographically diverse C. jejuni isolates were selected from a collection of cattle and chicken isolates. The BioNumerics software program was used for cluster analysis of the data from all 36 isolates for each of the three typing methods. Comparative genome indexing assigned a unique type to each isolate while PFGE and flaA SVR distinguished 29 and 35 different types, respectively. The four common types identified by PFGE were also closely related by CGI, and the overall similarity of the CGI results to those for PFGE indicates the value of CGI as a more informative alternative to PFGE. While flaA SVR was very discriminative, the isolates were all highly similar (>78%) resulting in finer distinctions between isolates and fewer genotypic relations to CGI or PFGE. Campylobacter jejuni is one of the most common causative agents of bacterial gastroenteritis in the world. The development of CGI as a molecular typing tool for C. jejuni offers a highly effective and informative means of further understanding the epidemiology of this ubiquitous pathogen.

A Quantitative Assessment of the Risk of Human Salmonellosis Arising from the Consumption of Almonds in the United States: The Impact of Preventive Treatment Levels

The presence of Salmonella on almonds continues to result in product-related outbreaks and recalls in the United States. In this study, the impact of microbial reduction treatment levels (1 to 5 log CFU) on the risk of human salmonellosis from the consumption of almond kernels in the United States was evaluated. An exposure model, including major steps in almond processing, was used to estimate prevalence and levels of contamination of Salmonella on almonds at the point of consumption. A Salmonella dose-response model and consumption data for almonds in the United States were used to assess risk of illness per serving and per year, quantifying variability and uncertainty separately. A 3-log reduction treatment resulted in a predicted mean risk of illness of two cases per year for almonds consumed as a core product not cooked at home (95% confidence interval [CI], one to four cases), one case per year for almonds consumed as an ingredient not cooked at home (95% CI, one to two cases), and less than one case per year for almonds consumed as an ingredient cooked at home (95% CI, 7 × 10-7 to 3 × 10-6 cases). A minimum 4-log reduction treatment resulted in an estimated mean risk of illness below one case per year in the United States. This study also includes an assessment of the risk of human salmonellosis as a result of an exceptional situation, which results in higher risk estimates compared with the baseline model. The exceptional situations modeled posttreatment resulted in estimates of mean risk that were not significantly affected by treatment level. Sensitivity analysis results showed initial Salmonella contamination level to be the factor with the most impact on risk per serving estimates, given a certain treatment level. The risk assessment also includes a simulation of the events that occurred in 2001. Treatment levels with a minimum 4-log microbial reduction would have been sufficient to prevent the outbreak cases. The uncertainty range in the estimates indicates that additional information is needed to make more precise predictions of this specific outbreak event.