Search

We sequenced 25 isolates of phenotypically multidrug-resistant Salmonella Indiana (n = 11), Typhimurium (n = 8), and Enteritidis (n = 6) using both MinION long-read [SQK-LSK109 and flow cell (R9.4.1)] and MiSeq short-read (Nextera XT and MiSeq Reagent Kit v2) sequencing technologies to determine the advantages of each approach in terms of the characteristics of genome structure, antimicrobial resistance (AMR), virulence potential, whole-genome phylogeny, and pan-genome. The MinION reads were base-called in real-time using MinKnow 3.4.8 integrated with Guppy 3.0.7. The long-read-only assembly, Illumina-only assembly, and hybrid assembly pipelines of Unicycler 0.4.8 were used to generate the MinION, MiSeq, and hybrid assemblies, respectively. The MinION assemblies were highly contiguous compared to the MiSeq assemblies but lacked accuracy, a deficiency that was mitigated by adding the MiSeq short reads through the Unicycler hybrid assembly which corrected erroneous single nucleotide polymorphisms (SNPs). The MinION assemblies provided similar predictions of AMR and virulence potential compared to the MiSeq and hybrid assemblies, although they produced more total false negatives of AMR genotypes, primarily due to failure in identifying tetracycline resistance genes in 11 of the 19 MinION assemblies of tetracycline-resistant isolates. The MinION assemblies displayed a large genetic distance from their corresponding MiSeq and hybrid assemblies on the whole-genome phylogenetic tree, indicating that the lower read accuracy of MinION sequencing caused incorrect clustering. The pan-genome of the MinION assemblies contained significantly more accessory genes and less core genes compared to the MiSeq and hybrid assemblies, suggesting that although these assemblies were more contiguous, their sequencing errors reduced accurate genome annotations. Our research demonstrates that MinION sequencing by itself provides an efficient assessment of the genome structure, antimicrobial resistance, and virulence potential of Salmonella; however, it is not sufficient for whole-genome phylogenetic and pan-genome analyses. MinION in combination with MiSeq facilitated the most accurate genomic analyses.

Escherichia coli serotype O157:H7 continues to pose a serious health threat to human beings. Cattle, a major reservoir of the pathogen, harbor E. coli O157:H7 in their gastrointestinal tract and shed variable concentrations of E. coli O157:H7 into the environment. Genetic characterization of cattle-shed E. coli O157 strains is of interest to the livestock industry, food business, and public health community. The present study applied whole genome shotgun sequencing (WGS) and single nucleotide variant (SNV) calling to characterize 279 cattle-shed E. coli O157:H7 strains isolated from a single feedlot located in southwestern region of the US. More than 4,000 SNVs were identified among the strains and the resultant phylogenomic tree revealed three major groups. Using the Sakai strain genome as reference, more than 2,000 SNVs were annotated and a detailed SNV map generated. Results clearly revealed highly polymorphic loci along the E. coli O157:H7 genome that aligned with the prophage regions and highly variant genes involved in processing bacterial genetic information. The WGS data were further profiled against a comprehensive virulence factor database (VFDB) for virulence gene identification. Among the total 285 virulence genes identified, only 132 were present in all the strains. There were six virulence genes unique to single isolates. Our findings suggested that the genome variations of the E. coli O157:H7 were mainly attributable to dynamics of certain phages, and the bacterial strains have variable virulence gene profiles, even though they came from a single cattle population, which may explain the differences in pathogenicity, host prevalence, and transmissibility by E. coli O157:H7.

