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

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.

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.

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.

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.

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.

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.

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.

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

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. Here we determine that S. enterica serovars Heidelberg and Minnesota were the most prevalent serovars in Brazilian poultry and in poultry products imported from Brazil into the UK. However, long-term surveillance data collected by the UK Health Security Agency 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. Salmonella Heidelberg and Salmonella Minnesota clades shared resistance to sulphonamides, tetracyclines and beta-lactams conferred by sul2, tetA and bla_CMY-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 bla_CMY-2 was either co-located with the sul2 and tetA genes on IncC plasmids or independently on IncI1 plasmids. 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.

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.

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.

A roadmap for the generation of benchmarking resources for antimicrobial resistance detection using next generation sequencing

Next Generation Sequencing technologies significantly impact the field of Antimicrobial Resistance (AMR) detection and monitoring, with immediate uses in diagnosis and risk assessment. For this application and in general, considerable challenges remain in demonstrating sufficient trust to act upon the meaningful information produced from raw data, partly because of the reliance on bioinformatics pipelines, which can produce different results and therefore lead to different interpretations. With the constant evolution of the field, it is difficult to identify, harmonise and recommend specific methods for large-scale implementations over time. In this article, we propose to address this challenge through establishing a transparent, performance-based, evaluation approach to provide flexibility in the bioinformatics tools of choice, while demonstrating proficiency in meeting common performance standards. The approach is two-fold: first, a community-driven effort to establish and maintain “live” (dynamic) benchmarking platforms to provide relevant performance metrics, based on different use-cases, that would evolve together with the AMR field; second, agreed and defined datasets to allow the pipelines’ implementation, validation, and quality-control over time. Following previous discussions on the main challenges linked to this approach, we provide concrete recommendations and future steps, related to different aspects of the design of benchmarks, such as the selection and the characteristics of the datasets (quality, choice of pathogens and resistances, etc.), the evaluation criteria of the pipelines, and the way these resources should be deployed in the community.

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.

Genomic adaptations of Campylobacter jejuni to long-term human colonization

BACKGROUND

Campylobacter is a genus of bacteria that has been isolated from the gastrointestinal tract of humans and animals, and the environments they inhabit around the world. Campylobacter adapt to new environments by changes in their gene content and expression, but little is known about how they adapt to long-term human colonization. In this study, the genomes of 31 isolates from a New Zealand patient and 22 isolates from a United Kingdom patient belonging to Campylobacter jejuni sequence type 45 (ST45) were compared with 209 ST45 genomes from other sources to identify the mechanisms by which Campylobacter adapts to long-term human colonization. In addition, the New Zealand patient had their microbiota investigated using 16S rRNA metabarcoding, and their level of inflammation and immunosuppression analyzed using biochemical tests, to determine how Campylobacter adapts to a changing gastrointestinal tract.

RESULTS

There was some evidence that long-term colonization led to genome degradation, but more evidence that Campylobacter adapted through the accumulation of non-synonymous single nucleotide polymorphisms (SNPs) and frameshifts in genes involved in cell motility, signal transduction and the major outer membrane protein (MOMP). The New Zealand patient also displayed considerable variation in their microbiome, inflammation and immunosuppression over five months, and the Campylobacter collected from this patient could be divided into two subpopulations, the proportion of which correlated with the amount of gastrointestinal inflammation.

CONCLUSIONS

This study demonstrates how genomics, phylogenetics, 16S rRNA metabarcoding and biochemical markers can provide insight into how Campylobacter adapts to changing environments within human hosts. This study also demonstrates that long-term human colonization selects for changes in Campylobacter genes involved in cell motility, signal transduction and the MOMP; and that genetically distinct subpopulations of Campylobacter evolve to adapt to the changing gastrointestinal environment.

Characterization of a pESI-like plasmid and analysis of multidrug-resistant Salmonella enterica Infantis isolates in England and Wales

Salmonella enterica serovar Infantis is the fifth most common Salmonella serovar isolated in England and Wales. Epidemiological, genotyping and antimicrobial-resistance data for S. enterica Infantis isolates were used to analyse English and Welsh demographics over a 5 year period. Travel cases associated with S. enterica Infantis were mainly from Asia, followed by cases from Europe and North America. Since 2000, increasing numbers of S. enterica Infantis had multidrug resistance determinants harboured on a large plasmid termed ‘plasmid of emerging S. enterica Infantis’ (pESI). Between 2013 and 2018, 42 S. enterica Infantis isolates were isolated from humans and food that harboured resistance determinants to multiple antimicrobial classes present on a pESI-like plasmid, including extended-spectrum β-lactamases (ESBLs; bla_CTX-M-65). Nanopore sequencing of an ESBL-producing human S. enterica Infantis isolate indicated the presence of two regions on an IncFIB pESI-like plasmid harbouring multiple resistance genes. Phylogenetic analysis of the English and Welsh S. enterica Infantis population indicated that the majority of multidrug-resistant isolates harbouring the pESI-like plasmid belonged to a single clade maintained within the population. The bla_CTX-M-65 ESBL isolates first isolated in 2013 comprise a lineage within this clade, which was mainly associated with South America. Our data, therefore, show the emergence of a stable resistant clone that has been in circulation for some time in the human population in England and Wales, highlighting the necessity of monitoring resistance in this serovar.

