Takehiro Tomita

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data.Source code is available from http://katholt.github.io/srst2/.

Prospective whole genome sequencing enhances national surveillance of Listeria monocytogenes

ABSTRACT Whole-genome sequencing (WGS) has emerged as a powerful tool for comparing bacterial isolates in outbreak detection and investigation. Here we demonstrate that WGS performed prospectively for national epidemiologic surveillance of Listeria monocytogenes has the capacity to be superior to our current approaches using pulsed-field gel electrophoresis (PFGE), multilocus sequence typing (MLST), multilocus variable-number tandem-repeat analysis (MLVA), binary typing, and serotyping. Initially 423 L. monocytogenes isolates underwent WGS, and comparisons uncovered a diverse genetic population structure derived from three distinct lineages. MLST, binary typing, and serotyping results inferred in silico from the WGS data were highly concordant (>99%) with laboratory typing performed in parallel. However, WGS was able to identify distinct nested clusters within groups of isolates that were otherwise indistinguishable using our current typing methods. Routine WGS was then used for prospective epidemiologic surveillance on a further 97 L. monocytogenes isolates over a 12-month period, which provided a greater level of discrimination than that of conventional typing for inferring linkage to point source outbreaks. A risk-based alert system based on WGS similarity was used to inform epidemiologists required to act on the data. Our experience shows that WGS can be adopted for prospective L. monocytogenes surveillance and investigated for other pathogens relevant to public health.

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

BACKGROUND Shigellosis is an acute, severe bacterial colitis that, in high-income countries, is typically associated with travel to high-risk regions (Africa, Asia, and Latin America). Since the 1970s, shigellosis has also been reported as a sexually transmitted infection in men who have sex with men (MSM), in whom transmission is an important component of shigellosis epidemiology in high-income nations. We aimed to use sophisticated subtyping and international sampling to determine factors driving shigellosis emergence in MSM linked to an outbreak in the UK. METHODS We did a large-scale, cross-sectional genomic epidemiological study of shigellosis cases collected from 29 countries between December, 1995, and June 8, 2014. Focusing on an ongoing epidemic in the UK, we collected and whole-genome sequenced clinical isolates of Shigella flexneri serotype 3a from high-risk and low-risk regions, including cases associated with travel and sex between men. We examined relationships between geographical, demographic, and clinical patient data with the isolate antimicrobial susceptibility, genetic data, and inferred evolutionary relationships. FINDINGS We obtained 331 clinical isolates of S flexneri serotype 3a, including 275 from low-risk regions (44 from individuals who travelled to high-risk regions), 52 from high-risk regions, and four outgroup samples (ie, closely related, but genetically distinct isolates used to determine the root of the phylogenetic tree). We identified a recently emerged lineage of S flexneri 3a that has spread intercontinentally in less than 20 years throughout regions traditionally at low risk for shigellosis via sexual transmission in MSM. The lineage had acquired multiple antimicrobial resistance determinants, and prevailing sublineages were strongly associated with resistance to the macrolide azithromycin. Eight (4%) of 206 isolates from the MSM-associated lineage were obtained from patients who had previously provided an isolate; these serial isolations indicated atypical infection patterns (eg, reinfection). INTERPRETATION We identified transmission-facilitating behaviours and atypical course(s) of infection as precipitating factors in shigellosis-affected MSM. The intercontinental spread of antimicrobial-resistant shigella through established transmission routes emphasises the need for new approaches to tackle the public health challenge of sexually transmitted infections in MSM. FUNDING Wellcome Trust (grant number 098051).

Failure of 500 mg of ceftriaxone to eradicate pharyngeal gonorrhoea, Australia

