Dynamics of microbiota and antimicrobial resistance in on-farm dairy processing plants using metagenomic and culture-dependent approaches.

On-farm dairy processing plants, which are situated close to farms and larger dairy processing facilities, face unique challenges in maintaining environmental hygiene. This can impact various stages of dairy processing. These plants operate on smaller scales and use Low-Temperature-Long-Time (LTLT) pasteurization, making them more susceptible to microbial contamination through direct and indirect contact. Antimicrobial-resistant bacteria found on dairy farms pose risks to human health by potentially transferring resistance via dairy products. Our study aimed to investigate microbial distribution and antimicrobial resistance at four key stages: the farm, pre-pasteurization, post-pasteurization, and processing environments. We assessed microbial distribution by quantifying indicator bacteria and conducting metagenomic analysis. Antimicrobial resistance was examined by identifying resistance phenotypes and detecting resistance genes in bacterial isolates and metagenomes. Our results showed that the indicator bacteria were detected at all stages of on-farm dairy processing. We observed a significant reduction in aerobic microbes and coliforms post-pasteurization. However, contamination of the final dairy products increased, suggesting potential cross-contamination during post-pasteurization. Metagenomic analysis revealed that Pseudomonas, a representative psychrotrophic bacterium, was predominant in both the farm (24.1 %) and pre-pasteurization (65.9 %) stages, indicating microbial transfer from the farms to the processing plants. Post-pasteurization, Pseudomonas and other psychrotrophs like Acinetobacter and Enterobacteriaceae remained dominant. Core microbiota analysis identified 74 genera in total, including 13 psychrotrophic bacteria, across all stages. Of the 59 strains isolated from these plants, 49 were psychrotrophic. Antimicrobial resistance analysis showed that 74.6 % (44/59) of isolates were resistant to at least one antibiotic, with cefoxitin-, ampicillin-, amoxicillin-, and ticarcillin-resistant bacteria present at all stages. Identical antimicrobial resistance patterns were observed in isolates from serial stages of the same farm and season, suggesting bacterial transmission across stages. Additionally, 27.1 % (16/59) of isolates carried plasmid-mediated resistance genes, which were also detected in the metagenomes of non-isolated samples, indicating potential antimicrobial resistance gene transmission and their presence in uncultured bacteria. These findings reveal the persistence of antimicrobial-resistant psychrotrophic bacteria in on-farm dairy processing plants, which pose potential health risks via dairy consumption. Our study underscores the importance of both culture-dependent and culture-independent methods to fully understand their distribution and impact.

Effect of biofilm formation by antimicrobial-resistant gram-negative bacteria in cold storage on survival in dairy processing lines.

Antimicrobial-resistant gram-negative bacteria in dairy products can transfer antimicrobial resistance to gut microbiota in humans and can adversely impact the product quality. In this study, we aimed to investigate their distribution in dairy processing lines and evaluate biofilm formation and heat tolerance under dairy processing line-like conditions. Additionally, we compared the relative expression of general and heat stress-related genes as well as spoilage-related gene between biofilm and planktonic cells under consecutive stresses, similar to those in dairy processing lines. Most species of gram-negative bacteria isolated from five different dairy processing plants were resistant to one or more antimicrobials. Biofilm formation by the bacteria at 5 °C increased with the increase in exposure time. Moreover, cells in biofilms remained viable under heat treatment, whereas all planktonic cells of the selected strains died. The expression of heat-shock-related genes significantly increased with heat treatment in the biofilms but mostly decreased in the planktonic cells. Thus, biofilm formation under raw milk storage conditions may improve the tolerance of antimicrobial-resistant gram-negative bacteria to pasteurization, thereby increasing their persistence in dairy processing lines and products. Furthermore, the difference in response to heat stress between biofilm and planktonic cells may be attributed to the differential expression of heat stress-related genes. Therefore, this study contributes to the understanding of how gram-negative bacteria persist under consecutive stresses in dairy processing procedures and the potential mechanism underlying heat tolerance in biofilms.

