Skip to main content
Advertisement
  • Loading metrics

Informal genomic surveillance of regional distribution of Salmonella Typhi genotypes and antimicrobial resistance via returning travellers

  • Danielle J. Ingle ,

    Roles Conceptualization, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review & editing

    danielle.ingle@anu.edu.au

    Affiliations National Centre for Epidemiology and Population Health, The Australian National University, Canberra, Australia, Microbiological Diagnostic Unit Public Health Laboratory, Department of Microbiology and Immunology, The University of Melbourne at The Peter Doherty Institute for Infection and Immunity, Melbourne, Australia

  • Satheesh Nair,

    Roles Conceptualization, Data curation, Investigation, Writing – original draft, Writing – review & editing

    Affiliation Gastrointestinal Bacteria Reference Unit, Bacteriology Reference Department, National Infection Service, Public Health England, London, United Kingdom

  • Hassan Hartman,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Gastrointestinal Bacteria Reference Unit, Bacteriology Reference Department, National Infection Service, Public Health England, London, United Kingdom

  • Philip M. Ashton,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Gastrointestinal Bacteria Reference Unit, Bacteriology Reference Department, National Infection Service, Public Health England, London, United Kingdom

  • Zoe A. Dyson,

    Roles Formal analysis, Investigation, Writing – review & editing

    Affiliations Department of Medicine, School of Clinical Medicine, University of Cambridge, Cambridge, United Kingdom, Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Victoria, Australia

  • Martin Day,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Gastrointestinal Bacteria Reference Unit, Bacteriology Reference Department, National Infection Service, Public Health England, London, United Kingdom

  • Joanne Freedman,

    Roles Data curation, Investigation

    Affiliation Travel Migrant and Health Department, Public Health England, 61 Colindale Avenue, London, United Kingdom

  • Marie A. Chattaway,

    Roles Conceptualization, Investigation, Resources, Writing – original draft, Writing – review & editing

    Affiliation Gastrointestinal Bacteria Reference Unit, Bacteriology Reference Department, National Infection Service, Public Health England, London, United Kingdom

  • Kathryn E. Holt ,

    Contributed equally to this work with: Kathryn E. Holt, Timothy J. Dallman

    Roles Conceptualization, Formal analysis, Investigation, Resources, Writing – original draft, Writing – review & editing

    Affiliations Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, University of Melbourne, Parkville, Victoria, Australia, London School of Hygiene and Tropical Medicine, Keppel Street, London, United Kingdom, Department of Infectious Diseases, Central Clinical School, Monash University, Melbourne, Victoria, Australia

  • Timothy J. Dallman

    Contributed equally to this work with: Kathryn E. Holt, Timothy J. Dallman

    Roles Conceptualization, Data curation, Investigation, Resources, Writing – original draft, Writing – review & editing

    Affiliation Gastrointestinal Bacteria Reference Unit, Bacteriology Reference Department, National Infection Service, Public Health England, London, United Kingdom

Abstract

Salmonella enterica serovar Typhi (S. Typhi) is the causative agent of typhoid fever, a systemic human infection with a burden exceeding 20 million cases each year that occurs disproportionately among children in low and middle income countries. Antimicrobial therapy is the mainstay for treatment, but resistance to multiple agents is common. Here we report genotypes and antimicrobial resistance (AMR) determinants detected from routine whole-genome sequencing (WGS) of 533 S. Typhi isolates referred to Public Health England between April 2014 and March 2017, 488 (92%) of which had accompanying patient travel information obtained via an enhanced surveillance questionnaire. The majority of cases involved S. Typhi 4.3.1 (H58) linked with travel to South Asia (59%). Travel to East and West Africa were associated with genotypes 4.3.1 and 3.3.1, respectively. Point mutations in the quinolone resistance determining region (QRDR), associated with reduced susceptibility to fluoroquinolones, were very common (85% of all cases) but the frequency varied significantly by region of travel: 95% in South Asia, 43% in East Africa, 27% in West Africa. QRDR triple mutants, resistant to ciprofloxacin, were restricted to 4.3.1 lineage II and associated with travel to India, accounting for 23% of cases reporting travel to the country. Overall 24% of isolates were MDR, however the frequency varied significantly by region and country of travel: 27% in West Africa, 52% in East Africa, 55% in Pakistan, 24% in Bangladesh, 3% in India. MDR determinants were plasmid-borne (IncHI1 PST2 plasmids) in S. Typhi 3.1.1 linked to West Africa, but in all other regions MDR was chromosomally integrated in 4.3.1 lineage I. We propose that routine WGS data from travel-associated cases in industrialised countries could serve as informal sentinel AMR genomic surveillance data for countries where WGS is not available or routinely performed.

Author summary

Our data demonstrate how routine WGS data produced by Public Health England can be further mined for informal passive surveillance of Salmonella Typhi circulating in different geographical regions where typhoid is endemic. We have shown the public health utility of a simplified approach to WGS reporting based on the GenoTyphi genotyping framework and nomenclature, which doesn’t require the generation of a phylogenetic tree or other phylogenetic analysis. These approaches yielded results consistent with previously reported antimicrobial resistance (AMR) patterns of S. Typhi, including prevalence of multi-drug resistant (MDR) and fluoroquinolone resistance in different regions in association with different pathogen variants. These data provide a rationale and framework for the extraction and reporting of geographically stratified genotype and AMR data from public health labs in non-endemic countries. Prospective analysis and reporting of such data could potentially detect shifts in regional S. Typhi populations, such as replacement or spread of different subclades and the emergence and dissemination of MDR, fluoroquinolone resistant and/or extensively drug resistant S. Typhi, providing valuable data to inform typhoid control measures in low and middle income countries that are still building their genomics capacity.

Introduction

Salmonella enterica serovar Typhi (S. Typhi) is the causative agent of typhoid fever[1], a systemic human infection responsible for an estimated 223,000 deaths each year resulting from more than 20 million infections[2,3]. The majority of the disease burden falls on children in low and middle income countries (LMICs)[2], however vaccination programmes are rare in endemic countries and antimicrobial therapy is considered crucial for the safe clearance of S. Typhi infections and the avoidance of clinical complications. Historically typhoid fever could be effectively treated using first-line drugs including chloramphenicol, ampillicin or co-trimoxazole. However the emergence of multi-drug resistant (MDR) S. Typhi, defined as displaying resistance to all three drugs, in the late 1980s and early 1990s resulted in changes to treatment guidelines, with fluoroquinolones becoming the recommended therapy[47]. Recent years have seen a rise in the proportion of S. Typhi disease isolates displaying reduced susceptibility to fluoroquinolones associated with point mutations in quinolone resistance determining region (QRDR) of gyrA and parC[2,8,9].