Whole genome sequencing reveals great diversity of Vibrio spp in prawns at retail

Consumption of prawns as a protein source has been on the rise worldwide with seafood identified as the predominant attributable source of human vibriosis. However, surveillance of non-cholera Vibrio is limited both in public health and in food. Using a population- and market share-weighted study design, 211 prawn samples were collected and cultured for Vibrio spp. Contamination was detected in 46 % of samples, and multiple diverse Vibrio isolates were obtained from 34 % of positive samples. Whole genome sequencing (WGS) and phylogenetic analysis illustrated a comprehensive view of Vibrio species diversity in prawns available at retail, with no known pathogenicity markers identified in Vibrio parahaemolyticus and V. cholerae. Antimicrobial resistance genes were found in 77 % of isolates, and 12 % carried genes conferring resistance to three or more drug classes. Resistance genes were found predominantly in V. parahaemolyticus, though multiple resistance genes were also identified in V. cholerae and V. vulnificus. This study highlights the large diversity in Vibrio derived from prawns at retail, even within a single sample. Although there was little evidence in this study that prawns are a major source of vibriosis in the UK, surveillance of non-cholera Vibrio is very limited. This study illustrates the value of expanding WGS surveillance efforts of non-cholera Vibrios in the food chain to identify critical control points for food safety through the production system and to determine the full extent of the public health impact.

Transmission of Multidrug-Resistant Salmonella enterica Subspecies enterica 4,[5],12:i:- Sequence Type 34 between Europe and the United States

Multidrug-resistant Salmonellaenterica subspecies enterica 4,[5],12:i:- sequence type 34 represents a worldwide public health risk. To determine its origin in the United States, we reconstructed a time-scaled phylogeny with a discrete trait geospatial model. The clone in the United States was introduced from Europe on multiple occasions in the early 2000s.

Large-scale sequencing of SARS-CoV-2 genomes from one region allows detailed epidemiology and enables local outbreak management

The COVID-19 pandemic has spread rapidly throughout the world. In the UK, the initial peak was in April 2020; in the county of Norfolk (UK) and surrounding areas, which has a stable, low-density population, over 3200 cases were reported between March and August 2020. As part of the activities of the national COVID-19 Genomics Consortium (COG-UK) we undertook whole genome sequencing of the SARS-CoV-2 genomes present in positive clinical samples from the Norfolk region. These samples were collected by four major hospitals, multiple minor hospitals, care facilities and community organizations within Norfolk and surrounding areas. We combined clinical metadata with the sequencing data from regional SARS-CoV-2 genomes to understand the origins, genetic variation, transmission and expansion (spread) of the virus within the region and provide context nationally. Data were fed back into the national effort for pandemic management, whilst simultaneously being used to assist local outbreak analyses. Overall, 1565 positive samples (172 per 100 000 population) from 1376 cases were evaluated; for 140 cases between two and six samples were available providing longitudinal data. This represented 42.6 % of all positive samples identified by hospital testing in the region and encompassed those with clinical need, and health and care workers and their families. In total, 1035 cases had genome sequences of sufficient quality to provide phylogenetic lineages. These genomes belonged to 26 distinct global lineages, indicating that there were multiple separate introductions into the region. Furthermore, 100 genetically distinct UK lineages were detected demonstrating local evolution, at a rate of ~2 SNPs per month, and multiple co-occurring lineages as the pandemic progressed. Our analysis: identified a discrete sublineage associated with six care facilities; found no evidence of reinfection in longitudinal samples; ruled out a nosocomial outbreak; identified 16 lineages in key workers which were not in patients, indicating infection control measures were effective; and found the D614G spike protein mutation which is linked to increased transmissibility dominates the samples and rapidly confirmed relatedness of cases in an outbreak at a food processing facility. The large-scale genome sequencing of SARS-CoV-2-positive samples has provided valuable additional data for public health epidemiology in the Norfolk region, and will continue to help identify and untangle hidden transmission chains as the pandemic evolves.

CoronaHiT: high-throughput sequencing of SARS-CoV-2 genomes

We present CoronaHiT, a platform and throughput flexible method for sequencing SARS-CoV-2 genomes (≤ 96 on MinION or > 96 on Illumina NextSeq) depending on changing requirements experienced during the pandemic. CoronaHiT uses transposase-based library preparation of ARTIC PCR products. Method performance was demonstrated by sequencing 2 plates containing 95 and 59 SARS-CoV-2 genomes on nanopore and Illumina platforms and comparing to the ARTIC LoCost nanopore method. Of the 154 samples sequenced using all 3 methods, ≥ 90% genome coverage was obtained for 64.3% using ARTIC LoCost, 71.4% using CoronaHiT-ONT and 76.6% using CoronaHiT-Illumina, with almost identical clustering on a maximum likelihood tree. This protocol will aid the rapid expansion of SARS-CoV-2 genome sequencing globally.

