key: cord-0012944-9zmnbjv7
authors: Guglielmino, Christine J. D.; Kakkanat, Asha; Forde, Brian M.; Rubenach, Sally; Merone, Lea; Stafford, Russell; Graham, Rikki M. A.; Beatson, Scott A.; Jennison, Amy V.
title: Outbreak of multi-drug-resistant (MDR) Shigella flexneri in northern Australia due to an endemic regional clone acquiring an IncFII plasmid
date: 2020-09-04
journal: Eur J Clin Microbiol Infect Dis
DOI: 10.1007/s10096-020-04029-w
sha: 74cce92c5b59b36bfeef98709a23fdeec6e94eeb
doc_id: 12944
cord_uid: 9zmnbjv7

Epidemiological surveillance of Shigella spp. in Australia is conducted to inform public health response. Multi-drug resistance has recently emerged as a contributing factor to sustained local transmission of Shigella spp. All data were collected as part of routine public health surveillance, and strains were whole-genome sequenced for further molecular characterisation. 108 patients with an endemic regional Shigella flexneri strain were identified between 2016 and 2019. The S. flexneri phylogroup 3 strain endemic to northern Australia acquired a multi-drug resistance conferring bla(DHA) plasmid, which has an IncFII plasmid backbone with virulence and resistance elements typically found in IncR plasmids. This is the first report of multi-drug resistance in Shigella sp. in Australia that is not associated with men who have sex with men. This strain caused an outbreak of multi-drug-resistant S. flexneri in northern Australia that disproportionality affects Aboriginal and Torres Strait Islander children. Community controlled public health action is recommended. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s10096-020-04029-w) contains supplementary material, which is available to authorized users.

Shigellosis in Australia is generally understood to occur as sporadic cases among returning travellers or amongst men who have sex with men (MSM) as locally acquired, sexually transmitted infections [1] . However, in northern Queensland and the Northern Territory, locally acquired strains of both Shigella sonnei and S. flexneri are persistently transmitted within other population groups [2] . These endemic strains have previously been mostly susceptible to commonly used antimicrobial agents including ampicillin, ciprofloxacin and co-trimoxazole [3] .

In 2019, both Australia and Queensland had a shigellosis notification rate of 12 per 100,000 head of population, while in contrast, the Northern Territory had a notification rate of 119 per 100,000 [4] , and in the northernmost health and hospital service districts of Queensland (Torres and Cape, Cairns and Hinterland, North West, and Townsville [5] ), the notification rate was 24.7 per 100,000. Northern Queensland and the Northern Territory cover an area of approximately 1.9 million km 2 , much of which is sparsely populated and geographically remote. The population of northern Australia has a large proportion of Aboriginal and Torres Strait Islander peoples who are uniquely vulnerable to diarrhoeal disease for a variety of reasons including significant social disadvantage and inadequate housing conditions [6, 7] . Children of remote communities in this region experience a disproportionate infectious disease burden when compared with the entire Australian population [8] , including for trachoma which, although is understood to be eliminated in Queensland, is still endemic in the Northern Territory and can also be managed with azithromycin [9, 10] .

The emergence of MDR strains in MSM was a driver for guideline changes in 2019; subsequently, Australian therapeutic guidelines now recommend treating shigellosis with ceftriaxone in immunocompromised patients or those with severe infection, or where required for those who are high risk, whilst laboratory phenotypic AMR testing is pending [11] . Prior to this recent change, the recommended agents were cotrimoxazole, ciprofloxacin or azithromycin. Notification to CARAlert, the national alert system for critical antimicrobial resistance (AMR), is currently required if phenotypic resistance is observed to any three of the following: ampicillin/ amoxicillin, ciprofloxacin/norfloxacin, co-trimoxazole, ceftriaxone/cefotaxime/ceftazidime (3GCR) or azithromycin with molecular methods for characterisation of the beta-lactamase gene/s as either AmpC-like (class C plasmid-mediated resistance) or ESBL (extended-spectrum beta-lactamase) is recommended but not required for notification. Prior to 2019, azithromycin was not one of the categories listed in the laboratory handbook [12] .

