# Genomic epidemiology of Shigella in the United Kingdom shows transmission of pathogen sublineages and determinants of antimicrobial resistance

Kate S Baker, Timothy Dallman, Nigel Field, Tristan Childs, Holly Mitchell, Martin Day, François-Xavier Weill, Sophie Lefevre, Mathieu Tourdjman, Gwenda Hughes, et al.

► To cite this version:

Kate S Baker, Timothy Dallman, Nigel Field, Tristan Childs, Holly Mitchell, et al.. Genomic epidemiology of Shigella in the United Kingdom shows transmission of pathogen sublineages and determinants of antimicrobial resistance. Scientific Reports, 2018, 8 (1), pp.7389. 10.1038/s41598-018-25764-3. pasteur-02006006

HAL Id: pasteur-02006006

https://pasteur.hal.science/pasteur-02006006

Submitted on 4 Feb 2019

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>
SCIENTIFIC REPORTS

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

OPEN

# Genomic epidemiology of *Shigella* in the United Kingdom shows transmission of pathogen sublineages and determinants of antimicrobial resistance

Received: 7 December 2017

Accepted: 27 April 2018

Published online: 09 May 2018

Kate S. Baker<sup>D 1</sup>, Timothy J. Dallman<sup>2</sup>, Nigel Field<sup>3</sup>, Tristan Childs<sup>4</sup>, Holly Mitchell<sup>3</sup>, Martin Day<sup>2</sup>, François-Xavier Weill<sup>5</sup>, Sophie Lefèvre<sup>5</sup>, Mathieu Tourdjman<sup>6</sup>, Gwenda Hughes<sup>7</sup>, Claire Jenkins<sup>2</sup> & Nicholas Thomson<sup>8,9</sup>

Shigella are globally important diarrhoeal pathogens that are endemic in low-to-middle income nations and also occur in high income nations, typically in travellers or community-based risk-groups. Shigella phylogenetics reveals population structures that are more reliable than those built with traditional typing methods, and has identified sublineages associated with specific geographical regions or patient groups. Genomic analyses reveal temporal increases in *Shigella* antimicrobial resistance (AMR) gene content, which is frequently encoded on mobile genetic elements. Here, we whole genome sequenced representative subsamples of *S. flexneri* 2a and *S. sonnei* (n = 366) from the United Kingdom from 2008 to 2014, and analysed these alongside publicly available data to make qualitative insights on the genomic epidemiology of shigellosis and its AMR within the broader global context. Combined phylogenetic, epidemiological and genomic anlayses revealed the presence of domestically-circulating sublineages in patient risk-groups and the importation of travel-related sublineages from both Africa and Asia, including ciprofloxacin-resistant sublineages of both species from Asia. Genomic analyses revealed common AMR determinants among travel-related and domestically-acquired isolates, and the evolution of mutations associated with reduced quinolone susceptibility in domestically-circulating sublineages. Collectively, this study provides unprecedented insights on the contribution and mobility of endemic and travel-imported sublineages and AMR determinants responsible for disease in a high-income nation.

*Shigella* are the most common bacterial cause of moderate to severe diarrhoea in children under 5 years old in Asia and Africa<fnref n="1" />. In high income nations, shigellosis is often related to travel but can also circulate as a sexually transmissible illness of men who have sex with men (MSM)²,<fnref n="3" />. Among the four species comprising the genus *Shigella*, *S. sonnei* and *S. flexneri* cause the greatest disease burden. Among the diverse (>50) serotypes of *S. flexneri*, the 2a serotype is historically the most prevalent worldwide and in the United Kingdom (UK)⁴,⁵. However, increased resolution of *Shigella* species population structures has been provided by Whole Genome Sequence

<footnote n="1">¹Institute for Integrative Biology, University of Liverpool, Liverpool, L69 7ZB, United Kingdom. ²Gastrointestinal Bacterial Reference Unit, National Infection Service, Public Health England, London, NW9 5HT, United Kingdom.</footnote>

<footnote n="3">³Centre for Molecular Epidemiology and Translational Research, Institute for Global Health, UCL, London, United Kingdom. ⁴Centre for Infectious Disease Surveillance and Control, National Infection Service, Public Health England, London, NW9 5HT, United Kingdom. ⁵Institut Pasteur, Unité des Bactéries Pathogènes Entériques, Paris, 75015, France. ⁶Santé Publique France, the French Public Health Agency, Saint-Maurice, France. ⁷Department of HIV and STIs, National Infection Service, Public Health England, London, NW9 5HT, United Kingdom. ⁸Wellcome Trust Sanger Institute, Hinxton, CB10 1SA, United Kingdom. ⁹London School of Hygiene and Tropical Medicine, London, WC1E 7HT, United Kingdom. Correspondence and requests for materials should be addressed to K.S.B. (email: kbaker@liverpool.ac.uk) or N.T. (email: nrt@sanger.ac.uk)</footnote>
Analysis (WGSA). *S. sonnei* is split into five WGSA subtypes called ‘lineages’, with the original four lineages⁶ being joined recently by a fifth reported in Latin America and Africa⁷. *S. flexneri* has been similarly divided into seven WGSA subtypes called ‘phylogroups’ (PGs)⁸. In addition to these species-defining studies, prevailing and emerging lower-order WGSA subtypes (usually called ‘sublineages’) of *S. sonnei* and *S. flexneri* have also been described in countries in Asia, the Middle East and Latin America⁷,⁹⁻¹². WGSA has also been used to demonstrate the international movement of *Shigella* sublineages where WGSA subtypes prevailing in one geographical region are found in another, isolated from patients with an epidemiological link (e.g. travel)¹³,¹⁴. For example, a specific sublineage of *S. sonnei* that dominates in Israel circulates globally among Orthodox Jewish communities¹⁰.