S. Typhi is genetically monomorphic, which has historically constrained molecular surveillance of S. Typhi before the era of high throughput whole genome sequencing (WGS)[5,10]. A phylogenetically informative genotyping scheme, GenoTyphi, was recently introduced to facilitate the interpretation of S. Typhi WGS data[11]. Application of the scheme to a global collection of isolates from >60 countries showed that the S. Typhi population is highly structured, comprising dozens of subclades associated with specific geographical regions[11,12]. This global genomic framework for S. Typhi revealed the majority of MDR S. Typhi infections worldwide are associated with a single phylogenetic lineage, designated as genotype 4.3.1 (H58 under the legacy scheme[5]), which has disseminated from South Asia since the 1990s, including into East Africa [11,12]. Contrastingly, in West Africa MDR S. Typhi is associated with a different genotype, 3.1.1[13,14]. The MDR phenotype in both S. Typhi 4.3.1 and 3.1.1 is encoded by a composite transposon carrying genes conferring resistance to five drug classes. These include chloramphenicol, penicillins and co-trimoxazole and the determinants are typically located on IncH1 plasmids (plasmid sequence type PST6 in 4.3.1, and PST2 in 3.1.1) [4,1217]. Recently, migration of the MDR transposon to the S. Typhi chromosome via IS1 transposition has been reported, with four different integration sites identified in the S. Typhi 4.3.1 strain background (fbp, yidA, STY4438 and the intergenic region between cyaA and cyaY)[12,1821].

Resistance to newer agents is also on the rise in S. Typhi. Reduced susceptibility to fluoroquinolones via QRDR mutations is commonly observed amongst S. Typhi 4.3.1[12], especially lineage II[16,22], but is rare in 3.1.1 and other genotypes[12,13]. S. Typhi 4.3.1 QRDR triple mutants (carrying two mutations in gyrA and one in parC) displaying complete resistance to ciprofloxacin are increasing in prevalence in South Asia[23]. WGS analysis showed that recent treatment failure with gatifloxacin in Nepal resulted from introduction of a 4.3.1 lineage II triple mutant from India, prompting a reconsideration of the reliance on fluoroquinolones in the region[8,22]. Subsequently an outbreak of S. Typhi 4.3.1 displaying resistance to ceftriaxone in addition to ciprofloxacin and all three first line drugs was reported in Pakistan[24,25]. The outbreak strain carried both the MDR composite transposon integrated in the chromosome at yidA, and an E. coli IncY plasmid harbouring the extended spectrum beta-lactamase (ESBL) gene blaCTX-M-15 (conferring resistance to ceftriaxone and other third generation cephalosporins) and the quinolone resistance gene qnrS (which combined with a gyrA-S83F mutation in the chromosome conferred resistance to ciprofloxacin)[24]. This strain has been designated XDR (extensively drug resistant, defined as resistance to chloramphenicol, penicillins, co-trimoxazole, ceftriaxone and ciprofloxacin) and severely limits treatment options, with azithromycin being the last remaining oral antibiotic (to which sporadic resistance has already been observed in South Asia)[7].

These dynamic trends highlight a need for prospective AMR surveillance in global S. Typhi populations, in order to inform empirical treatment options; genomic surveillance offers the added benefit of revealing resistance mechanisms and regional and international spread of emerging MDR and XDR strains[26]. Blood culture confirmation and isolate characterisation is a pre-requisite for such activities, but neither are routinely performed in laboratories in areas where typhoid fever is endemic. However, in many developed nations, S. Typhi infections are notifiable and the disease is primarily associated with travellers returning from high-risk regions[27], providing an opportunity for informal sentinel surveillance of those regions. In England, typhoid fever is notifiable and all S. Typhi isolates are sent to Public Health England (PHE) for confirmation and characterisation via WGS, and recent travel history is sought via enhanced surveillance questionnaires. Most cases identified in England (~80%) are associated with recent travel to typhoid endemic areas[27]. We recently reported that the geographic origin of travel-associated S. Typhi cases in London could be predicted by comparing WGS data to the global framework[11], and identified the first case of ESBL S. Typhi in England from a patient with travel to Pakistan[28], supporting the notion that WGS data on travel-associated cases has sentinel surveillance value. Here we analysed WGS data from S. Typhi isolates referred to PHE from three years of national surveillance between April 2014 and March 2017[29], and explored the distribution of lineages and AMR determinants in geographic regions frequented by UK travellers.

Methods

Bacterial isolates included in this study