been occurring since admission. On day 14, the patient developed leucocytosis and blood cultures were drawn, which grew MRSA. The MIC of vancomycin increased to 2 mg/L, with that of daptomycin remaining at 0.38 mg/L. Since the patient continued to be symptomatic, irrigation and debridement was scheduled, but due to surgeon availability issues this was not performed until hospital day 19. A sample was obtained during the procedure and cultured, and on day 22, after .2 weeks of treatment with daptomycin, MRSA grew with a daptomycin MIC of 3 mg/L [daptomycin-non-susceptible S. aureus (DNSSA)]. The primary cause of the daptomycin failure was most likely a delay in incision and drainage of the abscesses. Experimental osteomyelitis models have shown selective resistance with daptomycin monotherapy, 2 which possibly contributed to the resistance in this case. Daptomycin was discontinued and, unfortunately, vancomycin was restarted over the weekend. After the patient became bacteraemic on day 14, we recommended a TEE, which was ordered but not performed until day 23. Lesions were found on the left coronary cusp (LCC) and the tricuspid valve, consistent with endocarditis. Vancomycin was discontinued, since current expert opinion states alternative anti-microbials should be considered when the vancomycin MIC is ≥2 mg/L, due to treatment failures. 3,4 Considering this and the severity of infection in this patient, we turned to salvage therapy options and ceftaroline (600 mg intravenously every 12 h), which had an Etest MIC of 0.5 mg/L, was chosen. Ceftaroline is a fifth-generation cephalosporin with activity against MRSA. Ceftaroline was chosen for a number of reasons. First, it is rapidly bactericidal. Ho et al. 5 reported sterilization in 13 days in a case of endocarditis. Also, the incidence of thrombo-cytopenia with long-term use is low, which was particularly important since our patient had a history of thrombocytopenia secondary to hepatitis. Finally, a literature search revealed a series of case reports showing successful use of ceftaroline in the treatment of MRSA endocarditis 5 and an experimental model demonstrating potential use in osteomyelitis. 6 Subsequent blood cultures drawn on days 29 and 32 were negative. On day 32, a repeat irrigation and debridement was performed. After 37 consecutive days of ceftaroline treatment, a TEE showed resolution of echodensities at the base of the tri-cuspid valve and a stable, fibronodular lesion on the LCC. Ceftaro-line was continued for a total of 44 days, with no further episodes of bacteraemia, leucocytosis or fever, …

Identification of Streptococcus uberis Multilocus Sequence Types Highly Associated with Mastitis

ABSTRACT Multilocus sequence typing analysis of Streptococcus uberis has identified a cluster of isolates associated with clinical and subclinical mastitis and a cluster associated with cows with low somatic cell counts in their milk. Specific groups of genotypes (global clonal complex [GCC] sequence type 5s [ST5s] and GCC ST143s) were highly associated (P = 0.006) with clinical and subclinical mastitis and may represent a lineage of virulent isolates, whereas isolates belonging to GCC ST86 were associated with low-cell-count cows. This study has, for the first time, demonstrated the occurrence of identical sequence types (ST60 and ST184) between different continents (Australasia and Europe) and different countries (Australia and New Zealand). The standardized index of association and the empirical estimation of the rate of recombination showed substantial recombination within the S. uberis population in Australia, consistent with previous multilocus sequence type analyses.

Translating genomics into practice for real-time surveillance and response to carbapenemase-producing Enterobacteriaceae: evidence from a complex multi-institutional KPC outbreak

Background Until recently, Klebsiella pneumoniae carbapenemase (KPC)-producing Enterobacteriaceae were rarely identified in Australia. Following an increase in the number of incident cases across the state of Victoria, we undertook a real-time combined genomic and epidemiological investigation. The scope of this study included identifying risk factors and routes of transmission, and investigating the utility of genomics to enhance traditional field epidemiology for informing management of established widespread outbreaks. Methods All KPC-producing Enterobacteriaceae isolates referred to the state reference laboratory from 2012 onwards were included. Whole-genome sequencing was performed in parallel with a detailed descriptive epidemiological investigation of each case, using Illumina sequencing on each isolate. This was complemented with PacBio long-read sequencing on selected isolates to establish high-quality reference sequences and interrogate characteristics of KPC-encoding plasmids. Results Initial investigations indicated that the outbreak was widespread, with 86 KPC-producing Enterobacteriaceae isolates (K. pneumoniae 92%) identified from 35 different locations across metropolitan and rural Victoria between 2012 and 2015. Initial combined analyses of the epidemiological and genomic data resolved the outbreak into distinct nosocomial transmission networks, and identified healthcare facilities at the epicentre of KPC transmission. New cases were assigned to transmission networks in real-time, allowing focussed infection control efforts. PacBio sequencing confirmed a secondary transmission network arising from inter-species plasmid transmission. Insights from Bayesian transmission inference and analyses of within-host diversity informed the development of state-wide public health and infection control guidelines, including interventions such as an intensive approach to screening contacts following new case detection to minimise unrecognised colonisation. Conclusion A real-time combined epidemiological and genomic investigation proved critical to identifying and defining multiple transmission networks of KPC Enterobacteriaceae, while data from either investigation alone were inconclusive. The investigation was fundamental to informing infection control measures in real-time and the development of state-wide public health guidelines on carbapenemase-producing Enterobacteriaceae surveillance and management.

