Foodborne bacterial pathogens: genome-based approaches for enduring and emerging threats in a complex and changing world.

Foodborne illnesses pose a substantial health and economic burden, presenting challenges in prevention due to the diverse microbial hazards that can enter and spread within food systems. Various factors, including natural, political and commercial drivers, influence food production and distribution. The risks of foodborne illness will continue to evolve in step with these drivers and with changes to food systems. For example, climate impacts on water availability for agriculture, changes in food sustainability targets and evolving customer preferences can all have an impact on the ecology of foodborne pathogens and the agrifood niches that can carry microorganisms. Whole-genome and metagenome sequencing, combined with microbial surveillance schemes and insights from the food system, can provide authorities and businesses with transformative information to address risks and implement new food safety interventions across the food chain. In this Review, we describe how genome-based approaches have advanced our understanding of the evolution and spread of enduring bacterial foodborne hazards as well as their role in identifying emerging foodborne hazards. Furthermore, foodborne hazards exist in complex microbial communities across the entire food chain, and consideration of these co-existing organisms is essential to understanding the entire ecology supporting pathogen persistence and transmission in an evolving food system.

Global transmission of extended-spectrum cephalosporin resistance in Escherichia coli driven by epidemic plasmids.

BACKGROUND: Extended-spectrum cephalosporins (ESCs) are third and fourth generation cephalosporin antimicrobials used in humans and animals to treat infections due to multidrug-resistant (MDR) bacteria. Resistance to ESCs (ESC-R) in Enterobacterales is predominantly due to the production of extended-spectrum ß-lactamases (ESBLs) and plasmid-mediated AmpC ß-lactamases (AmpCs). The dynamics of ESBLs and AmpCs are changing across countries and host species, the result of global transmission of ESC-R genes. Plasmids are known to play a key role in this dissemination, but the relative importance of different types of plasmids is not fully understood. METHODS: In this study, Escherichia coli with the major ESC-R genes blaCTX-M-1, blaCTX-M-15, blaCTX-M-14 (ESBLs) and blaCMY-2 (AmpC), were selected from diverse host species and other sources across Canada, France and Germany, collected between 2003 and 2017. To examine in detail the vehicles of transmission of the ESC-R genes, long- and short-read sequences were generated to obtain complete contiguous chromosome and plasmid sequences (n = 192 ESC-R E. coli). The types, gene composition and genetic relatedness of these plasmids were investigated, along with association with isolate year, source and geographical origin, and put in context with publicly available plasmid sequences. FINDINGS: We identified five epidemic resistance plasmid subtypes with distinct genetic properties that are associated with the global dissemination of ESC-R genes across multiple E. coli lineages and host species. The IncI1 pST3 blaCTX-M-1 plasmid subtype was found in more diverse sources than the other main plasmid subtypes, whereas IncI1 pST12 blaCMY-2 was more frequent in Canadian and German human and chicken isolates. Clonal expansion also contributed to the dissemination of the IncI1 pST12 blaCMY-2 plasmid in ST131 and ST117 E. coli harbouring this plasmid. The IncI1 pST2 blaCMY-2 subtype was predominant in isolates from humans in France, while the IncF F31:A4:B1 blaCTX-M-15 and F2:A-:B- blaCTX-M-14 plasmid subtypes were frequent in human and cattle isolates across multiple countries. Beyond their epidemic nature with respect to ESC-R genes, in our collection almost all IncI1 pST3 blaCTX-M-1 and IncF F31:A4:B1 blaCTX-M-15 epidemic plasmids also carried multiple antimicrobial resistance (AMR) genes conferring resistance to other antimicrobial classes. Finally, we found genetic signatures in the regions surrounding specific ESC-R genes, identifying the predominant mechanisms of ESC-R gene movement, and using publicly available databases, we identified these epidemic plasmids from widespread bacterial species, host species, countries and continents. INTERPRETATION: We provide evidence that epidemic resistance plasmid subtypes contribute to the global dissemination of ESC-R genes, and in addition, some of these epidemic plasmids confer resistance to multiple other antimicrobial classes. The success of these plasmids suggests that they may have a fitness advantage over other plasmid types and subtypes. Identification and understanding of the vehicles of AMR transmission are crucial to develop and target strategies and interventions to reduce the spread of AMR. FUNDING: This project was supported by the Joint Programming Initiative on Antimicrobial Resistance (JPIAMR), through the Medical Research Council (MRC, MR/R000948/1), the Canadian Institutes of Health Research (CFC-150770), and the Genomics Research and Development Initiative (Government of Canada), the German Federal Ministry of Education and Research (BMBF) grant no. 01KI1709, the French Agency for food environmental and occupational health & safety (Anses), and the French National Reference Center (CNR) for antimicrobial resistance. Support was also provided by the Biotechnology and Biological Sciences Research Council (BBSRC) through the BBSRC Institute Strategic Programme Microbes in the Food ChainBB/R012504/1 and its constituent project BBS/E/F/000PR10348 (Theme 1, Epidemiology and Evolution of Pathogens in the Food Chain).