On top of this global population structure of shigellae WGSA subtypes, there is a horizontal and vertical transmission network of genes, including those encoding antimicrobial resistance (AMR). When WGSA subtyping is combined with epidemiological metadata and genome analyses (genomic epidemiology), the importance of AMR in shaping epidemic emergences and broader population structure in shigellae has been revealed. Studies over long time scales have shown the accumulation of AMR in *Shigella* over time⁸,¹⁵,¹⁶, and both local and global epidemic emergences of WGSA subtypes have been found to be associated with resistance to therapeutically-important antimicrobial classes⁶,⁹,¹²,¹³. AMR determinants in *Shigella* include Mobile Genetic Elements (MGEs), such as small and large resistance plasmids, and chromosomally-integrated islands, some of which have been found in multiple *Shigella* species. For example, the *Shigella* Resistance Locus Multiple Drug Resistance Element (SRL-MDRE) encoding resistance to ampicillin, chloramphenicol and tetracycline has been found in widely geographically-distributed *S. dysenteriae*, *S. flexneri* and *S. sonnei*⁷,⁸,¹⁵,¹⁷. In addition to these horizontally-acquired AMR determinants, chromosomal mutations in the Quinolone Resistance Determining Region (QRDR) (that confer resistance to ciprofloxacin, the recommended treatment for shigellosis)¹⁸ have also been reported in shigellae isolated from Asia, Latin America, and high-income nations⁷,¹¹,¹⁴,¹⁹. Ciprofloxacin-resistant *Shigella* were recently recognised by the World Health Organisation as one of the top dozen agents for which new antimicrobials are urgently needed. Common AMR determinants being found in geographically-separated and genetically-distinct *Shigella* *subtypes* point to the independent mobility of AMR determinants as well as convergent evolution.

With an increased resolution and volume of information on the population structure of *Shigella* and its AMR determinants available, these data can be used to enhance surveillance in public health laboratories that are using whole genome sequencing, such as the Gastrointestinal Bacteria Reference Unit (GBRU) at Public Health England (PHE) in the UK²⁰. In the UK, shigellosis has historically been associated with travel and only low-levels of domestic acquisition (i.e. illness in patients without a recent history of travel to a shigellosis-endemic area). Recently however, epidemics of domestically-acquired shigellosis (*S. sonnei*, *S. flexneri* 2a and *S. flexneri* 3a) have been linked with transmission among adult males²¹, and more specifically among men who have sex with men (MSM)⁵,²². As part of a recent study investigating AMR emergence among these parallel epidemics in MSM²³, we whole genome sequenced representative subsamples of *S. sonnei* (n = 187) and *S. flexneri* 2a (n = 179) from the UK sampled between 2008 and 2014. Here, in this allied study, we integrate this new cross sectional data from the UK with existing reference datasets of *Shigella* WGSA subtypes and AMR determinants from around the world to identify those responsible for shigellosis in the UK. By providing the context necessary for interpreting national surveillance patterns, we reveal important qualitative insights on the importation and domestic transmission of *Shigella* in the UK, as well as on the prevalence and mobility of AMR among these globally-important pathogens.

## Methods

**Nationally representative isolates.** To investigate an epidemic phase of *Shigella* transmission associated with MSM²¹,²³, isolates of *S. sonnei* and *S. flexneri* 2a submitted to the GBRU at PHE were representatively subsampled, with a random sample of isolates from each year proportional to the total number of annual cases (for each species) being taken from across 2008–2014 (Fig. 1). In doing so, we aimed to balance an economically-practical number of isolates for WGSA, with the need to have sufficient statistical power to detect emerging sublineages of domestic transmission (this calculation was based on data in¹³,²¹). For *S. flexneri* 2a, there were no selection criteria based on age, whereas for *S. sonnei* (owing to the larger number of cases that occur annually) selection was restricted to patients aged 16–60 years old. To focus on the epidemic domestic transmission, we sequenced a greater proportion of isolates from patients who did not report recent high-risk travel for both species (Fig. 1). The final proportions were as follows: for *S. sonnei* 3% (45 of 1539 cases) of travel-associated cases and 8% (142 of 1751 cases) of non-travel associated cases from patients aged 16–60 years old, and for *S. flexneri* 2a 9% (32 of 368 cases) of travel-associated cases and 18% (147 of 820 cases) of non-travel associated cases from patients of all ages. This resulted in a total of 187 *S. sonnei* and 179 *S. flexneri* 2a isolates (Fig. 1).