A total of 533 S. Typhi isolates from English cases received by PHE during the period of 1st April 2014 to 31st March 2017 were included in this study. Ethical approval for the detection of gastrointestinal bacterial pathogens from faecal specimens, or the identification, characteristation and typing of cultures of gastrointestinal pathogens, submitted to the Gastrointestinal Bacteria Reference Unit is not required as covered by Public Health England’s surveillance mandate. Patient travel information was available for 488/533 of the isolates and was obtained by PHE using an enhanced surveillance questionnaire (https://www.gov.uk/government/publications/typhoid-and-paratyphoid-enhanced-surveillance-questionnaire), this included questions pertaining to the destination of any foreign travel that occurred during the likely incubation period (28 days before onset of symptoms). No specific consent was required from the patients whose data were used in this analysis because Public Health England 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. For the remaining 45 isolates, no enhanced surveillance questionnaire was completed or the data collected was incomplete; these were spread across the four years (n = 15 for 2014, n = 24 for 2015, n = 4 for 2016 and n = 2 for 2017; mean 8.5%). Details of all 533 isolates are provided in S1 Table.

Whole genome sequencing

WGS was conducted as part of routine sequencing of all Salmonella isolates referred to the Gastrointestinal Bacteria Reference Unit at PHE, as previously described[29]. Briefly, DNA was fragmented and tagged for multiplexing with NexteraXT DNA Sample Preparation Kits (Illumina) and sequenced at PHE on a HiSeq 2500 yielding 100 bp paired end reads. FASTQ data is available from the NCBI Short Read Archive, BioProject accession PRJNA248792. Individual accession numbers for isolates analysed in this study are given in S1 Table.

Single nucleotide variant (SNV) analysis and in silico genotyping

Paired end Illumina reads were mapped to the CT18 reference genome (accession AL513382)[30], which is the standard reference for S. Typhi genomic studies, using RedDog (V1.beta10.3) available at (https://github.com/katholt/RedDog). Briefly, RedDog maps reads to the reference genome using Bowtie2 (v.2.2.9)[31], before using SAMtools (v1.3.1)[32] to identify high quality single nucleotide variant (SNV) calls as previously described[19]. A core SNV alignment was generated by concatenating alleles with high-quality consensus base calls (phred score >20), for all SNV sites that had such calls in >95% of genomes (this represents the 95% ‘soft’ core of the S. Typhi genomes); this alignment was filtered to exclude SNVs in phage regions and repetitive sequences (354 kb; ~7.4% of bases) in the CT18 reference chromosome as defined previously[10]; S2 Table) and recombinant regions identified by Gubbins[33]. The resulting alignment of 8053 SNVs was used as input to RAxML (v8.2.8) to infer a maximum likelihood (ML) phylogeny with a generalised time-reversible model and a Gamma distribution to model site-specific rate variation (GTR+ Γ substitution model; GTRGAMMA in RAxML), and 100 bootstrap pseudo-replicates to assess branch support.

Genotypes were inferred for all isolates by screening the Bowtie2 alignment (bam) files for SNVs used in the extended S. Typhi typing framework, GenoTyphi (code available at https://github.com/katholt/genotyphi)[11]. GenoTyphi uses 72 specific SNVs to assign isolates to one of four primary clusters; 16 clades; and 49 subclades[11], with the globally disseminated 4.3.1 (H58) subclade further delineated into lineages I and II (4.3.1.1 and 4.3.1.2).

In silico characterisation of AMR associated genes and mobile elements

S. Typhi genomes were screened for acquired AMR genes and plasmid replicons using PHE’s mapping-based Genefinder pipeline as previously described[34], which has been extensively validated in comparing phenotypic and WGS-derived AMR profiles on ~600 typhoidal Salmonella Typhi and Paratyphi [35]. Reference sequences in Genefinder were curated from those described in Comprehensive Antimicrobial Resistance Database[36] and PlasmidFinder[37]. Genes were called as present within a genome when detected with 100% coverage and >90% nucleotide identity to the reference gene. Genefinder was also used to detect point mutations in the QRDR of chromosomal genes gyrA (codons 83, 87) and parC (codons 79, 80, 84)[35,38]. Isolates were defined as being MDR if genes were detected by Genefinder in the Beta-Lactamases, Trimethoprim, Sulphonamides and Chloramphenicol classes.

ISmapper[39] (v2) was run with default parameters to screen all read sets for insertion sites of the transposase IS1 (accession X52534) relative to the CT18 reference chromosome sequence (accession AL513382). The binary data from ISmapper was processed in R using tidyverse (v1.2.1) (https://cran.r-project.org/package=tidyverse). For six isolates that carried MDR genes but no plasmid replicon genes, ISMapper was re-run with the minimum read depth threshold reduced to 1 read (from the default of 6 reads), in order to increase sensitivity to detect the associated IS1 site. This identified IS1 sites in five of the six isolates (SRR3049053, SRR5989319, SRR1967049, SRR5500465 and SRR5500455) but not in isolate SRR5500435.

The presence and subtypes of IncHI1 and IncN plasmids was further investigated using plasmid MLST (pMLST). Publicly available pMLST schemes for IncHI1[40] and IncN[41], available at https://pubmlst.org/plasmid/, were used to screen the relevant read sets using SRST2 (v0.1.8)[42]. The five isolates that were determined to have pST2 versions of the IncHI1 plasmid were then mapped to the reference sequence of plasmid pAKU1 (accession AM412236) with RedDog, using the same quality control thresholds as described for the chromosomal SNVs. The coverage level of the plasmid was addressed determining the ratio of average sequencing depth for the chromosome and the plasmid for all five isolates, with the ratio in the range 0.89–2.29 in all cases. A pST1 plasmid sequence (pUI1203_01 accession ERR340785) from S. Typhi strain UI1203 (genotype 3.2.1, isolated from Laos in 2001)[11] was included to provide an outgroup for tree rooting. The resulting SNV alignment was filtered to exclude non-backbone regions (66,458 bp) of the pAKU1 plasmid as previously defined [4] of the plasmid, and to exclude SNV sites present in <95% of plasmid sequences, and recombinant regions detected by Gubbins (v2.3.2), resulting in a final alignment of 17 SNVs for phylogenetic inference as above.

Genomes with less common AMR profiles were assembled using Unicycler v0.4.6 [43]. The assemblies were further interrogated for the presence of plasmid replicon genes using ABRicate (https://github.com/tseemann/abricate) and the PlasmidFinder database using a minimum coverage and minimum identity of 90. The location of AMR genes detected, IS1 sequences and IncHI1 plasmid backbone genes was visualised in Bandage v0.8.1 [44].

Trees annotated with MDR, IS1 insertion sites, plasmid replicons and QRDR point mutations were visualised using ggtree v1.8.1[45]. An interactive version of the S. Typhi ML phylogeny with associated metadata is found at microreact (https://microreact.org/project/r1njJ_qJ4).

Ciprofloxacin susceptibility phenotyping

The minimum inhibitory concentration (MIC) of ciprofloxacin was measured for 173 isolates (n = 117 from 2015 and n = 56 from 2016) as part of a previously reported study[35]. Briefly, MIC was determined by agar dilution using Mueller-Hinton agar, and S. Typhi were interpreted as displaying resistance (MIC ≥0.5 μg/mL) or susceptibility to ciprofloxacin (MIC <0.06 μg/mL), according to the EUCAST guidelines (v7.1)[35]. For the 173 phenotyped isolates, the total number of QRDR mutations per isolate was plotted against the reported MIC, in R using ggplot2 (v3.0.0)[46].

Results

Countries of travel

A total of 533 S. Typhi isolates were referred from English cases to PHE during the three-year period between April 2014 and March 2017 (n = 178, n = 176, n = 179 in each year running from April to March; listed in S1 Table). Of these, 449 cases (84.2%) reported recent foreign travel (within 28 days of onset of symptoms), spanning 26 countries. A further 39 cases (7.3%) reported no travel abroad. The proportion of cases reporting no recent travel, which may represent local transmission within the UK, displayed a non-significant increase from 5.6% in the first year to 8.9% in the third year (p = 0.5 using Chi-squared trend test). No information on travel was available for the remaining 45 cases (8.4%). The distribution of travel destinations for the cases with known travel history are given in Table 1 and Fig 1. The majority of cases (n = 387, 73%) reported travel to South Asia, particularly India (36%), Pakistan (29%) and Bangladesh (6%), reflecting travel and migration patterns between England and typhoid endemic areas. Annual case numbers from India and Bangladesh were constant (mean 64 and 11 per year, respectively; p = 0.2 in each case using Chi-squared trend test), but case numbers from Pakistan spiked significantly in the third year of the study, to n = 68 (38% of total cases) compared to n = 45 and n = 43 (25%) in the first two years (p = 0.01 using Chi-squared trend test). Small numbers of cases reported travel to South East Asia (n = 10, 2%), East Africa (n = 21, 4%), West Africa (n = 11, 2%) and the Middle East (n = 4, 1%; see Table 1). Five cases (1%) had reported travel to Europe or North America or South America only.

thumbnail
Fig 1. Summary of S. Typhi received at Public Health between April 2014 to March 2017.

a) Map highlighting countries that were reported as travel destinations for patients from which S. Typhi cultures were collected. Countries are coloured by major geographical regions. Blank map sourced from https://commons.wikimedia.org/wiki/File:BlankMap-FlatWorld6.svg b) Number of S. Typhi isolates reporting travel to each country over the three-year period of study, stratified by date of receipt at PHE.

https://doi.org/10.1371/journal.pntd.0007620.g001

S. Typhi genotypes by region

The 533 genomes were assigned to 31 unique S. Typhi genotypes using the GenoTyphi scheme (Fig 2). The majority (n = 391, 73.4%) belonged to subclade 4.3.1 (H58), which was identified in cases with travel to 14 different countries. Genotype 4.3.1 dominated amongst cases with reported travel to the South Asian countries associated with the greatest burden of travel-related cases (India, 79% 4.3.1; Pakistan, 92% 4.3.1; Bangladesh, 52% 4.3.1), and amongst cases with travel to East Africa (86% 4.3.1; see Table 1). Of the 4.3.1 isolates, 162 (41%) were further classified to lineage I (4.3.1.1), including 66% of the 4.3.1 isolates from cases with travel to Pakistan, 94% of those with travel to Bangladesh, and all 4.3.1 isolates from cases with travel to Mozambique, Tanzania and Zimbabwe (total n = 13). In contrast, lineage II (4.3.1.2) accounted for 156 (40%) of 4.3.1 genotypes, dominating in returning travellers to India (75% of 4.3.1 isolates) and accounting for the three 4.3.1 isolates associated with travel to Rwanda and Uganda (Fig 3). Notably, 4.3.1 was not detected in cases associated with travel to West Africa (0/11; see Table 1, Fig 3). Other common genotypes include clade 3.3 (6% of total isolates), found in cases with travel to India (9% of all isolates from this location), or Bangladesh (15%); clade 2.2 (3.1% of total isolates), found in cases with travel to India (2%) or Pakistan (4.5%); subclade 3.2.2 (2.6% of total isolates), found in 21% of cases with travel to Bangladesh; and subclade 3.1.1 (1.7% of total isolates), found in 7/11 (63%) of cases with travel to West Africa, including Nigeria (n = 4/6), Ghana (2/2) and Sierra Leone (1/1).

thumbnail
Fig 2. Membership of S. Typhi by genotypes and association with AMR.

a) Population framework tree of the S. Typhi genotypes sequenced at PHE. Background shading indicates the primary clades. Tree tips indicate individual genotypes assigned by GenoTyphi; tip sizes and colours indicate number of isolates belonging to each genotype, these values (n) are also printed to the right. b) Frequency of MDR for each genotype (defined based on presence of AMR genes encoding resistance to the first line drugs ampicillin, trimethoprim-sulfamethoxazole, or chloramphenicol). c) Frequency of QRDR mutations by genotype.

https://doi.org/10.1371/journal.pntd.0007620.g002

thumbnail
Fig 3. Geographic distribution of S. Typhi based on reported country of travel.

S. Typhi isolates with completed travel questionnaires and reporting travel to only one country are shown. The number in each circle indicates the total number of cases reporting travel to that country, circle graphs indicate the genotype composition of the corresponding S. Typhi isolates, stratified by MDR, coloured as per inset legend. Blank map sourced from https://commons.wikimedia.org/wiki/File:BlankMap-FlatWorld6.svg

https://doi.org/10.1371/journal.pntd.0007620.g003

Multi-drug resistance and QRDR mutations by region

The overall frequency of MDR, defined as carriage of genes associated with resistance to ampicillin (blaTEM), chloramphenicol (catA1) and co-trimoxazole (sul1 or sul2 plus a dfrA gene), was 129 (24%). MDR S. Typhi was associated with reported travel to 10 countries from three main regions, South Asia (26%), East Africa (52%) and West Africa (27%), with high frequencies amongst isolates whose cases report travel to Zimbabwe (8/10, 80%), Nigeria (2/6, 33%), Tanzania (2/4, 50%), Ghana (1/2, 50%), Bangladesh (8/33, 24%) and Pakistan (85/156, 55%) (Table 1, Fig 3). MDR was also present at low frequency amongst cases with reported travel to India (3%), and singleton MDR isolates were detected from Mozambique, United Arab Emirates and Afghanistan. Notably the majority of MDR isolates were linked with travel to Pakistan (n = 85, 66% of all MDR isolates) and other South Asian countries (n = 100, 78% of all MDR isolates; see Table 1). No ESBL or carbapenemase genes were detected.

All MDR isolates belonged to either 4.3.1 (n = 125) or 3.1.1 (n = 4, see Fig 2b). The frequency of MDR was highest in 4.3.1.1 (96%), then 3.1.1 (3%), and one isolate in 4.3.1 (Figs 2b and 3). At the regional level, MDR was common amongst cases associated with travel to East Africa (52%), West Africa (27%) and South Asia (26%), but not with travel elsewhere (with the exception of one isolate from UAE; see Table 1). However, there were significant differences between countries within these regions, associated with differences in the dominant S. Typhi clades (Fig 3). In East Africa and South Asia, MDR was detected amongst cases associated with travel to the countries dominated by MDR-associated lineage 4.3.1 lineage I (Mozambique, Tanzania and Zimbabwe in East Africa; and Bangladesh and Pakistan in South Asia; see Table 1, Fig 3). In West Africa, MDR was detected in three isolates (n = 2 Nigeria, n = 1 Ghana), all belonging to the region’s dominant subclade 3.1.1 (a fourth MDR 3.1.1 infection had no reported travel).

QRDR mutations were identified in 455 isolates (85.4%) originating from 18 countries (Table 1, S1 Table). QRDR mutants were most frequent in India (97%), Bangladesh (94%), Nepal (66%), Pakistan (94%), Myanmar (100%), Uganda (100%), and Nigeria (50%). Singleton QRDR mutants were also identified in China, Malaysia, Afghanistan, United Arab Emirates, Egypt, Rwanda, Peru, Greece and Zimbabwe. QRDR triple mutants were identified only in isolates associated with travel to South Asia: 33% of isolates from Nepal, 23% of those from India and 4% of those from Pakistan (Table 1). Overall, the presence of QRDR mutations was very common in cases associated with travel to South Asia (95% of all isolates linked to this region), and significantly less common (p < 2.2e-16 using Chi-squared trend test) in East Africa (43%,) and West Africa (27%; see Table 1).

Ciprofloxacin MICs have been previously reported for 173 of the S. Typhi isolates (S1 Table)[35], and a comparison of QRDR mutations with these phenotypes is shown in Fig 4 to facilitate the interpretation of QRDR mutations. These data showed that all isolates carrying a single QRDR mutation had ciprofloxacin MIC of at least 0.064 μg/mL (exceeding the EUCAST threshold for susceptibility); all those with two QRDR mutations had MIC of 0.125 μg/mL; and all those with three mutations had MIC ≥1 μg/mL (above the threshold for resistance) (Fig 4). The most common mutations were gyrA-S83F (n = 271, 51%) and gyrA-S83Y (n = 83, 16%), which were found in 16 genotypes in cases associated with travel to eleven and six countries respectively (Table 2). Double mutants (gyrA-S38F or -S83Y plus a mutation in parC), associated with elevation of ciprofloxacin MIC to ≥0.25 μg/mL (Fig 4), were found in 16 isolates (3.0%) associated with travel to India or Bangladesh. A total of 61 isolates were identified as QRDR triple mutants, which display resistance to ciprofloxacin (MIC >1 μg/mL, see Fig 4) and have been associated with fluoroquinolone treatment failure[8]. Most of these (n = 59) carried gyrA-S83F, gyrA -D87N and parC-S80I mutations, while the remaining two had a unique profile of gyrA-S83F, gyrA-D87G and parC-S80I and of gyrA-S83FY, parC-Y74X and parC-P98X (Table 2).

thumbnail
Table 2. Combinations of coding changes detected in the QRDR and their frequency by country of reported travel.

https://doi.org/10.1371/journal.pntd.0007620.t002

thumbnail
Fig 4. Minimum inhibitory concentration (MIC) for ciprofloxacin versus number of point mutations in QRDRs in two chromosomal genes.

The number of isolates with different combinations of MIC values and point mutations detected in the QRDRs of two chromosomal genes, gyrA and parC, are shown. The breakpoints of susceptible, reduced susceptibility and resistance are shown by the background gradient of the plot. The triangle indicates an isolate that carries qnrS-1.

https://doi.org/10.1371/journal.pntd.0007620.g004

All QRDR double and triple mutants belonged to genotype 4.3.1 and were associated with travel to South Asia (or no/unknown travel). Overall, QRDR mutations were significantly more common in 4.3.1 lineage II compared to lineage I (99% vs 91%, p = 0.006, two-sided test of difference in proportions). The 61 QRDR triple mutants were detected only in S. Typhi 4.3.1 lineage II isolates (Fig 2c) and were associated with either travel to India (n = 43, 22.5% of all isolates from this location, where lineage II is dominant); travel to neighbouring countries Pakistan (n = 6) or Nepal (n = 1); travel to multiple destinations including India (n = 5); no reported travel (n = 2); or no travel information available (n = 4) (Tables 1 and 2). Notably because all 4.3.1 MDR isolates belonged to lineage I, and all QRDR triple mutants belonged to lineage II (Fig 2), there were no isolates with both MDR and three QRDR mutations (Fig 2, S1 Table). There were however 119 cases with MDR plus 1–2 QRDR mutations (i.e. reduced susceptibility to fluoroquinolones, Fig 4). The vast majority of these were 4.3.1 lineage I (n = 115), most commonly associated with travel to Pakistan (n = 85, 54.4% of isolates from this country), Bangladesh (n = 8, 24%) and India (n = 6, 3%) but also occasional cases who reported travel to Zimbabwe (n = 2), Tanzania (n = 2) and United Arab Emirates (n = 1). Three 3.1.1 isolates were also MDR with one QRDR mutation (n = 2 travel to Nigeria, n = 1 with no recorded travel).

Plasmid vs chromosomal location of AMR genes in S. Typhi

All MDR isolates belonged to 4.3.1 (n = 125) or 3.1.1 (n = 4) and carried the typical S. Typhi MDR composite transposon comprising Tn6029 (encoding blaTEM-1, sul2, strAB) inserted in Tn21 (carrying a class I integron encoding dfrA alleles in the gene cassette and sul1 at the end), which is in turn inserted within Tn9 (encoding catA1)[47] (see Fig 5a). All 125 MDR 4.3.1 isolates (associated with South Asia and East Africa) carried dfrA7 in the integron cassette and no plasmid replicons (Fig 5a). In most of these (n = 123, 98%), we detected chromosomal IS1 insertions at sites previously associated with IS1-mediated integration of the MDR composite transposon (S1 Fig) (cyaA or yidA sites[12,18]). A putative IS1 insertion was detected in the novel site STY3168 in a single, genome (SRR5500440). Further, of the six isolates with four IS1 sites detected, five had recent travel to Pakistan with last isolate having no reported travel (S1 Fig). Notably, most (93%) of the MDR 4.3.1 isolates also carried a QRDR mutation, the most common being gyrA-S83F (Table 2, S1 Fig). The four MDR 3.1.1 isolates (associated with West Africa) carried IncHI1 PST2 plasmids with dfrA15 in the integron cassette (Fig 5b). An IncHI1 PST2 plasmid was also identified in a single non-MDR 2.3.1 isolate associated with travel to Nigeria. The plasmid backbone was very closely related to that of the 3.1.1 West African plasmids but carried dfrA1 in the integron cassette and lacked the chloramphenicol and ampicillin resistance genes catA1 and blaTEM (Fig 5a and 5b).

thumbnail
Fig 5. Mobilisation of MDR element in S. Typhi.

a) Structure of the MDR composite transposon in S. Typhi strain ERL12960 (4.3.1 lineage I) is shown (accession ERL12960). MDR genes are those encoding resistance to the first line drugs ampicillin, trimethoprim-sulfamethoxazole or chloramphenicol b) Phylogeny showing genetic relationships between the five IncHI1 PST2 S. Typhi isolates identified in this study, based on SNVs identified in the plasmid backbone sequence and rooted using the PST1 plasmid pUI1203_01 (accession ERR340785) as an outgroup. The presence of AMR genes detected in each genome is indicated in the heatmap, with alleles specified for dfrA and blaTEM genes.

