A novel ciprofloxacin-resistant subclade of H58 Salmonella Typhi is associated with fluoroquinolone treatment failure

  1. Duy Pham Thanh
  2. Abhilasha Karkey
  3. Sabina Dongol
  4. Nhan Ho Thi
  5. Corinne N Thompson
  6. Maia A Rabaa
  7. Amit Arjyal
  8. Kathryn E Holt
  9. Vanessa Wong
  10. Nga Tran Vu Thieu
  11. Phat Voong Vinh
  12. Tuyen Ha Thanh
  13. Ashish Pradhan
  14. Saroj Kumar Shrestha
  15. Damoder Gajurel
  16. Derek Pickard
  17. Christopher M Parry
  18. Gordon Dougan
  19. Marcel Wolbers
  20. Christiane Dolecek
  21. Guy E Thwaites
  22. Buddha Basnyat
  23. Stephen Baker  Is a corresponding author
  1. Oxford University Clinical Research Unit, Vietnam
  2. Patan Academy of Health Sciences, Nepal
  3. Oxford University, United Kingdom
  4. The London School of Hygiene and Tropical Medicine, United Kingdom
  5. The University of Melbourne, Australia
  6. The Wellcome Trust Sanger Institute, United Kingdom
  7. Civil Services Hospital, Nepal
  8. Nagasaki University, Japan
  9. University of Oxford, United Kingdom

Abstract

The interplay between bacterial antimicrobial susceptibility, phylogenetics and patient outcome is poorly understood. During a typhoid clinical treatment trial in Nepal, we observed several treatment failures and isolated highly fluoroquinolone-resistant Salmonella Typhi (S. Typhi). Seventy-eight S. Typhi isolates were genome sequenced and clinical observations, treatment failures and fever clearance times (FCTs) were stratified by lineage. Most fluoroquinolone-resistant S. Typhi belonged to a specific H58 subclade. Treatment failure with S. Typhi-H58 was significantly less frequent with ceftriaxone (3/31; 9.7%) than gatifloxacin (15/34; 44.1%)(Hazard Ratio 0.19, p=0.002). Further, for gatifloxacin-treated patients, those infected with fluoroquinolone-resistant organisms had significantly higher median FCTs (8.2 days) than those infected with susceptible (2.96) or intermediately resistant organisms (4.01)(p<0.001). H58 is the dominant S. Typhi clade internationally, but there are no data regarding disease outcome with this organism. We report an emergent new subclade of S. Typhi-H58 that is associated with fluoroquinolone treatment failure.

Clinical trial registration: ISRCTN63006567.

https://doi.org/10.7554/eLife.14003.001

eLife digest

People who ingest a type of bacteria called Salmonella Typhi can develop the symptoms of typhoid fever. This disease is common in low-income settings in Asia and Africa, and causes a high rate of death in people who are not treated with antimicrobial drugs.

During a study in Nepal, Thanh et al. tried to evaluate which of two antimicrobials was better for treating typhoid fever. One of the drugs – called gatifloxacin – did not work in some of the patients. To understand why this treatment failed, Thanh et al. decoded the entire DNA sequences of all the Salmonella Typhi bacteria isolated during the study. Comparing this genetic data to the clinical data of the patients identified a new variant of Salmonella Typhi. These bacteria have a specific combination of genetic mutations that render them resistant to the family of drugs that gatifloxacin belongs to – the fluoroquinolones.

Patients infected with the variant bacteria and treated with gatifloxacin were highly likely to completely fail treatment and have longer-lasting fevers. On further investigation Thanh et al. found these organisms were likely recently introduced into Nepal from India.

Fluoroquinolones are amongst the most effective and common antimicrobials used to treat typhoid fever and other bacterial infections. However, the presence of bacteria that are resistant to these compounds in South Asia means that they should no longer be the first choice of drug to treat typhoid fever in this location.

https://doi.org/10.7554/eLife.14003.002

Introduction

Enteric (typhoid) fever, a systemic infection caused predominantly by the bacterium Salmonella enterica subspecies enterica serovar Typhi (S. Typhi), remains one of the principal bacterial causes of febrile disease in low-income countries (Parry et al., 2002). S. Typhi is a distinct, monophyletic lineage of S. enterica that is exquisitely adapted to cause disease only in humans (Roumagnac et al., 2006), characterised by a non-specific fever with malaise and asymptomatic convalescent carriage (Parry et al., 2002). There are an estimated 20–30 million new cases of enteric fever per year globally (Crump and Mintz, 2010), with the majority occurring in Asia, but there is an increasingly recognised burden of disease across sub-Saharan Africa.

Antimicrobial resistance is a major global health challenge, and resistance against the most commonly used antimicrobials for treating enteric fever has evolved successively over the last 30 years. Ampicillin, chloramphenicol and trimethoprim-sulfamethoxazole were originally standard-of-care for enteric fever. However, multidrug resistance (MDR) against these agents began to emerge in the 1970s and 1980s (Olarte and Galindo, 1973; Wain et al., 2003). Consequently, third-generation cephalosporins and fluoroquinolones (FQs) became the most clinically reliable drugs for treating enteric fever (Kariuki et al., 2015), and were formally advocated by the World Health Organization (WHO) in 2003 (World Health Organization, 2003). S. Typhi isolates with acquired resistance against third-generation cephalosporins are rare (Hendriksen et al., 2015), but S. Typhi exhibiting reduced susceptibility to FQs, induced by sequential mutations in the gene encoding a target protein (gyrA), now dominate internationally (Emary et al., 2012; Kariuki et al., 2010). The global ascendency of S. Typhi strains with reduced susceptibility to FQs has been partly catalysed by the dissemination of a specific MDR lineage (H58) across Asia and Africa (Wong et al., 2015). These H58 strains are rapidly displacing other lineages, and strains with gyrA mutations may have a fitness advantage, even in the absence of antimicrobial exposure (Baker et al., 2013).

