Search

Volatile organic compounds (VOCs) are crucial precursors in the formation of ozone (O3). The sources of pollution are complex and significantly impact O3 generation. Long-term exposure to high-concentration O3 environments causes serious damage to organisms. High-altitude areas experience continuous high temperatures and strong solar radiation, which can easily produce O3 through photochemical reactions, making these areas prone to frequent air pollution. This study utilizes the National Positioning Station of the Yinchuan Urban Ecosystem in Ningxia to conduct field synchronous observations, gathering data on VOCs, meteorological variables, and O3. By employing machine learning algorithms, we examined the seasonal distribution characteristics of VOCs, their pollution sources, and their impact on O3. The results show that the seasonal distribution of VOCs areas is low in spring and autumn and high in summer and winter. The components with higher volume fraction contribution are m-tolualdehyde in spring and summer, ethane in autumn, and acetylene in winter. The main emission sources of VOCs pollution are hydrocarbon volatile emission in spring and winter, solvent volatilization emission in summer, and industrial sources in autumn. VOCs have a negative effect on O3, with a greater impact in winter and spring, evidenced by standardized effect values of −0.26 and −0.24. The key components affecting O3 in VOCs vary by season. Aromatic hydrocarbons are the key components in spring and winter, with contribution rates of 22 % and 21.3 %. Alkenes are the key components in summer and autumn (24.5 % and 26.8 %). Among meteorological variables, temperature is the key factor affecting O3, while wind speed is the key factor affecting O3 only in winter. This study aims to clarify the sources of VOCs pollution in different seasons and their impact on O3, providing a theoretical basis and technical support for the prevention and control of VOCs and O3 pollution in high-altitude areas.

ABSTRACT This study examined the diversity and persistence of Salmonella in the surface waters of agricultural regions of Brazil, Chile, and Mexico. Research groups (three in 2019–2020 and five in 2021–2022) conducted a long-term survey of surface water across 5–8 months annually (n = 30 monthly). On-site, each team filtered 10-L water samples with modified Moore Swabs to capture Salmonella, which were then isolated and identified using conventional microbiological techniques. Salmonella isolates were sequenced on Illumina platforms. Salmonella was present in 1,493/3,291 water samples (45.8%), with varying isolation rates across countries and years. Newport, Infantis, and Typhimurium were the most frequent among the 128 different serovars. Notably, 22 serovars were found in all three countries, representing almost half of the 1,911 different isolates collected. The resistome comprised 72 antimicrobial resistance (AMR) genes and six point mutations in three genes. At least one AMR determinant was observed in 33.8% (646/1,911) of the isolates, of which 47.4% (306/646) were potentially multidrug resistant. Phylogeny based on core genome multilocus sequence typing (cgMLST) showed that most isolates clustered according to sequence type and country of origin. Only 14 cgMLST multi-country clusters were detected among the 275 clusters. However, further analysis confirmed that close genetic relatedness occurred mostly among isolates from the same country, with three exceptions. Interestingly, isolates closely related phylogenetically were recovered over multiple years within the same country, indicating the persistence of certain Salmonella in those areas. In conclusion, surface waters in these regions are consistently contaminated with diverse Salmonella, including strains that persist over time.IMPORTANCESalmonella is a leading foodborne pathogen responsible for millions of illnesses, hospitalizations, and deaths annually. Although Salmonella-contaminated water has now been recognized as an important contamination source in the agrifood chain, there is a lack of knowledge on the global occurrence and diversity of Salmonella in surface water. Moreover, there has been insufficient research on Salmonella in surface waters from Latin American countries that are major producers and exporters of agricultural products. Incorporating genetic profiling of Salmonella isolates from underrepresented regions, such as Latin America, enhances our understanding of the pathogen’s ecology, evolution, antimicrobial resistance, and pathogenicity. Moreover, leveraging genomic data derived from pathogens isolated from diverse geographical areas is critical for assessing the potential public health risk posed by the pathogen and expediting investigations of foodborne outbreaks. Ultimately, global efforts contribute significantly to reducing the incidence of foodborne infections.

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.

ABSTRACTSurface 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.

Summary: Background: Multidrug-resistant (MDR) Salmonella Infantis has disseminated worldwide, mainly linked to the consumption of poultry products. Evidence shows dissemination of this pathogen in Chile; however, studies are primarily limited to phenotypic data or involve few isolates. As human cases of Salmonella Infantis infections have substantially increased in recent years, this study aimed to characterise the genomic epidemiology and antimicrobial-resistance profiles of isolates obtained from different sources, aiming to inform effective surveillance and control measures. Methods: We sequenced 396 Salmonella Infantis genomes and analysed them with all publicly available genomes of this pathogen from Chile (440 genomes in total), representing isolates from environmental, food, animal, and human sources obtained from 2009 to 2022. Based on bioinformatic and phenotypic methods, we assessed the population structure, dissemination among different niches, and antimicrobial resistance (AMR) profiles of Salmonella Infantis in the country. Findings: The genomic and phylogenetic analyses showed that Salmonella Infantis from Chile comprised several clusters of highly related isolates dominated by sequence type 32. The HC20_343 cluster grouped an important proportion of all isolates. This was the only cluster associated with pESI-like megaplasmids, and up to 12 acquired AMR genes/mutations predicted to result in an MDR phenotype. Accordingly, antimicrobial-susceptibility testing revealed a strong concordance between the AMR genetic determinants and their matching phenotypic expression, indicating that a significant proportion of HC20_343 isolates produce extended-spectrum β-lactamases and have intermediate fluoroquinolone resistance. HC20_343 Salmonella Infantis were spread among environmental, animal, food, and human niches, showing a close relationship between isolates from different years and sources, and a low intra-source genomic diversity. Interpretation: Our findings show a widespread dissemination of MDR Salmonella Infantis from the HC20_343 cluster in Chile. The high proportion of isolates with resistance to first-line antibiotics and the evidence of active transmission between the environment, animals, food, and humans highlight the urgency of improved surveillance and control measures in the country. As HC20_343 isolates predominate in the Americas, our results suggest a high prevalence of ESBL-producing Salmonella Infantis with intermediate fluoroquinolone resistance in the continent. Funding: Partially supported by the Food and Drug Administration (FDA) of the U.S. Department of Health and Human Services as part of an award, FDU001818, with 30% percent funded by FDA/HHS; and by Agencia de Investigación y Desarrollo de Chile (ANID) through FONDECYT de Postdoctorado Folio 3230796 and Folio 3210317, FONDECYT Regular Folio 1231082, and ANID—Millennium Science Initiative Program—ICN2021_044.