**Publically available data.** Various *S. sonnei* reference isolates (n = 40) from other WGSA studies were included to provide interpretative context. Representative contextual isolates and the associated studies were as follows: the original four global lineages of *S. sonnei*, including those classified into the multidrug resistant and globally-disseminated ‘Global III’ sublineage⁶; the newly identified global Lineage V and a regionally-epidemic Latin America sublineage, Latin America IIIa⁷; a sublineage associated with transmission in Israel and Orthodox Jewish communities elsewhere¹⁰; and a ciprofloxacin-resistant sublineage circulating in, and disseminating from, Asia, the ‘Central Asia III’ lineage¹⁴ (Supplementary Data). For *S. flexneri* 2a, isolates sequenced in a population-structure-defining study⁸ were included in analyses to anchor the positions of the *S. flexneri* 2a representative subsample, and investigate geographical associations. Specifically, this included representatives of all PGs (n = 9) and all PG3 *S. flexneri* 2a isolates (n = 69) (Supplementary Data).

**Bioinformatic analyses.** DNA was extracted and isolates whole genome sequenced as previously described¹³. Whole genome sequence phylogenetic trees were constructed through mapping of sequence data
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Representative sampling of isolates for whole genome sequencing. The figure shows a graph of the number of isolates sequenced per year by travel and species (upper graph), overlaid with the number of domestically-acquired (non-travel) and travel-associated cases of shigellosis (for all ages). Non-travel isolates were overrepresented (see Methods).</figcaption></figure>

(using smalt, http://www.sanger.ac.uk/science/tools/smalt-0) to reference genomes (*S. flexneri* 2a 2457 T, *S. sonnei* 53 G), followed by removal of mobile (plasmids, insertion sequences, other described mobile islands<sup>6,24)</sup> and repetitive elements (through removal of known repeat regions and repeat mapping discard), and sites of recombination (using gubbins<sup>25)</sup> followed by extraction of variable sites (using snpSites<sup>26)</sup>). This resulted in multiple sequence alignments of 19348, 11141, and 9790 variant sites in length for *S. flexneri* (all PGs), *S. flexneri* (PG3 only), and *S. sonnei* respectively. Multiple sequence alignments were used to infer maximum likelihood phylogenetic trees using RAxML with 100 bootstrap replicates<sup>27</sup>. Acquired AMR genes were identified using ResFinder<sup>28</sup> on draft genome assemblies (generated by<sup>29</sup>) and QRDR mutations were detected as previously described<sup>7</sup>. All sequence data are publicly available at the European Nucleotide Archive under project number PRJEB12097 (individual accessions in Supplementary Data).

**Ethical statement.** No individual patient consent was required or sought as PHE 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). The project was reviewed and approved by the PHE Research Support and Governance Office and Caldicott Panel deemed to comply with public health surveillance standards.

**Data availability statement.** All sequencing and genome assembly data for isolates in this study are available under project number PRJEB12097 in the European Nucleotide Archive. Further detail linking isolates and sample accession numbers within the project are available in Supplementary Data.

## Results

**Mobility of S. sonnei and S. flexneri 2a WGSA subtypes.** For both species, the whole genome phylogenies revealed a dominant subtype. Most *S. flexneri* 2a isolates (176 of 179) belonged to PG3 (Fig. 2A) and most (181 of 187) of the *S. sonnei* isolates belonged to Lineage III (Fig. 3). For *S. flexneri* 2a, this likely represents the strong association of the 2a serotype with the PG3 WGSA subtype. The association was not exclusive though as one PG2 and two PG1 isolates were among the *S. flexneri* 2a samples, reaffirming that serotype switching between PGs occurs<sup>8</sup>. For *S. sonnei*, the dominance of Lineage III reflects a greater burden of disease caused by this lineage, as previously shown<sup>6</sup>. Despite the Lineage III dominance however, further diversity of *S. sonnei* was detected in the UK samples. Specifically, we observed two Lineage V isolates from domestically-acquired cases (despite Lineage V only having been previously reported in Latin America or Africa)<sup>7</sup>, one Lineage I isolate from a patient recently-returned from Asia and three Lineage II isolates from patients without a history of recent travel (Fig. 3). In fact, the only lineage of *S. sonnei* not detected, Lineage IV, has yet to be seen outside of its original description<sup>6</sup>.