Scalable neighbour search and alignment with uvaia.

Despite millions of SARS-CoV-2 genomes being sequenced and shared globally, manipulating such data sets is still challenging, especially selecting sequences for focused phylogenetic analysis. We present a novel method, uvaia, which is based on partial and exact sequence similarity for quickly extracting database sequences similar to query sequences of interest. Many SARS-CoV-2 phylogenetic analyses rely on very low numbers of ambiguous sites as a measure of quality since ambiguous sites do not contribute to single nucleotide polymorphism (SNP) differences. Uvaia overcomes this limitation by using measures of sequence similarity which consider partially ambiguous sites, allowing for more ambiguous sequences to be included in the analysis if needed. Such fine-grained definition of similarity allows not only for better phylogenetic analyses, but could also lead to improved classification and biogeographical inferences. Uvaia works natively with compressed files, can use multiple cores and efficiently utilises memory, being able to analyse large data sets on a standard desktop.

A One Health Perspective on Salmonella enterica Serovar Infantis, an Emerging Human Multidrug-Resistant Pathogen.

Salmonella enterica serovar Infantis presents an ever-increasing threat to public health because of its spread throughout many countries and association with high levels of antimicrobial resistance (AMR). We analyzed whole-genome sequences of 5,284 Salmonella Infantis strains from 74 countries, isolated during 1989-2020 from a wide variety of human, animal, and food sources, to compare genetic phylogeny, AMR determinants, and plasmid presence. The global Salmonella Infantis population structure diverged into 3 clusters: a North American cluster, a European cluster, and a global cluster. The levels of AMR varied by Salmonella Infantis cluster and by isolation source; 73% of poultry isolates were multidrug resistant, compared with 35% of human isolates. This finding correlated with the presence of the pESI megaplasmid; 71% of poultry isolates contained pESI, compared with 32% of human isolates. This study provides key information for public health teams engaged in reducing the spread of this pathogen.

Yersinia enterocolitica biovar 1A: An underappreciated potential pathogen in the food chain.

Yersinia enterocolitica is an underreported cause of foodborne gastroenteritis. Little is known of the diversity of Y. enterocolitica isolated from food and which food commodities contribute to human disease. In this study, Y. enterocolitica was isolated from 37/50 raw chicken, 8/10 pork, 8/10 salmon and 1/10 leafy green samples collected at retail in the UK. Up to 10 presumptive Y. enterocolitica isolates per positive sample underwent whole genome sequencing (WGS) and were compared with publicly available genomes. In total, 207 Y. enterocolitica isolates were analyzed and belonged to 38 sequence types (STs). Up to five STs of Y. enterocolitica were isolated from individual food samples and isolates belonging to the same sample and ST differed by 0-74 single nucleotide polymorphisms (SNPs). Biotype was predicted for 205 (99 %) genomes that all belonged to biotype 1A, previously described as non-pathogenic. However, around half (51 %) of food samples contained isolates belonging to the same ST as previously isolated from UK human cases. The closest human-derived isolates shared between 17 and 7978 single nucleotide polymorphisms (SNPs) with the food isolates. Extensive food surveillance is required to determine what food sources are responsible for Y. enterocolitica infections and to re-examine the role of biotype 1A as a human pathogen.