We have previously shown that protracted fever clearance times (FCTs) are associated with organisms with higher Minimum Inhibitory Concentrations (MIC) against FQs in enteric fever patients treated with ciprofloxacin and ofloxacin (Parry et al., 2011). However, whilst the clinical efficacy of the older FQs in enteric fever is contentious, we have shown that the fourth-generation FQ, gatifloxacin, has remained efficacious for uncomplicated disease, even in patients infected with S. Typhi strains with reduced ciprofloxacin susceptibility (MIC ≥0.125 μg/mL) (Pandit et al., 2007; Koirala et al., 2013; Arjyal et al., 2011).

During a recent randomised controlled trial (RCT) comparing ceftriaxone and gatifloxacin, conducted in Nepal, we observed an increased number of treatment failures associated with FQ-resistant (ciprofloxacin MIC>32 μg/ml) S. Typhi, prompting the data safety and monitoring board to stop the trial (Arjyal et al., 2016). Aiming to assess the molecular epidemiology of the infecting isolates and investigate how genotype may be related to treatment outcome, we performed whole genome sequencing (WGS) on the S. Typhi isolated during this trial, and after stratifying by genotype, we assessed clinical presentation and outcome.

Results

Salmonella Typhi whole genome sequencing

We performed WGS on the 78 available S. Typhi isolates from patients in both RCT treatment arms (gatifloxacin and ceftriaxone) (Supplementary file 1). The resulting phylogeny, which incorporated reference sequence CT18, indicated that the majority of isolates (65/78; 83.3%) fell within the H58 lineage, while the remaining 13 (16.7%) represented eight different lineages (Figure 1). All but four of the H58 strains contained the common DNA gyrase (gyrA) mutation in codon 83 (S83F), which confers reduced susceptibility to FQs (ciprofloxacin MIC; 0.125–0 5 μg/ml) (Parry et al., 2010). Nested within the S83F H58 group, but separated from the rest of the group by a branch defined by 30 SNPs, was an H58 subclade comprised of 12 isolates containing the S83F gyrA mutation, a mutation in gyrA at codon 87 (D87N), and an additional mutation in the topoisomerase gene, parC (S80I) (H58 triple mutant). Notably, these H58 triple mutants shared high MICs against ciprofloxacin (≥24 μg/ml). Further, an additional two non-H58 RCT isolates with ciprofloxacin MIC≥24 μg/ml had the S83F gyrA mutation, an alternative mutation at codon 87 (D87V), the S80I parC mutation, and an A364V mutation in parE (Figure 1, Supplementary file 1). Notably, none of the sequenced isolates harboured plasmid-mediated quinolone resistance genes (PMQR) or contained additional antimicrobial resistance genes within the well-described S. Typhi-associated IncH1 family of plasmids.

The phylogenetic structure of 78 Nepali Salmonella Typhi isolated during a gatifloxacin versus ceftriaxone randomised controlled trial.

Maximum likelihood phylogeny based on core-genome SNPs of 78 Salmonella Typhi RCT isolates with the corresponding metadata, including the presence of mutations (dark grey) in gyrA (S83F, D87V and D87N), parC (S80I) and parC (A364V) and susceptibility to ciprofloxacin (susceptible, light blue; intermediate, mid-blue and non-susceptible, dark blue) by Minimum Inhibitory Concentration (MIC). The reference strain CT18 was used for context and highlighted by the black boxes. Red lines linking to metadata show isolates belonging to the Salmonella Typhi H58 lineage (with H58 triple mutants highlighted), other lineages (non-H58) are shown with black lines. The scale bar indicates the number of substitutions per variable site (see methods). Asterisks indicate ≥85% bootstrap support at nodes of interest.

https://doi.org/10.7554/eLife.14003.003

Clinical presentation of Salmonella Typhi infections

We stratified clinical data from the RCT by H58 status of the corresponding S. Typhi isolates (H58; N=65, non-H58; N=13) and compared baseline characteristics between these groups. We found no significant differences in demographics and no association between disease severity at presentation between those infected with an H58 S. Typhi isolate or a non-H58 isolate (Supplementary file 2A). Next, we compared the baseline characteristics of patients stratified by ciprofloxacin susceptibility (susceptible, intermediate and resistant), and found no differences in disease severity or demographics on presentation; the only exception being that FQ-resistant S. Typhi were more frequently isolated from adults (Supplementary file 2B). A significantly lower proportion of H58 S. Typhi (4/65; 6.2%) were susceptible to FQs compared to non-H58 isolates (6/13; 46%) (p=0.001) (Table 1) and, overall, H58 isolates had significantly higher (but not resistant) MICs against the majority of tested antimicrobials than non-H58 isolates (Table 1).

Table 1

Comparison of antimicrobial susceptibility by Salmonella Typhi lineage.