An important data gap in our understanding of the phyllosphere surrounds the origin of the many microbes described as phyllosphere communities. Most sampling in phyllosphere research has focused on the collection of microbiota without the use of a control, so the opportunity to determine which taxa are actually driven by the biology and physiology of plants as opposed to introduced by environmental forces has yet to be fully realized. To address this data gap, we used plastic plants as inanimate controls adjacent to live tomato plants (phyllosphere) in the field with the hope of distinguishing between bacterial microbiota that may be endemic to plants as opposed to introduced by environmental forces. Using 16S rRNA gene amplicons to study bacterial membership at four time points, we found that the vast majority of all species-level operational taxonomic units were shared at all time-points. Very few taxa were unique to phyllosphere samples. A higher taxonomic diversity was consistently observed in the control samples. The high level of shared taxonomy suggests that environmental forces likely play a very important role in the introduction of microbes to plant surfaces. The observation that very few taxa were unique to the plants compared to the number that were unique to controls was surprising and further suggests that a subset of environmentally introduced taxa thrive on plants. This finding has important implications for improving our approach to the description of core phytobiomes as well as potentially helping us better understand how foodborne pathogens may become associated with plant surfaces.

Establishing an association between possible food sources and clinical isolates requires discriminating the suspected pathogen from an environmental background, and distinguishing it from other closely-related foodborne pathogens. We used whole genome sequencing (WGS) to Salmonella subspecies enterica serotype Tennessee (S. Tennessee) to describe genomic diversity across the serovar as well as among and within outbreak clades of strains associated with contaminated peanut butter. We analyzed 71 isolates of S. Tennessee from disparate food, environmental, and clinical sources and 2 other closely-related Salmonella serovars as outgroups (S. Kentucky and S. Cubana), which were also shot-gun sequenced. A whole genome single nucleotide polymorphism (SNP) analysis was performed using a maximum likelihood approach to infer phylogenetic relationships. Several monophyletic lineages of S. Tennessee with limited SNP variability were identified that recapitulated several food contamination events. S. Tennessee clades were separated from outgroup salmonellae by more than sixteen thousand SNPs. Intra-serovar diversity of S. Tennessee was small compared to the chosen outgroups (1,153 SNPs), suggesting recent divergence of some S. Tennessee clades. Analysis of all 1,153 SNPs structuring an S. Tennessee peanut butter outbreak cluster revealed that isolates from several food, plant, and clinical isolates were very closely related, as they had only a few SNP differences between them. SNP-based cluster analyses linked specific food sources to several clinical S. Tennessee strains isolated in separate contamination events. Environmental and clinical isolates had very similar whole genome sequences; no markers were found that could be used to discriminate between these sources. Finally, we identified SNPs within variable S. Tennessee genes that may be useful markers for the development of rapid surveillance and typing methods, potentially aiding in traceback efforts during future outbreaks. Using WGS can delimit contamination sources for foodborne illnesses across multiple outbreaks and reveal otherwise undetected DNA sequence differences essential to the tracing of bacterial pathogens as they emerge.

Susceptibility of a particular species to prion disease is affected by small differences in the sequence of PrP and correlates with the propensity of its PrP to assume the β-state. A helix-cap motif in the β2-α2-loop of native α-helical rabbit PrP, a resistant species, contains sequence differences that influence intra- and interspecies transmission. To determine the effect the helix-cap motif on β-state refolding propensity, we mutated S170N, S174N, and S170N/S174N of the rabbit PrP helix-cap to resemble that of hamster PrP and conversely, N170S, N174S, and N170S/N174S of hamster PrP to resemble the helix-cap of rabbit PrP. High-resolution crystal structures (1.45-1.6 Å) revealed that these mutations ablate hydrogen-bonding interactions within the helix-cap motif in rabbit PrP(C). They also alter the β-state-misfolding propensity of PrP; the serine mutations in hamster PrP decrease the propensity up to 35%, whereas the asparagine mutations in rabbit PrP increase it up to 42%. Rapid dilution of rabbit and hamster into β-state buffer conditions causes quick conversion to β-state monomers. Kinetic monitoring using size-exclusion chromatography showed that the monomer population decreases exponentially mirrored by an increase in an octameric species. The monomer-octamer transition rates are faster for hamster than for rabbit PrP. The N170S/N174S mutant of hamster PrP has a smaller octamer component at the endpoint compared to the wild-type, whereas the kinetics of octamer formation in mutant and wild-type rabbit PrP are comparable. These findings demonstrate that the sequence of the β2-α2 helix-cap affects refolding to the β-state and subsequently, may influence susceptibility to prion disease.

