ABSTRACT
Oxford Nanopore long reads of simulated bacterial communities from fresh spinach and surface water were generated (R9.4.1+SQK-LSK109 and R10.4+SQK-LSK112; 0.5, one, and two million reads). Salmonella enterica serotype Heidelberg, Montevideo, or Typhimurium was included alone or in combination in the spinach community, while the water community harbored Pseudomonas aeruginosa.
ABSTRACT
Antimicrobial resistance (AMR) is a world-wide public health threat that is projected to lead to 10 million annual deaths globally by 2050. The AMR public health issue has led to the development of action plans to combat AMR, including improved antimicrobial stewardship, development of new antimicrobials, and advanced monitoring. The National Antimicrobial Resistance Monitoring System (NARMS) led by the United States (U.S) Food and Drug Administration along with the U.S. Centers for Disease Control and U.S. Department of Agriculture has monitored antimicrobial resistant bacteria in retail meats, humans, and food animals since the mid 1990's. NARMS is currently exploring an integrated One Health monitoring model recognizing that human, animal, plant, and environmental systems are linked to public health. Since 2020, the U.S. Environmental Protection Agency has led an interagency NARMS environmental working group (EWG) to implement a surface water AMR monitoring program (SWAM) at watershed and national scales. The NARMS EWG divided the development of the environmental monitoring effort into five areas: (i) defining objectives and questions, (ii) designing study/sampling design, (iii) selecting AMR indicators, (iv) establishing analytical methods, and (v) developing data management/analytics/metadata plans. For each of these areas, the consensus among the scientific community and literature was reviewed and carefully considered prior to the development of this environmental monitoring program. The data produced from the SWAM effort will help develop robust surface water monitoring programs with the goal of assessing public health risks associated with AMR pathogens in surface water (e.g., recreational water exposures), provide a comprehensive picture of how resistant strains are related spatially and temporally within a watershed, and help assess how anthropogenic drivers and intervention strategies impact the transmission of AMR within human, animal, and environmental systems.
ABSTRACT
A total of 55 food and clinical S. Schwarzengrund isolates were assayed for plasmid content, among which an IncFIB-IncFIC(FII) fusion plasmid, conferring streptomycin resistance, was detected in 17 isolates. Among the 17 isolates, 9 were food isolates primarily collected from poultry meat, and 8 clinical isolates collected from stool, urine, and gallbladder. SNP-based phylogenetic analyses showed that the isolates carrying the fusion plasmid formed a subclade indicating the plasmid was acquired and is now maintained by the lineage. Phylogenetic analysis of the plasmid suggested it is derived from avian pathogenic plasmids and might confer an adaptive advantage to the S. Schwarzengrund isolates within birds. IncFIB-IncFIC(FII) fusion plasmids from all food and three clinical isolates were self-conjugative and successfully transferred into E. coli J53 by conjugation. Food and clinical isolates had similar virulome profiles and were able to invade human Caco-2 cells. However, the IncFIB-IncFIC(FII) plasmid did not significantly add to their invasion and persistence potential in human Caco-2 cells.
ABSTRACT
Oxford Nanopore sequencing is one of the high-throughput sequencing technologies that facilitates the reconstruction of metagenome-assembled genomes (MAGs). This study aimed to assess the potential of long-read assembly algorithms in Oxford Nanopore sequencing to enhance the MAG-based identification of bacterial pathogens using both simulated and mock communities. Simulated communities were generated to mimic those on fresh spinach and in surface water. Long reads were produced using R9.4.1+SQK-LSK109 and R10.4 + SQK-LSK112, with 0.5, 1, and 2 million reads. The simulated bacterial communities included multidrug-resistant Salmonella enterica serotypes Heidelberg, Montevideo, and Typhimurium in the fresh spinach community individually or in combination, as well as multidrug-resistant Pseudomonas aeruginosa in the surface water community. Real data sets of the ZymoBIOMICS HMW DNA Standard were also studied. A bioinformatic pipeline (MAGenie, freely available at https://github.com/jackchen129/MAGenie) that combines metagenome assembly, taxonomic classification, and sequence extraction was developed to reconstruct draft MAGs from metagenome assemblies. Five assemblers were evaluated based on a series of genomic analyses. Overall, Flye outperformed the other assemblers, followed by Shasta, Raven, and Unicycler, while Canu performed least effectively. In some instances, the extracted sequences resulted in draft MAGs and provided the locations and structures of antimicrobial resistance genes and mobile genetic elements. Our study showcases the viability of utilizing the extracted sequences for precise phylogenetic inference, as demonstrated by the consistent alignment of phylogenetic topology between the reference genome and the extracted sequences. R9.4.1+SQK-LSK109 was more effective in most cases than R10.4+SQK-LSK112, and greater sequencing depths generally led to more accurate results.IMPORTANCEBy examining diverse bacterial communities, particularly those housing multiple Salmonella enterica serotypes, this study holds significance in uncovering the potential of long-read assembly algorithms to improve metagenome-assembled genome (MAG)-based pathogen identification through Oxford Nanopore sequencing. Our research demonstrates that long-read assembly stands out as a promising avenue for boosting precision in MAG-based pathogen identification, thus advancing the development of more robust surveillance measures. The findings also support ongoing endeavors to fine-tune a bioinformatic pipeline for accurate pathogen identification within complex metagenomic samples.
Subject(s)
Algorithms , Genome, Bacterial , High-Throughput Nucleotide Sequencing , Metagenome , Nanopore Sequencing , Nanopore Sequencing/methods , High-Throughput Nucleotide Sequencing/methods , Bacteria/genetics , Bacteria/classification , Bacteria/isolation & purification , Computational Biology/methods , Salmonella enterica/genetics , Salmonella enterica/classification , Salmonella enterica/isolation & purification , Metagenomics/methods , Pseudomonas aeruginosa/genetics , Pseudomonas aeruginosa/isolation & purification , Pseudomonas aeruginosa/classificationABSTRACT
Surface waters are considered ecological habitats where Salmonella enterica can persist and disseminate to fresh produce production systems. This study aimed to explore the genomic profiles of S. enterica serotypes Typhimurium, Newport, and Infantis from surface waters in Chile, Mexico, and Brazil collected between 2019 and 2022. We analyzed the whole genomes of 106 S. Typhimurium, 161 S. Newport, and 113 S. Infantis isolates. Our phylogenetic analysis exhibited distinct groupings of isolates by their respective countries except for a notable case involving a Chilean S. Newport isolate closely related to two Mexican isolates, showing 4 and 13 single nucleotide polymorphisms of difference, respectively. The patterns of the most frequently detected antimicrobial resistance genes varied across countries and serotypes. A strong correlation existed between integron carriage and genotypic multidrug resistance (MDR) across serotypes in Chile and Mexico (R > 0.90, P < 0.01), while integron(s) were not detected in any of the Brazilian isolates. By contrast, we did not identify any strong correlation between plasmid carriage and genotypic MDR across diverse countries and serotypes.IMPORTANCEUnveiling the genomic landscape of S. enterica in Latin American surface waters is pivotal for ensuring public health. This investigation sheds light on the intricate genomic diversity of S. enterica in surface waters across Chile, Mexico, and Brazil. Our research also addresses critical knowledge gaps, pioneering a comprehensive understanding of surface waters as a reservoir for multidrug-resistant S. enterica. By integrating our understanding of integron carriage as biomarkers into broader MDR control strategies, we can also work toward targeted interventions that mitigate the emergence and dissemination of MDR in S. enterica in surface waters. Given its potential implications for food safety, this study emphasizes the critical need for informed policies and collaborative initiatives to address the risks associated with S. enterica in surface waters.
Subject(s)
Drug Resistance, Multiple, Bacterial , Phylogeny , Salmonella enterica , Salmonella typhimurium , Serogroup , Salmonella enterica/genetics , Salmonella enterica/isolation & purification , Salmonella enterica/classification , Salmonella enterica/drug effects , Brazil , Drug Resistance, Multiple, Bacterial/genetics , Mexico , Salmonella typhimurium/genetics , Salmonella typhimurium/isolation & purification , Salmonella typhimurium/drug effects , Salmonella typhimurium/classification , Integrons/genetics , Genome, Bacterial , Chile , Genomics , Anti-Bacterial Agents/pharmacology , Latin America , Water Microbiology , Polymorphism, Single Nucleotide , Plasmids/genetics , Microbial Sensitivity TestsABSTRACT
Seventeen Salmonella enterica serovar Schwarzengrund isolates from chicken (n = 9) and clinical samples including stool (n = 6), urine (n = 1), and gallbladder (n = 1) were sequenced and found to carry an IncFIB-IncFIC (FII) fusion plasmid of approximately 145 Kb. This information provides reference genomic data for comparative studies of S. Schwarzengrund pathogenicity and plasmid genetics.
ABSTRACT
IMPORTANCE: In developed countries, the human diet is predominated by food commodities, which have been manufactured, processed, and stored in a food production facility. Little is known about the application of metagenomic sequencing approaches for detecting foodborne pathogens, such as L. monocytogenes, and characterizing microbial diversity in food production ecosystems. In this work, we investigated the utility of 16S rRNA amplicon and quasimetagenomic sequencing for the taxonomic and phylogenetic classification of Listeria culture enrichments of environmental swabs collected from dairy and seafood production facilities. We demonstrated that single-nucleotide polymorphism (SNP) analyses of L. monocytogenes metagenome-assembled genomes (MAGs) from quasimetagenomic data sets can achieve similar resolution as culture isolate whole-genome sequencing. To further understand the impact of genome coverage on MAG SNP cluster resolution, an in silico downsampling approach was employed to reduce the percentage of target pathogen sequence reads, providing an initial estimate of required MAG coverage for subtyping resolution of L. monocytogenes.
Subject(s)
Listeria monocytogenes , Humans , Listeria monocytogenes/genetics , Food Microbiology , Phylogeny , RNA, Ribosomal, 16S/genetics , Ecosystem , SeafoodABSTRACT
Most current Salmonella subtyping analyses rely on whole genome sequencing (WGS), which focuses on the high-resolution analysis of single genomes or multiple single genomes from the isolated colonies on microbiological agar plates. In this study, we introduce bioinformatics innovations for a metagenomic outbreak response workflow that accurately identifies multiple Salmonella serovars at the same time. bettercallsal is one of the first analysis tools to identify multiple Salmonella enterica serotypes from metagenomic or quasi-metagenomic datasets with high accuracy, allowing these isolate-independent methods to be incorporated into surveillance and root cause investigations. It was tested on an in silico benchmark dataset comprising 29 unique Salmonella serovars, 46 non-Salmonella bacterial genomes, and 10 viral genomes at varying read depths and on previously well-characterized and sequenced non-selective primary and selective enrichments of papaya and peach samples from separate outbreak investigations that resulted in the identification of multiple Salmonella serovars using traditional isolate culturing and WGS as well as nucleic acid assays. Analyses were also conducted on these datasets using a custom-built k-mer tool, SeqSero2, and Kallisto to compare serotype calling to bettercallsal. The in silico dataset analyzed with bettercallsal achieved the maximum precision, recall, and accuracy of 100, 83, and 94%, respectively. In the papaya outbreak samples, bettercallsal identified the presence of multiple serovars in agreement with the Luminex® xMAP assay results and also identified more serovars per sample, as evidenced by NCBI SNP clustering. In peach outbreak samples, bettercallsal identified two serovars in concordance with k-mer analysis and the Luminex xMAP assay. The genome hit reported by bettercallsal clustered with the chicken isolate genome, as reported by the FDA peach outbreak investigation from sequenced isolates (WGS). Overall, bettercallsal outperformed k-mer, Seqsero2, and Kallisto in identifying multiple serovars from enrichment cultures using shotgun metagenomic sequencing.
ABSTRACT
Aeromonas hydrophila and other closely related Aeromonas species cause motile aeromonad septicemia, a common fish disease. The disease affects many aquaculture sectors potentially requiring antimicrobial treatments. Therefore, researchers and laboratory diagnosticians need criteria called epidemiological cutoff values (ECVs) to determine whether a bacterial isolate has developed decreased susceptibility to an antimicrobial. To generate ECVs for this bacterium, we assembled a diverse collection of 245 isolates previously identified as A. hydrophila from fish. Using rpoD sequencing, we confirmed that 97 of the 245 isolates were A. hydrophila. We allocated the isolates among three laboratories and tested their susceptibility against eight antimicrobials using standard Clinical and Laboratory Standards Institute (CLSI) disk diffusion and broth microdilution methods. The resulting frequency distributions were statistically analyzed to determine wild-type cutoff estimates, which, along with scatterplots, were used to estimate potential ECVs. In collaboration with the CLSI, aquaculture working group, we proposed ECVs for six of the eight antimicrobials tested. Subsequently, the CLSI Subcommittee on Veterinary Antimicrobial Susceptibility Testing reviewed our data and approved the ECVs to be added to the 2020 edition of the VET04 performance standards for antimicrobial susceptibility testing of aquatic bacteria.
Subject(s)
Aeromonas , Anti-Infective Agents , Aeromonas/genetics , Aeromonas hydrophila , Animals , Anti-Bacterial Agents/pharmacology , Fishes , Microbial Sensitivity TestsABSTRACT
ABSTRACT: Cold-smoked salmon is a ready-to-eat seafood product of high commercial importance. The processing and storage steps facilitate the introduction, growth, and persistence of foodborne pathogens and spoilage bacteria. The growth of commensal bacteria during storage and once the product is opened also influence the quality and safety of cold-smoked salmon. Here we investigated the microbial community through targeted 16S rRNA gene and shotgun metagenomic sequencing as means to better understand the interactions among bacteria in cold-smoked salmon. Cold-smoked salmon samples were tested over 30 days of aerobic storage at 4°C and cultured at each time point in a buffered Listeria enrichment broth (BLEB) commonly used to detect Listeria in foods. The microbiomes were composed of Firmicutes and Proteobacteria, namely, Carnobacterium, Brochothrix, Pseudomonas, Serratia, and Psychrobacter. Pseudomonas species were the most diverse species, with 181 taxa identified. In addition, we identified potential homologs to 10 classes of bacteriocins in microbiomes of cold-smoked salmon stored at 4°C and corresponding BLEB culture enrichments. The findings presented here contribute to our understanding of microbiome population dynamics in cold-smoked salmon, including changes in bacterial taxa during aerobic cold storage and after culture enrichment. This may facilitate improvements to pathogen detection and quality preservation of this food.
Subject(s)
Listeria monocytogenes , Microbiota , Animals , Cold Temperature , Colony Count, Microbial , Food Microbiology , Food Preservation , Population Dynamics , RNA, Ribosomal, 16S , Salmon/microbiology , Seafood/microbiology , SmokeABSTRACT
ABSTRACT: Listeria monocytogenes (Lm) is one of the leading causes of death because of foodborne illness, affecting the elderly, pregnant women, neonates, and people who are immunocompromised. Serologically, Lm can be classified into 13 serotypes, although only 4 are typically linked with food contamination and illness. Since 2000, a shift in serotypes involved in listeriosis outbreaks has been observed, suggesting that tracking of serotypes could help identify emerging trends. A PCR method developed in 2004 allowed detection of the four major serotypes as molecular serogroups, corresponding to broad phylogenetic groups. In this study, a novel quantitative PCR (qPCR) method was developed that uses two multiplex qPCRs, one to confirm the Listeria genus and Lm species and the second for Lm molecular serogrouping. This method was compared with the U.S. Food and Drug Administration Bacteriological Analytical Manual (BAM) method for Lm and the seroagglutination method, using a 208-strain panel. Comparison of the genus and species qPCR assay with the BAM methods found an equal or slightly higher accuracy for the qPCR method (>98%), compared with the BAM protocol (>96%), when evaluated against independent characterization data. Molecular serogrouping using the qPCR method (96.6%) was more accurate than the seroagglutination assay (75.6%). The qPCR method identified Lm 4bV strains, which could not be resolved using seroagglutination. The qPCR could not identify lineage III and IV serotype 4b strains but did correctly identify 16 of 18 lineage III and IV strains. The qPCR method performed genus identification for the Listeria species Lm, L. innocua, L. welshimeri, L. ivanovii, and L. seeligeri. In addition, the method performed species identification for Lm and classified Lm into six molecular serogroups: 2A, 2B, 2C, 4B, NT, and 4bV. This method provided a rapid and accurate confirmation of Lm and serogroup determinations; furthermore, it could help identify otherwise unlinked strains by enabling whole genome sequencing analysis based on broad phylogeny, independent of other information.
Subject(s)
Listeria monocytogenes , Listeria , Listeriosis , Aged , Female , Humans , Infant, Newborn , Listeria monocytogenes/genetics , Phylogeny , Pregnancy , Serogroup , SerotypingABSTRACT
Cronobacter species are opportunistic pathogens capable of causing life-threatening infections in humans, with serious complications arising in neonates, infants, immuno-compromised individuals, and elderly adults. The genus is comprised of seven species: Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite a multiplicity of genomic data for the genus, little is known about likely transmission vectors. Using DNA microarray analysis, in parallel with whole genome sequencing, and targeted PCR analyses, the total gene content of two C. malonaticus, three C. turicensis, and 14 C. sakazaki isolated from various filth flies was assessed. Phylogenetic relatedness among these and other strains obtained during surveillance and outbreak investigations were comparatively assessed. Specifically, microarray analysis (MA) demonstrated its utility to cluster strains according to species-specific and sequence type (ST) phylogenetic relatedness, and that the fly strains clustered among strains obtained from clinical, food and environmental sources from United States, Europe, and Southeast Asia. This combinatorial approach was useful in data mining for virulence factor genes, and phage genes and gene clusters. In addition, results of plasmidotyping were in agreement with the species identity for each strain as determined by species-specific PCR assays, MA, and whole genome sequencing. Microarray and BLAST analyses of Cronobacter fly sequence datasets were corroborative and showed that the presence and absence of virulence factors followed species and ST evolutionary lines even though such genes were orthologous. Additionally, zebrafish infectivity studies showed that these pathotypes were as virulent to zebrafish embryos as other clinical strains. In summary, these findings support a striking phylogeny amongst fly, clinical, and surveillance strains isolated during 2010-2015, suggesting that flies are capable vectors for transmission of virulent Cronobacter spp.; they continue to circulate among United States and European populations, environments, and that this "pattern of circulation" has continued over decades.
ABSTRACT
Coccidiosis is one of the most prevalent diseases seen in the poultry industry leading to excessive economic losses. The aim of this study was to investigate the effect of butyric acid glycerol esters (BE) on the ileal and cecal microbiota in birds challenged with Eimeria maxima (EM). Ross 708 male broilers were fed a diet supplemented with 0 (control) or 0.25% BE from day 1. On day 21, half of the birds were infected with 103 EM oocysts. For determing microbiota, ileal and cecal contents and epithelial scrapings were collected at 7 and 10 D postinfection (PI). Alpha diversity of bacterial communities was mostly affected (P < 0.05) by time PI and EM infection. The richness of luminal bacterial populations in the ileum and ceca was affected (P < 0.05) by addition of BE and by time PI × EM × BE interaction, respectively. In the ileal and cecal luminal and mucosal bacterial communities, permutational multivariate analysis of variance (PERMANOVA, unweighted UniFrac) showed significant (P < 0.05) differences because of time PI and interaction between time PI, EM, and BE. Significant (P < 0.05) differences in taxonomic composition at the family level were observed in microbiota of luminal and mucosal populations of the ileum and ceca owing to time PI, EM, BE, and their interactions. The bacterial community present in the cecal lumen was characterized by the lowest number of differential bacteria, whereas the cecal mucosal community was characterized by the highest number of differentially abundant bacteria. In conclusion, our results show that EM infection and time PI has the biggest impact on microbial diversity in the chicken gut. The presence of BE in the diet had a limited effect on gut microbiota.
Subject(s)
Butyric Acid , Coccidiosis , Eimeria , Esters , Gastrointestinal Microbiome , Poultry Diseases , Animal Feed/analysis , Animals , Butyric Acid/pharmacology , Cecum/microbiology , Chickens , Coccidiosis/microbiology , Coccidiosis/veterinary , Diet/veterinary , Esters/pharmacology , Gastrointestinal Microbiome/drug effects , Glycerol/pharmacology , Ileum/microbiology , Intestinal Mucosa/microbiology , Intestinal Mucosa/parasitology , Male , Poultry Diseases/drug therapyABSTRACT
Salmonella enterica cause significant illnesses worldwide. There has been a marked increase in resistance to fluoroquinolones and ß-lactams/cephalosporins, antibiotics commonly used to treat salmonellosis. However, S. enterica serovars vary in their resistance to these and other antibiotics. The systemic virulence of some Salmonella serovars is due to a low copy number, IncF plasmid (65-100 kb) that contains the ADP-ribosylating toxin, SpvB. This virulence plasmid is present in only nine Salmonella serovars. It is possible that the spvB-virulence plasmid excludes other plasmids and may explain why antibiotic resistance is slow to develop in certain Salmonella serovars such as S. Enteritidis. The distribution of plasmid entry exclusion genes traS/traT and traY/excA are variable in Salmonella IncF and IncI plasmids, respectively and may account for differences in emergent antimicrobial resistance for some Salmonella serovars. The goal of this study is to determine the contribution of the Salmonella spvB-virulence plasmid in F-plasmid exclusion. From conjugation experiments, S. Typhimurium exhibited lower conjugation frequency with incFI and incFII plasmids when the spvB-virulence plasmid is present. Furthermore, introduction of cloned incFI traS into a "plasmidless" S. Typhimurium LT2 strain and Escherichia coli DH5α excluded incFI plasmid. However, deletion of the virulence plasmid traS did not affect plasmid exclusion significantly compared to a spvB control deletion. In addition, differences in F plasmid conjugation in natural Salmonella isolates did not correlate with IncF or SpvB-virulence plasmid genotype. There appear to be other plasmid or chromosomal genes at play in plasmid exclusion that may be responsible for the slow development of antibiotic resistance in certain serovars.
ABSTRACT
BACKGROUND: Salmonella enterica possess several iron acquisition systems, encoded on the chromosome and plasmids. Recently, we demonstrated that incompatibility group (Inc) FIB plasmid-encoded iron acquisition systems (Sit and aerobactin) likely play an important role in persistence of Salmonella in human intestinal epithelial cells (Caco-2). In this study, we sought to determine global transcriptome analyses of S. enterica in iron-rich (IR) and iron-depleted (ID) growth conditions. RESULTS: The number of differentially-expressed genes were substantially higher for recipient (SE819) (n = 966) and transconjugant (TC) (n = 945) compared to the wild type (WT) (SE163A) (n = 110) strain in ID as compared to IR growth conditions. Several virulence-associated factors including T3SS, flagellin, cold-shock protein (cspE), and regulatory genes were upregulated in TC in ID compared to IR conditions. Whereas, IS1 and acrR/tetR transposases located on the IncFIB plasmid, ferritin and several regulatory genes were downregulated in TC in ID conditions. Enterobactin transporter (entS), iron ABC transporter (fepCD), colicin transporter, IncFIB-encoded enolase, cyclic di-GMP regulator (cdgR) and other regulatory genes of the WT strain were upregulated in ID compared to IR conditions. Conversely, ferritin, ferrous iron transport protein A (feoA), IncFIB-encoded IS1 and acrR/tetR transposases and ArtA toxin of WT were downregulated in ID conditions. SDS-PAGE coupled with LC-MS/MS analyses revealed that siderophore receptor proteins such as chromosomally-encoded IroN and, IncFIB-encoded IutA were upregulated in WT and TC in ID growth conditions. Both chromosome and IncFIB plasmid-encoded SitA was overexpressed in WT, but not in TC or recipient in ID conditions. Increased expression of flagellin was detected in recipient and TC, but not in WT in ID conditions. CONCLUSION: Iron concentrations in growth media influenced differential gene expressions both at transcriptional and translational levels, including genes encoded on the IncFIB plasmid. Limited iron availability within the host may promote pathogenic Salmonella to differentially express subsets of genes encoded by chromosome and/or plasmids, facilitating establishment of successful infection.
Subject(s)
Culture Media/chemistry , Gene Expression Profiling , Iron/analysis , Iron/pharmacology , Salmonella enterica/growth & development , Salmonella enterica/genetics , Caco-2 Cells , Humans , Proteomics , Salmonella enterica/drug effectsABSTRACT
Food microbiome composition impacts food safety and quality. The resident microbiota of many food products is influenced throughout the farm to fork continuum by farming practices, environmental factors, and food manufacturing and processing procedures. Currently, most food microbiology studies rely on culture-dependent methods to identify bacteria. However, advances in high-throughput DNA sequencing technologies have enabled the use of targeted 16S rRNA gene sequencing to profile complex microbial communities including non-culturable members. In this study we used 16S rRNA gene sequencing to assess the microbiome profiles of plant and animal derived foods collected at two points in the manufacturing process; post-harvest/pre-retail (cilantro) and retail (cilantro, masala spice mixes, cucumbers, mung bean sprouts, and smoked salmon). Our findings revealed microbiome profiles, unique to each food, that were influenced by the moisture content (dry spices, fresh produce), packaging methods, such as modified atmospheric packaging (mung bean sprouts and smoked salmon), and manufacturing stage (cilantro prior to retail and at retail). The masala spice mixes and cucumbers were comprised mainly of Proteobacteria, Firmicutes, and Actinobacteria. Cilantro microbiome profiles consisted mainly of Proteobacteria, followed by Bacteroidetes, and low levels of Firmicutes and Actinobacteria. The two brands of mung bean sprouts and the three smoked salmon samples differed from one another in their microbiome composition, each predominated by either by Firmicutes or Proteobacteria. These data demonstrate diverse and highly variable resident microbial communities across food products, which is informative in the context of food safety, and spoilage where indigenous bacteria could hamper pathogen detection, and limit shelf life.
ABSTRACT
Recently, we have identified a link between four listeriosis incidents/outbreaks to a variant of Listeria monocytogenes (Lm) serotype 4b strains, 4bV. Although 4bV strains have been reported from clinical specimens as well as from foods, listeriosis outbreaks occurring in 2014-2016 were the first reported outbreaks involving 4bV in the USA. Since traditional typing methods do not detect members of this group, we undertook a systematic and retrospective analysis of all Lm in the NCBI WGS Sequence Read Archive database to investigate the burden of 4bV strains among all listeriosis cases. This analysis identified the presence of isolates causing sporadic cases as well as those associated with the aforementioned outbreaks, as determined by WGS and traditional epidemiology. In total, approximately 350 Lm 4bV strains were identified from multiple parts of the USA as well as from Australia and Chile, dating back to 2001. The genomic relatedness of these strains was compared using the CFSAN SNP Pipeline and multi-virulence-locus sequence typing (MVLST). Using the CFSAN Pipeline tool, the 4bV strains were found to group into seven clusters that were separate from 4b strains. All seven clades appeared to contain isolates from both clinical and non-clinical sources. Conversely, the MVLST analysis revealed that practically all of the strains belonged to a single clade, suggesting that 4bV strains from disparate geographic regions and sources are under varied selective pressure, restricting diversity across these six virulence loci while allowing more variability across the genome as a whole. Further evaluation of these 4bV strains identified genes potentially acquired from a lineage II source external to the lmo0733-lmo0739 region, as well as highly conserved SNPs unique to the 4bV strains when compared to those from other lineages. Taken together, these data suggest that 4bV strains have undergone adaptive responses to selective pressures that may enhance survival in the environment while maintaining the pathogenic potential of serotype 4b strains.
ABSTRACT
16S rRNA community profiling continues to be a useful tool to study microbiome composition and dynamics, in part due to advances in next generation sequencing technology that translate into reductions in cost. Reliable taxonomic identification to the species-level, however, remains difficult, especially for short-read sequencing platforms, due to incomplete coverage of the 16S rRNA gene. This is especially true for Salmonella enterica, which is often found as a low abundant member of the microbial community, and is often found in combination with several other closely related enteric species. Here, we report on the evaluation and application of Resphera Insight, an ultra-high resolution taxonomic assignment algorithm for 16S rRNA sequences to the species level. The analytical pipeline achieved 99.7% sensitivity to correctly identify S. enterica from WGS datasets extracted from the FDA GenomeTrakr Bioproject, while demonstrating 99.9% specificity over other Enterobacteriaceae members. From low-diversity and low-complexity samples, namely ice cream, the algorithm achieved 100% specificity and sensitivity for Salmonella detection. As demonstrated using cilantro and chili powder, for highly complex and diverse samples, especially those that contain closely related species, the detection threshold will likely have to be adjusted higher to account for misidentifications. We also demonstrate the utility of this approach to detect Salmonella in the clinical setting, in this case, bloodborne infections.
ABSTRACT
Salmonella Typhimurium is the leading cause of human non-typhoidal gastroenteritis in the US. S. Kentucky is one the most commonly recovered serovars from commercially processed poultry carcasses. This study compared the genotypic and phenotypic properties of two Salmonella enterica strains Typhimurium (ST221_31B) and Kentucky (SK222_32B) recovered from commercially processed chicken carcasses using whole genome sequencing, phenotype characterizations and an intracellular killing assay. Illumina MiSeq platform was used for sequencing of two Salmonella genomes. Phylogenetic analysis employing homologous alignment of a 1,185 non-duplicated protein-coding gene in the Salmonella core genome demonstrated fully resolved bifurcating patterns with varying levels of diversity that separated ST221_31B and SK222_32B genomes into distinct monophyletic serovar clades. Single nucleotide polymorphism (SNP) analysis identified 2,432 (ST19) SNPs within 13 Typhimurium genomes including ST221_31B representing Sequence Type ST19 and 650 (ST152) SNPs were detected within 13 Kentucky genomes including SK222_32B representing Sequence Type ST152. In addition to serovar-specific conserved coding sequences, the genomes of ST221_31B and SK222_32B harbor several genomic regions with significant genetic differences. These included phage and phage-like elements, carbon utilization or transport operons, fimbriae operons, putative membrane associated protein-encoding genes, antibiotic resistance genes, siderophore operons, and numerous hypothetical protein-encoding genes. Phenotype microarray results demonstrated that ST221_31B is capable of utilizing certain carbon compounds more efficiently as compared to SK222_3B; namely, 1,2-propanediol, M-inositol, L-threonine, α-D-lactose, D-tagatose, adonitol, formic acid, acetoacetic acid, and L-tartaric acid. ST221_31B survived for 48 h in macrophages, while SK222_32B was mostly eliminated. Further, a 3-fold growth of ST221_31B was observed at 24 hours post-infection in chicken granulosa cells while SK222_32B was unable to replicate in these cells. These results suggest that Salmonella Typhimurium can survive host defenses better and could be more invasive than Salmonella Kentucky and provide some insights into the genomic determinants responsible for these differences.
Subject(s)
Chickens/microbiology , Drug Resistance, Multiple, Bacterial/genetics , Genotype , Phenotype , Salmonella/drug effects , Animals , Phylogeny , Salmonella/geneticsABSTRACT
Four listeriosis incidences/outbreaks, spanning 19 months, have been linked to Listeria monocytogenes serotype 4b variant (4bV) strains. Three of these incidents can be linked to a defined geographical region, while the fourth is likely to be linked. In this study, whole genome sequencing (WGS) of strains from these incidents was used for genomic comparisons using two approached. The first was JSpecies tetramer, which analyzed tetranucleotide frequency to assess relatedness. The second, the CFSAN SNP Pipeline, was used to perform WGS SNP analyses against three different reference genomes to evaluate relatedness by SNP distances. In each case, unrelated strains were included as controls. The analyses showed that strains from these incidents form a highly related clade with SNP differences of ≤101 within the clade and >9000 against other strains. Multi-Virulence-Locus Sequence Typing, a third standardized approach for evaluation relatedness, was used to assess the genetic drift in six conserved, known virulence loci and showed a different clustering pattern indicating possible differences in selection pressure experienced by these genes. These data suggest a high degree of relatedness among these 4bV strains linked to a defined geographic region and also highlight the possibility of alterations related to adaptation and virulence.