Circular intermediate-mediated horizontal transfer of the chromosome-encoded cfr(C) gene in multi-drug resistant Campylobacter coli from swine sources.

Campylobacter is a major zoonotic pathogen that causes gastrointestinal and, rarely, immune diseases in humans. The antimicrobial-resistance gene cfr(C) carried by Campylobacter and is a cfr-like gene that targets bacterial 23S rRNA through A2503 methylation. cfr(C) confers cross-resistance to five antimicrobial classes (PhLOPSA), including lincosamide, streptogramin A, and pleuromutilin, which are classified as critically important antimicrobials to human by the World Health Organization. To elucidate the genetic variation and horizontal transfer mechanism of cfr(C), we analyzed the genetic background and horizontal transfer unit of Campylobacter-derived cfr(C) through comparative genomic analysis. We identified nine cfr(C)-positive C. coli strains of 157 strains isolated from swine sources. Three novel cfr(C) gene single nucleotide polymorphism (SNP) sites (19delA, 674C > A, and 890 T > C) were identified from nine cfr(C)-positive strains. Among six identified cfr(C) SNP variant types (SNP-I to -VI), five types of randomly inserted cfr(C)-cassettes on chromosome and one type of plasmid-like element were identified, their gene cassette composition differing depending on the cfr(C) variants. Three of six cfr(C) cassette types contained aminoglycoside-streptothricin resistance cluster "aphA3-sat4-aadE." The cfr(C) gene cassette with pcp gene (GC-1, GC-4, and GC-5) formed a pcp-mediated circular intermediate "pcp-hp-cfr(C)-aphA3," which has not been previously reported. Other two cfr(C) cassette-types with ISChh1 formed circular intermediate "ISChh1-aphA3-cfr(C)-lnu (G)-pnp-ant1-hp-ATPase" and "ISChh1-aphA3-cfr(C)-hp." In conjugation assay, the pcp-mediated circular intermediate was naturally transferred to the plasmid of recipient C. coli wild-type strain from swine source, and comparative genomic analysis revealed that cfr(C) encoded in pcp-mediated circular intermediate was inserted into the plasmid of recipient by homologous recombination with pcp and aphA3. This study revealed that novel multidrug resistance gene cfr(C) carried by C. coli from swine sources can be highly genetically diverse and transferable. Moreover, we suggest that the transferability of chromosomal cfr(C) may contribute to the global spread of multidrug resistance against clinically important antimicrobials.

Genomic characteristics and comparative genomics analysis of Salmonella enterica subsp. enterica serovar Thompson isolated from an outbreak in South Korea.

Salmonella infections represent an important public health problem. In 2018, a multistate outbreak of S. enterica subsp. enterica serovar Thompson infection associated with contaminated chocolate cakes in schools was reported in South Korea. In this study, we sequenced the 37 S. Thompson strains isolated from chocolate cakes, egg whites, preserves, and cookware associated with the outbreak. In addition, we analyze the genomic sequences of 61 S. Thompson strains (37 chocolate cake-related outbreak strains, 4 strains isolated from outbreaks in South Korea and 20 strains available in the National Center for Biotechnology Information) to assess the genomic characteristics of outbreak-related strains by comparative genomics and phylogenetic analysis. The results showed that identically classified clusters divided strains into two clusters, sub-clusters A & I (with strains from 2018 in South Korea) and sub-clusters B & II (with strains from 2014 to 2015 in South Korea). S. Thompson isolated from South Korea were accurately distinguished from publicly-available strains. Unlike other S. Thompson genomes, those of chocolate cake outbreak-related strains had three Salmonella phages (SEN8, vB SosS Oslo, and SI7) integrated into their chromosome. Comparative genomics revealed several genes responsible for the specific genomic features of chocolate cake outbreak-related strains and three bacteriophages that may contribute to the pathogenicity of other S. Thompson strains.

