Comparison of contemporary invasive and non-invasive Streptococcus pneumoniae isolates reveals new insights into circulating anti-microbial resistance determinants.

Streptococcus pneumoniae is a major human pathogen with a high burden of disease. Non-invasive isolates (those found in non-sterile sites) are thought to be a key source of invasive isolates (those found in sterile sites) and a reservoir of anti-microbial resistance (AMR) determinants. Despite this, pneumococcal surveillance has almost exclusively focused on invasive isolates. We aimed to compare contemporaneous invasive and non-invasive isolate populations to understand how they interact and identify differences in AMR gene distribution. We used a combination of whole-genome sequencing and phenotypic anti-microbial susceptibility testing and a data set of invasive (n = 1,288) and non-invasive (n = 186) pneumococcal isolates, collected in Victoria, Australia, between 2018 and 2022. The non-invasive population had increased levels of antibiotic resistance to multiple classes of antibiotics including beta-lactam antibiotics penicillin and ceftriaxone. We identified genomic intersections between the invasive and non-invasive populations and no distinct phylogenetic clustering of the two populations. However, this analysis revealed sub-populations overrepresented in each population. The sub-populations that had high levels of AMR were overrepresented in the non-invasive population. We determined that WamR-Pneumo was the most accurate in silico tool for predicting resistance to the antibiotics tested. This tool was then used to assess the allelic diversity of the penicillin-binding protein genes, which acquire mutations leading to beta-lactam antibiotic resistance, and found that they were highly conserved (≥80% shared) between the two populations. These findings show the potential of non-invasive isolates to serve as reservoirs of AMR determinants.

Genomic investigation of Staphylococcus aureus recovered from Gambian women and newborns following an oral dose of intra-partum azithromycin.

BACKGROUND: Oral azithromycin given during labour reduces carriage of bacteria responsible for neonatal sepsis, including Staphylococcus aureus. However, there is concern that this may promote drug resistance. OBJECTIVES: Here, we combine genomic and epidemiological data on S. aureus isolated from mothers and babies in a randomized intra-partum azithromycin trial (PregnAnZI) to describe bacterial population dynamics and resistance mechanisms. METHODS: Participants from both arms of the trial, who carried S. aureus in day 3 and day 28 samples post-intervention, were included. Sixty-six S. aureus isolates (from 7 mothers and 10 babies) underwent comparative genome analyses and the data were then combined with epidemiological data. Trial registration (main trial): ClinicalTrials.gov Identifier NCT01800942. RESULTS: Seven S. aureus STs were identified, with ST5 dominant (n = 40, 61.0%), followed by ST15 (n = 11, 17.0%). ST5 predominated in the placebo arm (73.0% versus 49.0%, P = 0.039) and ST15 in the azithromycin arm (27.0% versus 6.0%, P = 0.022). In azithromycin-resistant isolates, msr(A) was the main macrolide resistance gene (n = 36, 80%). Ten study participants, from both trial arms, acquired azithromycin-resistant S. aureus after initially harbouring a susceptible isolate. In nine (90%) of these cases, the acquired clone was an msr(A)-containing ST5 S. aureus. Long-read sequencing demonstrated that in ST5, msr(A) was found on an MDR plasmid. CONCLUSIONS: Our data reveal in this Gambian population the presence of a dominant clone of S. aureus harbouring plasmid-encoded azithromycin resistance, which was acquired by participants in both arms of the study. Understanding these resistance dynamics is crucial to defining the public health drug resistance impacts of azithromycin prophylaxis given during labour in Africa.

Longitudinal genomic analysis of Neisseria gonorrhoeae transmission dynamics in Australia.