https://doi.org/10.7554/eLife.14003.004
E testNon-H58 (N=13)H58 (N=65)p value*
MIC50MIC90GM (range)MIC50MIC90GM (range)
Amoxicillin0.510.77 (0.38–38)0.75>2561.43 (0.38–>256)0.0412
Chloramphenicol342.7 (1.5–8)4125.7 (2–>256)0.0147
Ceftriaxone0.060.060.06 (0.05–0.13)0.090.190.11 (0.03–0.64)0.0004
Gatifloxacin0.130.250.06 (0.01–2)0.1320.21 (0.01–3)0.1197
Nalidixic acid>256>25621.6 (1–>256)>256>256346.8 (1–>256)0.0004
Ofloxacin0.250.750.24 (0.03–>32)0.5>32321.09 (0.03–>32)0.0240
Trimethoprim sulphate0.020.050.03 (0.02–0.05)0.050.320.09 (0.01–>32)0.0016
Ciprofloxacin0.130.750.11 (0.01–>32)0.38>320.80 (0.02–>32)0.0051
Ciprofloxacin susceptibility group0.0008#
- Susceptible6 (46.2%)4 (6.2%)
- Intermediate4 (30.8%)48 (73.8%)
- Resistant3 (23.1%)13 (20.0%)
  1. *Comparisons between Salmonella Typhi lineage for MICs and ciprofloxacin susceptibility groups were based on the Wilcoxon rank sum test and Fisher’s exact test. respectively.

  2. MIC: minimum inhibitory concentration, measured in µg/ml

  3. #p value for comparison of susceptible vs. intermediate/resistant combined between groups by Fisher’s exact test is 0.001.

  4. GM: geometric mean, the upper range of the values was determined by multiplying the MIC by 2 if the result was >X (for example, >256 = 256*2 = 512).

Treatment failure and fever clearance times

The primary endpoint of the RCT in which these data were generated was a composite for treatment failure (see method and previous publication) (Arjyal et al., 2016). Treatment failure with H58 S. Typhi was significantly less common in the ceftriaxone group (3/31; 9.7%) than the gatifloxacin group (15/34; 44.1%) (Hazard Ratio (HR) of time to failure 0.19, 95%CI 0.05–0.56, p=0.002) (Table 2). Conversely, there was no significant difference in treatment failure between those infected with non-H58 isolates treated with gatifloxacin (0/6; 0%) or ceftriaxone (2/7; 28.6%) (p=0.32). Similarly, time to fever clearance differed significantly between the two treatment groups in H58 infections, with median FCTs of 5.03 days (interquartile range (IQR): 3.18–7.21) in the gatifloxacin group and 3.07 days (IQR: 1.89–4.52) in the ceftriaxone group (p<0.0006). Again, this trend was not mirrored in the non-H58 S. Typhi infections, with FCTs of 2.87 (IQR: 2.08–3.7) and 3.12 (IQR: 2.2–4.12) days for gatifloxacin and ceftriaxone, respectively (p=0.61) (Table 3). Moreover, in the gatifloxacin arm, H58 S. Typhi tended to be associated with a higher risk of treatment failure (p=0.06) and a longer fever clearance time (p=0.013) (Figure 2, Table 2 and Supplementary file 2C).

Table 2

Summary of time to treatment failure by Salmonella Typhi lineage and ciprofloxacin susceptibility.

https://doi.org/10.7554/eLife.14003.005
Time to treatment failureGatifloxacin (events/N)Ceftriaxone (events/N)Hazard ratio of time to failure (95%CI); p valueHeterogeneity test
(p value)
H58*0.020
- H5815/343/310.19 (0.05, 0.56); p=0.002
- Non-H580/62/73.87 (0.31, 534.24); p=0.32
Ciprofloxacin susceptibility group0.08
- Susceptible0/51/52.40 (0.13, 350.21); p=0.57
- Intermediate7/252/270.27 (0.05, 0.99); p=0.049
- Resistant8/102/60.27 (0.05, 1.01); p=0.052
  1. *Likelihood ratio test p=0.06 and 0.40 for comparison of time to treatment failure between H58 vs. non-H58 groups in gatifloxacin arm only and in all patients, respectively

  2. Likelihood ratio test p=0.007 for comparison of time to treatment failure between MIC groups in gatifloxacin arm only

The association of Salmonella Typhi lineage and ciprofloxacin susceptibility with treatment failure and fever clearance time in patients randomised to gatifloxacin.

(A) Kaplan-Meier curve for time to treatment failure by H58 and non-H58 Salmonella Typhi. (B) Kaplan-Meier curve for time to treatment failure by Salmonella Typhi susceptibility group (susceptible, intermediate, resistant to ciprofloxacin). (C) Non-parametric maximum likelihood estimators for interval-censored fever clearance time (see methods) by H58 and non-H58 Salmonella Typhi. (D) Non-parametric maximum likelihood estimators for interval-censored fever clearance time by Salmonella Typhi susceptibility group (susceptible, intermediate, resistant to ciprofloxacin).

https://doi.org/10.7554/eLife.14003.006

As we identified two non-H58 isolates that were also FQ-resistant (Figure 1), we additionally stratified outcome for the gatifloxacin arm (N=40 patients) by FQ susceptibility of the infecting organism. Those infected with FQ-resistant S. Typhi failed gatifloxacin treatment more frequently (8/10; 80%) than those infected with an intermediately resistant organism (7/25; 28%) or a susceptible organism (0/5; 0%) (p=0.007) (Figure 2 and Table 2). Furthermore, in the gatifloxacin arm, those infected with FQ-resistant organisms had significantly higher median FCTs than those infected with S. Typhi with alternative FQ susceptibility profiles (median FCTs (days): susceptible, 2.96 (IQR: 2.13–3.85), intermediate, 4.01 (IQR: 2.76–5.37) and resistant 8.2 (IQR: 5.99–10.5), respectively [p<0.0001]) (Table 3 and Supplementary file 2D). Comparatively, the median FCT for those infected with an FQ-resistant organism but randomised to ceftriaxone was 3.83 days (IQR: 2.96–4.7) (p<0.0001 for the between-treatment comparison).

Table 3