Characterization of Cronobacter sakazakii Strains Originating from Plant-Origin Foods Using Comparative Genomic Analyses and Zebrafish Infectivity Studies.

Cronobacter sakazakii continues to be isolated from ready-to-eat fresh and frozen produce, flours, dairy powders, cereals, nuts, and spices, in addition to the conventional sources of powdered infant formulae (PIF) and PIF production environments. To understand the sequence diversity, phylogenetic relationship, and virulence of C. sakazakii originating from plant-origin foods, comparative molecular and genomic analyses, and zebrafish infection (ZI) studies were applied to 88 strains. Whole genome sequences of the strains were generated for detailed bioinformatic analysis. PCR analysis showed that all strains possessed a pESA3-like virulence plasmid similar to reference C. sakazakii clinical strain BAA-894. Core genome analysis confirmed a shared genomic backbone with other C. sakazakii strains from food, clinical and environmental strains. Emerging nucleotide diversity in these plant-origin strains was highlighted using single nucleotide polymorphic alleles in 2000 core genes. DNA hybridization analyses using a pan-genomic microarray showed that these strains clustered according to sequence types (STs) identified by multi-locus sequence typing (MLST). PHASTER analysis identified 185 intact prophage gene clusters encompassing 22 different prophages, including three intact Cronobacter prophages: ENT47670, ENT39118, and phiES15. AMRFinderPlus analysis identified the CSA family class C ß-lactamase gene in all strains and a plasmid-borne mcr-9.1 gene was identified in three strains. ZI studies showed that some plant-origin C. sakazakii display virulence comparable to clinical strains. Finding virulent plant-origin C. sakazakii possessing significant genomic features of clinically relevant STs suggests that these foods can serve as potential transmission vehicles and supports widening the scope of continued surveillance for this important foodborne pathogen.

Prevalence, Characteristics, and Clonal Distribution of Escherichia coli Carrying Mobilized Colistin Resistance Gene mcr-1.1 in Swine Farms and Their Differences According to Swine Production Stages.

Global spread of Escherichia coli strains carrying the mobilized colistin resistance gene mcr-1.1 (MCR1-EC) poses serious threats to public health. Colistin has been generally prescribed for swine colibacillosis, having made swine farms as major reservoirs of MCR1-EC. The present study aimed to understand characteristic differences of MCR1-EC, including prevalence, antimicrobial resistance, and virulence, according to swine production stages. In addition, genetic relatedness was evaluated between MCR1-EC isolated from this study as well as pig-, human-, and chicken-derived strains published in the National Center for Biotechnology Information (NCBI), based on the multi-locus sequence types (MLSTs) and whole-genome sequences (WGS). Individual fecal samples (n = 331) were collected from asymptomatic weaning-piglets, growers, finishers, and sows from 10 farrow-to-finishing farms in South Korea between 2017 and 2019. The weighted prevalence of MCR1-EC was 11.6% (95% CI: 8.9%-15.0%, 55/331), with the highest prevalence at weaning stage. The 96.2% of MCR1-EC showed multi-drug resistance. Notably, weaning stage-derived MCR1-EC showed higher resistance rates (e.g., against extended-spectrum ß-lactams or quinolones) than those from other stages. MCR1-EC with virulence advantages (e.g., intestinal/extraintestinal pathogenic E. coli or robust biofilm formation) were identified from all pig stages, accounting for nearly half of the total strains. WGS-based in-depth characterization showed that intestinal pathogenic MCR1-EC harbored multi-drug resistance and multiple virulence factors, which were highly shared between strains isolated from pigs of different stages. The clonal distribution of MCR1-EC was shared within swine farms but rarely across farms. The major clonal type of MCR1-EC from swine farms and NCBI database was ST10-A. Core genomes of MCR1-EC isolated from individuals within closed environments (same farms or human hospitals) were highly shared (genetic distance < 0.01), suggesting a high probability of clonal expansion of MCR1-EC within closed environments such as livestock husbandry. To the best of our knowledge, this is the first study to analyze the differences in the characteristics and clonal distribution of MCR1-EC according to production stages in swine farms, an important reservoir of MCR1-EC. Our results highlight the need to establish MCR1-EC control plans in swine farms based on an in-depth understanding of MCR1-EC characteristics according to swine production stages, focusing especially on the weaning stages.