Plasmids play a pivotal role in the acquisition and dissemination AMR via mobile genetic elements (MGE), with IncF plasmids being predominant types among Enterobacterales [13] . A study into the global dissemination of MDR S. flexneri among MSM found that the IncFII 35 plasmid pKSR100 is associated with resistant lineages [8] . This plasmid carries mph(A) and ermB conferring resistance to azithromycin, as well as other AMR determinants. Related plasmids have subsequently been found and are mainly associated with MSM in other MDR lineages globally [14] , including IncFII 35 pKSR100-like plasmid described in MSM associated S. flexneri 2a in metropolitan south-east Australia [1] . The IncFII 2 R100 plasmid from S. flexneri from Japan in the 1950s is considered the classical IncFII plasmid [13] and contains tet(A) and catA AMR genes.

The association between virulence and AMR genes facilitates the persistence of pathogenic bacteria leading to the development of high-risk clones [15] . The phage shock protein (PSP) system mediated by the psp operon is an integral part of the stress response, enabling bacteria to compete for survival under energy-limited conditions. The IncR pKSP30 from Klebsiella pneumoniae carries the sap operon (plays a role in resistance to cationic peptides), psp operon and 3'CS extreme of a class 1 integron consisting of qacEΔ1 conferring resistance to quaternary ammonium compounds, followed by sul1 along with other MDR determinants [16, 17] . The presence of atypical class 1 integrons is commonly reported in Shigella spp. from Asia and Africa [18] [19] [20] .

Whole-genome sequencing (WGS) is established as the method of choice for identification and typing of Shigella spp. [21] . A species defining study has divided clonal complex 245 (CC245) S. flexneri (all S. flexneri serogroups other than O6) into seven phylogroups (PGs) by core single nucleotide polymorphism (SNP) typing analysis [22] . Genes responsible for determining S. flexneri serotypes are encoded on horizontally transmissible elements and can easily be lost or gained under selective pressure; thus, serotype results may not represent true phylogenetic relationships [23] . Queensland Health Forensic and Scientific Services Public Health Microbiology is the state reference laboratory for Shigella spp. isolated in Queensland, undertaking routine WGS surveillance of all Shigella spp. isolates received, which consists of approximately 40% of total Shigellosis notifications in the state. The remainder of notifications is PCR-only diagnoses, with no isolate available for epidemiological typing. This routine surveillance has identified the persistence of an endemic local strain of PG3 S. flexneri in northern Queensland and the Northern Territory among 108 different patients between 2016 and 2019.

All data was collected as part of routine health surveillance activities, outbreak support and reporting between 2016 and 2019 under the Queensland Public Health Act [24] , (supplementary Table 1 ). Rates are calculated as crude rates using population data sourced from the Australian Bureau of Statistics [25].

All isolates were serotyped with agglutinating sera (Denka Seiken, Japan) and, where antimicrobial susceptibility results were available, performed with Vitek 2 GN AST card and Etests for azithromycin (bioMérieux, France) using EUCAST [26] interpretations.

WGS was performed using the Nextera XT library preparation kit and the Nextseq 500 platform (Illumina, CA) and reads trimmed with trimmomatic [27] and assembled with Spades [28] . Core SNP analysis was performed using snippy [29] with S. flexneri serotype 2a PG2 strain 301 chromosome NC_004337.2 as a reference. The presence of AMR determinants was determined with resfinder and point finder [30] . One strain (M2901) was selected for promethION sequencing (Oxford Nanopore Technologies, UK). Nanopore sequencing libraries were prepared using the rapid barcode sequencing kit SQK-RBK004 and sequenced on a single flow cell (R9.4.1) for 48 h. Raw Nanopore sequence data was de novo assembled using Flye with default parameters [31] . The draft genome was polished with Medaka [32] , followed by polishing with both long and short reads using NextPolish v1.2.4 [33] . Annotation was performed with NCBI prokaryotic genome annotation pipeline [34] .

The rate of shigellosis notifications for Aboriginal and Torres Strait Island people living in north Queensland has been consistently higher than observed for non-Indigenous residents of this region. This is most evident in the three years from 2017 to 2019, with 76 notifications per 100,000 population in Aboriginal and Torres Strait Islander residents compared with 11 notifications per 100,000 in non-Indigenous residents, and even more striking for S. flexneri (Fig. 1) .

Of the 108 patients with strains belonging to the endemic northern Australia PG3 S. flexneri clone, 103 (95%) resided in either the Northern Territory or one of the four northernmost health and hospital service districts of Queensland. The median age of patients was 14 years which is significantly lower than the median age of 33 years for all Shigellosis notifications for the outbreak periods 2016 to 2019 from Queensland (Kruskal-Wallis H = 42.0; P < 0.001).