A roadmap for the generation of benchmarking resources for antimicrobial resistance detection using next generation sequencing

Next Generation Sequencing technologies significantly impact the field of Antimicrobial Resistance (AMR) detection and monitoring, with immediate uses in diagnosis and risk assessment. For this application and in general, considerable challenges remain in demonstrating sufficient trust to act upon the meaningful information produced from raw data, partly because of the reliance on bioinformatics pipelines, which can produce different results and therefore lead to different interpretations. With the constant evolution of the field, it is difficult to identify, harmonise and recommend specific methods for large-scale implementations over time. In this article, we propose to address this challenge through establishing a transparent, performance-based, evaluation approach to provide flexibility in the bioinformatics tools of choice, while demonstrating proficiency in meeting common performance standards. The approach is two-fold: first, a community-driven effort to establish and maintain “live” (dynamic) benchmarking platforms to provide relevant performance metrics, based on different use-cases, that would evolve together with the AMR field; second, agreed and defined datasets to allow the pipelines’ implementation, validation, and quality-control over time. Following previous discussions on the main challenges linked to this approach, we provide concrete recommendations and future steps, related to different aspects of the design of benchmarks, such as the selection and the characteristics of the datasets (quality, choice of pathogens and resistances, etc.), the evaluation criteria of the pipelines, and the way these resources should be deployed in the community.

The role of animals as a source of antimicrobial resistant nontyphoidal Salmonella causing invasive and non-invasive human disease in Vietnam

BACKGROUND

Nontyphoidal Salmonella (NTS) are associated with both diarrhea and bacteremia. Antimicrobial resistance (AMR) is common in NTS in low-middle income countries, but the major source(s) of AMR NTS in humans are not known. Here, we aimed to assess the role of animals as a source of AMR in human NTS infections in Vietnam. We retrospectively combined and analyzed 672 NTS human and animal isolates from four studies in southern Vietnam and compared serovars, sequence types (ST), and AMR profiles. We generated a population structure of circulating organisms and aimed to attribute sources of AMR in NTS causing invasive and noninvasive disease in humans using Bayesian multinomial mixture models.

RESULTS

Among 672 NTS isolates, 148 (22%) originated from human blood, 211 (31%) from human stool, and 313 (47%) from animal stool. The distribution of serovars, STs, and AMR profiles differed among sources; serovars Enteritidis, Typhimurium, and Weltevreden were the most common in human blood, human stool, and animals, respectively. We identified an association between the source of NTS and AMR profile; the majority of AMR isolates were isolated from human blood (p < 0.001). Modelling by ST-AMR profile found chickens and pigs were likely the major sources of AMR NTS in human blood and stool, respectively; but unsampled sources were found to be a major contributor.

CONCLUSIONS

Antimicrobial use in food animals is hypothesized to play role in the emergence of AMR in human pathogens. Our cross-sectional population-based approach suggests a significant overlap between AMR in NTS in animals and humans, but animal NTS does explain the full extent of AMR in human NTS infections in Vietnam.

Anti-inflammatories in Alzheimer's disease-potential therapy or spurious correlate?

Epidemiological evidence suggests non-steroidal anti-inflammatory drugs reduce the risk of Alzheimer’s disease. However, clinical trials have found no evidence of non-steroidal anti-inflammatory drug efficacy. This incongruence may be due to the wrong non-steroidal anti-inflammatory drugs being tested in robust clinical trials or the epidemiological findings being caused by confounding factors. Therefore, this study used logistic regression and the innovative approach of negative binomial generalized linear mixed modelling to investigate both prevalence and cognitive decline, respectively, in the Alzheimer’s Disease Neuroimaging dataset for each commonly used non-steroidal anti-inflammatory drug and paracetamol. Use of most non-steroidal anti-inflammatories was associated with reduced Alzheimer’s disease prevalence yet no effect on cognitive decline was observed. Paracetamol had a similar effect on prevalence to these non-steroidal anti-inflammatory drugs suggesting this association is independent of the anti-inflammatory effects and that previous results may be due to spurious associations. Interestingly, diclofenac use was significantly associated with both reduce incidence and slower cognitive decline warranting further research into the potential therapeutic effects of diclofenac in Alzheimer’s disease.

Genetic variation associated with infection and the environment in the accidental pathogen Burkholderia pseudomallei

The environmental bacterium Burkholderia pseudomallei causes melioidosis, an important endemic human disease in tropical and sub-tropical countries. This bacterium occupies broad ecological niches including soil, contaminated water, single-cell microbes, plants and infection in a range of animal species. Here, we performed genome-wide association studies for genetic determinants of environmental and human adaptation using a combined dataset of 1,010 whole genome sequences of B. pseudomallei from Northeast Thailand and Australia, representing two major disease hotspots. With these data, we identified 47 genes from 26 distinct loci associated with clinical or environmental isolates from Thailand and replicated 12 genes in an independent Australian cohort. We next outlined the selective pressures on the genetic loci (dN/dS) and the frequency at which they had been gained or lost throughout their evolutionary history, reflecting the bacterial adaptability to a wide range of ecological niches. Finally, we highlighted loci likely implicated in human disease.