N. gonorrhoeae, which causes the sexually transmissible infection gonorrhoea, remains a significant public health threat globally, with challenges posed by increasing transmission and antimicrobial resistance (AMR). The COVID-19 pandemic introduced exceptional circumstances into communicable disease control, impacting the transmission of gonorrhoea and other infectious diseases. Through phylogenomic and phylodynamic analysis of 5881 N. gonorrhoeae genomes from Australia, we investigated N. gonorrhoeae transmission over five years, including a time period during the COVID-19 pandemic. Using a novel cgMLST-based genetic threshold, we demonstrate persistence of large N. gonorrhoeae genomic clusters over several years, with some persistent clusters associated with heterosexual transmission. We observed a decline in both N. gonorrhoeae transmission and genomic diversity during the COVID-19 pandemic, suggestive of an evolutionary bottleneck. The longitudinal, occult transmission of N. gonorrhoeae over many years further highlights the urgent need for improved diagnostic, treatment, and prevention strategies for gonorrhoea.

Recombinational exchange of M-fibril and T-pilus genes generates extensive cell surface diversity in the global group A Streptococcus population.

Among genes present in all group A streptococci (GAS), those encoding M-fibril and T-pilus proteins display the highest levels of sequence diversity, giving rise to the two primary serological typing schemes historically used to define strain. A new genotyping scheme for the pilin adhesin and backbone genes is developed and, when combined with emm typing, provides an account of the global GAS strain population. Cluster analysis based on nucleotide sequence similarity assigns most T-serotypes to discrete pilin backbone sequence clusters, yet the established T-types correspond to only half the clusters. The major pilin adhesin and backbone sequence clusters yield 98 unique combinations, defined as "pilin types." Numerous horizontal transfer events that involve pilin or emm genes generate extensive antigenic and functional diversity on the bacterial cell surface and lead to the emergence of new strains. Inferred pilin genotypes applied to a meta-analysis of global population-based collections of pharyngitis and impetigo isolates reveal highly significant associations between pilin genotypes and GAS infection at distinct ecological niches, consistent with a role for pilin gene products in adaptive evolution. Integration of emm and pilin typing into open-access online tools (pubmlst.org) ensures broad utility for end-users wanting to determine the architecture of M-fibril and T-pilus genes from genome assemblies.IMPORTANCEPrecision in defining the variant forms of infectious agents is critical to understanding their population biology and the epidemiology of associated diseases. Group A Streptococcus (GAS) is a global pathogen that causes a wide range of diseases and displays a highly diverse cell surface due to the antigenic heterogeneity of M-fibril and T-pilus proteins which also act as virulence factors of varied functions. emm genotyping is well-established and highly utilized, but there is no counterpart for pilin genes. A global GAS collection provides the basis for a comprehensive pilin typing scheme, and online tools for determining emm and pilin genotypes are developed. Application of these tools reveals the expansion of structural-functional diversity among GAS via horizontal gene transfer, as evidenced by unique combinations of surface protein genes. Pilin and emm genotype correlations with superficial throat vs skin infection provide new insights on the molecular determinants underlying key ecological and epidemiological trends.

Population structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in Victoria, Australia.