In addition to this diversity at the highest level of nomenclature (i.e. lineage/PG), analysis at the sub-lineage level enabled further examination of the role of travel in the epidemiology of shigellosis in the UK. Travel-associated isolates were often phylogenetically proximate to reference isolates from that same region, differing by as few as 3 SNPs (Figs 2B and 3). For *S. flexneri* PG3, there was a sublineage related to travel from Asia,
<figure><img src="image_4.png" /><figcaption>**Figure 2.** *Shigella flexneri* 2a UK cross section in context. (A) The smaller tree shows the broad context of phylogroups (PG) 1–7 of *S. flexneri* and intercontinental MSM-associated (MSMA) *S. flexneri* 3a, with the relationships of three non-PG3 isolates from this study. PG3 is overlaid by a grey polygon in the smaller tree and is shown in full in the larger tree. (B) The larger tree is a mid-point rooted maximum likelihood tree and shows the detailed evolutionary relationships of UK cross section isolates (black branches) with reference isolates (coloured branches). Reference isolate branches are coloured by region of sample origin, and region of recent travel is overlaid for UK samples as coloured circles on tree tips (both coloured according to the inlaid region key). Two MSM-associated sublineages are highlighted by grey boxes, and the *S. flexneri* 2a Central Asia lineage is indicated at the defining internal node. AMR features are shown in adjacent tracks. The leftmost track shows the presence of QRDR mutations and one 'Other', which designates a quadruple mutant also carrying the *qnrS1* gene. Subsequent gray scale tracks indicate the presence of antimicrobial resistance genes (ARGs, labelled above) associated with known *Shigella* MGEs. Specifically (from left to right), the Tn7/Int2, SRL-MDRE, Small R-plasmids, Minor clade plasmid, pKSR100 (Table 1).</figcaption></figure>

containing Asian reference isolates as well as most (11 of the 17) of the patients recently returned from Asia (specific locations were India, Pakistan, Sri Lanka and Afghanistan). To reflect these origins, and mirror *S. sonnei* nomenclature, this sublineage was designated *S. flexneri* 2a Central Asia (Fig. 2B). The cognate *S. sonnei* Central Asia III sublineage (recently described as disseminating from Asia¹⁴) was similarly detectable (Fig. 3). Specifically, 12 of the 15 isolates of *S. sonnei* from patients who had recently travelled to Asia (China and India) belonged to this sublineage. In addition to these Asian-travel related sublineages, there was a similar *S. sonnei* sublineage related to travel from Africa (see African travel related, Fig. 3). Specifically, the majority (13 of 19) of isolates from patients who had recently travelled to Africa belonged to a sublineage within the Global III sublineage (Fig. 3). Twelve of these 13 patients had recently returned from Egypt, and one from Nigeria, so it is likely that this sublineage represents the repeated importation of a sublineage that is nationally- or regionally-dominant in Northern Africa. In addition to these Asian- and African-related travel lineages, an *S. sonnei* isolate belonging to the Latin America IIIa sublineage was detected in a patient who had recently returned from a Latin American country where this lineage is highly prevalent⁷.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** *Shigella sonnei* UK cross section in context. The midpoint-rooted maximum likelihood phylogenetic tree shows isolates from this study in the context of five global lineages of *S. sonnei*. Lineage III has further relevant sublineages indicated at the internal nodes and MSM-associated sublineages are highlighted in grey. Isolates from patients recently returned from endemic regions are indicated by circles overlying the tree tips which are coloured according to the inlaid key. Reference isolates from previously described WGSA subtypes (lineages/sublineages) are indicated in the leftmost adjacent track, and tree branches for references isolates are also coloured according to the inlaid key. AMR features are shown in subsequent tracks. The second leftmost track shows the presence of QRDR mutations and two isolates containing ‘Other’ quinolone resistance determinants (including an isolate with a single QRDR mutation and a *qepA* gene, and an isolate containing a *qnrB19* gene). Subsequent gray scale tracks indicate the presence of antimicrobial resistance genes (ARGs, labelled above) associated with known *Shigella* MGEs. Specifically (from left to right), the Tn7/Int2, SpA, MSM-associated sublineage 3 plasmid, pKSR100, SRL-MDRE (Table 1).</figcaption></figure>