Taxonomic resolution of the ribosomal RNA operon in bacteria: implications for its use with long-read sequencing

DNA barcoding through the use of amplified regions of the ribosomal operon, such as the 16S gene, is a routine method to gain an overview of the microbial taxonomic diversity within a sample without the need to isolate and culture the microbes present. However, bacterial cells usually have multiple copies of this ribosomal operon, and choosing the 'wrong' copy could provide a misleading species classification. While this presents less of a problem for well-characterized organisms with large sequence databases to interrogate, it is a significant challenge for lesser known organisms with unknown copy number and diversity. Using the entire length of the ribosomal operon, which encompasses the 16S, 23S, 5S and internal transcribed spacer regions, should provide greater taxonomic resolution but has not been well explored. Here, we use publicly available reference genomes and explore the theoretical boundaries when using concatenated genes and the full-length ribosomal operons, which has been made possible by the development and uptake of long-read sequencing technologies. We quantify the issues of both copy choice and operon length in a phylogenetic context to demonstrate that longer regions improve the phylogenetic signal while maintaining taxonomic accuracy.

Salmonella enterica Serotype 4,[5],12:i:- in Swine in the United States Midwest: An Emerging Multidrug-Resistant Clade

Background

Salmonella 4,[5],12:i:-, a worldwide emerging pathogen that causes many food-borne outbreaks mostly attributed to pig and pig products, is expanding in the United States.

Methods

Whole-genome sequencing was applied to conduct multiple comparisons of 659 S. 4,[5],12:i:- and 325 Salmonella Typhimurium from different sources and locations (ie, the United States and Europe) to assess their genetic heterogeneity, with a focus on strains recovered from swine in the US Midwest. In addition, the presence of resistance genes and other virulence factors was detected and the antimicrobial resistance phenotypes of 50 and 22 isolates of livestock and human origin, respectively, was determined.

Results

The S. 4,5,12:i:- strains formed two main clades regardless of their source and geographic origin. Most (84%) of the US isolates recovered in 2014-2016, including those (48 of 51) recovered from swine in the US Midwest, were part of an emerging clade. In this clade, multiple genotypic resistance determinants were predominant, including resistance against ampicillin, streptomycin, sulfonamides, and tetracyclines. Phenotypic resistance to enrofloxacin (11 of 50) and ceftiofur (9 of 50) was found in conjunction with the presence of plasmid-mediated resistance genes (qnrB19/qnrB2/qnrS1 and blaCMY-2/blaSHV-12, respectively). Higher similarity was also found between S. 4,[5],12:i:- from the emerging clade and S. Typhimurium from Europe than with S. Typhimurium from the United States.

Conclusions

Salmonella 4,[5],12:i:- currently circulating in swine in the US Midwest are likely to be part of an emerging multidrug-resistant clade first reported in Europe, and can carry plasmid-mediated resistance genes that may be transmitted horizontally to other bacteria, and thus may represent a public health concern.

Phylogenomic analysis of gastroenteritis-associated Clostridium perfringens in England and Wales over a 7-year period indicates distribution of clonal toxigenic strains in multiple outbreaks and extensive involvement of enterotoxin-encoding (CPE) plasmids

Clostridium perfringens is a major enteric pathogen known to cause gastroenteritis in human adults. Although major outbreak cases are frequently reported, only limited whole-genome sequencing (WGS) based studies have been performed to understand the genomic epidemiology and virulence gene content of outbreak-associated C. perfringens strains. We performed phylogenomic analysis on 109 C. perfringens isolates (human and food) obtained from disease cases in England and Wales between 2011 and 2017. Initial findings highlighted the enhanced discriminatory power of WGS in profiling outbreak C. perfringens strains, when compared to the current Public Health England referencing laboratory technique of fluorescent amplified fragment length polymorphism analysis. Further analysis identified that isogenic C. perfringens strains were associated with nine distinct care-home-associated outbreaks over the course of a 5-year interval, indicating a potential common source linked to these outbreaks or transmission over time and space. As expected, the enterotoxin cpe gene was encoded in all but 4 isolates (96.3 %; 105/109), with virulence plasmids encoding cpe (particularly pCPF5603 and pCPF4969 plasmids) extensively distributed (82.6 %; 90/109). Genes encoding accessory virulence factors, such as beta-2 toxin, were commonly detected (46.7 %; 51/109), and genes encoding phage proteins were also frequently identified. Overall, this large-scale genomic study of gastroenteritis-associated C. perfringens suggested that three major cpe-encoding (toxinotype F) genotypes underlie these outbreaks: strains carrying (1) pCPF5603 plasmid, (2) pCPF4969 plasmid and (3) chromosomal-cpe strains. Our findings substantially expanded our knowledge on type F C. perfringens involved in human-associated gastroenteritis, with further studies required to fully probe the dissemination and regional reservoirs of this enteric pathogen, which may help devise effective prevention strategies to reduce the food-poisoning disease burden in vulnerable patients, such as the elderly.