Comparative Genomics Shows That Mycobacterium ulcerans Migration and Expansion Preceded the Rise of Buruli Ulcer in Southeastern Australia

ABSTRACT Since 2000, cases of the neglected tropical disease Buruli ulcer, caused by infection with Mycobacterium ulcerans, have increased 100-fold around Melbourne (population 4.4 million), the capital of Victoria, in temperate southeastern Australia. The reasons for this increase are unclear. Here, we used whole-genome sequence comparisons of 178 M. ulcerans isolates obtained primarily from human clinical specimens, spanning 70 years, to model the population dynamics of this pathogen from this region. Using phylogeographic and advanced Bayesian phylogenetic approaches, we found that there has been a migration of the pathogen from the east end of the state, beginning in the 1980s, 300 km west to the major human population center around Melbourne. This move was then followed by a significant increase in M. ulcerans population size. These analyses inform our thinking around Buruli ulcer transmission and control, indicating that M. ulcerans is introduced to a new environment and then expands, rather than it being from the awakening of a quiescent pathogen reservoir. IMPORTANCE Buruli ulcer is a destructive skin and soft tissue infection caused by Mycobacterium ulcerans and is characterized by progressive skin ulceration, which can lead to permanent disfigurement and long-term disability. Despite the majority of disease burden occurring in regions of West and central Africa, Buruli ulcer is also becoming increasingly common in southeastern Australia. Major impediments to controlling disease spread are incomplete understandings of the environmental reservoirs and modes of transmission of M. ulcerans. The significance of our research is that we used genomics to assess the population structure of this pathogen at the Australian continental scale. We have then reconstructed a historical bacterial spread and modeled demographic dynamics to reveal bacterial population expansion across southeastern Australia. These findings provide explanations for the observed epidemiological trends with Buruli ulcer and suggest possible management to control disease spread.

Real-time genomic and epidemiological investigation of a multi-institution outbreak of KPC-producing Enterobacteriaceae: a translational study

Background Until recently, KPC-producing Enterobacteriaceae were rarely identified in Australia. Following an increase in the number of incident cases across the state of Victoria, we undertook a real-time combined genomic and epidemiological investigation. The scope of this study included identifying risk factors and routes of transmission, and investigating the utility of genomics to enhance traditional field epidemiology for informing management of established widespread outbreaks. Methods and Findings All KPC-producing Enterobacteriaceae isolates referred to the state reference laboratory from 2012 onwards were included. Whole-genome sequencing (WGS) was performed in parallel with a detailed descriptive epidemiological investigation of each case, using Illumina sequencing on each isolate. This was complemented with PacBio long-read sequencing on selected isolates to establish high-quality reference sequences and interrogate characteristics of KPC-encoding plasmids. Initial investigations indicated the outbreak was widespread, with 86 KPC-producing Enterobacteriaceae isolates (K. pneumoniae 92%) identified from 35 different locations across metropolitan and rural Victoria between 2012-2015. Initial combined analyses of the epidemiological and genomic data resolved the outbreak into distinct nosocomial transmission networks, and identified healthcare facilities at the epicentre of KPC transmission. New cases were assigned to transmission networks in real-time, allowing focussed infection control efforts. PacBio sequencing confirmed a secondary transmission network arising from inter-species plasmid transmission. Insights from Bayesian transmission inference and analyses of within-host diversity informed the development of state-wide public health and infection control guidelines, including interventions such as an intensive approach to screening contacts following new case detection to minimise unrecognised colonisation. Conclusions A real-time combined epidemiological and genomic investigation proved critical to identifying and defining multiple transmission networks of KPC Enterobacteriaceae, while data from either investigation alone were inconclusive. The investigation was fundamental to informing infection control measures in real-time and the development of state-wide public health guidelines on carbapenemase producing Enterobacteriaceae management.

Convergent Evolution Driven by Rifampin Exacerbates the Global Burden of Drug-Resistant Staphylococcus aureus