Abstract Multidrug-resistant (MDR) Salmonella has been a long-standing challenge in public health and food safety. The prevalence of MDR S. Enteritidis, especially isolated from humans, in China is significantly higher than those from the U.S. and other countries. A dataset of 197 S. Enteritidis genomes, including 16 sequenced clinical isolates from China and 181 downloaded genomes of human isolates from the U.S., Europe, and Africa, was analyzed for genomic diversity, virulence potential, and antimicrobial resistance (AMR). Phylogenomic analyses identified four major well-supported clades (I–IV). While AMR genotype in the majority of isolates in clades I and IV displayed as pan-susceptible, 81.8% (9/11) and 22.4% (13/58) of isolates in clades III and II were MDR, respectively. It is noted that 77% (10/13) of MDR isolates in clade II were from China. The most common antimicrobial resistance genes (ARGs) carried by the Chinese isolates were aph(3′)-IIa, bla CTX-M-55, and bla TEM-1B, whereas bla TEM-1B, sul1, sul2, drfA7, aph(3")-Ib/strA, and aph(6)-Id/strB were most often identified in those from Africa (clade III). Among the 14 plasmid types identified, IncX1 and IncFII(pHN7A8) were found exclusively in the Chinese MDR isolates, while IncQ1 was highly associated with the African MDR isolates. The spvRABCD virulence operon was present in 94.9% (187/197) of isolates tested and was highly associated with both the IncF (IncFII and IncFIB) plasmids. In addition, phylogenetic differences in distribution of Salmonella pathogenicity islands (SPIs), prophages and other accessory genes were also noted. Taken together, these findings provide new insights into the molecular mechanisms underpinning diversification of MDR S. Enteritidis.

ABSTRACT Phylogeny is a powerful tool that can be incorporated into quantitative descriptions of community diversity, yet its use has been limited largely due to the difficulty in constructing phylogenies which incorporate the wide genomic diversity of microbial communities. Here, we describe the development of a web portal, PhyloPlus, which enables users to generate customized phylogenies that may be applied to any bacterial or archaeal communities. We demonstrate the power of phylogeny by comparing metrics that employ phylogeny with those that do not when applied to data sets from two metagenomic studies (fermented food, n = 58; human microbiome, n = 60). This example shows how inclusion of all bacterial species identified by taxonomic classifiers (Kraken2 and Kaiju) made the phylogeny perfectly congruent to the corresponding classification outputs. Our phylogeny-based approach also enabled the construction of more constrained null models which (i) shed light into community structure and (ii) minimize potential inflation of type I errors. Construction of such null models allowed for the observation of under-dispersion in 44 (75.86%) food samples, with the metacommunity defined as bacteria that were found in different food matrices. We also observed that closely related species with high abundance and uneven distribution across different sites could potentially exaggerate the dissimilarity between phylogenetically similar communities if they were measured using traditional species-based metrics (Padj. = 0.003), whereas this effect was mitigated by incorporating phylogeny (Padj. = 1). In summary, our tool can provide additional insights into microbial communities of interest and facilitate the use of phylogeny-based approaches in metagenomic analyses. IMPORTANCE There has been an explosion of interest in how microbial diversity affects human health, food safety, and environmental functions among many other processes. Accurately measuring the diversity and structure of those communities is central to understanding their effects. Here, we describe the development of a freely available online tool, PhyloPlus, which allows users to generate custom phylogenies that may be applied to any data set, thereby removing a major obstacle to the application of phylogeny to metagenomic data analysis. We demonstrate that the genetic relatedness of the organisms within those communities is a critical feature of their overall diversity, and that using a phylogeny which captures and quantifies this diversity allows for much more accurate descriptions while preventing misleading conclusions based on estimates that ignore evolutionary relationships.

The increasing number of studies reporting the presence of Salmonella in environmental water sources suggests that it is beyond incidental findings originated from sparse fecal contamination events. However, there is no consensus on the occurrence of Salmonella as its relative serovar representation across non-recycled water sources. We conducted a meta-analysis of proportions by fitting a random-effects model using the restricted maximum-likelihood estimator to obtain the weighted average proportion and between-study variance associated with the occurrence of Salmonella in water sources. Moreover, meta-regression and non-parametric supervised machine learning method were performed to predict the effect of moderators on the frequency of Salmonella in non-recycled water sources. Three sequential steps (identification of information sources, screening and eligibility) were performed to obtain a preliminary selection from identified abstracts and article titles. Questions related to the frequency of Salmonella in aquatic environments, as well as putative differences in the relative frequencies of the reported Salmonella serovars and the role of potential variable moderators (sample source, country, and sample volume) were formulated according to the population, intervention, comparison, and outcome method (PICO). The results were reported according to the Preferred Reporting Items for Systematic Review and Meta-Analyzes statement (PRISMA). A total of 26 eligible papers reporting 148 different Salmonella serovars were retrieved. According to our model, the Salmonella frequency in non-recycled water sources was 0.19 [CI: 0.14; 0.25]. The source of water was identified as the most import variable affecting the frequency of Salmonella, estimated as 0.31 and 0.17% for surface and groundwater, respectively. There was a higher frequency of Salmonella in countries with lower human development index (HDI). Small volume samples of surface water resulted in lower detectable Salmonella frequencies both in high and low HDI regions. Relative frequencies of the 148 serovars were significantly affected only by HDI and volume. Considering that serovars representation can also be affected by water sample volume, efforts toward the standardization of water samplings for monitoring purposes should be considered. Further approaches such as metagenomics could provide more comprehensive insights about the microbial ecology of fresh water and its importance for the quality and safety of agricultural products.

Abstract Background We benchmarked the hybrid assembly approaches of MaSuRCA, SPAdes, and Unicycler for bacterial pathogens using Illumina and Oxford Nanopore sequencing by determining genome completeness and accuracy, antimicrobial resistance (AMR), virulence potential, multilocus sequence typing (MLST), phylogeny, and pan genome. Ten bacterial species (10 strains) were tested for simulated reads of both mediocre- and low-quality, whereas 11 bacterial species (12 strains) were tested for real reads. Results Unicycler performed the best for achieving contiguous genomes, closely followed by MaSuRCA, while all SPAdes assemblies were incomplete. MaSuRCA was less tolerant of low-quality long reads than SPAdes and Unicycler. The hybrid assemblies of five antimicrobial-resistant strains with simulated reads provided consistent AMR genotypes with the reference genomes. The MaSuRCA assembly of Staphylococcus aureus with real reads contained msr(A) and tet(K), while the reference genome and SPAdes and Unicycler assemblies harbored blaZ. The AMR genotypes of the reference genomes and hybrid assemblies were consistent for the other five antimicrobial-resistant strains with real reads. The numbers of virulence genes in all hybrid assemblies were similar to those of the reference genomes, irrespective of simulated or real reads. Only one exception existed that the reference genome and hybrid assemblies of Pseudomonas aeruginosa with mediocre-quality long reads carried 241 virulence genes, whereas 184 virulence genes were identified in the hybrid assemblies of low-quality long reads. The MaSuRCA assemblies of Escherichia coli O157:H7 and Salmonella Typhimurium with mediocre-quality long reads contained 126 and 118 virulence genes, respectively, while 110 and 107 virulence genes were detected in their MaSuRCA assemblies of low-quality long reads, respectively. All approaches performed well in our MLST and phylogenetic analyses. The pan genomes of the hybrid assemblies of S. Typhimurium with mediocre-quality long reads were similar to that of the reference genome, while SPAdes and Unicycler were more tolerant of low-quality long reads than MaSuRCA for the pan-genome analysis. All approaches functioned well in the pan-genome analysis of Campylobacter jejuni with real reads. Conclusions Our research demonstrates the hybrid assembly pipeline of Unicycler as a superior approach for genomic analyses of bacterial pathogens using Illumina and Oxford Nanopore sequencing.

Probiotics are recognized for outcompeting pathogenic bacteria by competitive receptor-mediated colonization and secretion of functional metabolites which are antimicrobial against certain microbes as well as improving host’s gut health and immunity. Recently, we have constructed a bioactive Lactobacillus casei (LC) strain, LC+mcra, by inserting mcra (myosin cross-reactive antigen) gene, which stimulates the conversion of conjugated linoleic acids. In this study, we evaluated the modulation of gut microbiome and protective roles of LC+mcra against pathogenic Salmonella enterica serovar Typhimurium (ST) and enterohemorrhagic E. coli (EHEC) infections in BALB/cJ mice. We observed that LC+mcra colonized efficiently in mice gut intestine and competitively reduced the infection with ST and EHEC in various locations of small and large intestine, specifically cecum, jejunum, and ileum (p

Metagenomics offers the highest level of strain discrimination of bacterial pathogens from complex food and water microbiota. With the rapid evolvement of assembly algorithms, defining an optimal assembler based on the performance in the metagenomic identification of foodborne and waterborne pathogens is warranted. We aimed to benchmark short-read assemblers for the metagenomic identification of foodborne and waterborne pathogens using simulated bacterial communities. Bacterial communities on fresh spinach and in surface water were simulated by generating paired-end short reads of Illumina HiSeq, MiSeq, and NovaSeq at different sequencing depths. Multidrug-resistant Salmonella Indiana SI43 and Pseudomonas aeruginosa PAO1 were included in the simulated communities on fresh spinach and in surface water, respectively. ABySS, IDBA-UD, MaSuRCA, MEGAHIT, metaSPAdes, and Ray Meta were benchmarked in terms of assembly quality, identifications of plasmids, virulence genes, Salmonella pathogenicity island, antimicrobial resistance genes, chromosomal point mutations, serotyping, multilocus sequence typing, and whole-genome phylogeny. Overall, MEGHIT, metaSPAdes, and Ray Meta were more effective for metagenomic identification. We did not obtain an optimal assembler when using the extracted reads classified as Salmonella or P. aeruginosa for downstream genomic analyses, but the extracted reads showed consistent phylogenetic topology with the reference genome when they were aligned with Salmonella or P. aeruginosa strains. In most cases, HiSeq, MiSeq, and NovaSeq were comparable at the same sequencing depth, while higher sequencing depths generally led to more accurate results. As assembly algorithms advance and mature, the evaluation of assemblers should be a continuous process.

Abstract Background Internalins are surface proteins that are utilized by Listeria monocytogenes to facilitate its invasion into human intestinal epithelial cells. The expression of a full-length InlA is one of essential virulence factors for L. monocytogenes to cross the intestinal barrier in order to invade epithelial cells. Results In this study, the gene sequences of inlA in 120 L. monocytogenes isolates from food (n = 107) and humans (n = 13) were analyzed. Premature stop codon (PMSC) mutations in inlA were identified in 51 isolates (50 from food and 1 from human). Six mutation types of PMSCs were identified. Among the 51 isolates with PMSCs in inlA, there were 44 serogroup 1/2c, 3c isolates from food, of which seven belonged to serogroups 1/2a, 3a. A total of 153,382 SNPs in 2247 core genes from 42 genomes were identified and used to construct a phylogenetic tree. Serotype 1/2c isolates with inlA PMSC mutations were grouped together. Cell culture studies on 21 isolates showed that the invasion to Caco-2 cells was significantly reduced among isolates with inlA PMSC mutations compared to those without PMSC mutations (P

Shiga toxin-producing Escherichia coli (STEC) causes foodborne outbreaks that can lead to complications such as hemolytic uremic syndrome. Their main reservoir is cattle, and ground beef has been frequently associated with disease and outbreaks. In this study, we attempted to understand the genetic relationship among STEC isolated in Chile from different sources, their relationship to STEC from the rest of the world, and to identify molecular markers of Chilean STEC. We sequenced 62 STEC isolated in Chile using MiSeq Illumina. In silico typing was determined using tools of the Center Genomic Epidemiology, Denmark University (CGE/DTU). Genomes of our local STEC collection were compared with 113 STEC isolated worldwide through a core genome MLST (cgMLST) approach, and we also searched for distinct genes to be used as molecular markers of Chilean isolates. Genomes in our local collection were grouped based on serogroup and sequence type, and clusters were formed within local STEC. In the worldwide STEC analysis, Chilean STEC did not cluster with genomes of the rest of the world suggesting that they are not phylogenetically related to previously described STEC. The pangenome of our STEC collection was 11,650 genes, but we did not identify distinct molecular markers of local STEC. Our results showed that there may be local emerging STEC with unique features, nevertheless, no molecular markers were detected. Therefore, there might be elements such as a syntenic organization that might explain differential clustering detected between local and worldwide STEC.

Listeria monocytogenes is a major foodborne pathogen that can contaminate food products and colonize food-producing facilities. Foodservice operations (FSOp) are frequently responsible for foodborne outbreaks due to food safety practices failures. We investigated the presence of and characterized L. monocytogenes from two FSOp (cafeterias) distributing ready-to-eat meals and verified FSOp’s compliance with good manufacturing practices (GMP). Two facilities (FSOp-A and FSOp-B) were visited three times each over 5 months. We sampled foods, ingredients, and surfaces for microbiological analysis, and L. monocytogenes isolates were characterized by phylogenetic analyses and phenotypic characteristics. GMP audits were performed in the first and third visits. A ready-to-eat salad (FSOp-A) and a frozen ingredient (FSOp-B) were contaminated with L. monocytogenes, which was also detected on Zone 3 surfaces (floor, drains, and a boot cover). The phylogenetic analysis demonstrated that FSOp-B had persistent L. monocytogenes strains, but environmental isolates were not closely related to food or ingredient isolates. GMP audits showed that both operations worked under “fair” conditions, and “facilities and equipment” was the section with the least compliances. The presence of L. monocytogenes in the environment and GMP failures could promote food contamination with this pathogen, presenting a risk to consumers.

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.

Water is vital to agriculture. It is essential that the water used for the production of fresh produce commodities be safe. Microbial pathogens are able to survive for extended periods of time in water. It is critical to understand their biology and ecology in this ecosystem in order to develop better mitigation strategies for farmers who grow these food crops. In this review the prevalence, persistence and ecology of four major foodborne pathogens, Shiga toxin-producing Escherichia coli (STEC), Salmonella, Campylobacter and closely related Arcobacter, and Listeria monocytogenes, in water are discussed. These pathogens have been linked to fresh produce outbreaks, some with devastating consequences, where, in a few cases, the contamination event has been traced to water used for crop production or post-harvest activities. In addition, antimicrobial resistance, methods improvements, including the role of genomics in aiding in the understanding of these pathogens, are discussed. Finally, global initiatives to improve our knowledge base of these pathogens around the world are touched upon.

The exposure of foodborne pathogens such as Salmonella enterica to a sub-lethal stress may protect bacterial cells against distinct stresses during the production of leafy greens, which can constitute potential health hazards to consumers. In this study, we evaluated how the prior exposure of S. enterica to sub-lethal food processing-related stresses influenced its subsequent persistence on baby spinach under cold (4 °C for 7 days) and temperature abuse (37 °C for 2 h + 4 °C for 7 days) conditions. We also compared the survival characteristics of pre-stressed S. enterica and Enterococcus faecium NRRL B-2354 as its surrogate on baby spinach. A cocktail of three S. enterica serovars, as well as S. Typhimurium ATCC 14028 wild type and its ΔrpoS mutant, and E. faecium NRRL B-2354, was first exposed to sub-lethal desiccation, oxidation, heat shock, and acid stresses. Afterward, baby spinach was inoculated with unstressed or pre-stressed cells at 7.0 log CFU/sample unit, followed by 7-day storage under cold and temperature abuse conditions. The unstressed S. enterica (fresh cells in sterile 0.85% saline) decreased rapidly within the first day and thereafter persisted around 5.5 log CFU/sample unit under both conditions. The desiccation-stressed S. enterica showed the highest bacterial counts (p < 0.05) compared to other conditions. The unstressed S. enterica survived better (p < 0.05) than the oxidation- and acid-stressed S. enterica, while there were no significant differences (p > 0.05) between the unstressed and heat-shocked S. enterica. Unlike the wild type, temperature abuse did not lead to the enhanced survival of the ΔrpoS mutant after exposure to desiccation stress, indicating that the rpoS gene could play a critical role in the persistence of desiccation-stressed S. enterica subjected to temperature abuse. E. faecium NRRL B-2354 was more persistent (p < 0.05) than the pre-stressed S. enterica under both conditions, suggesting its use as a suitable surrogate for pre-stressed S. enterica by providing a sufficient safety margin. Our results demonstrate the merit of considering the prior exposure of foodborne pathogens to sub-lethal stresses when validating the storage conditions for leafy greens.

Cattle are the main reservoir of Shiga toxin-producing Escherichia coli (STEC), one of the world’s most important foodborne pathogens. The pathogen causes severe human diseases and outbreaks. This study aimed to identify and characterize non-O157 STEC isolated from cattle feces from central and southern Chile. We analyzed 446 cattle fecal samples and isolated non-O157 STEC from 12.6% (56/446); a total of 93 different isolates were recovered. Most isolates displayed β-glucuronidase activity (96.8%; 90/93) and fermented sorbitol (86.0%; 80/93), whereas only 39.8% (37/93) were resistant to tellurite. A subgroup of 30 representative non-O157 STEC isolates was selected for whole-genome sequencing and bioinformatics analysis. In silico analysis showed that they grouped into 16 different sequence types and 17 serotypes; the serotypes most frequently identified were O116:H21 and O168:H8 (13% each). A single isolate of serotype O26:H11 was recovered. One isolate was resistant to tetracycline and carried resistance genes tet(A) and tet(R); no other isolate displayed antimicrobial resistance or carried antimicrobial resistance genes. The intimin gene (eae) was identified in 13.3% (4/30) of the genomes and 90% (27/30) carried the stx2 gene. A phylogenetic reconstruction demonstrated that the isolates clustered based on serotypes, independent of geographical origin. These results indicate that cattle in Chile carry a wide diversity of STEC potentially pathogenic for humans based on the presence of critical virulence genes.

Pickled ready-to-eat meat (PRTEM) is a meat product that is treated with various seasonings and then cooked. PRTEM is a popular food consumed mostly in China and some Asian countries. Since this food is considered ‘ready to eat’, once it is contaminated by foodborne pathogens such as Salmonella, the prospect for significant morbidity, mortality, and immeasurable economic losses can occur. Here we investigated the prevalence and concentration of Salmonella in 107 PRTEM samples collected from Shaanxi, China during 2015–2016. Furthermore, we analyzed the serotype, antibiotic susceptibility, and presence of antibiotic resistance genes and amino acid mutations in 219 Salmonella isolates, followed by subtyping of 115 representative isolates. The average detection rate of Salmonella-positive PRTEM was 58.9%, and the average most probable number (MPN) of Salmonella in positive samples was 2.27 logMPN per gram of sample (range: 2.10–2.43). Ten serotypes were identified from the 219 Salmonella isolates, with S. Thompson (37.9%) and S. Indiana (20.5%) being predominant. The remaining serotypes were S. Typhi (7.8%), S. Typhimurium (7.3%), S. Mbandaka (6.9%), S. Albany (6.4%), S. Blockley (5.5%), S. Infantis (4.1%), S. Escanaba (3.2%), and S. Dusseldorf (0.5%). All isolates were resistant to ceftiofur (100%), while most of them were resistant to ciprofloxacin (99.1%), amoxicillin-clavulanic acid (97.7%), trimethoprim-sulfamethoxazole (96.4%), ampicillin (92.3%), sulfisoxazole (92.2%), tetracyclines (90.4%), and nalidixic acid (90.4%), respectively. A single mutation of Ser83Phe (27.1%) and double mutations of Ser83Phe-Asp87Gly (25.9%) in GyrA were detected in 85 isolates, whereas mutations of Thr57Ser (63.9%) and Ser80Arg (36.1%) in ParC were detected in 122 isolates. qnrB, oqxAB, aac(6′)-Ib, and qnrA were present in 50 (22.8%), 48 (21.9%), 26 (11.9%), and 1 (0.5%) isolate(s), respectively. Pulse field gel electrophoresis results revealed that those isolates recovered from the same type of PRTEM or the same sampling place shared identical or similar DNA profiles, antibiotic resistance phenotypes, and even plasmid-mediated quinolone resistance encoding genes. The findings indicate that Salmonella is commonly prevalent in PRTEMs at high concentrations in Shaanxi, China. More attention should be paid to the processing and storage of this ready-to-eat food to prevent bacterial contamination and foodborne outbreaks.

Escherichia coli O145 serogroup is one of the big six non-O157 Shiga toxin producing E. coli (STEC) that causes foodborne illnesses in the United States and other countries. Cattle are a major reservoir of STEC, which harbor them in their hindgut and shed in the feces. Cattle feces is the main source of hide and subsequent carcass contaminations during harvest leading to foodborne illnesses in humans. The objective of our study was to determine the virulence potential of STEC O145 strains isolated from cattle feces and hide samples. A total of 71 STEC O145 strains isolated from cattle feces (n = 16), hide (n = 53), and human clinical samples (n = 2) were used in the study. The strains were subjected to whole genome sequencing using Illumina MiSeq platform. The average draft genome size of the fecal, hide, and human clinical strains were 5.41, 5.28, and 5.29 Mb, respectively. The average number of genes associated with mobile genetic elements was 260, 238, and 259, in cattle fecal, hide, and human clinical strains, respectively. All strains belonged to O145:H28 serotype and carried eae subtype γ. Shiga toxin 1a was the most common Shiga toxin gene subtype among the strains, followed by stx2a and stx2c. The strains also carried genes encoding type III secretory system proteins, nle, and plasmid-encoded virulence genes. Phylogenetic analysis revealed clustering of cattle fecal strains separately from hide strains, and the human clinical strains were more closely related to the hide strains. All the strains belonged to sequence type (ST)-32. The virulence gene profile of STEC O145 strains isolated from cattle sources was similar to that of human clinical strains and were phylogenetically closely related to human clinical strains. The genetic analysis suggests the potential of cattle STEC O145 strains to cause human illnesses. Inclusion of more strains from cattle and their environment in the analysis will help in further elucidation of the genetic diversity and virulence potential of cattle O145 strains.

ABSTRACT Using whole-genome sequence (WGS) data from the GenomeTrakr network, a globally distributed network of laboratories sequencing foodborne pathogens, we present a new phylogeny of Salmonella enterica comprising 445 isolates from 266 distinct serovars and originating from 52 countries. This phylogeny includes two previously unidentified S. enterica subsp. enterica clades. Serovar Typhi is shown to be nested within clade A. Our findings are supported by both phylogenetic support, based on a core genome alignment, and Bayesian approaches, based on single-nucleotide polymorphisms. Serovar assignments were refined by in silico analysis using SeqSero. More than 10% of serovars were either polyphyletic or paraphyletic. We found variable genetic content in these isolates relating to gene mobilization and virulence factors which have different distributions within clades. Gifsy-1- and Gifsy-2-like phages appear more prevalent in clade A; other viruses are more evenly distributed. Our analyses reveal IncFII is the predominant plasmid replicon in S. enterica. Few core or clade-defining virulence genes are observed, and their distributions appear probabilistic in nature. Together, these patterns demonstrate that genetic exchange within S. enterica is more extensive and frequent than previously realized, which significantly alters how we view the genetic structure of the bacterial species. IMPORTANCE Rapid improvements in nucleotide sequencing access and affordability have led to a drastic increase in availability of genetic information. This information will improve the accuracy of molecular descriptions, including serovars, within S. enterica. Although the concept of serovars continues to be useful, it may have more significant limitations than previously understood. Furthermore, the discrete absence or presence of specific genes can be an unstable indicator of phylogenetic identity. Whole-genome sequencing provides more rigorous tools for assessing the distributions of these genes. Our phylogenetic and genetic content analyses reveal how active genetic elements are dynamically distributed within a species, allowing us to better understand genetic reservoirs and underlying bacterial evolution.

Recently developed nanopore sequencing technologies offer a unique opportunity to rapidly close the genome and to identify complete sequences of mobile genetic elements (MGEs). In this study, 17 isolates of Listeria monocytogenes (Lm) epidemic clone II (ECII) from seven ready-to-eat meat or poultry processing facilities, not known to be associated with outbreaks, were shotgun sequenced, and among them, five isolates were further subjected to long-read sequencing. Additionally, 26 genomes of Lm ECII isolates associated with three listeriosis outbreaks in the U.S. and South Africa were obtained from the National Center for Biotechnology Information (NCBI) database and analyzed to evaluate if MGEs may be used as a high-resolution genetic marker for identifying and sourcing the origin of Lm. The analyses identified four comK prophages in 11 non-outbreak isolates from four facilities and three comK prophages in 20 isolates associated with two outbreaks that occurred in the U.S. In addition, three different plasmids were identified among 10 non-outbreak isolates and 14 outbreak isolates. Each comK prophage and plasmid was conserved among the isolates sharing it. Different prophages from different facilities or outbreaks had significant genetic variations, possibly due to horizontal gene transfer. Phylogenetic analysis showed that isolates from the same facility or the same outbreak always closely clustered. The time of most recent common ancestor of the Lm ECII isolates was estimated to be in March 1816 with the average nucleotide substitution rate of 3.1 × 10−7 substitutions per site per year. This study showed that complete MGE sequences provide a good signal to determine the genetic relatedness of Lm isolates, to identify persistence or repeated contamination that occurred within food processing environment, and to study the evolutionary history among closely related isolates.

To obtain full details of gut microbiota, including bacteria, fungi, bacteriophages, and helminths, in giant pandas (GPs), we created a comprehensive microbial genome database and used metagenomic sequences to align against the database. We delineated a detailed and different gut microbiota structures of GPs. A total of 680 species of bacteria, 198 fungi, 185 bacteriophages, and 45 helminths were found. Compared with 16S rRNA sequencing, the dominant bacterium phyla not only included Proteobacteria, Firmicutes, Bacteroidetes, and Actinobacteria but also Cyanobacteria and other eight phyla. Aside from Ascomycota, Basidiomycota, and Glomeromycota, Mucoromycota, and Microsporidia were the dominant fungi phyla. The bacteriophages were predominantly dsDNA Myoviridae, Siphoviridae, Podoviridae, ssDNA Inoviridae, and Microviridae. For helminths, phylum Nematoda was the dominant. In addition to previously described parasites, another 44 species of helminths were found in GPs. Also, differences in abundance of microbiota were found between the captive, semiwild, and wild GPs. A total of 1,739 genes encoding cellulase, β-glucosidase, and cellulose β-1,4-cellobiosidase were responsible for the metabolism of cellulose, and 128,707 putative glycoside hydrolase genes were found in bacteria/fungi. Taken together, the results indicated not only bacteria but also fungi, bacteriophages, and helminths were diverse in gut of giant pandas, which provided basis for the further identification of role of gut microbiota. Besides, metagenomics revealed that the bacteria/fungi in gut of GPs harbor the ability of cellulose and hemicellulose degradation.

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.

Escherichia coli O103, harbored in the hindgut and shed in the feces of cattle, can be enterohemorrhagic (EHEC), enteropathogenic (EPEC), or putative non-pathotype. The genetic diversity particularly that of virulence gene profiles within O103 serogroup is likely to be broad, considering the wide range in severity of illness. However, virulence descriptions of the E. coli O103 strains isolated from cattle feces have been primarily limited to major genes, such as Shiga toxin and intimin genes. Less is known about the frequency at which other virulence genes exist or about genes associated with the mobile genetic elements of E. coli O103 pathotypes. Our objective was to utilize whole genome sequencing (WGS) to identify and compare major and putative virulence genes of EHEC O103 (positive for Shiga toxin gene, stx1, and intimin gene, eae; n = 43), EPEC O103 (negative for stx1 and positive for eae; n = 13) and putative non-pathotype O103 strains (negative for stx and eae; n = 13) isolated from cattle feces. Six strains of EHEC O103 from human clinical cases were also included. All bovine EHEC strains (43/43) and a majority of EPEC (12/13) and putative non-pathotype strains (12/13) were O103:H2 serotype. Both bovine and human EHEC strains had significantly larger average genome sizes (P < 0.0001) and were positive for a higher number of adherence and toxin-based virulence genes and genes on mobile elements (prophages, transposable elements, and plasmids) than EPEC or putative non-pathotype strains. The genome size of the three pathotypes positively correlated (R2 = 0.7) with the number of genes carried on mobile genetic elements. Bovine strains clustered phylogenetically by pathotypes, which differed in several key virulence genes. The diversity of E. coli O103 pathotypes shed in cattle feces is likely reflective of the acquisition or loss of virulence genes carried on mobile genetic elements.

The DNA sequence of the O-antigen gene clusters of Escherichia coli serogroups O62, O68, O131, O140, O142, and O163 was determined, and primers based on the wzx (O-antigen flippase) and/or wzy (O-antigen polymerase) genes within the O-antigen gene clusters were designed and used in PCR assays to identify each serogroup. Specificity was tested with E. coli reference strains, field isolates belonging to the target serogroups, and non-E. coli bacteria. The PCR assays were highly specific for the respective serogroups; however, the PCR assay targeting the O62 wzx gene reacted positively with strains belonging to E. coli O68, which was determined by serotyping. Analysis of the O-antigen gene cluster sequences of serogroups O62 and O68 reference strains showed that they were 94% identical at the nucleotide level, although O62 contained an insertion sequence (IS) element located between the rmlA and rmlC genes within the O-antigen gene cluster. A PCR assay targeting the rmlA and rmlC genes flanking the IS element was used to differentiate O62 and O68 serogroups. The PCR assays developed in this study can be used for the detection and identification of E. coli O62/O68, O131, O140, O142, and O163 strains isolated from different sources.

Listeria monocytogenes is a significant foodborne pathogen causing severe systemic infections in humans with high mortality rates. The objectives of this work were to establish a phylogenetic framework of L. monocytogenes from China and to investigate sequence diversity among different serotypes. We selected 17 L. monocytogenes strains recovered from patients and foods in China representing serotypes 1/2a, 1/2b, and 1/2c. Draft genome sequences were determined using Illumina MiSeq technique and associated protocols. Open reading frames were assigned using prokaryotic genome annotation pipeline by NCBI. Twenty-four published genomes were included for comparative genomic and phylogenetic analysis. More than 154,000 single nucleotide polymorphisms (SNPs) were identified from multiple genome alignment and used to reconstruct maximum likelihood phylogenetic tree. The 41 genomes were differentiated into lineages I and II, which consisted of 4 and 11 subgroups, respectively. A clinical strain from China (SHL009) contained significant SNP differences compared to the rest genomes, whereas clinical strain SHL001 shared most recent common ancestor with strain SHL017 from food. Moreover, clinical strains SHL004 and SHL015 clustered together with two strains (08-5578 and 08-5923) recovered from an outbreak in Canada. Partial sequences of a plasmid found in the Canadian strain were also present in SHL004. We investigated the presence of various genes and gene clusters associated with virulence and subgroup-specific genes, including internalins, L. monocytogenes pathogenicity islands (LIPIs), L. monocytogenes genomic islands (LGIs), stress survival islet 1 (SSI-1), and clustered regularly interspaced short palindromic repeats (CRISPR)/cas system. A novel genomic island, denoted as LGI-2 was identified. Comparative sequence analysis revealed differences among the L. monocytogenes strains related to virulence, survival abilities, and attributes against foreign genetic elements. L. monocytogenes from China were genetically diverse. Strains from clinical specimens and food related closely suggesting foodborne transmission of human listeriosis.

To improve pulsed-field gel electrophoresis–based strain discrimination of 76 Salmonella Enteritidis strains, we evaluated 6 macro-restriction endonucleases, separately and in various combinations. One 3-enzyme subset, SfiI/PacI/NotI, was highly discriminatory. Five different indices, including the Simpson diversity index, supported this 3-enzyme combination for improved differentiation of S. Enteritidis.

Susceptibilities to fourteen antimicrobial agents important in clinical medicine and agriculture were determined for 752 Escherichia coli isolates of serotypes O26, O103, O111, O128, and O145. Strains of these serotypes may cause urinary tract and enteric infections in humans and have been implicated in infections with Shiga toxin–producing E. coli (STEC). Approximately 50% of the 137 isolates from humans were resistant to ampicillin, sulfamethoxazole, cephalothin, tetracycline, or streptomycin, and approximately 25% were resistant to chloramphenicol, trimethoprim-sulfamethoxazole, or amoxicillin-clavulanic acid. Approximately 50% of the 534 isolates from food animals were resistant to sulfamethoxazole, tetracycline, or streptomycin. Of 195 isolates with STEC-related virulence genes, approximately 40% were resistant to sulfamethoxazole, tetracycline, or streptomycin. Findings from this study suggest antimicrobial resistance is widespread among E. coli O26, O103, O111, O128, and O145 inhabiting humans and food animals.

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.

Overexpression of ramA has been implicated in resistance to multiple drugs in several enterobacterial pathogens. In the present study, Salmonella Typhimurium strain LTL with constitutive expression of ramA was compared to its ramA-deletion mutant by employing both DNA microarrays and phenotype microarrays (PM). The mutant strain with the disruption of ramA showed differential expression of at least 33 genes involved in 11 functional groups. The study confirmed at the transcriptional level that the constitutive expression of ramA was directly associated with increased expression of multidrug efflux pump AcrAB-TolC and decreased expression of porin protein OmpF, thereby conferring multiple drug resistance phenotype. Compared to the parent strain constitutively expressing ramA, the ramA mutant had increased susceptibility to over 70 antimicrobials and toxic compounds. The PM analysis also uncovered that the ramA mutant was better in utilization of 10 carbon sources and 5 phosphorus sources. This study suggested that the constitutive expression of ramA locus regulate not only multidrug efflux pump and accessory genes but also genes involved in carbon metabolic pathways.

Several studies are conducted to determine the presence of the Shiga toxin-producing Escherichiacoli (STEC) and other potentially diarrheagenic E. coli strains in retail meats. The results demonstrate the analyzed retail meats are mostly contaminated via different STEC strains.

The aim of this study was to investigate the toxin gene profile and antimicrobial resistance of Staphylococcus aureus isolates from raw chicken in the People's Republic of China. In total, 289 S. aureus isolates were characterized by antimicrobial susceptibility testing, and genes encoding enterotoxins, exfoliative toxins, Panton-Valentine leukocidin, and toxic shock syndrome toxin were revealed by PCR. Overall, 46.0% of the isolates were positive for one or more toxin genes. A high proportion of toxin genes were pvl (26.6%), followed by sej (12.5%), sea (9.0%), seh (8.3%), seb (6.9%), sec (6.9%), sed (4.8%), sei (3.1%), and see (2.4%). None of the isolates harbored seg, tsst-1, or exfoliative toxin genes. In total, 29 toxin gene profiles were obtained, and pvl (10.7%) was the most frequent genotype, followed by sea (5.9%), seb (4.8%), and sej (4.2%). Furthermore, 99.7% of the strains were resistant to at least one of the tested antimicrobial agents, and 87.2% of them displayed multidrug resistance. Resistance was most frequently observed to trimethoprim-sulfamethoxazole and erythromycin (86.2% for each), followed by tetracycline (69.9%), amoxicillin--clavulanic acid (45.0%), and ampicillin (42.6%). None of the strains were resistant to vancomycin. This study indicates that S. aureus isolates from raw chicken harbored multiple toxin genes and exhibited multiple antimicrobial resistance, which represents a potential health hazard for consumers. [ABSTRACT FROM AUTHOR]