Summary of fever clearance time by Salmonella Typhi lineage and ciprofloxacin susceptibility.

https://doi.org/10.7554/eLife.14003.007
Fever clearance timeGatifloxacin median (IQR) daysCeftriaxone median (IQR) daysAcceleration factor (95%CI); p valueHeterogeneity test
(p value)
H58¥0.07
- H585.03 (3.18, 7.21)3.07 (1.89 ,4.52)1.59 (1.22, 2.09); p=0.0006
- Non-H582.87 (2.08, 3.7)3.12 (2.2, 4.12)0.90 (0.59, 1.36); p=0.61
Ciprofloxacin susceptibility group0.015
- Susceptible2.96 (2.13, 3.85)4.78 (4.01, 5.5)0.71 (0.49, 1.02); p=0.07
- Intermediate4.01 (2.76, 5.37)2.63 (1.52, 4.05)1.31 (0.97, 1.76); p=0.07
- Resistant8.2 (5.99, 10.5)3.83 (2.96, 4.7)2.23 (1.57, 3.17); p<0.0001
  1. ¥p=0.013 and p=0.029 for comparison of interval censored time to fever clearance between H58 vs. non-H58 groups in gatifloxacin arm only and in all patients, respectively

  2. p<0.0001 for comparison of interval censored time to fever clearance between MIC groups in gatifloxacin arm only

The emergence of fluoroquinolone-resistant Salmonella Typhi

To measure the pattern of emergence of FQ-resistant S. Typhi in Nepal, we compiled FQ susceptibility data from 837 organisms isolated during enteric fever RCTs conducted at Patan Hospital between 2005 and 2014 (Figure 3) (Pandit et al., 2007; Koirala et al., 2013; Arjyal et al., 2011). MICs against FQs were generally higher for S. Paratyphi A than for S. Typhi. There was a significant temporal increase in S. Typhi MICs against both ciprofloxacin (p<0.0001) and gatifloxacin (p<0.0001), with a sharp increase from 2009. MICs against gatifloxacin in S. Paratyphi A also significantly increased with time (p<0.0001); however, MICs against ciprofloxacin showed only weak evidence of an upward trend over time (p=0.06).

Minimum Inhibitory Concentrations of Nepali Salmonella Typhi and Salmonella Paratyphi against ciprofloxacin and gatifloxacin over ten years.

Minimum Inhibitory Concentrations (μg/ml) for 568 Nepali Salmonella Typhi (blue) and 269 Nepali Salmonella Paratyphi A (grey) against (A) ciprofloxacin and (B) gatifloxacin collected from four randomised controlled trials conducted between 2005–2014 at Patan Hospital in Kathmandu, Nepal (Pandit et al., 2007; Koirala et al., 2013; Arjyal et al., 2011). The smoothed line derived from the generalized additive model showing a non-linear increase in Minimum Inhibitory Concentrations over time, with shading representing the 95% confidence interval. Lower and upper horizontal lines represent the current CLSI cut-offs for susceptible/intermediate and intermediate/resistant, respectively (CLSI, 2012).

https://doi.org/10.7554/eLife.14003.008

We hypothesised that the H58 triple mutants represented a contemporary importation into Nepal. To explore this, we compared the genomes of the 78 RCT S. Typhi isolates with those from 58 supplementary S. Typhi isolates from previous studies conducted between 2008 and 2013 in this setting (Figure 4, Supplementary file 1) (Wong et al., 2015). We found that the majority of the local H58 isolates (84/121; 69.4%) were closely related; these strains represented an 'endemic' Nepali H58 clade containing a single S83F gyrA mutation. Additionally, we identified a further five Nepali strains isolated in 2013 that belonged to the H58 triple mutant group, and had an MIC ≥24 μg/ml against ciprofloxacin. Incorporating additional genome sequences from a recent international study of the H58 lineage (Wong et al., 2015), we found that all the Nepali H58 triple mutants were very closely related (5 SNPs to nearest neighbour) to H58 triple mutants isolated previously in neighbouring India between 2008 and 2012 (Figure 4).

The phylogenetic structure of fluoroquinolone resistant Salmonella Typhi in a regional context.

Maximum likelihood phylogeny based on core-genome SNPs of 136 (78 from the RCT) Salmonella Typhi isolates from Nepal and neighbouring India (Supplementary file 1). Main tree shows the overall phylogenetic structure and the presence of specific combinations of mutations in gyrA (S83F, D87V and D87N), parC (S80I) and parE (A364V). The inset shows a magnified view of the fluoroquinolone-resistant Salmonella Typhi H58 triple mutants from Nepal and their close association with similarly fluoroquinolone-resistant Salmonella Typhi H58 triple mutants from India (Wong et al., 2015). The scale bar on the primary tree indicates the number of substitutions per variable site, while that in the inset indicates genetic distance in number of SNPs (see methods). Asterisks indicate ≥85% bootstrap support at nodes of interest.

https://doi.org/10.7554/eLife.14003.009

Discussion

Our study shows that a new FQ-resistant subclade of H58 S. Typhi has been introduced into Nepal and is associated with a lack of FQ efficacy. This subclade was associated with longer FCTs and treatment failure in patients treated with the FQ, gatifloxacin. For the first time, we can conclusively show how enteric fever patients respond to FQ treatment when infected with a specific subclade of H58, thereby linking organism genotype with a treatment phenotype. Given the international significance of FQs for the treatment of enteric fever and other bacterial infections, our findings have major global health implications for the long-term use and efficacy of this group of antimicrobials.