Increasing antibiotic resistance in the major human pathogen Staphylococcus aureus is threatening the ability to treat patients with these infections. Recent laboratory studies suggest that mutations in the gene commonly associated with rifampin resistance may also impact susceptibility to other last-line antibiotics in S. aureus; however, the overall frequency and clinical impact of these mutations are unknown. By mining a global collection of clinical S. aureus genomes and by mutagenesis experiments, this work reveals that common rifampin-induced rpoB mutations promote phenotypic plasticity that has led to the global emergence of stable, multidrug-resistant S. aureus lineages that are associated with increased risk of therapeutic failure through coresistance to other last-line antimicrobials. We recommend decreasing susceptibility breakpoints for rifampin to allow phenotypic detection of critical rpoB mutations conferring low resistance to rifampin and reconsidering the appropriate use of rifampin to reduce the fixation and spread of these deleterious mutations globally. ABSTRACT Mutations in the beta-subunit of bacterial RNA polymerase (RpoB) cause resistance to rifampin (Rifr), a critical antibiotic for treatment of multidrug-resistant Staphylococcus aureus. In vitro studies have shown that RpoB mutations confer decreased susceptibility to other antibiotics, but the clinical relevance is unknown. Here, by analyzing 7,099 S. aureus genomes, we demonstrate that the most prevalent RpoB mutations promote clinically relevant phenotypic plasticity resulting in the emergence of stable S. aureus lineages, associated with increased risk of therapeutic failure through generation of small-colony variants (SCVs) and coresistance to last-line antimicrobial agents. We found eight RpoB mutations that accounted for 93% (469/505) of the total number of Rifr mutations. The most frequently selected amino acid substitutions affecting residue 481 (H481N/Y) were associated with worldwide expansions of Rifr clones spanning decades. Recreating the H481N/Y mutations confirmed no impact on S. aureus growth, but the H481N mutation promoted the emergence of a subpopulation of stable Rifr SCVs with reduced susceptibility to vancomycin and daptomycin. Recreating the other frequent RpoB mutations showed similar impacts on resistance to these last-line agents. We found that 86% of all Rifr isolates in our global sample carried the mutations promoting cross-resistance to vancomycin and 52% to both vancomycin and daptomycin. As four of the most frequent RpoB mutations confer only low-level Rifr, equal to or below some international breakpoints, we recommend decreasing these breakpoints and reconsidering the appropriate use of rifampin to reduce the fixation and spread of these clinically deleterious mutations. IMPORTANCE Increasing antibiotic resistance in the major human pathogen Staphylococcus aureus is threatening the ability to treat patients with these infections. Recent laboratory studies suggest that mutations in the gene commonly associated with rifampin resistance may also impact susceptibility to other last-line antibiotics in S. aureus; however, the overall frequency and clinical impact of these mutations are unknown. By mining a global collection of clinical S. aureus genomes and by mutagenesis experiments, this work reveals that common rifampin-induced rpoB mutations promote phenotypic plasticity that has led to the global emergence of stable, multidrug-resistant S. aureus lineages that are associated with increased risk of therapeutic failure through coresistance to other last-line antimicrobials. We recommend decreasing susceptibility breakpoints for rifampin to allow phenotypic detection of critical rpoB mutations conferring low resistance to rifampin and reconsidering the appropriate use of rifampin to reduce the fixation and spread of these deleterious mutations globally.

Comprehensive antibiotic-linked mutation assessment by resistance mutation sequencing (RM-seq)

Mutation acquisition is a major mechanism of bacterial antibiotic resistance that remains insufficiently characterised. Here we present RM-seq, a new amplicon-based deep sequencing workflow based on a molecular barcoding technique adapted from Low Error Amplicon sequencing (LEA-seq). RM-seq allows detection and functional assessment of mutational resistance at high throughput from mixed bacterial populations. The sensitive detection of very low-frequency resistant sub-populations permits characterisation of antibiotic-linked mutational repertoires in vitro and detection of rare resistant populations during infections. Accurate quantification of resistance mutations enables phenotypic screening of mutations conferring pleiotropic phenotypes such as in vivo persistence, collateral sensitivity or cross-resistance. RM-seq will facilitate comprehensive detection, characterisation and surveillance of resistant bacterial populations (https://github.com/rguerillot/RM-seq).