The ability to detect a specific organism from a complex environment is vitally important to many fields of public health, including food safety. For example, tomatoes have been implicated numerous times as vehicles of foodborne outbreaks due to strains of Salmonella but few studies have ever recovered Salmonella from a tomato phyllosphere environment. Precision of culturing techniques that target agents associated with outbreaks depend on numerous factors. One important factor to better understand is which species co-enrich during enrichment procedures and how microbial dynamics may impede or enhance detection of target pathogens. We used a shotgun sequence approach to describe taxa associated with samples pre-enrichment and throughout the enrichment steps of the Bacteriological Analytical Manual's (BAM) protocol for detection of Salmonella from environmental tomato samples. Recent work has shown that during efforts to enrich Salmonella (Proteobacteria) from tomato field samples, Firmicute genera are also co-enriched and at least one co-enriching Firmicute genus (Paenibacillus sp.) can inhibit and even kills strains of Salmonella. Here we provide a baseline description of microflora that co-culture during detection efforts and the utility of a bioinformatic approach to detect specific taxa from metagenomic sequence data. We observed that uncultured samples clustered together with distinct taxonomic profiles relative to the three cultured treatments (Universal Pre-enrichment broth (UPB), Tetrathionate (TT), and Rappaport-Vassiliadis (RV)). There was little consistency among samples exposed to the same culturing medias, suggesting significant microbial differences in starting matrices or stochasticity associated with enrichment processes. Interestingly, Paenibacillus sp. (Salmonella inhibitor) was significantly enriched from uncultured to cultured (UPB) samples. Also of interest was the sequence based identification of a number of sequences as Salmonella despite indication by all media, that samples were culture negative for Salmonella. Our results substantiate the nascent utility of metagenomic methods to improve both biological and bioinformatic pathogen detection methods.

Salmonella Newport has ranked in the top three Salmonella serotypes associated with foodborne outbreaks from 1995 to 2011 in the United States. In the current study, we selected 26 S. Newport strains isolated from diverse sources and geographic locations and then conducted 454 shotgun pyrosequencing procedures to obtain 16-24 × coverage of high quality draft genomes for each strain. Comparative genomic analysis of 28 S. Newport strains (including 2 reference genomes) and 15 outgroup genomes identified more than 140,000 informative SNPs. A resulting phylogenetic tree consisted of four sublineages and indicated that S. Newport had a clear geographic structure. Strains from Asia were divergent from those from the Americas. Our findings demonstrated that analysis using whole genome sequencing data resulted in a more accurate picture of phylogeny compared to that using single genes or small sets of genes. We selected loci around the mutS gene of S. Newport to differentiate distinct lineages, including those between invH and mutS genes at the 3' end of Salmonella Pathogenicity Island 1 (SPI-1), ste fimbrial operon, and Clustered, Regularly Interspaced, Short Palindromic Repeats (CRISPR) associated-proteins (cas). These genes in the outgroup genomes held high similarity with either S. Newport Lineage II or III at the same loci. S. Newport Lineages II and III have different evolutionary histories in this region and our data demonstrated genetic flow and homologous recombination events around mutS. The findings suggested that S. Newport Lineages II and III diverged early in the serotype evolution and have evolved largely independently. Moreover, we identified genes that could delineate sublineages within the phylogenetic tree and that could be used as potential biomarkers for trace-back investigations during outbreaks. Thus, whole genome sequencing data enabled us to better understand the genetic background of pathogenicity and evolutionary history of S. Newport and also provided additional markers for epidemiological response.