Our data suggest these FQ-resistant S. Typhi strains circulating in Nepal most likely descended from a single ancestor carrying the triple gyrA/parC mutant, such as that isolated in Nepal in 2011 (Koirala et al., 2012). This isolate was also associated with treatment failure, although this organism was not genome sequenced and was assumed to be an isolated case. More significantly, several very closely related strains were genome sequenced during an international study of H58 S. Typhi (Wong et al., 2015). These organisms had the same combination of triple FQ resistance mutations as those described here; our analysis shows they belong to the same subclade of H58. These strains had equivalently high MICs against ciprofloxacin and were isolated in India between 2008 and 2012. However, there were no associated patient outcome data for these strains and other reports from India have been limited. Our data implies that this lineage was introduced into Nepal from India or elsewhere in South Asia within the last 4–5 years and has subsequently entered in an endemic transmission cycle in Kathmandu. Given the large extent of human movement between India and Nepal, we propose this is the most likely route of introduction. However, there is also a small possibility that multiple strains independently gained resistance against FQs through the same selective pressure.

Appropriate antimicrobial therapy is critical in the treatment of enteric fever, as effective drugs curtail symptoms and prevent life threatening complications. Our data has substantial repercussions for enteric fever treatment, and we advocate that FQs should no longer be used for empirical enteric fever therapy on the Indian subcontinent, as we predict these strains are likely to be widespread and are associated with poor outcomes with FQ therapy. Notably, in the RCT from which these data were derived we used the newer generation FQ, gatifloxacin, which binds to a different location on the DNA gyrase than the older FQs and is not as susceptible to the common resistance mutations (Lu et al., 1999). The isolates in this study were not generally resistant to gatifloxacin according to the current CLSI guidelines for Enterobacteriaceae (CLSI, 2012); we suggest that these guidelines be modified specifically for S. Typhi to reflect these new clinical data. We additionally propose that S. Typhi genotyping, mapping and susceptibility testing is performed routinely and rapidly in reference laboratories inside and outside of South Asia to monitor the international spread of these strains and ensure the provision of alternative efficacious therapies to returning travellers (Lee et al., 2013; García-Fernández et al., 2015). In cases of infection with these FQ-resistant isolates, we suggest that ceftriaxone and azithromycin be used as alternatives, and do not currently recommend a return to the use of first-line drugs without contemporary data on treatment outcome. Whilst none of the isolates in this study were MDR, we predict a rapid return of MDR strains if there is a hasty return to older first-line alternatives.

This study has limitations. First, the clinical data was collected from one study in a single location, thus limiting utility outside this setting. Second, the overall sample size (and the gatifloxacin group subsampling) of those with culture-positive S. Typhi-associated enteric fever was relatively small, and the analysis presented here was performed in a post hoc manner. Notwithstanding these limitations, we were able to show a highly significant association between disease outcome and susceptibility profile of the infecting organism. Further, by using WGS, we were able to pinpoint causative mutations, identify the subclade responsible for treatment failure and relate these strains to other isolates circulating outside Nepal in other parts of South Asia. The methodologies presented here, in which clinical outcome data are combined with genome sequences and antimicrobial susceptibility data, should become the gold standard for informing empiric treatment for all invasive bacterial infections and understanding the role of bacterial genotype and resistance profile on disease outcome for other bacterial infections. No other combination of methodologies would provide the granularity of data required to understand the epidemiology and clinical impact of this emergent strain in detail.

In conclusion, our data, for the first time, show a significant association between S. Typhi genotype, antimicrobial susceptibility and disease outcome for those treated with gatifloxacin in a cohort of Nepali enteric fever patients. A FQ-resistant variant of Typhi H58 has emerged in Nepal and is associated with the clinical failure of FQs. Our data suggest these isolates are likely widespread in the subcontinent and FQs should not be recommended for empirical enteric fever therapy in this setting.

Materials and methods

Study design and setting

Request a detailed protocol

The RCT from which the organisms and corresponding clinical data originated for these analyses was conducted at Patan Hospital and the Civil Hospital in the Lalitpur area of Kathmandu, Nepal, between 2011 and 2014, as described previously (Arjyal et al., 2016). The trial was registered at www.clinicaltrials.gov (ISRCTN63006567). Briefly, patients were randomly assigned to seven days of treatment with either oral gatifloxacin (400 mg tablets, Square Pharmaceuticals Limited, Bangladesh) at a dose of 10 mg/kg once daily or intravenous ceftriaxone (Powercef, 1000mg injection vial, Wockhardt Ltd, India), injected over 10 min at a dose of 60 mg/kg up to a maximum of two grams (aged 2 to 13 years) or two grams (≥14 years) once daily.

A detailed description of the RCT from which these data were generated has been previously published (Arjyal et al., 2016). The primary endpoint was a composite of treatment failure, defined as the occurrence of at least one of the following events: fever clearance time (FCT) (time from the first dose of a study drug until the temperature dropped to ≤37·5°C and remained there for at least two days) more than seven days post-treatment initiation; requirement for rescue treatment as judged by the treating physician; blood culture positivity for S. Typhi or S. Paratyphi on day eight of treatment (microbiological failure); culture-confirmed or syndromic enteric fever relapse within 28 days of initiation of treatment; and the development of any enteric fever-related complication (e.g. clinically significant bleeding, fall in the Glasgow Coma Score, perforation of the gastrointestinal tract and hospital admission) within 28 days after the initiation of treatment. Time to treatment failure was defined as the time from the first dose of treatment until the date of the earliest failure event. FCTs were calculated electronically using twice-daily recorded temperatures and treated as interval-censored outcomes. Patients without fever clearance or relapse, respectively, were censored at the time of their last follow-up visit (additional details regarding study procedures can be found in Arjyal at al. 2016 (Arjyal et al., 2016)).