Streptococcus pneumoniae is a major human pathogen and can cause a range of conditions from asymptomatic colonization to invasive pneumococcal disease (IPD). The epidemiology and distribution of IPD-causing serotypes in Australia has undergone large changes following the introduction of the 7-valent pneumococcal conjugate vaccine (PCV) in 2005 and the 13-valent PCV in 2011. In this study, to provide a contemporary understanding of the IPD causing population in Victoria, Australia, we aimed to examine the population structure and prevalence of antimicrobial resistance using whole-genome sequencing and comprehensive antimicrobial susceptibility data of 1288 isolates collected between 2018 and 2022. We observed high diversity among the isolates with 52 serotypes, 203 sequence types (STs) and 70 Global Pneumococcal Sequencing Project Clusters (GPSCs) identified. Serotypes contained in the 13v-PCV represented 35.3 % (n=405) of isolates. Antimicrobial resistance (AMR) to at least one antibiotic was identified in 23.8 % (n=358) of isolates with penicillin resistance the most prevalent (20.3 %, n=261 using meningitis breakpoints and 5.1 % n=65 using oral breakpoints). Of the AMR isolates, 28 % (n=101) were multidrug resistant (MDR) (resistant to three or more drug classes). Vaccination status of cases was determined for a subset of isolates with 34 cases classified as vaccine failure events (fully vaccinated IPD cases of vaccine serotype). However, no phylogenetic association with failure events was observed. Within the highly diverse IPD population, we identified six high-risk sub-populations of public health concern characterized by high prevalence, high rates of AMR and MDR, or serotype inclusion in vaccines. High-risk serotypes included serotypes 3, 19F, 19A, 14, 11A, 15A and serofamily 23. In addition, we present our data validating seroBA for in silico serotyping to facilitate ISO-accreditation of this test in routine use in a public health reference laboratory and have made this data set available. This study provides insights into the population dynamics, highlights non-vaccine serotypes of concern that are highly resistant, and provides a genomic framework for the ongoing surveillance of IPD in Australia which can inform next-generation IPD prevention strategies.

Optimising genomic approaches for identifying vancomycin-resistant Enterococcus faecium transmission in healthcare settings.

Vancomycin-resistant Enterococcus faecium (VREfm) is a major nosocomial pathogen. Identifying VREfm transmission dynamics permits targeted interventions, and while genomics is increasingly being utilised, methods are not yet standardised or optimised for accuracy. We aimed to develop a standardized genomic method for identifying putative VREfm transmission links. Using comprehensive genomic and epidemiological data from a cohort of 308 VREfm infection or colonization cases, we compared multiple approaches for quantifying genetic relatedness. We showed that clustering by core genome multilocus sequence type (cgMLST) was more informative of population structure than traditional MLST. Pairwise genome comparisons using split k-mer analysis (SKA) provided the high-level resolution needed to infer patient-to-patient transmission. The more common mapping to a reference genome was not sufficiently discriminatory, defining more than three times more genomic transmission events than SKA (3729 compared to 1079 events). Here, we show a standardized genomic framework for inferring VREfm transmission that can be the basis for global deployment of VREfm genomics into routine outbreak detection and investigation.

Multi-site implementation of whole genome sequencing for hospital infection control: A prospective genomic epidemiological analysis.

Background: Current microbiological methods lack the resolution to accurately identify multidrug-resistant organism (MDRO) transmission, however, whole genome sequencing can identify highly-related patient isolates providing opportunities for precision infection control interventions. We investigated the feasibility and potential impact of a prospective multi-centre genomics workflow for hospital infection control. Methods: We conducted a prospective genomics implementation study across eight Australian hospitals over 15 months (2017,2018), collecting all clinical and screening isolates from inpatients with vanA VRE, MRSA, ESBL Escherichia coli (ESBL-Ec), or ESBL Klebsiella pneumoniae (ESBL-Kp). Genomic and epidemiologic data were integrated to assess MDRO transmission. Findings: In total, 2275 isolates were included from 1970 patients, predominantly ESBL-Ec (40·8%) followed by MRSA (35·6%), vanA VRE (15·2%), and ESBL-Kp (8·3%).Overall, hospital and genomic epidemiology showed 607 patients (30·8%) acquired their MDRO in hospital, including the majority of vanA VRE (266 patients, 86·4%), with lower proportions of ESBL-Ec (186 patients, 23·0%), ESBL-Kp (42 patients, 26·3%), and MRSA (113 patients, 16·3%). Complex patient movements meant the majority of MDRO transmissions would remain undetected without genomic data.The genomics implementation had major impacts, identifying unexpected MDRO transmissions prompting new infection control interventions, and contributing to vanA VRE becoming a notifiable condition. We identified barriers to implementation and recommend strategies for mitigation. Interpretation: Implementation of a multi-centre genomics-informed infection control workflow is feasible and identifies many unrecognised MDRO transmissions. This provides critical opportunities for interventions to improve patient safety in hospitals. Funding: Melbourne Genomics Health Alliance (supported by State Government of Victoria, Australia), and National Health and Medical Research Council (Australia).

