Mashtree: a rapid comparison of whole genome sequence files

Novel oceanic cyanobacterium isolated from Bangaram island with profound acid neutralizing ability is proposed as Leptolyngbya iicbica sp. nov. strain LK

An acid-neutralizing, filamentous, non-heterocytous, marine cyanobacterium named 'LK' has been isolated from the seashore of Bangaram Island, an atoll of Lakshadweep, India, and is described here as a novel species. LK has been characterized using morphological, ecological, and genomic features. Based on 16S rRNA, whole-genome sequencing, and marker gene-based analysis, LK has been identified as a new species. LK clustered with Leptolyngbya-like strains belonging to the LPP group but diverged from Leptolyngbya sensu stricto, indicating the polyphyletic nature of the Leptolyngbya genus. Leptolyngbya sp. SIOISBB and Halomicronema sp. CCY15110 were identified as LK's two closest phylogenetic neighbors in various phylogenetic studies. The analysis of 16S rRNA, ITS secondary structures, and genome relatedness indices such as AAI, ANI, and gANI strongly support LK as a novel species of the Leptolyngbya genus. The mechanism behind acid neutralization in LK has been delineated, attributing it to a surface phenomenon most likely due to the presence of salts of calcium, magnesium, sodium, and potassium. We name LK as Leptolyngbya iicbica strain LK which is a novel species with prominent acidic pH-neutralizing properties.

Detection of Leptospira kirschneri in a short-beaked common dolphin (Delphinus delphis delphis) stranded off the coast of southern California, USA

BACKGROUND

Pathogenic Leptospira species are globally important zoonotic pathogens capable of infecting a wide range of host species. In marine mammals, reports of Leptospira have predominantly been in pinnipeds, with isolated reports of infections in cetaceans.

CASE PRESENTATION

On 28 June 2021, a 150.5 cm long female, short-beaked common dolphin (Delphinus delphis delphis) stranded alive on the coast of southern California and subsequently died. Gross necropsy revealed multifocal cortical pallor within the reniculi of the kidney, and lymphoplasmacytic tubulointerstitial nephritis was observed histologically. Immunohistochemistry confirmed Leptospira infection, and PCR followed by lfb1 gene amplicon sequencing suggested that the infecting organism was L.kirschneri. Leptospira DNA capture and enrichment allowed for whole-genome sequencing to be conducted. Phylogenetic analyses confirmed the causative agent was a previously undescribed, divergent lineage of L.kirschneri.

CONCLUSIONS

We report the first detection of pathogenic Leptospira in a short-beaked common dolphin, and the first detection in any cetacean in the northeastern Pacific Ocean. Renal lesions were consistent with leptospirosis in other host species, including marine mammals, and were the most significant lesions detected overall, suggesting leptospirosis as the likely cause of death. We identified the cause of the infection as L.kirschneri, a species detected only once before in a marine mammal - a northern elephant seal (Mirounga angustirostris) of the northeastern Pacific. These findings raise questions about the mechanism of transmission, given the obligate marine lifestyle of cetaceans (in contrast to pinnipeds, which spend time on land) and the commonly accepted view that Leptospira are quickly killed by salt water. They also raise important questions regarding the source of infection, and whether it arose from transmission among marine mammals or from terrestrial-to-marine spillover. Moving forward, surveillance and sampling must be expanded to better understand the extent to which Leptospira infections occur in the marine ecosystem and possible epidemiological linkages between and among marine and terrestrial host species. Generating Leptospira genomes from different host species will yield crucial information about possible transmission links, and our study highlights the power of new techniques such as DNA enrichment to illuminate the complex ecology of this important zoonotic pathogen.

Spread of antibiotic resistance genes in drinking water reservoirs: Insights from a deep metagenomic study using a curated database

The exploration of antibiotic resistance genes (ARGs) in drinking water reservoirs is an emerging field. Using a curated database, we enhanced the ARG detection and conducted a comprehensive analysis using 2.2 Tb of deep metagenomic sequencing data to determine the distribution of ARGs across 16 drinking water reservoirs and associated environments. Our findings reveal a greater diversity of ARGs in sediments than in water, underscoring the importance of extensive background surveys. Crucial ARG carriers-specifically Acinetobacter, Pseudomonas, and Mycobacterium were identified in drinking water reservoirs. Extensive analysis of the data uncovered a considerable concern for drinking water safety, particularly in regions reliant on river sources. Mobile genetic elements have been found to contribute markedly to the propagation of ARGs. The results of this research suggest that the establishment of drinking water reservoirs for supplying raw water may be an effective strategy for alleviating the spread of water-mediated ARGs.

Clinical, phenotypic, and genotypic characteristics of ESBL-producing Salmonella enterica bloodstream infections from Qatar

Background

Resistant Salmonella infections are a major global public health challenge particularly for multidrug-resistant (MDR) isolates manifesting as bloodstream infections (BSIs).

Objectives

To evaluate clinical, phenotypic, and genotypic characteristics of extended-spectrum beta-lactamase (ESBL) producing Salmonella enterica BSIs from Qatar.

Methods

Phenotypic ESBL Salmonella enterica from adult patients presenting with positive BSIs were collected between January 2019 to May 2020. Microbiological identification and characterization were performed using standard methods while genetic characteristics were examined through whole genome sequencing studies.

Results

Of 151 episodes of Salmonella enterica BSI, 15 (10%) phenotypic ESBL isolates were collected. Recent travel was recorded in most cases (80%) with recent exposure to antimicrobials (27%). High-level resistance to quinolines, aminoglycosides, and cephalosporins was recorded (80-100%) while meropenem, tigecycline and colistin demonstrated universal susceptibility. Genomic evaluation demonstrated dominance of serotype Salmonella Typhi sequence type 1 (93%) while antimicrobial resistance genes revealed dominance of aminoglycoside resistance (100%), qnrS1 quinolones resistance (80%), blaCTX-M-15 ESBLs (86.7%), and paucity of AmpC resistance genes (6.7%).

Conclusions

Invasive MDR Salmonella enterica is mainly imported, connected to patients from high prevalent regions with recent travel and antimicrobial use caused by specific resistant clones. In suspected cases of multidrug resistance, carbapenem therapy is recommended.

Accurate, scalable, and fully automated inference of species trees from raw genome assemblies using ROADIES

Inference of species trees plays a crucial role in advancing our understanding of evolutionary relationships and has immense significance for diverse biological and medical applications. Extensive genome sequencing efforts are currently in progress across a broad spectrum of life forms, holding the potential to unravel the intricate branching patterns within the tree of life. However, estimating species trees starting from raw genome sequences is quite challenging, and the current cutting-edge methodologies require a series of error-prone steps that are neither entirely automated nor standardized. In this paper, we present ROADIES, a novel pipeline for species tree inference from raw genome assemblies that is fully automated, easy to use, scalable, free from reference bias, and provides flexibility to adjust the tradeoff between accuracy and runtime. The ROADIES pipeline eliminates the need to align whole genomes, choose a single reference species, or pre-select loci such as functional genes found using cumbersome annotation steps. Moreover, it leverages recent advances in phylogenetic inference to allow multi-copy genes, eliminating the need to detect orthology. Using the genomic datasets released from large-scale sequencing consortia across three diverse life forms (placental mammals, pomace flies, and birds), we show that ROADIES infers species trees that are comparable in quality with the state-of-the-art approaches but in a fraction of the time. By incorporating optimal approaches and automating all steps from assembled genomes to species and gene trees, ROADIES is poised to improve the accuracy, scalability, and reproducibility of phylogenomic analyses.

Efficient and Robust Search of Microbial Genomes via Phylogenetic Compression

Comprehensive collections approaching millions of sequenced genomes have become central information sources in the life sciences. However, the rapid growth of these collections has made it effectively impossible to search these data using tools such as BLAST and its successors. Here, we present a technique called phylogenetic compression, which uses evolutionary history to guide compression and efficiently search large collections of microbial genomes using existing algorithms and data structures. We show that, when applied to modern diverse collections approaching millions of genomes, lossless phylogenetic compression improves the compression ratios of assemblies, de Bruijn graphs, and k -mer indexes by one to two orders of magnitude. Additionally, we develop a pipeline for a BLAST-like search over these phylogeny-compressed reference data, and demonstrate it can align genes, plasmids, or entire sequencing experiments against all sequenced bacteria until 2019 on ordinary desktop computers within a few hours. Phylogenetic compression has broad applications in computational biology and may provide a fundamental design principle for future genomics infrastructure.