Genomic characterization of Pseudomonas spp. on food: implications for spoilage, antimicrobial resistance and human infection.

BACKGROUND: Pseudomonas species are common on food, but their contribution to the antimicrobial resistance gene (ARG) burden within food or as a source of clinical infection is unknown. Pseudomonas aeruginosa is an opportunistic pathogen responsible for a wide range of infections and is often hard to treat due to intrinsic and acquired ARGs commonly carried by this species. This study aimed to understand the potential role of Pseudomonas on food as a reservoir of ARGs and to assess the presence of potentially clinically significant Pseudomonas aeruginosa strains on food. To achieve this, we assessed the genetic relatedness (using whole genome sequencing) and virulence of food-derived isolates to those collected from humans. RESULTS: A non-specific culturing approach for Pseudomonas recovered the bacterial genus from 28 of 32 (87.5%) retail food samples, although no P. aeruginosa was identified. The Pseudomonas species recovered were not clinically relevant, contained no ARGs and are likely associated with food spoilage. A specific culture method for P. aeruginosa resulted in the recovery of P. aeruginosa from 14 of 128 (11%) retail food samples; isolates contained between four and seven ARGs each and belonged to 16 sequence types (STs), four of which have been isolated from human infections. Food P. aeruginosa isolates from these STs demonstrated high similarity to human-derived isolates, differing by 41-312 single nucleotide polymorphisms (SNPs). There were diverse P. aeruginosa collected from the same food sample with distinct STs present on some samples and isolates belonging to the same ST differing by 19-67 SNPs. The Galleria mellonella infection model showed that 15 of 16 STs isolated from food displayed virulence between a low-virulence (PAO1) and a high virulence (PA14) control. CONCLUSION: The most frequent Pseudomonas recovered from food examined in this study carried no ARGs and are more likely to play a role in food spoilage rather than infection. P. aeruginosa isolates likely to be able to cause human infections and with multidrug resistant genotypes are present on a relatively small but still substantial proportions of retail foods examined. Given the frequency of exposure, the potential contribution of food to the burden of P. aeruginosa infections in humans should be evaluated more closely.

Multidrug resistance plasmids underlie clonal expansions and international spread of Salmonella enterica serotype 1,4,[5],12:i:- ST34 in Southeast Asia.

Salmonella enterica serotype 1,4,[5],12:i:- (Typhimurium monophasic variant) of sequence type (ST) 34 has emerged as the predominant pandemic genotype in recent decades. Despite increasing reports of resistance to antimicrobials in Southeast Asia, Salmonella ST34 population structure and evolution remained understudied in the region. Here we performed detailed genomic investigations on 454 ST34 genomes collected from Vietnam and diverse geographical sources to elucidate the pathogen's epidemiology, evolution and antimicrobial resistance. We showed that ST34 has been introduced into Vietnam in at least nine occasions since 2000, forming five co-circulating major clones responsible for paediatric diarrhoea and bloodstream infection. Most expansion events were associated with acquisitions of large multidrug resistance plasmids of IncHI2 or IncA/C2. Particularly, the self-conjugative IncA/C2 pST34VN2 (co-transferring blaCTX-M-55, mcr-3.1, and qnrS1) underlies local expansion and intercontinental spread in two separate ST34 clones. At the global scale, Southeast Asia was identified as a potential hub for the emergence and dissemination of multidrug resistant Salmonella ST34, and mutation analysis suggests of selection in antimicrobial responses and key virulence factors.

Genomic diversity and epidemiological significance of non-typhoidal Salmonella found in retail food collected in Norfolk, UK.