https://doi.org/10.1371/journal.pntd.0007620.g005

Plasmid replicon and AMR gene screening identified additional plasmid replicons and AMR genes at low frequency. Two of the MDR 4.3.1.1 isolates from Zimbabwe carried IncN (subtype PST3) plasmids. In addition to the genes typical of the MDR composite transposon (with dfrA7 in the integron), these isolates also carried qnrS, dfrA14 and tet(A). It was not possible to resolve the precise locations of the acquired AMR genes due to the limitations of short read assembly. As these isolates lack QRDR mutations, the presence of the qnrS gene is predicted to confer reduced susceptibility to fluoroquinolones but not full resistance; indeed, one of the isolates (SRR4063811) was phenotyped and displayed reduced susceptibility to ciprofloxacin, with MIC 0.25 μg/mL. IncN plasmids (subtype PST5) were found in three non-MDR 4.3.1.2 isolates that carried dfrA15, sul1 and tet(A). These isolates (two from India, one with no reported travel) were also QRDR triple mutants and thus predicted to be fully resistant to fluoroquinolones. The combination of dfrA14, sul2, blaTEM-1, strA and strB was detected in three 4.3.1.1 isolates with no QRDR mutations (two with travel to Tanzania, one with no travel information). All three carried sequences with similarity (100% nucleotide identity and 50% coverage) to a FIBK plasmid carrying dfrA14, sul2 and blaTEM-1 that was previously sequenced from a 2008 Tanzanian S. Typhi, strain 129–0238 (GenBank accession LT904889)[12]. The same combination of AMR genes (dfrA14, sul2, blaTEM-1, strA and strB) plus tet(A) were also found in a 3.1.1 gyrA-S83Y isolate from Nigeria that harboured an IncY plasmid replicon.