Sublineages transmitting domestically in patient communities were also detectable among the isolates. A separate analysis of this data incorporating patient age and gender identified the six sublineages associated with transmission in MSM (MSM-associated, Figs 2 and 3)²³. For S. flexneri 2a there was one large (47 isolates), low-diversity (maximum pair-wise distance 20 SNPs) MSM-associated sublineage, and one smaller (7 isolates), higher-diversity (maximum pair-wise distance 148 SNPs) MSM-associated sublineage (Fig. 2). Among the S. sonnei isolates, there were four MSM-associated sublineages with maximum pairwise distances of between 15 and 30 SNPs (Fig. 3). One MSM-associated sublineage of S. sonnei was nested within the Central Asia III sublineage, and three other MSM-associated S. sonnei sublineages belonged to the Global III sublineage (Fig. 3), which is thought to have disseminated globally after the acquisition of the Tn7/Int2 AMR determinant⁶. In addition to the MSM-associated sublineages, isolates belonging to the Orthodox Jewish community-associated sublineage (also within Global III sublineage) were detected from years spanning 2008–2014 (Fig. 3). This is consistent with a previous study that demonstrates that this Orthodox Jewish community-associated sublineage is circulating between Israel and the UK, where it causes domestically-transmitted outbreaks¹⁰,³⁰.
<table><caption>**Table 1.** Antimicrobial resistance determinants of *Shigella* and species distributions. The *strA* gene has been interrupted by the introduction of *dfra14* on this plasmid meaning no Streptomycin resistance would be conferred.</caption><thead><tr><th>Type</th><th>Name</th><th>Resistance genes</th><th>Antimicrobial classes</th><th>Species reports</th></tr></thead><tbody><tr><td rowspan="2">Chromosomal Islands</td><td>SRL-MDRE</td><td>tetB<br>bla<sub>OXA-1</sub><br>catA1</td><td>Tetracyclines<br>Ampicillin<br>Chloramphenicol</td><td>S. dysenteriae<sup>15</sup><br>S. flexneri<sup>8,17</sup><br>S. sonnei<sup>7</sup></td></tr><tr><td>Tn7/Int2</td><td>aadA1<br>dfrA1</td><td>Aminoglycosides<br>Trimethoprim</td><td>S. flexneri 2a (this study)<br>S. flexneri Xv<sup>10,33</sup><br>S. sonnei<sup>6,7</sup></td></tr><tr><td rowspan="3">Small resistance plasmids</td><td>pSS046_spA</td><td>sul2<br>strA/strB<br>tetA</td><td>Sulphonamides<br>Streptomycin<br>Tetracyclines</td><td>S. sonnei<sup>6,32</sup></td></tr><tr><td>pABC-3</td><td>dfrA14<br>sul2<br>strB<sup>a</sup></td><td>Trimethoprim<br>Sulphonamides</td><td>S. sonnei<sup>7,31</sup><br>S. flexneri 3a<sup>13</sup><br>S. flexneri 2a (this study)</td></tr><tr><td>pSFxv_3</td><td>strA/strB<br>sul2</td><td>Streptomycin<br>Sulphonamides</td><td>S. flexneri Xv<sup>12</sup><br>S. flexneri 2a (this study)</td></tr><tr><td rowspan="3">Large azithromycin resistance plasmids</td><td>pKSR100</td><td>mphA<br>ermB<br>bla<sub>TEM-1B</sub><br>aadA5<br>dfrA17<br>sul1</td><td>Macrolides<br>Macrolides<br>Ampicillin<br>Aminoglycosides<br>Trimethoprim<br>Sulphonamides</td><td>MSM-associated:<br>S. flexneri 3a<sup>23</sup><br>S. flexneri 2a<sup>23</sup><br>S. sonnei<sup>23</sup></td></tr><tr><td>Minor MSMA<br>sublineage plasmid</td><td>mphA<br>tetA<br>bla<sub>TEM-1A</sub></td><td>Macrolides<br>Tetracyclines<br>Ampicillin</td><td>S. flexneri 2a<sup>23</sup></td></tr><tr><td>MSMA sublineage<br>3 plasmid</td><td>mphA<br>dfrA1<br>sul1</td><td>Macrolides<br>Trimethoprim<br>Sulphonamides</td><td>S. sonnei<sup>23</sup></td></tr></tbody></table>

**Shared AMR determinants among Shigella species.** There were various horizontally-transmissible AMR genes among the isolates, often carried on known MGEs (Table 1). These Antimicrobial Resistance Genes (ARGs) encoded resistance against virtually all known antimicrobial classes. However, genes encoding extended spectrum beta lactamases and plasmid-encoded quinolone resistance genes were rare (10 and 3 of 366 isolates respectively, Supplementary Data). There was a chromosomally-encoded MGE encoding multiple-drug resistance in both *S. sonnei* and *S. flexneri* (i.e. the Tn7/Int2 determinant and the SRL-MDRE respectively, Table 1, Figs 2 and 3). Notably, the SRL-MDRE was also present in the only *S. sonnei* isolate belonging to the Latin America IIIa sublineage (Fig. 3), congruous with the known AMR profile of this sublineage in its endemic region<sup>7</sup>. This isolate also contained the small resistance plasmid pABC-3, typical of Latin America IIIa *S. sonnei*<sup>7,31</sup>, and also found in various MSM-associated *Shigella* sublineages (specifically the intercontinentally-distributed MSM-associated *S. flexneri* 3a<sup>13</sup> and the larger sublineage of MSM-associated *S. flexneri* 2a (Fig. 2)). Other closely-related small resistance plasmids, specifically the *S. sonnei* SpA plasmid, reported in the first sequenced genome<sup>32</sup>, and the resistance plasmid pSFxv_3, originally described in the Chinese epidemic *S. flexneri* Xv<sup>33</sup>, were also near-ubiquitous among the isolates (Figs 2 and 3, Table 1). Large resistance plasmids encoding azithromycin, and other, antimicrobial resistances were also found in isolates from MSM-associated sublineages (described more fully elsewhere)<sup>23</sup>.