Although most strains in the outbreak were serotype 2b, five were var X and four were 2a; however, all belonged to the same phylogroup (PG3) and clustered together by core SNP typing analysis. Figure 2 shows a phylogenetic tree built using SNP typing with the outbreak strains, compared with a publicly available PG3 S. flexneri 2a MDR strain AUSMDU00008332 associated with MSM in south-east Australia and globally, as well as the PG1 S. flexneri 1c strain AUSMDU00008355. The pKSR100like plasmid from strain AUSMDU00008332 does not possess drfA17, sul1 and aadA5; however, the one from AUSMDU00008355 does. Plasmid MLST [38] performed in silico via pubMLST revealed all 108 strains in the outbreak to have an IncFII 27 virulence plasmid carrying the invasion plasmid antigen IpaH. In this study, the outbreak strain appears to be established and endemic in the northern Australia region since 2016 but has acquired MDR determinants by gaining an additional bla DHA carrying IncFII 2 plasmid (pM2901) more recently in 2018. Prior to this, plasmid-based MGE associated with AMR has appeared in this clone; however, none has established themselves like the bla DHA plasmid. All strains were found to carry bla OXA-1 , tet(B), catA1 and dfrA1, plus at least one gene conferring aminoglycoside resistance while they expressed phenotypic resistance to ampicillin and/or amoxicillin. In addition to pM2901, some strains also carried an additional IncI1 14 plasmid carrying bla TEM1-C and sul2. All AMR conferring genes cited in this manuscript are listed with the class of AMR they are expected to confer resistance to in Table 1 .

Just two strains in the outbreak with a qnrS gene expressed phenotypic quinolone resistance. Further 27 strains carried the plasmid-mediated quinolone resistance gene qnrB [39] ; however, all 27 strains were classified as susceptible to quinolone antibiotics. No known quinolone resistance determining mutations were observed in any strain for gyrA or parC. Although all 108 strains had at least one trimethoprim dfrA resistance marker, with 26 of these also possessing sul1, only those strains also possessing the sul2 gene (n = 9) demonstrated co-trimoxazole resistance.

Many strains with multiple AMR markers did not meet the criteria for CARAlert notification owing to some not expressing phenotypic resistance to match their genotype. Moreover, azithromycin resistance is an additional test to those performed routinely, and if it is not tested for a Shigella sp., which is resistant to two other classes of antimicrobials, it will therefore not meet the criteria for notification even if it possesses a gene encoding for macrolide resistance. All four of the strains that showed phenotypic resistance to azithromycin when tested demonstrated the mph(A) macrolide resistance conferring gene. Although all but one outbreak strains had bla OXA-1 and were phenotypically resistant to ampicillin Fig. 2 Phylogenetic tree of core SNP typing built using snippy [29] and RAxML [35] , visualised with iTOL [36] , showing PG3 outbreak strains w i t h a c q u i r e d A M R g e n e s d e t e c t e d , i n c l u d i n g s t r a i n AUSMDU00008332 accession NZ_LR213455.1 from the global MDR MSM outbreak of PG3 S. flexneri 2a, and AUSMDU00008355 accession NZ_LR213452.1 from PG1 S. flexneri 1c for context. Strain M2901, for which the plasmid has been characterised in Fig. 3 , is marked. The year of isolation is noted. and/or amoxicillin, only those with an additional bla DHA/TEM/ CTX-M gene showed resistance to a third-generation cephalosporin (3GCR).

The pM2901 is an 81,732-bp plasmid with an average G + C content of 52.5% with transfer-encoding region. It has a backbone akin to pR100, but with an insertion of a psp operon and bla DHA , qacEΔ1, sul1 and qnrB genes similar to that of the IncR plasmid pKPS30 which carries bla OXA-1 , in addition to dfrA17 like the pKSR100-like plasmids found in non-PG3 Shigella spp. in MSM in south-east Australia. Unlike the pKSR100-like plasmids, aadA5 of pM2901 was truncated with a hypothetical protein. Additionally, the tet(A)/catA locus of pR100 was entirely absent from pM2901. The pM2901 plasmid is most similar to p7102_58-6, which has been described in Salmonella Worthington in vivo acquisition from Citrobacter amalonaticus [40] , and also pUB_DHA-1 from thirdgeneration cephalosporin-resistant (3GCR) Escherichia coli in England [41] . Figure 3 shows an alignment of the pM2901 against other plasmids. Sequence data for strain M2910 is available at the National Center for Biotechnology Information (NCBI) BioProject PRJNA599490, and assembly data at accessions CP058589-CP058593.