Discussion

Here we demonstrate the utility of using S. Typhi WGS data generated routinely at a public health laboratory in the UK to serve as informal surveillance for different geographical regions where typhoid fever is endemic. The PHE dataset encompasses a diverse collection of 533 S. Typhi isolates from multiple geographical regions collected over a three-year sampling period (Fig 1). As typhoid fever has not been endemic in England since the successful interventions of the major controlling measures, water sanitation and hygiene complemented with antimicrobial therapy[27,48], it is assumed that notified cases to PHE are associated with returned travellers and their contacts. The collection is biased towards isolates from South Asia and East Africa which reflects historical and contemporary ties between England and these regions, however these are also regions that experience a high burden of typhoid fever and could benefit from the AMR and WGS data obtained routinely by PHE. Notably, public health agencies in other countries receive S. Typhi isolates reflecting the distinct travel habits of their own populations (for example Institut Pasteur receives more S. Typhi from travellers visiting Francophone countries in Africa and former French colonies such as Vietnam), and the synthesis of these diverse collections could potentially provide more extensive sentinel surveillance coverage of typhoid endemic regions.

Typing the S. Typhi isolates using the GenoTyphi scheme enabled rapid classification of the 533 genomes into 31 distinct lineages (Fig 2). This simple tree-free approach showed clustering of subclades by geographical region of travel that was consistent with previous previously reported geographical patterns[11], providing support for the use of travel associated isolates as an indicator of local pathogen populations. Notable examples include the detection of MDR subclade 4.3.1 in East Africa and South Asia [8,15], and the complete absence of subclade 4.3.1 in the isolates from West Africa[14] (Fig 3). Genomes from cases reporting travel to West Africa were genotyped as 3.1.1, consistent with earlier studies where 3.1.1 was found to be the main S. Typhi lineage in the region[13,14]. Furthermore, for isolates with travel to more than one country, the genotype can help to discern the most likely origin of the infection; for example, for one case reporting travel to Nigeria and Turkey, the genome isolate (SRR558502) was identified as genotype 3.1.1, suggesting that the pathogen was most likely acquired in Nigeria. This further demonstrates the public health utility of WGS data on returning traveller isolates. Notably the GenoTyphi scheme provides a mechanism and nomenclature for such insights to be achieved simply and rapidly from individual genomes and by different laboratories working independently, without need for phylogenetic tree construction or other comparative analyses.