Owing to its importance for treatment, ciprofloxacin resistance was also examined among the isolates. Remarkably, for both *Shigella* species, triple QRDR mutations conferring resistance to ciprofloxacin were specifically-associated with the Central Asia sublineages (*S. flexneri* 2a Central Asia in Fig. 2, Central Asia III in Fig. 3). In addition to these triple-mutant sublineages that evolved in Asia and spread to the UK, isolates and sublineages with single QRDR mutations such as gyrA S83L were present in other parts of the phylogeny (Figs 2 and 3), including in one of the domestically-transmitting MSM-associated *S. sonnei* sublineages which contained a gyrA D87Y mutation (Fig. 3). These single mutations have previously shown to be step-wise toward the evolution of fully resistant mutants<sup>14</sup>, and might indicate the potential for ciprofloxacin resistance to evolve in sublineages transmitting locally in the UK.

**Discussion**

This study reveals important information about the epidemiology of shigellosis in a high-income nation where shigellosis is frequently related to travel. Indeed, this study confirms that importation of *Shigella* sublineages through travel is a driver of shigellosis in the UK. The most striking examples of this are the repeated importation of sublineages from both Africa and Asia for *S. sonnei*, and a similar Asian-travel related sublineage of *S. flexneri* 2a. These sublineages likely represent epidemiologically prevailing sublineages circulating in endemic regions that disseminate through travel, as previously shown for the Central Asia III sublineage of *S. sonnei*<sup>14</sup> and for African and Asian-travel related sublineages of *S. flexneri* 3a<sup>13</sup>. The importation of sublineages is concerning, particularly when they are resistant to key antimicrobial classes such as the ciprofloxacin-resistant Central Asia sublineages. As our study focused on domestically-acquired shigellosis, we were uniquely able to follow the fate of some of these imported sublineages. Isolates from domestically-acquired cases of shigellosis were present in the travel-related sublineages and likely represent onward transmission from travel-related cases (or alternatively could result from misreporting of travel histories). In our study, domestically-acquired cases of
travel-related sublineages often appeared to be dead-end transmissions (i.e. did not form larger clusters of isolates from domestically-acquired infections). This is consistent with recent comprehensive WGSA surveillance of S. sonnei in the UK where travel-related cases were typically in genetic clusters of only small numbers of isolates²⁰. In this study however, we did observe one antimicrobial resistant sublineage that was imported and then went on to contribute to a domestically-transmitting epidemic. Specifically, an MSM-associated sublineage of S. sonnei was nested within the Central Asia III sublineage, indicating that a newly imported AMR strain can cause domestically-propagated epidemics should they enter an appropriate transmission network)¹³,²²,²³.

Unfortunately, biases in our subsample (proportional sampling at different depths for travel and non-travel related isolates, as well as age-restrictions for S. sonnei) make it inappropriate for us to draw quantitative inferences on the relative contribution of this importation to the shigellosis burden in the UK. However, recent results from routine whole genome sequencing surveillance of S. sonnei in the UK for the period August 2015–January 2016 show that approximately 60% of genetic clusters are likely related to travel²⁰.

In addition to directly causing domestically-transmitting epidemics, imported strains can bring in new AMR determinants that might be transferred into locally-circulating sublineages. Encouragingly, there was little evidence for the mobility of mutated QRDRs from imported lineages through chromosomal recombination. However, this was not the case for ARGs carried on MGEs where there was evidence for the exchange of AMR determinants within the UK. For example, the small resistance plasmid carried in the imported Latin America IIIa S. sonnei (pABC-3) was also found in the domestically-transmitting MSM-associated sublineages of S. flexneri 2a (Fig. 2) and S. flexneri 3a¹³ (Table 1). This indicates the potential for imported AMR determinants to horizontally transmit to locally-circulating sublineages and contribute to epidemic emergences. The imported Latin American IIIa S. sonnei isolate also carried the SRL-MDRE, which is rare in S. sonnei outside of Latin America, but could mobilise to other pathogens capable of genetic exchange with S. sonnei. The genetic exchange of AMR determinants between imported and locally-circulating strains means that the public health impact of the importation of AMR pathogens may not be the disease caused by that specific pathogen, but rather the mobilisation of the AMR determinants into other already locally-adapted pathogens and/or patient microbiota, which would then act as a reservoir for AMR emergence. In this way, imported AMR pathogens enrich the UK AMR pool, and may only be noticed only when antimicrobial usage is sufficient to drive their spread to reportable pathogens. Indeed, the ability for horizontal transfer of an AMR plasmid to drive new population-level epidemics is shown in our allied study focused on the MSM-associated sublineages²³.

