Umeå University's logo

Salmonella enterica subspecies enterica serotype Minnesota (S. Minnesota) has recently emerged as a predominant serotype in poultry farming operations. Genomic surveillance efforts concentrated primarily in Europe have been used to evaluate food safety risks associated with S. Minnesota in imported poultry/poultry products. However, the burden imposed by S. Minnesota on consumers in sub-Saharan Africa is not understood. Here, we used whole-genome sequencing to characterize 36 S. Minnesota strains from raw poultry imported into South Africa, specifically (i) 11 strains isolated at port of entry and (ii) 25 strains from imported poultry in South African supermarkets. While all 36 strains belonged to the same sequence type (ST548), multiple ST548 lineages were present among poultry products. Comparison of the 36 strains sequenced here to all publicly available, high-quality ST548 genomes (n=460, from Enterobase) identified several public genomes differing by <30 core SNPs, including strains isolated previously from South American poultry imported into the UK. Notably, a cluster consisting of 14 highly similar genomes sequenced here (0 core SNPs) uniquely possessed blaCTX-M-8. A search of plasmids in public databases, alongside antimicrobial resistance (AMR) genes from >1.9 million bacterial genomes, revealed that this cluster harboured blaCTX-M-8 on an IncI1 plasmid-like region, which we hypothesize was acquired recently, from Escherichia coli. Overall, our study provides insight into the intercontinental dissemination of S. Minnesota and its associated AMR determinants via the global poultry trade.

The ability of Listeria to show reduced susceptibility to sanitizers commonly used in fresh produce packing and processing environments continues to be mentioned as a concern. We assessed the survival of 501 produce-associated Listeria isolates (328 Listeria monocytogenes [LM] and 173 Listeria spp. [LS]) after 30 s of exposure to benzalkonium chloride (BC, 300 ppm) and peroxyacetic acid (PAA, 80 ppm). A subset of 108 isolates was also exposed to sodium hypochlorite (NaOCl, 500 ppm) for 30 s. Isolates showed a range of log reductions, including 2.76–5.73 log for BC, 0.15–6.16 log for PAA, and 1.34–7.02 log for NaOCl; the variance of log reductions was significantly lower for BC compared to PAA and NaOCl. Cluster analysis on log reduction data identified four clusters, including one cluster of five LM isolates that showed reduced susceptibility to all three sanitizers. Log reductions of LS were significantly lower than LM after exposure to PAA, indicating reduced PAA susceptibility among LS. Whole genome sequence (WGS)-based characterization of all isolates revealed that the presence of known BC resistance genes (i.e., bcrABC, mdrL, and sugE1/2) was not significantly associated with log reductions to BC, and the presence of stress survival islet SSI-2 was not significantly associated with log reductions to PAA and NaOCl. Genome-wide association studies did not reveal any association of pangenome genes with phenotypic sanitizer susceptibility but identified several SNPs in core genes as associated with sanitizer susceptibility.