Identification of equine mares as reservoir hosts for pathogenic species of Leptospira

Equine leptospirosis can result in abortion, stillbirth, neonatal death, placentitis, and uveitis. Horses can also act as subclinical reservoir hosts of infection, which are characterized as asymptomatic carriers that persistently excrete leptospires and transmit disease. In this study, PCR and culture were used to assess urinary shedding of pathogenic Leptospira from 37 asymptomatic mares. Three asymptomatic mares, designated as H2, H8, and H9, were PCR-positive for lipL32, a gene specific for pathogenic species of Leptospira. One asymptomatic mare, H9, was culture-positive, and the recovered isolate was classified as L. kirschneri serogroup Australis serovar Rushan. DNA capture and enrichment of Leptospira genomic DNA from PCR-positive, culture-negative samples determined that asymptomatic mare H8 was also shedding L. kirschneri serogroup Australis, whereas asymptomatic mare H2 was shedding L. interrogans serogroup Icterohaemorrhagiae. Sera from all asymptomatic mares were tested by the microscopic agglutination test (MAT) and 35 of 37 (94.6%) were seropositive with titers ranging from 1:100 to 1:3200. In contrast to asymptomatic mares, mare H44 presented with acute spontaneous abortion and a serum MAT titer of 1:102,400 to L. interrogans serogroup Pomona serovar Pomona. Comparison of L. kirschneri serogroup Australis strain H9 with that of L. interrogans serogroup Pomona strain H44 in the hamster model of leptospirosis corroborated differences in virulence of strains. Since lipopolysaccharide (LPS) is a protective antigen in bacterin vaccines, the LPS of strain H9 (associated with subclinical carriage) was compared with strain H44 (associated with spontaneous abortion). This revealed different LPS profiles and immunoreactivity with reference antisera. It is essential to know what species and serovars of Leptospira are circulating in equine populations to design efficacious vaccines and diagnostic tests. Our results demonstrate that horses in the US can act as reservoir hosts of leptospirosis and shed diverse pathogenic Leptospira species via urine. This report also details the detection of L. kirschneri serogroup Australis serovar Rushan, a species and serotype of Leptospira, not previously reported in the US.

Optimized model architectures for deep learning on genomic data

The success of deep learning in various applications depends on task-specific architecture design choices, including the types, hyperparameters, and number of layers. In computational biology, there is no consensus on the optimal architecture design, and decisions are often made using insights from more well-established fields such as computer vision. These may not consider the domain-specific characteristics of genome sequences, potentially limiting performance. Here, we present GenomeNet-Architect, a neural architecture design framework that automatically optimizes deep learning models for genome sequence data. It optimizes the overall layout of the architecture, with a search space specifically designed for genomics. Additionally, it optimizes hyperparameters of individual layers and the model training procedure. On a viral classification task, GenomeNet-Architect reduced the read-level misclassification rate by 19%, with 67% faster inference and 83% fewer parameters, and achieved similar contig-level accuracy with ~100 times fewer parameters compared to the best-performing deep learning baselines.

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

In-depth characterization of multidrug-resistant NDM-1 and KPC-3 co-producing Klebsiella pneumoniae bloodstream isolates from Italian hospital patients