Non-typhoidal Salmonella (NTS) is a major cause of bacterial gastroenteritis. Although many countries have implemented whole genome sequencing (WGS) of NTS, there is limited knowledge on NTS diversity on food and its contribution to human disease. In this study, the aim was to characterise the NTS genomes from retail foods in a particular region of the UK and assess the contribution to human NTS infections. Raw food samples were collected at retail in a repeated cross-sectional design in Norfolk, UK, including chicken (n=311), leafy green (n=311), pork (n=311), prawn (n=279) and salmon (n=157) samples. Up to eight presumptive NTS isolates per positive sample underwent WGS and were compared to publicly available NTS genomes from UK human cases. NTS was isolated from chicken (9.6 %), prawn (2.9 %) and pork (1.3 %) samples and included 14 serovars, of which Salmonella Infantis and Salmonella Enteritidis were the most common. The S. Enteritidis isolates were only isolated from imported chicken. No antimicrobial resistance determinants were found in prawn isolates, whilst 5.1 % of chicken and 0.64 % of pork samples contained multi-drug resistant NTS. The maximum number of pairwise core non-recombinant single nucleotide polymorphisms (SNPs) amongst isolates from the same sample was used to measure diversity and most samples had a median of two SNPs (range: 0-251). NTS isolates that were within five SNPs to clinical UK isolates belonged to specific serovars: S. Enteritidis and S. Infantis (chicken), and S. I 4,[5],12:i- (pork and chicken). Most NTS isolates that were closely related to human-derived isolates were obtained from imported chicken, but further epidemiological data are required to assess definitively the probable source of the human cases. Continued WGS surveillance of Salmonella on retail food involving multiple isolates from each sample is necessary to capture the diversity of Salmonella and determine the relative importance of different sources of human disease.

Factors associated with plasmid antibiotic resistance gene carriage revealed using large-scale multivariable analysis.

Plasmids are major vectors of bacterial antibiotic resistance, but understanding of factors associated with plasmid antibiotic resistance gene (ARG) carriage is limited. We curated > 14,000 publicly available plasmid genomes and associated metadata. Duplicate and replicate plasmids were excluded; where possible, sample metadata was validated externally (BacDive database). Using Generalised Additive Models (GAMs) we assessed the influence of 12 biotic/abiotic factors (e.g. plasmid genetic factors, isolation source, collection date) on ARG carriage, modelled as a binary outcome. Separate GAMs were built for 10 major ARG types. Multivariable analysis indicated that plasmid ARG carriage patterns across time (collection years), isolation sources (human/livestock) and host bacterial taxa were consistent with antibiotic selection pressure as a driver of plasmid-mediated antibiotic resistance. Only 0.42% livestock plasmids carried carbapenem resistance (compared with 12% human plasmids); conversely, tetracycline resistance was enriched in livestock vs human plasmids, reflecting known prescribing practices. Interpreting results using a timeline of ARG type acquisition (determined by literature review) yielded additional novel insights. More recently acquired ARG types (e.g. colistin and carbapenem) showed increases in plasmid carriage during the date range analysed (1994-2019), potentially reflecting recent onset of selection pressure; they also co-occurred less commonly with ARGs of other types, and virulence genes. Overall, this suggests that following acquisition, plasmid ARGs tend to accumulate under antibiotic selection pressure and co-associate with other adaptive genes (other ARG types, virulence genes), potentially re-enforcing plasmid ARG carriage through co-selection.

Repeated cross-sectional study identifies differing risk factors associated with microbial contamination in common food products in the United Kingdom.

All foods carry microbes, many of which are harmless, but foods can also carry pathogens and/or microbial indicators of contamination. Limited information exists on the co-occurrence of microbes of food safety concern and the factors associated with their presence. Here, a population-based repeated cross-sectional design was used to determine the prevalence and co-occurrence of Escherichia coli, Klebsiella spp., Salmonella spp. and Vibrio spp. in key food commodities - chicken, pork, prawns, salmon and leafy greens. Prevalence in 1,369 food samples for these four target bacterial genera/species varied, while 25.6% of all samples had at least two of the target bacteria and eight different combinations of bacteria were observed as co-occurrence profiles in raw prawns. Imported frozen chicken was 6.4 times more likely to contain Salmonella than domestic chicken, and imported salmon was 5.5 times more likely to be contaminated with E. coli. Seasonality was significantly associated with E. coli and Klebsiella spp. contamination in leafy greens, with higher detection in summer and autumn. Moreover, the odds of Klebsiella spp. contamination were higher in summer in chicken and pork samples. These results provide insight on the bacterial species present on foods at retail, and identify factors associated with the presence of individual bacteria, which are highly relevant for food safety risk assessments and the design of surveillance programmes.