Key parameters for genomics-based real-time detection and tracking of multidrug-resistant bacteria: a systematic analysis.

BACKGROUND: Pairwise single nucleotide polymorphisms (SNPs) are a cornerstone of genomic approaches to the inference of transmission of multidrug-resistant (MDR) organisms in hospitals. However, the impact of many key analytical approaches on these inferences has not yet been systematically assessed. This study aims to make such a systematic assessment. METHODS: We conducted a 15-month prospective study (2-month pilot phase, 13-month implementation phase), across four hospital networks including eight hospitals in Melbourne, VIC, Australia. Patient clinical and screening samples containing one or more isolates of meticillin-resistant Staphylococcus aureus, vancomycin-resistant Enterococcus faecium, and extended-spectrum ß-lactamase-producing Escherichia coli and Klebsiella pneumoniae were collected and underwent whole genome sequencing. Using the genome data from the top four most numerous sequence types from each species, 16 in total, we systematically assessed the: (1) impact of sample and reference genome diversity through multiple core genome alignments using different data subsets and reference genomes, (2) effect of masking of prophage and regions of recombination in the core genome alignments by assessing SNP distances before and after masking, (3) differences between a cumulative versus a 3-month sliding-window approach to sample genome inclusion in the dataset over time, and (4) the comparative effects each of these approaches had when applying a previously defined SNP threshold for inferring likely transmission. FINDINGS: 2275 samples were collected (397 during the pilot phase from April 4 to June 18, 2017; 1878 during the implementation phase from Oct 30, 2017, to Nov 30, 2018) from 1870 patients. Of these 2275 samples, 1537 were identified as arising from the four most numerous sequence types from each of the four target species of MDR organisms in this dataset (16 sequence types in total: S aureus ST5, ST22, ST45, and ST93; E faecium ST80, ST203, ST1421, and ST1424; K pneumoniae ST15, ST17, ST307, and ST323; and E coli ST38, ST131, ST648, and ST1193). Across the species, using a reference genome of the same sequence type provided a greater degree of pairwise SNP resolution, compared with species and outgroup-reference alignments that mostly resulted in inflated SNP distances and the possibility of missed transmission events. Omitting prophage regions had minimal effect; however, omitting recombination regions had a highly variable effect, often inflating the number of closely related pairs. Estimated SNP distances between isolate pairs over time were more consistent using a sliding-window than a cumulative approach. INTERPRETATION: We propose that the use of a closely related reference genome, without masking of prophage or recombination regions, and of a sliding-window approach for isolate inclusion is best for accurate and consistent MDR organism transmission inference, when using core genome alignments and SNP thresholds. These approaches provide increased stability and resolution, so SNP thresholds can be more reliably applied for putative transmission inference among diverse MDR organisms, reducing the chance of incorrectly inferring the presence or absence of close genetic relatedness and, therefore, transmission. The establishment of a broadly applicable and standardised approach, as proposed here, is necessary to implement widespread prospective genomic surveillance for MDR organism transmission. FUNDING: Melbourne Genomics Health Alliance, and National Health and Medical Research Council of Australia.

Daptomycin Resistance Occurs Predominantly in vanA-Type Vancomycin-Resistant Enterococcus faecium in Australasia and Is Associated With Heterogeneous and Novel Mutations.