Determining antimicrobial susceptibility in Salmonella enterica serovar Typhimurium through whole genome sequencing: a comparison against multiple phenotypic susceptibility testing methods

Background

UK public health organisations perform routine antimicrobial susceptibility tests (ASTs) to characterise the potential for antimicrobial resistance in Salmonella enterica serovars. Genetic determinants of these resistance mechanisms are detectable by whole genome sequencing (WGS), however the viability of WGS-based genotyping as an alternative resistance screening tool remains uncertain. We compared WGS-based genotyping, disk diffusion and agar dilution to the broth microdilution reference AST for 102 Salmonella enterica serovar Typhimurium (S. Typhimurium) isolates across 11 antimicrobial compounds.

Results

Genotyping concordance, interpreted using epidemiological cut-offs (ECOFFs), was 89.8% (1007/1122) with 0.83 sensitivity and 0.96 specificity. For seven antimicrobials interpreted using Salmonella clinical breakpoints, genotyping produced 0.84 sensitivity and 0.88 specificity. Although less accurate than disk diffusion (0.94 sensitivity, 0.93 specificity) and agar dilution (0.83 sensitivity, 0.98 specificity), genotyping performance improved to 0.89 sensitivity and 0.97 specificity when two antimicrobials with relatively high very major error rates were excluded (streptomycin and sulfamethoxazole).

Conclusions

An 89.8% concordance from WGS-based AST predictions using ECOFF interpretations suggest that WGS would serve as an effective screening tool for the tracking of antimicrobial resistance mechanisms in S. Typhimurium. For use as a standalone clinical diagnostic screen, further work is required to reduce the error rates for specific antimicrobials.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1520-9) contains supplementary material, which is available to authorized users.

Genomic analysis of Klebsiella pneumoniae isolates from Malawi reveals acquisition of multiple ESBL determinants across diverse lineages

OBJECTIVES

ESBL-producing Klebsiella pneumoniae (KPN) pose a major threat to human health globally. We carried out a WGS study to understand the genetic background of ESBL-producing KPN in Malawi and place them in the context of other global isolates.

METHODS

We sequenced genomes of 72 invasive and carriage KPN isolates collected from patients admitted to Queen Elizabeth Central Hospital, Blantyre, Malawi. We performed phylogenetic and population structure analyses on these and previously published genomes from Kenya (n = 66) and from outside sub-Saharan Africa (n = 67). We screened for presence of antimicrobial resistance (AMR) genetic determinants and carried out association analyses by genomic sequence cluster, AMR phenotype and time.

RESULTS

Malawian isolates fit within the global population structure of KPN, clustering into the major lineages of KpI, KpII and KpIII. KpI isolates from Malawi were more related to those from Kenya, with both collections exhibiting more clonality than isolates from the rest of the world. We identified multiple ESBL genes, including blaCTX-M-15, several blaSHV, blaTEM-63 and blaOXA-10, and other AMR genes, across diverse lineages of the KPN isolates from Malawi. No carbapenem resistance genes were detected; however, we detected IncFII and IncFIB plasmids that were similar to the carbapenem resistance-associated plasmid pNDM-mar.

CONCLUSIONS

There are multiple ESBL genes across diverse KPN lineages in Malawi and plasmids in circulation that are capable of carrying carbapenem resistance. Unless appropriate interventions are rapidly put in place, these may lead to a high burden of locally untreatable infection in vulnerable populations.

Antimicrobial Resistance Diversity Suggestive of Distinct Salmonella Typhimurium Sources or Selective Pressures in Food-Production Animals

Salmonella enterica subsp. enterica serovar Typhimurium is a common cause of enterocolitis in humans globally, with multidrug resistant (MDR) strains posing an enhanced threat. S. Typhimurium is also a pathogen in food-production animals, and these populations can act as reservoirs of the bacterium. Therefore, surveillance and control measures within food-production animal populations are of importance both to animal and human health and have the potential to be enhanced though improved understanding of the epidemiology of S. Typhimurium within and between food-production animal populations. Here, data from Scotland and national surveillance England and Wales data for isolates from cattle (n = 1115), chickens (n = 248) and pigs (n = 2174) collected between 2003 and 2014 were analyzed. Ecological diversity analyses and rarefaction curves were used to compare the diversity of observed antimicrobial resistance (AMR) profiles between the host species, and within host species populations. Higher AMR profile diversity was observed in isolates from pigs compared to chickens across diversity measures and isolates from cattle for three of four diversity measures. Variation in AMR profile diversity between production sectors was noted, with higher AMR diversity of isolates from broiler compared to layer chickens, breeder compared to rearer and finisher pigs and beef compared to dairy cattle. Findings indicate variation in AMR profile diversity both within and between food-production animal host species. These observations suggest alternate sources of AMR bacteria and/or variation in selective evolutionary pressures within and between food-production animal host species populations.

Genomic Analysis of Salmonella enterica Serovar Typhimurium DT160 Associated with a 14-Year Outbreak, New Zealand, 1998-2012