Data from this study revealed that reduced susceptibility to fluoroquinolones was common amongst S. Typhi associated with diverse geographic sites and genotypes. However, the ciprofloxacin resistant QRDR triple mutants were all from cases belonging to subclade 4.3.1.2 (Fig 2, Table 2), the majority of which had reported travel to South Asia, mainly India. None of these QRDR triple mutants were also MDR, which was associated with cases belonging to subclade 4.3.1.1. This is in line with previous observations, and others have speculated there may be a fitness cost to the carriage of QRDR triple mutants particularly in the MDR background [19]. These data align with previous findings[8,19], and support the hypothesis that high levels of fluoroquinolone exposure in India through healthcare and the environment are driving the emergence of resistant S. Typhi and other pathogens in the region[49] and contributing to treatment failure for typhoid fever[8]. Further, in addition to confirming that QRDR triple mutants are resistant to ciprofloxacin, the MIC data clearly showed that even a single QRDR mutation is associated with reduced ciprofloxacin susceptibility (MIC 0.06–0.25 μg/mL) (Fig 4), which has also been associated with clinical failure[35,50]. Single QRDR mutations are not associated with a fitness cost[51], and are present in MDR isolates, suggesting that this may further limit treatment options. While phenotypic MICs are still advised for clinical treatment, inferring susceptibility to antimicrobials from WGS data is appropriate for surveillance purposes [35].

Strikingly in this data set all 4.3.1 MDR isolates carried the composite transposon integrated in the chromosome, suggesting an enormous shift in the burden of MDR typhoid from plasmid-borne resistance. This is of grave concern as it means that there is likely to be very little fitness cost associated with carriage of the MDR transposon. In the late 1990s the increase in MDR typhoid in Asia prompted a switch to fluoroquinolones for treatment, which in Nepal and other regions was followed by almost complete loss of the MDR plasmid from the S. Typhi population, suggesting that the fitness cost of the plasmid leads to plasmid loss in the absence of selection from the first-line drugs. However, the integration of the MDR transposon into the S. Typhi chromosome likely alleviates any fitness cost, making it more likely that MDR will be maintained even in the absence of selection for the specific resistances encoded.

A key limitation of our MDR element analysis is that the routine genomic surveillance data generated at PHE is restricted to short-read Illumina sequences, which prohibited the full resolution of complex MDR loci and their precise locations in the S. Typhi genome. Nevertheless, using short-read mapping approaches, we were able to confirm that in nearly all MDR isolates the MDR loci could be localised to plasmids or chromosomal integration sites that have been previously fully resolved in S. Typhi using long reads [52,18]. Notably, 98% of plasmid-free MDR isolates showed evidence of chromosomal integration of the MDR element at the known sites yidA or upstream of cyaA in the 4.3.1.1 lineage. The other MDR isolates may harbour the MDR transposon in novel genetic contexts, which could potentially be resolved using long-read sequencing. PHE has demonstrated the utility of this approach previously, including in S. Typhi [28], but does not currently conduct long-read sequencing as a routine part of Salmonella genomic surveillance.

Of particular note was the increase in number of IS1 insertions in from cases in subclade 4.3.1.1 with reported travel to Pakistan (S1 Fig). We may hypothesise that the S. Typhi from this region may be more likely to acquire novel mechanisms in response to local selective pressures. Indeed, the recent acquisition of an IncY plasmid harbouring blaCTX-M-15 and qnrS genes has resulted in the emergence of an XDR lineage of S. Typhi from Pakistan[24] and an IncI1 plasmid encoding blaCTX-M-15 in an S. Typhi from Bangladesh[52] provides evidence for this hypothesis, highlighting the importance of ongoing surveillance of these regions that experience a high burden of typhoid fever.

Currently there are three critical AMR threats posed by S. Typhi, namely the dissemination of mobile AMR genes mediating MDR profiles, the evolution of point mutations in gyrA and parC, two core housekeeping genes, that confer differing levels of fluoroquinolone resistance (Figs 2 and 4), and the recent emergence of XDR S. Typhi. The WGS data presented here provide insight into changing AMR dynamics within S. Typhi. Importantly, the concordance of genomic and phenotypic AMR data for ciprofloxacin resistance in this study (Fig 4) which have been extensively characterised for S. Typhi for multiple drugs previously[35], demonstrates the utility of WGS for robustly characterising AMR profiles. Here, two of the AMR threats where characterised within the S. Typhi collection, reflecting geographical differences in AMR profiles. While no XDR S. Typhi had been detected in the PHE collection between April 2014 and March 2017, the first XDR S. Typhi isolate in a returned traveller with recent travel to Pakistan shortly after XDR S. Typhi were reported from Pakistan has been identified[24,28]. This highlights the value in using these data as informal surveillance of S. Typhi as we hypothesise that resistance to second-line drugs such as azithromycin will arise under continued selective pressure as has occurred with previous drugs.

