Microbiology & Immunology - Research Publications

Permanent URI for this collection

Search Results

Now showing 1 - 10 of 17
  • Item
    No Preview Available
    High-performance enrichment-based genome sequencing to support the investigation of hepatitis A virus outbreaks
    Zufan, SE ; Mercoulia, K ; Kwong, JC ; Judd, LM ; Howden, BP ; Seemann, T ; Stinear, TP ; Rantsiou, K (American Society for Microbiology, 2024-01)
    This proof-of-concept study introduces a hybrid capture oligo panel for whole-genome sequencing of all six human pathogenic hepatitis A virus (HAV) subgenotypes, exhibiting a higher sensitivity than some conventional genotyping assays. The ability of hybrid capture to enrich multiple targets allows for a single, streamlined workflow, thus facilitating the potential harmonization of molecular surveillance of HAV with other enteric viruses. Even challenging sample matrices can be accommodated, making them suitable for broad implementation in clinical and public health laboratories. This innovative approach has significant implications for enhancing multijurisdictional outbreak investigations as well as our understanding of the global diversity and transmission dynamics of HAV.
  • Item
    Thumbnail Image
    Optimising genomic approaches for identifying vancomycin-resistant Enterococcus faecium transmission in healthcare settings
    Higgs, C ; Sherry, NL ; Seemann, T ; Horan, K ; Walpola, H ; Kinsella, P ; Bond, K ; Williamson, DA ; Marshall, C ; Kwong, JC ; Grayson, ML ; Stinear, TP ; Gorrie, CL ; Howden, BP (NATURE PORTFOLIO, 2022-01-26)
    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.
  • Item
    Thumbnail Image
    The changing landscape of vancomycin-resistant Enterococcus faecium in Australia: a population-level genomic study
    Lee, RS ; da Silva, AG ; Baines, SL ; Strachan, J ; Ballard, S ; Carter, GP ; Kwong, JC ; Schultz, MB ; Bulach, DM ; Seemann, T ; Stinear, TP ; Howden, BP (OXFORD UNIV PRESS, 2018-12)
    BACKGROUND: Vancomycin-resistant Enterococcus faecium (VREfm) represent a major source of nosocomial infection worldwide. In Australia, there has been a recent concerning increase in bacteraemia associated with the vanA genotype, prompting investigation into the genomic epidemiology of VREfm. METHODS: A population-level study of VREfm (10 November-9 December 2015) was conducted. A total of 321 VREfm isolates (from 286 patients) across Victoria State were collected and sequenced with Illumina NextSeq. SNPs were used to assess relatedness. STs and genes associated with resistance and virulence were identified. The vanA-harbouring plasmid from an isolate from each ST was assembled using long-read data. Illumina reads from remaining isolates were then mapped to these assemblies to identify their probable vanA-harbouring plasmid. RESULTS: vanA-VREfm comprised 17.8% of isolates. ST203, ST80 and a pstS(-) clade, ST1421, predominated (30.5%, 30.5% and 37.2%, respectively). Most vanB-VREfm were ST796 (77.7%). vanA-VREfm were more closely related within hospitals versus between them [core SNPs 10 (IQR 1-357) versus 356 (179-416), respectively], suggesting discrete introductions of vanA-VREfm, with subsequent intra-hospital transmission. In contrast, vanB-VREfm had similar core SNP distributions within versus between hospitals, due to widespread dissemination of ST796. Different vanA-harbouring plasmids were found across STs. With the exception of ST78 and ST796, Tn1546 transposons also varied. Phylogenetic analysis revealed Australian strains were often interspersed with those from other countries, suggesting ongoing cross-continental transmission. CONCLUSIONS: Emerging vanA-VREfm in Australia is polyclonal, indicating repeat introductions of vanA-VREfm into hospitals and subsequent dissemination. The close relationship to global strains reinforces the need for ongoing screening and control of VREfm in Australia and abroad.
  • Item
    Thumbnail Image
    Systematic analysis of key parameters for genomics-based real-time detection and tracking of multidrug-resistant bacteria
    Gorrie, C ; Da Silva, AG ; Ingle, D ; Higgs, C ; Seemann, T ; Stinear, T ; Williamson, D ; Kwong, J ; Grayson, L ; Sherry, N ; Howden, B ( 2020-09-25)
    Background: Pairwise single nucleotide polymorphisms (SNPs) are a cornerstone for genomic approaches to multidrug-resistant organisms (MDROs) transmission inference in hospitals. However, the impact of key analysis parameters on these inferences has not been systematically analysed. Methods: We conducted a multi-hospital 15-month prospective study, sequencing 1537 MDRO genomes for comparison; methicillin-resistant Staphylococcus aureus , vancomycin-resistant Enterococcus faecium , and extended-spectrum beta-lactamase-producing Escherichia coli and Klebsiella pneumoniae . We systematically assessed the impact of sample and reference genome diversity, masking of prophage and regions of recombination, cumulative genome analysis compared to a three-month sliding-window, and the comparative effects each of these had when applying a SNP threshold for inferring likely transmission (≤15 SNPs for S. aureus , ≤25 for other species). Findings: Across the species, using a reference genome of the same sequence type provided a greater degree of pairwise SNP resolution, compared to species and outgroup-reference alignments that typically resulted in inflated SNP distances and the possibility of missed transmission events. Omitting prophage regions had minimal impacts, however, omitting recombination regions a highly variable effect, often inflating the number of closely related pairs. Estimating pairwise SNP distances was more consistent using a sliding-window than a cumulative approach. Interpretation: The use of a closely-related reference genome, without masking of prophage or recombination regions, and a sliding-window for isolate inclusion is best for accurate and consistent MDRO transmission inference. The increased stability and resolution provided by these approaches means SNP thresholds for putative transmission inference can be more reliably applied among diverse MDROs. Funding: This work was supported by the Melbourne Genomics Health Alliance (funded by the State Government of Victoria, Department of Health and Human Services, and the ten member organizations); an National Health and Medical Research Council (Australia) Partnership grant (GNT1149991) and individual grants from National Health and Medical Research Council (Australia) to NLS (GNT1093468), JCK (GNT1008549) and BPH (GNT1105905).
  • Item
    Thumbnail Image
    Key parameters for genomics-based real-time detection and tracking of multidrug-resistant bacteria: a systematic analysis
    Gorrie, CL ; Da Silva, AG ; Ingle, DJ ; Higgs, C ; Seemann, T ; Stinear, TP ; Williamson, DA ; Kwong, JC ; Grayson, ML ; Sherry, NL ; Howden, BP (ELSEVIER, 2021-11)
    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.
  • Item
    Thumbnail Image
    Comprehensive Genomic Investigation of Adaptive Mutations Driving the Low-Level Oxacillin Resistance Phenotype in Staphylococcus aureus
    Giulieri, SG ; Guerillot, R ; Kwong, JC ; Monk, IR ; Hayes, AS ; Daniel, D ; Baines, S ; Sherry, NL ; Holmes, NE ; Ward, P ; Gao, W ; Seemann, T ; Stinear, TP ; Howden, BP ; Pettigrew, MM (AMER SOC MICROBIOLOGY, 2020-12-08)
    Antistaphylococcal penicillins such as oxacillin are the key antibiotics in the treatment of invasive methicillin-susceptible Staphylococcus aureus (MSSA) infections; however, mec gene-independent resistance adaptation can cause treatment failure. Despite its clinical relevance, the basis of this phenomenon remains poorly understood. Here, we investigated the genomic adaptation to oxacillin at an unprecedented scale using a large collection of 503 clinical mec-negative isolates and 30 in vitro-adapted isolates from independent oxacillin exposures. By combining comparative genomics, evolutionary convergence, and genome-wide association analysis, we found 21 genetic loci associated with low-level oxacillin resistance, underscoring the polygenic nature of this phenotype. Evidence of adaptation was particularly strong for the c-di-AMP signal transduction pathways (gdpP and dacA) and in the clpXP chaperone-protease complex. The role of mutations in gdpP in conferring low-level oxacillin resistance was confirmed by allele-swapping experiments. We found that resistance to oxacillin emerges at high frequency in vitro (median, 2.9 × 10-6; interquartile range [IQR], 1.9 × 10-6 to 3.9 × 10-6), which is consistent with a recurrent minimum inhibitory concentration (MIC) increase across the global phylogeny of clinical isolates. Nevertheless, adaptation in clinical isolates appears sporadically, with no stably adapted lineages, suggesting a high fitness cost of resistance, confirmed by growth assessment of mutants in rich media. Our data provide a broader understanding of the emergence and dynamics of oxacillin resistance adaptation in S. aureus and a framework for future surveillance of this clinically important phenomenon.IMPORTANCE The majority of Staphylococcus aureus strains causing human disease are methicillin-susceptible (MSSA) and can be treated with antistaphylococcal penicillins (such as oxacillin). While acquisition of the mec gene represents the main resistance mechanism to oxacillin, S. aureus can acquire low-level resistance through adaptive mutations in other genes. In this study, we used genomic approaches to understand the basis of S. aureus adaption to oxacillin and its dynamic at the population level. By combining a genome analysis of clinical isolates from persistent MSSA infections, in vitro selection of oxacillin resistance, and genome-wide association analysis on a large collection of isolates, we identified 21 genes linked to secondary oxacillin resistance. Adaptive mutations in these genes were easy to select when S. aureus was exposed to oxacillin, but they also came at a substantial cost in terms of bacterial fitness, suggesting that this phenotype emerges preferentially in the setting of sustained antibiotic exposure.
  • Item
    Thumbnail Image
    Whole-genome sequencing reveals transmission of gonococcal antibiotic resistance among men who have sex with men: an observational study
    Kwong, JC ; Chow, EPF ; Stevens, K ; Stinear, TP ; Seemann, T ; Fairley, CK ; Chen, MY ; Howden, BP (BMJ PUBLISHING GROUP, 2018-03)
    OBJECTIVES: Drug-resistant Neisseria gonorrhoeae are now a global public health threat. Direct transmission of antibiotic-resistant gonococci between individuals has been proposed as a driver for the increased transmission of resistance, but direct evidence of such transmission is limited. Whole-genome sequencing (WGS) has superior resolution to investigate outbreaks and disease transmission compared with traditional molecular typing methods such as multilocus sequence typing (MLST) and N. gonorrhoeae multiantigen sequence (NG-MAST). We therefore aimed to systematically investigate the transmission of N. gonorrhoeae between men in sexual partnerships using WGS to compare isolates and their resistance to antibiotics at a genome level. METHODS: 458 couples from a large prospective cohort of men who have sex with men (MSM) tested for gonorrhoea together between 2005 and 2014 were included, and WGS was conducted on all isolates from couples where both men were culture-positive for N. gonorrhoeae. Resistance-determining sequences were identified from genome assemblies, and comparison of isolates between and within individuals was performed by pairwise single nucleotide polymorphism and pangenome comparisons, and in silico predictions of NG-MAST and MLST. RESULTS: For 33 of 34 (97%; 95% CI 85% to 100%) couples where both partners were positive for gonorrhoea, the resistance-determining genes and mutations were identical in isolates from each partner (94 isolates in total). Resistance determinants in isolates from 23 of 23 (100%; 95% CI 86% to 100%) men with multisite infections were also identical within an individual. These partner and within-host isolates were indistinguishable by NG-MAST, MLST and whole genomic comparisons. CONCLUSIONS: These data support the transmission of antibiotic-resistant strains between sexual partners as a key driver of resistance rates in gonorrhoea among MSM. This improved understanding of the transmission dynamics of N. gonorrhoeae between sexual partners will inform treatment and prevention guidelines.
  • Item
    Thumbnail Image
    Rapid Emergence and Evolution of Staphylococcus aureus Clones Harboring fusC-Containing Staphylococcal Cassette Chromosome Elements
    Baines, SL ; Howden, BP ; Heffernan, H ; Stinear, TP ; Carter, GP ; Seemann, T ; Kwong, JC ; Ritchie, SR ; Williamson, DA (AMER SOC MICROBIOLOGY, 2016-04)
    The prevalence of fusidic acid (FA) resistance amongStaphylococcus aureusstrains in New Zealand (NZ) is among the highest reported globally, with a recent study describing a resistance rate of approximately 28%. Three FA-resistantS. aureusclones (ST5 MRSA, ST1 MSSA, and ST1 MRSA) have emerged over the past decade and now predominate in NZ, and in all three clones FA resistance is mediated by thefusCgene. In particular, ST5 MRSA has rapidly become the dominant MRSA clone in NZ, although the origin of FA-resistant ST5 MRSA has not been explored, and the genetic context offusCin FA-resistant NZ isolates is unknown. To better understand the rapid emergence of FA-resistantS. aureus, we used population-based comparative genomics to characterize a collection of FA-resistant and FA-susceptible isolates from NZ. FA-resistant NZ ST5 MRSA displayed minimal genetic diversity and represented a phylogenetically distinct clade within a global population model of clonal complex 5 (CC5)S. aureus In all lineages,fusCwas invariably located within staphylococcal cassette chromosome (SCC) elements, suggesting that SCC-mediated horizontal transfer is the primary mechanism offusCdissemination. The genotypic association offusCwithmecAhas important implications for the emergence of MRSA clones in populations with high usage of fusidic acid. In addition, we found thatfusCwas colocated with a recently described virulence factor (tirS) in dominant NZS. aureusclones, suggesting a fitness advantage. This study points to the likely molecular mechanisms responsible for the successful emergence and spread of FA-resistantS. aureus.
  • Item
    Thumbnail Image
    A phylogenomic framework for assessing the global emergence and evolution of clonal complex 398 methicillin-resistant Staphylococcus aureus
    da Silva, AG ; Baines, SL ; Carter, GP ; Heffernan, H ; French, NP ; Ren, X ; Seemann, T ; Bulach, D ; Kwong, J ; Stinear, TP ; Howden, BP ; Williamson, DA (MICROBIOLOGY SOC, 2017-01)
    Distinct clones of methicillin-resistant Staphylococcus aureus (MRSA) have emerged as important causes of infection in individuals who have exposure to livestock (livestock-associated MRSA; LA-MRSA). Clonal complex 398 (CC398) is the most prevalent LA-MRSA clone, and has been reported from several geographical settings, including Europe, the Americas and Asia. To understand the factors contributing to the global dissemination of this clone, we analysed CC398 MRSA isolates from New Zealand (NZ), a geographically isolated country with an economy strongly dependent on livestock farming. We supplemented the NZ CC398 MRSA collection with global datasets of CC398 MRSA and CC398 methicillin-susceptible S. aureus. Here, we demonstrate multiple sporadic incursions of CC398 MRSA into NZ, as well as recent importation and spread of a swine-associated clade related to the European LA-MRSA lineage. Within a larger global phylogenomic framework, Bayesian modelling suggested that this NZ clade emerged in the late 2000s, with a probable origin in swine from Western Europe. Elucidating the factors responsible for the incursion and spread of LA-MRSA in geographically distant regions, such as NZ, provides important insights into global pathways of S. aureus transmission, and will inform strategies to control importation and spread.
  • Item
    Thumbnail Image
    Genomic epidemiology and antimicrobial resistance of Neisseria gonorrhoeae in New Zealand
    Lee, RS ; Seemann, T ; Heffernan, H ; Kwong, JC ; da Silva, AG ; Carter, GP ; Woodhouse, R ; Dyet, KH ; Bulach, DM ; Stinear, TP ; Howden, BP ; Williamson, DA (OXFORD UNIV PRESS, 2018-02)
    BACKGROUND: Antimicrobial-resistant Neisseria gonorrhoeae is a major threat to public health. No studies to date have examined the genomic epidemiology of gonorrhoea in the Western Pacific Region, where the incidence of gonorrhoea is particularly high. METHODS: A population-level study of N. gonorrhoeae in New Zealand (October 2014 to May 2015). Comprehensive susceptibility testing and WGS data were obtained for 398 isolates. Relatedness was inferred using phylogenetic trees, and pairwise core SNPs. Mutations and genes known to be associated with resistance were identified, and correlated with phenotype. RESULTS: Eleven clusters were identified. In six of these clusters, >25% of isolates were from females, while in eight of them, >15% of isolates were from females. Drug resistance was common; 98%, 32% and 68% of isolates were non-susceptible to penicillin, ciprofloxacin and tetracycline, respectively. Elevated MICs to extended-spectrum cephalosporins (ESCs) were observed in 3.5% of isolates (cefixime MICs ≥ 0.12 mg/L, ceftriaxone MICs ≥ 0.06 mg/L). Only nine isolates had penA XXXIV genotypes, three of which had decreased susceptibility to ESCs (MIC = 0.12 mg/L). Azithromycin non-susceptibility was identified in 43 isolates (10.8%); two of these isolates had 23S mutations (C2611T, 4/4 alleles), while all had mutations in mtrR or its promoter. CONCLUSIONS: The high proportion of females in clusters suggests transmission is not exclusively among MSM in New Zealand; re-assessment of risk factors for transmission may be warranted in this context. As elevated MICs of ESCs and/or azithromycin were found in closely related strains, targeted public health interventions to halt transmission are urgently needed.