During 1998-2012, an extended outbreak of Salmonella enterica serovar Typhimurium definitive type 160 (DT160) affected >3,000 humans and killed wild birds in New Zealand. However, the relationship between DT160 within these 2 host groups and the origin of the outbreak are unknown. Whole-genome sequencing was used to compare 109 Salmonella Typhimurium DT160 isolates from sources throughout New Zealand. We provide evidence that DT160 was introduced into New Zealand around 1997 and rapidly propagated throughout the country, becoming more genetically diverse over time. The genetic heterogeneity was evenly distributed across multiple predicted functional protein groups, and we found no evidence of host group differentiation between isolates collected from human, poultry, bovid, and wild bird sources, indicating ongoing transmission between these host groups. Our findings demonstrate how a comparative genomic approach can be used to gain insight into outbreaks, disease transmission, and the evolution of a multihost pathogen after a probable point-source introduction.

Integrated view of Vibrio cholerae in the Americas

Latin America has experienced two of the largest cholera epidemics in modern history; one in 1991 and the other in 2010. However, confusion still surrounds the relationships between globally circulating pandemic Vibrio cholerae clones and local bacterial populations. We used whole-genome sequencing to characterize cholera across the Americas over a 40-year time span. We found that both epidemics were the result of intercontinental introductions of seventh pandemic El Tor V. cholerae and that at least seven lineages local to the Americas are associated with disease that differs epidemiologically from epidemic cholera. Our results consolidate historical accounts of pandemic cholera with data to show the importance of local lineages, presenting an integrated view of cholera that is important to the design of future disease control strategies.

Genomic landscape of extended-spectrum β-lactamase resistance in Escherichia coli from an urban African setting

Objectives

Efforts to treat Escherichia coli infections are increasingly being compromised by the rapid, global spread of antimicrobial resistance (AMR). Whilst AMR in E. coli has been extensively investigated in resource-rich settings, in sub-Saharan Africa molecular patterns of AMR are not well described. In this study, we have begun to explore the population structure and molecular determinants of AMR amongst E. coli isolates from Malawi.

Methods

Ninety-four E. coli isolates from patients admitted to Queen's Hospital, Malawi, were whole-genome sequenced. The isolates were selected on the basis of diversity of phenotypic resistance profiles and clinical source of isolation (blood, CSF and rectal swab). Sequence data were analysed using comparative genomics and phylogenetics.

Results

Our results revealed the presence of five clades, which were strongly associated with E. coli phylogroups A, B1, B2, D and F. We identified 43 multilocus STs, of which ST131 (14.9%) and ST12 (9.6%) were the most common. We identified 25 AMR genes. The most common ESBL gene was bla CTX-M-15 and it was present in all five phylogroups and 11 STs, and most commonly detected in ST391 (4/4 isolates), ST648 (3/3 isolates) and ST131 [3/14 (21.4%) isolates].

Conclusions

This study has revealed a high diversity of lineages associated with AMR, including ESBL and fluoroquinolone resistance, in Malawi. The data highlight the value of longitudinal bacteraemia surveillance coupled with detailed molecular epidemiology in all settings, including low-income settings, in describing the global epidemiology of ESBL resistance.

The Microevolution and Epidemiology of Staphylococcus aureus Colonization during Atopic Eczema Disease Flare

Staphylococcus aureus is an opportunistic pathogen and variable component of the human microbiota. A characteristic of atopic eczema (AE) is colonization by S. aureus, with exacerbations associated with an increased bacterial burden of the organism. Despite this, the origins and genetic diversity of S. aureus colonizing individual patients during AE disease flares is poorly understood. To examine the microevolution of S. aureus colonization, we deep sequenced S. aureus populations from nine children with moderate to severe AE and 18 non-atopic children asymptomatically carrying S. aureus nasally. Colonization by clonal S. aureus populations was observed in both AE patients and control participants, with all but one of the individuals carrying colonies belonging to a single sequence type. Phylogenetic analysis showed that disease flares were associated with the clonal expansion of the S. aureus population, occurring over a period of weeks to months. There was a significant difference in the genetic backgrounds of S. aureus colonizing AE cases versus controls (Fisher exact test, P = 0.03). Examination of intra-host genetic heterogeneity of the colonizing S. aureus populations identified evidence of within-host selection in the AE patients, with AE variants being potentially selectively advantageous for intracellular persistence and treatment resistance.

ARIBA: rapid antimicrobial resistance genotyping directly from sequencing reads

Antimicrobial resistance (AMR) is one of the major threats to human and animal health worldwide, yet few high-throughput tools exist to analyse and predict the resistance of a bacterial isolate from sequencing data. Here we present a new tool, ARIBA, that identifies AMR-associated genes and single nucleotide polymorphisms directly from short reads, and generates detailed and customizable output. The accuracy and advantages of ARIBA over other tools are demonstrated on three datasets from Gram-positive and Gram-negative bacteria, with ARIBA outperforming existing methods.