Determination and quantification of microbial communities and antimicrobial resistance on food through host DNA-depleted metagenomics.

Food products carry bacteria unless specifically sterilised. These bacteria can be pathogenic, commensal or associated with food spoilage, and may also be resistant to antimicrobials. Current methods for detecting bacteria on food rely on culturing for specific bacteria, a time-consuming process, or 16S rRNA metabarcoding that can identify different taxa but not their genetic content. Directly sequencing metagenomes of food is inefficient as its own DNA vastly outnumbers the bacterial DNA present. We optimised host DNA depletion enabling efficient sequencing of food microbiota, thereby increasing the proportion of non-host DNA sequenced 13-fold (mean; range: 1.3-40-fold) compared to untreated samples. The method performed best on chicken, pork and leafy green samples which had high mean prokaryotic read proportions post-depletion (0.64, 0.74 and 0.74, respectively), with lower mean prokaryotic read proportions in salmon (0.50) and prawn samples (0.19). We show that bacterial compositions and concentrations of antimicrobial resistance (AMR) genes differed by food type, and that salmon metagenomes were influenced by the production/harvesting method. The approach described in this study is an efficient and effective method of identifying and quantifying the predominant bacteria and AMR genes on food.

Comparative genomics of Campylobacter jejuni from clinical campylobacteriosis stool specimens.

BACKGROUND: Campylobacter jejuni is a pervasive pathogen of major public health concern with a complex ecology requiring accurate and informative approaches to define pathogen diversity during outbreak investigations. Source attribution analysis may be confounded if the genetic diversity of a C. jejuni population is not adequately captured in a single specimen. The aim of this study was to determine the genomic diversity of C. jejuni within individual stool specimens from four campylobacteriosis patients. Direct plating and pre-culture filtration of one stool specimen per patient was used to culture multiple isolates per stool specimen. Whole genome sequencing and pangenome level analysis were used to investigate genomic diversity of C. jejuni within a patient. RESULTS: A total 92 C. jejuni isolates were recovered from four patients presenting with gastroenteritis. The number of isolates ranged from 13 to 30 per patient stool. Three patients yielded a single C. jejuni multilocus sequence type: ST-21 (n = 26, patient 4), ST-61 (n = 30, patient 1) and ST-2066 (n = 23, patient 2). Patient 3 was infected with two different sequence types [ST-51 (n = 12) and ST-354 (n = 1)]. Isolates belonging to the same sequence type from the same patient specimen shared 12-43 core non-recombinant SNPs and 0-20 frameshifts with each other, and the pangenomes of each sequence type consisted of 1406-1491 core genes and 231-264 accessory genes. However, neither the mutation nor the accessory genes were connected to a specific functional gene category. CONCLUSIONS: Our findings show that the C. jejuni population recovered from an individual patient's stool are genetically diverse even within the same ST and may have shared common ancestors before specimens were obtained. The population is unlikely to have evolved from a single isolate at the time point of initial patient infection, leading us to conclude that patients were likely infected with a heterogeneous C. jejuni population. The diversity of the C. jejuni population found within individual stool specimens can inform future methodological approaches to attribution and outbreak investigations.

Dynamics of extended-spectrum cephalosporin resistance genes in Escherichia coli from Europe and North America.

Extended-spectrum cephalosporins (ESCs) are critically important antimicrobial agents for human and veterinary medicine. ESC resistance (ESC-R) genes have spread worldwide through plasmids and clonal expansion, yet the distribution and dynamics of ESC-R genes in different ecological compartments are poorly understood. Here we use whole genome sequence data of Enterobacterales isolates of human and animal origin from Europe and North America and identify contrasting temporal dynamics. AmpC ß-lactamases were initially more dominant in North America in humans and farm animals, only later emerging in Europe. In contrast, specific extended-spectrum ß-lactamases (ESBLs) were initially common in animals from Europe and later emerged in North America. This study identifies differences in the relative importance of plasmids and clonal expansion across different compartments for the spread of different ESC-R genes. Understanding the mechanisms of transmission will be critical in the design of interventions to reduce the spread of antimicrobial resistance.