Hyper-aerotolerant Campylobacter coli, an emerging foodborne pathogen, shows differential expressions of oxidative stress-related genes.

Although Campylobacter, an obligate microaerophilic foodborne pathogen, is susceptible to oxygen, aerotolerant/hyper-aerotolerant (HAT) Campylobacter can survive under aerobic conditions. Here, we aimed to reveal what affects the enhanced aerotolerance in HAT Campylobacter coli at genome and gene expression levels. We compared the whole genomes between HAT and oxygen-sensitive (OS) C. coli isolates from swine and analyzed the relative expressions of oxidative stress-related (sodB, ahpC, katA, and trxB) and iron transport/uptake-related (cfbpA, ceuE, feuB, and feoB) genes. The comparative genomics showed no relation between the clustering of the strains and aerotolerance levels. The reactive oxygen species-related factors involved in respiration, stress response, and iron acquisition/uptake were similar among the strains, regardless of their aerotolerance levels. However, the expressions of the oxidative stress-related genes under aerobic conditions compared to that of microaerobic conditions increased in the HAT strains, while decreased in the OS strains. Our findings suggest that what influences differences in aerotolerance between HAT and OS C. coli may be due to the differential expressions of oxidative stress-related genes despite the similarities in genomic structure. This study provides insights into the genetic basis of aerotolerance in C. coli. Therefore, it could assist in managing HAT C. coli that has the potential to be easily transmitted to humans through the food chain.

Emergence of Transferable mcr-9 Gene-Carrying Colistin-Resistant Salmonella enterica Dessau ST14 Isolated from Retail Chicken Meat in Korea.