The average nucleotide identity (ANI) metric has become the gold standard for prokaryotic species delineation in the genomics era. The most popular ANI algorithms are available as command-line tools and/or web applications, making it inconvenient to incorporate them into bioinformatic workflows, which utilize the popular Python programming language. Here, we present PyOrthoANI, PyFastANI, and Pyskani, Python libraries for three popular ANI computation methods. ANI values produced by PyOrthoANI, PyFastANI, and Pyskani are virtually identical to those produced by OrthoANI, FastANI, and skani, respectively (adjusted R2>0.999). Compared to OrthoANI, PyOrthoANI is, on average, 3× faster per genome, while PyFastANI has multithreading support for single queries. All three libraries integrate seamlessly with BioPython, making it easy and convenient to use, compare, and benchmark popular ANI algorithms within Python-based bioinformatic workflows, software programs, and notebooks. Each library is available as part of the Python Package Index repository under the open-source MIT license, with source code available via GitHub (PyOrthoANI, https://github.com/althonos/orthoani; PyFastANI, https://github.com/althonos/pyfastani; Pyskani, https://github.com/althonos/pyskani).

Single-cell RNA-seq methods can be used to delineate cell types and states at unprecedented resolution but do little to explain why certain genes are expressed. Single-cell ATAC-seq and multiome (ATAC + RNA) have emerged to give a complementary view of the cell state. It is however unclear what additional information can be extracted from ATAC-seq data besides transcription factor binding sites. Here, we show that ATAC-seq telomere-like reads counter-inituively cannot be used to infer telomere length, as they mostly originate from the subtelomere, but can be used as a biomarker for chromatin condensation. Using long-read sequencing, we further show that modern hyperactive Tn5 does not duplicate 9 bp of its target sequence, contrary to common belief. We provide a new tool, Telomemore, which can quantify nonaligning subtelomeric reads. By analyzing several public datasets and generating new multiome fibroblast and B-cell atlases, we show how this new readout can aid single-cell data interpretation. We show how drivers of condensation processes can be inferred, and how it complements common RNA-seq-based cell cycle inference, which fails for monocytes. Telomemore-based analysis of the condensation state is thus a valuable complement to the single-cell analysis toolbox.

Listeria monocytogenes and Salmonella enterica Typhimurium are leading causes of foodborne illness in the United States and frequently implicated in produce outbreaks. Conventional decontamination methods, such as cold-water washes with chlorine, have limited antibacterial efficacy and environmental sustainability. Power ultrasound has emerged as a promising non-thermal alternative, but the molecular mechanisms remain insufficiently elucidated. This study evaluated transcriptomic responses of L. monocytogenes and S. enterica Typhimurium to (i) ultrasound (20 kHz), (ii) chlorine (50 ppm), and (iii) combined ultrasound + chlorine treatments. RNA-seq analysis identified differentially expressed genes, as well as enriched Gene Ontology and KEGG terms. Results showed that ultrasound and chlorine triggered distinct transcriptomic responses. L. monocytogenes exhibited broad transcriptional shifts under ultrasound, including significant upregulation of phosphotransferase system components and central metabolism. Chlorine alone induced a narrower response, with fewer differentially expressed genes clustering into limited functional categories. In contrast, the combined ultrasound + chlorine treatment elicited the strongest response in S. enterica Typhimurium, with enrichment of multiple energy- and metabolism-related pathways, including the citrate cycle, carbon metabolism, and microbial metabolism in diverse environments. These findings provide new insights into ultrasound-triggered responses in foodborne pathogens and may inform development of optimized ultrasound-based hurdle sanitization strategies for produce safety.

This paper presents the annotated genomes of Salmonella Haifa, Salmonella Bangkok, and Salmonella Reading, which are uncommonly isolated from meat in South Africa. Despite their rarity in South Africa, these serotypes have been linked to several high-profile outbreaks in other parts of the world.

This study addresses the limited tools available for assessing food safety risks from cytotoxic Bacillus cereus group strains in contaminated food. We quantified the growth, in skim milk broth, of 17 cytotoxic B. cereus strains across 6 phylogenetic groups with various virulence gene profiles. The strains did not grow in HTST milk at 4 or 6°C. At 10°C, 15 strains exhibited growth; at 8°C, one strain grew; and all strains grew at temperatures ≥14°C. Using growth data from 16 strains, we developed linear secondary growth models and an exposure assessment model. This model, simulating a 5-stage HTST milk supply chain and up to 35 d of consumer storage with an initial contamination of 100 cfu/mL, estimated that 2.81 ± 0.66% and 4.13 ± 2.53% of milk containers would surpass 105 cfu/mL of B. cereus by d 21 and 35, respectively. A sensitivity analysis identified the initial physiological state of cells as the most influential variable affecting predictions for specific isolates. What-if scenarios indicated that increases in mean and variability of consumer storage temperatures significantly affected the predicted B. cereus concentrations in milk. This model serves as an initial tool for risk-based food safety decision-making regarding low-level B. cereus contamination.

The Bacillus cereus sensu stricto (s.s.) species comprises strains of biovar Thuringiensis (Bt) known for their bioinsecticidal activity, as well as strains with foodborne pathogenic potential. Bt strains are identified (i) based on the production of insecticidal crystal proteins, also known as Bt toxins, or (ii) based on the presence of cry, cyt, and vip genes, which encode Bt toxins. Multiple bioinformatics tools have been developed for the detection of crystal protein-encoding genes based on whole-genome sequencing (WGS) data. However, the performance of these tools is yet to be evaluated using phenotypic data. Thus, the goal of this study was to assess the performance of four bioinformatics tools for the detection of crystal protein-encoding genes. The accuracy of sequence-based identification of Bt was determined in reference to phenotypic microscope-based screening for the production of crystal proteins. A total of 58 diverse B. cereus sensu lato strains isolated from clinical, food, environmental, and commercial biopesticide products underwent WGS. Isolates were examined for crystal protein production using phase contrast microscopy. Crystal protein-encoding genes were detected using BtToxin_Digger, BTyper3, IDOPS (identification of pesticidal sequences), and Cry_processor. Out of 58 isolates, the phenotypic production of crystal proteins was confirmed for 18 isolates. Specificity and sensitivity of Bt identification based on sequences were 0.85 and 0.94 for BtToxin_Digger, 0.97 and 0.89 for BTyper3, 0.95 and 0.94 for IDOPS, and 0.88 and 1.00 for Cry_processor, respectively. Cry_processor predicted crystal protein production with the highest specificity, and BtToxin_Digger and IDOPS predicted crystal protein production with the highest sensitivity. Three out of four tested bioinformatics tools performed well overall, with IDOPS achieving high sensitivity and specificity (>0.90).IMPORTANCEStrains of Bacillus cereus sensu stricto (s.s.) biovar Thuringiensis (Bt) are used as organic biopesticides. Bt is differentiated from the foodborne pathogen Bacillus cereus s.s. by the production of insecticidal crystal proteins. Thus, reliable genomic identification of biovar Thuringiensis is necessary to ensure food safety and facilitate risk assessment. This study assessed the accuracy of whole-genome sequencing (WGS)-based identification of Bt compared to phenotypic microscopy-based screening for crystal protein production. Multiple bioinformatics tools were compared to assess their performance in predicting crystal protein production. Among them, identification of pesticidal sequences performed best overall at WGS-based Bt identification.

Salmonella enterica subsp. enterica serotype Typhimurium definitive type 104 (DT104) can infect both humans and animals and is often multidrug-resistant (MDR). Previous studies have indicated that, unlike most S. Typhimurium, the overwhelming majority of DT104 strains produce pertussis-like toxin ArtAB via prophage-encoded genes artAB. However, DT104 that lack artAB have been described on occasion. Here, we identify an MDR DT104 complex lineage circulating among humans and cattle in the USA, which lacks artAB (i.e. the ‘U.S. artAB-negative major clade’; n=42 genomes). Unlike most other bovine- and human-associated DT104 complex strains from the USA (n=230 total genomes), which harbour artAB on prophage Gifsy-1 (n=177), members of the U.S. artAB-negative major clade lack Gifsy-1, as well as anti-inflammatory effector gogB. The U.S. artAB-negative major clade encompasses human- and cattle-associated strains isolated from ≥11 USA states over a 20-year period. The clade was predicted to have lost artAB, Gifsy-1 and gogB circa 1985–1987 (95 % highest posterior density interval 1979.0–1992.1). When compared to DT104 genomes from other regions of the world (n=752 total genomes), several additional, sporadic artAB, Gifsy-1 and/or gogB loss events among clades encompassing five or fewer genomes were observed. Using phenotypic assays that simulate conditions encountered during human and/or bovine digestion, members of the U.S. artAB-negative major clade did not differ from closely related Gifsy-1/artAB/gogB-harbouring U.S. DT104 complex strains (ANOVA raw P>0.05); thus, future research is needed to elucidate the roles that artAB, gogB and Gifsy-1 play in DT104 virulence in humans and animals.

Introduction: Macrococcus species have been isolated from a range of mammals and mammal-derived food products. While they are largely considered to be animal commensals, Macrococcus spp. can be opportunistic pathogens in both veterinary and human clinical settings. This study aimed to provide insight into the evolution, population structure, and functional potential of the Macrococcus genus, with an emphasis on antimicrobial resistance (AMR) and virulence potential.

Methods: All high-quality, publicly available Macrococcus genomes (n = 104, accessed 27 August 2022), plus six South African genomes sequenced here (two strains from bovine clinical mastitis cases and four strains from beef products), underwent taxonomic assignment (using four different approaches), AMR determinant detection (via AMRFinderPlus), and virulence factor detection (using DIAMOND and the core Virulence Factor Database).

Results: Overall, the 110 Macrococcus genomes were of animal commensal, veterinary clinical, food-associated (including food spoilage), and environmental origins; five genomes (4.5%) originated from human clinical cases. Notably, none of the taxonomic assignment methods produced identical results, highlighting the potential for Macrococcus species misidentifications. The most common predicted antimicrobial classes associated with AMR determinants identified across Macrococcus included macrolides, beta-lactams, and aminoglycosides (n = 81, 61, and 44 of 110 genomes; 73.6, 55.5, and 40.0%, respectively). Genes showing homology to Staphylococcus aureus exoenzyme aureolysin were detected across multiple species (using 90% coverage, n = 40 and 77 genomes harboring aureolysin-like genes at 60 and 40% amino acid [AA] identity, respectively). S. aureus Panton-Valentine leucocidin toxin-associated lukF-PV and lukS-PV homologs were identified in eight M. canis genomes (≥40% AA identity, >85% coverage). Using a method that delineates populations using recent gene flow (PopCOGenT), two species (M. caseolyticus and M. armenti) were composed of multiple within-species populations. Notably, M. armenti was partitioned into two populations, which differed in functional potential (e.g., one harbored beta-lactamase family, type II toxin-antitoxin system, and stress response proteins, while the other possessed a Type VII secretion system; PopCOGenT p < 0.05).

Discussion: Overall, this study leverages all publicly available Macrococcus genomes in addition to newly sequenced genomes from South Africa to identify genomic elements associated with AMR or virulence potential, which can be queried in future experiments.