The mobility of AMR determinants among *Shigella* sp. is more broadly evidenced by the detection of the same AMR determinants among the diverse global data available. Analysis of previous reports and data from this study show that different *Shigella* species across disparate geographical locations share common AMR MGEs (Table 1). That this mere handful of AMR determinants dominate the genetic landscape of *Shigella* is remarkable considering the breadth of what shigellae (globally ubiquitous, genetically-diverse pathogens capable of lateral gene transfer) would have the opportunity to acquire in the global enteric milieu. The presence of epidemiologically-dominant sublineages with stable AMR determinant associations is commonly interpreted to result from epidemiological success of the sublineage consequent to the conferred AMR phenotype(s). However, in addition to that effect, the small number of recurrent AMR determinants among enormously diverse shigellae might suggest that the success of such AMR determinant-sublineage combinations may also be attributable to some measure of fitness-matching of the AMR determinant for that sublineage/species. It is picking apart these finer details of what limits and drives AMR flow among pathogenic organisms such as *Shigella* that will allow us to understand and perhaps eventually predict new emergences to prevent their occurrence.

## References

1. Liu, J. et al. Use of quantitative molecular diagnostic methods to identify causes of diarrhoea in children: a reanalysis of the GEMS case-control study. *Lancet* **388**, 1291–1301, https://doi.org/10.1016/S0140-6736(16)31529-X (2016).

2. Kotloff, K. L. et al. Global burden of Shigella infections: implications for vaccine development and implementation of control strategies. *Bulletin of the World Health Organization* **77**, 651–666 (1999).

3. Dritz, S. K. & Back, A. F. Letter: Shigella enteritis venereally transmitted. *The New England journal of medicine* **291**, 1194 (1974).

4. Livio, S. et al. Shigella isolates from the global enteric multicenter study inform vaccine development. *Clinical infectious diseases: an official publication of the Infectious Diseases Society of America* **59**, 933–941, https://doi.org/10.1093/cid/ciu468 (2014).

5. Borg, M. L. et al. Ongoing outbreak of Shigella flexneri serotype 3a in men who have sex with men in England and Wales, data from 2009–2011. *Euro surveillance: bulletin europeen sur les maladies transmissibles = European communicable disease bulletin* **17** (2012).

6. Holt, K. E. et al. Shigella sonnei genome sequencing and phylogenetic analysis indicate recent global dissemination from Europe. *Nature genetics* **44**, 1056–1059, https://doi.org/10.1038/ng.2369 (2012).

7. Baker, K. S. et al. Whole genome sequencing of Shigella sonnei through PulseNet Latin America and Caribbean: advancing global surveillance of foodborne illnesses. *Clinical microbiology and infection: the official publication of the European Society of Clinical Microbiology and Infectious Diseases*, https://doi.org/10.1016/j.cmi.2017.03.021 (2017).

8. Connor, T. R. et al. Species-wide whole genome sequencing reveals historical global spread and recent local persistence in Shigella flexneri. *Elife* **4**, e07335, https://doi.org/10.7554/eLife.07335 (2015).

9. Holt, K. E. et al. Tracking the establishment of local endemic populations of an emergent enteric pathogen. *Proceedings of the National Academy of Sciences of the United States of America* **110**, 17522–17527, https://doi.org/10.1073/pnas.1308632110 (2013).

10. Baker, K. S. et al. Travel- and Community-Based Transmission of Multidrug-Resistant Shigella sonnei Lineage among International Orthodox Jewish Communities. *Emerging infectious diseases* **22**, 1545–1553, https://doi.org/10.3201/eid2209.151953 (2016).

11. Chung The, H. et al. Introduction and establishment of fluoroquinolone-resistant Shigella sonnei into Bhutan. *Microb Genom* **1**, e000042, https://doi.org/10.1099/mgen.0.000042 (2015).

12. Zhang, N. et al. Genomic Portrait of the Evolution and Epidemic Spread of a Recently Emerged Multidrug-Resistant Shigella flexneri Clone in China. *Journal of clinical microbiology* **52**, 1119–1126, https://doi.org/10.1128/JCM.02669-13 (2014).

13. Baker, K. S. et al. Intercontinental dissemination of azithromycin-resistant shigellosis through sexual transmission: a cross-sectional study. *The Lancet infectious diseases* **15**, 913–921, https://doi.org/10.1016/S1473-3099(15)00002-X (2015).

14. Chung The, H. et al. South Asia as a Reservoir for the Global Spread of Ciprofloxacin-Resistant Shigella sonnei: A Cross-Sectional Study. *PLoS medicine* **13**, e1002055, https://doi.org/10.1371/journal.pmed.1002055 (2016).
15. Njamkepo, E. et al. Global phylogeography and evolutionary history of Shigella dysenteriae type 1. *Nat Microbiol* 1, 16027, https://doi.org/10.1038/nmicrobiol.2016.27 (2016).