Genetic characterization of Salmonella Infantis from South Africa, 2004-2016.

Salmonella Infantis is presenting an increasing risk to public health. Of particular concern are the reports of pESI, a multidrug resistance (MDR) encoding megaplasmid, in isolates from multiple countries, but little is known about its presence or diversity in South Africa. Whole genome sequences of 387 S. Infantis isolates from South Africa (2004-2020) were analysed for genetic phylogeny, recombination frequency, antimicrobial resistance (AMR) determinants, plasmid presence and overall gene content. The population structure of South African S. Infantis was substantially different to S. Infantis reported elsewhere; only two thirds of isolates belonged to eBG31, while the remainder were identified as eBG297, a much rarer group globally. Significantly higher levels of recombination were observed in the eBG297 isolates, which was associated with the presence of prophages. The majority of isolates were putatively susceptible to antimicrobials (335/387) and lacked any plasmids (311/387); the megaplasmid pESI was present in just one isolate. A larger proportion of eBG31 isolates, 19% (49/263), contained at least one AMR determinant, compared to eBG297 at 2% (3/124). Comparison of the pan-genomes of isolates from either eBG identified 943 genes significantly associated with eBG, with 43 found exclusively in eBG31 isolates and 34 in eBG297 isolates. This, along with the single nucleotide polymorphism distance and difference in resistance profiles, suggests that eBG31 and eBG297 isolates occupy different niches within South Africa. If antibiotic-resistant S. Infantis emerges in South Africa, probably through the spread of the pESI plasmid, treatment of this infection would be compromised.

Contrasting long-term dynamics of antimicrobial resistance and virulence plasmids in Salmonella Typhimurium from animals.

Plasmids are mobile elements that can carry genes encoding traits of clinical concern, including antimicrobial resistance (AMR) and virulence. Population-level studies of Enterobacterales, including Escherichia coli, Shigella and Klebsiella, indicate that plasmids are important drivers of lineage expansions and dissemination of AMR genes. Salmonella Typhimurium is the second most common cause of salmonellosis in humans and livestock in the UK and Europe. The long-term dynamics of plasmids between S. Typhimurium were investigated using isolates collected through national surveillance of animals in England and Wales over a 25-year period. The population structure of S. Typhimurium and its virulence plasmid (where present) were inferred through phylogenetic analyses using whole-genome sequence data for 496 isolates. Antimicrobial resistance genes and plasmid markers were detected in silico. Phenotypic plasmid characterization, using the Kado and Liu method, was used to confirm the number and size of plasmids. The differences in AMR and plasmids between clades were striking, with livestock clades more likely to carry one or more AMR plasmid and be multi-drug-resistant compared to clades associated with wildlife and companion animals. Multiple small non-AMR plasmids were distributed across clades. However, all hybrid AMR-virulence plasmids and most AMR plasmids were highly clade-associated and persisted over decades, with minimal evidence of horizontal transfer between clades. This contrasts with the role of plasmids in the short-term dissemination of AMR between diverse strains in other Enterobacterales in high-antimicrobial-use settings, with implications for predicting plasmid dissemination amongst S. Typhimurium.

Dynamics of Salmonella enterica and antimicrobial resistance in the Brazilian poultry industry and global impacts on public health.