Colistin is an important antibiotic currently used to manage infections caused by multidrug-resistant pathogens in both humans and livestock animals. A new mobile colistin-resistance (mcr-9) gene was recently discovered; this discovery highlighted the need for rigorous monitoring of bacterial resistance against colistin. Salmonella is one of the major pathogens responsible for foodborne illnesses; however, there is minimal information regarding the presence of mcr genes in foodborne Salmonella strains. The aim of this study was to investigate the presence of mcr genes among 178 Salmonella strains isolated from chicken meat in Korea. Antimicrobial susceptibility was measured using the broth microdilution method. Bioinformatics characterization of colistin-resistant strains and genetic environment of the mcr-9 gene were analyzed using next-generation sequencing. Transferability of the mcr-9 carrying colistin-resistant Salmonella strain was tested using broth-mating conjugation. Thirteen of the 178 Salmonella isolates showed colistin resistance, but only one strain, Salmonella Dessau ST14 (KUFSE-SAL043) from a traditional chicken market in Korea, carried an mcr family gene, mcr-9. This strain also carried other acquired antimicrobial resistance genes such as blaTEM-1B, qnrS1, and aac(6')-Iaa. Only the IncX1 plasmid replicon type was detected in this strain. In the strain KUFSE-SAL043, the mcr-9 gene was located between two insertion sequences, IS903B and IS26, followed by the downstream regulatory genes qseB-like and qseC-like, which were located between IS1R and ΔIS1R. Conjugation tests revealed that the mcr-9 gene was successfully transferred to Escherichia coli J53 at a mean frequency of 2.03 × 10-7. This is the first report of a transferable mcr-9 gene in Salmonella isolated from chicken meat in Korea, highlighting the possibility of transfer of colistin resistance. Therefore, the wide use of colistin should be reconsidered, and a One Health perspective should be adopted to monitor the antimicrobial resistance of Enterobacteriaceae strains in humans, livestock, and the environment.

The Secretion of Toxins and Other Exoproteins of Cronobacter: Role in Virulence, Adaption, and Persistence.

: Cronobacter species are considered an opportunistic group of foodborne pathogenic bacteria capable of causing both intestinal and systemic human disease. This review describes common virulence themes shared among the seven Cronobacter species and describes multiple exoproteins secreted by Cronobacter, many of which are bacterial toxins that may play a role in human disease. The review will particularly concentrate on the virulence factors secreted by C. sakazakii, C. malonaticus, and C. turicensis, which are the primary human pathogens of interest. It has been discovered that various species-specific virulence factors adversely affect a wide range of eukaryotic cell processes including protein synthesis, cell division, and ion secretion. Many of these factors are toxins which have been shown to also modulate the host immune response. These factors are encoded on a variety of mobile genetic elements such as plasmids and transposons; this genomic plasticity implies ongoing re-assortment of virulence factor genes which has complicated our efforts to categorize Cronobacter into sharply defined genomic pathotypes.

Prevalence, Distribution, and Phylogeny of Type Two Toxin-Antitoxin Genes Possessed by Cronobacter Species where C. sakazakii Homologs Follow Sequence Type Lineages.

Cronobacter species are a group of foodborne pathogenic bacteria that cause both intestinal and systemic human disease in individuals of all age groups. Little is known about the mechanisms that Cronobacter employ to survive and persist in foods and other environments. Toxin-antitoxin (TA) genes are thought to play a role in bacterial stress physiology, as well as in the stabilization of horizontally-acquired re-combinatorial elements such as plasmids, phage, and transposons. TA systems have been implicated in the formation of a persistence phenotype in some bacterial species including Escherichia coli and Salmonella. This project's goal was to understand the phylogenetic relatedness among TA genes present in Cronobacter. Preliminary studies showed that two typical toxin genes, fic and hipA followed species evolutionary lines. A local database of 22 TA homologs was created for Cronobacter sakazakii and a Python version 3 shell script was generated to extract TA FASTA sequences present in 234 C. sakazakii genomes previously sequenced as part of Center for Food Safety and Applied Nutrition's (CFSAN) GenomeTrakr project. BLAST analysis showed that not every C. sakazakii strain possessed all twenty-two TA loci. Interestingly, some strains contained either a toxin or an antitoxin component, but not both. Five common toxin genes: ESA_00258 (parDE toxin-antitoxin family), ESA_00804 (relBE family), ESA_01887 (relBE family), ESA_03838 (relBE family), and ESA_04273 (YhfG-Fic family) were selected for PCR analysis and the primers were designed to detect these genes. PCR analysis showed that 55 of 63 strains possessed three of these genes Sequence analysis identified homologs of the target genes and some of the strains were PCR-negative for one or more of the genes, pointing to potential nucleotide polymorphisms in those loci or that these toxin genes were absent. Phylogenetic studies using a Cronobacter pan genomic microarray showed that for the most part TAs follow species evolutionary lines except for a few toxin genes possessed by some C. malonaticus and C. universalis strains; this demonstrates that some TA orthologues share a common phylogeny. Within the C. sakazakii strains, the prevalence and distribution of these TA homologs by C. sakazakii strain BAA-894 (a powdered infant formula isolate) followed sequence-type evolutionary lineages. Understanding the phylogeny of TAs among the Cronobacter species is essential to design future studies to realize the physiological mechanisms and roles for TAs in stress adaptation and persistence of Cronobacter within food matrices and food processing environments.

Draft genomes of Cronobacter sakazakii strains isolated from dried spices bring unique insights into the diversity of plant-associated strains.

Cronobacter sakazakii is a Gram-negative opportunistic pathogen that causes life- threatening infantile infections, such as meningitis, septicemia, and necrotizing enterocolitis, as well as pneumonia, septicemia, and urinary tract and wound infections in adults. Here, we report 26 draft genome sequences of C. sakazakii, which were obtained from dried spices from the USA, the Middle East, China, and the Republic of Korea. The average genome size of the C. sakazakii genomes was 4393 kb, with an average of 4055 protein coding genes, and an average genome G + C content of 56.9%. The genomes contained genes related to carbohydrate transport and metabolism, amino acid transport and metabolism, and cell wall/membrane biogenesis. In addition, we identified genes encoding proteins involved in osmotic responses such as DnaJ, Aquaproin Z, ProQ, and TreF, as well as virulence-related and heat shock-related proteins. Interestingly, a metabolic island comprised of a variably-sized xylose utilization operon was found within the spice-associated C. sakazakii genomes, which supports the hypothesis that plants may serve as transmission vectors or alternative hosts for Cronobacter species. The presence of the genes identified in this study can support the remarkable phenotypic traits of C. sakazakii such as the organism's capabilities of adaptation and survival in response to adverse growth environmental conditions (e.g. osmotic and desiccative stresses). Accordingly, the genome analyses provided insights into many aspects of physiology and evolutionary history of this important foodborne pathogen.

Whole-Genome Sequences of Cronobacter sakazakii Isolates Obtained from Foods of Plant Origin and Dried-Food Manufacturing Environments.

Here, we present draft genome sequences of 29 Cronobacter sakazakii isolates obtained from foods of plant origin and dried-food manufacturing facilities. Assemblies and annotations resulted in genome sizes ranging from 4.3 to 4.5 Mb and 3,977 to 4,256 gene-coding sequences with G+C contents of ∼57.0%.

Genomic characterization of malonate positive Cronobacter sakazakii serotype O:2, sequence type 64 strains, isolated from clinical, food, and environment samples.

BACKGROUND: Malonate utilization, an important differential trait, well recognized as being possessed by six of the seven Cronobacter species is thought to be largely absent in Cronobacter sakazakii (Csak). The current study provides experimental evidence that confirms the presence of a malonate utilization operon in 24 strains of sequence type (ST) 64, obtained from Europe, Middle East, China, and USA; it offers explanations regarding the genomic diversity and phylogenetic relatedness among these strains, and that of other C. sakazakii strains. RESULTS: In this study, the presence of a malonate utilization operon in these strains was initially identified by DNA microarray analysis (MA) out of a pool of 347 strains obtained from various surveillance studies involving clinical, spices, milk powder sources and powdered infant formula production facilities in Ireland and Germany, and dried dairy powder manufacturing facilities in the USA. All ST64 C. sakazakii strains tested could utilize malonate. Zebrafish embryo infection studies showed that C. sakazakii ST64 strains are as virulent as other Cronobacter species. Parallel whole genome sequencing (WGS) and MA showed that the strains phylogenetically grouped as a separate clade among the Csak species cluster. Additionally, these strains possessed the Csak O:2 serotype. The nine-gene, ~ 7.7 kbp malonate utilization operon was located in these strains between two conserved flanking genes, gyrB and katG. Plasmidotyping results showed that these strains possessed the virulence plasmid pESA3, but in contrast to the USA ST64 Csak strains, ST64 Csak strains isolated from sources in Europe and the Middle East, did not possess the type six secretion system effector vgrG gene. CONCLUSIONS: Until this investigation, the presence of malonate-positive Csak strains, which are associated with foods and clinical cases, was under appreciated. If this trait was used solely to identify Cronobacter strains, many strains would likely be misidentified. Parallel WGS and MA were useful in characterizing the total genome content of these Csak O:2, ST64, malonate-positive strains and further provides an understanding of their phylogenetic relatedness among other virulent C. sakazakii strains.