16. Baker, K. S. et al. The extant World War 1 dysentery bacillus NCTC1: a genomic analysis. *Lancet* **384**, 1691–1697, https://doi.org/10.1016/S0140-6736(14)61789-X (2014).

17. Turner, S. A., Luck, S. N., Sakellaris, H., Rajakumar, K. & Adler, B. Nested deletions of the SRL pathogenicity island of Shigella flexneri 2a. *Journal of bacteriology* **183**, 5535–5543, https://doi.org/10.1128/JB.183.19.5535-5543.2001 (2001).

18. World Helath Organisation, Guidelines for the control of shigellosis, including epidemics due to Shigella dysenteriae type 1. 70 (WHO Document Production Services, Geneva, Switzerland, 2005).

19. Gaudreau, C. et al. Ciprofloxacin-resistant Shigella sonnei among men who have sex with men, Canada, 2010. *Emerging infectious diseases* **17**, 1747–1750, https://doi.org/10.3201/eid1709.102034 (2011).

20. Dallman, T. J. et al. Use of whole-genome sequencing for the public health surveillance of Shigella sonnei in England and Wales, 2015. *Journal of medical microbiology* **65**, 882–884, https://doi.org/10.1099/jmm.0.000296 (2016).

21. 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. *Euro surveillance: bulletin europeen sur les maladies transmissibles = European communicable disease bulletin* **20** (2015).

22. Gilbart, V. L. et al. Sex, drugs and smart phone applications: findings from semistructured interviews with men who have sex with men diagnosed with Shigella flexneri 3a in England and Wales. *Sex Transm Infect* **91**, 598–602, https://doi.org/10.1136/sextrans-2015-052014 (2015).

23. Baker, K. S. et al. Horizontal antimicrobial resistance transfer drives epidemics of multiple Shigella species. *Nature communications* In press (2018).

24. Wei, J. et al. Complete genome sequence and comparative genomics of Shigella flexneri serotype 2a strain 2457T. *Infection and immunity* **71**, 2775–2786 (2003).

25. Croucher, N. J. et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. *Nucleic acids research* **43**, e15, https://doi.org/10.1093/nar/gku1196 (2015).

26. Page, A. J. et al. SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. *Microb Genom* **2**, e000056, https://doi.org/10.1099/mgen.0.000056 (2016).

27. Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. *Bioinformatics* **22**, 2688–2690, https://doi.org/10.1093/bioinformatics/btl446 (2006).

28. Zankari, E. et al. Identification of acquired antimicrobial resistance genes. *The Journal of antimicrobial chemotherapy* **67**, 2640–2644, https://doi.org/10.1093/jac/dks261 (2012).

29. Page, A. J. et al. Robust high-throughput prokaryote de novo assembly and improvement pipeline for Illumina data. *Microbial Genomics* **2**, https://doi.org/10.1099/mgen.0.000083 (2016).

30. Rew, V. et al. Whole genome sequencing revealed concurrent outbreaks of shigellosis in the English Orthodox Jewish Community caused by multiple importations of Shigella sonnei from Israel. *Microbial Genomics* **4**, https://doi.org/10.1099/mgen.0.000170 (2018).

31. Miranda, A. et al. Emergence of Plasmid-Borne dfrA14 Trimethoprim Resistance Gene in Shigella sonnei. *Front Cell Infect Microbiol* **6**, 77, https://doi.org/10.3389/fcimb.2016.00077 (2016).

32. Yang, F. et al. Genome dynamics and diversity of Shigella species, the etiologic agents of bacillary dysentery. *Nucleic acids research* **33**, 6445–6458, https://doi.org/10.1093/nar/gki954 (2005).

33. Ye, C. et al. Emergence of a new multidrug-resistant serotype X variant in an epidemic clone of Shigella flexneri. *Journal of clinical microbiology* **48**, 419–426, https://doi.org/10.1128/JCM.00614-09 (2010).

## Acknowledgements

The authors are grateful to the administrative, sequencing and pathogen informatics teams at the Wellcome Trust Sanger Institute and helpful scientists who commented on the work in progress during formal and informal presentations. The authors also thank John Were for his contribution to managing epidemiological data. This work was supported by Wellcome Trust grant number 206194, and KSB is supported by a Wellcome Trust Clinical Research Career Development Fellowship (106690/A/14/Z).

## Author Contributions

K.S.B., T.D., N.F., G.H., C.J., N.R.T. conceived and developed study. K.S.B. and N.R.T. secured funding. K.S.B. directed research, performed analyses and wrote the manuscript drafts. T.D., N.F., T.C., G.H. and C.J. subsampled prepared English isolates for sequencing. F.X.W. and S.L. selected, susceptibility tested and prepared French isolates for sequencing. N.F., H.M., G.H. identified UK patients with repeat infections. M.T. conducted epidemiological interviews of French patients. M.D. conducted susceptibility testing of English isolates. All authors contributed to and approved the final manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-25764-3.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**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) 2018