Non-typhoidal Salmonella enterica is a common cause of diarrhoeal disease; in humans, consumption of contaminated poultry meat is believed to be a major source. Brazil is the world's largest exporter of chicken meat globally, and previous studies have indicated the introduction of Salmonella serovars through imported food products from Brazil. Here we provide an in-depth genomic characterisation and evolutionary analysis to investigate the most prevalent serovars and antimicrobial resistance (AMR) in Brazilian chickens and assess the impact to public health of products contaminated with S. enterica imported into the United Kingdom from Brazil. To do so, we examine 183 Salmonella genomes from chickens in Brazil and 357 genomes from humans, domestic poultry and imported Brazilian poultry products isolated in the United Kingdom. S. enterica serovars Heidelberg and Minnesota were the most prevalent serovars in Brazil and in meat products imported from Brazil into the UK. We extended our analysis to include 1,259 publicly available Salmonella Heidelberg and Salmonella Minnesota genomes for context. The Brazil genomes form clades distinct from global isolates, with temporal analysis suggesting emergence of these Salmonella Heidelberg and Salmonella Minnesota clades in the early 2000s, around the time of the 2003 introduction of the Enteritidis vaccine in Brazilian poultry. Analysis showed genomes within the Salmonella Heidelberg and Salmonella Minnesota clades shared resistance to sulphonamides, tetracyclines and beta-lactams conferred by sul2, tetA and blaCMY-2 genes, not widely observed in other co-circulating serovars despite similar selection pressures. The sul2 and tetA genes were concomitantly carried on IncC plasmids, whereas blaCMY-2 was either co-located with the sul2 and tetA genes on IncC plasmids or independently on IncI1 plasmids. Long-term surveillance data collected in the UK showed no increase in the incidence of Salmonella Heidelberg or Salmonella Minnesota in human cases of clinical disease in the UK following the increase of these two serovars in Brazilian poultry. In addition, almost all of the small number of UK-derived genomes which cluster with the Brazilian poultry-derived sequences could either be attributed to human cases with a recent history of foreign travel or were from imported Brazilian food products. These findings indicate that even should Salmonella from imported Brazilian poultry products reach UK consumers, they are very unlikely to be causing disease. No evidence of the Brazilian strains of Salmonella Heidelberg or Salmonella Minnesota were observed in UK domestic chickens. These findings suggest that introduction of the Salmonella Enteritidis vaccine, in addition to increasing antimicrobial use, could have resulted in replacement of salmonellae in Brazilian poultry flocks with serovars that are more drug resistant, but less associated with disease in humans in the UK. The plasmids conferring resistance to beta-lactams, sulphonamides and tetracyclines likely conferred a competitive advantage to the Salmonella Minnesota and Salmonella Heidelberg serovars in this setting of high antimicrobial use, but the apparent lack of transfer to other serovars present in the same setting suggests barriers to horizontal gene transfer that could be exploited in intervention strategies to reduce AMR. The insights obtained reinforce the importance of One Health genomic surveillance.

Mobility of antimicrobial resistance across serovars and disease presentations in non-typhoidal Salmonella from animals and humans in Vietnam.

Non-typhoidal Salmonella (NTS) is a major cause of bacterial enterocolitis globally but also causes invasive bloodstream infections. Antimicrobial resistance (AMR) hampers the treatment of these infections and understanding how AMR spreads between NTS may help in developing effective strategies. We investigated NTS isolates associated with invasive disease, diarrhoeal disease and asymptomatic carriage in animals and humans from Vietnam. Isolates included multiple serovars and both common and rare phenotypic AMR profiles; long- and short-read sequencing was used to investigate the genetic mechanisms and genomic backgrounds associated with phenotypic AMR profiles. We demonstrate concordance between most AMR genotypes and phenotypes but identified large genotypic diversity in clinically relevant phenotypes and the high mobility potential of AMR genes (ARGs) in this setting. We found that 84 % of ARGs identified were located on plasmids, most commonly those containing IncHI1A_1 and IncHI1B(R27)_1_R27 replicons (33%), and those containing IncHI2_1 and IncHI2A_1 replicons (31%). The vast majority (95%) of ARGS were found within 10 kbp of IS6/IS26 elements, which provide plasmids with a mechanism to exchange ARGs between plasmids and other parts of the genome. Whole genome sequencing with targeted long-read sequencing applied in a One Health context identified a comparatively limited number of insertion sequences and plasmid replicons associated with AMR. Therefore, in the context of NTS from Vietnam and likely for other settings as well, the mechanisms by which ARGs move contribute to a more successful AMR profile than the specific ARGs, facilitating the adaptation of bacteria to different environments or selection pressures.

Characterization of Salmonella enterica Contamination in Pork and Poultry Meat from São Paulo/Brazil: Serotypes, Genotypes and Antimicrobial Resistance Profiles.

