# The evolution and international spread of extensively drug resistant Shigella sonnei

Lewis Mason, David Greig, Lauren Cowley, Sally Partridge, Elena Martinez, Grace Blackwell, Charlotte Chong, P. Malaka de Silva, Rebecca Bengtsson, Jenny Draper, et al.

► To cite this version:

Lewis Mason, David Greig, Lauren Cowley, Sally Partridge, Elena Martinez, et al.. The evolution and international spread of extensively drug resistant Shigella sonnei. Nature Communications, 2023, 14, pp.1983. 10.1038/s41467-023-37672-w. [https://hal.science/hal-04613474](https://hal.science/hal-04613474)

HAL Id: hal-04613474

https://hal.science/hal-04613474

Submitted on 17 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# The evolution and international spread of extensively drug resistant Shigella sonnei

Lewis Mason, David Greig, Lauren Cowley, Sally Partridge, Elena Martinez, Grace Blackwell, Charlotte Chong, P. Malaka de Silva, Rebecca Bengtsson, Jenny Draper, et al.

► To cite this version:

Lewis Mason, David Greig, Lauren Cowley, Sally Partridge, Elena Martinez, et al.. The evolution and international spread of extensively drug resistant Shigella sonnei. Nature Communications, 2023, 14, pp.1983. 10.1038/s41467-023-37672-w. [https://hal.science/hal-04613474v1](https://hal.science/hal-04613474v1)

HAL Id: hal-04613474

Submitted on 17 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# The evolution and international spread of extensively drug resistant *Shigella sonnei*

Received: 12 September 2022

Accepted: 24 March 2023

Published online: 08 April 2023

<figure><img src="image_3.png" /></figure>

Check for updates

Lewis C. E. Mason®<fnref n="1" />,², David R. Greig³, Lauren A. Cowley<fnref n="4" />, Sally R. Partridge®⁵,⁶,⁷,⁸, Elena Martinez⁷,⁹, Grace A. Blackwell⁷,⁹, Charlotte E. Chong², P. Malaka De Silva², Rebecca J. Bengtsson², Jenny L. Draper®⁷,⁹, Andrew N. Ginn⁷,⁸,⁹,¹⁰, Indy Sandaradura®⁶,⁷,⁹, Eby M. Sim⁵,⁶,⁸, Jonathan R. Iredell⁵,⁶,⁷,⁸, Vitali Sintchenko®⁵,⁶,⁷,⁸,⁹,¹¹, Danielle J. Ingle®¹², Benjamin P. Howden®¹², Sophie Lefèvre®¹³, Elisabeth Njamkepo®¹³, François-Xavier Weill®¹³, Pieter-Jan Ceyssens¹⁴, Claire Jenkins³ & Kate S. Baker®<fnref n="1" />,²✉

*Shigella sonnei* causes shigellosis, a severe gastrointestinal illness that is sexually transmissible among men who have sex with men (MSM). Multidrug resistance in *S. sonnei* is common including against World Health Organisation recommended treatment options, azithromycin, and ciprofloxacin. Recently, an MSM-associated outbreak of extended-spectrum β-lactamase producing, extensively drug resistant *S. sonnei* was reported in the United Kingdom. Here, we aimed to identify the genetic basis, evolutionary history, and international dissemination of the outbreak strain. Our genomic epidemiological analyses of 3,304 isolates from the United Kingdom, Australia, Belgium, France, and the United States of America revealed an internationally connected outbreak with a most recent common ancestor in 2018 carrying a low-fitness cost resistance plasmid, previously observed in travel associated sublineages of *S. flexneri*. Our results highlight the persistent threat of horizontally transmitted antimicrobial resistance and the value of continuing to work towards early and open international sharing of genomic surveillance data.

Bacteria belonging to the genus *Shigella* cause shigellosis, a gastrointestinal infection characterised by a combination of profuse, often bloody, diarrhoea with abdominal pain, fever, nausea and/or vomiting<fnref n="1" />. In middle- and higher-income nations, although symptoms

can be severe, the duration of illness is usually short-lived². However, in lower to middle income countries (LMIC), shigellosis is associated with higher levels of morbidity and mortality³. *Shigella* infection has no substantial non-human reservoir and is transmitted from person to

<footnote n="1">¹NIHR HPRU in Gastrointestinal Infections at University of Liverpool, Liverpool, UK. ²Department of Clinical Infection, Microbiology, and Immunology; Institute for Infection, Veterinary and Ecological Sciences, Liverpool, UK. ³Gastro and Food Safety (One Health) Division, UK Health Security Agency, London, UK.</footnote>

<footnote n="4">⁴Milner Centre for Evolution, University of Bath, Bath, UK. ⁵Centre for Infectious Diseases and Microbiology, The Westmead Institute for Medical Research, Westmead, NSW, Australia. ⁶Western Sydney Local Health District, Westmead, NSW, Australia. ⁷Faculty of Medicine and Health, University of Sydney, Sydney, NSW, Australia. ⁸Sydney Infectious Diseases Institute, University of Sydney, Sydney, NSW, Australia. ⁹New South Wales Health Pathology, Dee Why, NSW, Australia. ¹⁰Douglass Hanly Moir Pathology, Macquarie Park, NSW, Australia. ¹¹Centre for Infectious Diseases and Microbiology – Public Health, Institute for Clinical Pathology and Microbiology Research, Westmead Hospital, Westmead, NSW, Australia. ¹²Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia. ¹³Institut Pasteur, Université Paris Cité, Unité des Bactéries pathogènes entériques, Centre National de Référence des Escherichia coli, Shigella et Salmonella, Paris F-75015, France. ¹⁴Division of Human Bacterial Diseases, Sciensano, Ixelles, Belgium. ✉e-mail: kbaker@liverpool.ac.uk</footnote>
<figure><img src="image_2.png" /></figure>

# The evolution and international spread of extensively drug resistant *Shigella sonnei*

Received: 12 September 2022

Accepted: 24 March 2023

Published online: 08 April 2023

<figure><img src="image_3.png" /></figure>

Check for updates

Lewis C. E. Mason®<fnref n="1" />,², David R. Greig³, Lauren A. Cowley<fnref n="4" />, Sally R. Partridge®⁵,⁶,⁷,⁸, Elena Martinez⁷,⁹, Grace A. Blackwell⁷,⁹, Charlotte E. Chong², P. Malaka De Silva², Rebecca J. Bengtsson², Jenny L. Draper®⁷,⁹, Andrew N. Ginn⁷,⁸,⁹,¹⁰, Indy Sandaradura®⁶,⁷,⁹, Eby M. Sim⁵,⁶,⁸, Jonathan R. Iredell⁵,⁶,⁷,⁸, Vitali Sintchenko®⁵,⁶,⁷,⁸,⁹,¹¹, Danielle J. Ingle®¹², Benjamin P. Howden®¹², Sophie Lefèvre®¹³, Elisabeth Njamkepo®¹³, François-Xavier Weill®¹³, Pieter-Jan Ceyssens¹⁴, Claire Jenkins³ & Kate S. Baker®<fnref n="1" />,²✉

*Shigella sonnei* causes shigellosis, a severe gastrointestinal illness that is sexually transmissible among men who have sex with men (MSM). Multidrug resistance in *S. sonnei* is common including against World Health Organisation recommended treatment options, azithromycin, and ciprofloxacin. Recently, an MSM-associated outbreak of extended-spectrum β-lactamase producing, extensively drug resistant *S. sonnei* was reported in the United Kingdom. Here, we aimed to identify the genetic basis, evolutionary history, and international dissemination of the outbreak strain. Our genomic epidemiological analyses of 3,304 isolates from the United Kingdom, Australia, Belgium, France, and the United States of America revealed an internationally connected outbreak with a most recent common ancestor in 2018 carrying a low-fitness cost resistance plasmid, previously observed in travel associated sublineages of *S. flexneri*. Our results highlight the persistent threat of horizontally transmitted antimicrobial resistance and the value of continuing to work towards early and open international sharing of genomic surveillance data.

Bacteria belonging to the genus *Shigella* cause shigellosis, a gastrointestinal infection characterised by a combination of profuse, often bloody, diarrhoea with abdominal pain, fever, nausea and/or vomiting<fnref n="1" />. In middle- and higher-income nations, although symptoms

can be severe, the duration of illness is usually short-lived². However, in lower to middle income countries (LMIC), shigellosis is associated with higher levels of morbidity and mortality³. *Shigella* infection has no substantial non-human reservoir and is transmitted from person to

<footnote n="1">¹NIHR HPRU in Gastrointestinal Infections at University of Liverpool, Liverpool, UK. ²Department of Clinical Infection, Microbiology, and Immunology; Institute for Infection, Veterinary and Ecological Sciences, Liverpool, UK. ³Gastro and Food Safety (One Health) Division, UK Health Security Agency, London, UK.</footnote>

<footnote n="4">⁴Milner Centre for Evolution, University of Bath, Bath, UK. ⁵Centre for Infectious Diseases and Microbiology, The Westmead Institute for Medical Research, Westmead, NSW, Australia. ⁶Western Sydney Local Health District, Westmead, NSW, Australia. ⁷Faculty of Medicine and Health, University of Sydney, Sydney, NSW, Australia. ⁸Sydney Infectious Diseases Institute, University of Sydney, Sydney, NSW, Australia. ⁹New South Wales Health Pathology, Dee Why, NSW, Australia. ¹⁰Douglass Hanly Moir Pathology, Macquarie Park, NSW, Australia. ¹¹Centre for Infectious Diseases and Microbiology – Public Health, Institute for Clinical Pathology and Microbiology Research, Westmead Hospital, Westmead, NSW, Australia. ¹²Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia. ¹³Institut Pasteur, Université Paris Cité, Unité des Bactéries pathogènes entériques, Centre National de Référence des Escherichia coli, Shigella et Salmonella, Paris F-75015, France. ¹⁴Division of Human Bacterial Diseases, Sciensano, Ixelles, Belgium. ✉e-mail: kbaker@liverpool.ac.uk</footnote>
person by contact with an infected individual¹, or by the consumption of food and/or water contaminated by human faeces⁴,⁵. The four species of *Shigella* (*S. sonnei*, *S. flexneri*, *S. boydii* and *S. dysenteriae*) are responsible for approximately 212,400 deaths per year globally, 30% of which occur in young children aged under five years old³. In the UK, approximately 1650 cases of shigellosis are notified to UK Health Security Agency (UKHSA) each year, and *S. sonnei* is the most frequently detected species⁶. Prior to the 1980s, outbreaks of *S. sonnei* in the UK and other higher-income nations were associated with schools and other institutional settings⁷. However, in recent decades more infections were imported by travellers returning from low and middle income countries (LMIC)⁸, sometimes followed by on-going transmission within close knit communities⁹. Since 2012, sexual transmission of *S. sonnei* between men who have sex with men (MSM) has contributed to a high proportion of diagnoses in the UK¹⁰.

*S. sonnei* already exhibited multidrug resistance (MDR), to sulphonamides, ampicillin, streptomycin, and tetracycline, throughout the 1960s, conferred by different combinations of antimicrobial resistance (AMR) genes¹¹,¹². As a single serotype species, the splitting of *S. sonnei* into four genomic lineages¹³, and further clades and subclades¹⁴, has facilitated the study of its relationship with AMR and surveillance of global transmission. For example, a study of the introduction of *S. sonnei* into Vietnam highlighted the ability of *Shigella* to refine its accessory genome repertoire whereby introduced subclades were observed to acquire AMR genes, and genes encoding colicins that facilitate killing of other bacteria¹⁵. These characteristics contributed to its ability to become endemic in the human population and evolve over time towards resistance to an increasing number of antimicrobials. Furthermore, a recently defined subclade of the globally dominant genomic subtype of *S. sonnei*, Lineage 3¹³, disseminating from Southeast Asia, has acquired ciprofloxacin resistance, conferred by three mutations in the quinolone resistance determining region (QRDR) of the *gyrA* and *parC* genes¹⁶. More recently, ciprofloxacin resistant *S. sonnei* has acquired resistance to azithromycin, one of the second line agents, particularly in communities where antibiotic use is high¹⁷. The first global transmission of azithromycin resistant shigellosis was described among MSM, driven by the acquisition of a low-fitness cost AMR plasmid, and frequent use of antimicrobials for other sexually transmitted illnesses (STIs) leading to bystander resistance¹⁸,¹⁹. The acquisition of resistance to ciprofloxacin, recommended by the World Health Organisation (WHO) for empirical treatment of shigellosis²⁰, and increasing MDR in *S. sonnei* has earned it a place as a WHO priority organism for which new antimicrobials are urgently needed²¹.

The extensive global dissemination of MDR *S. sonnei* subsequently led to periodic outbreaks of XDR *S. sonnei*, defined here as being resistant to the only first-line antimicrobial (ciprofloxacin) and to all but one of the second-line antimicrobials (pivmecillinam, ceftriaxone and azithromycin) recommended by the WHO for the treatment of shigellosis²⁰. These outbreak isolates have resistance to ciprofloxacin, azithromycin and extended-spectrum β-lactam antibiotics (e.g., ceftriaxone)²². The latter resistance is typically associated with acquisition of an extended-spectrum β-lactamase (ESBL) gene, such as those in the *bla*<sub>CTX-M</sub> family. Sporadic, *bla*<sub>CTX-M</sub>-containing XDR *S. sonnei* isolates and local outbreaks have been previously described in various countries¹⁴,²³, but there has been little evidence of widespread dissemination of these subtypes, unlike their ciprofloxacin and ciprofloxacin-azithromycin resistant predecessors. Since 2015, transmission of MSM-associated shigellosis in the UK has been dominated by *S. sonnei* resistant to both azithromycin and ciprofloxacin⁶, until the COVID-19 pandemic led to a sharp reduction in case numbers²⁴. However, during the resurgence of shigellosis in the last quarter of 2021, an outbreak of *S. sonnei* occurred that likely involved multiple countries²²,²⁵,²⁶. To further explore the strain responsible for the newly sustained domestic transmission of XDR *S. sonnei* in the UK and explore its potential links to strains internationally (only possible with

genomics in this single serotype pathogen), we use genomic epidemiology to study the evolution of this strain from its MDR predecessor, the genetic basis for the XDR genotype, and the extent of dissemination across Europe, Australia, and the United States of America (USA).

## Results

### Evolutionary history of the UK XDR outbreak

To explore the evolutionary history and global context of *S. sonnei* belonging to CipR.MSM5 and specifically those isolates associated with the XDR outbreak in the United Kingdom, we constructed a cgMLST tree of a subset of 2820 clinical isolates of *S. sonnei* processed by UKHSA between 2016 and 2021, alongside 120 subclade representatives of the global genotyping framework (Fig. 1)¹⁴. The majority (*n* = 2730/2820, 97%) of UK isolates belonged to the globally dominant *S. sonnei* Lineage 3. Among these isolates, overlaying the year of isolation revealed a pattern of clonal replacement whereby isolates collected in later years tended to cluster together in more distant parts of the phylogeny and the sustained domestic transmission of a ciprofloxacin resistant sublineage (Fig. 1, Supplementary Fig. 1). This is consistent with previous reports of CipR.MSM5 (otherwise known as genotype 3.6.1.1.2) replacing earlier circulating sublineages of *S. sonnei* in the UK⁶,¹⁷. Among the UK isolates, predicted ESBL activity was solely associated with the presence of a *bla*<sub>CTX-M</sub> genes (*n* = 567, Supplementary Data 1). Historically, these genes had been seen only sporadically in conjunction with *mphA*, conferring azithromycin resistance in the ciprofloxacin resistant lineage (Fig. 1). More recently however, a sub-clade (a 10 SNP single linkage cluster, called t10.377) acquired *bla*<sub>CTX-M-27</sub>, to give rise to the XDR outbreak cluster (approximated by BAPS5 and indicated by the presence of *bla*<sub>CTX-M-27</sub> at the extremity of the tree) such that 90% (*n* = 35/39) of *S. sonnei* isolates from 2021 belonged to the XDR-outbreak (Fig. 1, see Supplementary Fig. 1 for full genotyping).

The patient demography of isolates in the outbreak linkage cluster supported epidemiological investigations that confirmed this was an MSM-associated outbreak²⁶. Specifically, isolates from male patients aged 16 to 60 years who had not reported recent international travel were overrepresented among the outbreak isolates, with these patients contributing 84% (405/482) of the UK isolates in t10.377 compared with 20% (475/2412) of UK isolates in the remainder of the phylogeny (*p*-value < 0.0001). We also found topological clustering of travel-associated isolates, i.e., isolates from patients recently returned from Asia co-located in the phylogeny, as did those from Africa (Fig. 1), reflecting the different bacterial populations endemic to these regions; a previously-described phenomenon¹⁷,¹⁸,²⁷. The MSM-associated outbreak lineage (t10.377) was almost exclusively (99%, 478/482) comprised of isolates from patients who did not report travel to Asia, Africa, or Latin America, compared to 69% (1661/2412) in the rest of the UK phylogeny (*p*-value < 0.00001).

### International linkage and accessory genome evolution of the outbreak

We explored potential international links and obtained a more detailed evolutionary history of the UK outbreak cluster (t10.377) by integrating and analysing datasets from international partners (Fig. 2). The subset of UK isolates belonging to t10.377 (*n* = 468) were analysed alongside 475 CipR.MSM5 isolates from Australia (*n* = 231), France (*n* = 101), Belgium (*n* = 112), and the USA (*n* = 31) (as these two genotypes t10.377 and CipR.MSM5/3.6.1.1.2 were 96% concordant (see Methods). These 943 isolates were constructed into a SNP-based phylogeny with improved resolution and the population was further divided into BAPS clusters (Fig. 2). This analysis supports the suggestion of clonal replacement of the CipR.MSM5 subclade over time, and geographic admixing within the BAPS clusters suggests that the evolution of MSM-associated *S. sonnei* was structured temporally rather than geographically (Fig. 2). Temporal analysis of the outbreak cluster showed
person by contact with an infected individual¹, or by the consumption of food and/or water contaminated by human faeces⁴,⁵. The four species of *Shigella* (*S. sonnei*, *S. flexneri*, *S. boydii* and *S. dysenteriae*) are responsible for approximately 212,400 deaths per year globally, 30% of which occur in young children aged under five years old³. In the UK, approximately 1650 cases of shigellosis are notified to UK Health Security Agency (UKHSA) each year, and *S. sonnei* is the most frequently detected species⁶. Prior to the 1980s, outbreaks of *S. sonnei* in the UK and other higher-income nations were associated with schools and other institutional settings⁷. However, in recent decades more infections were imported by travellers returning from low and middle income countries (LMIC)⁸, sometimes followed by on-going transmission within close knit communities⁹. Since 2012, sexual transmission of *S. sonnei* between men who have sex with men (MSM) has contributed to a high proportion of diagnoses in the UK¹⁰.

*S. sonnei* already exhibited multidrug resistance (MDR), to sulphonamides, ampicillin, streptomycin, and tetracycline, throughout the 1960s, conferred by different combinations of antimicrobial resistance (AMR) genes¹¹,¹². As a single serotype species, the splitting of *S. sonnei* into four genomic lineages¹³, and further clades and subclades¹⁴, has facilitated the study of its relationship with AMR and surveillance of global transmission. For example, a study of the introduction of *S. sonnei* into Vietnam highlighted the ability of *Shigella* to refine its accessory genome repertoire whereby introduced subclades were observed to acquire AMR genes, and genes encoding colicins that facilitate killing of other bacteria¹⁵. These characteristics contributed to its ability to become endemic in the human population and evolve over time towards resistance to an increasing number of antimicrobials. Furthermore, a recently defined subclade of the globally dominant genomic subtype of *S. sonnei*, Lineage 3¹³, disseminating from Southeast Asia, has acquired ciprofloxacin resistance, conferred by three mutations in the quinolone resistance determining region (QRDR) of the *gyrA* and *parC* genes¹⁶. More recently, ciprofloxacin resistant *S. sonnei* has acquired resistance to azithromycin, one of the second line agents, particularly in communities where antibiotic use is high¹⁷. The first global transmission of azithromycin resistant shigellosis was described among MSM, driven by the acquisition of a low-fitness cost AMR plasmid, and frequent use of antimicrobials for other sexually transmitted illnesses (STIs) leading to bystander resistance¹⁸,¹⁹. The acquisition of resistance to ciprofloxacin, recommended by the World Health Organisation (WHO) for empirical treatment of shigellosis²⁰, and increasing MDR in *S. sonnei* has earned it a place as a WHO priority organism for which new antimicrobials are urgently needed²¹.

The extensive global dissemination of MDR *S. sonnei* subsequently led to periodic outbreaks of XDR *S. sonnei*, defined here as being resistant to the only first-line antimicrobial (ciprofloxacin) and to all but one of the second-line antimicrobials (pivmecillinam, ceftriaxone and azithromycin) recommended by the WHO for the treatment of shigellosis²⁰. These outbreak isolates have resistance to ciprofloxacin, azithromycin and extended-spectrum β-lactam antibiotics (e.g., ceftriaxone)²². The latter resistance is typically associated with acquisition of an extended-spectrum β-lactamase (ESBL) gene, such as those in the *bla*<sub>CTX-M</sub> family. Sporadic, *bla*<sub>CTX-M</sub>-containing XDR *S. sonnei* isolates and local outbreaks have been previously described in various countries¹⁴,²³, but there has been little evidence of widespread dissemination of these subtypes, unlike their ciprofloxacin and ciprofloxacin-azithromycin resistant predecessors. Since 2015, transmission of MSM-associated shigellosis in the UK has been dominated by *S. sonnei* resistant to both azithromycin and ciprofloxacin⁶, until the COVID-19 pandemic led to a sharp reduction in case numbers²⁴. However, during the resurgence of shigellosis in the last quarter of 2021, an outbreak of *S. sonnei* occurred that likely involved multiple countries²²,²⁵,²⁶. To further explore the strain responsible for the newly sustained domestic transmission of XDR *S. sonnei* in the UK and explore its potential links to strains internationally (only possible with

genomics in this single serotype pathogen), we use genomic epidemiology to study the evolution of this strain from its MDR predecessor, the genetic basis for the XDR genotype, and the extent of dissemination across Europe, Australia, and the United States of America (USA).

## Results

### Evolutionary history of the UK XDR outbreak

To explore the evolutionary history and global context of *S. sonnei* belonging to CipR.MSM5 and specifically those isolates associated with the XDR outbreak in the United Kingdom, we constructed a cgMLST tree of a subset of 2820 clinical isolates of *S. sonnei* processed by UKHSA between 2016 and 2021, alongside 120 subclade representatives of the global genotyping framework (Fig. 1)¹⁴. The majority (*n* = 2730/2820, 97%) of UK isolates belonged to the globally dominant *S. sonnei* Lineage 3. Among these isolates, overlaying the year of isolation revealed a pattern of clonal replacement whereby isolates collected in later years tended to cluster together in more distant parts of the phylogeny and the sustained domestic transmission of a ciprofloxacin resistant sublineage (Fig. 1, Supplementary Fig. 1). This is consistent with previous reports of CipR.MSM5 (otherwise known as genotype 3.6.1.1.2) replacing earlier circulating sublineages of *S. sonnei* in the UK⁶,¹⁷. Among the UK isolates, predicted ESBL activity was solely associated with the presence of a *bla*<sub>CTX-M</sub> genes (*n* = 567, Supplementary Data 1). Historically, these genes had been seen only sporadically in conjunction with *mphA*, conferring azithromycin resistance in the ciprofloxacin resistant lineage (Fig. 1). More recently however, a sub-clade (a 10 SNP single linkage cluster, called t10.377) acquired *bla*<sub>CTX-M-27</sub>, to give rise to the XDR outbreak cluster (approximated by BAPS5 and indicated by the presence of *bla*<sub>CTX-M-27</sub> at the extremity of the tree) such that 90% (*n* = 35/39) of *S. sonnei* isolates from 2021 belonged to the XDR-outbreak (Fig. 1, see Supplementary Fig. 1 for full genotyping).

The patient demography of isolates in the outbreak linkage cluster supported epidemiological investigations that confirmed this was an MSM-associated outbreak²⁶. Specifically, isolates from male patients aged 16 to 60 years who had not reported recent international travel were overrepresented among the outbreak isolates, with these patients contributing 84% (405/482) of the UK isolates in t10.377 compared with 20% (475/2412) of UK isolates in the remainder of the phylogeny (*p*-value < 0.0001). We also found topological clustering of travel-associated isolates, i.e., isolates from patients recently returned from Asia co-located in the phylogeny, as did those from Africa (Fig. 1), reflecting the different bacterial populations endemic to these regions; a previously-described phenomenon¹⁷,¹⁸,²⁷. The MSM-associated outbreak lineage (t10.377) was almost exclusively (99%, 478/482) comprised of isolates from patients who did not report travel to Asia, Africa, or Latin America, compared to 69% (1661/2412) in the rest of the UK phylogeny (*p*-value < 0.00001).

### International linkage and accessory genome evolution of the outbreak

We explored potential international links and obtained a more detailed evolutionary history of the UK outbreak cluster (t10.377) by integrating and analysing datasets from international partners (Fig. 2). The subset of UK isolates belonging to t10.377 (*n* = 468) were analysed alongside 475 CipR.MSM5 isolates from Australia (*n* = 231), France (*n* = 101), Belgium (*n* = 112), and the USA (*n* = 31) (as these two genotypes t10.377 and CipR.MSM5/3.6.1.1.2 were 96% concordant (see Methods). These 943 isolates were constructed into a SNP-based phylogeny with improved resolution and the population was further divided into BAPS clusters (Fig. 2). This analysis supports the suggestion of clonal replacement of the CipR.MSM5 subclade over time, and geographic admixing within the BAPS clusters suggests that the evolution of MSM-associated *S. sonnei* was structured temporally rather than geographically (Fig. 2). Temporal analysis of the outbreak cluster showed
<figure><img src="image_4.png" /></figure>

**Fig. 1 |** *The emergence of an XDR *S. sonnei* outbreak in the United Kingdom. A* cgMLST dendrogram (midpoint rooted) of clinical isolates from the UK (n = 2,820) and genomic subtype references (n = 120) with the scale bar indicating distance in cgMLST alleles. Metadata tracks show patient and genomic features for isolates coloured according to the inlaid keys. Specifically, from inner to outer, the patient data comprises: year of isolation, patient sex, age group, and travel history (for UK isolates only), with missing/unavailable/not determined data shown as white.

Genomic features then show: isolate lineage, BAPS clusters (for isolates belonging to the t10.377 cluster only), and the presence of mutations in the Quinolone Resistance Determining Region (QRDR*3 denotes all three canonical mutations; *gyrA_D87G*, *gyrA_S83L* and *parC_S80I*) and *mph(A)*, and *bla*<sub>CTX-M</sub> genes, where white indicates absence of the gene.

that the CipR.MSM5 genotype likely emerged in 2014, and that the BAPS 5 cluster emerged in 2018 (Supplementary Fig. 4).

The clustering of CipR.MSM5 isolates carrying *bla*<sub>CTX-M-27</sub> from all countries in BAPS5 indicates that this recently emerged monophyletic lineage of *S. sonnei* was circulating intercontinentally across regions historically considered low-risk for shigellosis. To investigate the evolution of the XDR phenotype in the BAPS5 outbreak cluster, AMR gene profiles were determined for isolates in BAPS1-5 from the UK, Belgium, and NSW, as sets of isolates from other locations included only those with a *bla*<sub>CTX-M</sub> gene, (Table 1). This was complemented by phenotypic confirmation for a subset of 14 UK isolates from the CipR.MSM5 subclade (Supplementary Table 1). All isolates were MDR and acquisition of *bla*<sub>CTX-M-27</sub> in the BAPS cluster 5 correlated with phenotypic ceftriaxone resistance. BAPS5 isolates lacked *bla*<sub>TEM-1</sub> compared with other BAPS, but as expected this did not affect their resistance profile against the antibiotics tested (Supplementary Table 1).

The resistance profiles against the recommended therapeutic antimicrobials (azithromycin, ciprofloxacin, and third-generation cephalosporins) were broadly consistent across BAPS clusters 1–4. Although sharing similar levels of azithromycin resistance with these clusters, the BAPS5 outbreak cluster distinguished itself early in its evolution by acquiring *qnrB19*, and losing *bla*<sub>TEM-1</sub> (Fig. 2, Table 1). Notably, *qnrB19* (found on a 2579 bp plasmid) in BAPS5 seems redundant as all isolates in the broader CipR.MSM5 subclade have QRDR mutations that conferred full resistance to ciprofloxacin (MIC ≥ 4 µg/mL), with no increase in MIC associated with the presence of *qnrB19* (Supplementary Table 1). Following the acquisition of this

*qnrB19* plasmid, the BAPS5 cluster subsequently gained *bla*<sub>CTX-M-27</sub> which confers resistance to ceftriaxone (MICs ≥12 µg/mL, Supplementary Table 1).

## The *bla*<sub>CTX-M-27</sub> containing plasmid, p893816

Analysis of hybrid assemblies from eight outbreak isolates from the UK revealed that all contained a near identical ~83.4 kbp FII plasmid that carries *bla*<sub>CTX-M-27</sub>. Plasmids from single isolates from France and NSW also match the plasmid from the UK and a BLASTn search revealed that all of these plasmids are also near identical to p893816, reported in MSM-associated CipR.MSM5 *S. sonnei* in London in 2020²⁸. This plasmid carries several other AMR genes (Table 1, Supplementary Figs. 2, 3) interspersed with transposon fragments and insertion sequences that may facilitate their movement (Fig. 3a). The plasmid backbone encodes conjugation and stabilisation machinery and carries *psiB* (plasmid SOS inhibition), recently shown to inhibit the disruptive host SOS response to subminimal inhibitory concentrations of ciprofloxacin¹⁹.

Mapping of Illumina reads from BAPS 5 isolates from all countries (n = 233) revealed that coverage was 96.1–98.7% for most with *bla*<sub>CTX-M-27</sub>-containing isolates (n = 161/184, 87.5%), including those with assembled plasmids virtually identical to p893816 (Table 2, Supplementary Fig. 5) coverage. Coverage for the remaining BAPS 5 isolates with *bla*<sub>CTX-M-27</sub> and for those without this gene (n = 49) was distinctly lower (81–93%; Supplementary Fig. 5). This suggests that a p893916-like plasmid carries *bla*<sub>CTX-M-27</sub> in the majority of BAPS 5 isolates, but that less related plasmid(s) are present in a minority.
<figure><img src="image_4.png" /></figure>

**Fig. 1 |** *The emergence of an XDR *S. sonnei* outbreak in the United Kingdom.* A cgMLST dendrogram (midpoint rooted) of clinical isolates from the UK (n = 2,820) and genomic subtype references (n = 120) with the scale bar indicating distance in cgMLST alleles. Metadata tracks show patient and genomic features for isolates coloured according to the inlaid keys. Specifically, from inner to outer, the patient data comprises: year of isolation, patient sex, age group, and travel history (for UK isolates only), with missing/unavailable/not determined data shown as white.

Genomic features then show: isolate lineage, BAPS clusters (for isolates belonging to the t10.377 cluster only), and the presence of mutations in the Quinolone Resistance Determining Region (QRDR*3 denotes all three canonical mutations; *gyrA_D87G*, *gyrA_S83L* and *parC_S80I*) and *mph(A)*, and *bla*<sub>CTX-M</sub> genes, where white indicates absence of the gene.

that the CipR.MSM5 genotype likely emerged in 2014, and that the BAPS 5 cluster emerged in 2018 (Supplementary Fig. 4).

The clustering of CipR.MSM5 isolates carrying *bla*<sub>CTX-M-27</sub> from all countries in BAPS5 indicates that this recently emerged monophyletic lineage of *S. sonnei* was circulating intercontinentally across regions historically considered low-risk for shigellosis. To investigate the evolution of the XDR phenotype in the BAPS5 outbreak cluster, AMR gene profiles were determined for isolates in BAPS1-5 from the UK, Belgium, and NSW, as sets of isolates from other locations included only those with a *bla*<sub>CTX-M</sub> gene, (Table 1). This was complemented by phenotypic confirmation for a subset of 14 UK isolates from the CipR.MSM5 subclade (Supplementary Table 1). All isolates were MDR and acquisition of *bla*<sub>CTX-M-27</sub> in the BAPS cluster 5 correlated with phenotypic ceftriaxone resistance. BAPS5 isolates lacked *bla*<sub>TEM-1</sub> compared with other BAPS, but as expected this did not affect their resistance profile against the antibiotics tested (Supplementary Table 1).

The resistance profiles against the recommended therapeutic antimicrobials (azithromycin, ciprofloxacin, and third-generation cephalosporins) were broadly consistent across BAPS clusters 1-4. Although sharing similar levels of azithromycin resistance with these clusters, the BAPS5 outbreak cluster distinguished itself early in its evolution by acquiring *qnrB19*, and losing *bla*<sub>TEM-1</sub> (Fig. 2, Table 1). Notably, *qnrB19* (found on a 2579 bp plasmid) in BAPS5 seems redundant as all isolates in the broader CipR.MSM5 subclade have QRDR mutations that conferred full resistance to ciprofloxacin (MIC ≥ 4 µg/mL), with no increase in MIC associated with the presence of *qnrB19* (Supplementary Table 1). Following the acquisition of this

*qnrB19* plasmid, the BAPS5 cluster subsequently gained *bla*<sub>CTX-M-27</sub> which confers resistance to ceftriaxone (MICs ≥12 µg/mL, Supplementary Table 1).

## The *bla*<sub>CTX-M-27</sub> containing plasmid, p893816

Analysis of hybrid assemblies from eight outbreak isolates from the UK revealed that all contained a near identical ~83.4 kbp FII plasmid that carries *bla*<sub>CTX-M-27</sub>. Plasmids from single isolates from France and NSW also match the plasmid from the UK and a BLASTn search revealed that all of these plasmids are also near identical to p893816, reported in MSM-associated CipR.MSM5 *S. sonnei* in London in 2020²⁸. This plasmid carries several other AMR genes (Table 1, Supplementary Figs. 2, 3) interspersed with transposon fragments and insertion sequences that may facilitate their movement (Fig. 3a). The plasmid backbone encodes conjugation and stabilisation machinery and carries *psiB* (plasmid SOS inhibition), recently shown to inhibit the disruptive host SOS response to subminimal inhibitory concentrations of ciprofloxacin¹⁹.

Mapping of Illumina reads from BAPS 5 isolates from all countries (n = 233) revealed that coverage was 96.1–98.7% for most with *bla*<sub>CTX-M-27</sub>-containing isolates (n = 161/184, 87.5%), including those with assembled plasmids virtually identical to p893816 (Table 2, Supplementary Fig. 5) coverage. Coverage for the remaining BAPS 5 isolates with *bla*<sub>CTX-M-27</sub> and for those without this gene (n = 49) was distinctly lower (81–93%; Supplementary Fig. 5). This suggests that a p893916-like plasmid carries *bla*<sub>CTX-M-27</sub> in the majority of BAPS 5 isolates, but that less related plasmid(s) are present in a minority.
<figure><img src="image_5.png" /><figcaption>**Fig. 2 |** The evolution and international spread of MSM-associated XDR *S. sonnei*. A midpoint rooted maximum likelihood phylogenetic tree shows the distribution of UK isolates (belonging to both CipR.MSM5 and the t10.377 outbreak cluster, *n* = 468) and relevant related international isolates belonging to CipR.MSM5 (Supplementary Fig. 1, *n* = 475). Metadata tracks show year and country</figcaption></figure>

(area) of isolation, BAPS subtype, and the presence of selected AMR genes according to the inlaid keys. The scalebar is provided by IQTree, and represents expected number of substitutions per site across a 1717 bp alignment. Bold branches represent a bootstrap value of ≥ 70 out of 100.

<table><caption>**Table 1 |** Proportion of isolates from the UK, Belgium and NSW carrying different AMR genes by BAPS cluster</caption><thead><tr><th rowspan="2">Antimicrobial class</th><th rowspan="2">Gene<sup>a</sup></th><th colspan="5">Proportion of BAPS cluster containing gene (%)</th></tr><tr><th>1 (n = 240)</th><th>2 (n = 53)</th><th>3 (n = 249)</th><th>4 (n = 148)</th><th>5 (n = 118)</th></tr></thead><tbody><tr><td rowspan="3">Aminoglycosides</td><td>aadA5</td><td>97.08</td><td>100</td><td>97.99</td><td>76.35</td><td>100</td></tr><tr><td>aph(3')-lb/strA</td><td>87.92</td><td>98.11</td><td>85.94</td><td>79.05</td><td>83.05</td></tr><tr><td>aph(6)-ld/strB</td><td>87.92</td><td>98.11</td><td>85.94</td><td>79.05</td><td>83.05</td></tr><tr><td>β-lactams</td><td>bla<sub>TEM-1</sub></td><td>97.08</td><td>94.34</td><td>97.19</td><td>76.35</td><td>36.44</td></tr><tr><td>Cephalosporins</td><td>bla<sub>CTX-M-27</sub></td><td>0</td><td>1.89</td><td>0</td><td>0</td><td>58.47</td></tr><tr><td>Fluroquinolones<sup>b</sup></td><td>qnrB19</td><td>0</td><td>0</td><td>0</td><td>0</td><td>94.92</td></tr><tr><td rowspan="2">Macrolides</td><td>erm(B)</td><td>97.50</td><td>100</td><td>99.60</td><td>93.92</td><td>100</td></tr><tr><td>mph(A)</td><td>97.08</td><td>100</td><td>97.99</td><td>93.92</td><td>100</td></tr><tr><td rowspan="2">Sulphonamides</td><td>sul1</td><td>89.17</td><td>98.11</td><td>85.94</td><td>79.05</td><td>84.75</td></tr><tr><td>sul2</td><td>87.92</td><td>96.23</td><td>85.94</td><td>79.73</td><td>83.05</td></tr><tr><td>Tetracyclines</td><td>tet(A)</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td></tr><tr><td>Trimethoprim</td><td>dfrA17</td><td>97.08</td><td>100</td><td>97.99</td><td>81.08</td><td>100</td></tr></tbody></table>

These plasmids also match Plasmid 3 from an Australian *S. flexneri* 3a isolated even earlier in 2018 (Table 2). Finding the same plasmid in another species suggests it is mobilising among *Shigella*, consistent with previous work showing that p893816 had close relatives in another circulating genotype of MSM-associated *S. sonnei*, MSM2/

VN.KH1<sup>29</sup>, and echoing the global expansion of the azithromycin resistance pKSR100-like plasmids among different *Shigella* species<sup>17,19</sup>.

To determine the extent and spread of close relatives of p893816 among bacteria, we compared the plasmid sequence with publicly available data. Firstly, comparisons (using BLAST) of the

<footnote><sup>a</sup>All isolates also carry *sat2* and *dfrA1* and the intrinsic chromosomal *ampC* (annotated as *blaEC*), *acrF* and *emrD* genes.</footnote>

<footnote><sup>b</sup>All isolates have *gyrA_D87G*, *gyrA_S83L* and *parC_S80I* QRDR mutations associated with quinolone resistance.</footnote>
<figure><img src="image_5.png" /><figcaption>**Fig. 2 |** The evolution and international spread of MSM-associated XDR *S. sonnei*. A midpoint rooted maximum likelihood phylogenetic tree shows the distribution of UK isolates (belonging to both CipR.MSM5 and the t10.377 outbreak cluster, *n* = 468) and relevant related international isolates belonging to CipR.MSM5 (Supplementary Fig. 1, *n* = 475). Metadata tracks show year and country</figcaption></figure>

(area) of isolation, BAPS subtype, and the presence of selected AMR genes according to the inlaid keys. The scalebar is provided by IQTree, and represents expected number of substitutions per site across a 1717 bp alignment. Bold branches represent a bootstrap value of ≥ 70 out of 100.

<table><caption>**Table 1 |** Proportion of isolates from the UK, Belgium and NSW carrying different AMR genes by BAPS cluster</caption><thead><tr><th rowspan="2">Antimicrobial class</th><th rowspan="2">Gene<sup>a</sup></th><th colspan="5">Proportion of BAPS cluster containing gene (%)</th></tr><tr><th>1 (n = 240)</th><th>2 (n = 53)</th><th>3 (n = 249)</th><th>4 (n = 148)</th><th>5 (n = 118)</th></tr></thead><tbody><tr><td rowspan="3">Aminoglycosides</td><td>aadA5</td><td>97.08</td><td>100</td><td>97.99</td><td>76.35</td><td>100</td></tr><tr><td>aph(3')-lb/strA</td><td>87.92</td><td>98.11</td><td>85.94</td><td>79.05</td><td>83.05</td></tr><tr><td>aph(6)-ld/strB</td><td>87.92</td><td>98.11</td><td>85.94</td><td>79.05</td><td>83.05</td></tr><tr><td>β-lactams</td><td>bla<sub>TEM-1</sub></td><td>97.08</td><td>94.34</td><td>97.19</td><td>76.35</td><td>36.44</td></tr><tr><td>Cephalosporins</td><td>bla<sub>CTX-M-27</sub></td><td>0</td><td>1.89</td><td>0</td><td>0</td><td>58.47</td></tr><tr><td>Fluroquinolones<sup>b</sup></td><td>qnrB19</td><td>0</td><td>0</td><td>0</td><td>0</td><td>94.92</td></tr><tr><td rowspan="2">Macrolides</td><td>erm(B)</td><td>97.50</td><td>100</td><td>99.60</td><td>93.92</td><td>100</td></tr><tr><td>mph(A)</td><td>97.08</td><td>100</td><td>97.99</td><td>93.92</td><td>100</td></tr><tr><td rowspan="2">Sulphonamides</td><td>sul1</td><td>89.17</td><td>98.11</td><td>85.94</td><td>79.05</td><td>84.75</td></tr><tr><td>sul2</td><td>87.92</td><td>96.23</td><td>85.94</td><td>79.73</td><td>83.05</td></tr><tr><td>Tetracyclines</td><td>tet(A)</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td></tr><tr><td>Trimethoprim</td><td>dfrA17</td><td>97.08</td><td>100</td><td>97.99</td><td>81.08</td><td>100</td></tr></tbody></table>

These plasmids also match Plasmid 3 from an Australian *S. flexneri* 3a isolated even earlier in 2018 (Table 2). Finding the same plasmid in another species suggests it is mobilising among *Shigella*, consistent with previous work showing that p893816 had close relatives in another circulating genotype of MSM-associated *S. sonnei*, MSM2/

VN.KH1<sup>29</sup>, and echoing the global expansion of the azithromycin resistance pKSR100-like plasmids among different *Shigella* species<sup>17,19</sup>.

To determine the extent and spread of close relatives of p893816 among bacteria, we compared the plasmid sequence with publicly available data. Firstly, comparisons (using BLAST) of the

<footnote><sup>a</sup>All isolates also carry sat2 and dfrA1 and the intrinsic chromosomal *ampC* (annotated as *blaEC*), *acrF* and *emrD* genes.</footnote>

<footnote><sup>b</sup>All isolates have gyrA_D87G, gyrA_S83L and parC_S80I QRDR mutations associated with quinolone resistance.</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** | Organisation and distribution of p893816-like plasmids. **a** Organisation of p893816-like plasmids. The lengths and directions of transcription of various genes are shown as arrows. The upper track shows the resistance region. AMR genes indicated by solid red arrows, with resistance-related genes outlined in red. The boxes around *dfrA17* and *aadA5* represent the gene cassettes that carry these genes. Block arrows in white indicate various insertion sequences (IS), named/ numbered as labelled, and the direction of their transposase gene. Fragments of different transposons (Tn) are shown as interrupted blocks in different colours (blue, purple, orange) and labelled by name, and a class 1 integron is also shown in grey and labelled. The lower track shows the plasmid backbone, with the AMR</figcaption></figure>

region insertion site indicated by dashed lines. Replication genes are in pink, conjugation genes in purple (tra) and grey (trb), toxin-antitoxin systems in blue, other genes encoding known plasmid functions in cyan and genes encoding hypothetical proteins in white. **b** The distribution of p893816-like plasmids. Number of draft genome assemblies in the 661 K COBS data structure showing k-mer similarity (≥0.80) to p893916, coloured by species and sequence type according to the inlaid keys. ST that are represented 10 or more times are named, while rare sequence types are collapsed into the category of 'other' within the respective species.

p893816 sequence with the nucleotide collection at the National Centre for Biotechnology Information (NCBI) revealed multiple similar sequences (≥95% identity, Supplementary Fig. 3). Secondly, a broader search for bacterial isolates containing elements of p893816 was performed using an indexed data structure created from 661,405 uniformly assembled bacterial genomes that were present in the European Nucleotide Archive in November of 2018 (see methods). The index was queried with the constituent k-mers (short sequences of length k) of p893816, identifying datasets that contained ≥80% of the query k-mers; or ≥80% k-mer similarity. It is important to note that

k-mer similarity drops more rapidly than nucleotide identity, as for each SNP difference between the query and the nearest sequence in the index, a window of overlapping k-mers is absent. This search confirmed that some *Shigella* and *E. coli* genomes had sequences of high k-mer similarity and that some more distantly related *Enterobacteriaceae* harboured more distant p893816 relatives (Fig. 3b). A limitation of this analysis is that the index was built from a snapshot in 2018, when this plasmid appears to have emerged, and will miss relationships to more recently deposited sequences. This highlights the value of dynamic ongoing surveillance and how rapid analysis of the
<figure><img src="image_6.png" /><figcaption>**Fig. 3** | Organisation and distribution of p893816-like plasmids. **a** Organisation of p893816-like plasmids. The lengths and directions of transcription of various genes are shown as arrows. The upper track shows the resistance region. AMR genes indicated by solid red arrows, with resistance-related genes outlined in red. The boxes around *dfrA17* and *aadA5* represent the gene cassettes that carry these genes. Block arrows in white indicate various insertion sequences (IS), named/ numbered as labelled, and the direction of their transposase gene. Fragments of different transposons (Tn) are shown as interrupted blocks in different colours (blue, purple, orange) and labelled by name, and a class 1 integron is also shown in grey and labelled. The lower track shows the plasmid backbone, with the AMR</figcaption></figure>

region insertion site indicated by dashed lines. Replication genes are in pink, conjugation genes in purple (tra) and grey (trb), toxin-antitoxin systems in blue, other genes encoding known plasmid functions in cyan and genes encoding hypothetical proteins in white. **b** The distribution of p893816-like plasmids. Number of draft genome assemblies in the 661 K COBS data structure showing k-mer similarity (≥0.80) to p893916, coloured by species and sequence type according to the inlaid keys. ST that are represented 10 or more times are named, while rare sequence types are collapsed into the category of 'other' within the respective species.

p893816 sequence with the nucleotide collection at the National Centre for Biotechnology Information (NCBI) revealed multiple similar sequences (≥95% identity, Supplementary Fig. 3). Secondly, a broader search for bacterial isolates containing elements of p893816 was performed using an indexed data structure created from 661,405 uniformly assembled bacterial genomes that were present in the European Nucleotide Archive in November of 2018 (see methods). The index was queried with the constituent k-mers (short sequences of length k) of p893816, identifying datasets that contained ≥80% of the query k-mers; or ≥80% k-mer similarity. It is important to note that

k-mer similarity drops more rapidly than nucleotide identity, as for each SNP difference between the query and the nearest sequence in the index, a window of overlapping k-mers is absent. This search confirmed that some *Shigella* and *E. coli* genomes had sequences of high k-mer similarity and that some more distantly related *Enterobacteriaceae* harboured more distant p893816 relatives (Fig. 3b). A limitation of this analysis is that the index was built from a snapshot in 2018, when this plasmid appears to have emerged, and will miss relationships to more recently deposited sequences. This highlights the value of dynamic ongoing surveillance and how rapid analysis of the
<table><caption>Table 2 | Details of selected p893816-like plasmids compared in this study (and see Supplementary)</caption><thead><tr><th>Plasmid</th><th>p893816</th><th>p1538171_3</th><th>p202008564-6</th><th>p20-001-0088_F</th><th>Plasmid-3</th></tr></thead><tbody><tr><td>Accession</td><td>MW396858.1</td><td>CP104412</td><td>OP038290</td><td>CP115395</td><td>LR861790</td></tr><tr><td>Strain name</td><td>893816</td><td>1538171</td><td>202008564</td><td>20-001-0088</td><td>AUSMDU0022017</td></tr><tr><td>Length (bp)</td><td>83,397</td><td>83,416</td><td>83,397</td><td>83,389</td><td>83,425</td></tr><tr><td>Species</td><td>S. sonnei</td><td>S. sonnei</td><td>S. sonnei</td><td>S. sonnei</td><td>S. flexneri</td></tr><tr><td>Subtype</td><td>CipR.MSM5</td><td>ND</td><td>CipR.MSM5, BAPS 5</td><td>CipR.MSM5, BAPS 5</td><td>Serotype 3a</td></tr><tr><td>Year</td><td>2020</td><td>2021</td><td>2020</td><td>2019</td><td>2018</td></tr><tr><td>Country</td><td>UK</td><td>UK</td><td>France</td><td>Australia</td><td>Australia</td></tr><tr><td>Reference</td><td>Locke et al.<sup>28</sup></td><td>This study</td><td>This study</td><td>This study</td><td>Publicly available data</td></tr></tbody></table>

distribution of mobile genetic elements would further our understanding of transmissible AMR.

As our previous investigations suggested that pKSR100 (GenBank accession CP090162), an epidemiologically successful *Shigella* AMR plasmid, conferred a lower fitness cost on the bacterial host than another less successful co-circulating AMR plasmid, pAPR100 (GenBank accession CP090161)¹⁹, we explored the fitness cost of the p893816-like plasmid p1538171_3 in a naïve *E. coli* host. Analogous experiments here based on area under the curve analysis of growth curves showed that p1538171_3 had a similarly low fitness cost as pKSR100 (Fig. 4). Minimum inhibitory concentration experiments confirmed functional expression of the AMR genes present on the p1538171_3 plasmid (Supplementary Table 2). Unfortunately, this portends that this MDR-plasmid may readily spread throughout MSM-associated *Shigella* globally, perhaps in further MDR strains that have already acquired ciprofloxacin resistance.

Discussion

Antimicrobial use for treatment of STIs among the at-risk MSM community³⁰, the possibility of prolonged infection in this group, and the propensity of *Shigella* to acquire genes through horizontal gene transfer from the commensal gut microflora³¹,³², creates conditions where these pathogens can hone their AMR repertoire. MDR and XDR strains have evolved to become problematic within the MSM community, creating a significant burden with respect to potential treatment failure, and both persistent and reoccurring infections³³. In 2004, the first UK outbreak of sexually transmitted shigellosis was caused by *S. sonnei*, though relatively few of the isolates carried azithromycin resistance genes such as *mph*(A)³⁴. A spate of similar MSM-associated outbreaks of MDR *S. sonnei* occurred across other countries at this time, and perhaps represented an early example of international dissemination of *Shigella* throughout the global MSM community, which has been subsequently observed for MDR subtypes, including azithromycin resistant *Shigella*¹⁸ and now the BAPS5 lineage described in this study. The acquisition of resistance against recommended treatments for shigellosis and other STIs is a global public health concern, and potentially herald an increased risk of treatment failure³⁵, and eventual emergence of pan-drug resistant *S. sonnei*.

The drivers and trajectory of this emergent XDR lineage, which is additionally ceftriaxone resistant compared with previously described internationally disseminated MSM-associated *Shigella*, are likely to be similar to previously identified phenomena. Specifically, evolution of resistance to third-generation cephalosporins is likely being driven by bystander resistance resulting from shifting treatment regimens for other sexually transmissible illnesses. While historical use of azithromycin for treating many STIs led to the acquisition of *mph*(A) and *erm*(B) among MSM associated *Shigella*¹⁸, the current first-line treatment recommendation for gonorrhoea in the USA, UK, France, and Belgium is now third-generation cephalosporins (e.g. ceftriaxone)³⁶⁻³⁹. Treatment for ongoing high rates of gonorrhoea in the UK⁴⁰ thus likely

<figure><img src="image_7.png" /><figcaption>**Fig. 4 | Relative fitness of different *Shigella* MDR plasmids in *E. coli*.** Relative fitness of *E. coli* MG1655 carrying a previously characterised azithromycin resistance plasmid (pAPR100, green; pKSR100, purple), or p1538171_3 (navy) compared with a plasmid-free strain (grey). Each dataset is represented by a box and whisker plot indicating maximum and minimum values (upper and lower limits, respectively), and interquartile ranges (25th to 75th percentiles) of the original data, overlaid with black dots, indicating individual data points. The centre emboldened line on each box represents the median value. The results of pairwise statistical comparisons are shown above and with two-tailed *p*-values (two sample t-tests) indicated above black lines joining the pair comparators. The sample size for statistical purposes is *n* = 3 biologically independent samples (see methods).</figcaption></figure>
<table><caption>Table 2 | Details of selected p893816-like plasmids compared in this study (and see Supplementary)</caption><thead><tr><th>Plasmid</th><th>p893816</th><th>p1538171_3</th><th>p202008564-6</th><th>p20-001-0088_F</th><th>Plasmid-3</th></tr></thead><tbody><tr><td>Accession</td><td>MW396858.1</td><td>CP104412</td><td>OP038290</td><td>CP115395</td><td>LR861790</td></tr><tr><td>Strain name</td><td>893816</td><td>1538171</td><td>202008564</td><td>20-001-0088</td><td>AUSMDU0022017</td></tr><tr><td>Length (bp)</td><td>83,397</td><td>83,416</td><td>83,397</td><td>83,389</td><td>83,425</td></tr><tr><td>Species</td><td>S. sonnei</td><td>S. sonnei</td><td>S. sonnei</td><td>S. sonnei</td><td>S. flexneri</td></tr><tr><td>Subtype</td><td>CipR.MSM5</td><td>ND</td><td>CipR.MSM5, BAPS 5</td><td>CipR.MSM5, BAPS 5</td><td>Serotype 3a</td></tr><tr><td>Year</td><td>2020</td><td>2021</td><td>2020</td><td>2019</td><td>2018</td></tr><tr><td>Country</td><td>UK</td><td>UK</td><td>France</td><td>Australia</td><td>Australia</td></tr><tr><td>Reference</td><td>Locke et al.<sup>28</sup></td><td>This study</td><td>This study</td><td>This study</td><td>Publicly available data</td></tr></tbody></table>

distribution of mobile genetic elements would further our understanding of transmissible AMR.

As our previous investigations suggested that pKSR100 (GenBank accession CP090162), an epidemiologically successful *Shigella* AMR plasmid, conferred a lower fitness cost on the bacterial host than another less successful co-circulating AMR plasmid, pAPR100 (GenBank accession CP090161)¹⁹, we explored the fitness cost of the p893816-like plasmid p1538171_3 in a naïve *E. coli* host. Analogous experiments here based on area under the curve analysis of growth curves showed that p1538171_3 had a similarly low fitness cost as pKSR100 (Fig. 4). Minimum inhibitory concentration experiments confirmed functional expression of the AMR genes present on the p1538171_3 plasmid (Supplementary Table 2). Unfortunately, this tends that this MDR-plasmid may readily spread throughout MSM-associated *Shigella* globally, perhaps in further MDR strains that have already acquired ciprofloxacin resistance.

Discussion

Antimicrobial use for treatment of STIs among the at-risk MSM community³⁰, the possibility of prolonged infection in this group, and the propensity of *Shigella* to acquire genes through horizontal gene transfer from the commensal gut microflora³¹,³², creates conditions where these pathogens can hone their AMR repertoire. MDR and XDR strains have evolved to become problematic within the MSM community, creating a significant burden with respect to potential treatment failure, and both persistent and reoccurring infections³³. In 2004, the first UK outbreak of sexually transmitted shigellosis was caused by *S. sonnei*, though relatively few of the isolates carried azithromycin resistance genes such as *mph*(A)³⁴. A spate of similar MSM-associated outbreaks of MDR *S. sonnei* occurred across other countries at this time, and perhaps represented an early example of international dissemination of *Shigella* throughout the global MSM community, which has been subsequently observed for MDR subtypes, including azithromycin resistant *Shigella*¹⁸ and now the BAPS5 lineage described in this study. The acquisition of resistance against recommended treatments for shigellosis and other STIs is a global public health concern, and potentially herald an increased risk of treatment failure³⁵, and eventual emergence of pan-drug resistant *S. sonnei*.

The drivers and trajectory of this emergent XDR lineage, which is additionally ceftriaxone resistant compared with previously described internationally disseminated MSM-associated *Shigella*, are likely to be similar to previously identified phenomena. Specifically, evolution of resistance to third-generation cephalosporins is likely being driven by bystander resistance resulting from shifting treatment regimens for other sexually transmissible illnesses. While historical use of azithromycin for treating many STIs led to the acquisition of *mph*(A) and *erm*(B) among MSM associated *Shigella*¹⁸, the current first-line treatment recommendation for gonorrhoea in the USA, UK, France, and Belgium is now third-generation cephalosporins (e.g. ceftriaxone)³⁶⁻³⁹. Treatment for ongoing high rates of gonorrhoea in the UK⁴⁰ thus likely

<figure><img src="image_7.png" /><figcaption>**Fig. 4 | Relative fitness of different *Shigella* MDR plasmids in *E. coli*.** Relative fitness of *E. coli* MG1655 carrying a previously characterised azithromycin resistance plasmid (pAPR100, green; pKSR100, purple), or p1538171_3 (navy) compared with a plasmid-free strain (grey). Each dataset is represented by a box and whisker plot indicating maximum and minimum values (upper and lower limits, respectively), and interquartile ranges (25th to 75th percentiles) of the original data, overlaid with black dots, indicating individual data points. The centre emboldened line on each box represents the median value. The results of pairwise statistical comparisons are shown above and with two-tailed *p*-values (two sample t-tests) indicated above black lines joining the pair comparators. The sample size for statistical purposes is *n* = 3 biologically independent samples (see methods).</figcaption></figure>
contributed to the acquisition of resistance to third-generation cephalosporins into this circulating ciprofloxacin resistant lineage of *S. sonnei*. It is also possible that plasmid fitness cost played a role in the emergence of this XDR strain. Given the global transmission of the low fitness cost MDR pKSR100-type plasmids among various serotypes of *Shigella*¹⁹, it is likely that the ESBL producing p983816-like plasmids will continue to spread via horizontal gene transfer to other species and serotypes of *Shigella* in MSM, and early evidence of this is emerging in the United Kingdom. (K. Thorley, personal communication 2023).

The ongoing emergence of highly MDR and XDR *Shigella* sub-lineages conferred by low-fitness plasmids among MSM highlights important future directions for the genomic surveillance of AMR and sexually transmissible enteric infections (STEIs) in this patient group, and beyond. Firstly, given the evidence for antimicrobial selection pressures driving the evolution of AMR in STEIs, the impact of STI treatments on AMR could be monitored more holistically through active surveillance of the microbiota in this risk group (e.g. through detection of AMR genes in stool and/or other sentinel organisms). Active surveillance like this could feed into evidence-based updates of treatment guidelines for STIs to also monitor and mitigate the further development of XDR, and PDR strains of STEIs. Secondly, collaborative genomic epidemiology studies such as ours highlight that the emer-gence of a concerning resistance phenotype in MSM associated *Shi-gella* in one country may lead to rapid global dissemination within a short time frame (i.e., a decade) and efforts to develop, coordinate, and harmonise, global genomic surveillance in line with the new WHO strategy in this area⁴¹ are imperative to counter the public health threat posed by AMR STEIs. The final lesson from our study stems from this XDR outbreak lineage evolving from a clade that dominated MSM-associated shigellosis for nearly a decade, and a plasmid detected as early as 2018, highlighting a window of opportunity to have pre-empted this new AMR threat. To ensure global genomic surveillance efforts can eventually transition to a prospective, predictive frame-work where preventative action can be taken, and intervention efforts focused, we must work toward identifying early indicators for AMR threats with a high potential for global spread.

## Methods

### Reference isolates

Representatives of all *S. sonnei* subclades, including 120 reference isolates representative of a global genotyping framework (accessions in Supplementary Data 1)¹⁴ were included in the core genome Multi-Locus Sequence Typing (cgMLST) tree in this study to contextualise the UK outbreak.

### Isolates and data from the United Kingdom Health Security Agency (UKHSA)

All *S. sonnei* isolates with an assigned Single Nucleotide Polymorphism (SNP) address in England (97%, *n* = 2895 of 2982 submitted) received by the United Kingdom national reference laboratory (Gastrointestinal Bacteria Reference Unit) between 01/01/2016 and 30/12/2021 were included in this study. The SNP address is a linkage-clustering system in which reference-mapping based pairwise distances at 250, 100, 50, 25, 10, 5, and 0 SNP are numbered and listed in this order, joined by periods⁴².

Data on patient age was available for all isolates, data on patient sex for 99% (n = 2875/2895), and explicitly declared information about travel for approximately half (n = 1487/2895, 51%). In total, 429 isolates were from patients who reported recent travel to Asia, 289 to Africa, 73 to Europe, 121 to North America, 6 to Oceania, 38 to South America, and 73 to an unknown location. There were 458 isolates from patients where 'no travel' was explicitly indicated and for 1408 isolates information from patients was missing. In the absence of individual patient risk factors on sexual activity being available, possible MSM associated

isolates were classified by proxy as deriving from men aged 16 years or older, without a history of recent international travel, in line with a previous validation of this approach in this setting¹⁴.

All isolates were routinely sequenced (Illumina) according to previously described protocols⁴⁴. Isolates were assigned SNP addres-ses, according to in house subtyping methods⁴². The SNP address of relevance here is 1.1.1.1.377. This lineage became XDR and the outbreak was investigated at the level of the 10-SNP threshold, so the 1.1.1.1.377 SNP linkage cluster is referred to here as t10.377 for consistency with the initial report²⁶. Among the 2895 *S. sonnei* isolates in this study, 483 belonged to SNP linkage cluster t10.377.

To enable exploring international links to the UK outbreak, we employed a recently described global genotyping framework¹⁴. This indicated that the UK-based SNP linkage cluster t10.377 roughly approximated the international 3.6.1.1.2 genotype of *S. sonnei* (96% concordance statistic, Supplementary Data 1). The human-readable nomenclature for this genotype is CipR.MSM5, due to its known ciprofloxacin resistance and previous description in MSM as 'Clade 5', the prevailing *S. sonnei* subtype in the UK⁶,¹⁴.

Eight plasmid sequences from UK isolates (see accessions in Supplementary Data 1) generated from hybrid assemblies with Oxford Nanopore technology generated as previously⁴⁵, were used in plasmid comparisons.

### International outbreak associated isolates

CipR.MSM5 *S. sonnei* genotype isolates, with or without *bla*<sub>CTX-M-27</sub>, from public health surveillance organisations in other countries were integrated into the analysis (based on data volume and sharing arrangements). A brief overview for each country and aggregate data on patient age and gender, which were consistent with isolates deriv-ing from MSM-associated transmission, is provided here.

**Belgium.** All *S. sonnei* confirmed as CipR.MSM5 isolated at Sciensano's reference laboratory between January 2017 and March 2022 and with sequence data were included. All (n = 112) isolates came from patients aged up to 84 years who were mostly male (n = 101/112, 90%). Sequence data was obtained using the following protocol: isolates were cultured overnight in BHI broth (BD) at 37 °C. DNA was extracted using an MgC Bacterial DNA Kit™ with a 60 µL elution volume (Atrida, Amersfoort, The Netherlands), following the manufacturer's instructions. Sequen-cing libraries were constructed using the Illumina Nextera XT DNA sample preparation kit and sequenced on an Illumina MiSeq instru-ment with a 250 bp paired-end protocol (MiSeq v3 chemistry), according to the manufacturer's instructions (Illumina, San Diego, CA, USA).

**France.** All CipR.MSM5 isolates that contained *bla*<sub>CTX-M-27</sub> (n = 101) isolated at the Enteric Bacterial Pathogens Unit at the Institut Pasteur between 01/09/2020 and 15/02/2022 were included in this study. The isolates came from patients from 13 to 68 years and the majority were male (n = 98/101, 97%). These were sequenced using the following protocol: Genomic DNA was extracted with the MagNA Pure DNA isolation kit (Roche Molecular Systems, Indianapolis, IN, USA), in accordance with the manufacturer's instructions. Whole-genome sequencing was performed with the Mutualized Platform for Micro-biology (P2M) at Institut Pasteur, Paris. The libraries were prepared with the Nextera XT kit (Illumina, San Diego, CA, USA) and sequencing was performed with the NextSeq 500 system.

The sequence of a single plasmid (p202008564-6) from isolate 202008564 was obtained. Genomic DNA was prepared as follows: the isolates were cultured overnight at 37 °C in alkaline nutrient agar (20 g casein meat peptone E2 from Organotechnie; 5 g sodium chloride from Sigma; 15 g Bacto agar from Difco; distilled water to 1 L; adjusted to pH 8.4; autoclaved at 121 °C for 15 min). A few isolated colonies from the overnight culture were used to inoculate 20 mL of Brain-Heart-Infusion
contributed to the acquisition of resistance to third-generation cephalosporins into this circulating ciprofloxacin resistant lineage of S. sonnei. It is also possible that plasmid fitness cost played a role in the emergence of this XDR strain. Given the global transmission of the low fitness cost MDR pKSR100-type plasmids among various serotypes of Shigella¹⁹, it is likely that the ESBL producing p983816-like plasmids will continue to spread via horizontal gene transfer to other species and serotypes of Shigella in MSM, and early evidence of this is emerging in the United Kingdom. (K. Thorley, personal communication 2023).

The ongoing emergence of highly MDR and XDR *Shigella* sub-lineages conferred by low-fitness plasmids among MSM highlights important future directions for the genomic surveillance of AMR and sexually transmissible enteric infections (STEIs) in this patient group, and beyond. Firstly, given the evidence for antimicrobial selection pressures driving the evolution of AMR in STEIs, the impact of STI treatments on AMR could be monitored more holistically through active surveillance of the microbiota in this risk group (e.g. through detection of AMR genes in stool and/or other sentinel organisms). Active surveillance like this could feed into evidence-based updates of treatment guidelines for STIs to also monitor and mitigate the further development of XDR, and PDR strains of STEIs. Secondly, collaborative genomic epidemiology studies such as ours highlight that the emer-gence of a concerning resistance phenotype in MSM associated *Shi-gella* in one country may lead to rapid global dissemination within a short time frame (i.e., a decade) and efforts to develop, coordinate, and harmonise, global genomic surveillance in line with the new WHO strategy in this area⁴¹ are imperative to counter the public health threat posed by AMR STEIs. The final lesson from our study stems from this XDR outbreak lineage evolving from a clade that dominated MSM-associated shigellosis for nearly a decade, and a plasmid detected as early as 2018, highlighting a window of opportunity to have pre-empted this new AMR threat. To ensure global genomic surveillance efforts can eventually transition to a prospective, predictive frame-work where preventative action can be taken, and intervention efforts focused, we must work toward identifying early indicators for AMR threats with a high potential for global spread.

## Methods

### Reference isolates

Representatives of all S. sonnei subclades, including 120 reference isolates representative of a global genotyping framework (accessions in Supplementary Data 1)¹⁴ were included in the core genome Multi-Locus Sequence Typing (cgMLST) tree in this study to contextualise the UK outbreak.

### Isolates and data from the United Kingdom Health Security Agency (UKHSA)

All S. sonnei isolates with an assigned Single Nucleotide Polymorphism (SNP) address in England (97%, n = 2895 of 2982 submitted) received by the United Kingdom national reference laboratory (Gastrointestinal Bacteria Reference Unit) between 01/01/2016 and 30/12/2021 were included in this study. The SNP address is a linkage-clustering system in which reference-mapping based pairwise distances at 250, 100, 50, 25, 10, 5, and 0 SNP are numbered and listed in this order, joined by periods⁴².

Data on patient age was available for all isolates, data on patient sex for 99% (n = 2875/2895), and explicitly declared information about travel for approximately half (n = 1487/2895, 51%). In total, 429 isolates were from patients who reported recent travel to Asia, 289 to Africa, 73 to Europe, 121 to North America, 6 to Oceania, 38 to South America, and 73 to an unknown location. There were 458 isolates from patients where 'no travel' was explicitly indicated and for 1408 isolates information from patients was missing. In the absence of individual patient risk factors on sexual activity being available, possible MSM associated

isolates were classified by proxy as deriving from men aged 16 years or older, without a history of recent international travel, in line with a previous validation of this approach in this setting¹³.

All isolates were routinely sequenced (Illumina) according to previously described protocols⁴⁴. Isolates were assigned SNP addres-ses, according to in house subtyping methods⁴². The SNP address of relevance here is 1.1.1.1.377. This lineage became XDR and the outbreak was investigated at the level of the 10-SNP threshold, so the 1.1.1.1.377 SNP linkage cluster is referred to here as t10.377 for consistency with the initial report²⁶. Among the 2895 S. sonnei isolates in this study, 483 belonged to SNP linkage cluster t10.377.

To enable exploring international links to the UK outbreak, we employed a recently described global genotyping framework¹⁴. This indicated that the UK-based SNP linkage cluster t10.377 roughly approximated the international 3.6.1.1.2 genotype of S. sonnei (96% concordance statistic, Supplementary Data 1). The human-readable nomenclature for this genotype is CipR.MSM5, due to its known ciprofloxacin resistance and previous description in MSM as 'Clade 5', the prevailing S. sonnei subtype in the UK⁶,¹⁴.

Eight plasmid sequences from UK isolates (see accessions in Supplementary Data 1) generated from hybrid assemblies with Oxford Nanopore technology generated as previously⁴⁵, were used in plasmid comparisons.

### International outbreak associated isolates

CipR.MSM5 S. sonnei genotype isolates, with or without bla<sub>CTX-M-27</sub>, from public health surveillance organisations in other countries were integrated into the analysis (based on data volume and sharing arrangements). A brief overview for each country and aggregate data on patient age and gender, which were consistent with isolates deriv-ing from MSM-associated transmission, is provided here.

**Belgium.** All S. sonnei confirmed as CipR.MSM5 isolated at Sciensano's reference laboratory between January 2017 and March 2022 and with sequence data were included. All (n = 112) isolates came from patients aged up to 84 years who were mostly male (n = 101/112, 90%). Sequence data was obtained using the following protocol: isolates were cultured overnight in BHI broth (BD) at 37 °C. DNA was extracted using an MgC Bacterial DNA Kit™ with a 60 µL elution volume (Atrida, Amersfoort, The Netherlands), following the manufacturer's instructions. Sequen-cing libraries were constructed using the Illumina Nextera XT DNA sample preparation kit and sequenced on an Illumina MiSeq instru-ment with a 250 bp paired-end protocol (MiSeq v3 chemistry), according to the manufacturer's instructions (Illumina, San Diego, CA, USA).

**France.** All CipR.MSM5 isolates that contained bla<sub>CTX-M-27</sub> (n = 101) isolated at the Enteric Bacterial Pathogens Unit at the Institut Pasteur between 01/09/2020 and 15/02/2022 were included in this study. The isolates came from patients from 13 to 68 years and the majority were male (n = 98/101, 97%). These were sequenced using the following protocol: Genomic DNA was extracted with the MagNA Pure DNA isolation kit (Roche Molecular Systems, Indianapolis, IN, USA), in accordance with the manufacturer's instructions. Whole-genome sequencing was performed with the Mutualized Platform for Micro-biology (P2M) at Institut Pasteur, Paris. The libraries were prepared with the Nextera XT kit (Illumina, San Diego, CA, USA) and sequencing was performed with the NextSeq 500 system.

The sequence of a single plasmid (p202008564-6) from isolate 202008564 was obtained. Genomic DNA was prepared as follows: the isolates were cultured overnight at 37 °C in alkaline nutrient agar (20 g casein meat peptone E2 from Organotechnie; 5 g sodium chloride from Sigma; 15 g Bacto agar from Difco; distilled water to 1 L; adjusted to pH 8.4; autoclaved at 121 °C for 15 min). A few isolated colonies from the overnight culture were used to inoculate 20 mL of Brain-Heart-Infusion
(BHI) broth and were cultured until a final OD600 of 0.8 was reached at 37 °C with shaking (200 rpm–Thermo Scientific MaxQ 6800). Bacterial cells were harvested by centrifugation and the DNA extraction was performed by two different methods. For isolate 202008118, we followed the protocol described by von Mentzer et al.⁴⁶, except that MaXtract High Density columns (Qiagen) were used (instead of phase lock tubes) and DNA was resuspended in molecular biology grade water (instead of 10 mM Tris pH 8.0). We used a Genomic-tip 100/G column (Qiagen) according to the manufacturer's protocol. The library was prepared according to the instructions of the Native barcoding genomic DNA (with EXP-NBD104, EXP-NBD114, and SQK-LSK109) procedure provided by Oxford Nanopore Technology. Sequencing was then performed on a MinION Mk1C apparatus (Oxford Nanopore Technologies). The genomic sequences of the isolates were assembled from long and short reads, with a hybrid approach and Unicycer (v 0.4.8)⁴⁷. A polishing step was performed with Pilon (v 1.23), to generate a high-quality plasmid sequence⁴⁸, p202008564-6 (GenBank accession no. OP038290).

**United States of America.** Data for all publicly available CipR.MSM5 isolates which contained a gene in the *bla*<sub>CTX-M</sub> family (n = 31) and associated age and gender data were kindly provided by PulseNet, Centre for Diseases Control. Aggregate metadata showed that these isolates came from patients aged between 20 and 79 years that were predominately male (n = 30/31, 97%).

**Australia.** From the state of Victoria, sequences of CipR.MSM5 *bla*<sub>CTX-M-27</sub> isolates (n = 3) from routine health surveillance (Microbiological Diagnostic Unit, Public Health Laboratory, Doherty Institute) were analysed, two of which have been described previously (AUSMDU00040532 and AUSMDU00040564)²³. All three were from male patients aged between 16 and 60 years. From the state of New South Wales, all CipR.MSM5 isolates collected during routine surveillance October 2017-November 2020 (n = 232) were included. Isolates were from patients between 19 and 83 years (mean = 40 years) and mainly male (n = 224/232, 97%). Genomic DNA was extracted from pure cultures using a QIAGEN DNeasy Blood and Tissue Mini Kit (QIAGEN). Sequencing libraries were prepared using Nextera XT DNA Library Prep Kit (Illumina) and sequenced on a NextSeq500 at the Microbial Genomics Reference Laboratory, Westmead Hospital. Isolate 20-001-0088 was used for long-read sequencing. DNA was extracted from growth on Horse Blood agar (Thermofisher) using DNeasy® Ultra-Clean® Microbial Kit (Qiagen) with mechanical lysis reduced to 2 min. DNA quality and quantity were assessed on a Nano-300 Micro-Spectrophotometer (Allsheng) and a Qubit™ 2 Fluorometer (Life Technologies) while integrity was assessed by 0.6% (w/v) agarose gel electrophoresis. A library was prepared using the SQK-RBK004 Rapid barcoding kit (Oxford Nanopore Technologies) according to manufacturer's instructions, loaded into a R4.9.1 flow cell and sequenced on a MinION Mk 1B (Oxford Nanopore Technologies). Sequencing was performed for up to 24 h and base-calling was performed post sequencing using Guppy (v 3.4.5 + flfbbf)⁴⁹. Hybrid assembly was performed using Unicycer (v 0.4.8)⁴⁷. The assembly was further short-read polished on Pilon (v 1.23)⁴⁸, using BAM files generated by MiniMap2 (v 2.17-r941)⁵⁰, until no further changes could be made. The assembly was automatically annotated using Prokka (v 1.14.6)⁵¹, running on metagenome mode and using genetic code 11.

## Bioinformatic analyses

**Clustering dendrogram of UK isolates.** A dendrogram was constructed for the UK isolates available on Enterobase (n = 2820/2895) and the global context representatives (n = 120) using the NINJA NJ tree algorithm (v 1.0) on the cgMLST hierarchical clustering results from the algorithm implemented in Enterobase (v 1.1.3)⁵²,⁵³. The tree was visualised using the Interactive Tree of Life (iTOL, v 5)⁵⁴.

**Phylogenetic reconstruction of the international outbreak.** To explore international relationships, we constructed a tree from international CipR.MSM5 isolates (n = 475) alongside isolates belonging to the UK outbreak cluster (n = 468 of 482 which generated pseudogenomes of uniform length, t10.377)²⁶. Illumina paired end reads for the isolates from the UK were downloaded from the Sequence Read Archive (SRA) using fastq-dump (SRA toolkit v 2.11.0)⁵⁵ with default settings except the '--split 3' parameter was used to automatically separate paired end read data into two separate FASTQ files. Nucleotide sequences were trimmed using Trimmomatic (v 0.39)⁵⁶ with the parameters: Leading: 20, Trailing: 20, SlidingWindow:4:20, MinLen:40, and quality checked with MultiQC (v 1.12)⁵⁷. Burrows-Wheeler Aligner (BWA v 0.7.17) was used with default parameters to align the S. sonnei sequences to *Shigella sonnei* 53G as the reference genome (NCBI accession number: GCA_00283715.1:ASM28371vi)⁵⁸. PICARD (v 2.27.2)⁵⁹ was used to mark and remove artificial nucleotide duplicates from the genomes of the S. sonnei isolates using default parameters. SAMTOOLS (v 1.11) was used to index the files⁶⁰. BCFTOOLS (v 1.9)⁶¹, was used to call and filter variants, using default parameters and '-c --ploidy 1 -O -o' for variant calling and '-SnpGap 15 --IndelGap 15' for filtering⁶². A chromosomal pseudogenome was created for all isolates and prophages and plasmids were masked using sed (v 4.2.2)⁶³. Gubbins (v 3.2.1)⁶⁴ was then used to remove remaining regions of recombination using default parameters except: '-c 30, f 60'. Multiple sequence alignments were then used to 1) impute a phylogenetic tree using IQTREE (v 2.2.0.3)⁶⁵ with the default parameters except: '-ntmax 25 -bb 1000 -m GTR+I+R+ASC, and 2) determine population clusters using Rhier Bayesian Analysis of Genetic Population Structure (RhierBAPS v 1.01)⁶⁶,⁶⁷.

**Temporal signal analyses.** Newick tree files from iTOL (v 5)⁵⁴, derived from the IQTREE (v 2.2.0.3)⁵⁶ output described above for isolates where the full dates were known (day, month and year, UK, n = 555, NSW n = 232), were imported into TempEst⁶⁸, (only the year of isolate collection is listed in Supplementary Data 1, for patient confidentiality reasons). Default parameters were used, with the function set to 'heuristic residual mean squared', and with 'best-fitting root' activated. Dates were parsed as a calendar date (DD/MM/YYYY). Dates were specified as 'years' and 'since some time in the past' and lineage 'root-to-tip' data used to infer the date of emergence of the most recent common ancestors.

**Genome assembly and identifying genes of interest.** Draft genomes were assembled for isolates from all locations using Unicycer (v 0.5.0)⁴⁷, using the 'forward paired', 'reverse paired', 'forward unpaired', 'reverse unpaired' trimmed and quality controlled FASTQ files. Draft genome assemblies were checked for quality using QUAST (v 5.02)⁶⁹. Default parameters were used except *min fasta length 200, -vcf* were used instead. Draft genomes were then analysed for AMR genes using NCBI AMRFinder Plus (v 3.10.24)⁷⁰ with the organism set as *Escherichia*, and the '--plus' parameter to obtain information on stress genes. SonneiTyping Script (v 20210201)⁷¹ embedded in Mykrobe (v 0.11.0)⁷² was used for genotyping and retrieval of QRDR mutations in *gyrA* and *parC* to infer ciprofloxacin resistance⁷³.

**Statistical testing.** Standard proportions were calculated, and statistical support for phylogenetic analyses evaluated using chi-square testing. Plasmid fitness experiments were evaluated by two sample *t* test.

**In silico plasmid analyses.** The eight UK IncFII plasmids were reorientated to start with the *repA* gene using circulator Fixstart (v 1.5.5)⁷⁴, followed by annotation using Prokka (v 1.14.6),⁷¹. Plasmids (GBK files) were compared (>95% nucleotide identity) and visualised using Clinker (v 0.0.21)⁷⁵. For a broader search, the sequence of p893916 was queried against a Compact Bit-sliced Signature (COBS)
(BHI) broth and were cultured until a final OD₆₀₀ of 0.8 was reached at 37 °C with shaking (200 rpm–Thermo Scientific MaxQ 6800). Bacterial cells were harvested by centrifugation and the DNA extraction was performed by two different methods. For isolate 202008118, we followed the protocol described by von Mentzer et al.⁴⁶, except that MaXtract High Density columns (Qiagen) were used (instead of phase lock tubes) and DNA was resuspended in molecular biology grade water (instead of 10 mM Tris pH 8.0). We used a Genomic-tip 100/G column (Qiagen) according to the manufacturer's protocol. The library was prepared according to the instructions of the Native barcoding genomic DNA (with EXP-NBD104, EXP-NBD114, and SQK-LSK109) procedure provided by Oxford Nanopore Technology. Sequencing was then performed on a MinION Mk1C apparatus (Oxford Nanopore Technologies). The genomic sequences of the isolates were assembled from long and short reads, with a hybrid approach and Unicycer (v 0.4.8)⁴⁷. A polishing step was performed with Pilon (v 1.23), to generate a high-quality plasmid sequence⁴⁸, p202008564-6 (GenBank accession no. OP038290).

**United States of America.** Data for all publicly available CipR.MSM5 isolates which contained a gene in the *bla*<sub>CTX-M</sub> family (n = 31) and associated age and gender data were kindly provided by PulseNet, Centre for Diseases Control. Aggregate metadata showed that these isolates came from patients aged between 20 and 79 years that were predominately male (n = 30/31, 97%).

**Australia.** From the state of Victoria, sequences of CipR.MSM5 *bla*<sub>CTX-M-27</sub> isolates (n = 3) from routine health surveillance (Microbiological Diagnostic Unit, Public Health Laboratory, Doherty Institute) were analysed, two of which have been described previously (AUSMDU00040532 and AUSMDU00040564)²³. All three were from male patients aged between 16 and 60 years. From the state of New South Wales, all CipR.MSM5 isolates collected during routine surveillance October 2017-November 2020 (n = 232) were included. Isolates were from patients between 19 and 83 years (mean = 40 years) and mainly male (n = 224/232, 97%). Genomic DNA was extracted from pure cultures using a QIAGEN DNeasy Blood and Tissue Mini Kit (QIAGEN). Sequencing libraries were prepared using Nextera XT DNA Library Prep Kit (Illumina) and sequenced on a NextSeq500 at the Microbial Genomics Reference Laboratory, Westmead Hospital. Isolate 20-001-0088 was used for long-read sequencing. DNA was extracted from growth on Horse Blood agar (Thermofisher) using DNeasy® Ultra-Clean® Microbial Kit (Qiagen) with mechanical lysis reduced to 2 min. DNA quality and quantity were assessed on a Nano-300 Micro-Spectrophotometer (Allsheng) and a Qubit™ 2 Fluorometer (Life Technologies) while integrity was assessed by 0.6% (w/v) agarose gel electrophoresis. A library was prepared using the SQK-RBK004 Rapid barcoding kit (Oxford Nanopore Technologies) according to manufacturer's instructions, loaded into a R4.9.1 flow cell and sequenced on a MinION Mk 1B (Oxford Nanopore Technologies). Sequencing was performed for up to 24 h and base-calling was performed post sequencing using Guppy (v 3.4.5 + flfbbf)⁴⁹. Hybrid assembly was performed using Unicycer (v 0.4.8)⁴⁷. The assembly was further short-read polished on Pilon (v 1.23)⁴⁸, using BAM files generated by MiniMap2 (v 2.17-r941)⁵⁰, until no further changes could be made. The assembly was automatically annotated using Prokka (v 1.14.6)⁵¹, running on metagenome mode and using genetic code 11.

## Bioinformatic analyses

**Clustering dendrogram of UK isolates.** A dendrogram was constructed for the UK isolates available on Enterobase (n = 2820/2895) and the global context representatives (n = 120) using the NINJA NJ tree algorithm (v 1.0) on the cgMLST hierarchical clustering results from the algorithm implemented in Enterobase (v 1.1.3)⁵²,⁵³. The tree was visualised using the Interactive Tree of Life (iTOL, v 5)⁵⁴.

**Phylogenetic reconstruction of the international outbreak.** To explore international relationships, we constructed a tree from international CipR.MSM5 isolates (n = 475) alongside isolates belonging to the UK outbreak cluster (n = 468 of 482 which generated pseudogenomes of uniform length, t10.377)²⁶. Illumina paired end reads for the isolates from the UK were downloaded from the Sequence Read Archive (SRA) using fastq-dump (SRA toolkit v 2.11.0)⁵⁵ with default settings except the '--split 3' parameter was used to automatically separate paired end read data into two separate FASTQ files. Nucleotide sequences were trimmed using Trimmomatic (v 0.39)⁵⁶ with the parameters: Leading: 20, Trailing: 20, SlidingWindow:4:20, MinLen:40, and quality checked with MultiQC (v 1.12)⁵⁷. Burrows-Wheeler Aligner (BWA v 0.7.17) was used with default parameters to align the S. sonnei sequences to *Shigella sonnei* 53G as the reference genome (NCBI accession number: GCA_00283715.1:ASM28371vi)⁵⁸. PICARD (v 2.27.2)⁵⁹ was used to mark and remove artificial nucleotide duplicates from the genomes of the S. sonnei isolates using default parameters. SAMTOOLS (v 1.11) was used to index the files⁶⁰. BCFTOOLS (v 1.9)⁶¹, was used to call and filter variants, using default parameters and '-c --ploidy 1 -O -o' for variant calling and '-SnpGap 15 --IndelGap 15' for filtering⁶². A chromosomal pseudogenome was created for all isolates and prophages and plasmids were masked using sed (v 4.2.2)⁶³. Gubbins (v 3.2.1)⁶⁴ was then used to remove remaining regions of recombination using default parameters except: '-c 30, f 60'. Multiple sequence alignments were then used to 1) impute a phylogenetic tree using IQTREE (v 2.2.0.3)⁶⁵ with the default parameters except: '-ntmax 25 -bb 1000 -m GTR+I+R+ASC, and 2) determine population clusters using Rhier Bayesian Analysis of Genetic Population Structure (RhierBAPS v 1.01)⁶⁶,⁶⁷.

**Temporal signal analyses.** Newick tree files from iTOL (v 5)⁵⁴, derived from the IQTREE (v 2.2.0.3)⁵⁶ output described above for isolates where the full dates were known (day, month and year, UK, n = 555, NSW n = 232), were imported into TempEst⁶⁸, (only the year of isolate collection is listed in Supplementary Data 1, for patient confidentiality reasons). Default parameters were used, with the function set to 'heuristic residual mean squared', and with 'best-fitting root' activated. Dates were parsed as a calendar date (DD/MM/YYYY). Dates were specified as 'years' and 'since some time in the past' and lineage 'root-to-tip' data used to infer the date of emergence of the most recent common ancestors.

**Genome assembly and identifying genes of interest.** Draft genomes were assembled for isolates from all locations using Unicycer (v 0.5.0)⁴⁷, using the 'forward paired', 'reverse paired', 'forward unpaired', 'reverse unpaired' trimmed and quality controlled FASTQ files. Draft genome assemblies were checked for quality using QUAST (v 5.02)⁶⁹. Default parameters were used except *'min fasta length 200, -vcf'* were used instead. Draft genomes were then analysed for AMR genes using NCBI AMRFinder Plus (v 3.10.24)⁷⁰ with the organism set as *Escherichia*, and the '--plus' parameter to obtain information on stress genes. SonneiTyping Script (v 20210201)⁷¹ embedded in Mykrobe (v 0.11.0)⁷² was used for genotyping and retrieval of QRDR mutations in *gyrA* and *parC* to infer ciprofloxacin resistance⁷³.

**Statistical testing.** Standard proportions were calculated, and statistical support for phylogenetic analyses evaluated using chi-square testing. Plasmid fitness experiments were evaluated by two sample *t* test.

**In silico plasmid analyses.** The eight UK IncFII plasmids were reorientated to start with the *repA* gene using circulator Fixstart (v 1.5.5)⁷⁴, followed by annotation using Prokka (v 1.14.6),⁷¹. Plasmids (GBK files) were compared (>95% nucleotide identity) and visualised using Clinker (v 0.0.21)⁷⁵. For a broader search, the sequence of p893916 was queried against a Compact Bit-sliced Signature (COBS)
index<sup>76</sup>, of 661,405 curated draft genomes<sup>77</sup>, with a k-mer similarity cut-off of 80%. Genomes that were high quality and had a species ID of 0.8 or higher (n = 1505) were included in the distribution figure. Species and sequence types (ST) were extracted from Supplementary file File4_QC_characterisation_611K<sup>77</sup>. But due to the difficulty in distinguishing *E. coli* and *Shigella* spp., some species initially identified as *E. coli* were updated to the relevant *Shigella* species based on the assigned ST. A BLAST Ring Image Generator (BRIG, v 0.95)<sup>78</sup>, plot was constructed using default settings with p893816 (GenBank accession MW396858) as the reference plasmid genome. Illumina reads for all BAPS 5 isolates were mapped against p893816 using Burrow-Wheeler Aligner (v 0.7.17-r118)<sup>78</sup>. Percentage coverage at 1X sequencing depth calculated by Qualimap (v 2.2.2)<sup>79</sup> with default parameters. Isolates known to have plasmids virtually identical to p893816 (Table 2, Fig. S2) gave values of 96-97%.

**In vitro plasmid analyses.** The p893816-like plasmid from *S. sonnei* clinical isolate 1538717 (p1538717_3, CP104412) was conjugated into *E. coli* MG1655 and used in antimicrobial susceptibility testing and relative fitness assays, all as previously described<sup>19</sup>. *E. coli* MG1655 carrying pKSR100 or pAPR100 were generated from isolates obtained from the UKHSA, as previously described<sup>19</sup>. *E. coli* MG1655 and derivatives that each carry a different plasmid were inoculated into M9 minimal media supplemented with thiamine (10 µg/ml) and grown overnight with shaking at 220 rpm at 37 °C. These pre-cultures were adjusted to OD<sub>600</sub> 0.05 in a total volume of 1 ml and 150 µl aliquots inoculated into individual wells of a 96-well plate (Greiner Bio One, UK). The 96-well plate was then incubated with shaking at 37 °C in a Synergy H1 multi-mode plate reader (BioTek Instruments) taking optical density measurements at 600 nm every 15 min. Data were analysed using the R package Growthcurver and Area Under the Curve (AUC) output was used to calculate the relative fitness of each strain carrying the plasmid compared with the plasmid-free strain<sup>80</sup>. Each experiment consisted of three technical replicates and was repeated three times (i.e. three biological replicates). Details of all strains and plasmids used in this study can be found in Supplementary Table 3.

**Antimicrobial susceptibility testing.** To confirm resistances inferred from sequence data, 14 UK isolates with representative genotypic AMR profiles, and *E. coli* MG1655 with and without p1538717_3 underwent antimicrobial susceptibility testing using Liofilchem® gradient MIC Test Strips (MTS) according to manufacturer's instructions to determine the minimum inhibitory concentrations (MICs). The following antimicrobial and MIC ranges (all in µg/mL) were used: ciprofloxacin (CIP, 0.002–32), ertapenem (ERT, 0.002–32), mecillinam (MEC, 0.016–256), azithromycin (AZM, 0.016–256), ceftriaxone (CRO, 0.002–32), trimethoprim sulfamethoxazole (SXT, 0.002–32), gentamicin (GEN, 0.016–256). Results were interpreted as defined by the European Committee on Antimicrobial Susceptibility Testing (EUCAST)<sup>81</sup>.

**Bioinformatics packages and services used.** A summary of all bioinformatics packages and services used in this study, alongside information about commands and versions used, can be found in Supplementary Table 4 and the Reporting Summary document.

**Ethics statement**

No individual patient consent was required or sought as UKHSA has authority to handle patient data for public health monitoring and infection control under section 251 of the UK National Health Service Act of 2006 (previously section 60 of the Health and Social Care Act of 2001).

**Reporting summary**

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

**Data availability**

Sequence data relating to all isolates has been deposited in the Sequence Read Archive under BioProject numbers PRJNA315192 (United Kingdom), PRJEB44801 (France), PRJEB40097 (Belgium), and PRJNA613115 (Australia, New South Wales). Individual isolate accession numbers and isolate metadata are available in the Supplementary Data 1 as indicated in the text. Correspondence regarding isolates from the United Kingdom and the overall study should be directed to Dr. Claire Jenkins and Professor Kate Baker. Correspondence regarding isolates: from France should be directed to Professor François-Xavier Weill; from Belgium should be directed to Pieter-Jan Ceyssens; from the United States of America should be directed to Kaitlin Tagg or PulseNet; and from Australia to Ben Howden and Danielle Ingle (Victoria) or Elena Martinez or Jonathan Iredell (New South Wales).

**Code availability**

No custom code was written to conduct this study. Details of bioinformatic packages and commands used can be found in Supplementary Table 4.

**References**

1. Kotloff, K. L., Riddle, M. S., Platts-Mills, J. A., Pavlinac, P. & Zaidi, A. K. *M. Shigellosis*. *Lancet* **391**, 801–812 (2018).

2. Baker, S. & The, H. C. Recent insights into *Shigella*. *Curr. Opin. Infect. Dis.* **31**, 449–454 (2018).

3. Khalil, I., Troeger, C. E., Blacker, B. F. & Reiner, R. C. Capturing the true burden of *Shigella* and ETEC: The way forward. *Vaccine* **37**, 4784–4786 (2019).

4. Warren, B. R., Parish, M. E. & Schneider, K. R. *Shigella* as a foodborne pathogen and current methods for detection in food. *Crit. Rev. Food Sci. Nutr.* **46**, 551–567 (2006).

5. He, F. et al. *Shigellosis outbreak associated with contaminated well water in a rural elementary school: Sichuan Province, China, June 7–16, 2009*. *PLoS ONE* **7**, e47239 (2012).

6. Bardsley, M. et al. Persistent transmission of Shigellosis in England is associated with a recently emerged multidrug-resistant strain of *Shigella sonnei*. *J. Clin. Microbiol.* **58**, e01692–19 (2020).

7. Baine, W. B. et al. Waterborne shigellosis at a public school. *Am. J. Epidemiol.* **101**, 323–332 (1975).

8. Ekdahl, K. & Andersson, Y. The epidemiology of travel-associated shigellosis—regional risks, seasonality and serogroups. *J. Infect.* **51**, 222–229 (2005).

9. Baker, K. S. et al. Travel- and community-based transmission of multidrug-resistant *Shigella sonnei* lineage among International Orthodox Jewish Communities. *Emerg. Infect. Dis.* **22**, 1545–1553 (2016).

10. Simms, I. et al. Intensified shigellosis epidemic associated with sexual transmission in men who have sex with men—*Shigella flexneri* and *S. sonnei* in England, 2004 to end of February 2015. *Eurosurveillance* **20**, 21097 (2015).

11. Farrar, W. E. & Eidson, M. Antibiotic resistance in "*Shigella*" mediated by R factors. *J. Infect. Dis.* **123**, 477–484 (1971).

12. Davies, J. R., Farrant, W. N. & Tomlinson, A. J. Further studies on the antibiotic resistance of *Shigella sonnei*. *II. The acquisition of transferable antibiotic resistance in vivo*. *J. Hyg., Lond.* **66**, 479–487 (1968).

13. Holt, K. E. et al. *Shigella sonnei* genome sequencing and phylogenetic analysis indicate recent global dissemination from Europe. *Nat. Genet.* **44**, 1056–1059 (2012).

14. Hawkey, J. et al. Global population structure and genotyping framework for genomic surveillance of the major dysentery pathogen, *Shigella sonnei*. *Nat. Commun.* **12**, 2684 (2021).

15. Holt, K. E. et al. Tracking the establishment of local endemic populations of an emergent enteric pathogen. *Proc. Natl Acad. Sci.* **110**, 17522–17527 (2013).
index<sup>76</sup>, of 661,405 curated draft genomes<sup>77</sup>, with a k-mer similarity cut-off of 80%. Genomes that were high quality and had a species ID of 0.8 or higher (n = 1505) were included in the distribution figure. Species and sequence types (ST) were extracted from Supplementary file File4_QC_characterisation_611K<sup>77</sup>. But due to the difficulty in distinguishing *E. coli* and *Shigella* spp., some species initially identified as *E. coli* were updated to the relevant *Shigella* species based on the assigned ST. A BLAST Ring Image Generator (BRIG, v 0.95)<sup>78</sup>, plot was constructed using default settings with p893816 (GenBank accession MW396858) as the reference plasmid genome. Illumina reads for all BAPS 5 isolates were mapped against p893816 using Burrow-Wheeler Aligner (v 0.7.17-r118)<sup>78</sup>. Percentage coverage at 1X sequencing depth calculated by Qualimap (v 2.2.2)<sup>79</sup> with default parameters. Isolates known to have plasmids virtually identical to p893816 (Table 2, Fig. S2) gave values of 96-97%.

**In vitro plasmid analyses.** The p893816-like plasmid from *S. sonnei* clinical isolate 1538717 (p1538717_3, CP104412) was conjugated into *E. coli* MG1655 and used in antimicrobial susceptibility testing and relative fitness assays, all as previously described<sup>19</sup>. *E. coli* MG1655 carrying pKSR100 or pAPR100 were generated from isolates obtained from the UKHSA, as previously described<sup>19</sup>. *E. coli* MG1655 and derivatives that each carry a different plasmid were inoculated into M9 minimal media supplemented with thiamine (10 µg/ml) and grown overnight with shaking at 220 rpm at 37 °C. These pre-cultures were adjusted to OD<sub>600</sub> 0.05 in a total volume of 1 ml and 150 µl aliquots inoculated into individual wells of a 96-well plate (Greiner Bio One, UK). The 96-well plate was then incubated with shaking at 37 °C in a Synergy H1 multi-mode plate reader (BioTek Instruments) taking optical density measurements at 600 nm every 15 min. Data were analysed using the R package Growthcurver and Area Under the Curve (AUC) output was used to calculate the relative fitness of each strain carrying the plasmid compared with the plasmid-free strain<sup>80</sup>. Each experiment consisted of three technical replicates and was repeated three times (i.e. three biological replicates). Details of all strains and plasmids used in this study can be found in Supplementary Table 3.

**Antimicrobial susceptibility testing.** To confirm resistances inferred from sequence data, 14 UK isolates with representative genotypic AMR profiles, and *E. coli* MG1655 with and without p1538717_3 underwent antimicrobial susceptibility testing using Liofilchem® gradient MIC Test Strips (MTS) according to manufacturer's instructions to determine the minimum inhibitory concentrations (MICs). The following antimicrobial and MIC ranges (all in µg/mL) were used: ciprofloxacin (CIP, 0.002–32), ertapenem (ERT, 0.002–32), mecillinam (MEC, 0.016–256), azithromycin (AZM, 0.016–256), ceftriaxone (CRO, 0.002–32), trimethoprim sulfamethoxazole (SXT, 0.002–32), gentamicin (GEN, 0.016–256). Results were interpreted as defined by the European Committee on Antimicrobial Susceptibility Testing (EUCAST)<sup>81</sup>.

**Bioinformatics packages and services used.** A summary of all bioinformatics packages and services used in this study, alongside information about commands and versions used, can be found in Supplementary Table 4 and the Reporting Summary document.

**Ethics statement**

No individual patient consent was required or sought as UKHSA has authority to handle patient data for public health monitoring and infection control under section 251 of the UK National Health Service Act of 2006 (previously section 60 of the Health and Social Care Act of 2001).

**Reporting summary**

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

**Data availability**

Sequence data relating to all isolates has been deposited in the Sequence Read Archive under BioProject numbers PRJNA315192 (United Kingdom), PRJEB44801 (France), PRJEB40097 (Belgium), and PRJNA613115 (Australia, New South Wales). Individual isolate accession numbers and isolate metadata are available in the Supplementary Data 1 as indicated in the text. Correspondence regarding isolates from the United Kingdom and the overall study should be directed to Dr. Claire Jenkins and Professor Kate Baker. Correspondence regarding isolates: from France should be directed to Professor François-Xavier Weill; from Belgium should be directed to Pieter-Jan Ceyssens; from the United States of America should be directed to Kaitlin Tagg or PulseNet; and from Australia to Ben Howden and Danielle Ingle (Victoria) or Elena Martinez or Jonathan Iredell (New South Wales).

**Code availability**

No custom code was written to conduct this study. Details of bioinformatic packages and commands used can be found in Supplementary Table 4.

**References**

1. Kotloff, K. L., Riddle, M. S., Platts-Mills, J. A., Pavlinac, P. & Zaidi, A. K. *M. Shigellosis*. *Lancet* **391**, 801–812 (2018).

2. Baker, S. & The, H. C. Recent insights into *Shigella*. *Curr. Opin. Infect. Dis.* **31**, 449–454 (2018).

3. Khalil, I., Troeger, C. E., Blacker, B. F. & Reiner, R. C. Capturing the true burden of *Shigella* and ETEC: The way forward. *Vaccine* **37**, 4784–4786 (2019).

4. Warren, B. R., Parish, M. E. & Schneider, K. R. *Shigella* as a foodborne pathogen and current methods for detection in food. *Crit. Rev. Food Sci. Nutr.* **46**, 551–567 (2006).

5. He, F. et al. *Shigellosis outbreak associated with contaminated well water in a rural elementary school: Sichuan Province, China, June 7–16, 2009*. *PLoS ONE* **7**, e47239 (2012).

6. Bardsley, M. et al. Persistent transmission of Shigellosis in England is associated with a recently emerged multidrug-resistant strain of *Shigella sonnei*. *J. Clin. Microbiol.* **58**, e01692–19 (2020).

7. Baine, W. B. et al. Waterborne shigellosis at a public school. *Am. J. Epidemiol.* **101**, 323–332 (1975).

8. Ekdahl, K. & Andersson, Y. The epidemiology of travel-associated shigellosis—regional risks, seasonality and serogroups. *J. Infect.* **51**, 222–229 (2005).

9. Baker, K. S. et al. Travel- and community-based transmission of multidrug-resistant *Shigella sonnei* lineage among International Orthodox Jewish Communities. *Emerg. Infect. Dis.* **22**, 1545–1553 (2016).

10. Simms, I. et al. Intensified shigellosis epidemic associated with sexual transmission in men who have sex with men—*Shigella flexneri* and *S. sonnei* in England, 2004 to end of February 2015. *Eurosurveillance* **20**, 21097 (2015).

11. Farrar, W. E. & Eidson, M. Antibiotic resistance in "*Shigella*" mediated by R factors. *J. Infect. Dis.* **123**, 477–484 (1971).

12. Davies, J. R., Farrant, W. N. & Tomlinson, A. J. Further studies on the antibiotic resistance of *Shigella sonnei*. *II. The acquisition of transferable antibiotic resistance in vivo*. *J. Hyg., Lond.* **66**, 479–487 (1968).

13. Holt, K. E. et al. *Shigella sonnei* genome sequencing and phylogenetic analysis indicate recent global dissemination from Europe. *Nat. Genet.* **44**, 1056–1059 (2012).

14. Hawkey, J. et al. Global population structure and genotyping framework for genomic surveillance of the major dysentery pathogen, *Shigella sonnei*. *Nat. Commun.* **12**, 2684 (2021).

15. Holt, K. E. et al. Tracking the establishment of local endemic populations of an emergent enteric pathogen. *Proc. Natl Acad. Sci.* **110**, 17522–17527 (2013).
16. The, H. C. et al. South Asia as a Reservoir for the Global Spread of Ciprofloxacin-Resistant *Shigella sonnei*: A Cross-Sectional Study. *PLOS Med.* **13**, e1002055 (2016).

17. Baker, K. S. et al. Horizontal antimicrobial resistance transfer drives epidemics of multiple *Shigella* species. *Nat. Commun.* **9**, 1462 (2018).

18. Baker, K. S. et al. Intercontinental dissemination of azithromycin-resistant shigellosis through sexual transmission: a cross-sectional study. *Lancet Infect. Dis.* **15**, 913–921 (2015).

19. Malaka De Silva, P. et al. A tale of two plasmids: contributions of plasmid associated phenotypes to epidemiological success among *Shigella*. *Proc. R. Soc. B: Biol. Sci.* **289**, 20220581 (2022).

20. World Health Organisation (WHO). Guidelines for the Control of Shigellosis, Including Epidemics Due to Shigella Dysenteriae Type 1. https://apps.who.int/iris/handle/10665/43252 (Geneva, 2005).

21. World Health Organisation (WHO). WHO Publishes List of Bacteria for Which New Antibiotics Are Urgently Needed. https://www.who.int/news/item/27-02-2017-who-publishes-list-of-bacteria-for-which-new-antibiotics-are-urgently-needed (Geneva, 2017).

22. European Centre for Disease Prevention and Control. Increase in Extensively-drug Resistant Shigella sonnei Infections in Men Who Have Sex with Men in the EU/EEA and the UK, https://www.ecdc.europa.eu/en/publications-data/rapid-risk-assessment-increase-extensively-drug-resistant-shigella-sonnei (Stockholm, 2022).

23. Ingle, D. J. et al. Prolonged outbreak of multidrug-resistant *Shigella sonnei* Harboring *bla*<sub>CTX-M-27</sub> in Victoria, Australia. *Antimicrobial Agents Chemother.* **64**, e01518–20 (2020).

24. Love, N. K. et al. Impact of the COVID-19 pandemic on gastro-intestinal infection trends in England, February–July 2020. *BMJ Open* **12**, e050469 (2022).

25. Public Health England (PHE). Multi-drug Resistant Shigella sonnei Cluster (CTX-M-27) Probably Associated with MSM. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/769163/SHGLL-SONN_2018275.pdf (London, 2019).

26. Charles, H. et al. Outbreak of sexually transmitted, extensively drug-resistant *Shigella sonnei* in the UK, 2021-22: a descriptive epidemiological study. *Lancet Infect. Dis.* **22**, 1503–1510 (2022).

27. Baker, K. S. et al. Genomic epidemiology of *Shigella* in the United Kingdom shows transmission of pathogen sublineages and determinants of antimicrobial resistance. *Sci. Rep.* **8**, 7389 (2018).

28. Locke, R. K., Greig, D. R., Jenkins, C., Dallman, T. J. & Cowley, L. A. Acquisition and loss of CTX-M plasmids in *Shigella* species associated with MSM transmission in the UK. *Microb. Genomics* **7**, 000644 (2021).

29. Holt, K. E. et al. Tracking the establishment of local endemic populations of an emergent enteric pathogen. *Proc. Natl Acad. Sci. USA* **110**, 17522–17527 (2013).

30. Rowlinson, E., Berzkalns, A., Thibault, C., Golden, M. & Barbee, L. O14.3 High incidence of antimicrobial use and overuse in cisgender men who have sex with men at risk of bacterial STIs. *Sexually Transmitted Infect.* **97**, A48–A49 (2021).

31. Bengtsson Rebecca, J. et al. Accessory genome dynamics and structural variation of *Shigella* from persistent infections. *mBio* **12**, e00254–00221 (2021).

32. Thanh Duy, P. et al. Commensal *Escherichia coli* are a reservoir for the transfer of XDR plasmids into epidemic fluoroquinolone-resistant *Shigella sonnei*. *Nat. Microbiol.* **5**, 256–264 (2020).

33. Allen, H. et al. Evidence for re-infection and persistent carriage of *Shigella* species in adult males reporting domestically acquired infection in England. *Clin. Microbiol. Infect.* **27**, 126.e127–126.e113 (2021).

34. Mason, L., Dallman, T., Jenkins, C. & Baker, K. Genomic epidemiology of the first London outbreak of antimicrobial resistant sexually

transmitted shigellosis. Access *Microbiology* **4**, https://doi.org/10.1099/acmi.ac2021-po0344 (2022).

35. Collins, J. P. et al. Evidence of failure of oral third-generation cephalosporin treatment for *Shigella sonnei* infection. *Open Forum Infect. Dis.* **7**, ofaa113 (2020).

36. St Cyr, S. et al. Update to CDC's treatment guidelines for gonococcal infection, 2020. *MMWR Morbidity Mortal. Wkly. Rep.* **69**, 1911–1916 (2020).

37. Fifer, H., Saunders, J., Soni, S., Sadiq, S. T. & FitzGerald, M. 2018 UK national guideline for the management of infection with *Neisseria gonorrhoeae*. *Int. J. STD AIDS* **31**, 4–15 (2020).

38. Fouéré, S. et al. Update on French recommendations for the treatment of uncomplicated *Neisseria gonorrhoeae* infections. *Int. J. STD AIDS* **32**, 1081–1083 (2021).

39. Jespers, V. et al. Diagnosis and treatment of gonorrhoea: 2019 Belgian National guideline for primary care. *Acta Clin. Belgica* **77**, 186–194 (2022).

40. Migchelsen, S. S. et al. Sexually Transmitted Infections and Screening for Chlamydia in England, 2021 (UK Health Security Agency (UKHSA), 2022).

41. Carter, L. L. et al. Global genomic surveillance strategy for pathogens with pandemic and epidemic potential 2022–2032. *Bull. World Health Organ.* **100**, 239–239a (2022).

42. Dallman, T. J. et al. Use of whole-genome sequencing for the public health surveillance of *Shigella sonnei* in England and Wales, 2015. *J. Med. Microbiol.* **65**, 882–884 (2016).

43. Mitchell, H. D. et al. Use of whole-genome sequencing to identify clusters of *Shigella flexneri* associated with sexual transmission in men who have sex with men in England: a validation study using linked behavioural data. *Microb. Genomics* **5**, e000311 (2019).

44. Dallman, T. J. et al. Emergence of novel strains of *Shigella flexneri* associated with sexual transmission in adult men in England, 2019-2020. *J. Med. Microbiol.* **70**, 001437 (2021).

45. Yara, D. A., Greig, D. R., Gally, D. L., Dallman, T. J. & Jenkins, C. Comparison of Shiga toxin-encoding bacteriophages in highly pathogenic strains of Shiga toxin-producing *Escherichia coli* O157:H7 in the UK. *Microb. Genomics* **6**, e000334 (2020).

46. von Mentzer, A. et al. Long-read-sequenced reference genomes of the seven major lineages of enterotoxigenic *Escherichia coli* (ETEC) circulating in modern time. *Sci. Rep.* **11**, 9256 (2021).

47. Wick, R. R., Judd, L. M., Gorrie, C. L. & Holt, K. E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. *PLoS Comput. Biol.* **13**, e1005595 (2017).

48. Walker, B. J. et al. Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement. *PLoS ONE* **9**, e112963 (2014).

49. Nanopore Technologies. Guppy, https://nanoporetech.com/ (Oxford, 2022).

50. Li, H. Minimap2: pairwise alignment for nucleotide sequences. *Bioinformatics* **34**, 3094–3100 (2018).

51. Seemann, T. Prokka: rapid prokaryotic genome annotation. *Bioinformatics* **30**, 2068–2069 (2014).

52. Zhou, Z., Alikhan, N. F., Mohamed, K., Fan, Y. & Achtman, M. The EnteroBase user's guide, with case studies on *Salmonella* transmissions, *Yersinia pestis* phylogeny, and *Escherichia* core genomic diversity. *Genome Res.* **30**, 138–152 (2020).

53. Zhou, Z. et al. GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens. *Genome Res.* **28**, 1395–1404 (2018).

54. Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. *Nucleic Acids Res.* **49**, W293–W296 (2021).

55. SRAToolkitDevelopmentTeam. SRA Toolkit https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=software, (GitHub Repository, 2022).
16. The, H. C. et al. South Asia as a Reservoir for the Global Spread of Ciprofloxacin-Resistant *Shigella sonnei*: A Cross-Sectional Study. *PLOS Med.* **13**, e1002055 (2016).

17. Baker, K. S. et al. Horizontal antimicrobial resistance transfer drives epidemics of multiple *Shigella* species. *Nat. Commun.* **9**, 1462 (2018).

18. Baker, K. S. et al. Intercontinental dissemination of azithromycin-resistant shigellosis through sexual transmission: a cross-sectional study. *Lancet Infect. Dis.* **15**, 913–921 (2015).

19. Malaka De Silva, P. et al. A tale of two plasmids: contributions of plasmid associated phenotypes to epidemiological success among *Shigella*. *Proc. R. Soc. B: Biol. Sci.* **289**, 20220581 (2022).

20. World Health Organisation (WHO). Guidelines for the Control of Shigellosis, Including Epidemics Due to Shigella Dysenteriae Type 1. https://apps.who.int/iris/handle/10665/43252 (Geneva, 2005).

21. World Health Organisation (WHO). WHO Publishes List of Bacteria for Which New Antibiotics Are Urgently Needed. https://www.who.int/news/item/27-02-2017-who-publishes-list-of-bacteria-for-which-new-antibiotics-are-urgently-needed (Geneva, 2017).

22. European Centre for Disease Prevention and Control. Increase in Extensively-drug Resistant Shigella sonnei Infections in Men Who Have Sex with Men in the EU/EEA and the UK, https://www.ecdc.europa.eu/en/publications-data/rapid-risk-assessment-increase-extensively-drug-resistant-shigella-sonnei (Stockholm, 2022).

23. Ingle, D. J. et al. Prolonged outbreak of multidrug-resistant *Shigella sonnei* Harboring *bla*<sub>CTX-M-27</sub> in Victoria, Australia. *Antimicrobial Agents Chemother.* **64**, e01518–20 (2020).

24. Love, N. K. et al. Impact of the COVID-19 pandemic on gastro-intestinal infection trends in England, February–July 2020. *BMJ Open* **12**, e050469 (2022).

25. Public Health England (PHE). Multi-drug Resistant *Shigella sonnei* Cluster (CTX-M-27) Probably Associated with MSM. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/769163/SHGLL-SONN_2018275.pdf (London, 2019).

26. Charles, H. et al. Outbreak of sexually transmitted, extensively drug-resistant *Shigella sonnei* in the UK, 2021-22: a descriptive epidemiological study. *Lancet Infect. Dis.* **22**, 1503–1510 (2022).

27. Baker, K. S. et al. Genomic epidemiology of *Shigella* in the United Kingdom shows transmission of pathogen sublineages and determinants of antimicrobial resistance. *Sci. Rep.* **8**, 7389 (2018).

28. Locke, R. K., Greig, D. R., Jenkins, C., Dallman, T. J. & Cowley, L. A. Acquisition and loss of CTX-M plasmids in *Shigella* species associated with MSM transmission in the UK. *Microb. Genomics* **7**, 000644 (2021).

29. Holt, K. E. et al. Tracking the establishment of local endemic populations of an emergent enteric pathogen. *Proc. Natl Acad. Sci. USA* **110**, 17522–17527 (2013).

30. Rowlinson, E., Berzkalns, A., Thibault, C., Golden, M. & Barbee, L. O14.3 High incidence of antimicrobial use and overuse in cisgender men who have sex with men at risk of bacterial STIs. *Sexually Transmitted Infect.* **97**, A48–A49 (2021).

31. Bengtsson Rebecca, J. et al. Accessory genome dynamics and structural variation of *Shigella* from persistent infections. *mBio* **12**, e00254–00221 (2021).

32. Thanh Duy, P. et al. Commensal *Escherichia coli* are a reservoir for the transfer of XDR plasmids into epidemic fluoroquinolone-resistant *Shigella sonnei*. *Nat. Microbiol.* **5**, 256–264 (2020).

33. Allen, H. et al. Evidence for re-infection and persistent carriage of *Shigella* species in adult males reporting domestically acquired infection in England. *Clin. Microbiol. Infect.* **27**, 126.e127–126.e113 (2021).

34. Mason, L., Dallman, T., Jenkins, C. & Baker, K. Genomic epidemiology of the first London outbreak of antimicrobial resistant sexually

transmitted shigellosis. Access *Microbiology* **4**, https://doi.org/10.1099/acmi.ac2021-po0344 (2022).

35. Collins, J. P. et al. Evidence of failure of oral third-generation cephalosporin treatment for *Shigella sonnei* infection. *Open Forum Infect. Dis.* **7**, ofaa113 (2020).

36. St Cyr, S. et al. Update to CDC's treatment guidelines for gonococcal infection, 2020. *MMWR Morbidity Mortal. Wkly. Rep.* **69**, 1911–1916 (2020).

37. Fifer, H., Saunders, J., Soni, S., Sadiq, S. T. & FitzGerald, M. 2018 UK national guideline for the management of infection with *Neisseria gonorrhoeae*. *Int. J. STD AIDS* **31**, 4–15 (2020).

38. Fouéré, S. et al. Update on French recommendations for the treatment of uncomplicated *Neisseria gonorrhoeae* infections. *Int. J. STD AIDS* **32**, 1081–1083 (2021).

39. Jespers, V. et al. Diagnosis and treatment of gonorrhoea: 2019 Belgian National guideline for primary care. *Acta Clin. Belgica* **77**, 186–194 (2022).

40. Migchelsen, S. S. et al. Sexually Transmitted Infections and Screening for Chlamydia in England, 2021 (UK Health Security Agency (UKHSA), 2022).

41. Carter, L. L. et al. Global genomic surveillance strategy for pathogens with pandemic and epidemic potential 2022–2032. *Bull. World Health Organ.* **100**, 239–239a (2022).

42. Dallman, T. J. et al. Use of whole-genome sequencing for the public health surveillance of *Shigella sonnei* in England and Wales, 2015. *J. Med. Microbiol.* **65**, 882–884 (2016).

43. Mitchell, H. D. et al. Use of whole-genome sequencing to identify clusters of *Shigella flexneri* associated with sexual transmission in men who have sex with men in England: a validation study using linked behavioural data. *Microb. Genomics* **5**, e000311 (2019).

44. Dallman, T. J. et al. Emergence of novel strains of *Shigella flexneri* associated with sexual transmission in adult men in England, 2019-2020. *J. Med. Microbiol.* **70**, 001437 (2021).

45. Yara, D. A., Greig, D. R., Gally, D. L., Dallman, T. J. & Jenkins, C. Comparison of Shiga toxin-encoding bacteriophages in highly pathogenic strains of Shiga toxin-producing *Escherichia coli* O157:H7 in the UK. *Microb. Genomics* **6**, e000334 (2020).

46. von Mentzer, A. et al. Long-read-sequenced reference genomes of the seven major lineages of enterotoxigenic *Escherichia coli* (ETEC) circulating in modern time. *Sci. Rep.* **11**, 9256 (2021).

47. Wick, R. R., Judd, L. M., Gorrie, C. L. & Holt, K. E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. *PLoS Comput. Biol.* **13**, e1005595 (2017).

48. Walker, B. J. et al. Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement. *PLoS ONE* **9**, e112963 (2014).

49. Nanopore Technologies. Guppy, https://nanoporetech.com/ (Oxford, 2022).

50. Li, H. Minimap2: pairwise alignment for nucleotide sequences. *Bioinformatics* **34**, 3094–3100 (2018).

51. Seemann, T. Prokka: rapid prokaryotic genome annotation. *Bioinformatics* **30**, 2068–2069 (2014).

52. Zhou, Z., Alikhan, N. F., Mohamed, K., Fan, Y. & Achtman, M. The EnteroBase user's guide, with case studies on *Salmonella* transmissions, *Yersinia pestis* phylogeny, and *Escherichia* core genomic diversity. *Genome Res.* **30**, 138–152 (2020).

53. Zhou, Z. et al. GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens. *Genome Res.* **28**, 1395–1404 (2018).

54. Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. *Nucleic Acids Res.* **49**, W293–W296 (2021).

55. SRAToolkitDevelopmentTeam. SRA Toolkit https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=software, (GitHub Repository, 2022).
56. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* **30**, 2114–2120 (2014).

57. Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. *Bioinformatics* **32**, 3047–3048 (2016).

58. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* **25**, 1754–1760 (2009).

59. Broad Institute. Picard Toolkit v. 2.27.2 https://broadinstitute.github.io/picard/ (GitHub Repository, 2019).

60. Li, H. et al. The Sequence Alignment/Map format and SAMtools. *Bioinformatics* **25**, 2078–2079 (2009).

61. Danecek, P. et al. Twelve years of SAMtools and BCFtools. *Giga-science* **10**, giab008 (2021).

62. Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. *Bioinformatics* **27**, 2987–2993 (2011).

63. Free Software Foundation. GNU sed https://www.gnu.org/software/sed/ (2013).

64. Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. *Nucleic Acids Res.* **43**, e15 (2015).

65. Nguyen, L. T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. *Mol. Biol. Evol.* **32**, 268–274 (2015).

66. Tonkin-Hill, G., Lees, J., Bentley, S., Frost, S. & Corander, J. Rhier-BAPS: An R implementation of the population clustering algorithm hierBAPS. *Wellcome Open Res.* **3**, 93 (2018).

67. Cheng, L., Connor, T. R., Sirén, J., Aanensen, D. M. & Corander, J. Hierarchical and spatially explicit clustering of DNA sequences with BAPS software. *Mol. Biol. Evol.* **30**, 1224–1228 (2013).

68. Rambaut, A., Lam, T. T., Max Carvalho, L. & Pybus, O. G. Exploring the temporal structure of heterochronous sequences using Tem-pEst (formerly Path-O-Gen). *Virus Evol.* **2**, vew007 (2016).

69. Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. *Bioinformatics* **29**, 1072–1075 (2013).

70. Feldgarden, M. et al. AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among anti-microbial resistance, stress response, and virulence. *Sci. Rep.* **11**, 12728 (2021).

71. Holt, K., Hawkey, J., Paranagama, K. & Wright, S. Sonneityping script parses Mykrobe predict results for *Shigella sonnei* https://github.com/katholt/sonneityping (2021).

72. Hunt, M. et al. Antibiotic resistance prediction for *Mycobacterium tuberculosis* from genome sequence data with Mykrobe. *Wellcome Open Res.* **4**, 191 (2019).

73. Bradley, P. et al. Rapid antibiotic-resistance predictions from genome sequence data for *Staphylococcus aureus* and *Mycobacterium tuberculosis*. *Nat. Commun.* **6**, 10063 (2015).

74. Hunt, M. et al. Circlator: automated circularization of genome assemblies using long sequencing reads. *Genome Biol.* **16**, 294 (2015).

75. Gilchrist, C. L. M. & Chooi, Y. H. Clinker & clusterman.js: Automatic generation of gene cluster comparison figures. *Bioinformatics*, https://doi.org/10.1093/bioinformatics/btab007 (2021).

76. Bingmann, T., Bradley, P., Gauger, F. & Iqbal, Z. COBS: A Compact Bit-Sliced Signature Index Preprint at https://arxiv.org/abs/1905.09624 (2019).

77. Blackwell, G. A. et al. Exploring bacterial diversity via a curated and searchable snapshot of archived DNA sequences. *PLoS Biol.* **19**, e3001421 (2021).

78. Alikhan, N.-F., Petty, N. K., Ben Zakour, N. L. & Beatson, S. A. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. *BMC Genomics* **12**, 402 (2011).

79. García-Alcalde, F. et al. Qualimap: evaluating next-generation sequencing alignment data. *Bioinformatics* **28**, 2678–2679 (2012).

80. Sprouffske, K. & Wagner, A. Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves. *BMC Bioinform.* **17**, 172 (2016).

81. European Committee on Antimicrobial Susceptibility Testing (EUCAST). *Clinical Breakpoints and Dosing of Antibiotics—Bacteria (V 12.0)*, https://www.eucast.org/clinical_breakpoints/ (Vaxjo, 2022).

## Acknowledgements

This study was funded by the National Institute for Health and Care Research (NIHR) Health Protection Research Unit in Gastrointestinal Infections, a partnership between the UK Health Security Agency, the University of Liverpool, and the University of Warwick. The views expressed are those of the author(s) and not necessarily those of the NIHR, the UK Health Security Agency or the Department of Health and Social Care. The authors thank Hattie Webb, Kaitlin Tagg, Jason Folster, Jean Whichard, and other staff at the United States Centres for Disease Control for helpful discussions to guide access to genome data from PulseNet. The authors are also grateful to pathology laboratories that provided isolates for genomic surveillance. This work was also supported by MRC and BBSRC grants held by KSB (MR/R020787/1 and BB/V009184/1, respectively).

## Author contributions

L.C.E.M., C.J., and K.B. were responsible for the writing of the original draft. All authors contributed to the review and editing of the manuscript. L.C.E.M., D.G., L.C., S.P., E.M., G.B., C.C., M.D.S., J.D., E.S., E.N., C.J., and K.S.B. were responsible for formal analysis. L.C.E.M., D.G., L.C., S.P., E.M., G.B., C.C., M.D.S., J.D., A.G., I.S., E.S., V.S., J.I., S.L., C.J., and K.S.B. were responsible for the investigation. L.C.E.M., D.G., S.P., E.M., G.B., J.D., A.G., E.S., C.J., K.S.B. were responsible for the methodology. L.C.E.M., F.X.W., P.J.C., C.J., K.S.B. were responsible for conceptualisation. D.G., E.M., J.D., A.G., D.I., S.L., E.N., and C.J. were responsible for data curation. D.G., I.S., V.S., J.I., B.P.H., S.L., E.N., F.X.W., C.J., and K.S.B. were responsible for providing resources. M.D.S. was responsible for validation. V.S., J.I., B.P.H., F.X.W., C.J., K.S.B. were responsible for funding acquisition. V.S., J.I., C.J. and K.S.B. were responsible for project administration. R.J.B., V.S., J.I., C.J., and K.S.B. were responsible for supervision. L.C.E.M., D.G., S.P., G.B., C.C., and M.D.S. were responsible for visualisation.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-37672-w.

**Correspondence** and requests for materials should be addressed to Kate S. Baker.

**Peer review information** Nature Communications thanks Jay Worley and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
56. Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* **30**, 2114–2120 (2014).

57. Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. *Bioinformatics* **32**, 3047–3048 (2016).

58. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* **25**, 1754–1760 (2009).

59. Broad Institute. Picard Toolkit v. 2.27.2 https://broadinstitute.github.io/picard/ (GitHub Repository, 2019).

60. Li, H. et al. The Sequence Alignment/Map format and SAMtools. *Bioinformatics* **25**, 2078–2079 (2009).

61. Danecek, P. et al. Twelve years of SAMtools and BCFtools. *Giga-science* **10**, giab008 (2021).

62. Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. *Bioinformatics* **27**, 2987–2993 (2011).

63. Free Software Foundation. GNU sed https://www.gnu.org/software/sed/ (2013).

64. Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. *Nucleic Acids Res.* **43**, e15 (2015).

65. Nguyen, L. T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. *Mol. Biol. Evol.* **32**, 268–274 (2015).

66. Tonkin-Hill, G., Lees, J., Bentley, S., Frost, S. & Corander, J. Rhier-BAPS: An R implementation of the population clustering algorithm hierBAPS. *Wellcome Open Res.* **3**, 93 (2018).

67. Cheng, L., Connor, T. R., Sirén, J., Aanensen, D. M. & Corander, J. Hierarchical and spatially explicit clustering of DNA sequences with BAPS software. *Mol. Biol. Evol.* **30**, 1224–1228 (2013).

68. Rambaut, A., Lam, T. T., Max Carvalho, L. & Pybus, O. G. Exploring the temporal structure of heterochronous sequences using Tem-pEst (formerly Path-O-Gen). *Virus Evol.* **2**, vew007 (2016).

69. Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. *Bioinformatics* **29**, 1072–1075 (2013).

70. Feldgarden, M. et al. AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among anti-microbial resistance, stress response, and virulence. *Sci. Rep.* **11**, 12728 (2021).

71. Holt, K., Hawkey, J., Paranagama, K. & Wright, S. Sonneityping script parses Mykrobe predict results for *Shigella sonnei* https://github.com/katholt/sonneityping (2021).

72. Hunt, M. et al. Antibiotic resistance prediction for *Mycobacterium tuberculosis* from genome sequence data with Mykrobe. *Wellcome Open Res.* **4**, 191 (2019).

73. Bradley, P. et al. Rapid antibiotic-resistance predictions from genome sequence data for *Staphylococcus aureus* and *Mycobacterium tuberculosis*. *Nat. Commun.* **6**, 10063 (2015).

74. Hunt, M. et al. Circlator: automated circularization of genome assemblies using long sequencing reads. *Genome Biol.* **16**, 294 (2015).

75. Gilchrist, C. L. M. & Chooi, Y. H. Clinker & clusterman.js: Automatic generation of gene cluster comparison figures. *Bioinformatics*, https://doi.org/10.1093/bioinformatics/btab007 (2021).

76. Bingmann, T., Bradley, P., Gauger, F. & Iqbal, Z. COBS: A Compact Bit-Sliced Signature Index Preprint at https://arxiv.org/abs/1905.09624 (2019).

77. Blackwell, G. A. et al. Exploring bacterial diversity via a curated and searchable snapshot of archived DNA sequences. *PLoS Biol.* **19**, e3001421 (2021).

78. Alikhan, N.-F., Petty, N. K., Ben Zakour, N. L. & Beatson, S. A. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. *BMC Genomics* **12**, 402 (2011).

79. García-Alcalde, F. et al. Qualmap: evaluating next-generation sequencing alignment data. *Bioinformatics* **28**, 2678–2679 (2012).

80. Sprouffske, K. & Wagner, A. Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves. *BMC Bioinform.* **17**, 172 (2016).

81. European Committee on Antimicrobial Susceptibility Testing (EUCAST). *Clinical Breakpoints and Dosing of Antibiotics—Bacteria (V 12.0)*, https://www.eucast.org/clinical_breakpoints/ (Vaxjo, 2022).

## Acknowledgements

This study was funded by the National Institute for Health and Care Research (NIHR) Health Protection Research Unit in Gastrointestinal Infections, a partnership between the UK Health Security Agency, the University of Liverpool, and the University of Warwick. The views expressed are those of the author(s) and not necessarily those of the NIHR, the UK Health Security Agency or the Department of Health and Social Care. The authors thank Hattie Webb, Kaitlin Tagg, Jason Folster, Jean Whichard, and other staff at the United States Centres for Disease Control for helpful discussions to guide access to genome data from PulseNet. The authors are also grateful to pathology laboratories that provided isolates for genomic surveillance. This work was also supported by MRC and BBSRC grants held by KSB (MR/R020787/1 and BB/V009184/1, respectively).

## Author contributions

L.C.E.M., C.J., and K.B. were responsible for the writing of the original draft. All authors contributed to the review and editing of the manuscript. L.C.E.M., D.G., L.C., S.P., E.M., G.B., C.C., M.D.S., J.D., E.S., E.N., C.J., and K.S.B. were responsible for formal analysis. L.C.E.M., D.G., L.C., S.P., E.M., G.B., C.C., M.D.S., J.D., A.G., I.S., E.S., V.S., J.I., S.L., C.J., and K.S.B. were responsible for the investigation. L.C.E.M., D.G., S.P., E.M., G.B., J.D., A.G., E.S., C.J., K.S.B. were responsible for the methodology. L.C.E.M., F.X.W., P.J.C., C.J., K.S.B. were responsible for conceptualisation. D.G., E.M., J.D., A.G., D.I., S.L., E.N., and C.J. were responsible for data curation. D.G., I.S., V.S., J.I., B.P.H., S.L., E.N., F.X.W., C.J., and K.S.B. were responsible for providing resources. M.D.S. was responsible for validation. V.S., J.I., B.P.H., F.X.W., C.J., K.S.B. were responsible for funding acquisition. V.S., J.I., C.J. and K.S.B. were responsible for project administration. R.J.B., V.S., J.I., C.J., and K.S.B. were responsible for supervision. L.C.E.M., D.G., S.P., G.B., C.C., and M.D.S. were responsible for visualisation.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-37672-w.

**Correspondence** and requests for materials should be addressed to Kate S. Baker.

**Peer review information** *Nature Communications* thanks Jay Worley and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023, corrected publication 2023
**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023, corrected publication 2023