Blood (3 ml if aged <14 years; 8 ml if aged ≥14 years) was taken from all patients for bacterial culture on enrolment. Adult blood samples were inoculated into media containing tryptone soya broth and sodium polyanethol sulphonate, up to a total volume of 50 mL. Bactec Peds Plus culture bottles (Becton Dickinson, New Jersey, USA) were used for paediatric blood samples. Culture results were reported for up to seven days, positive bottles were subcultured onto blood, chocolate and MacConkey agar and presumptive Salmonella colonies were identified using standard biochemical tests and serotype-specific antisera (Murex Biotech, Dartford, England). Antimicrobial susceptibility testing was performed by the modified Bauer-Kirby disc diffusion method with zone size interpretation based on CLSI guidelines (CLSI, 2012). Etests were used to determine MICs, following the manufacturer's recommendations (bioMérieux, France). Ciprofloxacin MICs were used to categorise S. Typhi isolates as susceptible (≤0.06 μg/mL), intermediate (0.12–0.5 μg/mL) and resistant (≥1 μg/mL) following CLSI guidelines (CLSI, 2012).

Whole genome sequencing and analysis

Request a detailed protocol

Genomic DNA from Nepali S. Typhi organisms originating from this RCT (78 isolates) was extracted using the Wizard Genomic DNA Extraction Kit (Promega, Wisconsin, USA) (Supplementary file 1) (Karkey et al., 2013). Two μg of genomic DNA was subjected to WGS on an Illumina Miseq platform, following the manufacturer’s recommendations to generate 250bp/100bp paired-end reads. All reads were mapped to the reference sequence of S. Typhi CT18 (accession no: AL515582) using SMALT (version 0.7.4). Candidate single nucleotide polymorphisms (SNPs) were called against the reference sequence using SAMtools (Li et al., 2009) and filtered with a minimal phred quality of 30 and a quality cut-off of 0.75. The allele at each locus in each isolate was determined by reference to the consensus base in that genome, using samtools mpileup and removing low confidence alleles with consensus base quality ≤20, read depth ≤5 or a heterozygous base call. SNPs called in phage regions, repetitive sequences or recombinant regions were excluded, (Wong et al., 2015) resulting in a final set of 1,607 chromosomal SNPs. Strains belonging to haplotype H58 were defined by the SNP glpA-C1047T (position 2348902 in S. Typhi CT18, BiP33) (Emary et al., 2012; Holt et al., 2008; Parkhill et al., 2001).

A maximum likelihood (ML) phylogeny was estimated using a 1440 SNP alignment of the 78 RCT isolates in RAxML (version 7.8.6) with the generalized time-reversible substitution model (GTR) and a gamma distribution, with support for the phylogeny assessed via 1000 bootstrap replicates. The alignment was then compared to a global S. Typhi sequence database, with a particular focus on identifying sequences with a mutational profile suggestive of shared ancestry with a divergent H58 clade identified in the previous phylogeny. A secondary ML phylogenetic tree was then inferred from the SNP alignment of the 136 Nepali Typhi along with 19 recently described Typhi H58 with the aforementioned mutational profile, using the same parameters as above (1642 SNPs; Supplementary file 1) (Wong et al., 2015). Raw sequence data are available in the European Nucleotide Archive (ENA) (Supplementary file 1).

Statistical analysis

Request a detailed protocol

Comparison of baseline characteristics within patient groups, stratified by the H58 status or susceptibility category of their corresponding S. Typhi isolates was performed using the Kruskal Wallis test for continuous variables and Fisher’s exact test for categorical variables. Time to treatment failure was analysed using Firth’s penalized maximum likelihood bias reduction method for Cox regression as a solution for the non-convergence of likelihood function in the case of zero event counts in subgroups (Firth, 1993). For comparisons between treatment arms, H58 status, or ciprofloxacin susceptibility group, the model included treatment arm, H58 status, or susceptibility group as a single covariate. Confidence intervals (CI) and p-values were calculated by profile-penalized likelihood. FCT was analysed as an interval-censored outcome, i.e. as the time interval from the last febrile temperature assessment until the first afebrile assessment, using parametric Weibull accelerated failure time models (Kalbfleisch and Prentice, 2002). Median and inter-quartile range (IQR) FCT calculations for subgroups were based on models for each subgroup separately. Acceleration factors were based on models that included treatment arm as the only covariate. The non-parametric maximum likelihood estimator (NPMLE) was used to visualize the distribution of FCT between groups. Heterogeneity between subgroups was tested with models that included an interaction between treatment arm and the sub-grouping variable. To study the emergence of FQ resistance, data from previous enteric fever trials from 2005–2014 (Pandit et al., 2007; Koirala et al., 2013; Arjyal et al., 2011) was pooled and generalized additive models (GAM) were used to examine potential non-linear trends of ciprofloxacin and gatifloxacin MICs over time. All analyses were performed using R software version 3.2.2 (Team, 2012).

Data availability

The following data sets were generated

References

    1. CLSI
    (2012)
    Clinical and laboratory standards institute: performance standards for antimicrobial susceptibility testing
    Twentieth Informational Supplement pp. M100–S22.
    1. Lu T
    2. Zhao X
    3. Drlica K
    (1999)
    Gatifloxacin activity against quinolone-resistant gyrase: allele-specific enhancement of bacteriostatic and bactericidal activities by the C-8-methoxy group
    Antimicrobial Agents and Chemotherapy 43:2969–2974.
    1. Team R
    (2012) R Foundation for Statistical Computing
    R: a language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, Open access available http//cran r-project org.

Article and author information