Salmonellosis is a zoonosis of major relevance to global public health. Here we present the assessment of Salmonella enterica contamination in pork and poultry meat sold at retail markets in São Paulo, Brazil. A total of 780 meat samples (386 poultry meat and 394 pork samples) were collected from 132 markets. From these, 57 samples (7.3%) were positive for S. enterica isolation, including 32 (8.3%) poultry meat and 25 (6.3%) pork samples. S. enterica isolates were further characterized for serotyping, antimicrobial resistance and genotyping by amplified fragment length polymorphism and pulsed field gel electrophoresis. Antimicrobial resistance analysis demonstrated two main profiles: pork isolates were more resistant to macrolides, ß-lactams, tetracycline, phenicols, and fluoroquinolones, and poultry meat isolates presented higher resistance to fluoroquinolones, sulfonamides, tetracycline, and ß-lactams. A total of 72.4% of poultry meat isolates were identified as S. Heidelberg, while most of pork isolates were S. Typhimurium (31.7%) and S. Give (16.7%). Genotyping resulted in most clusters consisting exclusively of pork or poultry meat, no cross-contamination was detected, and a tendency to differentiate isolates according to their serotypes and markets of origin. High resistance rates to critically important antimicrobials reinforce the importance of controlling Salmonella contamination in meat production chains.

Diversity of blaCTX-M-1-carrying plasmids recovered from Escherichia coli isolated from Canadian domestic animals.

Conserved IncI1 and IncHI1 plasmids carrying blaCTX-M-1 have been found circulating in chickens and horses from continental Europe, respectively. In Canada, blaCTX-M-1 is overwhelmingly the most common blaCTX-M variant found in Escherichia coli from chicken and horses and can be recovered at lower frequencies in swine, cattle, and dogs. Whole-genome sequencing has identified a large genetic diversity of isolates carrying this variant, warranting further investigations into the plasmids carrying this gene. Therefore, the objective of this study was to describe the genetic profiles of blaCTX-M-1 plasmids circulating in E. coli from Canadian domestic animals and compare them to those recovered in animals in Europe. Fifty-one blaCTX-M-1 positive E. coli isolates from chicken (n = 14), horses (racetrack horses n = 11; community horses n = 3), swine (n = 7), turkey (n = 6), dogs (n = 5), beef cattle (n = 3), and dairy cattle (n = 2) were selected for plasmid characterization. Sequences were obtained through both Illumina and Oxford Nanopore technologies. Genomes were assembled using either Unicycler hybrid assembly or Flye with polishing performed using Pilon. blaCTX-M-1 was found residing on a plasmid in 45 isolates and chromosomally located in six isolates. A conserved IncI1/ST3 plasmid was identified among chicken (n = 12), turkey (n = 4), swine (n = 6), dog (n = 2), and beef cattle (n = 2) isolates. When compared against publicly available data, these plasmids showed a high degree of similarity to those identified in isolates from poultry and swine in Europe. These results suggest that an epidemic IncI1/ST3 plasmid similar to the one found in Europe is contributing to the spread of blaCTX-M-1 in Canada. A conserved IncHI1/FIA(HI1)/ST2 plasmid was also recovered from nearly all racetrack horse isolates (n = 10). Although IncHI1/ST2 plasmids have been reported among European horse isolates, IncHI1/ST9 plasmids appear to be more widespread. Further studies are necessary to understand the factors contributing to these plasmids' success in their respective populations.

Tatajuba: exploring the distribution of homopolymer tracts.

Length variation of homopolymeric tracts, which induces phase variation, is known to regulate gene expression leading to phenotypic variation in a wide range of bacterial species. There is no specialized bioinformatics software which can, at scale, exhaustively explore and describe these features from sequencing data. Identifying these is non-trivial as sequencing and bioinformatics methods are prone to introducing artefacts when presented with homopolymeric tracts due to the decreased base diversity. We present tatajuba, which can automatically identify potential homopolymeric tracts and help predict their putative phenotypic impact, allowing for rapid investigation. We use it to detect all tracts in two separate datasets, one of Campylobacter jejuni and one of three Bordetella species, and to highlight those tracts that are polymorphic across samples. With this we confirm homopolymer tract variation with phenotypic impact found in previous studies and additionally find many more with potential variability. The software is written in C and is available under the open source licence GNU GPLv3.