Supporting information

S1 Fig. Overview of S. Typhi PHE collection.

a) The phylogeny of the 533 S. Typhi isolates in the PHE collection is shown on the left with the geographical region of reported travel. b) The presence of IncHI1 or IncN plasmid replicons is shown with the different colours indicating PST. c) The presence of MDR profiles is shown in grey. The presence of point mutations in QRDR genes is shown with the different gradients of blue indicating different mutations. d) Detected IS1 sites of insertion are shown in grey. The total number of IS detected in each of the S. Typhi isolates is shown on the far right.

https://doi.org/10.1371/journal.pntd.0007620.s001

(PDF)

S1 Table. Data of 533 S. Typhi isolates received and sequenced at Public Health England between April 2014 and March 2017.

https://doi.org/10.1371/journal.pntd.0007620.s002

(XLSX)

S2 Table. Excluded repeat and phage regions in CT18 reference.

https://doi.org/10.1371/journal.pntd.0007620.s003

(XLSX)

References

  1. 1. Parry CM, Hien TT, Dougan G, White NJ, Farrar JJ. Typhoid Fever. N Engl J Med. 2002;347: 1770–1782. pmid:12456854
  2. 2. Crump JA, Mintz ED. Global Trends in Typhoid and Paratyphoid Fever. Clin Infect Dis. 2010;50: 241–246. pmid:20014951
  3. 3. Mogasale V, Maskery B, Ochiai RL, Lee JS, Mogasale VV, Ramani E, et al. Burden of typhoid fever in low-income and middle-income countries: a systematic, literature-based update with risk-factor adjustment. The Lancet Global Health. 2014;2: e570–e580. pmid:25304633
  4. 4. Holt KE, Phan MD, Baker S, Duy PT, Nga TVT, Nair S, et al. Emergence of a globally dominant IncHI1 plasmid type associated with multiple drug resistant Typhoid. PLoS Negl Trop Dis. 2011;5: 1–13.
  5. 5. Roumagnac P, Weill F-X, Dolecek C, Baker S, Brisse S, Chinh NT, et al. Evolutionary history of Salmonella Typhi. Science. 2006;314: 1301–1304. pmid:17124322
  6. 6. Akhtar S, Sarker MR, Jabeen K, Sattar A, Qamar A, Fasih N. Antimicrobial resistance in Salmonella enterica serovar Typhi and Paratyphi in South Asia-current status, issues and prospects. Critical Reviews in Microbiology. 2014;41: 536–545. pmid:24645636
  7. 7. Levine MM, Simon R. The gathering storm: Is untreatable Typhoid Fever on the way. mBio. 2018;9: e00482–18. pmid:29559573
  8. 8. Pham Thanh D, Karkey A, Dongol S, Ho Thi N, Thompson CN, Rabaa MA, et al. A novel ciprofloxacin-resistant subclade of H58 Salmonella Typhi is associated with fluoroquinolone treatment failure. eLIFE. 2016;5: e14003. pmid:26974227
  9. 9. Chiou C-S, Lauderdale T-L, Phung DC, Watanabe H, Kuo J-C, Wang P-J, et al. Antimicrobial resistance in Salmonella enterica serovar Typhi isolates from Bangladesh, Indonesia, Taiwan, and Vietnam. Antimicrob Agents Chemother. 2014;58: 6501–6507. pmid:25136011
  10. 10. Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill F-X, Goodhead I, et al. High-throughput sequencing provides insights into genome variation and evolution in Salmonella Typhi. Nat Genet. 2008;40: 987–993. pmid:18660809
  11. 11. Wong VK, Baker S, Connor TR, Pickard D, Page AJ, Dave J, et al. An extended genotyping framework for Salmonella enterica serovar Typhi, the cause of human typhoid. Nat Comms. 2016;7: 1–11. pmid:27703135
  12. 12. Wong VK, Baker S, Pickard DJ, Parkhill J, Page AJ, Feasey NA, et al. Phylogeographical analysis of the dominant multidrug-resistant H58 clade of Salmonella Typhi identifies inter- and intracontinental transmission events. Nat Genet. 2015;47: 632–639. pmid:25961941
  13. 13. International Typhoid Consortium, Wong VK, Holt KE, Okoro C, Baker S, Pickard DJ, et al. Molecular surveillance identifies multiple transmissions of Typhoid in West Africa. PLoS Negl Trop Dis. 2016;10: e0004781–22. pmid:27657909
  14. 14. Park SE, Pham DT, Boinett C, Wong VK, Pak GD, Panzner U, et al. The phylogeography and incidence of multi-drug resistant typhoid fever in sub-Saharan Africa. Nat Comms. 2018;9: 5094. pmid:30504848
  15. 15. Dyson ZA, Thanh DP, Bodhidatta L, Mason CJ, Srijan A, Rabaa MA, et al. Whole Genome sequence analysis of Salmonella Typhi isolated in Thailand before and after the introduction of a National Immunization Program. PLoS Negl Trop Dis. 2017;11: e0005274–15. pmid:28060810
  16. 16. Kariuki S, Revathi G, Kiiru J, Mengo DM, Mwituria J, Muyodi J, et al. Typhoid in Kenya is associated with a dominant multidrug-resistant Salmonella enterica serovar Typhi haplotype that is also widespread in Southeast Asia. J Clin Microbiol. 2010;48: 2171–2176. pmid:20392916
  17. 17. Phan MD, Kidgell C, Nair S, Holt KE, Turner AK, Hinds J, et al. Variation in Salmonella enterica serovar Typhi IncHI1 plasmids during the global spread of resistant Typhoid Fever. Antimicrob Agents Chemother. 2009;53: 716–727. pmid:19015365
  18. 18. Ashton PM, Nair S, Dallman T, Rubino S, Rabsch W, Mwaigwisya S, et al. MinION nanopore sequencing identifies the position and structure of a bacterial antibiotic resistance island. Nat Biotechnol. 2014;33: 296–300. pmid:25485618
  19. 19. Britto CD, Dyson ZA, Duchêne S, Carter MJ, Gurung M, Kelly DF, et al. Laboratory and molecular surveillance of paediatric typhoidal Salmonella in Nepal: Antimicrobial resistance and implications for vaccine policy. PLoS Negl Trop Dis. 2018;12: e0006408–19. pmid:29684021
  20. 20. Feasey NA, Gaskell K, Wong V, Msefula C, Selemani G, Kumwenda S, et al. Rapid emergence of multidrug resistant, H58-Lineage Salmonella Typhi in Blantyre, Malawi. PLoS Negl Trop Dis. 2015;9: e0003748–13. pmid:25909750
  21. 21. Hendriksen RS, Leekitcharoenphon P, Lukjancenko O, Lukwesa-Musyani C, Tambatamba B, Mwaba J, et al. Genomic signature of multidrug-resistant Salmonella enterica serovar Typhi isolates related to a massive outbreak in Zambia between 2010 and 2012. J Clin Microbiol. 2015;53: 262–272. pmid:25392358
  22. 22. Britto C, Pollard AJ, Voysey M, Blohmke CJ. An appraisal of the clinical features of pediatric enteric fever: Systematic review and meta-analysis of the age-stratified disease occurrence. Clin Infect Dis. 2017;64: 1604–1611. pmid:28369224
  23. 23. Cuypers WL, Jacobs J, Wong V, Klemm EJ, Deborggraeve S, Van Puyvelde S. Fluoroquinolone resistance in Salmonella: insights by whole-genome sequencing. Microb Genom. 2018;4: e000195. pmid:29975627
  24. 24. Klemm EJ, Shakoor S, Page AJ, Qamar FN, Judge K, Saeed DK, et al. Emergence of an extensively drug-resistant Salmonella enterica serovar Typhi clone harboring a promiscuous plasmid encoding resistance to Fluoroquinolones and Third-Generation Cephalosporins. mBio. 2018;9: e00105–18. pmid:29463654
  25. 25. Munir T, Lodhi M, Ansari JK, Andleeb S, Ahmed M. Extended spectrum Beta lactamase producing cephalosporin resistant Salmonella Typhi, reported from Rawalpindi Pakistan. J Pak Med Assoc. 2016;66: 1035–1036. pmid:27524545
  26. 26. Baker S, Thomson N, Weill F-X, Holt KE, kathryn . Genomic insights into the emergence and spread of antimicrobial-resistant bacterial pathogens. Science. 2018;360: 733–738. pmid:29773743
  27. 27. Bhan MK, Bahl R, Bhatnagar S. Typhoid and paratyphoid fever. The Lancet. 2005;366: 749–762.
  28. 28. Godbole GS, Day MR, Murthy S, Chattaway MA, Nair S. First report of CTX-M-15 Salmonella Typhi from England. Clin Infect Dis. 2018;66: 1976–1977. pmid:29471386
  29. 29. Ashton PM, Nair S, Peters TM, Bale JA, Powell DG, Painset A, et al. Identification of Salmonella for public health surveillance using whole genome sequencing. PeerJ. 2016;4: e1752–18. pmid:27069781
  30. 30. Parkhill J, Dougan G, James KD, Thomson NR, Pickard D, Wain J, et al. Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature. 2001;413: 848–852. pmid:11677608
  31. 31. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Meth. 2012;9: 357–359. pmid:22388286
  32. 32. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25: 2078–2079. pmid:19505943
  33. 33. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43: e15–e15. pmid:25414349
  34. 34. Day M, Doumith M, Jenkins C, Dallman TJ, Hopkins KL, Elson R, et al. Antimicrobial resistance in Shiga toxin-producing Escherichia coli serogroups O157 and O26 isolated from human cases of diarrhoeal disease in England, 2015. J Antimicrob Chemother. 2017;72: 145–152. pmid:27678285
  35. 35. Day MR, Doumith M, Do Nascimento V, Nair S, Ashton PM, Jenkins C, et al. Comparison of phenotypic and WGS-derived antimicrobial resistance profiles of Salmonella enterica serovars Typhi and Paratyphi. J Antimicrob Chemother. 2018;73: 365–372. pmid:29216342
  36. 36. McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother.; 2013;57: 3348–3357. pmid:23650175
  37. 37. Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, Villa L, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58: 3895–3903. pmid:24777092
  38. 38. Eaves DJ, Randall L, Gray DT, Buckley A, Woodward MJ, White AP, et al. Prevalence of mutations within the quinolone resistance-determining region of gyrA, gyrB, parC, and parE and association with antibiotic resistance in quinolone-resistant Salmonella enterica. Antimicrob Agents Chemother. 2004;48: 4012–4015. pmid:15388468
  39. 39. Hawkey J, Hamidian M, Wick RR, Edwards DJ, Billman-Jacobe H, Hall RM, et al. ISMapper: identifying transposase insertion sites in bacterial genomes from short read sequence data. BMC Genomics. 2015;16: 1–11.
  40. 40. Phan MD, Kidgell C, Nair S, Holt KE, Turner AK, Hinds J, et al. Variation in Salmonella enterica serovar Typhi IncHI1 plasmids during the global spread of resistant typhoid fever. Antimicrob Agents Chemother. 2009;53: 716–727. pmid:19015365
  41. 41. García-Fernández A, Villa L, Moodley A, Hasman H, Miriagou V, Guardabassi L, et al. Multilocus sequence typing of IncN plasmids. 2011;66: 1987–1991. pmid:21653604
  42. 42. Inouye M, Dashnow H, Raven L, Schultz MB, Pope BJ, Tomita T, et al. SRST2: Rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 2014;6: 1–16.
  43. 43. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13: e1005595–22. pmid:28594827
  44. 44. Wick RR, Schultz MB, Zobel J, Holt KE. Bandage: interactive visualization of de novo genome assemblies. Bioinformatics. 2015;31: 3550–3552.
  45. 45. Yu G, Smith DK, Zhu H, Guan Y, Lam TT-Y. ggtree: an rpackage for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2016;8: 28–36.
  46. 46. Wickham H. ggplot2: Elegant Graphics for Data Analysis. Cham: Springer International Publishing; 2016. https://doi.org/10.1007/978-3-319-24277-4
  47. 47. Holt KE, Thomson NR, Wain J, Phan MD, Nair S, Hasan R, et al. Multidrug-resistant Salmonella enterica serovar Paratyphi A harbors IncHI1 plasmids similar to those found in serovar Typhi. J Bacteriol. 2007;189: 4257–4264. pmid:17384186
  48. 48. Crump JA, Sjölund-Karlsson M, Gordon MA, Parry CM. Epidemiology, Clinical Presentation, laboratory diagnosis, antimicrobial resistance, and antimicrobial management of invasive Salmonella Infections. Clin Microbiol Rev. 2015;28: 901–937. pmid:26180063
  49. 49. The HC, Rabaa MA, Pham Thanh D, De Lappe N, Cormican M, Valcanis M, et al. South Asia as a reservoir for the global spread of ciprofloxacin-resistant Shigella sonnei: A cross-sectional study. PLoS Med. 2016;13: e1002055–12. pmid:27483136
  50. 50. Threlfall EJ, Ward LR, Skinner JA, Smith HR, Lacey S. Ciprofloxacin-resistant Salmonella Typhi and treatment failure. The Lancet. 1999;353: 1590–1591.
  51. 51. Baker S, Duy PT, Nga TVT, Dung TTN, Phat VV, Chau TT, et al. Fitness benefits in fluoroquinolone-resistant Salmonella Typhi in the absence of antimicrobial pressure. eLIFE. 2013;2: e01229. pmid:24327559
  52. 52. Djeghout B, Saha S, Sajib MSI, Tanmoy AM, Islam M, Kay GL, et al. Ceftriaxone-resistant Salmonella Typhi carries an IncI1-ST31 plasmid encoding CTX-M-15. J Med Microbiol. 2018;67: 620–627. pmid:29616895