This is the first report of locally acquired MDR or azithromycin-resistant Shigella spp. in Australia that is not associated with MSM; this endemic strain disproportionately affects Aboriginal and Torres Strait Islander communities of remote Northern Territory and far-north Queensland.

The broad host range of the main plasmid involved in this outbreak, even though restricted to Enterobacterales [42] , heightens its public health impact having previously been reported in Salmonella sp. and E. coli., bla DHA-1 has traditionally been found on narrow-host range IncFII plasmids; however, Mata et al. [43] have reported the co-localisation of qnrB with bla DHA-1 on broad host range L/M plasmids and an IncN plasmid. This is the first report of the co-localisation of these genes in an IncFII plasmid in Shigella sp. Wang et al. recently reported a bla IMP plasmid with a psp operon clustered with bla DHA and qnrB in an Enterobacter cloacae IncHI2 plasmid [44] . pK245 (NCBI accession NC_010886) from Klebsiella pneumoniae is considered the protype plasmid of IncR [45] and does not contain the psp operon. Although IncR plasmids are thought to be non-transferable and non-mobilisable due to Fig. 3 Outbreak associated plasmid pM2901 aligned using BLASTn, as implemented in Easyfig [37] , against (A) pR100 IncFII 2 accession AP000342 and pKSP30 IncR accession KF793937. (B) pKSR100 IncFII 35 accession LN624486 and pKSR100-like IncFII 35 accession NZ_LR213453.1 from strain AUSMDU00008355 lacking a transfer system and a relaxase [45] , this demonstrates the transposition/plasmid recombination of IncR resistance and virulence elements to an IncF plasmid (containing a transfer system and relaxase), demonstrating the high plasticity of bacterial plasmids as anticipated by Compain et al. [46] . pHN7A8 from China JN232517 is another contemporary IncFII plasmid with a resistance region derived from R100 [47] . The dual presence of IncFII 2 and IncR in Shigella sp. strains in this study demonstrates the potential for strains to carry multiple plasmid types with different incompatibility types, which was historically thought to not be possible; however, IncFII plasmids have been shown previously to diverge and overcome the incompatibility barrier with other plasmids [48] . This study did not examine the host microbiomes of patients to determine if the plasmid is hosted in organisms other than Shigella spp., limiting the insight into highresolution plasmid relationships.

This outbreak consisted of multiple serotypes, which highlights a challenge for vaccine development where immunity to Shigella spp., is considered to be serotype-specific [49, 50] . There was no evidence in this outbreak of patients being reinfected with either the same or different serotypes or species of Shigella spp. The inclusion of different serotypes in the single cluster confirms phylogrouping as a more relevant method for public health surveillance and outbreak investigation than serotyping, as these cases would have been excluded from the investigation based on serotype results if WGS were not performed. We maintain the position that molecular genotyping of AMR is more suitable for public health surveillance purposes, while phenotyping is relevant for individual patient management as AMR genes may be present but not be expressed until selective pressure is induced by exposure to antimicrobial agents.

Without public health intervention, continued transmission of this now MDR strain of S. flexneri is likely to persist in this vulnerable population and has the potential to spread to other geographical regions. The emergence of azithromycin resistance markers in a population potentially treated with azithromycin for trachoma and/or other conditions supports the theory of Baker et al. [14] that selective pressure from azithromycin administration has contributed to the emergence of resistance in this organism. This theory is supported by the co-emergence of low-level azithromycin resistance in Neisseria gonorrhoeae (data not shown) and macrolide resistance genetic markers in Mycoplasma genitalium [51] that are currently occurring via different mechanisms in this same geographical area, and which is likely to also emerge in other infectious agents with continued social disadvantage. The recent inclusion of azithromycin as a separate category of AMR class contributing towards CARAlert notification is likely to facilitate more notifications than the previous criteria and therefore improve public health outcomes. This was evidenced by a 218% increase is CARAlert notifications for MDR Shigella spp. in Australia in 2019 from 2018, although this is primarily attributed to MDR outbreaks among MSM in other states [52] . Azithromycin testing is now required much more frequently than it was previously, for Shigella spp. from Australian patients in both urban and remote regions.

This report highlights the need for continuing AMR notification and surveillance which includes azithromycin for various organisms of public health significance among populations potentially subjected to community azithromycin distribution as a control measure for trachoma [53] . Addressing social inequalities is however the most fundamental aspect of reducing disease burden. Community controlled responses including health promotion initiatives to reduce the burden of diarrhoeal diseases, as well as rapid progression of healthy housing initiatives to reduce overcrowding and sanitation in remote communities are recommended. This is vital in improving housing, sanitation and socio-economic status and in continuing the important work of closing the gap in health and life expectancy between Aboriginal and Torres Strait Islander peoples and non-Indigenous Australians [54] .