Author details

  1. Duy Pham Thanh

    The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    Contribution
    DPT, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  2. Abhilasha Karkey

    Oxford University Clinical Research Unit, Patan Academy of Health Sciences, Kathmandu, Nepal
    Contribution
    AK, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  3. Sabina Dongol

    Oxford University Clinical Research Unit, Patan Academy of Health Sciences, Kathmandu, Nepal
    Contribution
    SD, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  4. Nhan Ho Thi

    The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    Contribution
    NHT, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  5. Corinne N Thompson

    1. The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    2. Centre for Tropical Medicine and Global Health, Oxford University, Oxford, United Kingdom
    3. The London School of Hygiene and Tropical Medicine, London, United Kingdom
    Contribution
    CNT, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  6. Maia A Rabaa

    1. The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    2. Centre for Tropical Medicine and Global Health, Oxford University, Oxford, United Kingdom
    Contribution
    MAR, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0529-2228
  7. Amit Arjyal

    Oxford University Clinical Research Unit, Patan Academy of Health Sciences, Kathmandu, Nepal
    Contribution
    AA, Conception and design, Acquisition of data
    Competing interests
    The authors declare that no competing interests exist.
  8. Kathryn E Holt

    Department of Biochemistry and Molecular Biology, Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Melbourne, Australia
    Contribution
    KEH, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  9. Vanessa Wong

    The Wellcome Trust Sanger Institute, Cambridge, United Kingdom
    Contribution
    VW, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  10. Nga Tran Vu Thieu

    The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    Contribution
    NTVT, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  11. Phat Voong Vinh

    The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    Contribution
    PVV, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  12. Tuyen Ha Thanh

    The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    Contribution
    THT, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  13. Ashish Pradhan

    Civil Services Hospital, Kathmandu, Nepal
    Contribution
    AP, Recruited patients into study and collected clinical data, Acquisition of data
    Competing interests
    The authors declare that no competing interests exist.
  14. Saroj Kumar Shrestha

    Civil Services Hospital, Kathmandu, Nepal
    Contribution
    SKS, Recruited patients into study and collected clinical data, Acquisition of data
    Competing interests
    The authors declare that no competing interests exist.
  15. Damoder Gajurel

    Civil Services Hospital, Kathmandu, Nepal
    Contribution
    DG, Recruited patients into study and collected clinical data, Acquisition of data
    Competing interests
    The authors declare that no competing interests exist.
  16. Derek Pickard

    The Wellcome Trust Sanger Institute, Cambridge, United Kingdom
    Contribution
    DP, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  17. Christopher M Parry

    1. The London School of Hygiene and Tropical Medicine, London, United Kingdom
    2. School of Tropical Medicine and Global Health, Department of Clinical Research, Nagasaki University, Nagasaki, Japan
    Contribution
    CMP, Conception and design, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  18. Gordon Dougan

    The Wellcome Trust Sanger Institute, Cambridge, United Kingdom
    Contribution
    GD, Conception and design, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  19. Marcel Wolbers

    1. The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    2. Centre for Tropical Medicine and Global Health, Oxford University, Oxford, United Kingdom
    Contribution
    MW, Analysis and interpretation of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  20. Christiane Dolecek

    1. The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    2. Centre for Tropical Medicine and Global Health, Oxford University, Oxford, United Kingdom
    3. Department of Zoology, University of Oxford, Oxford, United Kingdom
    Contribution
    CD, Conception and design, Acquisition of data, Contributed unpublished essential data or reagents
    Competing interests
    The authors declare that no competing interests exist.
  21. Guy E Thwaites

    1. The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    2. Centre for Tropical Medicine and Global Health, Oxford University, Oxford, United Kingdom
    Contribution
    GET, Conception and design, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  22. Buddha Basnyat

    Oxford University Clinical Research Unit, Patan Academy of Health Sciences, Kathmandu, Nepal
    Contribution
    BB, Conception and design, Acquisition of data
    Competing interests
    The authors declare that no competing interests exist.
  23. Stephen Baker

    1. The Hospital for Tropical Diseases, Wellcome Trust Major Overseas Programme, Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam
    2. Centre for Tropical Medicine and Global Health, Oxford University, Oxford, United Kingdom
    3. The London School of Hygiene and Tropical Medicine, London, United Kingdom
    Contribution
    SB, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    sbaker@oucru.org
    Competing interests
    The authors declare that no competing interests exist.

Funding

Wellcome Trust (100087/Z/12/Z)

  • Stephen Baker

Royal Society (100087/Z/12/Z)

  • Stephen Baker

National Health and Medical Research Council (1061409)

  • Kathryn E Holt

Li Ka Shing Foundation

  • Christiane Dolecek

Wellcome Trust (106158/Z/14/Z)

  • Kathryn E Holt
  • Gordon Dougan
  • Christiane Dolecek
  • Buddha Basnyat
  • Stephen Baker

The Oak Foundation (OCAY-15-547)

  • Duy Pham Thanh
  • Abhilasha Karkey
  • Stephen Baker

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This project was funded by the Wellcome Trust of Great Britain (106158/Z/14/Z). SB is a Sir Henry Dale Fellow, jointly funded by the Wellcome Trust and the Royal Society (100087/Z/12/Z). KEH is supported by fellowship #1061409 from the NHMRC of Australia. CD was funded by the Li Ka Shing Foundation Global Health Programme at the University of Oxford. DTP and AK are funded as leadership fellows through the Oak Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We wish to acknowledge Mangal Rawal, Sumi Munankarmi, Bibek Karki, Radheshyam KC, Sudeep Dhoj Thapa, Rabin Gautam, Priyanka Tiwari, Manisha Risal, Surendra Shrestha, Balmukunda Neupane, Nabraj Regmi, Krishna Prajapati, Bimal Thapa and the trial monitors Nguyen Thi Phuong Dung and Nguyen Thi Thanh Thuy for their assistance in conducting this trial.