Evolution of mobile genetic element composition in an epidemic methicillin-resistant Staphylococcus aureus: temporal changes correlated with frequent loss and gain events

BACKGROUND

Horizontal transfer of mobile genetic elements (MGEs) that carry virulence and antimicrobial resistance genes mediates the evolution of methicillin-resistant Staphylococcus aureus, and the emergence of new MRSA clones. Most MRSA lineages show an association with specific MGEs and the evolution of MGE composition following clonal expansion has not been widely studied.

RESULTS

We investigated the genomes of 1193 S. aureus bloodstream isolates, 1169 of which were MRSA, collected in the UK and the Republic of Ireland between 2001 and 2010. The majority of isolates belonged to clonal complex (CC)22 (n = 923), which contained diverse MGEs including elements that were found in other MRSA lineages. Several MGEs showed variable distribution across the CC22 phylogeny, including two antimicrobial resistance plasmids (pWBG751-like and SAP078A-like, carrying erythromycin and heavy metal resistance genes, respectively), a pathogenicity island carrying the enterotoxin C gene and two phage types Sa1int and Sa6int. Multiple gains and losses of these five MGEs were identified in the CC22 phylogeny using ancestral state reconstruction. Analysis of the temporal distribution of the five MGEs between 2001 and 2010 revealed an unexpected reduction in prevalence of the two plasmids and the pathogenicity island, and an increase in the two phage types. This occurred across the lineage and was not correlated with changes in the relative prevalence of CC22, or of any sub-lineages within in.

CONCLUSIONS

Ancestral state reconstruction coupled with temporal trend analysis demonstrated that epidemic MRSA CC22 has an evolving MGE composition, and indicates that this important MRSA lineage has continued to adapt to changing selective pressure since its emergence.

Patterns of antimicrobial resistance in Streptococcus suis isolates from pigs with or without streptococcal disease in England between 2009 and 2014

Antimicrobial resistance in Streptococcus suis, a global zoonotic pathogen of pigs, has been mostly studied only in diseased animals using surveys that have not evaluated changes over time. We compared patterns of resistance between S. suis isolates from clinical cases of disease (CC) and non-clinical case (NCC) pigs in England, collected over two discrete periods, 2009-2011 and 2013-2014. Minimum inhibitory concentrations (MIC) of 17 antimicrobials (nine classes) were determined on 405 S. suis isolates categorised by sampling period and disease association to assess changes in resistance over time and association with disease. First, isolates were characterized as resistant or susceptible using published clinical breakpoints. Second, epidemiological cut-offs (ECOFF) were derived from MIC values, and isolates classified as wild type (WT) below the ECOFF and non-wild type (NWT) above the ECOFF. Finally, isolate subsets were analysed for shifts in MIC distribution. NCC isolates were more resistant than CC isolates to cephalosporins, penams, pleuromutilins, potentiated sulphonamides and tetracyclines in both study periods. Resistance levels among CC isolates increased in 2013-2014 relative to 2009-2011 for antimicrobials including aminoglycosides, cephalosporins, fluoroquinolones, pleuromutilins, potentiated sulphonamides and tetracyclines. The prevalence of isolates categorised as NWT for five or more classes of antimicrobials was greater among NCC than CC isolates for both time periods, and increased with time. This study used standardised methods to identify significant shifts in antimicrobial resistance phenotypes of S. suis isolated from pigs in England, not only over time but also between isolates from known clinical cases or disease-free pigs.

Whole genome sequencing of Shigella sonnei through PulseNet Latin America and Caribbean: advancing global surveillance of foodborne illnesses

Objectives

Shigella sonnei is a globally important diarrhoeal pathogen tracked through the surveillance network PulseNet Latin America and Caribbean (PNLA&C), which participates in PulseNet International. PNLA&C laboratories use common molecular techniques to track pathogens causing foodborne illness. We aimed to demonstrate the possibility and advantages of transitioning to whole genome sequencing (WGS) for surveillance within existing networks across a continent where S. sonnei is endemic.

Methods

We applied WGS to representative archive isolates of S. sonnei (n = 323) from laboratories in nine PNLA&C countries to generate a regional phylogenomic reference for S. sonnei and put this in the global context. We used this reference to contextualise 16 S. sonnei from three Argentinian outbreaks, using locally generated sequence data. Assembled genome sequences were used to predict antimicrobial resistance (AMR) phenotypes and identify AMR determinants.

Results

S. sonnei isolates clustered in five Latin American sublineages in the global phylogeny, with many (46%, 149 of 323) belonging to previously undescribed sublineages. Predicted multidrug resistance was common (77%, 249 of 323), and clinically relevant differences in AMR were found among sublineages. The regional overview showed that Argentinian outbreak isolates belonged to distinct sublineages and had different epidemiologic origins.

Conclusions

Latin America contains novel genetic diversity of S. sonnei that is relevant on a global scale and commonly exhibits multidrug resistance. Retrospective passive surveillance with WGS has utility for informing treatment, identifying regionally epidemic sublineages and providing a framework for interpretation of prospective, locally sequenced outbreaks.