Co-circulation of multidrug-resistant Shigella among men who have sex with men

Shigella flexneri 2b in the northern territory in 2017

Antimicrobial sensitivity trends and virulence genes in Shigella spp. from the Oceania region

Australian Government Department of Health

Hospital and Health Service maps (2020) Queensland Health

Inequalities in pediatric avoidable hospitalizations between Aboriginal and non-Aboriginal children in Australia: a population data linkage study

Housing conditions and health in Indigenous Australian communities: current status and recent trends

Disease burden and health-care clinic attendances for young children in remote aboriginal communities of northern Australia

CDNA National Guidelines for the Public Health Management of Trachoma (2014) Communicable diseases network Australia 11

Australia Commission of Safety and Quality in Health Care

Mobile genetic elements associated with antimicrobial resistance

Intercontinental dissemination of azithromycin-resistant shigellosis through sexual transmission: a cross-sectional study

Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world?

The phage-shock-protein response

The 3' conserved segment of integrons contains a gene associated with multidrug resistance to antiseptics and disinfectants

Dissemination of multidrug-resistant Shigella flexneri and Shigella sonnei with class 1, class 2, and atypical class 1 integrons in China

Prevalence, resistance patterns, and characterization of integrons of Shigella flexneri isolated from Jiangsu Province in China

Molecular characterization of Shigella spp. from patients in Gabon

Identification of Escherichia coli and Shigella species from whole genome sequences

Species-wide whole genome sequencing reveals historical global spread and recent local persistence in Shigella flexneri

Serotype-converting bacteriophages and O-antigen modification in Shigella flexneri

Queensland Government 25. Regional Population by Age and Sex

Breakpoint tables for interpretation of MICs and zone diameters (2019) The European Committee on Antimicrobial Susceptibility Testing

Trimmomatic: a flexible trimmer for Illumina sequence data

SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing

Snippy: rapid haploid variant calling and core SNP phylogeny

Identification of acquired antimicrobial resistance genes

Assembly of long, error-prone reads using repeat graphs

NextPolish: a fast and efficient genome polishing tool for long-read assembly

NCBI prokaryotic genome annotation pipeline

RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies

Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees

Easyfig: a genome comparison visualizer

Replicon sequence typing of IncF plasmids carrying virulence and resistance determinants

) qnrB, another plasmid-mediated gene for quinolone resistance

First clinical case of in vivo acquisition of DHA-1 plasmid-mediated AmpC in a Salmonella enterica subsp. enterica isolate

Characterization of cefotaxime-resistant urinary Escherichia coli from primary care in South-West England 2017-18

Resistance plasmid families in Enterobacteriaceae

Association of bla DHA-1 and qnrB genes carried by broadhost-range plasmids among isolates of Enterobacteriaceae at a Spanish hospital

A multidrug resistance plasmid pIMP26, carrying bla IMP-26 , fosA5, bla DHA-1 , and qnrB4 in Enterobacter cloacae

A degenerate primer MOB typing (DPMT) method to classify gammaproteobacterial plasmids in clinical and environmental settings

Complete nucleotide sequence of two multidrugresistant IncR plasmids from Klebsiella pneumoniae

Complete nucleotide sequence of pHN7A8, an F33: A-: B-type epidemic plasmid carrying bla CTX-M-65 , fosA3 and rmtB from China

Mosaic plasmids and mosaic replicons: evolutionary lessons from the analysis of genetic diversity in IncFII-related replicons

Prospective study of the association between serum antibodies to lipopolysaccharide O antigen and the attack rate of shigellosis

Clinical trials of Shigella vaccines: two steps forward and one step back on a long, hard road

Levels of Mycoplasma genitalium antimicrobial resistance differ by both region and gender in the state of Queensland, Australia: implications for treatment guidelines

CARAlert annual report: 2019 (2020) Australian Commission on Safety and Quality in Health Care

Antimicrobial resistance following mass azithromycin distribution for trachoma: a systematic review

Closing the gap and indigenous housing

Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

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

Acknowledgements This publication made use of the plasmid MLST website (https://pubmlst.org/plasmid/) sited at the University of Oxford [55] . The development of this site has been funded by the Wellcome Trust.Data availability All data generated or analysed during this study are included in this published article and its supplementary information files.