Bloodstream infection (BSI) caused by carbapenem-resistant Klebsiella pneumoniae (KP) poses significant challenges, particularly when the infecting isolate carries multiple antimicrobial resistance (AMR) genes/determinants. This study, employing short- and long-read whole-genome sequencing, characterizes six New Delhi metallo-β-lactamase (NDM) 1 and KP carbapenemase (KPC) 3 co-producing KP isolates, the largest cohort investigated in Europe to date. Five [sequence type (ST) 512] and one (ST11) isolates were recovered from patients who developed BSI from February to August 2022 or February 2023 at two different hospitals in Rome, Italy. Phylogenetic analysis revealed two distinct clusters among ST512 isolates and a separate cluster for the ST11 isolate. Beyond blaNDM-1 and blaKPC-3, various AMR genes, indicative of a multidrug resistance phenotype, including colistin resistance, were found. Each cluster-representative ST512 isolate harbored a blaNDM-1 plasmid (IncC) and a blaKPC-3 plasmid [IncFIB(pQil)/IncFII(K)], while the ST11 isolate harbored a blaNDM-1 plasmid [IncFII(pKPX1)] and a blaKPC-3 plasmid [IncFIB(K)/IncFII(K)]. The blaNDM-1 plasmids carried genes conferring resistance to clinically relevant antimicrobial agents, and the aminoglycoside resistance gene aac(6')-Ib was found on different plasmids. Colistin resistance-associated mgrB/pmrB gene mutations were present in all isolates, and the yersiniabactin-encoding ybt gene was unique to the ST11 isolate. In conclusion, our findings provide insights into the genomic context of blaNDM-1/blaKPC-3 carbapenemase-producing KP isolates.IMPORTANCEThis study underscores the critical role of genomic surveillance as a proactive measure to restrict the spread of carbapenemase-producing KP isolates, especially when key antimicrobial resistance genes, such as blaNDM-1/blaKPC-3, are plasmid borne. In-depth characterization of these isolates may help identify plasmid similarities contributing to their intra-hospital/inter-hospital adaptation and transmission. Despite the lack of data on patient movements, it is possible that carbapenem-resistant isolates were selected to co-produce KP carbapenemase and New Delhi metallo-β-lactamase via plasmid acquisition. Studies employing long-read whole-genome sequencing should be encouraged to address the emergence of KP clones with converging phenotypes of virulence and resistance to last-resort antimicrobial agents.

Nanopore and Illumina sequencing reveal different viral populations from human gut samples

The advent of viral metagenomics, or viromics, has improved our knowledge and understanding of global viral diversity. High-throughput sequencing technologies enable explorations of the ecological roles, contributions to host metabolism, and the influence of viruses in various environments, including the human intestinal microbiome. However, bacterial metagenomic studies frequently have the advantage. The adoption of advanced technologies like long-read sequencing has the potential to be transformative in refining viromics and metagenomics. Here, we examined the effectiveness of long-read and hybrid sequencing by comparing Illumina short-read and Oxford Nanopore Technology (ONT) long-read sequencing technologies and different assembly strategies on recovering viral genomes from human faecal samples. Our findings showed that if a single sequencing technology is to be chosen for virome analysis, Illumina is preferable due to its superior ability to recover fully resolved viral genomes and minimise erroneous genomes. While ONT assemblies were effective in recovering viral diversity, the challenges related to input requirements and the necessity for amplification made it less ideal as a standalone solution. However, using a combined, hybrid approach enabled a more authentic representation of viral diversity to be obtained within samples.

Emergence of OXA-48-like producing Citrobacter species, Germany, 2011 to 2022

BackgroundCarbapenemase-producing Enterobacterales are a public health threat worldwide and OXA-48 is the most prevalent carbapenemase in Germany and western Europe. However, the molecular epidemiology of OXA-48 in species other than Escherichia coli and Klebsiella pneumoniae remains poorly understood.AimTo analyse the molecular epidemiology of OXA-48 and OXA-48-like carbapenemases in Citrobacter species (spp.) in Germany between 2011 and 2022.MethodsData of 26,822 Enterobacterales isolates sent to the National Reference Centre (NRC) for Gram-negative bacteria were evaluated. Ninety-one Citrobacter isolates from 40 German hospitals harbouring bla OXA-48/OXA-48‑like were analysed by whole genome sequencing and conjugation experiments.ResultsThe frequency of OXA-48 in Citrobacter freundii (CF) has increased steadily since 2011 and is now the most prevalent carbapenemase in this species in Germany. Among 91 in-depth analysed Citrobacter spp. isolates, CF (n = 73) and C. koseri (n = 8) were the most common species and OXA-48 was the most common variant (n = 77), followed by OXA-162 (n = 11) and OXA‑181 (n = 3). Forty percent of the isolates belonged to only two sequence types (ST19 and ST22), while most other STs were singletons. The plasmids harbouring bla OXA‑48 and bla OXA-162 belonged to the plasmid types IncL (n = 85) or IncF (n = 3), and plasmids harbouring bla OXA‑181 to IncX3 (n = 3). Three IncL plasmid clusters (57/85 IncL plasmids) were identified, which were highly transferable in contrast to sporadic plasmids.ConclusionIn CF in Germany, OXA-48 is the predominant carbapenemase. Dissemination is likely due to distinct highly transmissible plasmids harbouring bla OXA‑48 or bla OXA-48-like and the spread of the high-risk clonal lineages ST19 and ST22.

Towards personalised anti-microbial and immune approaches to infections in acute care. Can real-time genomic-informed diagnosis of pathogens, and immune-focused therapies improve outcomes for patients? An observational, experimental study protocol

INTRODUCTION

Infection causes a vast burden of disease, with significant mortality, morbidity and costs to health-care systems. However, identifying the pathogen causative infection can be challenging, resulting in high use of broad-spectrum antibiotics, much of which may be inappropriate. Novel metagenomic methods have potential to rapidly identify pathogens, however their clinical utility for many infections is currently unclear. Outcome from infection is also impacted by the effectiveness of immune responses, which can be impaired by age, co-morbidity and the infection itself. The aims of this study are twofold: To compare diversity of organisms identified and time-to-result using metagenomic methods versus traditional culture -based techniques, to explore the potential clinical role of metagenomic approaches to pathogen identification in a range of infections.To characterise the ex vivo function of immune cells from patients with acute infection, exploring host and pathogen-specific factors which may affect immune function and overall outcomes.

METHODS

This is a prospective observational study of patients with acute infection. Patients with symptoms suggestive of an acute infection will be recruited, and blood and bodily fluid relevant to the site of infection collected (for example, sputum and naso-oropharyngeal swabs for respiratory tract infections, or urine for a suspected urinary tract infection). Metagenomic analysis of samples will be compared to traditional microbiology, alongside the antimicrobials received. Blood and respiratory samples such as bronchoalveolar lavage will be used to isolate immune cells and interrogate immune cell function. Where possible, similar samples will be collected from matched participants without a suspected infection to determine the impact of infection on both microbiome and immune cell function.

Pseudomonas aeruginosa MipA-MipB envelope proteins act as new sensors of polymyxins

Due to the rising incidence of antibiotic-resistant infections, the last-line antibiotics, polymyxins, have resurged in the clinics in parallel with new bacterial strategies of escape. The Gram-negative opportunistic pathogen Pseudomonas aeruginosa develops resistance to colistin/polymyxin B by distinct molecular mechanisms, mostly through modification of the lipid A component of the LPS by proteins encoded within the arnBCDATEF-ugd (arn) operon. In this work, we characterized a polymyxin-induced operon named mipBA, present in P. aeruginosa strains devoid of the arn operon. We showed that mipBA is activated by the ParR/ParS two-component regulatory system in response to polymyxins. Structural modeling revealed that MipA folds as an outer-membrane β-barrel, harboring an internal negatively charged channel, able to host a polymyxin molecule, while the lipoprotein MipB adopts a β-lactamase fold with two additional C-terminal domains. Experimental work confirmed that MipA and MipB localize to the bacterial envelope, and they co-purify in vitro. Nano differential scanning fluorimetry showed that polymyxins stabilized MipA in a specific and dose-dependent manner. Mass spectrometry-based quantitative proteomics on P. aeruginosa membranes demonstrated that ∆mipBA synthesized fourfold less MexXY-OprA proteins in response to polymyxin B compared to the wild-type strain. The decrease was a direct consequence of impaired transcriptional activation of the mex operon operated by ParR/ParS. We propose MipA/MipB to act as membrane (co)sensors working in concert to activate ParS histidine kinase and help the bacterium to cope with polymyxin-mediated envelope stress through synthesis of the efflux pump, MexXY-OprA.IMPORTANCEDue to the emergence of multidrug-resistant isolates, antibiotic options may be limited to polymyxins to eradicate Gram-negative infections. Pseudomonas aeruginosa, a leading opportunistic pathogen, has the ability to develop resistance to these cationic lipopeptides by modifying its lipopolysaccharide through proteins encoded within the arn operon. Herein, we describe a sub-group of P. aeruginosa strains lacking the arn operon yet exhibiting adaptability to polymyxins. Exposition to sub-lethal polymyxin concentrations induced the expression and production of two envelope-associated proteins. Among those, MipA, an outer-membrane barrel, is able to specifically bind polymyxins with an affinity in the 10-µM range. Using membrane proteomics and phenotypic assays, we showed that MipA and MipB participate in the adaptive response to polymyxins via ParR/ParS regulatory signaling. We propose a new model wherein the MipA-MipB module functions as a novel polymyxin sensing mechanism.

Genetic characterization of multidrug-resistant Escherichia coli harboring colistin-resistant gene isolated from food animals in food supply chain

Colistin is widely used for the prophylaxis and treatment of infectious disease in humans and livestock. However, the global food chain may actively promote the dissemination of colistin-resistant bacteria in the world. Mobile colistin-resistant (mcr) genes have spread globally, in both communities and hospitals. This study sought to genomically characterize mcr-mediated colistin resistance in 16 Escherichia coli strains isolated from retail meat samples using whole genome sequencing with short-read and long-read platforms. To assess colistin resistance and the transferability of mcr genes, antimicrobial susceptibility testing and conjugation experiments were conducted. Among the 16 isolates, 11 contained mcr-1, whereas three carried mcr-3 and two contained mcr-1 and mcr-3. All isolates had minimum inhibitory concentration (MIC) for colistin in the range 1-64 μg/mL. Notably, 15 out of the 16 isolates demonstrated successful transfer of mcr genes via conjugation, indicative of their presence on plasmids. In contrast, the KK3 strain did not exhibit such transferability. Replicon types of mcr-1-containing plasmids included IncI2 and IncX4, while IncFIB, IncFII, and IncP1 contained mcr-3. Another single strain carried mcr-1.1 on IncX4 and mcr-3.5 on IncP1. Notably, one isolate contained mcr-1.1 located on a chromosome and carrying mcr-3.1 on the IncFIB plasmid. The chromosomal location of the mcr gene may ensure a steady spread of resistance in the absence of selective pressure. Retail meat products may act as critical reservoirs of plasmid-mediated colistin resistance that has been transmitted to humans.

Genomic surveillance of extended-spectrum cephalosporin-resistant Escherichia coli isolated from poultry in the UK from 2016 to 2020

Introduction

Surveillance is vital for monitoring the increasing risk of antimicrobial resistance (AMR) in bacteria leading to failures in humans and animals to treat infections. In a One Health context, AMR bacteria from livestock and food can transfer through the food chain to humans, and vice versa, which can be characterized in detail through genomics. We investigated the critical aspects of AMR and the dynamics of AMR in poultry in the UK.

Methods

In this study, we performed whole genome sequencing for genomic characterization of 761 extended-spectrum cephalosporinases (ESCs) harboring Escherichia coli isolated from poultry caeca and meat through EU harmonized monitoring of AMR in zoonotic and commensal bacteria from 2016 and 2018 and UK national monitoring in 2020.

Results

The most common ESC in 2016 and 2018 was blaCTX-M-1; however, 2020 had a greater diversity of ESCs with blaCTX-M-55 dominant in chickens and blaCTX-M-15 more prevalent in turkeys. Co-resistance to sulphonamides, tetracycline, and trimethoprim was widespread, and there were several positive correlations between the sequence types (STs) and ESC genes. We identified certain AMR genotypes and STs that were frequent each year but not as successful in subsequent years, e.g., ST350 harboring blaCTX-M-1, sul2, and tetA-v4.Phylogenetic comparison of isolates prevalent in our panel with global ones from the same STs available in public databases showed that isolates from the UK generally clustered together, suggesting greater within-country than between-country transmission.

Discussion

We conclude that future genomic surveillance of indicator organisms will be invaluable as it will enable detailed comparisons of AMR between and within neighboring countries, potentially identifying the most successful sequence types, plasmids, or emerging threats.

Genomic analysis of Citrobacter from Australian wastewater and silver gulls reveals novel sequence types carrying critically important antibiotic resistance genes

Antimicrobial resistance (AMR) is a major public health concern, and environmental bacteria have been recognized as important reservoirs of antimicrobial resistance genes (ARGs). Citrobacter, a common environmental bacterium and opportunistic pathogen in humans and other animals, has been largely understudied in terms of its diversity and AMR potential. Whole-genome (short-read) sequencing on a total of 77 Citrobacter isolates obtained from Australian silver gull (Chroicocephalus novaehollandiae) (n = 17) and influent wastewater samples (n = 60) was performed, revealing a diverse Citrobacter population, with seven different species and 33 sequence types, 17 of which were novel. From silver gull using non-selective media we isolated a broader range of species with little to no mobilised ARG carriage. Wastewater isolates (selected using Carbapenem- Resistant Enterobacterales (CRE) selective media) carried a heavy burden of ARGs (up to 21 ARGs, conferring resistance to nine classes of antibiotics), with several novel multidrug-resistant (MDR) lineages identified, including C. braakii ST1110, which carried ARGs conferring resistance to eight to nine classes of antibiotics, and C. freundii ST1105, which carried two carbapenemase genes, blaIMP-4 in class 1 integron structure, and blaKPC-2. Additionally, we identified an MDR C. portucalensis isolate carrying blaNDM-1, blaSHV-12, and mcr-9. We identified IncC, IncM2, and IncP6 plasmids as the likely vectors for many of the critically important mobilised ARGs. Phylogenetic analyses were performed to assess any epidemiological linkages between isolation sources, demonstrating low relatedness across sources beyond the ST level. However, these analyses did reveal some closer relationships between strains from disparate wastewater sources despite their collection some 13,000 km apart. These findings support the need for future surveillance of Citrobacter populations in wastewater and wildlife populations to monitor for potential opportunistic human pathogens.

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

Characterization of nontyphoidal Salmonella strains from a tertiary hospital in China: serotype diversity, multidrug resistance, and genetic insights

Objective

Nontyphoidal Salmonella is a significant public health concern due to its ability to cause foodborne illnesses worldwide. This study aims to characterize the nontyphoidal Salmonella strains isolated from patients in China.

Methods

A total of 19 nontyphoidal Salmonella strains were characterized through serovar identification, antimicrobial susceptibility testing (AST), biofilm formation assessment. Genetic relatedness was determined using pulsed-field gel electrophoresis (PFGE). WGS was employed to decipher the resistance mechanism and to contextualize the S. serovar Mbandaka strains among previously sequenced isolates in China. The biofilm associated mrkA gene was examined by PCR.

Results

The predominant serovar identified was S. Enteritidis, followed by S. Mbandaka, S. Thompson, S. Livingston, S. Alachua, and S. Infantis. PFGE analysis indicated a notable genetic similarity among the S. Mbandaka isolates. Phylogenetic analysis suggested that these strains were likely derived from a single source that had persisted in China for over five years. One multidrug resistance (MDR) S. Enteritidis isolate carried a highly transferable IncB/O/K/Z plasmid with bla CTX-M-15. One S. Thompson strain, harboring the mrkABCDF operon in an IncX1 plasmid, isolated from cutaneous lesions, demonstrated robust biofilm formation. However, no mrkABCDF loci were detected in other strains.

Conclusion

Our study emphasizes the importance of persisted surveillance and prompt response to Salmonella infections to protect public health. The dissemination of bla CTX-M-15-harboring IncB/O/K/Z plasmid and the spread of virulent mrkABCDF operon among Salmonella in China and other global regions warrant close monitoring.

GWAS reveals a rapidly evolving candidate avirulence effector in the Cercospora leaf spot pathogen

The major resistance gene BvCR4 recently bred into sugar beet hybrids provides a high level of resistance to Cercospora leaf spot caused by the fungal pathogen Cercospora beticola. The occurrence of pathogen strains that overcome BvCR4 was studied using field trials in Switzerland conducted under natural disease pressure. Virulence of a subset of these strains was evaluated in a field trial conducted under elevated artificial disease pressure. We created a new C. beticola reference genome and mapped whole genome sequences of 256 isolates collected in Switzerland and Germany. These were combined with virulence phenotypes to conduct three separate genome-wide association studies (GWAS) to identify candidate avirulence genes. We identified a locus associated with avirulence containing a putative avirulence effector gene named AvrCR4. All virulent isolates either lacked AvrCR4 or had nonsynonymous mutations within the gene. AvrCR4 was present in all 74 isolates from non-BvCR4 hybrids, whereas 33 of 89 isolates from BvCR4 hybrids carried a deletion. We also mapped genomic data from 190 publicly available US isolates to our new reference genome. The AvrCR4 deletion was found in only one of 95 unique isolates from non-BvCR4 hybrids in the United States. AvrCR4 presents a unique example of an avirulence effector in which virulent alleles have only recently emerged. Most likely these were selected out of standing genetic variation after deployment of BvCR4. Identification of AvrCR4 will enable real-time screening of C. beticola populations for the emergence and spread of virulent isolates.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Antimicrobial resistance level and conjugation permissiveness shape plasmid distribution in clinical enterobacteria

Conjugative plasmids play a key role in the dissemination of antimicrobial resistance (AMR) genes across bacterial pathogens. AMR plasmids are widespread in clinical settings, but their distribution is not random, and certain associations between plasmids and bacterial clones are particularly successful. For example, the globally spread carbapenem resistance plasmid pOXA-48 can use a wide range of enterobacterial species as hosts, but it is usually associated with a small number of specific Klebsiella pneumoniae clones. These successful associations represent an important threat for hospitalized patients. However, knowledge remains limited about the factors determining AMR plasmid distribution in clinically relevant bacteria. Here, we combined in vitro and in vivo experimental approaches to analyze pOXA-48-associated AMR levels and conjugation dynamics in a collection of wild-type enterobacterial strains isolated from hospitalized patients. Our results revealed significant variability in these traits across different bacterial hosts, with Klebsiella spp. strains showing higher pOXA-48-mediated AMR and conjugation frequencies than Escherichia coli strains. Using experimentally determined parameters, we developed a simple mathematical model to interrogate the contribution of AMR levels and conjugation permissiveness to plasmid distribution in bacterial communities. The simulations revealed that a small subset of clones, combining high AMR levels and conjugation permissiveness, play a critical role in stabilizing the plasmid in different polyclonal microbial communities. These results help to explain the preferential association of plasmid pOXA-48 with K. pneumoniae clones in clinical settings. More generally, our study reveals that species- and strain-specific variability in plasmid-associated phenotypes shape AMR evolution in clinically relevant bacterial communities.

Klebsiella pneumoniae clinical isolates with features of both multidrug-resistance and hypervirulence have unexpectedly low virulence

Klebsiella pneumoniae has been classified into two types, classical K. pneumoniae (cKP) and hypervirulent K. pneumoniae (hvKP). cKP isolates are highly diverse and important causes of nosocomial infections; they include globally disseminated antibiotic-resistant clones. hvKP isolates are sensitive to most antibiotics but are highly virulent, causing community-acquired infections in healthy individuals. The virulence phenotype of hvKP is associated with pathogenicity loci responsible for siderophore and hypermucoid capsule production. Recently, convergent strains of K. pneumoniae, which possess features of both cKP and hvKP, have emerged and are cause of much concern. Here, we screen the genomes of 2,608 multidrug-resistant K. pneumoniae isolates from the United States and identify 47 convergent isolates. We perform phenotypic and genomic characterization of 12 representative isolates. These 12 convergent isolates contain a variety of antimicrobial resistance plasmids and virulence plasmids. Most convergent isolates contain aerobactin biosynthesis genes and produce more siderophores than cKP isolates but not more capsule. Unexpectedly, only 1 of the 12 tested convergent isolates has a level of virulence consistent with hvKP isolates in a murine pneumonia model. These findings suggest that additional studies should be performed to clarify whether convergent strains are indeed more virulent than cKP in mouse and human infections.

Diversity, Taxonomic Novelty, and Encoded Functions of Salar de Ascotán Microbiota, as Revealed by Metagenome-Assembled Genomes

Salar de Ascotán is a high-altitude arsenic-rich salt flat exposed to high ultraviolet radiation in the Atacama Desert, Chile. It hosts unique endemic flora and fauna and is an essential habitat for migratory birds, making it an important site for conservation and protection. However, there is limited information on the resident microbiota's diversity, genomic features, metabolic potential, and molecular mechanisms that enable it to thrive in this extreme environment. We used long- and short-read metagenomics to investigate the microbial communities in Ascotán's water, sediment, and soil. Bacteria predominated, mainly Pseudomonadota, Acidobacteriota, and Bacteroidota, with a remarkable diversity of archaea in the soil. Following hybrid assembly, we recovered high-quality bacterial (101) and archaeal (6) metagenome-assembled genomes (MAGs), including representatives of two putative novel families of Patescibacteria and Pseudomonadota and two novel orders from the archaeal classes Halobacteriota and Thermoplasmata. We found different metabolic capabilities across distinct lineages and a widespread presence of genes related to stress response, DNA repair, and resistance to arsenic and other metals. These results highlight the remarkable diversity and taxonomic novelty of the Salar de Ascotán microbiota and its rich functional repertoire, making it able to resist different harsh conditions. The highly complete MAGs described here could serve future studies and bioprospection efforts focused on salt flat extremophiles, and contribute to enriching databases with microbial genome data from underrepresented regions of our planet.

Comparative Genomics Reveals Three Genetic Groups of the Whitefly Obligate Endosymbiont Candidatus Portiera aleyrodidarum

Maternally inherited obligate endosymbionts codiverge with their invertebrate hosts and reflect their host's evolutionary history. Whiteflies (Hemiptera: Aleyrodidae) harbor one obligate endosymbiont, Candidatus Portiera aleyrodidarum (hereafter Portiera). Portiera was anciently acquired by whitefly and has been coevolving with its host ever since. Uncovering the divergence of endosymbionts provides a fundamental basis for inspecting the coevolutionary processes between the bacteria and their hosts. To illustrate the divergence of Portiera lineages across different whitefly species, we sequenced the Portiera genome from Aleyrodes shizuokensis and conducted a comparative analysis on the basic features and gene evolution with bacterial genomes from five whitefly genera, namely Aleurodicus, Aleyrodes, Bemisia, Pealius, and Trialeurodes. The results indicated that Portiera from Bemisia possessed significantly larger genomes, fewer coding sequences (CDSs), and a lower coding density. Their gene arrangement differed notably from those of other genera. The phylogeny of the nine Portiera lineages resembled that of their hosts. Moreover, the lineages were classified into three distinct genetic groups based on the genetic distance, one from Aleurodicus (Aleurodicinae), one from Bemisia (Aleyrodinae), and another from Aleyrodes, Pealius, and Trialeurrodes (Aleyrodinae). Synonymous and nonsynonymous rate analyses, parity rule 2 plot analyses, neutrality plot analyses, and effective number of codons analyses supported the distinction of the three genetic groups. Our results indicated that Portiera from distant hosts exhibit distinct genomic contents, implying codivergence between hosts and their endosymbionts. This work will enhance our understanding of coevolution between hosts and their endosymbionts.

Phylogenetic analysis of Shewanella spp. isolated from fish

Different species of Shewanella spp. widely inhabit freshwater and marine environments. Some of them are opportunistic fish pathogens. The application of high-throughput sequencing enabled the characterization and taxonomic reclassification of many Shewanella spp. species. Still, some strains collected from fish need to be better recognized. The aim of the present study was to classify and determine the phylogenetic relationships of Shewanella spp. collected from fish. The complete genomes of 94 strains of Shewanella spp. from different fish species were sequenced using Illumina platform (MiSeq). The 16S rRNA gene, genomic features and whole-genome relationships of those bacteria were comprehensively analysed in comparison to reference strains. Whole-genome analysis showed that the tested Shewanella spp. strains were clustered into six groups similar to reference strains of S. xiamenensis, S. oneidensis, S. glacialipiscicola, S. hafniensis, S. baltica and S. oncorhynchi. Our study indicates that the whole-genome sequence analysis enabled taxonomic classification and assessment of the diversity of the Shewanella spp. strains, as opposed to recently the gold standard method of 16S rRNA amplicon sequencing. The high genetic diversity and low similarity to the reference genome of S. oneidensis indicate that the group of strains may be a subspecies or even new species. Furthermore, we showed that the most frequent Shewanella spp. species occurring in freshwater fish in our study is the recently described species S. oncorhynchi.

Genomic Analysis of Group B Streptococcus Carriage Isolates From Botswana Reveals Distinct Local Epidemiology and Identifies Novel Strains

In pregnant people colonized with group B Streptococcus (GBS) in Botswana, we report the presence/expansion of sequence types 223 and 109, a low rate of erythromycin resistance, and 3 novel sequence types. These data highlight the importance of local epidemiologic studies of GBS, a significant source of neonatal disease.

simpiTB - a pipeline designed to extract meaningful information from whole genome sequencing data of Mycobacterium tuberculosis complex, allows to combine genomic, phylogenetic and clustering analyses in existing SITVIT databases

Data obtained from new sequencing technologies are evolving rapidly, leading to the development of specific bioinformatic tools, pipelines and softwares. Several algorithms and tools are today available allowing a better identification and description of Mycobacterium tuberculosis complex (MTBC) isolates worldwide. Our approach consists in applying existing methods to analyze DNA sequencing data (from FASTA or FASTQ files), and tentatively extract meaningful information that would facilitate identification as well as a better understanding and management of MTBC isolates (taking into account whole genome sequencing and classical genotyping data). The aim of this study is to propose a pipeline analysis allowing to potentially simplify MTBC data analysis by providing different ways to interpret genomic or genotyping information based on existing tools. Furthermore, we propose a "reconciledTB" list making a link with results directly obtained from whole genome sequencing (WGS) data and results obtained from classical genotyping analysis (data inferred from SpoTyping and MIRUReader). Data visualization graphics and trees generated provide additional elements to better understand and confer associations among information overlap analyses. Additionally, comparison between data entered in an international genotyping database (SITVITEXTEND) and ensuing data obtained from the pipeline not only provide meaningful information, but further suggest that simpiTB could also be suitable for new data integration in specific TB genotyping databases.

Complete De Novo Assembly of Wolbachia Endosymbiont of Frankliniella intonsa

As an endosymbiont, Wolbachia exerts significant effects on the host, including on reproduction, immunity, and metabolism. However, the study of Wolbachia in Thysanopteran insects, such as flower thrips Frankliniella intonsa, remains limited. Here, we assembled a gap-free looped genome assembly of Wolbachia strain wFI in a length of 1,463,884 bp (GC content 33.80%), using Nanopore long reads and Illumina short reads. The annotation of wFI identified a total of 1838 protein-coding genes (including 85 pseudogenes), 3 ribosomal RNAs (rRNAs), 35 transfer RNAs (tRNAs), and 1 transfer-messenger RNA (tmRNA). Beyond this basic description, we identified mobile genetic elements, such as prophage and insertion sequences (ISs), which make up 17% of the entire wFI genome, as well as genes involved in riboflavin and biotin synthesis and metabolism. This research lays the foundation for understanding the nutritional mutualism between Wolbachia and flower thrips. It also serves as a valuable resource for future studies delving into the intricate interactions between Wolbachia and its host.

Genomic dissection of endemic carbapenem resistance reveals metallo-beta-lactamase dissemination through clonal, plasmid and integron transfer

Infections caused by metallo-beta-lactamase-producing organisms (MBLs) are a global health threat. Our understanding of transmission dynamics and how MBLs establish endemicity remains limited. We analysed two decades of blaIMP-4 evolution in a hospital using sequence data from 270 clinical and environmental isolates (including 169 completed genomes) and identified the blaIMP-4 gene across 7 Gram-negative genera, 68 bacterial strains and 7 distinct plasmid types. We showed how an initial multi-species outbreak of conserved IncC plasmids (95 genomes across 37 strains) allowed endemicity to be established through the ability of blaIMP-4 to disseminate in successful strain-genetic setting pairs we termed propagators, in particular Serratia marcescens and Enterobacter hormaechei. From this reservoir, blaIMP-4 persisted through diversification of genetic settings that resulted from transfer of blaIMP-4 plasmids between bacterial hosts and of the integron carrying blaIMP-4 between plasmids. Our findings provide a framework for understanding endemicity and spread of MBLs and may have broader applicability to other carbapenemase-producing organisms.

The conserved iol gene cluster in Pseudomonas is involved in rhizosphere competence

The Pseudomonas genus has shown great potential as a sustainable solution to support agriculture through its plant-growth-promoting and biocontrol activities. However, their efficacy as bioinoculants is limited by unpredictable colonization in natural conditions. Our study identifies the iol locus, a gene cluster in Pseudomonas involved in inositol catabolism, as a feature enriched among superior root colonizers in natural soil. Further characterization revealed that the iol locus increases competitiveness, potentially caused by an observed induction of swimming motility and the production of fluorescent siderophore in response to inositol, a plant-derived compound. Public data analyses indicate that the iol locus is broadly conserved in the Pseudomonas genus and linked to diverse host-microbe interactions. Together, our findings suggest the iol locus as a potential target for developing more effective bioinoculants for sustainable agriculture.

TheiaEuk: a species-agnostic bioinformatics workflow for fungal genomic characterization

Introduction

The clinical incidence of antimicrobial-resistant fungal infections has dramatically increased in recent years. Certain fungal pathogens colonize various body cavities, leading to life-threatening bloodstream infections. However, the identification and characterization of fungal isolates in laboratories remain a significant diagnostic challenge in medicine and public health. Whole-genome sequencing provides an unbiased and uniform identification pipeline for fungal pathogens but most bioinformatic analysis pipelines focus on prokaryotic species. To this end, TheiaEuk_Illumina_PE_PHB (TheiaEuk) was designed to focus on genomic analysis specialized to fungal pathogens.

Methods

TheiaEuk was designed using containerized components and written in the workflow description language (WDL) to facilitate deployment on the cloud-based open bioinformatics platform Terra. This species-agnostic workflow enables the analysis of fungal genomes without requiring coding, thereby reducing the entry barrier for laboratory scientists. To demonstrate the usefulness of this pipeline, an ongoing outbreak of C. auris in southern Nevada was investigated. We performed whole-genome sequence analysis of 752 new C. auris isolates from this outbreak. Furthermore, TheiaEuk was utilized to observe the accumulation of mutations in the FKS1 gene over the course of the outbreak, highlighting the utility of TheiaEuk as a monitor of emerging public health threats when combined with whole-genome sequencing surveillance of fungal pathogens.

Results

A primary result of this work is a curated fungal database containing 5,667 unique genomes representing 245 species. TheiaEuk also incorporates taxon-specific submodules for specific species, including clade-typing for Candida auris (C. auris). In addition, for several fungal species, it performs dynamic reference genome selection and variant calling, reporting mutations found in genes currently associated with antifungal resistance (FKS1, ERG11, FUR1). Using genome assemblies from the ATCC Mycology collection, the taxonomic identification module used by TheiaEuk correctly assigned genomes to the species level in 126/135 (93.3%) instances and to the genus level in 131/135 (97%) of instances, and provided zero false calls. Application of TheiaEuk to actual specimens obtained in the course of work at a local public health laboratory resulted in 13/15 (86.7%) correct calls at the species level, with 2/15 called at the genus level. It made zero incorrect calls. TheiaEuk accurately assessed clade type of Candida auris in 297/302 (98.3%) of instances.

Discussion

TheiaEuk demonstrated effectiveness in identifying fungal species from whole genome sequence. It further showed accuracy in both clade-typing of C. auris and in the identification of mutations known to associate with drug resistance in that organism.

The distribution of fitness effects of plasmid pOXA-48 in clinical enterobacteria

Antimicrobial resistance (AMR) in bacteria is a major public health problem. The main route for AMR acquisition in clinically important bacteria is the horizontal transfer of plasmids carrying resistance genes. AMR plasmids allow bacteria to survive antibiotics, but they also entail physiological alterations in the host cell. Multiple studies over the last few years have indicated that these alterations can translate into a fitness cost when antibiotics are absent. However, due to technical limitations, most of these studies are based on analysing new associations between plasmids and bacteria generated in vitro, and we know very little about the effects of plasmids in their native bacterial hosts. In this study, we used a CRISPR-Cas9-tool to selectively cure plasmids from clinical enterobacteria to overcome this limitation. Using this approach, we were able to study the fitness effects of the carbapenem resistance plasmid pOXA-48 in 35 pOXA-48-carrying isolates recovered from hospitalized patients. Our results revealed that pOXA-48 produces variable effects across the collection of wild-type enterobacterial strains naturally carrying the plasmid, ranging from fitness costs to fitness benefits. Importantly, the plasmid was only associated with a significant fitness reduction in four out of 35 clones, and produced no significant changes in fitness in the great majority of isolates. Our results suggest that plasmids produce neutral fitness effects in most native bacterial hosts, helping to explain the great prevalence of plasmids in natural microbial communities.

Particular genomic and virulence traits associated with preterm infant-derived toxigenic Clostridium perfringens strains

Clostridium perfringens is an anaerobic toxin-producing bacterium associated with intestinal diseases, particularly in neonatal humans and animals. Infant gut microbiome studies have recently indicated a link between C. perfringens and the preterm infant disease necrotizing enterocolitis (NEC), with specific NEC cases associated with overabundant C. perfringens termed C. perfringens-associated NEC (CPA-NEC). In the present study, we carried out whole-genome sequencing of 272 C. perfringens isolates from 70 infants across 5 hospitals in the United Kingdom. In this retrospective analysis, we performed in-depth genomic analyses (virulence profiling, strain tracking and plasmid analysis) and experimentally characterized pathogenic traits of 31 strains, including 4 from CPA-NEC patients. We found that the gene encoding toxin perfringolysin O, pfoA, was largely deficient in a human-derived hypovirulent lineage, as well as certain colonization factors, in contrast to typical pfoA-encoding virulent lineages. We determined that infant-associated pfoA+ strains caused significantly more cellular damage than pfoA- strains in vitro, and further confirmed this virulence trait in vivo using an oral-challenge C57BL/6 murine model. These findings suggest both the importance of pfoA+ C. perfringens as a gut pathogen in preterm infants and areas for further investigation, including potential intervention and therapeutic strategies.

Genetic and enzymatic characterization of two novel blaNDM-36, -37 variants in Escherichia coli strains

The widespread of different NDM variants in clinical Enterobacterales isolates poses a serious public health concern, which requires continuous monitoring. In this study, three E. coli strains carrying two novel blaNDM variants of blaNDM-36, -37 were identified from a patient with refractory urinary tract infection (UTI) in China. We conducted antimicrobial susceptibility testing (AST), enzyme kinetics analysis, conjugation experiment, whole-genome sequencing (WGS), and bioinformatics analysis to characterize the blaNDM-36, -37 enzymes and their carrying strains. The blaNDM-36, -37 harboring E. coli isolates belonged to ST227, O9:H10 serotype and exhibited intermediate or resistance to all β-lactams tested except aztreonam and aztreonam/avibactam. The genes of blaNDM-36, -37 were located on a conjugative IncHI2-type plasmid. NDM-37 differed from NDM-5 by a single amino acid substitution (His261Tyr). NDM-36 differed from NDM-37 by an additional missense mutation (Ala233Val). NDM-36 had increased hydrolytic activity toward ampicillin and cefotaxime relative to NDM-37 and NDM-5, while NDM-37 and NDM-36 had lower catalytic activity toward imipenem but higher activity against meropenem in comparison to NDM-5. This is the first report of co-occurrence of two novel blaNDM variants in E. coli isolated from the same patient. The work provides insights into the enzymatic function and demonstrates the ongoing evolution of NDM enzymes.

Integrative omics identifies conserved and pathogen-specific responses of sepsis-causing bacteria

Even in the setting of optimal resuscitation in high-income countries severe sepsis and septic shock have a mortality of 20-40%, with antibiotic resistance dramatically increasing this mortality risk. To develop a reference dataset enabling the identification of common bacterial targets for therapeutic intervention, we applied a standardized genomic, transcriptomic, proteomic and metabolomic technological framework to multiple clinical isolates of four sepsis-causing pathogens: Escherichia coli, Klebsiella pneumoniae species complex, Staphylococcus aureus and Streptococcus pyogenes. Exposure to human serum generated a sepsis molecular signature containing global increases in fatty acid and lipid biosynthesis and metabolism, consistent with cell envelope remodelling and nutrient adaptation for osmoprotection. In addition, acquisition of cholesterol was identified across the bacterial species. This detailed reference dataset has been established as an open resource to support discovery and translational research.

Genomic profiling of dioecious Amaranthus species provides novel insights into species relatedness and sex genes

BACKGROUND

Amaranthus L. is a diverse genus consisting of domesticated, weedy, and non-invasive species distributed around the world. Nine species are dioecious, of which Amaranthus palmeri S. Watson and Amaranthus tuberculatus (Moq.) J.D. Sauer are troublesome weeds of agronomic crops in the USA and elsewhere. Shallow relationships among the dioecious Amaranthus species and the conservation of candidate genes within previously identified A. palmeri and A. tuberculatus male-specific regions of the Y (MSYs) in other dioecious species are poorly understood. In this study, seven genomes of dioecious amaranths were obtained by paired-end short-read sequencing and combined with short reads of seventeen species in the family Amaranthaceae from NCBI database. The species were phylogenomically analyzed to understand their relatedness. Genome characteristics for the dioecious species were evaluated and coverage analysis was used to investigate the conservation of sequences within the MSY regions.

RESULTS

We provide genome size, heterozygosity, and ploidy level inference for seven newly sequenced dioecious Amaranthus species and two additional dioecious species from the NCBI database. We report a pattern of transposable element proliferation in the species, in which seven species had more Ty3 elements than copia elements while A. palmeri and A. watsonii had more copia elements than Ty3 elements, similar to the TE pattern in some monoecious amaranths. Using a Mash-based phylogenomic analysis, we accurately recovered taxonomic relationships among the dioecious Amaranthus species that were previously identified based on comparative morphology. Coverage analysis revealed eleven candidate gene models within the A. palmeri MSY region with male-enriched coverages, as well as regions on scaffold 19 with female-enriched coverage, based on A. watsonii read alignments. A previously reported FLOWERING LOCUS T (FT) within A. tuberculatus MSY contig was also found to exhibit male-enriched coverages for three species closely related to A. tuberculatus but not for A. watsonii reads. Additional characterization of the A. palmeri MSY region revealed that 78% of the region is made of repetitive elements, typical of a sex determination region with reduced recombination.

CONCLUSIONS

The results of this study further increase our understanding of the relationships among the dioecious species of the Amaranthus genus as well as revealed genes with potential roles in sex function in the species.

Contribution bias of parental genomes to the hybrid lineages of black Amur bream and topmouth culter revealed by low-coverage whole-genome sequencing

Interspecific hybridization has the potential to increase animal genetic diversity, enable the introgression of advantageous genetic variation, and even create hybrid lineages. Using the low-coverage whole-genome sequencing data of black Amur bream (Megalobrama terminalis [SJ]), topmouth culter (Culter alburnus [QZ]), and their hybrid progeny, we focused on the sequence variation and contribution bias to the genomic composition of the hybrid lineages of SJ♀ × QZ♂ [ZJ] and QZ♀ × SJ♂ [FJ]. Both mash-based genetic distance analysis and variant analysis based on various reference genomes revealed that the hybrid lineages' genomes are more prone to QZ composition. Furthermore, we assess the degree of subgenomic intermixing between parents QZ and SJ in the FJ and ZJ hybrid lineages. In the Jaccard index-zero genomic area, 89.41% of the genomic region in ZJ and 90.23% in FJ are entirely from one of the parental genomes. This genomic area is primarily from QZ, regardless of FJ or ZJ. In FJ, the proportion is 55.70%; in ZJ, the proportion is 54.60%. Using the reference-free variant discovery tool, we identified two InDel markers (Both were SSR-type InDel markers) and one (T/A)-type SNP marker that displayed differences in agarose electrophoresis. Our current findings suggest that parental genomic contribution bias leads to an imbalance in the genomic composition of hybrid lineages. The genomic composition bias is connected to species but not maternal effects, giving insights into the development of vertebrate genomes shortly after hybridization.

PorinPredict: In Silico Identification of OprD Loss from WGS Data for Improved Genotype-Phenotype Predictions of P. aeruginosa Carbapenem Resistance

The increasing integration of genomics into routine clinical diagnostics requires reliable computational tools to identify determinants of antimicrobial resistance (AMR) from whole-genome sequencing data. Here, we developed PorinPredict, a bioinformatic tool that predicts defects of the Pseudomonas aeruginosa outer membrane porin OprD, which are strongly associated with reduced carbapenem susceptibility. PorinPredict relies on a database of intact OprD variants and reports inactivating mutations in the coding or promoter region. PorinPredict was validated against 987 carbapenemase-negative P. aeruginosa genomes, of which OprD loss was predicted for 454 out of 522 (87.0%) meropenem-nonsusceptible and 46 out of 465 (9.9%) meropenem-susceptible isolates. OprD loss was also found to be common among carbapenemase-producing isolates, resulting in even further increased MICs. Chromosomal mutations in quinolone resistance-determining regions and OprD loss commonly co-occurred, likely reflecting the restricted use of carbapenems for multidrug-resistant infections as recommended in antimicrobial stewardship programs. In combination with available AMR gene detection tools, PorinPredict provides a robust and standardized approach to link P. aeruginosa phenotypes to genotypes. IMPORTANCE Pseudomonas aeruginosa is a major cause of multidrug-resistant nosocomial infections. The emergence and spread of clones exhibiting resistance to carbapenems, a class of critical last-line antibiotics, is therefore closely monitored. Carbapenem resistance is frequently mediated by chromosomal mutations that lead to a defective outer membrane porin OprD. Here, we determined the genetic diversity of OprD variants across the P. aeruginosa population and developed PorinPredict, a bioinformatic tool that enables the prediction of OprD loss from whole-genome sequencing data. We show a high correlation between predicted OprD loss and meropenem nonsusceptibility irrespective of the presence of carbapenemases, which are a second widespread determinant of carbapenem resistance. Isolates with resistance determinants to other antibiotics were disproportionally affected by OprD loss, possibly due to an increased exposure to carbapenems. Integration of PorinPredict into genomic surveillance platforms will facilitate a better understanding of the clinical impact of OprD modifications and transmission dynamics of resistant clones.

Antifungal susceptibility and molecular characteristics of Cryptococcus spp. based on whole-genome sequencing in Zhejiang Province, China

Cryptococcus spp. is a complex species that often causes cryptococcosis, which is one of the most common opportunistic infections in adults living with HIV and has very high morbidity and mortality rates. This study aimed to investigate the antifungal susceptibility profiles and epidemiological characteristics of the Cryptococcus neoformans species complex (CNSC) and the Cryptococcus gattii species complex (CGSC) in Zhejiang Province, China. A total of 177 CNSC and 3 CGSC isolates were collected, and antifungal susceptibility was tested by FUNGUS 3 and verified with an E-test. Moreover, multiple classification methods and genomic analyses were performed. The majority of the isolates (96.11%) were C. neoformans (formerly C. neoformans var. grubii) (ST5-VNI-A-α). Our study highlights that most of the patients with cryptococcosis were non-HIV patients in China, and nearly half of them did not have underlying diseases that led to immune insufficiency. Most of the Cryptococcus spp. isolates in this study were sensitive to common antifungal drugs. Two 5-flucytosine (5-FC)-resistant strains were identified, and FUR1 mutation was detected in the 5-FC-resistant isolates. Typing based on whole-genome sequencing (WGS) showed better discrimination than that achieved with multilocus sequence typing (MLST) and indicated a clear population structure. A phylogenetic analysis based on WGS included more genomic information than traditional classification methods.

The hazard of carbapenemase (OXA-181)-producing Escherichia coli spreading in pig and veal calf holdings in Italy in the genomics era: Risk of spill over and spill back between humans and animals

Carbapenemase-producing Enterobacterales (CPE) are considered a major public health issue. In the frame of the EU Harmonized AMR Monitoring program conducted in Italy in 2021, 21 epidemiological units of fattening pigs (6.98%; 95% CI 4.37-10.47%; 21/301) and four epidemiological units of bovines <12 months (1.29%; 95% CI 0.35-3.27%, 4/310) resulted positive to OXA-48-like-producing E. coli (n = 24 OXA-181, n = 1 OXA-48). Whole Genome Sequencing (WGS) for in-depth characterization, genomics and cluster analysis of OXA-181-(and one OXA-48) producing E. coli isolated, was performed. Tracing-back activities at: (a) the fattening holding of origin of one positive slaughter batch, (b) the breeding holding, and (c) one epidemiologically related dairy cattle holding, allowed detection of OXA-48-like-producing E. coli in different units and comparison of further human isolates from fecal samples of farm workers. The OXA-181-producing isolates were multidrug resistant (MDR), belonged to different Sequence Types (STs), harbored the IncX and IncF plasmid replicons and multiple virulence genes. Bioinformatics analysis of combined Oxford Nanopore Technologies (ONT) long reads and Illumina short reads identified bla OXA-181 as part of a transposon in IncX1, IncX3, and IncFII fully resolved plasmids from 16 selected E. coli, mostly belonging to ST5229, isolated during the survey at slaughter and tracing-back activities. Although human source could be the most likely cause for the introduction of the bla OXA-181-carrying IncX1 plasmid in the breeding holding, concerns arise from carbapenemase OXA-48-like-producing E. coli spreading in 2021 in Italian fattening pigs and, to a lesser extent, in veal calf holdings.

Plasmidome in mcr-1 harboring carbapenem-resistant enterobacterales isolates from human in Thailand

The emergence of the mobile colistin-resistance genes mcr-1 has attracted significant attention worldwide. This study aimed to investigate the genetic features of mcr-1-carrying plasmid among carbapenem-resistant Enterobacterales (CRE) isolates and the potential genetic basis governing transmission. Seventeen mcr-harboring isolates were analyzed based on whole genome sequencing using short-read and long-read platforms. All the mcr-1-carrying isolates could be conjugatively transferred into a recipient Escherichia coli UB1637. Among these 17 isolates, mcr-1 was located on diverse plasmid Inc types, consisting of IncX4 (11/17; 64.7%), IncI2 (4/17; 23.53%), and IncHI/IncN (2/17; 11.76%). Each of these exhibited remarkable similarity in the backbone set that is responsible for plasmid replication, maintenance, and transfer, with differences being in the upstream and downstream regions containing mcr-1. The IncHI/IncN type also carried other resistance genes (bla_TEM-1B or bla_TEM-135). The mcr-1-harboring IncX4 plasmids were carried in E. coli ST410 (7/11; 63.6%) and ST10 (1/11; 9.1%) and Klebsiella pneumoniae ST15 (1/11; 9.1%), ST336 (1/11; 9.1%), and ST340 (1/11; 9.1%). The IncI2-type plasmid was harbored in E. coli ST3052 (1/4; 25%) and ST1287 (1/4; 25%) and in K. pneumoniae ST336 (2/4; 50%), whereas IncHI/IncN were carried in E. coli ST6721 (1/2; 50%) and new ST (1/2; 50%). The diverse promiscuous plasmids may facilitate the spread of mcr-1 among commensal E. coli or K. pneumoniae strains in patients. These results can provide information for a surveillance system and infection control for dynamic tracing.

Genetic Structure of the Liriope muscari Polyploid Complex and the Possibility of Its Genetic Disturbance in Japan

Anthropogenic activities, such as the movement of plants through greening, can result in genetic disturbance that can interfere with local adaptation in wild populations. Although research is underway to prevent genetic disturbance associated with greening, genetic disturbance of intraspecific polyploidy, which is estimated to be present in 24% of vascular plants, has not been well studied. Liriope muscari is a polyploid complex with known diploid (2n = 36), tetraploid (2n = 72), and hexaploid (2n = 108) forms. The plants of this species tolerate dry and hot conditions and are therefore frequently used for greening and gardening. However, the distribution of this polyploid in Japan, its genetic structure, and genetic disturbance are not known. In this study, we investigated the polyploidy distribution and genetic structure in naturally distributed L. muscari in Japan using chloroplast DNA (cpDNA) haplotypes and nuclear DNA (nDNA). Commercially produced individuals were also studied and compared with natural populations to assess any genetic disturbance of the ploidy complex in this species. Chromosome counts, cpDNA, and nDNA results showed three genetically and cytologically distinct groups in Japan: first, a tetraploid group in mainland Japan; second, a hexaploid group in the Ryukyu Islands; and third, a diploid and tetraploid group in the Ryukyu Islands. Significant isolation by distance was also detected within the three groups (p = 0.001). Genetic disturbance due to greening and gardening should be avoided among the three groups. Genetic disturbance can be reduced by using individuals derived from natural populations that are close to the sites used for greening and gardening. For commercially produced individuals, genetic disturbance is unlikely in the Kanto region, an area of high usage, while genetic disturbance is thought possible in the Ryukyu Islands.

Membrane Clustering of Coronavirus Variants Using Document Similarity

Currently, as an effect of the COVID-19 pandemic, bioinformatics, genomics, and biological computations are gaining increased attention. Genomes of viruses can be represented by character strings based on their nucleobases. Document similarity metrics can be applied to these strings to measure their similarities. Clustering algorithms can be applied to the results of their document similarities to cluster them. P systems or membrane systems are computation models inspired by the flow of information in the membrane cells. These can be used for various purposes, one of them being data clustering. This paper studies a novel and versatile clustering method for genomes and the utilization of such membrane clustering models using document similarity metrics, which is not yet a well-studied use of membrane clustering models.

Global diversity and distribution of prophages are lineage-specific within the Ralstonia solanacearum species complex

Background

Ralstonia solanacearum species complex (RSSC) strains are destructive plant pathogenic bacteria and the causative agents of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, their virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape RSSC virulence, the global diversity, distribution, and wider functional gene content of RSSC prophages are unknown. In this study, prophages were identified in a diverse collection of 192 RSSC draft genome assemblies originating from six continents.

Results

Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distributions were characterised using metadata on RSSC strain geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 313 intact prophages were identified, forming ten genetically distinct clusters. These included six prophage clusters with similarity to the Inoviridae, Myoviridae, and Siphoviridae phage families, and four uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distributions, being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained many auxiliary genes involved in metabolism and virulence of both phage and bacteria.

Conclusions

Our results show that while RSSC prophages are highly diverse globally, they make lineage-specific contributions to the RSSC accessory genome, which could have resulted from shared coevolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08909-7.

Genomic epidemiology of Mycobacterium abscessus in a Canadian cystic fibrosis centre

The Mycobacterium abscessus complex causes significant morbidity and mortality among patients with Cystic Fibrosis (CF). It has been hypothesized that these organisms are transmitted from patient to patient based on genomics. However, few studies incorporate epidemiologic data to confirm this hypothesis. We longitudinally sampled 27 CF and 7 non-CF patients attending a metropolitan hospital in Ontario, Canada from 2013 to 2018. Whole genome sequencing along with epidemiological data was used to evaluate the likelihood of transmission. Overall, the genetic diversity of M. abscessus was large, with a median pairwise distance (IQR) of 1,279 (143-134) SNVs between all Ontario M. abscessus isolates and 2,908 (21-3,204) single nucleotide variants (SNVs) between M. massiliense isolates. This reflects the global diversity of this pathogen, with Ontario isolates widely dispersed throughout global phylogenetic trees of each subspecies. Using a maximum distance of 25 SNVs as a threshold to identify possible transmission, we identified 23 (of 276 total) pairs of closely-related isolates. However, transmission was probable for only one pair based on both genomic and epidemiological data. This suggests that person-to-person transmission of M. abscessus among CF patients is indeed rare and reinforces the critical importance of epidemiological data for inferences of transmission.

Delineating Mycobacterium abscessus population structure and transmission employing high-resolution core genome multilocus sequence typing

Mycobacterium abscessus is an emerging multidrug-resistant non-tuberculous mycobacterium that causes a wide spectrum of infections and has caused several local outbreaks worldwide. To facilitate standardized prospective molecular surveillance, we established a novel core genome multilocus sequence typing (cgMLST) scheme. Whole genome sequencing data of 1991 isolates were employed to validate the scheme, re-analyze global population structure and set genetic distance thresholds for cluster detection and taxonomic identification. We confirmed and amended the nomenclature of the main dominant circulating clones and found that these also correlate well with traditional 7-loci MLST. Dominant circulating clones could be linked to a corresponding reference genome with less than 250 alleles while 99% of pairwise comparisons between epidemiologically linked isolates were below 25 alleles and 90% below 10 alleles. These thresholds can be used to guide further epidemiological investigations. Overall, the scheme will help to unravel the apparent global spread of certain clonal complexes and as yet undiscovered transmission routes.

Convergence of virulence and antimicrobial resistance in increasingly prevalent Escherichia coli ST131 papGII+ sublineages

Escherichia coli lineage ST131 is an important cause of urinary tract and bloodstream infections worldwide and is highly resistant to antimicrobials. Specific ST131 lineages carrying invasiveness-associated papGII pathogenicity islands (PAIs) were previously described, but it is unknown how invasiveness relates to the acquisition of antimicrobial resistance (AMR). In this study, we analysed 1638 ST131 genomes and found that papGII+ isolates carry significantly more AMR genes than papGII-negative isolates, suggesting a convergence of virulence and AMR. The prevalence of papGII+ isolates among human clinical ST131 isolates increased dramatically since 2005, accounting for half of the recent E. coli bloodstream isolates. Emerging papGII+ lineages within clade C2 were characterized by a chromosomally integrated blaCTX-M-15 and the loss and replacement of F2:A1:B- plasmids. Convergence of virulence and AMR is worrying, and further dissemination of papGII+ ST131 lineages may lead to a rise in severe and difficult-to-treat extraintestinal infections.

Genomic analyses indicate that the ST131 lineage of E. coli, responsible for urinary and bloodstream infections globally, is evolving towards both increased virulence and increased resistance to antimicrobials.

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.