Global and regional dissemination and evolution of Burkholderia pseudomallei

The environmental bacterium Burkholderia pseudomallei causes an estimated 165,000 cases of human melioidosis per year worldwide and is also classified as a biothreat agent. We used whole genome sequences of 469 B. pseudomallei isolates from 30 countries collected over 79 years to explore its geographic transmission. Our data point to Australia as an early reservoir, with transmission to Southeast Asia followed by onward transmission to South Asia and East Asia. Repeated reintroductions were observed within the Malay Peninsula and between countries bordered by the Mekong River. Our data support an African origin of the Central and South American isolates with introduction of B. pseudomallei into the Americas between 1650 and 1850, providing a temporal link with the slave trade. We also identified geographically distinct genes/variants in Australasian or Southeast Asian isolates alone, with virulence-associated genes being among those over-represented. This provides a potential explanation for clinical manifestations of melioidosis that are geographically restricted.

Microevolution of Monophasic Salmonella Typhimurium during Epidemic, United Kingdom, 2005-2010

Microevolution resulted in considerable genotypic variation.

Microevolution associated with emergence and expansion of new epidemic clones of bacterial pathogens holds the key to epidemiologic success. To determine microevolution associated with monophasic Salmonella Typhimurium during an epidemic, we performed comparative whole-genome sequencing and phylogenomic analysis of isolates from the United Kingdom and Italy during 2005–2012. These isolates formed a single clade distinct from recent monophasic epidemic clones previously described from North America and Spain. The UK monophasic epidemic clones showed a novel genomic island encoding resistance to heavy metals and a composite transposon encoding antimicrobial drug resistance genes not present in other Salmonella Typhimurium isolates, which may have contributed to epidemiologic success. A remarkable amount of genotypic variation accumulated during clonal expansion that occurred during the epidemic, including multiple independent acquisitions of a novel prophage carrying the sopE gene and multiple deletion events affecting the phase II flagellin locus. This high level of microevolution may affect antigenicity, pathogenicity, and transmission.

Genomic variations leading to alterations in cell morphology of Campylobacter spp

Campylobacter jejuni, the most common cause of bacterial diarrhoeal disease, is normally helical. However, it can also adopt straight rod, elongated helical and coccoid forms. Studying how helical morphology is generated, and how it switches between its different forms, is an important objective for understanding this pathogen. Here, we aimed to determine the genetic factors involved in generating the helical shape of Campylobacter. A C. jejuni transposon (Tn) mutant library was screened for non-helical mutants with inconsistent results. Whole genome sequence variation and morphological trends within this Tn library, and in various C. jejuni wild type strains, were compared and correlated to detect genomic elements associated with helical and rod morphologies. All rod-shaped C. jejuni Tn mutants and all rod-shaped laboratory, clinical and environmental C. jejuni and Campylobacter coli contained genetic changes within the pgp1 or pgp2 genes, which encode peptidoglycan modifying enzymes. We therefore confirm the importance of Pgp1 and Pgp2 in the maintenance of helical shape and extended this to a wide range of C. jejuni and C. coli isolates. Genome sequence analysis revealed variation in the sequence and length of homopolymeric tracts found within these genes, providing a potential mechanism of phase variation of cell shape.

Factors determining human-to-human transmissibility of zoonotic pathogens via contact

•

There are several modes of contact transmission of pathogens amongst humans.

•

Factors promoting contact transmission act at the pathogen, host or environmental levels.

•

Common pathogen factors are immune evasion, high viral load and low infectious dose.

•

Common host factors are crowding, promiscuity and the presence of co-infections.

The pandemic potential of zoonotic pathogens lies in their ability to become efficiently transmissible amongst humans. Here, we focus on contact-transmitted pathogens and discuss the factors, at the pathogen, host and environmental levels that promote or hinder their human-to-human transmissibility via the following modes of contact transmission: skin contact, sexual contact, respiratory contact and multiple route contact. Factors common to several modes of transmission were immune evasion, high viral load, low infectious dose, crowding, promiscuity, and co-infections; other factors were specific for a pathogen or mode of contact transmission. The identification of such factors will lead to a better understanding of the requirements for human-to-human spread of pathogens, as well as improving risk assessment of newly emerging pathogens.

Molecular Approaches to Understanding Transmission and Source Attribution in Nontyphoidal Salmonella and Their Application in Africa

Nontyphoidal Salmonella (NTS) is a frequent cause of diarrhea around the world, yet in many African countries it is more commonly associated with invasive bacterial disease. Various source attribution models have been developed that utilize microbial subtyping data to assign cases of human NTS infection to different animal populations and foods of animal origin. Advances in molecular microbial subtyping approaches, in particular whole-genome sequencing, provide higher resolution data with which to investigate these sources. In this review, we provide updates on the source attribution models developed for Salmonella, and examine the application of whole-genome sequencing data combined with evolutionary modeling to investigate the putative sources and transmission pathways of NTS, with a focus on the epidemiology of NTS in Africa. This is essential information to decide where, what, and how control strategies might be applied most effectively.