Healthcare associated infections caused by vancomycin-resistant Enterococcus faecium (VREfm) have a major impact on health outcomes. VREfm is difficult to treat because of intrinsic and acquired resistance to many clinically used antimicrobials, with daptomycin being one of the few last line therapeutic options for treating multidrug-resistant VREfm. The emergence of daptomycin-resistant VREfm is therefore of serious clinical concern. Despite this, the impact that daptomycin-resistant VREfm have on patient health outcomes is not clearly defined and knowledge on the mechanisms and genetic signatures linked with daptomycin resistance in VREfm remains incomplete. To address these knowledge gaps, phenotypic daptomycin susceptibility testing was undertaken on 324 E. faecium isolates from Australia and New Zealand. Approximately 15% of study isolates were phenotypically resistant to daptomycin. Whole genome sequencing revealed a strong association between vanA-VREfm and daptomycin resistance, with 95% of daptomycin-resistant study isolates harbouring vanA. Genomic analyses showed that daptomycin-resistant VREfm isolates were polyclonal and carried several previously characterised mutations in the liaR and liaS genes as well as several novel mutations within the rpoB, rpoC, and dltC genes. Overall, 70% of daptomycin-resistant study isolates were found to carry mutations within the liaR, rpoB, rpoC, or dltC genes. Finally, in a mouse model of VREfm bacteraemia, infection with the locally dominant daptomycin-resistant clone led to reduced daptomycin treatment efficacy in comparison to daptomycin-susceptible E. faecium. These findings have important implications for ongoing VREfm surveillance activities and the treatment of VREfm infections.

Pilot study of a combined genomic and epidemiologic surveillance program for hospital-acquired multidrug-resistant pathogens across multiple hospital networks in Australia.

OBJECTIVES: To conduct a pilot study implementing combined genomic and epidemiologic surveillance for hospital-acquired multidrug-resistant organisms (MDROs) to predict transmission between patients and to estimate the local burden of MDRO transmission. DESIGN: Pilot prospective multicenter surveillance study. SETTING: The study was conducted in 8 university hospitals (2,800 beds total) in Melbourne, Australia (population 4.8 million), including 4 acute-care, 1 specialist cancer care, and 3 subacute-care hospitals. METHODS: All clinical and screening isolates from hospital inpatients (April 24 to June 18, 2017) were collected for 6 MDROs: vanA VRE, MRSA, ESBL Escherichia coli (ESBL-Ec) and Klebsiella pneumoniae (ESBL-Kp), and carbapenem-resistant Pseudomonas aeruginosa (CRPa) and Acinetobacter baumannii (CRAb). Isolates were analyzed and reported as routine by hospital laboratories, underwent whole-genome sequencing at the central laboratory, and were analyzed using open-source bioinformatic tools. MDRO burden and transmission were assessed using combined genomic and epidemiologic data. RESULTS: In total, 408 isolates were collected from 358 patients; 47.5% were screening isolates. ESBL-Ec was most common (52.5%), then MRSA (21.6%), vanA VRE (15.7%), and ESBL-Kp (7.6%). Most MDROs (88.3%) were isolated from patients with recent healthcare exposure.Combining genomics and epidemiology identified that at least 27.1% of MDROs were likely acquired in a hospital; most of these transmission events would not have been detected without genomics. The highest proportion of transmission occurred with vanA VRE (88.4% of patients). CONCLUSIONS: Genomic and epidemiologic data from multiple institutions can feasibly be combined prospectively, providing substantial insights into the burden and distribution of MDROs, including in-hospital transmission. This analysis enables infection control teams to target interventions more effectively.

Genomics of vancomycin-resistant Enterococcus faecium.

Vancomycin-resistant Enterococcus faecium (VREfm) is a globally significant public health threat and was listed on the World Health Organization's 2017 list of high-priority pathogens for which new treatments are urgently needed. Treatment options for invasive VREfm infections are very limited, and outcomes are often poor. Whole-genome sequencing is providing important new insights into VREfm evolution, drug resistance and hospital adaptation, and is increasingly being used to track VREfm transmission within hospitals to detect outbreaks and inform infection control practices. This mini-review provides an overview of recent data on the use of genomics to understand and respond to the global problem of VREfm.