Ethics

Clinical trial Registry: ISRCTN.

Registration ID: ISRCTN63006567.

Human subjects: This study was performed following the principles of the declaration of Helsinki. Written informed consent to participate in all studies from Nepal contributing data for this analysis was required from all patients. For those aged <18 years, written informed consent was obtained from a parent or an adult guardian. The protocol was reviewed and approved by the Ethics Committee of the Nepal Health Research Council (NHRC) and the Oxford Tropical Research Ethics Committee (OxTREC) UK.

Version history

  1. Received: December 22, 2015
  2. Accepted: February 19, 2016
  3. Version of Record published: March 14, 2016 (version 1)

Copyright

© 2016, Pham Thanh et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 2,418
    Page views
  • 455
    Downloads
  • 79
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Duy Pham Thanh
  2. Abhilasha Karkey
  3. Sabina Dongol
  4. Nhan Ho Thi
  5. Corinne N Thompson
  6. Maia A Rabaa
  7. Amit Arjyal
  8. Kathryn E Holt
  9. Vanessa Wong
  10. Nga Tran Vu Thieu
  11. Phat Voong Vinh
  12. Tuyen Ha Thanh
  13. Ashish Pradhan
  14. Saroj Kumar Shrestha
  15. Damoder Gajurel
  16. Derek Pickard
  17. Christopher M Parry
  18. Gordon Dougan
  19. Marcel Wolbers
  20. Christiane Dolecek
  21. Guy E Thwaites
  22. Buddha Basnyat
  23. Stephen Baker
(2016)
A novel ciprofloxacin-resistant subclade of H58 Salmonella Typhi is associated with fluoroquinolone treatment failure
eLife 5:e14003.
https://doi.org/10.7554/eLife.14003

Share this article

https://doi.org/10.7554/eLife.14003

Further reading

    1. Ecology
    2. Epidemiology and Global Health
    Aleksandra Kovacevic, David RM Smith ... Lulla Opatowski
    Research Article

    Non-pharmaceutical interventions implemented to block SARS-CoV-2 transmission in early 2020 led to global reductions in the incidence of invasive pneumococcal disease (IPD). By contrast, most European countries reported an increase in antibiotic resistance among invasive Streptococcus pneumoniae isolates from 2019 to 2020, while an increasing number of studies reported stable pneumococcal carriage prevalence over the same period. To disentangle the impacts of the COVID-19 pandemic on pneumococcal epidemiology in the community setting, we propose a mathematical model formalizing simultaneous transmission of SARS-CoV-2 and antibiotic-sensitive and -resistant strains of S. pneumoniae. To test hypotheses underlying these trends five mechanisms were built into the model and examined: (1) a population-wide reduction of antibiotic prescriptions in the community, (2) lockdown effect on pneumococcal transmission, (3) a reduced risk of developing an IPD due to the absence of common respiratory viruses, (4) community azithromycin use in COVID-19 infected individuals, (5) and a longer carriage duration of antibiotic-resistant pneumococcal strains. Among 31 possible pandemic scenarios involving mechanisms individually or in combination, model simulations surprisingly identified only two scenarios that reproduced the reported trends in the general population. They included factors (1), (3), and (4). These scenarios replicated a nearly 50% reduction in annual IPD, and an increase in antibiotic resistance from 20% to 22%, all while maintaining a relatively stable pneumococcal carriage. Exploring further, higher SARS-CoV-2 R0 values and synergistic within-host virus-bacteria interaction mechanisms could have additionally contributed to the observed antibiotic resistance increase. Our work demonstrates the utility of the mathematical modeling approach in unraveling the complex effects of the COVID-19 pandemic responses on AMR dynamics.

    1. Epidemiology and Global Health
    Olivera Djuric, Elisabetta Larosa ... The Reggio Emilia Covid-19 Working Group
    Research Article

    Background:

    The aim of our study was to test the hypothesis that the community contact tracing strategy of testing contacts in households immediately instead of at the end of quarantine had an impact on the transmission of SARS-CoV-2 in schools in Reggio Emilia Province.

    Methods:

    We analysed surveillance data on notification of COVID-19 cases in schools between 1 September 2020 and 4 April 2021. We have applied a mediation analysis that allows for interaction between the intervention (before/after period) and the mediator.

    Results:

    Median tracing delay decreased from 7 to 3.1 days and the percentage of the known infection source increased from 34–54.8% (incident rate ratio-IRR 1.61 1.40–1.86). Implementation of prompt contact tracing was associated with a 10% decrease in the number of secondary cases (excess relative risk –0.1 95% CI –0.35–0.15). Knowing the source of infection of the index case led to a decrease in secondary transmission (IRR 0.75 95% CI 0.63–0.91) while the decrease in tracing delay was associated with decreased risk of secondary cases (1/IRR 0.97 95% CI 0.94–1.01 per one day of delay). The direct effect of the intervention accounted for the 29% decrease in the number of secondary cases (excess relative risk –0.29 95%–0.61 to 0.03).

    Conclusions:

    Prompt contact testing in the community reduces the time of contact tracing and increases the ability to identify the source of infection in school outbreaks. Although there are strong reasons for thinking it is a causal link, observed differences can be also due to differences in the force of infection and to other control measures put in place.

    Funding:

    This project was carried out with the technical and financial support of the Italian Ministry of Health – CCM 2020 and Ricerca Corrente Annual Program 2023.